A Chebyshev spectral element approximation of acoustic propagation problems based on linearized Euler equations is introduced, and the numerical approach is based on spectral elements in space with first-order Clayton–Engquist–Majda absorbing boundary conditions and implicit Newmark method in time. An initial perturbation problem has been solved to test the accuracy and stability of the numerical method. Then the sound propagation by source terms is also studied, including the radiation of a monopole and dipolar source in both stationary medium and uniform mean flow. The numerical simulation leads to good results in both accuracy and stability. Compared with the analytical solutions, the numerical results show the advantages in spectral accuracy even with relatively fewer grid points. Moreover, the implicit Newmark method in time marching also presents its superiority in stability. Finally, a problem of sound propagation in pipes is simulated as well.
SPECTRAL ELEMENT APPROXIMATION OF FREDHOLM INTEGRAL EIGENVALUE PROBLEMS
