次Bjerknes力作用下气泡的体积振动和散射声场

利用Lagrange方程得到了次Bjerknes力作用下气泡的体积振动方程,并探讨了次Bjerknes力作用下不同参数对气泡体积振动振幅和振动初相位的影响,研究了振动初相位差为π和0的气泡对在液体中形成的散射声场特征.结果表明:次Bjerknes作用力下,相邻气泡半径、气泡间距、多方指数均能影响气泡的体积振动振幅,气泡对的均衡半径、气泡间距和驱动频率则对气泡振动初相位产生明显影响;相距很近、相位相差为π的两个气泡的散射声压与气泡体积振动振幅、气泡间距、驱动频率和振动初相位有关,随声场距离成反比减小,与声场位置有关,其平均散射声功率是单个孤立气泡的1/6(kd_(12))~2半径相同、相距很近、相位相同的两个6气泡的散射声压与气泡振动初相位、体积振动振幅、气泡间距、驱动频率有关,随声场距离成反比减小,其平均散射声功率是单个孤立气泡的4倍.     

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

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

含少量气泡流体饱和孔隙介质中的弹性波

考虑孔隙流体中含有少量气泡,且气泡在声波作用下线性振动,研究声波在这种孔隙介质中的传播特性.本文先由流体质量守恒方程和孔隙度微分与流体压力微分的关系推导出了含有气泡形式的渗流连续性方程;在处理渗流连续性方程中的气体体积分数时间导数时,应用Commander气泡线性振动理论导出气体体积分数时间导数与流体压强时间导数的关系,进而得到了修正的Biot形式的渗流连续性方程;最后结合Biot动力学方程求得了含气泡形式的位移场方程,便可得到两类纵波及一类横波的声学特性.通过对快、慢纵波的频散、衰减及两类波引起的流体位移与固体位移关系的考察,发现少量气泡的存在对快纵波和慢纵波的传播特性影响较大.     

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

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

声场作用下两空化泡相互作用的研究

建立了声场作用下两空化泡泡壁的运动方程,得出了双空化泡的共振频率,振动半径及空化噪声声压．由频率方程,振动半径和声压方程可以看出两气泡的运动情况与单气泡的运动情况有着明显的不同．共振频率,共振振幅及声压与两气泡之间的间距有关．在一定的简化条件下,运用MATLAB语言对共振频率,共振振幅及空化噪声声压进行了数值求解,发现共振频率和共振振幅随空泡间距的增大而增大,空化噪声声压随距离增大先增大后减小． 关键词： 超声 空化 频率 声压     

振动对水平管内液氢两相流影响的数值模拟

液氢加注过程中,管路的振动会影响液氢加注的安全性与效率。为了研究低温管路振动对液氢加注过程的影响,从而使得加注过程快速安全进行,建立了液氢管路的三维模型,采用流体体积函数(VOF)耦合水平集(Level-Set)方法,针对不同管路振动频率及振幅条件,模拟并分析了管内液氢两相流的流动及传热特性。研究结果表明,管路振动在频率和振幅较大时增强了液氢与管路之间的换热,使得管内气泡增多;管内的压降随着振幅频率的增加而增大,管路的振动破坏了管内液氢气液两相流的稳定边界,并且促进了气泡的分离与聚合,进而影响了低温输运管路的稳定性。通过数值模拟为液氢的加注提供了理论指导。     

非理想流体中Rayleigh-Taylor和Richtmyer-Meshkov不稳定性气泡速度研究

在随气泡顶端运动的坐标系中, 通过将理想流体模型推广到非理想流体的情况, 研究了流体黏性和表面张力对Rayleigh-Taylor (RT)和Richtmyer-Meshkov (RM)不稳定性气泡速度的影响. 首先得到了RT和RM不稳定性气泡运动的控制方程 (自洽的微分方程组); 其次给出了二维平面坐标和三维柱坐标中气泡速度的数值解和渐近解, 并定量分析了流体黏性和表面张力对RT和RM气泡速度和振幅的影响. 结果表明: 从线性阶段到非线性阶段的全过程中, 非理想流体中的气泡速度和振幅小于理想流体中的气泡速度和振幅. 也就是说, 流体黏性和表面张力对RT和RM不稳定性的发展都具有致稳作用. 关键词： Rayleigh-Taylor不稳定性 Richtmyer-Meshkov不稳定性 气泡速度 非理想流体     

