次声波在非均匀运动大气中非线性传播特性的研究

基于有限时域差分方法将大气中近似到二阶的非线性波动方程进行离散化,得到了数值模拟所采用的差分方程. 在此基础上,对线阵列辐射的脉冲声波在非均匀运动大气中的垂直和斜向传播进行了二维数值模拟,模拟了武汉地区(114:20°E, 30:37°N)在夏季和冬季UT=29000 s时开始传播的脉冲声波在不同时刻的声压分布. 模拟时通过采用Msise00和HWM93 两个大气模型,考虑了由于大气温度和密度变化以及大气风场存在所引起的大气不均匀性和运动性. 通过研究上述两季有风与无风条件下的声压差值p_r,可以发现:风场对次声波在传播中声压分布的影响较大;由于不同季节和不同传播距离上"有效声速"的不同,导致了两季p_r分布波形存在差异;风场对声波非线性传播的影响要远大于其对线性传播的影响. 关键词： 次声波传播 非均匀运动大气 有效声速     

含气泡液体中的非线性声传播

考虑到分布在液体中的气泡是声波在含气泡液体中传播时引起非线性的一个很重要的因素,本文研究了声波在含气泡液体中的非线性传播.将气体含量的影响引入到声波在液体中传播的方程中,从而得到声波在气液混合物中传播的数学模型.通过对该模型进行数值模拟发现,气体含量、驱动声场声压幅值及驱动声场作用时间均会影响到气液混合物中的声场分布及声压幅值大小.液体中的气泡会"阻滞"液体中声场的传播并将能量"聚集"在声源附近.对于连续大功率的驱动声场来说,液体中的气泡会"阻滞"气液混合物中声场及其能量的传播.     

非线性声场中的参量激发

本文基于求解单频声波方程近似解的方法,得到了非线性声场中谐波的声压与介质性质、初始声压幅值及频率之间的定量关系.并对两列相对声压幅值和相对频率不同情况下的声场分布进行了研究.通过分析单、双频声源辐射场中的谐波分布和传播规律发现:在非线性声场中会不断地出现新的谐波,激发的各阶谐波随着声波传播距离的增大逐渐增强而后减弱.在声源的附近,谐波的声压随基波声压振幅的增大而增大;但在基波的频率增大时反而会减小.在输入总声能相同的情况下,与单频声场相比双频声源辐射场的声能量分布较均匀,声的传播距离较大,远场中的谐波含量较大.结果表明,基波的频率越高,衰减得越快,谐波的积累越缓慢;声压的极值越大,基波声能量转移得越多,产生的谐波越多,基波的衰减越快,声压对远场声能的负效应增大;如果改用多频声源,并适当地控制输入声波的组成成分,可以达到改善声场分布均匀性、增大声辐射距离的效果.     

格子Boltzmann方法应用于气动声学研究

应用格子Boltzmann方法(LBM)对不同类型的气动声学问题进行数值研究.通过模拟一维平面声波和二维点源声波的传播,得到沿传播方向的声压脉动,其振荡幅值和衰减趋势与理论值相吻合.其次进行声波衍射和干涉现象的数值模拟.最后,模拟处于流场中运动声源辐射声场的多普勒效应.模拟结果说明LBM方法能较好地模拟低马赫数下的声学问题,包括声压脉动的传播,声波的波动特性以及流动与声波间的相互作用.     

超声在液体中的非线性传播及反常衰减

本文从广义的Navier-Stokes流体方程出发,考虑到流体介质的黏滞性和存在的热传导,导出了更接近实际流体的三维非线性声波动方程.鉴于声传播所涉及的空间和时间尺度的复杂性和多样性,文中针对一维情形下的非线性波动方程进行了求解和分析.由方程的二级近似解可以看出,声压振幅的衰减遵循几何级数规律,而且驱动声波的频率越高声压的衰减就越快.在满足条件ωb《ρ0c_0~2时,基波的衰减系数与驱动频率的平方及耗散系数的乘积成正比;二次谐波的衰减规律更加复杂,与频率的更高次幂相关.对声衰减系数及声压的分布进行数值计算发现,声压的分布还与初始的声压幅值及频率有关,初始的声压与频率越高衰减得越快.另外,当声压高于液体的空化阈值时,液体中就会出现大量的空化泡,文中模拟了单个空化泡的运动,发现随着声压的增大空化泡的振动越剧烈、空化泡所受的黏滞力变大,随着声波作用时间的增大黏滞力的幅值迅速增大并与驱动声压值同阶,因而空化泡的非线性径向运动引起的声衰减不容忽视.结果表明,驱动声压越高在空化区域附近引起的声衰减越快、输出的声压越低.     

