Schwarz method for earthquake source dynamics

Dynamic faulting under slip-dependent friction in a linear elastic domain (in-plane and 3D configurations) is considered. The use of an implicit time-stepping scheme (Newmark method) allows much larger values of the time step than the critical CFL time step, and higher accuracy to handle the non-smoothness of the interface constitutive law (slip weakening friction). The finite element form of the quasi-variational inequality is solved by a Schwarz domain decomposition method, by separating the inner nodes of the domain from the nodes on the fault. In this way, the quasi-variational inequality splits into two subproblems. The first one is a large linear system of equations, and its unknowns are related to the mesh nodes of the first subdomain (i.e. lying inside the domain). The unknowns of the second subproblem are the degrees of freedom of the mesh nodes of the second subdomain (i.e. lying on the domain boundary where the conditions of contact and friction are imposed). This nonlinear subproblem is solved by the same Schwarz algorithm, leading to some local nonlinear subproblems of a very small size. Numerical experiments are performed to illustrate convergence in time and space, instability capturing, energy dissipation and the influence of normal stress variations. We have used the proposed numerical method to compute source dynamics phenomena on complex and realistic 2D fault models (branched fault systems).

[1]  C. Dascalu,et al.  SLIP-WEAKENING FRICTION INSTABILITIES: EIGENVALUE ANALYSIS , 2004 .

[2]  O. Scotti,et al.  Dynamic transfer of rupture across differently oriented segments in a complex 3‐D fault system , 2005 .

[3]  A. Toselli.  Domain Decomposition Methods , 2004 .

[4]  L. Badea,et al.  Domain decomposition method for dynamic faulting under slip-dependent friction , 2004 .

[5]  Harsha S. Bhat,et al.  DYNAMIC SLIP TRANSFER FROM THE DENALI TO TOTSCHUNDA FAULTS, ALASKA: TESTING THEORY FOR FAULT BRANCHING , 2004 .

[6]  Brad Thomas Aagaard.  Finite-element simulations of earthquakes , 2000 .

[7]  Jean Virieux,et al.  Dynamic rupture simulation of non-planar faults with a finite-difference approach , 2004 .

[8]  Rolf Krause,et al.  Solving Frictional Contact Problems with Multigrid Efficiency , 2007 .

[9]  Jaroslav Haslinger,et al.  Numerical methods for unilateral problems in solid mechanics , 1996 .

[10]  Hideo Aochi,et al.  3D finite-difference dynamic-rupture modeling along nonplanar faults , 2007 .

[11]  T. Laursen.  Computational Contact and Impact Mechanics , 2003 .

[12]  Andrea Toselli,et al.  Domain decomposition methods : algorithms and theory , 2005 .

[13]  B. Duan,et al.  Nonuniform prestress from prior earthquakes and the effect on dynamics of branched fault systems , 2007 .

[14]  É. Oudet,et al.  Critical friction for wedged configurations , 2007 .

[15]  W. Ellsworth,et al.  Seismic Evidence for an Earthquake Nucleation Phase , 1995, Science.

[16]  R. Kornhuber.  Adaptive monotone multigrid methods for nonlinear variational problems , 1997 .

[17]  Jean-Pierre Vilotte,et al.  Influence of the rupture initiation on the intersonic transition: Crack‐like versus pulse‐like modes , 2006 .

[18]  Robert Burridge,et al.  Admissible Speeds for Plane-Strain Self-Similar Shear Cracks with Friction but Lacking Cohesion , 1973 .

[19]  J. Oden,et al.  Contact problems in elasticity , 1988 .

[20]  Y. Iio.  Slow initial phase of the P‐wave velocity pulse generated by microearthquakes , 1992 .

[21]  Steven M. Day,et al.  The 1999 Hector Mine Earthquake: The Dynamics of a Branched Fault System , 2001 .

[22]  Z. Zhong.  Finite Element Procedures for Contact-Impact Problems , 1993 .

[23]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[24]  J. C. Simo,et al.  A continuum-based finite element formulation for the implicit solution of multibody, large deformation-frictional contact problems , 1993 .

[25]  Barbara I. Wohlmuth,et al.  A Mortar Finite Element Method Using Dual Spaces for the Lagrange Multiplier , 2000, SIAM J. Numer. Anal..

[26]  Mitiyasu Ohnaka,et al.  Constitutive relations between dynamic physical parameters near a tip of the propagating slip zone during stick-slip shear failure , 1987 .

[27]  D. Dreger,et al.  Inverse Kinematic and Forward Dynamic Models of the 2002 Denali Fault Earthquake, Alaska , 2004 .

[28]  M. Campillo,et al.  Influence of the shape of the friction law and fault finiteness on the duration of initiation , 1999 .

[29]  David D. Oglesby,et al.  Earthquake dynamics on dip-slip faults , 1999 .

[30]  D. Oglesby,et al.  Erratum: ``Dynamics of dip-slip faulting: Explorations in two dimensions'' , 2000 .

[31]  Ralph J. Archuleta,et al.  Three-dimensional numerical simulations of dynamic faulting in a half-space , 1978, Bulletin of the Seismological Society of America.

[32]  T. Yamashita,et al.  The paradox of smooth and abrupt bends in two-dimensional in-plane shear-crack mechanics , 1996 .

[33]  D. J. Andrews,et al.  Rupture velocity of plane strain shear cracks , 1976 .

[34]  Barry F. Smith,et al.  Domain Decomposition: Parallel Multilevel Methods for Elliptic Partial Differential Equations , 1996 .

[35]  Jean Virieux,et al.  Dynamic non-planar crack rupture by a finite volume method , 2006 .

[36]  Sonia Fliss,et al.  Fault branching and rupture directivity , 2002 .

[37]  Michel Campillo,et al.  Initiation of in-plane shear instability under slip-dependent friction , 1997, Bulletin of the Seismological Society of America.

[38]  R. Glowinski.  Lectures on Numerical Methods for Non-Linear Variational Problems , 1981 .

[39]  H. Aochi,et al.  Spontaneous Rupture Propagation on a Non-planar Fault in 3-D Elastic Medium , 2000 .

[40]  Michel Campillo,et al.  Fault finiteness and initiation of dynamic shear instability , 2000 .

[41]  Peter Wriggers,et al.  Computational Contact Mechanics , 2002 .

[42]  John R. Rice,et al.  Universal nucleation length for slip-weakening rupture instability under nonuniform fault loading , 2003 .

[43]  T. Fung.  Unconditionally stable higher-order Newmark methods by sub-stepping procedure , 1997 .

[44]  S. Day,et al.  The effect of fault geometry on the 1999 Chi‐Chi (Taiwan) Earthquake , 2001 .

[45]  Nobuki Kame,et al.  Effects of prestress state and rupture velocity on dynamic fault branching , 2002 .

[46]  Barry Smith,et al.  Domain Decomposition Methods for Partial Differential Equations , 1997 .

[47]  M. Campillo,et al.  Initiation of shear instability in three-dimensional elastodynamics , 2002 .