Robust characteristics method for modelling multiphase visco-elasto-plastic thermo-mechanical problems
- 1 August 2007
- journal article
- Published by Elsevier BV in Physics of the Earth and Planetary Interiors
- Vol. 163 (1-4), 83-105
- https://doi.org/10.1016/j.pepi.2007.04.015
Abstract
We have extended our previous 2D method [Gerya, T.V., Yuen, D.A., 2003. Characteristics-based marker-in-cell method with conservative finite-differences schemes for modeling geological flows with strongly variable transport properties. Phys. Earth Planet. Interiors 140, 295–320], which is a combination of conservative finite-differences with marker-in-cell techniques to include the effects of visco-elasto-plastic rheology, self-gravitation and a self-consistently derived evolving curvilinear planetary surface. This code is called I2ELVIS and can solve a new class of computationally challenging problems in geodynamics, such as shear localization with large strains, crustal intrusion emplacement of magmas, bending of realistic visco-elasto-plastic plates and core-formation by vigorous shell tectonics activities related to a global Rayleigh–Taylor instability of a metal layer formed around silicate-rich lower density (primordial) core during planetary accretion. We discuss in detail the computational strategy required the rheological constraints to be satisfied at each time step and spatial location. We show analytical benchmarks and examples drawn from comparing between numerical and analogue experiments in structural geology, subducting slab bending with a visco-elasto-plastic rheology and equilibrium spherical configurations from self-gravitation. We have also tested possibilities of future applications by addressing 3D geometries based on multigrid method and including inertial effects in the momentum equation with tracers in order to simulate meteoritic impact events and eventually earthquake instabilities.Keywords
This publication has 33 references indexed in Scilit:
- Intrusion of ultramafic magmatic bodies into the continental crust: Numerical simulationPhysics of the Earth and Planetary Interiors, 2007
- Seismic implications of mantle wedge plumesPhysics of the Earth and Planetary Interiors, 2006
- Large‐scale rigid‐body rotation in the mantle wedge and its implications for seismic tomographyGeochemistry, Geophysics, Geosystems, 2006
- Relationship between bend‐faulting at trenches and intermediate‐depth seismicityGeochemistry, Geophysics, Geosystems, 2005
- Visualization of multi-scale dynamics of hydrous cold plumes at subduction zonesVisual Geosciences, 2004
- Catastrophic initiation of subduction following forced convergence across fracture zonesEarth and Planetary Science Letters, 2003
- Thermal evolution of a growing MarsJournal of Geophysical Research, 2002
- Thermal-mechanical effects of low-temperature plasticity (the Peierls mechanism) on the deformation of a viscoelastic shear zoneEarth and Planetary Science Letters, 1999
- Self‐organized criticality of plastic shear bands in rocksGeophysical Research Letters, 1994
- Polydiapirs: multiwavelength gravity structuresJournal of Structural Geology, 1992