水平变化波导中的全局矩阵耦合简正波方法

提出了一种新的水平变化波导中声场的耦合简正波求解方法,该方法能够处理二维点源和线源问题,提供声场的双向解。该方法利用全局矩阵(DGM)一次性求解耦合模式的系数,消除了传播矩阵递推求解中存在的误差累积问题;此外,改善了现有模型中对距离函数的归一化方法,从而避免了泄露模式指数增长导致的数值溢出问题。本文还给出了绝对软海底理想波导中耦合矩阵的闭合表达式,并分析了单个阶梯下简正波耦合现象。此外,本文还计算了理想楔形波导中的声传播问题(ASA标准问题),并与解析解及COUPLE07计算结果进行了比较,结果表明该方法是一种稳定、精确的水平变化波导中的声场计算方法。      

一种水平变化可穿透波导中声传播问题的耦合简正波方法

通过利用标准简正波程序KRAKEN计算本地简正波解及耦合矩阵, 进一步发展了求解水平变化波导中声场的全局矩阵耦合简正波方法(Luo et al., "A numerically stable coupled-mode formulation for acoustic propagation in range-dependent waveguides," Sci. China-Phys. Mech. Astron. 55, 572 (2012)), 使得该方法可以处理具有可穿透海底及随深度变化声速剖面等实际问题, 并提供声场的完全双向解. 本文还给出了双层波导中耦合矩阵的解析表达式, 并利用其验证了本方法中耦合矩阵数值算法的精度. 最后, 利用改善后的全局矩阵耦合简正波模型(DGMCM)计算了美国声学学会(ASA)提出的可穿透楔形波导标准问题, 将所得数值解与参考解比较, 结果表明DGMCM方法可以精确处理水平变化波导中声传播实际问题. 关键词： 耦合简正波理论 全局矩阵方法 可穿透楔形波导     

水平变化波导中声传播问题的高效稳定耦合简正波解

提出了一种高效、精确而且数值稳定的耦合简正波方法,用于求解水平变化波导中的声传播问题.与现有的数值模型相比,本方法采用合理的归一化距离解,极大地提高了数值稳定性;通过结合一个前向传播过程和一个后向传播过程,该方法极大地提高了计算精度,可用于精确求解水平变化波导中的声传播问题,尤其适用于大角度斜坡以及海底与海水声阻抗差别...     

水平变化双层波导声传播问题的耦合简正波解

为了验证现有模型的精度,导出了全反射下边界双层波导中简正波耦合矩阵的解析表达式,并将其应用到全局矩阵耦合简正波模型(Direct Global Matrix Coupled-Mode)中,使得该模型可以提供水平变化双层波导问题的标准解。文中首先利用COUPLE的简正波及耦合矩阵数值解验证了该简正波及耦合矩阵解析表达式的正确性;其次,采用改进的DGMCM模型求解了双层波导海山声传播损失,结果表明,改进后的DGMCM模型可以非常精确地求解水平变化双层波导问题,可作为求解此类问题的标准模型使用。     

水平变化波导中一种声场计算方法的数值实现

全局矩阵方法计算水平变化环境中的声场具有稳定性好、速度快、精度高等优点,在数值实现中如何快速、准确的求解大规模矩阵是该方法的一个关键问题。本文针对全局矩阵的特点,分别利用两种矩阵求解器,PARDISO和LAPACK,求解该问题。经过比较和讨论,得出结论:LAPACK对于全局矩阵的求解更有优势,求解速度快,而且可求解问题规模相对较大。使用不同的数值实现方法计算了楔形波导的传播损失(ASA标准问题),与解析解比较,证明全局矩阵方法计算精度很高。     

非均匀波导中的声聚焦

理论研究了声波在非均匀波导中的空间聚焦问题,利用多模态导纳法构建波导内任意位置处声压与入射声压在模态域的映射关系,计算使声波聚焦于空间某位置时的最佳入射波,并画出了相应的聚焦声场.研究了三种非均匀情况:水平变截面波导、含散射体波导以及声速垂直变化波导.结果表明,当输入最佳入射波时,在非均匀波导中可以产生良好的单点或多点声聚焦效果,声波的聚焦过程充分地利用了波导结构及介质非均匀性对声波的散射作用.     

