High order finite difference methods for wave propagation in discontinuous media

High order finite difference approximations are derived for the second order wave equation with discontinuous coefficients, on rectangular geometries. The discontinuity is treated by splitting the domain at the discontinuities in a multi block fashion. Each sub-domain is discretized with compact second derivative summation by parts operators and the blocks are patched together to a global domain using the projection method. This guarantees a conservative, strictly stable and high order accurate scheme. The analysis is verified by numerical simulations in one and two spatial dimensions.

[1]  Jing Gong,et al.  A stable hybrid method for hyperbolic problems , 2006, J. Comput. Phys..

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

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

[4]  Tobin A. Driscoll,et al.  Block Pseudospectral Methods for Maxwell's Equations II: Two-Dimensional, Discontinuous-Coefficient Case , 1999, SIAM J. Sci. Comput..

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

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

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

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

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

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

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

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

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

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

[15]  Jan Nordström,et al.  A Stable and Efficient Hybrid Method for Aeroacoustic Sound Generation and Propagation , 2005 .

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

[17]  D. Gottlieb,et al.  The Stability of Numerical Boundary Treatments for Compact High-Order Finite-Difference Schemes , 1993 .

[18]  Jan Nordström,et al.  High-order finite difference methods, multidimensional linear problems, and curvilinear coordinates , 2001 .

[19]  Oscar Reula,et al.  Summation by parts and dissipation for domains with excised regions , 2003, gr-qc/0308007.

[20]  J. Hesthaven,et al.  Convergent Cartesian Grid Methods for Maxwell's Equations in Complex Geometries , 2001 .

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

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

[23]  Bo Strand,et al.  High-Order Difference Approximations for Hyperbolic Initial Boundary Value Problems , 1996 .

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

[25]  S. Abarbanel,et al.  Asymptotically Stable Fourth-Order Accurate Schemes for the Diffusion Equation on Complex Shapes , 1997 .

[26]  H. Kreiss.  Stability theory for difference approximations of mixed initial boundary value problems. I , 1968 .

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

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

[29]  Chi Hou Chan,et al.  An explicit fourth-order orthogonal curvilinear staggered-grid FDTD method for Maxwell's equations , 2002 .

[30]  David Neilsen,et al.  The discrete energy method in numerical relativity: Towards long-term stability , 2004, gr-qc/0406116.