A unified first order hyperbolic model for nonlinear dynamic rupture processes in diffuse fracture zones

Abstract

Earthquake fault zones are more complex, both geometrically and rheologically, than an idealised infinitely thin plane embedded in linear elastic material. To incorporate nonlinear material behaviour, natural complexities, and multi-physics coupling within and outside of fault zones, here we present a first-order hyperbolic and thermodynamically compatible mathematical model for a continuum in a gravitational field which provides a unified description of nonlinear elasto-plasticity, material damage and of viscous Newtonian flows with phase transition between solid and liquid phases. The fault geometry and secondary cracks are described via a scalar function $\xi \in [0,1]$ that indicates the local level of material damage. The model also permits the representation of arbitrarily complex geometries via a diffuse interface approach based on the solid volume fraction function $\alpha \in [0,1]$. Neither of the two scalar fields $\xi$ and $\alpha$ needs to be mesh-aligned, allowing thus faults and cracks with complex topology and the use of adaptive Cartesian meshes (AMR). The model shares common features with phase-field approaches but substantially extends them. We show a wide range of numerical applications that are relevant for dynamic earthquake rupture in fault zones, including the co-seismic generation of secondary off-fault shear cracks, tensile rock fracture in the Brazilian disc test, as well as a natural convection problem in molten rock-like material.

Supplementary Materials are at the end of the main manuscript

References in corpus (16)

• A Posteriori Subcell Limiting of the Discontinuous Galerkin Finite Element Method for Hyperbolic Conservation Laws
• Dynamics, radiation and overall energy budget of earthquake rupture with coseismic off-fault damage
• High order ADER schemes for a unified first order hyperbolic formulation of continuum mechanics: viscous heat-conducting fluids and elastic solids
• Phase field modeling of crack propagation
• Space-time adaptive ADER discontinuous Galerkin finite element schemes with a posteriori sub-cell finite volume limiting
• Fracture in Mode I using a Conserved Phase-Field Model
• High order ADER schemes for continuum mechanics
• Continuum Mechanics and Thermodynamics in the Hamilton and the Godunov-type Formulations
• ExaHyPE: An Engine for Parallel Dynamically Adaptive Simulations of Wave Problems
• Efficient high order accurate staggered semi-implicit discontinuous Galerkin methods for natural convection problems
• High order ADER schemes and GLM curl cleaning for a first order hyperbolic formulation of compressible flow with surface tension
• Space-time adaptive ADER discontinuous Galerkin schemes for nonlinear hyperelasticity with material failure
• Continuum mechanics with torsion
• A simple diffuse interface approach on adaptive Cartesian grids for the linear elastic wave equations with complex topography
• Modeling wavefields in saturated elastic porous media based on thermodynamically compatible system theory for multiphase mixtures
• A Numerical Scheme for Non-Newtonian Fluids and Plastic Solids under the GPR Model

Cited by in corpus (7)

• Simulation of non-Newtonian viscoplastic flows with a unified first order hyperbolic model and a structure-preserving semi-implicit scheme
• Exact and numerical solutions of the Riemann problem for a conservative model of compressible two-phase flows
• A phase-field model for quasi-dynamic nucleation, growth, and propagation of rate-and-state faults
• Two-phase hyperbolic model for porous media saturated with a viscous fluid and its application to wavefields simulation
• High order well-balanced Arbitrary-Lagrangian-Eulerian ADER discontinuous Galerkin schemes on general polygonal moving meshes
• A unified multi-phase and multi-material formulation for combustion modelling
• Thermodynamically consistent model for poroelastic rocks towards tectonic and volcanic processes and earthquakes