Optimized translation of microbubbles driven by acoustic fields

The problem of a single acoustically driven bubble translating unsteadily in a fluid is considered. The investigation of the translation equation identifies the inverse Reynolds number as a small perturbation parameter. The objective is to obtain a closed-form, leading order solution for the translation of the bubble, assuming nonlinear radial oscillations and a pressure field as the forcing term. In a second part, the periodic attractor of the Rayleigh-Plesset equation serves as basis for an optimal acoustic forcing designed to achieve maximized bubble translation over one dimensionless period. At near-resonant or super-resonant driving frequencies, it seems one cannot improve much on sinusoidal forcing. However at moderate acoustic intensity and sub-resonant frequencies, acoustic wave forms that enhance bubble collapse lead to displacement many times larger than the case of purely sinusoidal forcing. The survey covers a wide spectrum of driving ratios and bubble diameters including those relevant to biomedical applications. Shape stability issues are considered. Together, these results suggest new ways to predict some of the direct and indirect effects of the acoustic radiation force in applications such as targeted drug delivery, selective bubble driving, and accumulation.     

Nonlinear oscillation and acoustic scattering of bubbles

The scattered acoustic pressure and scattered cross section of bubbles is studied using the scattered theory of bubbles. The nonlinear oscillations of bubbles and the scattering acoustic fields of a spherical bubble cluster are numerically simulated based on the bubble dynamic and fluid dynamic. The influences of the interaction between bubbles on scattering acoustic field of bubbles are researched. The results of numerical simulation show that the oscillation phases of bubbles are delayed to a certain extent at different positions in the bubble cluster, but the radii of bubbles during oscillation do not differ too much at different positions. Furthermore, directivity of the acoustic scattering of bubbles is obvious. The scattered acoustic pressures of bubbles are different at the different positions inside and outside of the bubble cluster. The scattering acoustic fields of a spherical bubble cluster depend on the driving pressure amplitude, driving frequency, the equilibrium radii of bubbles, bubble number and the radius of the spherical bubble cluster. These theoretical predictions provide a further understanding of physics behind ultrasonic technique and should be useful for guiding ultrasonic application.     

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

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

声场中空化气泡的耦合振动及形状不稳定性的研究

计算了两个具有非球形扰动的气泡所组成系统的能量,并基于Lagrange方程得到了有声相互作用的非球形气泡的动力学方程和形状稳定性方程,研究了声场中非球形气泡间相互作用力对非球形气泡的形状不稳定性和气泡形状模态振幅的影响.研究结果表明声场中具有非球形扰动的气泡之间的耦合方式有两种:形状耦合模式和径向耦合模式,气泡之间的耦合方式取决于气泡形状扰动模态.由形状耦合及径向耦合产生的气泡之间的相互作用力能够改变单个气泡的形状不稳定及形状模态振幅,具体影响因素取决于声场驱动条件、气泡形状模态、相邻气泡的初始半径.     

A two-frequency acoustic technique for bubble resonant oscillation studies

A two-frequency acoustic apparatus has been developed to study the dynamics of a single gas or vapor bubble in water. An advantage of the apparatus is its capability of trapping a bubble by an ultrasonic standing wave while independently driving it into oscillations by a second lower frequency acoustic wave. For a preliminary application, the apparatus is used to study resonant oscillations. First, near-resonant coupling between the volume and the n = 3 shape oscillation modes of air bubbles at room temperature is studied, where n is the mode number. The stability boundary, amplitude versus frequency, of the volume oscillation forms a wedge centered at the resonant frequency, which qualitatively agrees with a theoretical prediction based on a phase-space analysis. Next, the resonant volume oscillations of vapor bubbles are studied. The resonant radius of vapor bubbles at 80 degrees C driven at 1682 Hz is determined to be 0.7 mm, in agreement with a prediction obtained by numerical simulation.     

含气泡液体中气泡振动的研究

