用物理声学方法计算界面附近目标的回波

用物理声学（KIRCHHOFF近似）方法研究界面附近目标的声散射。导出了计算形态函数或目标强度的积分表示式．所产生的散射场由三部分组成：（1）不存在界面时目标的散射场;（2）入射或散射波之一经受界面反射时的散射场,如同双站散射场;（3）目标相对于界面的镜象产生的散射场,其中入射和散射波都经受界面反射．对于一个自由界面下的刚硬球,数值比较说明本文的方法与严格的解析方法符合良好．本文提供了一种计算沉浸在海面下或躺在海底上的水下目标的回波计算方法．     

适用于任意形状封闭空间的声场预测相干模型

针对现有几何声学的方法对封闭空间内声场进行预测时在中低频段出现较大误差的问题,该文提出一种近似圆锥声束追踪法和相干反射场理论相结合的声场预测新模型。在近似圆锥声束追踪法基础上,考虑声束轴线在边界多次反射时声压和相位的改变,最后计算不同声波之间的干涉效应,建立一种适用于任意形状封闭空间的声场预测相干模型。利用该模型对某一矩形封闭空间进行声场预测,通过对边界元法、Raynoise软件相干和非相干算法的预测结果和本模型的数值模拟结果对比。结果表明,文中提出的方法和边界元法的计算结果在中低频段非常吻合,两者的计算结果平均绝对误差为1.1 d B。本模型在中低频率下与同样考虑了相位的Raynoise相干算法相比有更好的准确性,在较高频率上,本模型计算结果与Raynoise相干算法计算结果非常吻合。     

水下环形凹槽圆柱体散射声场空间指向性调控

本文提出一种具有深度梯度的环形凹槽结构,可用于调控水中有限长刚性圆柱体散射声场空间指向性.基于声学相位阵列理论分析了环形凹槽圆柱声散射空间指向性改变的机理,研究表明:凹槽深度方向相位延迟和凹槽间Bragg散射的相互作用使得平面声波垂直于圆柱方向入射其正横方向散射声波发生偏转.采用有限元方法讨论了凹槽结构参数如占空比、梯度等对圆柱散射声场空间分布特征的影响规律.多个不同深度梯度环形凹槽单元组合圆柱体散射声场数值计算和实验结果显示:具有环形凹槽结构圆柱体正横方向散射声波均匀偏转到预定的空间范围内,使得圆柱体声散射场空间指向性均衡化,改变了圆柱整体的散射特征,这为水下目标声隐身设计和声波定向传播提供了新的方法.     

利用声辐射模态重构任意目标的散射声场

水下目标散射声场的重构可以作为水下目标散射特性的研究基础。本文主要利用声辐射模态对水下目标进行散射声场重构研究。首先,在借助声传递矩阵给出的任意结构声辐射模态的流体域求解方法基础上,通过理论证明了目标的散射声压与声辐射模态具有函数关系。其次,借助声场分布模态的概念,同时考虑到声场分布模态病态及声压测量易受噪声污染,提出基于声辐射模态的正则化散射声场重构算法。仿真结果表明,波数越低,重构所需声辐射模态阶数越少,在较高波数时仅需总模态数的大约20%即可对声场进行重构。与基于边界元的声场重构算法相比,计算量减小了至少80%,且克服了赫姆霍兹积分方程最小二乘法仅对球壳结构的重构效果较好而不适用于长条形结构重构的缺陷。     

水下掩埋目标的散射声场计算与实验

水下掩埋目标声散射问题是识别和探测掩埋目标的理论基础, 是声散射研究领域的热点问题. 本文基于射线声学推导了掩埋情况下目标声散射计算的格林函数近似式, 并在此基础上进一步给出了相应的远场积分公式. 在有限元方法的基础上, 将推导得到的公式写入有限元仿真软件, 对软件功能进行拓展, 构建二维轴对称目标的声散射模型, 并计算掩埋情况下弹性实心球在不同条件下的目标强度, 获得了其散射声场随频率、掩埋深度、沙层吸收系数等参数的变化规律. 开展实心球的自由空间和浅掩埋条件下水池声散射实验, 利用共振隔离技术处理实验数据, 提取目标声散射的纯弹性共振特征进行分析, 结果表明可将其用于掩埋目标识别和探测. 最后利用总散射声场与理论计算结果进行对比, 验证了理论仿真的正确性.     

水下低频振荡涡流场声散射调制机理与特性研究

