The numerical modelling of plate subduction requires solving a coupled thermo-mechanical highly-nonlinear transient problem. The mechanical description of the phenomenon results in a multiphase quasi-static Stokes flow, where the inertia terms are neglected. The transient thermal problem is dominated by the advection term. Here, the representation and evolution of the different phases are described using level sets. The phase tracking is carried out transporting the level set along with the material, using a pure advective model. The gradient discontinuities induced by the viscosity jump across the interface are resolved numerically by enriching the solution using a partition of unity method in a eXtended Finite Element Method (X-FEM) context. These numerical tools are used to simulate plate subduction with different parameters and to derive useful correlations between relevant geophysical factors.

The PDF file did not load properly or your web browser does not support viewing PDF files. Download directly to your device: Download PDF document
Back to Top

Document information

Published on 28/10/19

DOI: 10.1016/j.cma.2007.04.006
Licence: CC BY-NC-SA license

Document Score


Views 0
Recommendations 0

Share this document

claim authorship

Are you one of the authors of this document?