研究了含气泡液体中单个气泡在驱动声场一定情况下的振动过程. 让每次驱动声场作用的时间特别短, 使气泡半径发生微小变化后再将其变化反馈到气泡群对驱动声场的散射作用中去, 从而可以得到某单个气泡周围受气泡散射影响后的声场, 接着再让气泡在该声场作用下做短时振动, 如此反复. 通过这样的方法, 研究了液体中单个气泡的振动情况并对其半径变化进行了数值模拟, 结果发现, 在液体中含有大量气泡的情况下, 某单个气泡的振动过程明显区别于液体中只有一个气泡的情况. 由于大量气泡和驱动声场的相互作用, 使气泡半径的变化存在多种不同的振动情况, 在不同的气泡大小和含量的情况下, 半径变化过程分别表现为: 在平衡位置附近振荡的过程; 周期性的空化过程; 一次空化过程后保持某一大小振荡的过程; 增长后维持某一大小振荡的过程等. 所以, 对于含气泡液体中气泡振动的研究, 在驱动声场一定的情况下, 必须考虑气泡含量的因素. 关键词： 含气泡液体 超声空化 散射 数值模拟     

声场与电场作用下空化泡的动力学特性

以液体为工作介质, 利用空化泡的RP控制方程, 模拟分析了无量纲化的电场频率、场强的幅值以及无量纲化的声波频率、 声压幅值的变化对空化泡运动特性的影响. 结果表明: 声场和电场联合作用时, 空化泡运动处于混沌区域范围远高于两者单独作用下空化泡的混沌区域范围. 这不仅对声空化的进一步研究具有重要的理论意义, 而且对于提高和改进空化降解有机污染物的技术也具有指导意义.     

两种气泡混合的声空化

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

Effect of an entrained air bubble on the acoustics of an ink channel

Piezo-driven inkjet systems are very sensitive to air entrapment. The entrapped air bubbles grow by rectified diffusion in the ink channel and finally result in nozzle failure. Experimental results on the dynamics of fully grown air bubbles are presented. It is found that the bubble counteracts the pressure buildup necessary for the droplet formation. The channel acoustics and the air bubble dynamics are modeled. For good agreement with the experimental data it is crucial to include the confined geometry into the model: The air bubble acts back on the acoustic field in the channel and thus on its own dynamics. This two-way coupling limits further bubble growth and thus determines the saturation size of the bubble.     

非球形效应对强声场中次Bjerknes力的影响

考虑了非球形气泡在声场中的形状振动,推导了非球形气泡和球形气泡之间的次Bjerknes力方程,数值模拟了声场中非球形气泡和球形气泡之间的次Bjerknes力和两个球形气泡之间的次Bjerknes力,并对非球形气泡和球形气泡之间的次Bjerknes力的影响因素进行了分析讨论.研究结果表明:当驱动声压振幅大于非球形气泡的Black阈值且又能使得非球形气泡稳定振动时,在第一个声驱动周期内,非球形气泡和球形气泡之间的次Bjerknes力和两个球形气泡的次Bjerknes力方向差异较大,在大小上是两个球形气泡次Bjerkens力的数倍,且有着更长的作用距离.非球形气泡和球形气泡之间的次Bjerknes力取决于非球形气泡的形状模态、两个气泡初始半径的比值、驱动声压振幅、气泡间距和两个气泡的相对位置.     

Numerical simulations of stable cavitation bubble generation and primary Bjerknes forces in a three-dimensional nonlinear phased array focused ultrasound field

