Dispersion-relation-preserving finit difference schemes for computational acoustics

• Published in: Journal of computational physics Vol. 107:2
• Main Authors: Tam, C.K.W., Webb, J.C.
• Format: Journal Article
• Language: English
• Published: United States 01.08.1993
• Subjects: 990200 > Mathematics & Computers; ACOUSTICS; ASYMPTOTIC SOLUTIONS; CALCULATION METHODS; CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; COMPUTERIZED SIMULATION; FINITE DIFFERENCE METHOD; GENERAL AND MISCELLANEOUS//MATHEMATICS, COMPUTING, AND INFORMATION SCIENCE; ITERATIVE METHODS; MESH GENERATION; NUMERICAL SOLUTION; SIMULATION 661300 > Other Aspects of Physical Science > (1992-); STABILITY
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1006/jcph.1993.1142

Acoustics problems are governed by the linearized Euler equations. According to wave propagation theory, the number of wave modes and their wave propagation characteristics are all encoded in the dispersion relation of the governing equations. Thus one is assured that the numerical solutions of high order finite difference scheme will have the same number of wave modes (namely, the acoustic, vorticity, and entropy waves), the same wave propagation characteristics (namely, nondispersive, nondissipative, and isotropic) and the same wave speeds as those of the solutions of the Euler equations if both systems of equations have the same dispersion relations. Finite difference schemes which have the same dispersion relations as the original partial differential equations are referred to as dispersion-relation-preserving (DRP) schemes. A way to construct time marching DRP schemes by optimizing the finite difference approximations of the space and time derivatives in the wave number and frequency space is proposed. The stability of these schemes is analyzed and a sufficient condition for numerical stability is established. A set of radiation and outflow boundary conditions compatible with the DRP schemes is constructed. These conditions are derived from the asymptotic solutions of the governing equations. The asymptotic solutions are found by the use of Fourier-Laplace transforms and the method of stationary phase. A sequence of numerical simulations has been carried out. These simulation are designed to test the effectiveness of the DRP schemes and the radiation and outflow boundary conditions. The computed solutions agree very favorably with the exact solutions. The radiation boundary conditions perform satisfactorily causing little acoustic reflections. The outflow boundary conditions are found to be quite transparent to outgoing disturbances even when the disturbances are made up of a combination of acoustic, vorticity, and entropy waves. 26 refs., 14 figs.