A three‐dimensional hybrid finite element — spectral boundary integral method for modeling earthquakes in complex unbounded domains

We present a 3D hybrid method which combines the Finite Element Method (FEM) and the Spectral Boundary Integral method (SBIM) to model nonlinear problems in unbounded domains. The flexibility of FEM is used to model the complex, heterogeneous, and nonlinear part -- such as the dynamic rupture along a fault with near fault plasticity -- and the high accuracy and computational efficiency of SBIM is used to simulate the exterior half spaces perfectly truncating all incident waves. The exact truncation allows us to greatly reduce the domain of spatial discretization compared to a traditional FEM approach, leading to considerable savings in computational cost and memory requirements. The coupling of FEM and SBIM is achieved by the exchange of traction and displacement boundary conditions at the computationally defined boundary. The method is suited to implementation on massively parallel computers. We validate the developed method by means of a benchmark problem. Three more complex examples with a low velocity fault zone, low velocity off-fault inclusion, and interaction of multiple faults, respectively, demonstrate the capability of the hybrid scheme in solving problems of very large sizes. Finally, we discuss potential applications of the hybrid method for problems in geophysics and engineering.

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

[2]  Yijun Liu,et al.  Hypersingular boundary integral equations for radiation and scattering of elastic waves in three dimensions , 1993 .

[3]  J. Virieux,et al.  Dynamic faulting studied by a finite difference method : Bull seismol soc am, V72, N2, April 1982, P345–369 , 1982 .

[4]  J. Rice,et al.  Earthquake ruptures with thermal weakening and the operation of major faults at low overall stress levels , 2009 .

[5]  D. Kammer,et al.  Off‐fault heterogeneities promote supershear transition of dynamic mode II cracks , 2017 .

[6]  Luis A. Dalguer,et al.  Staggered-grid split-node method for spontaneous rupture simulation , 2007 .

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

[8]  J. Rice,et al.  Elastodynamic analysis for slow tectonic loading with spontaneous rupture episodes on faults with rate‐ and state‐dependent friction , 2000 .

[9]  Raul Madariaga,et al.  Dynamic faulting under rate-dependent friction , 1994 .

[10]  S. Day.  Three-dimensional simulation of spontaneous rupture: The effect of nonuniform prestress , 1982, Bulletin of the Seismological Society of America.

[11]  D. Boore,et al.  A two-dimensional moving dislocation model for a strike-slip fault , 1971 .

[12]  Jean-Paul Ampuero.  Étude physique et numérique de la nucléation des séismes , 2002 .

[13]  Shuo Ma,et al.  Radiated seismic energy based on dynamic rupture models of faulting , 2006 .

[14]  J. Ampuero,et al.  Three‐dimensional dynamic rupture simulation with a high‐order discontinuous Galerkin method on unstructured tetrahedral meshes , 2012 .

[15]  J. Ampuero,et al.  Spectral element modeling of spontaneous earthquake rupture on rate and state faults: Effect of velocity‐strengthening friction at shallow depths , 2008 .

[16]  Jean-Paul Ampuero,et al.  Pulse-like ruptures induced by low-velocity fault zones , 2011 .

[17]  Jeremy E. Kozdon,et al.  Simulation of Dynamic Earthquake Ruptures in Complex Geometries Using High-Order Finite Difference Methods , 2013, J. Sci. Comput..

[18]  R. Courant,et al.  Über die partiellen Differenzengleichungen der mathematischen Physik , 1928 .

[19]  A. Elbanna,et al.  A new hybrid numerical scheme for modelling elastodynamics in unbounded media with near-source heterogeneities , 2017, 1702.04802.

[20]  L. B. Freund,et al.  The mechanics of dynamic shear crack propagation , 1979 .

[21]  Numerical analysis of dynamic debonding under 2D in-plane and 3D loading , 1998 .

[22]  Francisco J. Sánchez-Sesma,et al.  A 3D hp‐adaptive discontinuous Galerkin method for modeling earthquake dynamics , 2012 .

[23]  R. Taylor,et al.  Lagrange constraints for transient finite element surface contact , 1991 .

[24]  B. V. Kostrov,et al.  An investigation of the complexity of the earthquake source time function using dynamic faulting models , 1988 .

[25]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

[26]  James H. Dieterich,et al.  Modeling of rock friction: 2. Simulation of preseismic slip , 1979 .

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

[28]  S. Day,et al.  Comparison of finite difference and boundary integral solutions to three‐dimensional spontaneous rupture , 2005 .