We present a model developed for studying the generation of stable cavitation bubbles and their motion in a three-dimensional volume of liquid with axial symmetry under the effect of finite-amplitude phased array focused ultrasound. The density of bubbles per unit volume is determined by a nonlinear law which is a threshold-dependent function of the negative acoustic pressure reached in the liquid, in which nuclei are initially distributed. The nonlinear mutual interaction of ultrasound and bubble oscillations is modeled by a nonlinear coupled differential system formed by the wave and a Rayleigh-Plesset equations, for which both the pressure and the bubble oscillation variables are unknown. The system, which accounts for nonlinearity, dispersion, and attenuation due to the bubbles, is solved by numerical approximations. The nonlinear acoustic pressure field is then used to evaluate the primary Bjerknes force field and to predict the subsequent motion of bubbles. In order to illustrate the procedure, a medium-high and a low ultrasonic frequency configurations are assumed. Simulation results show where bubbles are generated, the nonlinear effects they have on ultrasound, and where they are relocated. Despite many physical restrictions and thanks to its particularities (two nonlinear coupled fields, bubble generation, bubble motion), the numerical model used in this work gives results that show qualitative coherence with data observed experimentally in the framework of stable cavitation and suggest their usefulness in some application contexts.     

Numerical investigation of acoustic vaporization threshold of microdroplets

A numerical model is presented for the acoustic vaporization threshold of a dodecafluoropentane (or perfluoropentane) microdroplet. The model is based on the Rayleigh-Plesset equation and is improved by properly treating the supercritical state that occurs when a bubble collapses rapidly and by employing the van der Waals equation of state to consider the supercritical state. The present computations demonstrate that the microdroplet vaporization behavior depends intricately on bubble compressibility, liquid inertia and phase-change heat transfer under acoustic excitation conditions. We present acoustic pressure-frequency diagrams for bubble growth regimes and the ADV threshold conditions. The effects of acoustic parameters, fluid properties and the droplet radius on the ADV threshold are investigated.     

Dynamics of a bubble cluster in an acoustic field

A mathematical model describing the dynamics of clustered gas bubbles under the effect of an acoustic field is presented. The proposed model is used as the basis for an analytical study of small bubble oscillations in monodisperse and polydisperse clusters and for a numerical study of nonlinear bubble oscillations under high-amplitude external pressures. The following effects are found to occur in a polydisperse cluster: a synchronization of the collapse phases of bubbles with different radii and a collapse intensification for bubbles of one size in the presence of bubbles of another size. These effects are explained by the interaction between bubbles of different radii in the cluster.     

The dependence of ultrasound contrast agents backscatter on acoustic pressure: theory versus experiment

Experimental investigations have not fully explored the interaction between ultrasound beams and microbubble contrast agents. Moreover theoretical investigations have not solved the problem of the microbubble oscillation. A simple in-vitro system based on a commercial scanner (ATL UM9) was used to insonate (3 MHz transmission) diluted contrast suspensions of Definity and Quantison at different acoustic pressures (0.27-1.52 MPa). The experimental data were referred to a blood mimicking fluid in order to extract an estimate of their scattering cross-section. The results were compared with the solutions of the three main bubble oscillatidn models, Rayleigh-Plesset, Herring and Gilmore. Non-linear solutions of the above models were produced numerically using the Mathematica Package Software. The experiments showed that both agents provided a linear increase in scattering cross-section with increasing acoustic pressure. The thick shelled Quantison provided an increasing number of scatterers with increasing acoustic pressure, which proved that free bubbles leaked out of the shell. At high acoustic pressures both Quantison and Definity scattering cross-sections were almost identical, and were probably that of a free bubble. The Rayleigh-Plesset model provided a scattering cross-section almost independent of acoustic pressure. On the contrary the scattering cross-sections calculated by the Herring and Gilmore models solutions displayed a definite dependence on acoustic pressure of an order higher than one, which is slightly higher than the order of dependence exhibited by the experimental data. However, the increase of the experimentally measured scattering cross-section with acoustic pressure was sharper than the calculated one by the above two models. This is most probably due to the fact that the models simulated damped and not free bubble oscillations. In conclusion the Rayleigh-Plesset model was inadequate in describing the bubble oscillations even at small diagnostic acoustic pressures. The Herring and Gilmore models could simulate the dependence of the scattering cross-section of encapsulated microbubbles on acoustic pressure. However the contribution of free bubble oscillations has still to be modelled.     

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

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