超声场下刚性界面附近溃灭空化气泡的速度分析

为了揭示刚性界面附近气泡空化参数与微射流的相互关系, 从两气泡控制方程出发, 利用镜像原理, 建立了考虑刚性壁面作用的空化泡动力学模型. 数值对比了刚性界面与自由界面下气泡的运动特性, 并分析了气泡初始半径、气泡到固壁面的距离、声压幅值和超声频率对气泡溃灭的影响. 在此基础上, 建立了气泡溃灭速度和微射流的相互关系. 结果表明: 刚性界面对气泡振动主要起到抑制作用; 气泡溃灭的剧烈程度随气泡初始半径和超声频率的增加而降低, 随着气泡到固壁面距离的增加而增加; 声压幅值存在最优值, 固壁面附近的气泡在该最优值下气泡溃灭最为剧烈; 通过研究气泡溃灭速度和微射流的关系发现, 调节气泡溃灭速度可以达到间接控制微射流的目的.     

利用阻尼振动测量气垫导轨阻尼系数的新思路

利用气垫导轨上两根弹簧连接滑块做阻尼振动的模型,通过求解阻尼振动方程推算得到了滑块运动的最大速度衰减式,从而求得气垫导轨的阻尼系数.此外,从品质因数的角度推算阻尼系数的计算式,在阻尼很小时,该方法可视为最大速度衰减法的近似.通过气垫导轨阻尼振动实验,测量振幅、周期和光电门遮光时间等物理量,利用最大速度衰减和品质因数法得出阻尼系数的值非常接近.进一步对做阻尼振动的滑块拍摄视频,通过慢放视频得到其振动过程中振幅的衰减值,利用最大振幅衰减法计算出阻尼系数,所得结果与上述两种方法的数值符合较好,从另一个角度证明了本文提出的最大速度衰减法和品质因数法测阻尼系数的可靠性,该类方法在实验中易于操作、实践性强、所得结论可靠,为气垫导轨阻尼系数的测量提供了一种新的思路.     

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

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

Nonlinear effects in the dynamics of clouds of bubbles

This paper presents a spectral analysis of the response of a fluid containing bubbles to the motions of a wall oscillating normal to itself. First, a Fourier analysis of the Rayleigh-Plesset equation is used to obtain an approximate solution for the nonlinear effects in the oscillation of a single bubble in an infinite fluid. This is used in the approximate solution of the oscillating wall problem, and the resulting expressions are evaluated numerically in order to examine the nonlinear effects. Harmonic generation results from the nonlinearity. It is observed that the bubble natural frequency remains the dominant natural frequency in the volume oscillations of the bubbles near the wall. On the other hand, the pressure perturbations near the wall are dominated by the first and second harmonics present at twice the natural frequency while the pressure perturbation at the natural frequency of the bubble is inhibited. The response at the forcing frequency and its harmonics is explored along with the variation with amplitude of wall oscillation, void fraction, and viscous and surface tension effects. Splitting and cancellation of frequencies of maximum and minimum response due to enhanced nonlinear effects are also observed.     

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.     

Dynamics of double bubbles under the driving of burst ultrasound

This paper investigates the pulsations and translation of bubbles in a double-bubble system driven by burst ultrasound. Results illustrate that for two identical bubbles, decreasing the frequency of burst or increasing its amplitude can enhance the pulsations and improve the translation velocities of bubbles. In a certain scope, large bubble brings about fast translation velocity, but the velocity will fall down for too large bubble, such as the bubble with ambient radius over about its resonance radius. When the ambient radii of two bubbles are different, translation of the large bubble is smaller than that of the small bubble. In addition, the effect of initial distance between bubbles is described as well. If burst serials are used, shortening the time interval between each burst and improving the acoustic amplitude of bursts are beneficial for the translations of bubbles.     