海底爆破辐射声场的理论及数值研究

对固体基底浅海中声波传播的色散特性作了数值模拟与分析. 利用液体及固体中声传播理论导出了固体基底内爆破源在浅海中所激发声场的解析解, 在此基础上, 通过数值计算, 得出了浅海中不同空间点处的瞬态声压时域波形, 并在瞬态声压时域波形的平滑伪Wigner-Ville分布基础上, 对固体基底内爆破源在浅海中所激发瞬态声场的特性作了详细分析. 关键词： 海底爆破 声场 色散 瞬态声压     

除垢超声波传播影响因素的理论研究

从一维平面波理论入手分析了超声波声压分布特性。依据多普勒频移原理,在声场的运动方程,连续性方程,波动方程的基础上,建立一个超声波在流动的液体中传播的控制方程。根据轴对称模型的实际特点,简化了所得方程,并求出解析解。结果表明流动液体可以产生声波的衰减。液体的黏滞性是产生超声波衰减的重要原因。超声波的频率较高,液体的黏滞性对超声波衰减影响明显。依据黏滞力与速度梯度的关系,建立一个超声波在黏滞液体传播的控制方程,并依据边界条件求出解析解,反映了媒质黏滞性对超声波传播尤其是衰减特性的影响。     

次声波在平流层波导中传播的数值模拟

对次声波在大气中传播进行了建模。通过结合保色散关系空间差分格式和Runge-Kutta时间格式的数值方法,建立了次声波传播模型。应用该次声波模型,研究了在耗散的重力分层大气中次声波的平流层导行传播。数值模拟结果表明,当次声波波包在平流层高度上被反射时,反射区域存在焦散现象,在声波的声压下降的同时,声波的能量得到聚焦。通过数值模拟结果与射线计算结果的对比表明,大气中声波传播的轨迹的精确描述需要应用全波解。      

驻波法测定声速实验的matlab模拟

在驻波法测定声速实验中,根据发射器和接收器之间合成声场的位移波动方程和声压波动方程,用matlab程序模拟了声场和声压的振幅和初相位分布,并对模拟结果进行了分析和讨论。     

声波在固体板中的多径传播及其时间反转声场

采用无穷长条形阵元在固体板自由表面激发,研究了固体板中声波(纵波和横波)传播的多径效应,在射线近似下深入分析了柱面纵波和横波在固体板中的多次反射及波型转换特性,给出了简明解析表达式;将时间反转法应用于固体板中声波的传播,在理论和实验上分析了时间反转声场规律,证明了从不同阵元发出的沿不同路径传播的声波时间反转后同时同相到达原接收点,表明时间反转法能自动补偿固体板中由于多径效应造成的波形畸变;还通过聚焦增益和主副瓣比定量地分析了时间反转声场的自适应聚焦过程,考查了焦点位置和换能器阵列孔径对聚焦效果的影响,得到了理论和实验相符的结果。     

Nonlinear propagation and control of acoustic waves in phononic superlattices

The propagation of intense acoustic waves in a one-dimensional phononic crystal is studied. The medium consists in a structured fluid, formed by a periodic array of fluid layers with alternating linear acoustic properties and quadratic nonlinearity coefficient. The spacing between layers is of the order of the wavelength, therefore Bragg effects such as band gaps appear. We show that the interplay between strong dispersion and nonlinearity leads to new scenarios of wave propagation. The classical waveform distortion process typical of intense acoustic waves in homogeneous media can be strongly altered when nonlinearly generated harmonics lie inside or close to band gaps. This allows the possibility of engineer a medium in order to get a particular waveform. Examples of this include the design of media with effective (e.g., cubic) nonlinearities, or extremely linear media (where distortion can be canceled). The presented ideas open a way towards the control of acoustic wave propagation in nonlinear regime.     

