Comparison of High-Accuracy Finite-Difference Methods for Linear Wave Propagation

This paper analyzes a number of high-order and optimized finite-difference methods for numerically simulating the propagation and scattering of linear waves, such as electromagnetic, acoustic, and elastic waves. The spatial operators analyzed include compact schemes, noncompact schemes, schemes on staggered grids, and schemes which are optimized to produce specific characteristics. The time-marching methods include Runge--Kutta methods, Adams--Bashforth methods, and the leapfrog method. In addition, the following fully-discrete finite-difference methods are studied: a one-step implicit scheme with a three-point spatial stencil, a one-step explicit scheme with a five-point spatial stencil, and a two-step explicit scheme with a five-point spatial stencil. For each method, the number of grid points per wavelength required for accurate simulation of wave propagation over large distances is presented. The results provide a clear understanding of the relative merits of the methods compared, especially the trade-offs associated with the use of optimized methods. A numerical example is given which shows that the benefits of an optimized scheme can be small if the waveform has broad spectral content.

[1]  James Lighthill.  The Final Panel Discussion , 1993 .

[2]  Robert Vichnevetsky,et al.  A frequency analysis of finite difference and finite element methods for initial value problems , 1975 .

[3]  K. Marfurt.  Accuracy of finite-difference and finite-element modeling of the scalar and elastic wave equations , 1984 .

[4]  David W. Zingg,et al.  Runge-Kutta methods for linear ordinary differential equations , 1999 .

[5]  Richard W. Ziolkowski,et al.  A Three-Dimensional Modified Finite Volume Technique for Maxwell's Equations , 1990 .

[6]  Patrick Joly,et al.  Fourth order schemes for the heterogeneous acoustics equation , 1990 .

[7]  David W. Zingg,et al.  An optimized finite-difference scheme for wave propagation problems , 1993 .

[8]  B. Gustafsson.  The convergence rate for difference approximations to mixed initial boundary value problems , 1975 .

[9]  H. Lomax,et al.  Finite-Difference Schemes on Regular Triangular Grids , 1993 .

[10]  Jeffrey L. Young,et al.  Toward the construction of a fourth-order difference scheme for transient EM wave simulation: staggered grid approach , 1997 .

[11]  P. Sguazzero,et al.  Dispersion-bounded numerical integration of the elastodynamic equations with cost-effective staggered schemes , 1990 .

[12]  S. Gedney.  An anisotropic perfectly matched layer-absorbing medium for the truncation of FDTD lattices , 1996 .

[13]  Finite-surface method for the Maxwell equations with corner singularities , 1994 .

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

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

[16]  Shlomo Ta'asan,et al.  Finite difference schemes for long-time integration , 1994 .

[17]  Allen Taflove,et al.  Finite-difference time-domain modeling of curved surfaces (EM scattering) , 1992 .

[18]  Alain Sei.  A Family of Numerical Schemes for the Computation of Elastic Waves , 1995, SIAM J. Sci. Comput..

[19]  Finite-surface method for the Maxwell equations in generalized coordinates , 1993 .

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

[21]  V. Shankar,et al.  A Time-Domain, Finite-Volume Treatment for the Maxwell Equations , 1990 .

[22]  J. S. Shang.  A fractional-step method for solving 3D, time-domain Maxwell equations , 1995 .

[23]  D. Gottlieb,et al.  Regular Article: Well-posed Perfectly Matched Layers for Advective Acoustics , 1999 .

[24]  Yen Liu,et al.  Fourier analysis of numerical algorithms for the Maxwell equations , 1993 .

[25]  R. M. Alford,et al.  ACCURACY OF FINITE‐DIFFERENCE MODELING OF THE ACOUSTIC WAVE EQUATION , 1974 .

[26]  J. Bowles,et al.  Fourier Analysis of Numerical Approximations of Hyperbolic Equations , 1987 .

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

[28]  J. Lambert.  Computational Methods in Ordinary Differential Equations , 1973 .

[29]  M. A. Dablain,et al.  The application of high-order differencing to the scalar wave equation , 1986 .

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

[31]  Andreas C. Cangellaris,et al.  A Reflectionless Sponge Layer Absorbing Boundary Condition for the Solution of Maxwell's Equations with High-Order Staggered Finite Difference Schemes , 1998 .

[32]  Richard W. Ziolkowski,et al.  Time-derivative Lorentz material model-based absorbing boundary condition , 1997 .

[33]  Rosemary A. Renaut,et al.  COMPUTING AERODYNAMICALLY GENERATED NOISE , 1997 .

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

[35]  Philip Roe.  Linear Bicharacteristic Schemes Without Dissipation , 1998, SIAM J. Sci. Comput..

[36]  Olav Holberg,et al.  COMPUTATIONAL ASPECTS OF THE CHOICE OF OPERATOR AND SAMPLING INTERVAL FOR NUMERICAL DIFFERENTIATION IN LARGE-SCALE SIMULATION OF WAVE PHENOMENA* , 1987 .

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

[38]  M.Y. Hussaini,et al.  Low-Dissipation and Low-Dispersion Runge-Kutta Schemes for Computational Acoustics , 1994 .

[39]  M. Fusco,et al.  FDTD algorithm in curvilinear coordinates (EM scattering) , 1990 .

[40]  Christopher K. W. Tam,et al.  Computational aeroacoustics - Issues and methods , 1995 .

[41]  Andreas C. Cangellaris,et al.  Point-matched time domain finite element methods for electromagnetic radiation and scattering , 1987 .

[42]  J. Kim,et al.  Optimized Compact Finite Difference Schemes with Maximum Resolution , 1996 .

[43]  Sanford S. Davis,et al.  Low‐dispersion finite difference methods for acoustic waves in a pipe , 1991 .

[44]  David W. Zingg,et al.  Numerical Solution of the Time-Domain Maxwell Equations Using High-Accuracy Finite-Difference Methods , 2000, SIAM J. Sci. Comput..

[45]  David W. Zingg,et al.  High-Accuracy Finite-Difference Schemes for Linear Wave Propagation , 1996, SIAM J. Sci. Comput..

[46]  D. Gottlieb,et al.  Stable and accurate boundary treatments for compact, high-order finite-difference schemes , 1993 .

[47]  Peter G. Petropoulos,et al.  Phase error control for FD-TD methods of second and fourth order accuracy , 1994 .

[48]  C. Tam,et al.  Dispersion-relation-preserving finite difference schemes for computational acoustics , 1993 .

[49]  Allen Taflove,et al.  Computational Electrodynamics the Finite-Difference Time-Domain Method , 1995 .

[50]  David P. Lockard,et al.  High accuracy algorithms for computational aeroacoustics , 1994 .

[51]  E. Turkel,et al.  Fourth order method for Maxwell equations on a staggered mesh , 1997, IEEE Antennas and Propagation Society International Symposium 1997. Digest.

[52]  A. Taflove.  The Finite-Difference Time-Domain Method , 1995 .

[53]  Datta V. Gaitonde,et al.  On high resolution schemes for time-dependent Maxwell equations , 1996 .