水下涡流场对声波的散射问题是声波在复杂流场中传播的基本问题,在水下目标探测和流场声成像领域具有重要意义.针对水下低频振荡涡流场声散射调制问题建立了理论分析模型与数值计算方法,探究了其声散射调制声场的产生机理与时空频特性.首先,基于运动介质的波动方程,通过引入势函数将波动方程分解为流声耦合项和非耦合项,并对流声耦合项进行频域分析处理,揭示了水下振荡涡流场的声散射调制机理;其次,采用间断伽辽金数值方法对水下低频振荡涡流场中声传播过程进行了数值模拟,分析了低马赫数条件下,不同入射声波频率、涡流场的振荡频率和涡核尺度对涡流场声散射调制声场时空频特性的影响规律,并结合理论分析模型对其特性进行了解释.研究表明:低马赫数下,振荡涡流场对声波的散射可产生包含涡流场振荡频率双边带调制谐波的散射调制声场,且随着入射声波频率、涡核尺度的增大,散射调制声场强度增强,总散射声场空间分布具有对称性和明显主瓣,且主瓣方位角趋近于入射波传播方向;在频率比远大于1条件下,涡流场振荡频率对散射调制声场强度影响较小.     

水下等离子体声源的聚束特性研究

分析水下等离子体声源的反射聚束特性,对研究和应用水下定向辐射技术具有重要的指导意义。对不同曲面障板的反射规律做了分析,利用Euler方程建立了水下等离子体声源的聚束声场分布模型,利用迎风型WENO格式对方程进行求解,仿真计算了两种曲面障板的聚束过程和声场分布特性,并设计相关实验对仿真结果进行了验证。研究结果表明,仿真结果与实验测量结果基本吻合,得到的声场云图较为直观地反映了声源的反射聚束特性,证明了仿真模型的有效性,为进一步研究水下定向辐射技术、提高声源级奠定了基础。     

水下双层加肋圆柱壳全空间收发分置散射声场快速预报方法

针对水下双层加肋圆柱壳的全空间收发分置散射声场求解,提出了一种快速预报方法。该方法将散射声场表示为声散射传递函数与声源密度函数的乘积,以目标表面网格信息、少量的仿真或测试多基地散射声压数据作为已知信息,借助数值积分、矩阵理论、最小二乘法对其他收发分置散射声场进行预报。分别以有限元仿真和试验测试的散射声压数据作为输入,对水下双层加肋圆柱壳的多基地散射声场进行了计算,并与完全采用有限元方法的计算结果进行了对比。结果表明:该方法在目标表面结构和部分散射声场数据已知条件下能对目标的全空间收发分置散射声场进行预报;已知散射声压数据量越多,计算频率越低,预报精度越高。     

瞬态波位移场计算方法在相控阵声场模拟中的实验验证

通过实验研究了声束偏转不同角度和声场半径变化时线性相控阵的声场特性,并与基于瞬态波位移解析解的相控阵声场计算方法得到的理论结果进行了对比。首先研究了基于半圆圆弧离面瞬态位移的声场计算方法,并进行了数值计算可靠性的优化。然后建立了一套相控阵瞬态波位移场的测量系统,以半圆形钢板为试件,纵波传感器为接收传感器,测量了偏转角度为-10°,-30°,-50°,-70°,及声场半径为50 mm,75 mm,100 mm和125 mm时线性超声相控阵的声束指向性;同时以基于瞬态波位移解析解的相控阵声场计算方法,计算了相同条件下该相控阵的声束指向性。最后将实验结果与理论结果进行了对比,取得了较好的一致。说明基于瞬态波位移的相控阵声场计算方法可以较好地计算线阵的声场,对于声束偏转较大角度和声场半径变化时也同样有效。      

逆向傅里叶衍射定理快速预报水下目标三维声散射指向分布

提出逆向运用傅里叶衍射定理预报水下弱散射目标三维声散射指向分布的快速计算方法。依据目标形状、周围介质的密度和声速构建三维声场图像模型,建立散射远场积分结果与图像频域幅值的关系式,提取频域中半径为水中波数k;的球型表面上的幅值,获得精细化的宽带、全方位散射声压指向特性。数值计算表明:将傅里叶衍射定理逆向运用于解决声学正问题,适用于分层的、不均匀的、非规则及多体弱散射目标散射声场的求解。通过插值提取频域样本获得远场声压的方法,避免了有限元法(3D-FEM)所必须的大规模的网格划分和迭代运算,可以有效地减少计算成本并拓展散射频率响应的带宽。在水池中完成两种具有不同声学参数和形状目标指向性测试实验,得到散射声压指向性幅度函数与理论预报相一致。      

Bessel驻波对椭球形液滴的声辐射转矩

