一种可稳定计算Pekeris波导中声场的波数积分方法

提出一种可稳定计算Pekeris波导中声场的波数积分方法,并在此基础上开发出一个数值模型,可用于提供Pekeris波导中声场的精确、稳定的数值解。在这个方法中,由于与深度有关的波动方程齐次解中所有的上行波与下行波均采用了合理的归一化表示,从而得到的系统方程是无条件稳定的。在简正波方法中,割线积分一般只对近场有显著影响。因此,传统的简正波模型一般都忽略割线积分对声场的贡献。但是,如果某号简正波离割线非常近,则割线积分对非常远距离的声场仍可能有显著影响。在这种情况下,传统的简正波模型由于忽略割线积分的贡献,从而得到的声场结果是不准确的。本文通过数值算例比较本文提出的波数积分模型与传统的简正波模型。数值结果表明,本文提出的模型可以提供精确、稳定的Pekeris波导中声场的数值解,而在某些情况下传统的简正波模型得到的声场结果是不准确的。因此,本文提出的模型可以作为Pekeris波导中声传播问题的标准模型使用。      

Propagation in an elastic wedge using the virtual source technique

The virtual source technique, which is based on the boundary integral method, provides the means to impose boundary conditions on arbitrarily shaped boundaries by replacing them by a collection of sources whose amplitudes are determined from the boundary conditions. In this paper the virtual source technique is used to model propagation of waves in a range-dependent ocean overlying an elastic bottom with arbitrarily shaped ocean-bottom interface. The method is applied to propagation in an elastic Pekeris waveguide, an acoustic wedge, and an elastic wedge. In the case of propagation in an elastic Pekeris waveguide, the results agree very well with those obtained from the wavenumber integral technique, as they do with the solution of the parabolic equation (PE) technique in the case of propagation in an acoustic wedge. The results for propagation in an elastic wedge qualitatively agree with those obtained from an elastic PE solution.     

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

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

On the acoustic field in a Pekeris waveguide with attenuation in the bottom half-space

The acoustic field in a Pekeris channel with an attenuating basement is critically examined, based on contour integrations of the wave number integrals for the two domains. In both regions, the field consists of a finite sum of proper (square integrable) normal modes plus a branch line integral around a hyperbolic cut. For low bottom attenuation, only "trapped" modes exist but as the loss increases additional "dissipation" modes contribute to the mode sum. A Newton-Raphson iterative procedure is introduced for finding the complex eigenvalues of the modes and a new expression is derived which shows that the total number of proper (trapped plus dissipation) modes supported by the waveguide increases essentially linearly with rising bottom attenuation. Approximations for the complex shape functions of the modes in the water column and the basement are developed and compared with the exact shape functions determined from the Newton-Raphson procedure. An expression derived for the modal attenuation coefficient scales in proportion to the square of the mode number and inversely with the square of the frequency. Stationary-phase approximations for the branch line integrals for both domains are developed, which serve to illustrate the asymptotic range dependence of the lateral wave but otherwise have little utility.     

三维浅海下弹性结构声辐射预报的有限元-抛物方程法

利用多物理场耦合有限元法对结构和流体适应性强、抛物方程声场计算高效准确的特点,提出了三维浅海波导下弹性结构声振特性研究的有限元-抛物方程法.该方法采用多物理场耦合有限元理论建立浅海下结构近场声辐射模型,计算局域波导下结构声振信息,并提取深度方向上复声压值作为抛物方程初始值;然后采用隐式差分法求解抛物方程以步进计算结构辐射声场.重点介绍了该方法对浅海下结构声辐射计算的准确性、高效性以及快速收敛性后,对Pekeris波导中有限长弹性圆柱壳的声振特性进行了分析.研究得出,当圆柱壳靠近海面(海底)时,其耦合频率比自由场下的要高(低),当潜深达到一定范围时,与自由场耦合频率基本趋于一致;在低频远场,结构辐射场与同强度点源声场具有一定的等效性,且等效距离随着频率增加而增加;由于辐射声场受结构振动模态、几何尺寸和简正波模式影响,结构辐射场传播的衰减规律按近场声影响区、球面波衰减区、介于球面波和柱面波衰减区、柱面波衰减区四个扩展区依次进行.     

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.     

A method of images for a penetrable acoustic waveguide

In this paper the complex-image approximation to the reflection coefficient for water over a seabed half-space is used to generate an image representation for a bounded acoustic waveguide with an underlying layered seabed. The images are true point sources; they have constant amplitudes which are raypath independent and, in the case of a Pekeris waveguide, frequency-independent. This image representation is ideal for constructing the Green's function kernel of the boundary integral equation method for target scattering in a waveguide. The singular behavior of the Green's function for an infinitesimal source/receiver separation, possibly with the target adjacent to one of the interfaces, is modeled correctly and the image expansion has a simple analytic form which can be analytically differentiated. The method is also accurate for significant source/receiver separations, which means that it can be used in the modeling of scattering from large-sized objects and can also be used as an efficient and accurate short-range propagation model for harmonic and broadband propagation in a penetrable waveguide.     

一种基于简正波模态消频散变换的声源距离深度估计方法

