Stable and accurate wave-propagation in discontinuous media

A time stable discretization is derived for the second-order wave equation with discontinuous coefficients. The discontinuity corresponds to inhomogeneity in the underlying medium and is treated by splitting the domain. Each (homogeneous) sub domain is discretized using narrow-diagonal summation by parts operators and, then, patched to its neighbors by using a penalty method, leading to fully explicit time integration. This discretization yields a time stable and efficient scheme. The analysis is verified by numerical simulations in one-dimension using high-order finite difference discretizations, and in three-dimensions using an unstructured finite volume discretization.

[1]  Patrick Joly,et al.  Construction and Analysis of Fourth-Order Finite Difference Schemes for the Acoustic Wave Equation in Nonhomogeneous Media , 1996 .

[2]  Eli Turkel,et al.  A fourth-order accurate finite-difference scheme for the computation of elastic waves , 1986 .

[3]  Roland Glowinski,et al.  A Domain Decomposition Method for the Acoustic Wave Equation with Discontinuous Coefficients and Grid Change , 1997 .

[4]  H. Kreiss,et al.  Modeling the black hole excision problem , 2004, gr-qc/0412101.

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

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

[7]  Jan Nordström,et al.  Boundary conditions for a divergence free velocity-pressure formulation of the Navier-Stokes equations , 2007, J. Comput. Phys..

[8]  Gioel Calabrese.  Finite differencing second order systems describing black hole spacetimes , 2004, gr-qc/0410062.

[9]  T. Hagstrom.  Radiation boundary conditions for the numerical simulation of waves , 1999, Acta Numerica.

[10]  K. Yee.  Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media , 1966 .

[11]  H. Kreiss,et al.  Stability Theory of Difference Approximations for Mixed Initial Boundary Value Problems. II , 1972 .

[12]  Magnus Svärd,et al.  Stable and accurate schemes for the compressible Navier-Stokes equations , 2008, J. Comput. Phys..

[13]  Jonas Nycander,et al.  Tidal generation of internal waves from a periodic array of steep ridges , 2006, Journal of Fluid Mechanics.

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

[15]  Marcus J. Grote,et al.  Discontinuous Galerkin Finite Element Method for the Wave Equation , 2006, SIAM J. Numer. Anal..

[16]  Bertil Gustafsson,et al.  Time Compact High Order Difference Methods for Wave Propagation , 2004, SIAM J. Sci. Comput..

[17]  D. Schötzau,et al.  Interior penalty discontinuous Galerkin method for Maxwell's equations , 2007 .

[18]  Jan Nordström,et al.  High order finite difference methods for wave propagation in discontinuous media , 2006, J. Comput. Phys..

[19]  John B. Bell,et al.  A modified equation approach to constructing fourth order methods for acoustic wave propagation , 1987 .

[20]  Heinz-Otto Kreiss,et al.  Difference Approximations for the Second Order Wave Equation , 2002, SIAM J. Numer. Anal..

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

[22]  J. Nordström,et al.  Summation by Parts Operators for Finite Difference Approximations of Second-Derivatives with Variable Coefficients , 2004, Journal of Scientific Computing.

[23]  S. Tsynkov.  Numerical solution of problems on unbounded domains. a review , 1998 .

[24]  S. Lele.  Compact finite difference schemes with spectral-like resolution , 1992 .

[25]  Heinz-Otto Kreiss,et al.  An Embedded Boundary Method for the Wave Equation with Discontinuous Coefficients , 2005, SIAM J. Sci. Comput..

[26]  J. Virieux.  P-SV wave propagation in heterogeneous media: Velocity‐stress finite‐difference method , 1986 .

[27]  Jan Nordström,et al.  Boundary Conditions for a Divergence Free Velocity-Pressure Formulation of the Incompressible Navier-Stokes Equations , 2006 .

[28]  Jean Virieux,et al.  SH-wave propagation in heterogeneous media: velocity-stress finite-difference method , 1984 .

[29]  K. R. Kelly,et al.  SYNTHETIC SEISMOGRAMS: A FINITE ‐DIFFERENCE APPROACH , 1976 .

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

[31]  Bertil Gustafsson,et al.  Time Compact Difference Methods for Wave Propagation in Discontinuous Media , 2004, SIAM J. Sci. Comput..

[32]  Heinz-Otto Kreiss,et al.  Difference Approximations of the Neumann Problem for the Second Order Wave Equation , 2004, SIAM J. Numer. Anal..

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

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

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

[36]  Björn Engquist,et al.  High order difference methods for wave propagation in discontinuous media , 2006 .

[37]  HEINZ-OTTO KREISS,et al.  A Second Order Accurate Embedded Boundary Method for the Wave Equation with Dirichlet Data , 2005, SIAM J. Sci. Comput..

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

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

[40]  Parviz Moin,et al.  Towards time-stable and accurate LES on unstructured grids , 2007 .