An effective algorithm for simulating acoustical wave propagation. H.M. Sun, J.B. Wang, J. Pan.

PROGRAM SUMMARY
Title of program: AWP.f90, v 1.0
Catalogue identifier: ADRE
Ref. in CPC: 151(2003)241
Distribution format: gzip file
Number of lines in distributed program, including test data, etc: 777
Keywords: Acoustical wave propagator, Sound wave propagation, Sound pressure in duct, Computational methods, Differential equations.
Programming language used: Fortran

Nature of physical problem:
In this paper, the acoustical wave propagator scheme is implemented in Fortran for predicting sound propagation in a one-dimensional duct. Example calculations are performed for a semi-infinite duct and a duct with a solid blockage. Numerical accuracy of our results is examined and compared with the finite-difference time-domain method. This scheme is found to be highly accurate and computationally effective for describing the time-domain evolution of acoustic waves. Multiple reflections within the solid blockage and phase changes of the transmitting wave from solid back into air are illustrated through the implementation of this scheme.

Method of solution:
The acoustical wave propagator is used to describe sound wave propagation in a duct, which maps the spatially distributed sound pressure at t0 to that at t>t0. The propagator is expanded by real Chebyshev polynomials and spatial derivatives are evaluated by a Fast Fourier Transformation scheme.