Ultrasonic scattering cross sections of shell-encapsulated gas bubbles immersed in a viscoelastic liquid: first and second harmonics

The oscillations of gas bubbles, without shell, immersed in viscoelastic liquids and driven by an acoustic wave have been the subject of several investigations. They demonstrate that the viscosity coefficient and the spring constant of the liquid have significant influence on the scattering cross section of the gas bubble. For shell-encapsulated gas bubbles, the investigations have been concentrated to bubbles immersed in a pure viscous liquid. This present work computes the ultrasonic scattering cross section, first and second harmonics, of shell-encapsulated gas bubbles immersed in a viscoelastic liquid. The theoretical model of the bubble oscillation is based on the generalized Rayleigh-Plesset equation of motion of a spherical cavity immersed in a viscoelastic liquid represented by a three-parameter linear Oldroyd model. The scattering cross section is computed for Albunex type of bubble (shell thickness=15 nm, shell shear viscosity=1.77 Pas, shell modulus of rigidity=88.8 MPa) irradiated by a 3.5 MHz ultrasonic pressure wave with an amplitude of 30 kPa. The results demonstrate that encapsulated bubbles respond independently of the surrounding liquid being pure viscous or viscoelastic as long as the surrounding liquid shear viscosity is as low as 10(-3) Pas. Nevertheless, for higher shear viscosities, the bubble responds differently if the surrounding liquid is pure viscous or viscoelastic. In general, the scattering cross sections of first and second harmonics are larger for the viscoelastic liquid.     

Acoustic vibrational resonance in a Rayleigh-Plesset bubble oscillator

The phenomenon of vibrational resonance (VR) has been investigated in a Rayleigh-Plesset oscillator for a gas bubble oscillating in an incompressible liquid while driven by a dual-frequency force consisting of high-frequency, amplitude-modulated, weak, acoustic waves. The complex equation of the Rayleigh-Plesset bubble oscillator model was expressed as the dynamics of a classical particle in a potential well of the Liénard type, thus allowing us to use both numerical and analytic approaches to investigate the occurrence of VR. We provide clear evidence that an acoustically-driven bubble oscillates in a time-dependent single or double-well potential whose properties are determined by the density of the liquid and its surface tension. We show both theoretically and numerically that, besides the VR effect facilitated by the variation of the parameters on which the high-frequency depends, amplitude modulation, the properties of the liquid in which the gas bubble oscillates contribute significantly to the occurrence of VR. In addition, we discuss the observation of multiple resonances and their origin for the double-well case, as well as their connection to the low frequency, weak, acoustic force field.     

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.     

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

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

Oscillating bubble concentration and its size distribution using acoustic emission spectra

New method has been proposed for the estimation of size and number density distribution of oscillating bubbles in a sonochemical reactor using acoustic emission spectra measurements. Bubble size distribution has been determined using Minnaert's equation [M. Minnaert, On musical air bubbles and sound of running water, Philanthr. Mag. 16 (1933) 235], i.e., size of oscillating bubble is inversely related to the frequency of its volume oscillations. Decomposition of the pressure signal measured by the hydrophone in frequency domain of FFT spectrum and then inverse FFT reconstruction of the signal at each frequency level has been carried out to get the information about each of the bubble/cavity oscillation event. The number mean radius of the bubble size is calculated to be in the range of 50-80mum and it was not found to vary much with the spatial distribution of acoustic field strength of the ultrasound processor used in the work. However, the number density of the oscillating bubbles and the nature of the distribution were found to vary in different horizontal planes away from the driving transducer surface in the ultrasonic bath. A separate set of experiments on erosion assessment studies were carried out using a thin aluminium foil, revealing a phenomena of active region of oscillating bubbles at antinodal points of the stationary waves, identical to the information provided by the acoustic emission spectra at the same location in the ultrasonic bath.