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

为了探讨含气泡液体对声波传播的影响, 研究了声波在含气泡液体中的线性传播. 在建立含气泡液体的声学模型时引入气泡含量的影响,建立气泡模型时引用 Keller的气泡振动模型并同时考虑气泡间的声相互作用,得到了经过修正的气泡振动方程. 通过对含气泡液体的声传播方程和气泡振动方程联立并线性化求解,在满足 (ω R₀)/c << 1 的前提下,得到了描述含气泡液体对声波传播的衰减系数和传播速度. 通过数值分析发现,在驱动声场频率一定的情况下,气泡含量的增加及气泡的变小均会导致衰减系数增加和声速减小;气泡的体积分数和大小一定时, 驱动声场频率在远小于气泡谐振频率的情况下,声速会随驱动频率的增加而减小; 气泡间的声相互作用对声波传播速度及含气泡液体衰减系数的影响不明显.最终认为气泡的大小、 数量和驱动声场频率是影响声波在含气泡液体中线性传播的主要因素. 关键词： 含气泡液体 线性声波 声衰减系数 声速     

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

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

含混合气泡液体中声波共振传播的抑制效应

当声波在含气泡的液体中传播时会出现共振传播现象,即在气泡的共振频率附近声衰减和声速会显著地增大,这是声空化领域的一个重要现象.以往的研究一般假设液体中只存在单一种类的气泡,因此忽略了声波共振传播的某些重要信息.本文研究了含混合气泡液体中声波的共振传播,混合气泡是指液体中包含多种静态半径不同的气泡.结果显示:在这种系统中存在声波共振传播的抑制效应,即与含单一种类气泡的系统相比,在含混合气泡的系统中声波的共振衰减和共振声速会明显变小.对于两种气泡混合、多种气泡混合以及气泡满足某种连续分布的系统,研究了抑制效应的本质和主要特征,此外还探究了黏性和空化率等对抑制效应的影响.本文的研究结果是对该领域现有知识的必要补充.     

气泡线性振动对含气泡水饱和多孔介质声传播的影响

为了研究孔隙水含少量气泡时多孔介质中波的传播,本文在Biot模型的基础上,将孔隙水中气泡的体积振动融合到多孔介质的孔隙流体渗流连续性方程中,从而得到了考虑气泡体积振动的孔隙流体渗流连续性方程.在此基础上,根据气泡线性振动下气泡瞬时半径和介质背景压力的关系,以及多孔介质运动方程和流体介质运动方程,导出了受气泡影响下多孔介质位移矢量波动方程,建立了非水饱和多孔介质声速频散和衰减预报模型.气泡的存在增大了孔隙水的压缩率,导致含气泡水饱和多孔介质声速的降低.当声波频率等于气泡的共振频率时,在声波激励下,介质呈现高频散,且孔隙水中的气泡产生共振,吸收截面达到最大,使得多孔介质的声衰减也达到最大.文中数值分析验证了上述结论,表明了气泡含量、大小和驱动声场频率是影响声波在含少量气泡的水饱和多孔介质中传播的主要因素.     

超声波作用下泡群的共振声响应

从球状泡群气泡动力学方程出发, 考虑泡群间次级声辐射的影响, 得到了声场中两泡群共同存在时气泡振动的动力学方程, 并以此为基础探讨声波驱动下双泡群振动系统的共振响应特征. 由于泡群间气泡间的相互作用, 系统存在低频共振和高频共振现象, 两不同共振频率的数值与泡群内气泡的本征频率相关. 泡群内气泡的本征频率又受到初始半径、泡群大小和泡群内气泡数量的影响. 气泡自由振动和驱动声波的耦合激起泡群内气泡的受迫振动, 气泡初始半径、气泡数密度和驱动声波频率等都会影响泡群内气泡的振动幅值和初相位. 关键词： 气泡群 共振 声响应 超声空化     

球状泡群内气泡的耦合振动

振动气泡形成辐射场影响其他气泡的运动, 故多气泡体系中气泡处于耦合振动状态. 本文在气泡群振动模型的基础上, 考虑气泡间耦合振动的影响, 得到了均匀球状泡群内振动气泡的动力学方程, 以此为基础分析了气泡的非线性声响应特征. 气泡间的耦合振动增加了系统对每个气泡的约束, 降低了气泡的自然共振频率, 增强了气泡的非线性声响应. 随着气泡数密度的增加, 振动气泡受到的抑制增强; 增加液体静压力同样可抑制泡群内气泡的振动, 且存在静压力敏感区(1–2 atm, 1 atm=1.01325×10⁵ Pa); 驱动声波对气泡振动影响很大, 随着声波频率的增加, 能够形成空化影响的气泡尺度范围变窄. 在同样的声条件、泡群尺寸以及气泡内外环境下, 初始半径小于5 μm 的气泡具有较强的声响应. 气泡耦合振动会削弱单个气泡的空化影响, 但可延长多气泡系统空化泡崩溃发生的时间间隔和增大作用范围, 整体空化效应增强.     

