High-Accuracy Finite-Difference Schemes for Linear Wave Propagation

Two high-accuracy finite-difference schemes for simulating long-range linear wave propagation are presented: a maximum-order scheme and an optimized scheme. The schemes combine a seven-point spatial operator and an explicit six-stage low-storage time-march method of Runge--Kutta type. The maximum-order scheme can accurately simulate the propagation of waves over distances greater than five hundred wavelengths with a grid resolution of less than twenty points per wavelength. The optimized scheme is found by minimizing the maximum phase and amplitude errors for waves which are resolved with at least ten points per wavelength, based on Fourier error analysis. It is intended for simulations in which waves travel under three hundred wavelengths. For such cases, good accuracy is obtained with roughly ten points per wavelength.

[1]  A. Majda,et al.  Absorbing boundary conditions for the numerical simulation of waves , 1977 .

[2]  V. Shankar,et al.  A three-dimensional Maxwell's equation solver for computation of scattering from layered media , 1989 .

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

[4]  M. Fusco,et al.  A three-dimensional FDTD algorithm in curvilinear coordinates (EM scattering) , 1991 .

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

[6]  Ray Hixon,et al.  Evaluation of boundary conditions for computational aeroacoustics , 1995 .

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

[8]  Robert L. Higdon,et al.  Numerical absorbing boundary conditions for the wave equation , 1987 .

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

[10]  K. Yee,et al.  A subgridding method for the time-domain finite-difference method to solve Maxwell's equations , 1991 .

[11]  D. Gaitonde,et al.  Scattering simulations of computational electromagnetics , 1996 .

[12]  K. Thompson.  Time-dependent boundary conditions for hyperbolic systems, II , 1990 .

[13]  D. Katz,et al.  Validation and extension to three dimensions of the Berenger PML absorbing boundary condition for FD-TD meshes , 1994, IEEE Microwave and Guided Wave Letters.

[14]  Roger F. Harrington,et al.  Field computation by moment methods , 1968 .