No abstract

In this paper, the Chebyshev spectral method is used to solve the normal mode and parabolic equation models of underwater acoustic propagation, and the results of the Chebyshev spectral method and the traditional finite difference method are compared for an ideal fluid waveguide with a constant sound velocity and an ideal fluid waveguide with a deep-sea Munk speed profile. The research shows that, compared with the finite difference method, the Chebyshev spectral method has the advantages of a high computational accuracy and short computational time in underwater acoustic propagation.

No abstract

No abstract

The normal mode model is important in computational atmospheric acoustics. It is often used to compute the atmospheric acoustic field under a time-independent single-frequency sound source. Its solution consists of a set of discrete modes radiating into the upper atmosphere, usually related to the continuous spectrum. In this article, we present two spectral methods, the Chebyshev-Tau and Chebyshev-Collocation methods, to solve for the atmospheric acoustic normal modes, and corresponding programs are developed. The two spectral methods successfully transform the problem of searching for the modal wavenumbers in the complex plane into a simple dense matrix eigenvalue problem by projecting the governing equation onto a set of orthogonal bases, which can be easily solved through linear algebra methods. After the eigenvalues and eigenvectors are obtained, the horizontal wavenumbers and their corresponding modes can be obtained with simple processing. Numerical experiments were examined for both downwind and upwind conditions to verify the effectiveness of the methods. The running time data indicated that both spectral methods proposed in this article are faster than the Legendre-Galerkin spectral method proposed previously.

Solving an acoustic wave equation using a parabolic approximation is a popular approach for many existing ocean acoustic models. Commonly used parabolic equation (PE) model programs, such as the range-dependent acoustic model (RAM), are discretized by the finite difference method (FDM). Considering the idea and theory of the wide-angle rational approximation, a discrete PE model using the Chebyshev spectral method (CSM) is derived, and the code is developed. This method is currently suitable only for range-independent waveguides. Taking three ideal fluid waveguides as examples, the correctness of using the CSM discrete PE model in solving the underwater acoustic propagation problem is verified. The test results show that compared with the RAM, the method proposed in this paper can achieve higher accuracy in computational underwater acoustics and requires fewer discrete grid points. After optimization, this method is more advantageous than the FDM in terms of speed. Thus, the CSM provides high-precision reference standards for benchmark examples of the range-independent PE model.

Sound propagation in a range-dependent ocean environment has long been a matter of widespread concern in ocean acoustics. Stepwise coupled modes is one of the main methods to solve range-dependent acoustic propagation problems. Underwater sound propagation satisfies a Helmholtz equation, the solution of which represents the core of computational ocean acoustics. Due to its high accuracy in solving differential equations, the spectral method has been introduced into computational ocean acoustics in recent years and has achieved remarkable results. However, the existing underwater acoustic propagation algorithms based on the spectral method can calculate only range-independent ocean acoustic waveguides, which excludes applications in more general range-dependent environments. In this paper, a complete and efficient algorithm is designed using an improved global matrix of coupled modes to solve the range dependence of the ocean environment and using the Chebyshev-Tau spectral method to precisely solve the eigenmodes in a stepped range-independent stair. Based on this algorithm, a complete and efficient numerical program is developed, and the numerical simulation results verify that this algorithm is extremely computationally fast and accurate for various range dependence and seabed environments.

No abstract