Numerical simulation of the nonlinear ultrasonic pressure wave propagation in a cavitating bubbly liquid inside a sonochemical reactor

We investigate the acoustic wave propagation in bubbly liquid inside a pilot sonochemical reactor which aims to produce antibacterial medical textile fabrics by coating the textile with ZnO or CuO nanoparticles. Computational models on acoustic propagation are developed in order to aid the design procedures. The acoustic pressure wave propagation in the sonoreactor is simulated by solving the Helmholtz equation using a meshless numerical method. The paper implements both the state-of-the-art linear model and a nonlinear wave propagation model recently introduced by Louisnard (2012), and presents a novel iterative solution procedure for the nonlinear propagation model which can be implemented using any numerical method and/or programming tool. Comparative results regarding both the linear and the nonlinear wave propagation are shown. Effects of bubble size distribution and bubble volume fraction on the acoustic wave propagation are discussed in detail. The simulations demonstrate that the nonlinear model successfully captures the realistic spatial distribution of the cavitation zones and the associated acoustic pressure amplitudes.     

A derating method for therapeutic applications of high intensity focused ultrasound

Current methods of determining high intensity focused ultrasound (HIFU) fields in tissue rely on extrapolation of measurements in water assuming linear wave propagation both in water and in tissue. Neglecting nonlinear propagation effects in the derating process can result in significant errors. A new method based on scaling the source amplitude is introduced to estimate focal parameters of nonlinear HIFU fields in tissue. Focal values of acoustic field parameters in absorptive tissue are obtained from a numerical solution to a KZK-type equation and are compared to those simulated for propagation in water. Focal wave-forms, peak pressures, and intensities are calculated over a wide range of source outputs and linear focusing gains. Our modeling indicates, that for the high gain sources which are typically used in therapeutic medical applications, the focal field parameters derated with our method agree well with numerical simulation in tissue. The feasibility of the derating method is demonstrated experimentally in excised bovine liver tissue.     

Linear and nonlinear excitations in complex plasmas with nonadiabatic dust charge fluctuation and dust size distribution

Both linear and nonlinear excitation in dusty plasmas have been investigated including the nonadiabatic dust charge fluctuation and Gaussian size distribution dust particles. A linear dispersion relation and a Korteweg-de Vries-Burgers equation governing the dust acoustic shock waves are obtained. The relevance of the instability of wave and the wave evolution to the dust size distribution and nonadiabatic dust charge fluctuation is illustrated both analytically and numerically. The numerical results show that the Gaussian size distribution of dust particles and the nonadiabatic dust charge fluctuation have strong common influence on the propagation of both linear and nonlinear excitations.     

高频聚焦超声声场和温度场的仿真研究

为探究临床常用的7 MHz高频聚焦超声在多层生物组织中的声传播以及毫秒级时间内的生物传热规律问题,基于Westervelt方程和Pennes传热方程,使用有限元方法建立高频聚焦超声辐照多层组织的非线性热黏性声传播及传热模型。首先分析了线性模型和非线性模型之间的差异,然后在非线性模型下探究换能器的参数对声场和温度场的影响。仿真结果显示：在7 MHz频率下,当换能器输出声功率超过5 W时,声波传播的非线性效应不可忽视(p <0.05);当声功率从5 W增大到15 W时,非线性模型与线性模型预测的温度偏差从20%增加到34.703%;高频聚焦超声波的非线性行为比低频更加显著,基频能量向高次谐波转移的程度增大,声功率为10 W和15 W时4次谐波与基波之比分别达到7.33%和12.12%;高频换能器参数的改变对组织中声场和温度场分布的影响较大,换能器焦距从12 mm减小到11.2 mm,焦点处最高温度增加了77%。结果表明,7 MHz聚焦超声的非线性声传播需要考虑到4次谐波的影响。该文提出的多层组织非线性仿真模型可为高频聚焦超声换能器参数优化及制定安全、有效的术前治疗方案提供理论参考。     

非热离子对非均匀碰撞热尘埃等离子体中三维孤波的影响