含气泡液体中声传播的解析解及其强非线性声特性

声波在含气泡的液体中传播时,气泡的受迫振动会引起强的声散射,并且由于振动的非线性,使得气泡产生的次级波不仅含有基频成分,而且还会有高次谐波。本文从理论上描述了气泡个数随尺并给出了含气泡液体的等效非线性声参数Ｂ／Ａ的计算公式理论与已有的实验观测符合较好,文中对含气泡水的声速和声衰减等特性也进行了讨论。     

低频超声空化场中柱状泡群内气泡的声响应

本文在气泡群振动模型的基础上,考虑气泡间耦合振动的影响,得到了均匀柱状泡群内振动气泡的动力学方程,以此为基础分析了低频超声空化场中柱形气泡聚集区内气泡的非线性声响应特征.气泡间的耦合振动增加了系统对每个气泡的约束,降低了气泡的自然频率,增强了气泡的非线性声响应.随着气泡数密度的增加,气泡的自然共振频率降低,受迫振动气泡受到的抑制增强.数值分析结果表明:1)驱动声波频率越低,气泡的初始半径越小,气泡数密度变化对气泡最大半径变化幅度的影响越大;2)气泡振动幅值响应存在不稳定区,不稳定区域分布与气泡初始半径、驱动声波压力幅值、驱动声波频率等因素有关.在低频超声波作用下,对初始半径处在1—10μm之间的空化气泡而言,气泡初始半径越小,气泡最大半径不稳定区分布范围越大,表明小气泡具有更强的非线性特征.因此,气泡初始半径越小,声环境变化对空化泡声响应稳定性影响越显著.     

声波在水-含气沉积物界面的反射

含气泡海洋沉积物的声学特性是海底探测的重要问题。为了研究气泡存在对水-含气沉积物界面声反射系数的影响,本文基于气泡振动修正的Biot波动方程推导了气泡存在修正的Biot弹性模量,并结合水-沉积物界面的“开孔”边界条件推导了声波从水入射到水-含气沉积物界面的反射系数。数值分析表明气泡的振动导致反射系数呈现显著的频率特性。在气泡共振频率附近,由于气泡的共振引发的强散射和强衰减,使得反射系数很大,无论以何种角度入射,声波都很难进入含气泡的沉积物。本文研究结果表明,气泡半径、含量、声波频率以及入射角度都是影响水-含气沉积物界面反射系数的主要因素。      

气泡体积分数对沙质沉积物低频声学特性的影响

由于有机物质分解等原因,实际的海底沉积物中存在气泡,气泡的存在会显著影响沉积物低频段的声学特性,因此研究气泡对沉积物低频段声速的影响机理具有重要意义.考虑到外场环境的不可控性,在室内水池中搭建了大尺度含气非饱和沙质沉积物声学特性获取平台,在有界空间中应用多水听器反演方法首次获取了含气非饱和沙质沉积物300—3000 Hz频段内的声速数据(79—142 m/s),并同时利用双水听器法获取了同一频段的数据(112—121 m/s).在声波频率远低于沉积物中最大气泡的共振频率时,根据等效介质理论,将孔隙水和气泡等效为一种均匀流体,改进了水饱和等效密度流体近似模型.模型揭示了气泡对沉积物低频段声学特性的影响规律,理论上解释了沉积物中声速的降低.通过分析模型预报声速对模型参数的敏感性,根据测量得到的声速分布反演得到了沉积物不同区域的气泡体积分数,气泡体积分数从1.07%变化到2.81%.改进的模型为沉积物中气泡体积分数估计提供了一种新方法.     

Weakly nonlinear focused ultrasound in viscoelastic media containing multiple bubbles

To facilitate practical medical applications such as cancer treatment utilizing focused ultrasound and bubbles, a mathematical model that can describe the soft viscoelasticity of human body, the nonlinear propagation of focused ultrasound, and the nonlinear oscillations of multiple bubbles is theoretically derived and numerically solved. The Zener viscoelastic model and Keller–Miksis bubble equation, which have been used for analyses of single or few bubbles in viscoelastic liquid, are used to model the liquid containing multiple bubbles. From the theoretical analysis based on the perturbation expansion with the multiple-scales method, the Khokhlov–Zabolotskaya–Kuznetsov (KZK) equation, which has been used as a mathematical model of weakly nonlinear propagation in single phase liquid, is extended to viscoelastic liquid containing multiple bubbles. The results show that liquid elasticity decreases the magnitudes of the nonlinearity, dissipation, and dispersion of ultrasound and increases the phase velocity of the ultrasound and linear natural frequency of the bubble oscillation. From the numerical calculation of resultant KZK equation, the spatial distribution of the liquid pressure fluctuation for the focused ultrasound is obtained for cases in which the liquid is water or liver tissue. In addition, frequency analysis is carried out using the fast Fourier transform, and the generation of higher harmonic components is compared for water and liver tissue. The elasticity suppresses the generation of higher harmonic components and promotes the remnant of the fundamental frequency components. This indicates that the elasticity of liquid suppresses shock wave formation in practical applications.     

