Interface error analysis for numerical wave propagation

The numerical error associated with finite-difference simulation of wave propagation in discontinuous media consists of two components. The first component is a higher-order error that leads to grid dispersion; it can be controlled by higher-order methods. The second component results from misalignment between numerical grids and material interfaces. We provide an explicit estimate of the interface misalignment error for the second order in time and space staggered finite-difference scheme applied to the acoustic wave equation. Our analysis, confirmed by numerical experiments, demonstrates that the interface error results in a first-order time shift proportional to the distance between the interface and computational grids. A 2D experiment shows that the interface error cannot be suppressed by higher-order methods and indicates that our 1D analysis gives a good prediction about the behavior of the numerical solution in higher dimensions.

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

[2]  R. LeVeque,et al.  The immersed interface method for acoustic wave equations with discontinuous coefficients , 1997 .

[3]  Qing Huo Liu,et al.  Higher-Order Numerical Methods for Transient Wave Equations , 2003 .

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

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

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

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

[8]  M. B. Dobrin.  Introduction to Geophysical Prospecting , 1976 .

[9]  J. Douglas Faires,et al.  Numerical Analysis , 1981 .

[10]  David L. Brown.  A note on the numerical solution of the wave equation with piecewise smooth coefficients , 1984 .

[11]  A. T. Walden,et al.  The nature of the non-Gaussianity of primary reflection coefficients and its significance for deconvolution , 1986 .

[12]  Gary Cohen.  Higher-Order Numerical Methods for Transient Wave Equations , 2001 .

[13]  Olivier Coussy,et al.  Acoustics of Porous Media , 1988 .

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