非均匀波导中的最大声能流透射及鲁棒性分析

以变截面含可穿透散射体波导为模型,理论研究声波在非均匀波导中的最大透射问题.通过耦合简正波理论构建模态域内透射矩阵和水平波数矩阵,推导透射波能流的具体表达式,分析任意入射波的能流透射率随频率的变化,进而讨论任意给定频率下能够产生最大能流透射率的最佳入射波,并给出数组全透射声场算例.最佳入射波仅由可传播模态决定,与衰逝模态无关.利用衰逝模态不携带能流的特性,讨论衰逝模态对产生能流最大透射声场的影响,并分析最大能流透射的鲁棒性.在频率满足一定条件时,全透射声场可能表现出完美鲁棒性.文中所述方法可延伸至多种非均匀波导以分析其中的能流最大透射问题.     

复等效深度耦合简正波模型

为了考虑海底地形随距离变化的非水平分层介质中割线积分对声场的贡献,提出了复等效深度耦合简正波模型。该耦合简正波模型由介质运动方程和连续性方程推导得到了耦合微分方程组,此方程组满足海底地形随距离变化情况下的边界条件且仅包含一个耦合矩阵,并通过引入复等效深度理论处理连续谱和离散谱之间的相互耦合。仿真计算表明,复等效深度耦合简正波模型提高了波导简正波本征值位于割线枝点附近情况下声传播损失的计算精度,充分考虑了波导简正波、非波导简正波和割线积分对声场的贡献,可快速而准确地计算非水平分层介质中的声场。      

水平变化波导中两种耦合简正波方法的比较与分析

分析比较了处理水平变化波导中声传播问题的两种耦合简正波方法:DGMCM(Direct-Global-Matrix Coupled-Mode Method)和CCMM(Consistent Coupled-Mode Method)。首先,两种方法都提供声场的双向解,具有很高的计算精度。其次,DGMCM和CCMM中本地垂直模式序列均具有较快的收敛速度,但DGMCM比CCMM需要较少的水平分段数。再次,两种方法通过求解不同的耦合简正波系统得到声场解,但求解过程中所需参数的计算量基本相同。另外,DGMCM能够处理某些CCMM不容易解决的问题,如海底形状不平滑,线源在斜坡海底上方,以及多声源问题。在DGMCM方法中给出了双层波导问题的耦合矩阵解析表达式,还推导更新了CCMM模型,使其能够处理二维线源问题。     

水平变化波导中两种耦合简正波方法的比较与分析

分析比较了处理水平变化波导中声传播问题的两种耦合简正波方法:DGMCM(Direct-Global-Matrix Coupled-Mode Method)和CCMM(Consistent Coupled-Mode Method)。首先,两种方法都提供声场的双向解,具有很高的计算精度。其次,DGMCM和CCMM中本地垂直模式序列均具有较快的收敛速度,但DGMCM比CCMM需要较少的水平分段数。再次,两种方法通过求解不同的耦合简正波系统得到声场解,但求解过程中所需参数的计算量基本相同。另外,DGMCM能够处理某些CCMM不容易解决的问题,如海底形状不平滑,线源在斜坡海底上方,以及多声源问题。在DGMCM方法中给出了双层波导问题的耦合矩阵解析表达式,还推导更新了CCMM模型,使其能够处理二维线源问题。      

Sound propagation in inhomogeneous waveguides with sound-speed profiles using the multimodal admittance method