Weakly nonlinear propagation of focused ultrasound in bubbly liquids with a thermal effect: Derivation of two cases of Khokolov–Zabolotskaya–Kuznetsoz equations

A physico-mathematical model composed of a single equation that consistently describes nonlinear focused ultrasound, bubble oscillations, and temperature fluctuations is theoretically proposed for microbubble-enhanced medical applications. The Khokhlov–Zabolotskaya–Kuznetsov (KZK) equation that has been widely used as a simplified model for nonlinear propagation of focused ultrasound in pure liquid is extended to that in liquid containing many spherical microbubbles, by applying the method of multiple scales to the volumetric averaged basic equations for bubbly liquids. As a result, for two-dimensional and three-dimensional cases, KZK equations composed of the linear combination of nonlinear, dissipation, dispersion, and focusing terms are derived. Especially, the dissipation term depends on three factors, i.e., interfacial liquid viscosity, liquid compressibility, and thermal conductivity of gas inside bubbles; the thermal conduction is evaluated by using four types of temperature gradient models. Finally, we numerically solve the derived KZK equation and show a moderate temperature rise appropriate to medical applications.     

A simple model of ultrasound propagation in a cavitating liquid. Part I: Theory, nonlinear attenuation and traveling wave generation

The bubbles involved in sonochemistry and other applications of cavitation oscillate inertially. A correct estimation of the wave attenuation in such bubbly media requires a realistic estimation of the power dissipated by the oscillation of each bubble, by thermal diffusion in the gas and viscous friction in the liquid. Both quantities and calculated numerically for a single inertial bubble driven at 20 kHz, and are found to be several orders of magnitude larger than the linear prediction. Viscous dissipation is found to be the predominant cause of energy loss for bubbles small enough. Then, the classical nonlinear Caflish equations describing the propagation of acoustic waves in a bubbly liquid are recast and simplified conveniently. The main harmonic part of the sound field is found to fulfill a nonlinear Helmholtz equation, where the imaginary part of the squared wave number is directly correlated with the energy lost by a single bubble. For low acoustic driving, linear theory is recovered, but for larger drivings, namely above the Blake threshold, the attenuation coefficient is found to be more than 3 orders of magnitude larger then the linear prediction. A huge attenuation of the wave is thus expected in regions where inertial bubbles are present, which is confirmed by numerical simulations of the nonlinear Helmholtz equation in a 1D standing wave configuration. The expected strong attenuation is not only observed but furthermore, the examination of the phase between the pressure field and its gradient clearly demonstrates that a traveling wave appears in the medium.     

Stationary soliton waves in a liquid with heat-conducting gas bubbles

Formation of an undamped soliton wave in the process of propagation of sound perturbations in a liquid with uniformly distributed gas bubbles is first revealed. The stationary soliton wave exists only in the presence of two pairs of competing factors, one of which is a balance between the nonlinearity of the medium and the linear wave dispersion, and another is a balance between the heat inflow into and outflow from the gas phase. The thermodynamic convertibility of dissipation of gas bubbles in the process of soliton wave propagation is revealed. The Korteweg-de Vries equation is derived for adiabatic bubbles. It differs radically from the Korteweg-de Vries equation obtained by Van Wijngaarden. __________ Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 6, pp. 84–88, June, 2006.     

丛聚的含气泡水对线性声传播的影响

含气泡水内气泡的空间分布会对线性声传播产生影响,导致实验结论与理论预测存在较大偏差.为解决这一问题,将准晶体近似引入到自洽方法中,导出了考虑空间分布时多分散含气泡水的等效声波波数.考虑到含气泡水内,气泡间存在小范围的聚集趋势(简称丛聚现象),在此基础上引入Neyman-Scott点过程描述了含气泡水内气泡的丛聚现象.分析发现,丛聚时,声速、声衰减的峰值将受到抑制,并向低频偏移,且抑制和频偏现象会随丛聚加剧而变强;随频率远离峰值段,丛聚对声传播的影响逐渐减弱.此外,考虑到空间分布的统计信息提取对相关研究的精确与否起到重要作用,引入了一种比例无偏估计,通过该方法获得了仿真环境下丛聚含气泡水模型的相速度及衰减系数,该建模及统计方法也可为相关实验工作提供理论基础.     

Numerical simulations of three-dimensional nonlinear acoustic waves in bubbly liquids

This paper presents three-dimensional simulations of nonlinear propagation of ultrasonic waves through bubbly liquids, which represent the continuity of our previous works included in the numerical tool SNOW-BL. The behavior of three-dimensional nonlinear acoustic waves in bubbly liquids is analyzed by means of numerical predictions. Nonlinearity, attenuation, and dispersion due to the presence of bubbles in the liquid are taken into account. The numerical solution to the differential problem is obtained by means of a finite-difference scheme. The simulations we present here consider a homogeneous distribution of bubbles in the liquid. Results compare high and low-amplitude waves to detect the nonlinear effects of the bubbles. Results are shown for radiation and enclosure problems.