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

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

空气中快速运动声源水下声场的波数积分模型

本文采用二维波数积分，对空气中高速运动声源激发的水下声场进行建模。针对二维波数积分计算声场时域解计算量大的问题，提出一种快速计算方法。用本文提出的方法，对深海和浅海情况下，空气中高速运动单频点声源激发的水下声场进行了计算和仿真。计算结果表明：在深海，水下接收信号的幅度和瞬时频率随时间发生变化；接收器深度、接收器与声源运动轨迹的最小距离对接收信号的变化快慢有较大影响，而声源高度的影响较小；在浅海中，接收信号呈现快速的幅度起伏，明显的多普勒频移和大的频率展宽效应。与简正波方法相比，本文方法主要适用于近场计算，而简正波方法适用于远场。另外，当声源频率较高时，二维波数积分方法的计算量将迅速增大。     

管内声场计算的一种新的差分方法

本文提出一种新的计算管内声传播的差分方法。它采用MacCormack时间分裂显式格式,直接求解声压和声学速度的基本方程组,通过合适处理声学边界条件,成功地计算了三维矩形管内声场。文中还给出了全部计算格点的差分格式和稳定性条件。这种差分方法的特点是内存少,计算比较稳定,程序简单,调试方便。特别适合于三维声场计算。     

平面光波导放大器中光场分布的一种计算方法

用交变隐式差分方向的时域有限差分法(ADI-FDTD)计算了掺Er³⁺玻璃平面光波导放大器中光场的分布.其算法简单,结果无条件稳定,能完全模拟波导中的光场特性.     

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

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

快速计算阵列波导光栅波导耦合系数的修正重叠积分方法

提出了一种修正的重叠积分方法用以计算阵列波导光栅(AWG)波分复用器件中光从自由传输区域到阵列波导的耦合系数,并和光束传播法(BPM)数值方法计算得到的结果做了比较通过比较分析,得出结论:当波导中心距不是太小时,用修正重叠积分这一快速方法是合适的.     

基于声粒子分布积分的无网格声场计算方法

小尺度封闭空间内部声场的数值计算是声学设计、噪声控制等领域的关键技术。由于波动声学及几何声学方法计算频率上的限制,中频段声场计算问题一直是个难点。本文以声学无网格法为基础,提出了一种基于声粒子分布积分的无网格声场数值计算方法。文中利用声线跟踪理论计算声场中的声粒子分布,并以某个时间点上的声粒子作为蒙特卡罗法中的积分点,将其应用于无网格法中,从而获得声场中的节点声压。利用该方法对一个矩形封闭空间的中低频声场进行了计算,并与模态叠加法、商用声场计算软件、经典无网格法的结果进行了对比,证明基于声粒子分布积分的无网格声场数值计算方法在中低频段相较于传统基于网格的方法具有更高的精度。     

Pekeris波导中垂直发射阵的阵增益快速估计

为了利用单换能器的声强级估计垂直发射阵的声强级,研究了Pekeris波导中快速估计发射阵增益的方法。首先,基于等效深度简正波理论,计算了Pekeris波导中垂直发射阵的平均声强,进而得到用模态阶数求和表示的发射阵增益的表达式。然后对模态阶数用积分代替求和,并对引入的误差进行修,得到3/2次方衰减场区发射阵增益的快速估计表达式;有效声影区声场由于只剩1阶简正模态,可直接简化得到发射阵增益的快速估计表达式,在海底参数未知时可用海深近似代替等效深度来估计有效声影区的发射阵增益。计算机仿真和海上实验得到的阵增益快速估计误差均小于1 dB,验证了快速估计方法的有效性,这一发射阵增益快速估计方法可以为工程应用提供参考。     

一种计算固体中超声脉冲散射声场的边界元法

本文分析了用边界元方法计算固体中超声的散射,提出了一种时间域超声散射计算方法,这种方法能动态显示超声散射的全过程,并和动态光弹照片进行了对比。     

Pekeris波导中简正波声强流及其互谱信号处理

探讨了Pekeris波导中声压场和质点振速场的联合描述,尤其关注垂直声强流的分析。研究表明,由于各阶简正波的干涉作用,水平和垂直声强流均既有有功分量,亦有无功分量。低频声场的垂直声强流的无功分量虽不参与声能的输运,当单个声矢量传感器做适当放置后却可用以判断声源的特定深度。这对于矢量信号处理是有意义的。给出了声压和质点振速联合互谱处理进行目标判别的算法。     

