Two-dimensional mantle convection simulations using an internal state variable model: The role of a history dependent rheology on mantle convection

Abstract

We apply the Bammann inelastic internal state variable model (BIISV) to a mantle convection code TERRA2D to investigate the influence of a history dependent solid mechanics model on mantle convection. We compare and contrast the general purpose BIISV model to the commonly used power-law model. We implemented the BIISV model using a radial return algorithm and tested it against previously published mantle convection simulation results for verification. Model constants for the BIISV are used based on experimental stress-strain behaviour found in the literature. After implementation we give illustrative simulation examples were the BIISV produces hardened areas on the cold thermal boundary layer that the power-law model cannot produce. The hardened boundary layers divert material downward giving a plausible reason for the current subduction zones that are present on the Earth. © 2011 RAS.