Influence of rigid wall on the nonlinear pulsation of nearby bubble

This paper mainly focuses on the nonlinear pulsation of a bubble near the rigid wall. Dynamics of near-wall bubble and free bubble are discussed and compared in details. Investigation reveals as the driving acoustic pressure amplitude increases, nonlinear pulsation of bubble becomes intense gradually. Besides, decreasing the viscosity of host liquid is advantageous for the nonlinear pulsation of bubble. Bifurcation diagrams of bubble radius show acoustic reflection of the rigid wall makes the initial bifurcation appear at low driving acoustic amplitude and on bubble with small ambient radius, and makes the bifurcation still exist for bubble in high-viscosity liquids. That indicates the rigid wall will produce enhancement on the nonlinearity of nearby bubble. As the bubble approaches the wall, the enhancement becomes strong. Moreover, research on the influence of driving frequency shows the rigid wall makes the frequency band corresponding to chaos around the resonant frequency of free bubble shift downward.     

A new pressure formulation for gas-compressibility dampening in bubble dynamics models

We formulated a pressure equation for bubbles performing nonlinear radial oscillations under ultrasonic high pressure amplitudes. The proposed equation corrects the gas pressure at the gas–liquid interface on inertial bubbles. This pressure formulation, expressed in terms of gas-Mach number, accounts for dampening due to gas compressibility during the violent collapse of cavitation bubbles and during subsequent rebounds. We refer to this as inhomogeneous pressure, where the gas pressure at the gas–liquid interface can differ to the pressure at the centre of the bubble, in contrast to homogenous pressure formulations that consider that pressure inside the bubble is spatially uniform from the wall to the centre. The pressure correction was applied to two bubble dynamic models: the incompressible Rayleigh–Plesset equation and the compressible Keller and Miksis equation. This improved the predictions of the nonlinear radial motion of the bubble vs time obtained with both models. Those simulations were also compared with other bubble dynamics models that account for liquid and gas compressibility effects. It was found that our corrected models are in closer agreement with experimental data than alternative models. It was concluded that the Rayleigh–Plesset family of equations improve accuracy by using our proposed pressure correction.     

超声珩磨区实际气体的单空泡动力学分析

为进一步揭示功率超声振动的珩磨机理,以珩磨液为工作介质,研究了功率超声珩磨环境中实际气体的单空泡动力学特性。基于Rayleigh-Plesset方程,应用实际气体绝热方程和范德瓦尔斯方程对其进行了修正,建立了功率超声珩磨环境中实际气体的单空泡动力学方程以及实际气体单空泡共振频率方程。并运用4~5阶RungeKutta法模拟了不同超声条件(声压幅值、空泡初始半径、振动频率)对泡壁的运动以及运动速度的影响。结果表明:较高的声压幅值,空泡理论共振半径R'0与初始半径R0的比值为102数量级以及较低的超声频率有利于超声珩磨磨削区空化效应的发生。     

两种气泡混合的声空化

将非线性声波方程和改进的Rayleigh-Plesset方程联立可以描述空化环境中的声场及相应的气泡动力学特征. 用时域有限差分方法模拟了圆柱形容器内两种气泡相互混合时的空化情况. 在烧杯内的稳态背景声场形成过程中, 瓶壁耗散吸收扮演了重要的角色. 在稳态背景声场的基础上, 分析了混合气泡与声场的相互作用、气泡之间的相互作用、混合情况下的频谱特性. 结果表明: 两种气泡平衡半径都不太大时, 气泡与声场的相互作用不强, 声场及气泡的行为也比较规律; 相反, 当其中一种气泡平衡半径相对比较大时, 声场与气泡具有较强的非线性相互作用, 声场及气泡的行为表现出复杂的特性.