海洋声学中三维抛物方程非均匀网格模型

在海洋声学中,三维抛物方程模型可以有效考虑三维空间的声传播效应。然而,采用三维抛物方程模型分析三维空间内的声传播问题时,计算时间较长,并且需要消耗较大的计算机内存,因此给远距离声场的快速精确计算带来了很大困难。为此,将非均匀网格Galerkin离散化方法用于三维直角坐标系下的水声抛物方程模型中,深度算子和水平算子Galerkin离散方式由均匀网格变为非均匀网格。仿真结果表明,三维直角坐标系下非均匀网格离散的抛物方程模型,在保持计算精度、提高计算速度的同时,可以实现远距离声场的快速预报。另外,针对远距离局部海底地形与距离有关的三维声传播问题,给出了声场快速计算方法;在海底保持水平的区域,采用经典Kraken模型,重构抛物方程算法的初始场,随后依次递推求解地形与距离有关海底下的三维声场。采用改进模型,证明了远距离楔形波导声强增强效应。     

Developments of parabolic equation method in the period of 2000–2016

Parabolic equation(PE) method is an efficient tool for modelling underwater sound propagation, particularly for problems involving range dependence. Since the PE method was first introduced into the field of underwater acoustics,it has been about 40 years, during which contributions to extending its capability has been continuously made. The most recent review paper surveyed the contributions made before 1999. In the period of 2000–2016, the development of PE method basically focuses on seismo-acoustic problems, three-dimensional problems, and realistic applications. In this paper, a review covering the contribution from 2000 to 2016 is given, and what should be done in future work is also discussed.     

Wave mode computing method using the step-split Padé parabolic equation

Models based on a parabolic equation (PE) can accurately predict sound propagation problems in range-dependent ocean waveguides. Consequently, this method has developed rapidly in recent years. Compared with normal mode theory, PE focuses on numerical calculation, which is difficult to use in the mode domain analysis of sound propagation, such as the calculation of mode phase velocity and group velocity. To broaden the capability of PE models in analyzing the underwater sound field, a wave mode calculation method based on PE is proposed in this study. Step-split Padé PE recursive matrix equations are combined to obtain a propagation matrix. Then, the eigenvalue decomposition technique is applied to the matrix to extract sound mode eigenvalues and eigenfunctions. Numerical experiments on some typical waveguides are performed to test the accuracy and flexibility of the new method. Discussions on different orders of Padé approximant demonstrate angle limitations in PE and the missing root problem is also discussed to prove the advantage of the new method. The PE mode method can be expanded in the future to solve smooth wave modes in ocean waveguides, including fluctuating boundaries and sound speed profiles.     

采用能量守恒和高阶Padé近似的三维水声抛物方程模型

为了充分考虑海底地形随三维空间变化的海洋环境中水平方位角耦合效应对声传播的影响,建立了一种三维柱坐标系下流体高阶抛物方程算法。该算法采用泰勒近似将二维方根算子分裂成一维方根算子,并采用分裂步进的高阶Pade近似将一维方根算子写成微分算子有理分式连乘的形式,进而应用Galerkin离散化方法来处理微分算子,最终将微分方程写成矩阵方程的形式;采用能量守恒近似来处理海底边界,以考虑复杂海底对于声传播的影响;采用交替方向隐式格式,实现了三维声场的步进计算。楔形和海底山等典型海域声场仿真计算表明,相比于已有的声场计算模型,三维柱坐标系下高阶抛物方程模型可以更加精确地计算楔形海域和海底山区域的三维声场,实现水平方位全空间声场计算。