[29]  J. Fineberg,et al.  The equation of motion for supershear frictional rupture fronts , 2018, Science Advances.

[30]  J. Molinari,et al.  The existence of a critical length scale in regularised friction , 2014, 1402.3043.

[31]  Xiao Ma,et al.  Effect of off-fault low-velocity elastic inclusions on supershear rupture dynamics , 2015 .

[32]  D. Andrews.  Dynamic growth of mixed-mode shear cracks , 1994, Bulletin of the Seismological Society of America.

[33]  A. N. Stokes,et al.  An Improved Method for Numerical Inversion of Laplace Transforms , 1982 .

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

[35]  M. Dumbser,et al.  An arbitrary high-order discontinuous Galerkin method for elastic waves on unstructured meshes - I. The two-dimensional isotropic case with external source terms , 2006 .

[36]  The Seismic Radiation from a 3‐D Dynamic Model of a Complex Rupture Process. Part I: Confined Ruptures , 2013 .

[37]  Yoshiaki Ida,et al.  Cohesive force across the tip of a longitudinal‐shear crack and Griffith's specific surface energy , 1972 .

[38]  Martin Käser,et al.  Dynamic rupture modeling on unstructured meshes using a discontinuous Galerkin method , 2009 .

[39]  John R. Rice,et al.  Fault rupture between dissimilar materials: Ill-posedness, regularization, and slip-pulse response , 2000 .

[40]  Ronaldo I. Borja,et al.  An extended finite element framework for slow‐rate frictional faulting with bulk plasticity and variable friction , 2009 .

[41]  A. Ruina.  Slip instability and state variable friction laws , 1983 .

[42]  Shamita Das,et al.  A numerical method for determination of source time functions for general three‐dimensional rupture propagation , 1980 .

[43]  Steven M. Day,et al.  DYNAMIC RUPTURE IN A LAYERED MEDIUM: THE 1966 PARKFIELD EARTHQUAKE , 1980 .

[44]  M. Bonnet,et al.  Regularization of the Displacement and Traction BIE for 3D Elastodynamics Using Indirect Methods , 1993 .

[45]  Jean-François Molinari,et al.  Implementation of a parallel finite-element library , 2015 .

[46]  J. Lysmer,et al.  Finite Dynamic Model for Infinite Media , 1969 .

[47]  Marc Bonnet,et al.  Modelling of dynamical crack propagation using time-domain boundary integral equations , 1992 .

[48]  Keiiti Aki,et al.  A numerical study of two-dimensional spontaneous rupture propagation , 1977 .

[49]  D. J. Andrews,et al.  Bulletin of the Seismological Society of America , 1985 .

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

[51]  S. Day,et al.  Bimaterial effects in an earthquake cycle model using rate‐and‐state friction , 2016 .

[52]  Xiao Ma,et al.  A hybrid finite element‐spectral boundary integral approach: Applications to dynamic rupture modeling in unbounded domains , 2018, International Journal for Numerical and Analytical Methods in Geomechanics.

[53]  D. Helmberger,et al.  Earthquake ruptures modulated by waves in damaged fault zones , 2014 .

[54]  Wing Kam Liu,et al.  Nonlinear Finite Elements for Continua and Structures , 2000 .

[55]  J. Rice,et al.  A spectral method for three-dimensional elastodynamic fracture problems , 1995 .

[56]  D. Helmberger,et al.  The potential for supershear earthquakes in damaged fault zones – theory and observations , 2016 .

[57]  Xiao Ma,et al.  A Novel Hybrid Finite Element‐Spectral Boundary Integral Scheme for Modeling Earthquake Cycles: Application to Rate and State Faults With Low‐Velocity Zones , 2019, Journal of Geophysical Research: Solid Earth.

[58]  Prasanta K. Banerjee,et al.  Two-dimensional transient wave-propagation problems by time-domain BEM , 1990 .

[59]  Gabriele Albertini,et al.  UGUCA: A spectral-boundary-integral method for modeling fracture and friction , 2021, SoftwareX.

[60]  D. Komatitsch,et al.  Introduction to the spectral element method for three-dimensional seismic wave propagation , 1999 .

[61]  S. Das,et al.  On the Numerical Boundary Integral Equation Method for Three-Dimensional Dynamic Shear Crack Problems , 1987 .

[62]  Xiao Ma,et al.  Dynamic rupture propagation on fault planes with explicit representation of short branches , 2019, Earth and Planetary Science Letters.