The multimodal admittance method and its improvement are presented to deal with various aspects in underwater acoustics, mostly for the sound propagation in inhomogeneous waveguides with sound-speed profiles, arbitrary-shaped liquid-like scatterers, and range-dependent environments. In all cases, the propagation problem governed by the Helmholtz equation is transformed into initial value problems of two coupled first-order evolution equations with respect to the modal components of field quantities(sound pressure and its derivative), by projecting the Helmholtz equation on a constructed orthogonal and complete local basis. The admittance matrix, which is the modal representation of Direchlet-to-Neumann operator, is introduced to compute the first-order evolution equations with no numerical instability caused by evanescent modes. The fourth-order Magnus scheme is used for the numerical integration of differential equations in the numerical implementation. The numerical experiments of sound field in underwater inhomogeneous waveguides generated by point sources are performed. Besides, the numerical results computed by simulation software COMSOL Multiphysics are given to validate the correction of the multimodal admittance method. It is shown that the multimodal admittance method is an efficient and stable numerical method to solve the wave propagation problem in inhomogeneous underwater waveguides with sound-speed profiles, liquid-like scatterers, and range-dependent environments. The extension of the method to more complicated waveguides such as horizontally stratified waveguides is available.     

A coupled-mode sound propagation model with complex effective depth

A coupled-mode sound propagation model with complex effective depth is presented,in order to involve the effect of branch line integral for acoustic field in a range-dependent waveguide.The equations of motion and continuity are used to obtain the coupled equations,which satisfy boundary conditions in the waveguide with varying topography and contain one coupling matrix.Meanwhile,the couplings between discrete and continuous spectrum are dealt with based on complex effective depth theory.Numerical simulations show that the accuracy of transmission loss is improved by the coupled mode model when eigenvalues of trapped modes are located near the branch point.The acoustic field in a non-horizontally stratified waveguide can be calculated efficiently and accurately by this model,and the energy corresponding to trapped modes,leaky modes and branch line integral can be considered adequately.     

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.     

水平变化波导中的简正波耦合与能量转移

针对海底地形水平变化对声场能量传播和声场干涉结构的影响,对简正波之间的耦合和能量转移进行了研究.建立了一种二维大步长格式的耦合简正波模型和三维楔形波导耦合简正波模型,以便快速有效地分析简正波之间的耦合和能量转移.基于耦合简正波模型,阐述了前向声场能量在水平变化波导中传播时的转移过程.并根据射线简正波理论,解释了海底地形变化对声场能量分布的影响机理.水平变化波导中声场的仿真计算表明,当本征值虚部发生剧烈变化时声场存在着较强的简正波耦合和能量转移,且海底地形变化将导致声场能量的水平传播方向偏转至海水深度增加的方向.在声场能量转移和传播方向变化中,声场的能量趋于保留在波导中而不向海底泄漏.同时,声场能量分布受到类似于压缩或稀疏的作用,从而形成椭圆状的干涉结构.     

Propagation characteristics of highly elliptical core optical waveguides: a perturbation approach

A simple first-order perturbation approach to study the propagation characteristics of highly elliptical core waveguides is presented. The unperturbed index profile is taken to be a pseudo-rectangular core waveguide for which exact scalar wave modes are known. The method is used to obtain the various propagation characteristics of elliptical core, single mode waveguides. Results for various propagation characteristics obtained by the present analysis agree very well with those reported by other authors obtained by accurate numerical techniques.     

The contribution of a lateral wave in simulating low-frequency sound fields in an irregular waveguide with a liquid bottom

A further development of a previously proposed approach to calculating the sound field in an arbitrarily irregular ocean is presented. The approach is based on solving the first-order causal mode equations, which are equivalent to the boundary-value problem for acoustic wave equations in terms of the cross-section method. For the mode functions depending on the horizontal coordinate, additional terms are introduced in the cross-section equations to allow for the multilayer structure of the medium. A numerical solution to the causal equations is sought using the fundamental matrix equation. For the modes of the discrete spectrum and two fixed low frequencies, calculations are performed for an irregular two-layer waveguide model with fluid sediments, which is close to the actual conditions of low-frequency sound propagation in the coastal zone of the oceanic shelf. The calculated propagation loss curves are used as an example for comparison with results that can be obtained for the given waveguide model with the use of adiabatic and one-way propagation approximations.