Analytical solution based on the wavenumber integration method for the acoustic field in a Pekeris waveguide

An exact solution based on the wavenumber integration method is proposed and implemented in a numerical model for the acoustic field in a Pekeris waveguide excited by either a point source in cylindrical geometry or a line source in plane geometry. Besides, an unconditionally stable numerical solution is also presented, which entirely resolves the stability problem in previous methods. Generally the branch line integral contributes to the total field only at short ranges, and hence is usually ignored in traditional normal mode models. However, for the special case where a mode lies near the branch cut, the branch line integral can contribute to the total field significantly at all ranges. The wavenumber integration method is well-suited for such problems. Numerical results are also provided, which show that the present model can serve as a benchmark for sound propagation in a Pekeris waveguide.     

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.     

Piston acoustic radiator in an impedance screen of a Pekeris waveguide

A model problem of seismic hydroacoustics is studied for a piston radiator inserted in an impedance screen that coincides with the lower boundary of a Pekeris waveguide. The radiation resistance of the piston is numerically calculated as a function of the screen type and parameters of the radiator and the waveguide.     

Pekeris waveguide comparisons of methods for predicting acoustic field amplitude uncertainty caused by a spatially uniform environmental uncertainty (L)

Acoustic field calculations in underwater environments are often uncertain because the environmental parameters required for such calculations are uncertain. This letter compares the accuracy of direct simulations, the field shifting approximation, and polynomial chaos expansions for predicting acoustic amplitude uncertainty in 100-m-deep Pekeris waveguides having spatially uniform uncertain water-column sound speed. When this sound speed is Gaussian-distributed with a standard deviation of 1 m/s, direct simulations and polynomial chaos expansions, based on 21 field calculations, are more accurate than the field shifting approximation, based on two field calculations. This ranking reverses as the sound-speed standard deviation increases to 20 m/s.     

Waveguide invariant estimation in elastic Pekeris waveguide

Focus on sound intensity,particle velocity,complex sound intensity and waveguide impedance,the formation mechanism of interference structures in range-frequency domain for the above four kinds of physical quantities in Pekeris waveguides is analyzed based on Normal Mode Theory.Also,a sea-trial with broadband radiated source is conducted to verify the analysis results.Both the simulation results and the analysis of sea-trial data indicate that,in range-frequency domain,the four kinds of physical quantities will exhibit a stable interference structure which can be expressed with waveguide invariant β.Among these quantities,the coherent components of complex sound intensity can better reflect the interference characteristics of the waveguide.Finally,a multi-scale linear filter is introduced to deal with the sea-trial LOFAR(Low Frequency Analysis Recording) spectrum,the processing results show that the proposed filter can effectively enhance the interference characteristics in images,and conductively detect and extract striations' information from interference structures in LOFAR spectrum.     

Relating waveguide invariant and bottom reflection phase-shift parameter P in a Pekeris waveguide

The waveguide invariant β is affected by the shallow-water environment. The effect due to bottom sediment on β is investigated in this paper. It is found that the effect of sediment bottom can be concentrated on one parameter P-the bottom reflection phase-shift parameter. For a Pekeris waveguide, under Wentzel-Kramers-Brillouin (WKB) approximation, a very simple analytic relation is given: β ≈ 1 + P/(k(0)H(eff)), where H(eff) is the "effective depth," and H(eff) = H + P/2 k(0). The value of β related to different high-speed sediments (including layered sediment) ranges from 1.0 to 1.5. Some numerical examples including the layered sediment case are conducted to verify this result. Good agreement between the results calculated by KRAKEN and by WKB with parameter P has been found. Hence, the application of parameter P provides a model-free platform to investigate the bottom effect on the waveguide invariant β in shallow-water.     

Power emitted by a vertical compensated linear array in a Pekeris waveguide

An expression is derived for the acoustic power emitted by a vertical compensated discrete linear array in a Pekeris waveguide. The sound field is represented by the sum of a discrete spectrum and a continuous one. We consider the dependence of the power emitted by the array on both the number of array elements and the array compensation angle.