在实际的声操控中,由于声辐射力、表面张力和重力的共同作用,液滴往往呈现出椭球的形状,在螺旋声场中会受到力矩的作用而发生转动。从声波的散射理论出发,根据部分波展开法求解得到了椭球形液滴在Bessel驻波场中的声散射系数,并给出了其受到的声辐射转矩的解析式。在此基础上,对椭球形不可压缩液滴和椭球形可压缩液滴分别进行数值计算。仿真结果表明,不可压缩液滴的声辐射转矩与声束半锥角的关系更密切,而可压缩液滴则更依赖于特定的频率;提升Bessel驻波场的阶数有利于增强声辐射转矩的峰值,但在中低频处较难对可压缩液滴产生明显的力矩。该研究结果预期对利用螺旋声场进行液滴的操控具有理论指导作用。      

Stochastic differential equation analysis for sound scattering by random internal waves in the ocean

The scattering of a weakly divergent narrow sound beam by random inhomogeneities of a fluctuating ocean is considered in the coupled-mode approximation. The random index of sound refraction is described using the Garrett-Munk internal wave spectrum. The problem is solved using the stochastic differential equations for the first-and second-order statistical moments of the acoustic field. The equations are formulated according to the cumulant expansion method. The existence of weakly divergent narrow sound beams in long-range sound propagation was one of the last discoveries of L.M. Brekhovskikh, to which he attached much importance. The concentration of sound into narrow beams away from the axis of the underwater sound channel was first observed experimentally and then explained by Brekhovskikh and his former students Goncharov, Kurtepov, and Petukhov. In the present paper, the scattered field intensity of a sound beam is calculated for different frequencies and source depths. Analytical expressions are obtained for the coefficients of the differential equation. The intermode energy transfer that accompanies the long-range propagation of a weakly divergent sound beam is analyzed. A comparison with the conventionally used Monte Carlo simulation in the parabolic equation approximation is performed.     

深海海底反射区声场的相干结构及其有源声呐估距

有源声呐在探测深海海底反射区的目标时,由于声线大掠射角弯曲且声速沿声线传播路径不断变化,造成了常规估距方法产生较大的误差。有效声速法是减小常规估距方法误差的有效途径,但由于需要预先计算空间每一位置点的"声线时延有效声速"对,复杂度高,实时性差。针对有效声速法的实时性问题,本文基于深海海底反射区声场的相干结构,提出一种改进的有效声速估距方法.首先指出深海声场能量沿声源出射角强弱相间变化及其引起的海底反射区离散声呐可探测区现象,并利用深海近水面声源的声线干涉效应解释了该现象的物理机理,建立了声呐可探测区与高能量声线的量化关系。在此基础上,计算声呐可探测区边界位置的"声线时延-有效声速"对,并线性拟合出可探测区所有位置点对应的值。经仿真验证,该方法与传统的有效声速法均可实现对常规估距方法估距误差的有效校正。虽然该方法估距精度较传统的有效声速法略有增大,但计算复杂度和计算时间显著减小,实时性好,具有良好的工程应用前景。      

超声吸收体物理参数对红外热成像测量高强度聚焦超声声场的影响

超声吸收体的物理参数对利用水听器和红外热成像技术的高强度聚焦超声(High-Intensity Focused Ultrasound,HIFU)声场测量结果具有重要影响。为了探索超声吸收体的物理参数(密度、声速、衰减系数、热扩散系数、定压比热容)对测量结果的影响规律,本文根据层状模型计算出相同声源功率输出时不同物理参数对应的超声吸收体内部声场和热场,利用有限差分法计算出超声吸收体表面在辐照过程中的温度变化;利用基于水听器和红外热成像技术的聚焦声场测量方法测量出焦域内不同位置上声场特征值(轴线声强和-6 dB声束宽度),与通过理论计算得到的声源在纯水中声场特征值进行比较,分析了不同物理参数对测量结果的影响。超声吸收体的声、热学参数中除了声速外,其它物理参数的变化引起声场特征值的测量造成最大相对差异率小于15%。因为声波传播速度的改变会导致超声吸收体内部热场分布变化,使测量结果与理论计算值有较大偏差,其中-6 dB声束宽度和轴线声强最大相对差异率为95.37%和69.97%。因此在选择超声吸收体的声、热学参数时应重点关注声波在吸收体内声速的影响。超声吸收体的声学参数与水的声学参数相近时,可以在焦域内获得较好的测量结果。      

Simulation of sound field in a tissue medium generated by a concave spherically annular transducer

The concave spherically annular transducer is regarded as a negative and a positive concave spherical transducer, and the spheroidal beam equation is used to simulate the linear and nonlinear sound field in a tissue medium generated by this transducer. It is found that the acoustic focus of the ring does not coincide with the acoustic focus of its central part. If the width of the ring increases, its acoustic focus will move toward the geometric focus and the amplitudes of nonlinear harmonics will increase obviously. If there are several coaxial rings placed on the concave spherical surface, more than one peak will appear along the axial direction for the fundamental, and high harmonics focus better. The distribution of sound field will change with the number and the excited signals of rings, so it maybe is a potential approach to treat locally big tumors.     

