An improved acoustical wave propagator method and its application to a duct structure

The pseudospectral time-domain method has long been used to describe the acoustical wave propagation. However, due to the limitation and difficulties of the fast Fourier transform (FFT) in dealing with nonperiodic problems, the dispersion error is inevitable and the numerical accuracy greatly decreases after the waves arrive at the boundary. To resolve this problem, the Lagrange-Chebyshev interpolation polynomials were used to replace the previous FFT, which, however, brings in an additional restriction on the time step. In this paper, a mapped Chebyshev method is introduced, providing the dual benefit of preserving the spectral accuracy and overcoming the time step restriction at the same time. Three main issues are addressed to assess the proposed technique: (a) Spatial derivatives in the system operator and the boundary treatment; (b) parameter selections; and (c) the maximum time step in the temporal operator. Furthermore, a numerical example involving the time-domain evolution of wave propagation in a duct structure is carried out, with comparisons to those obtained by Euler method, the fourth-order Runge-Kutta method, and the exact analytical solution, to demonstrate the numerical performance of the proposed technique.