从理论上研究了非热离子、外部磁场、碰撞对非均匀热尘埃等离子体中三维非线性尘埃声孤波的影响。运用约化摄动法得到描述三维非线性尘埃声孤波的非标准的变系数Korteweg-de Vries(KdV)方程。然后把非标准KdV方程变为标准的变系数KdV方程,并且得到了标准的变系数KdV方程的近似解析解。由此解析解可以看出,非热离子的数目、碰撞、非均匀性、波的斜向传播、尘埃颗粒和非热离子的温度对三维非线性尘埃声孤波的振幅和宽度有很大的影响。外部磁场对三维非线性尘埃声孤波的宽度有影响,而对其振幅没有影响。此外,波的相速度与非热离子、波的斜向传播、尘埃颗粒的温度和非均匀性有关。     

一维非线性声波传播特性

针对一维非线性声波的传播问题进行了有限元仿真和实验研究. 首先推导了一维非线性声波方程的有限元形式, 含有高阶矩阵的非线性项导致声波具有波形畸变、谐波滋生、基频信号能量向高次谐波传递等非线性特性. 编制有限元程序对一维非线性声波进行了计算并对仿真得到的畸变非线性声波信号进行处理, 分析其传播性质和物理意义. 为验证有限元计算结果, 开展了水中的非线性声波传播的实验研究, 得到了不同输入信号幅度激励下和不同传播距离的畸变非线性声波信号. 然后对基波和二次谐波的传播性质进行详细讨论, 分析了二次谐波幅度与传播距离和输入信号幅度的变化关系及其意义, 拟合出二次谐波幅度随传播距离变化的方程并阐述了拟合方程的物理意义. 结果表明, 数值仿真信号及其频谱均与实验结果有较好的一致性, 证实计算方法和结果的正确性, 并提出了具有一定物理意义的二次谐波随传播距离变化的简单数学关系. 最后还对固体中的非线性声波传播性质进行了初步探讨. 本研究工作可为流体介质中的非线性声传播问题提供理论和实验依据.     

Numerical simulation of nonlinear propagation of sound waves in a finite horn

Nonlinear acoustic propagation generated by a piston in a finite horn is numerically studied.A quasi-one-dimensional nonlinear model with varying cross-section uses high-order low-dispersion numerical schemes to solve the governing equation.Because of the nonlinear wave distortion and reflected sound waves at the mouth,broadband time-domain impedance boundary conditions are employed.The impedance approximation can be optimized to identify the complex-conjugate pole-residue pairs of the impedance functions,which can be calculated by fast and efficient recursive convolution.The numerical results agree very well with experimental data in the situations of weak nonlinear wave propagation in an exponential horn,it is shown that the model can describe the broadband characteristics caused by nonlinear distortion.Moreover,finite-amplitude acoustic propagation in types of horns is simulated,including hyperbolic,conical,exponential and sinusoidal horns.It is found that sound pressure levels at the horn mouth are strongly affected by the horn profiles,the driving velocity and frequency of the piston.The paper also discusses the influence of the horn geometry on nonlinear waveform distortion.     

Effect of acoustic nonlinearity on heating of biological tissue by high-intensity focused ultrasound

Effect of strong acoustic nonlinearity on the efficiency of heating of a biological tissue by high-intensity focused ultrasound in the modes of operation used in real clinical setups is studied. The spatial distributions of thermal sources and the corresponding temperature increments caused by ultrasonic absorption are analyzed. Numerical algorithms are developed for simulating the nonlinear focusing of ultrasound in the calculations of both the heat sources on the basis of the Khokhlov-Zabolotskaya-Kuznetsov-type equations and the temperature field in a tissue on the basis of an inhomogeneous thermal conduction equation with a relaxation term. It is demonstrated that in the mode of operation typical of acoustic surgery, the nonlinearity improves the locality of heating and leads to an increase in the power of thermal sources in the focus by approximately an order of magnitude. The diffusion phenomena in the tissue lead to a smoothing of the spatial temperature distributions, as compared to the distributions of thermal sources. In the case of one-second exposure in the nonlinear mode of focusing, the maximal temperature in the focus exceeds the values obtained in the approximation of linear wave propagation by a factor of three.