Stable, high order accurate adaptive schemes for long time, highly intermittent geophysics problems

Many geophysical phenomena are characterized by properties that evolve over a wide range of scales which introduce difficulties when attempting to model these features in one computational method. We have developed a high-order finite difference method for the elastic wave equation that is able to efficiently handle varying temporal and spatial scales in a single, stand-alone framework. We apply this method to earthquake cycle models characterized by extremely long interseismic periods interspersed with abrupt, short periods of dynamic rupture. Through the use of summation-by-parts operators and weak enforcement of boundary conditions we derive a provably stable discretization. Time stepping is achieved through the implicit θθ-method which allows us to take large time steps during the intermittent period between earthquakes and adapts appropriately to fully resolve rupture.

[1]  D. Gottlieb,et al.  A Stable and Conservative Interface Treatment of Arbitrary Spatial Accuracy , 1999 .

[2]  B. Strand.  Summation by parts for finite difference approximations for d/dx , 1994 .

[3]  P. Olsson.  Summation by parts, projections, and stability. II , 1995 .

[4]  E. Dunham,et al.  Earthquake Ruptures with Strongly Rate-Weakening Friction and Off-Fault Plasticity , Part 2 : Nonplanar Faults by , 2011 .

[5]  S. Day,et al.  Rupture dynamics and ground motion from 3‐D rough‐fault simulations , 2013 .

[6]  J. Dieterich.  Modeling of rock friction: 1. Experimental results and constitutive equations , 1979 .

[7]  N. Anders Petersson,et al.  Stable Grid Refinement and Singular Source Discretization for Seismic Wave Simulations , 2009 .

[8]  B. Shibazaki,et al.  Spontaneous processes for nucleation, dynamic propagation, and stop of earthquake rupture , 1992 .

[9]  J. Nordström,et al.  A flexible boundary procedure for hyperbolic problems : Multiple penalty terms applied in a domain , 2014 .

[10]  Jeremy E. Kozdon,et al.  Rupture to the Trench: Dynamic Rupture Simulations of the 11 March 2011 Tohoku Earthquake , 2013 .

[11]  Ken Mattsson,et al.  Boundary Procedures for Summation-by-Parts Operators , 2003, J. Sci. Comput..

[12]  H. Kreiss,et al.  Time-Dependent Problems and Difference Methods , 1996 .

[13]  S. Roux,et al.  Velocity weakening friction: A renormalization approach , 1996 .

[14]  Frank Chorlton.  Summation by Parts , 1998 .

[15]  E. Dunham,et al.  Earthquake Ruptures with Strongly Rate-Weakening Friction and Off-Fault Plasticity, Part 1: Planar Faults , 2011 .

[16]  Olsson,et al.  SUMMATION BY PARTS, PROJECTIONS, AND STABILITY. I , 2010 .

[17]  Jing Gong,et al.  Interface procedures for finite difference approximations of the advection-diffusion equation , 2011, J. Comput. Appl. Math..

[18]  Jan Nordström,et al.  Conjugate heat transfer for the unsteady compressible Navier–Stokes equations using a multi-block coupling , 2013 .

[19]  Roger B. Nelsen,et al.  Summation by Parts , 1992 .

[20]  A computational study of vortex–airfoil interaction using high-order finite difference methods , 2010 .

[21]  Magnus Svärd,et al.  Well-Posed Boundary Conditions for the Navier-Stokes Equations , 2005, SIAM J. Numer. Anal..

[22]  C. Marone.  LABORATORY-DERIVED FRICTION LAWS AND THEIR APPLICATION TO SEISMIC FAULTING , 1998 .

[23]  Timothy G. Trucano,et al.  Verification and validation. , 2005 .

[24]  P. Okubo.  Dynamic rupture modeling with laboratory‐derived constitutive relations , 1989 .

[25]  Jeremy E. Kozdon,et al.  Interaction of Waves with Frictional Interfaces Using Summation-by-Parts Difference Operators: Weak Enforcement of Nonlinear Boundary Conditions , 2012, J. Sci. Comput..

[26]  Magnus Svärd,et al.  On the order of accuracy for difference approximations of initial-boundary value problems , 2006, J. Comput. Phys..

[27]  Daniel J. Bodony,et al.  Accuracy of the Simultaneous-Approximation-Term Boundary Condition for Time-Dependent Problems , 2010, J. Sci. Comput..

[28]  G. Bergeles,et al.  Notes on Numerical Fluid Mechanics and Multidisciplinary Design , 2012 .

[29]  H. Kreiss,et al.  Finite Element and Finite Difference Methods for Hyperbolic Partial Differential Equations , 1974 .

[30]  E. Dunham,et al.  An efficient numerical method for earthquake cycles in heterogeneous media: Alternating subbasin and surface‐rupturing events on faults crossing a sedimentary basin , 2014 .

[31]  Jing Gong,et al.  A stable and conservative high order multi-block method for the compressible Navier-Stokes equations , 2009, J. Comput. Phys..

[32]  J. Dieterich.  Time-dependent friction and the mechanics of stick-slip , 1978 .

[33]  Jan Nordström,et al.  High Order Finite Difference Approximations of Electromagnetic Wave Propagation Close to Material Discontinuities , 2003, J. Sci. Comput..

[34]  Jan Nordström,et al.  High-order accurate difference schemes for the Hodgkin-Huxley equations , 2012, J. Comput. Phys..

[35]  Magnus Svärd,et al.  Review of summation-by-parts schemes for initial-boundary-value problems , 2013, J. Comput. Phys..

[36]  Jan Nordström,et al.  Revisiting and Extending Interface Penalties for Multi-domain Summation-by-Parts Operators , 2010, J. Sci. Comput..

[37]  Patrick J. Roache,et al.  Verification and Validation in Computational Science and Engineering , 1998 .

[38]  Jan Nordström,et al.  Boundary and Interface Conditions for High-Order Finite-Difference Methods Applied to the Euler and Navier-Stokes Equations , 1999 .

[39]  J. Dieterich,et al.  Direct observation of frictional contacts: New insights for state-dependent properties , 1994 .

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

[41]  David W. Zingg,et al.  Comparison of High-Accuracy Finite-Difference Methods for Linear Wave Propagation , 2000, SIAM J. Sci. Comput..

[42]  D. Gottlieb,et al.  Time-stable boundary conditions for finite-difference schemes solving hyperbolic systems: methodology and application to high-order compact schemes , 1994 .

[43]  J. Nordström.  Linear and Nonlinear Boundary Conditions for Wave Propagation Problems , 2013 .

[44]  H. Kreiss,et al.  Comparison of accurate methods for the integration of hyperbolic equations , 1972 .

[45]  Uri M. Ascher,et al.  Computer methods for ordinary differential equations and differential-algebraic equations , 1998 .

[46]  N. SIAMJ.,et al.  WELL-POSED BOUNDARY CONDITIONS FOR THE NAVIER – STOKES EQUATIONS , 2005 .