针对浅海环境中传播的低频宽带水声脉冲信号,基于简正波水平波数差和波导不变量之间的关系,本文提出了一种利用距离-频散参数二维平面聚焦测距与匹配模态能量定深的目标声源定位方法.首先,通过将由频散参数和波导不变量表示的前几阶模态相速度与由环境模型计算的相速度进行对比分析,从而估计出前几阶模态的频散参数和环境的波导不变量.其次,利用估计出的频散参数值和波导不变量对接收信号进行消频散变换处理,只有当接收信号的距离参数等于目标声源距离时,各号简正波的幅度均达到最大值,在距离-频散参数二维平面上,出现声压聚焦的现象,利用此现象可以估计目标声源的距离.不仅如此,消频散变换后的接收信号,前几阶模态在时域上明显地分离开来,可以准确地估计出前几阶模态的能量,采用多模态能量匹配的方式,可以估计出目标声源的深度.最后,通过对仿真和冬季获得的气枪信号数据处理结果验证了本文方法的有效性.     

Scattering in a Pekeris waveguide from a rough bottom using a two-way coupled mode approach

Scattering from a rough ocean bottom is described numerically with a two-way coupled-mode formalism that contains scattering effects to all orders and provides an exact solution to the wave equation. Both scattered field and direct blast components are computed within the formalism framework. A comparison of the scattered component solution from the coupled mode with the Born approximation (BA) solution for scattering from a rough bottom Pekeris waveguide shows that the BA predicts correctly the scattered field levels but not detailed structure. The transition from direct blast to scattered field dominance is identified in the total field time series.     

基于双水听器的目标深度辨识

随着水下目标越来越安静,探测并区分水下目标变得越来越困难。现在对目标的深度辨识多采用匹配场的方法,但其需要大量的先验知识,由于海洋情况多变复杂,这种方法并不适合进行工程实践。本文旨在探索出一种适用于浅海波导的目标深度辨识方法,并减轻工程实践的复杂度。此方法基于Pekeris波导模型,应用了数值仿真技术对垂直双水听器声压互谱的有功分量进行了分析,发现当声源深度变化时互谱有功分量符号随之变化,这个变化规律可以用于目标深度辨识。本文总结了简正波阶数以及水听器深度变化对目标深度辨识结果的影响。通过实验验证,此方法可以较好的对浅海目标进行深度辨识,方法较为可靠,有一定的应用前景。     

Piston-type radiator in a soft screen in the pekeris waveguide

A model problem is considered for a radiator in the form of a circular disk with a given pressure jump at its surface. The radiator is inserted in a soft screen coinciding with the upper boundary of the Pekeris waveguide. A series expansion of the sound field in normal modes is obtained. A numerical analysis of the radiation impedance and its components that are responsible for the radiation into the waveguide and into the halfspace is carried out.     

浅海波导中目标散射的简正波-Kirchhoff近似混合方法

提出了计算浅海波导中复杂目标散射的数值方法:简正波-Kirchhoff近似混合方法。通过把目标散射的Kirchhoff近似方法和简正波声传播模型相结合,可对大尺寸复杂目标在浅海波导中的散射声场进行计算。以浅海波导中刚性球体散射的解析解为标准解验证了本方法,说明简正波-Kirchhoff近似混合方法是有一定计算精度的工程预报方法。数值计算Pekeris浅海波导中球体目标散射声场在深度-频率平面和距离-频率平面上的二维干涉结构及其与自由空间中的差异。进而通过FFT获得目标时域回波随深度的分布图,分析浅海波导中目标姿态、声速剖面对Benchmark潜艇目标散射的影响。      

Depolarization of the electromagnetic field scattered by a three-dimensional large-scale irregularity of the lower ionosphere

This paper develops analytical and numerical methods for the solution of three-dimensional problems of radio wave propagation. We consider a three-dimensional vector problem for the electromagnetic field of a vertical electric dipole in a planar Earth-ionosphere waveguide with a largescale local irregularity of negative characteristics at the anisotropic ionospheric boundary. The field components at the boundary surfaces obey the Leontovich boundary conditions. The problem is reduced to a system of two-dimensional integral equations taking into account the overexcitation and depolarization of the field scattered by the irregularity. Using asymptotic (with respect to the parameter kr≫1, where r is the distance from the source or receiver to the nearest point of the irregularity, k=2π/λ, and λ is the radio wavelength) integration over the direction perpendicular to the ray path, we transform this system to one-dimensional integral equations where integration contours represent the geometric contour of the irregularity. The system is numerically solved in the diagonal approximation, combining direct inversion of the Volterra integral operator and subsequent iterations. The proposed numerical algorithm reduces the computer time required for the solution of this problem and is applicable for studying both small-scale and large-scale irregularities. We obtained novel estimates for the field components that are not excited by the source but result entirely from scattering by the sample three-dimensional ionospheric irregularity.     

Spectral and modal formulations for the Doppler-shifted field scattered by an object moving in a stratified medium