Experiments with narrow movable beams of pulsed sound waves at 35 GHz

Only the study of absorption and angular dependent scattering of GHz and THz sound can give us detailed information about interaction processes of acoustical phonons. Therefore, the necessary techniques of sound generation, detection and undisturbed transmission to the target must be improved or to a large extent must still be developed. Corresponding attempts at 35 GHz are described. Sound-beam topography allowed to study crystal defects and also showed us where transmission without scattering occurred. For the piezoelectric generation a reflex klystron was used feeding a planar Hertzian resonator. The needle-like sound beams appeared as a double beam, each showing a minimum diameter of 35 μm for longitudinal polarized sound at an acoustic wavelength of 165 nm in X-cut quartz and 320 nm in a-cut sapphire. To obtain the smallest possible beam diameter for high resolution work, model calculations are presented for an excited fundamental Gaussian sound beam in an isotropic medium.     

Acoustooptic Bragg diffraction without overemodulation with a phased-array transducer

Acoustooptic Bragg diffraction without overemodulation (i.e., with the extended high-efficiency region on the diffraction efficiency dependence vs. sound amplitude) is studied theoretically and experimentally. It is found that this effect appears in the case of a symmetrically nonuniform acoustic field and is due to the equality of additional opposite phase shifts of light beams passing through symmetric regions of the acoustic field. The situation is considered when an acoustic field is excited by a three-section phased-array transducer. The conditions are determined, in which reverse optical power transfer from the diffracted beam to transmitted beam (overemodulation) in the case of a high (close to 100%) diffraction efficiency is considerably suppressed. In the case of a phased array, the effect weakly depends on the frequency of sound and the size of transducer sections, which makes it possible to observe it in a wide range of acoustooptic interaction parameters.     

Acoustic beam scattering and excitation of sphere resonance: Bessel beam example

The exact partial wave series for the scattering by a sphere centered on an ideal Bessel beam was recently given by Marston ["Scattering of a Bessel beam by a sphere," J. Acoust. Soc. Am. 121, 753-758 (2007)]. That series is applied here to solid elastic spheres in water and to an empty spherical shell in water. The examples are selected to illustrate the effect of varying the beam's conical angle so as to modify the coupling to specific resonances in the response of each type of sphere considered. The backscattering may be reduced or increased depending on properties of the resonance and of the specular contribution. Changing the conical angle is equivalent to changing the beamwidth. Some applications of the Van de Hulst localization principle to the interpretation of the partial wave series and to the interpretation of the scattering dependence on the beam's conical angle are discussed. Some potential applications to the analysis of the scattering by spheres of more general axisymmetric beams are noted.     

Acousto-optic interaction of diverging Gaussian beams in anisotropic media

A strong acousto-optic interaction of bounded light and sound beams of a Gaussian shape is considered for different geometries permitting long-term interaction. The anisotropic spreading of an acoustic beam in the course of its propagation is taken into account. The spectra of light beams are described by a Gaussian-polynomial expression, and a set of differential equations is derived which allows one to describe various geometries of acousto-optic diffraction. It is demonstrated that the transmission function of an acousto-optic cell essentially depends on the ratio between the dimensions of the light and sound beams and on the angle between the propagation direction of the diffracted beam and that of the incident beam. The degree of spreading of the acoustic beam noticeably affects the suppression of the transmission side lobes.     

Acoustic beam interaction with a rigid sphere: The case of a first-order non-diffracting Bessel trigonometric beam

Mathematical expressions for the acoustic scattering, instantaneous (linear), and time-averaged (nonlinear) forces resulting from the interaction of a new type of Bessel beam, termed here a first-order non-diffracting Bessel trigonometric beam (FOBTB) with a sphere, are derived. The beam is termed “trigonometric” because of the dependence of its phase on the cosine function. The FOBTB is regarded as a superposition of two equi-amplitude first-order Bessel vortex (helicoidal) beams having a unit positive and negative order (known also as topological charge), respectively. The FOBTB is non-diffracting, possesses an axial null, a geometric phase, and has an azimuthal phase that depends on cos(?±?₀), where ?₀ is an initial arbitrary phase angle. Beam rotation around its wave propagation axis can be achieved by varying ?₀. The 3D directivity patterns are computed, and the resulting modifications of the scattering are illustrated for a rigid sphere centered on the beam's axis and immersed in water. Moreover, the backward and forward acoustic scattering by a sphere vanish for all frequencies. The present paper will shed light on the novel scattering properties of an acoustical FOBTB by a sphere that may be useful in particle manipulation and entrapment, non-destructive/medical imaging, and may be extended to other potentially useful applications in optics and electromagnetism.