Spectral and normal mode formulations for the three-dimensional field scattered by an object moving in a stratified medium are derived using full-field wave theory. The derivations are based on Green's theorem for the time-domain scalar wave equation and account for Doppler effects induced by target motion as well as source and receiver motion. The formulations are valid when multiple scattering between the object and waveguide boundaries can be neglected, and the scattered field can be expressed as a linear function of the object's plane wave scattering function. The advantage of the spectral formulation is that it incorporates the entire wave number spectrum, including evanescent waves, and therefore can potentially be used at much closer ranges to the target than the modal formulation. The normal mode formulation is more computationally efficient but is limited to longer ranges. For a monochromatic source that excites N incident modes in the waveguide, there will be roughly N2 distinct harmonic components in the scattered field. The Doppler shifts in the scattered field are highly dependent upon the waveguide environment, target shape, and measurement geometry. The Doppler effects are illustrated through a number of canonical examples.     

Asymptotic Analysis of Vertical Branch-Cut Integral of Shear Waves in a Fluid-Filled Borehole Utilizing the Steepest-Descent Method

We obtain an asymptotic solution to the vertical branch-cut integral of shear waves excited by an impulsive pressure point source in a fluid-filled borehole, by taking the effect of the infinite singularity of the Hankel functions related to shear waves in the integrand at the shear branch point into account and using the method of steepest-descent to expand the vertical branch-cut integral of shear waves. It is theoretically proven that the saddle point of the integrand is located at ks-i/z, where ks and z are the shear branch point and the offset. The continuous and smooth amplitude spectra and the resonant peaks of shear waves are numerically calculated from the asymptotic solution. These asymptotic results are generally in agreement with the numerical integral results. It is also found by the comparison and analysis of two results that the resonant factor and the effect of the normal and leaking mode poles around the shear branch point lead to the two-peak characteristics of the amplitude spectra of shear waves in the resonant peak zones from the numerical integral calculations.     

Scattering of plane guided waves obliquely incident on a straight feature with uniform cross-section

A frequency-domain finite element (FE) method is presented for modeling the scattering of plane guided waves incident on an infinitely-long, straight feature with uniform cross-section in a planar host waveguide. The method utilizes a mesh of 2-dimensional finite elements with harmonic shape functions in the perpendicular direction. The model domain comprises a cross-section through the feature and short lengths of the adjoining host waveguide. A spatial frequency equal to the wavenumber of the desired incident mode multiplied by the sine of the desired incidence angle is prescribed for the element shape functions. An integral representation of the incident mode is used to determine a suitable system of harmonic forces to uniquely excite that mode. These are applied at nodes through the thickness of the host waveguide on one side of the feature. The displacement field is measured at nodes through the thickness of the host waveguide on either side of the feature and decomposed into reflected and transmitted modes. The cases of guided wave transmission in a featureless waveguide and the reflection of guided waves from a free-edge are examined as validation cases. Finally, the results for transmission at an adhesively-bonded stiffener are presented and compared with experimental measurements.     

A method for multi-frequency calculation of boundary integral equation in acoustics based on series expansion

The paper presents a method to solve the problem of multi-frequency calculation of Helmholtz boundary integral equation in acoustics. Based on series expansion, system matrices are independent of wavenumber and become the matrix power series of wavenumber. As a result, all matrices in the matrix power series are only dependent on the structure geometry. In addition, an element transform method to calculate the singular integral and Cauchy singular integral is also discussed because the singular integral need to be solved using the method. The convergence of the series expansion method is also proved in this paper. The effectiveness of the method is confirmed by two numerical examples.     

Scattering of sound by an elastic spherical shell immersed in a waveguide with a fluid bottom

The paper is devoted to modeling an acoustic field scattered by an elastic spherical shell. The shell is immersed in a waveguide with a fluid attenuating bottom. Modal analysis is applied to calculations of both the point source field in a free waveguide and the scattering coefficients of the shell. The occurring integrals along a branch cut are expressed via the probability integral. Cases are analyzed when the source frequency differs from the critical frequencies of normal waves and when it coincides with one of them. The method is applied to estimate the effect of the field scattered by a thin elastic shell filled with air with transmission loss of the total acoustic field with distance and its depth distribution.     

Pekeris波导中简正波的复声强及其应用

在Pekeris波导模型下,关注了简正波的矢量场,讨论了简正波水平复声强和垂直复声强的表述,并分析了其特征.单阶简正波在水平方向是行波,相应的水平复声强仅为有功的;在垂直方向为驻波,相应的垂直复声强仅为无功的.而多阶简正波相互干涉,因此总声场的复声强既有有功分量,也有无功分量,其中只有有功分量参与声能的输运,但无功分量是反映声场信息的重要组成部分.通过对垂直(交互)复声强无功分量和水平交互复声强有功分量的数值分析,对于甚低频率的点源声场,发现当声源深度变化时,上述声场分量的正负号呈有规变化,当接收传感器置 关键词： 目标深度分类 复声强 矢量场 Pekeris波导     

Exact numerical solution for the diffused channel waveguide

We present an exact numerical solution based on a vector integral equation to investigate the diffused channel waveguide. Various examples are given and compared with other numerical and approximate methods. An approximate numerical solution is also given using the effective index method and the integral equation solution for the inhomogeneous slab.