Effects of medium viscoelasticity on bubble collapse strength of interacting polydisperse bubbles

Due to its physical and/or chemical effects, acoustic cavitation plays a crucial role in various emerging applications ranging from advanced materials to biomedicine. The cavitation bubbles usually undergo oscillatory dynamics and violent collapse within a viscoelastic medium, which are closely related to the cavitation-associated effects. However, the role of medium viscoelasticity on the cavitation dynamics has received little attention, especially for the bubble collapse strength during multi-bubble cavitation with the complex interactions between size polydisperse bubbles. In this study, modified Gilmore equations accounting for inter-bubble interactions were coupled with the Zener viscoelastic model to simulate the dynamics of multi-bubble cavitation in viscoelastic media. Results showed that the cavitation dynamics (e.g., acoustic resonant response, nonlinear oscillation behavior and bubble collapse strength) of differently-sized bubbles depend differently on the medium viscoelasticity and each bubble is affected by its neighboring bubbles to a different degree. More specifically, increasing medium viscosity drastically dampens the bubble dynamics and weakens the bubble collapse strength, while medium elasticity mainly affects the bubble resonance at which the bubble collapse strength is maximum. Differently-sized bubbles can achieve resonances and even subharmonic resonances at high driving acoustic pressures as the elasticity changes to certain values, and the resonance frequency of each bubble increases with the elasticity increasing. For the interactions between the size polydisperse bubbles, it indicated that the largest bubble generally has a dominant effect on the dynamics of smaller ones while in turn it is almost unaffected, exhibiting a pattern of destructive and constructive interactions. This study provides a valuable insight into the acoustic cavitation dynamics of multiple interacting polydisperse bubbles in viscoelastic media, which may offer a potential of controlling the medium viscoelasticity to appropriately manipulate the dynamics of multi-bubble cavitation for achieving proper cavitation effects according to the desired application.     

The acoustic emissions of cavitation bubbles in stretched vortices

Pairs of unequal strength, counter-rotating vortices were produced in order to examine the inception, dynamics, and acoustic emission of cavitation bubbles in rapidly stretching vortices. The acoustic signatures of these cavitation bubbles were characterized during their inception, growth, and collapse. Growing and collapsing bubbles often produced a sharp, broadband, pop sound. The spectrum of these bubbles, and the peak resonant frequency can generally be related to quiescent flow bubble dynamics and corresponding resonant frequencies. However, some elongated cavitation bubbles produced a short tonal burst, or chirp, with frequencies on the order of a few kilohertz. Theses frequencies are too low to be related to resonant frequencies of a bubble in a quiescent flow. Instead, the frequency content of the acoustic signal during bubble inception and growth is related to the volumetric oscillations of the bubble while it interacted with vortical flow that surrounds the bubble (i.e., the resonant frequency of the vortex-bubble system). A relationship was determined between the observed peak frequency of the oscillations, the highly stretched vortex properties, and the water nuclei content. It was found that different cavitation spectra could relate to different flow and fluid properties and therefore would not scale in the same manner.     

Theoretical estimation of sonochemical yield in bubble cluster in acoustic field

In order to learn more about the physical phenomena occurring in cloud cavitation, the nonlinear dynamics of a spherical cluster of cavitation bubbles and cavitation bubbles in cluster in an acoustic field excited by a square pressure wave are numerically investigated by considering viscosity, surface tension, and the weak compressibility of the liquid.The theoretical prediction of the yield of oxidants produced inside bubbles during the strong collapse stage of cavitation bubbles is also investigated. The effects of acoustic frequency, acoustic pressure amplitude, and the number of bubbles in cluster on bubble temperature and the quantity of oxidants produced inside bubbles are analyzed. The results show that the change of acoustic frequency, acoustic pressure amplitude, and the number of bubbles in cluster have an effect not only on temperature and the quantity of oxidants inside the bubble, but also on the degradation types of pollutants, which provides a guidance in improving the sonochemical degradation of organic pollutants.     

An Euler–Lagrange method considering bubble radial dynamics for modeling sonochemical reactors

Unsteady numerical computations are performed to investigate the flow field, wave propagation and the structure of bubbles in sonochemical reactors. The turbulent flow field is simulated using a two-equation Reynolds-Averaged Navier–Stokes (RANS) model. The distribution of the acoustic pressure is solved based on the Helmholtz equation using a finite volume method (FVM). The radial dynamics of a single bubble are considered by applying the Keller–Miksis equation to consider the compressibility of the liquid to the first order of acoustical Mach number. To investigate the structure of bubbles, a one-way coupling Euler–Lagrange approach is used to simulate the bulk medium and the bubbles as the dispersed phase. Drag, gravity, buoyancy, added mass, volume change and first Bjerknes forces are considered and their orders of magnitude are compared. To verify the implemented numerical algorithms, results for one- and two-dimensional simplified test cases are compared with analytical solutions. The results show good agreement with experimental results for the relationship between the acoustic pressure amplitude and the volume fraction of the bubbles. The two-dimensional axi-symmetric results are in good agreement with experimentally observed structure of bubbles close to sonotrode.     

Nonlinear dynamics and acoustic emissions of interacting cavitation bubbles in viscoelastic tissues

The cavitation-mediated bioeffects are primarily associated with the dynamic behaviors of bubbles in viscoelastic tissues, which involves complex interactions of cavitation bubbles with surrounding bubbles and tissues. The radial and translational motions, as well as the resultant acoustic emissions of two interacting cavitation bubbles in viscoelastic tissues were numerically investigated. Due to the bubble–bubble interactions, a remarkable suppression effect on the small bubble, whereas a slight enhancement effect on the large one were observed within the acoustic exposure parameters and the initial radii of the bubbles examined in this paper. Moreover, as the initial distance between bubbles increases, the strong suppression effect is reduced gradually and it could effectively enhance the nonlinear dynamics of bubbles, exactly as the bifurcation diagrams exhibit a similar mode of successive period doubling to chaos. Correspondingly, the resultant acoustic emissions present a progressive evolution of harmonics, subharmonics, ultraharmonics and broadband components in the frequency spectra. In addition, with the elasticity and/or viscosity of the surrounding medium increasing, both the nonlinear dynamics and translational motions of bubbles were reduced prominently. This study provides a comprehensive insight into the nonlinear behaviors and acoustic emissions of two interacting cavitation bubbles in viscoelastic media, it may contribute to optimizing and monitoring the cavitation-mediated biomedical applications.     

方波驱动下双气泡的动力学行为*

探索方波驱动下双气泡的脉动规律,能够促进方波在声空化工程中的实际应用。本文通过数值求解双气泡耦合方程组,研究了方波驱动下双气泡的动力学行为,得到了多种条件下不同时刻两个气泡半径的数值,并以此计算出气泡间的次Bjerknes力。研究表明,增大驱动频率会使得两个气泡膨胀时能达到的最大半径和次Bjerknes力减小。当两个气泡的平衡半径不同时,其中一个气泡的剧烈收缩会使得另一个气泡产生一个振动方向相反的声脉冲。随着两个气泡平衡半径差距的增加,气泡收缩的时间间隔增大。此外,当驱动声压幅值逐渐增大时,气泡脉动规律也会发生很大的变化。     

Sonoluminescence and dynamics of cavitation bubble populations in sulfuric acid

The detailed link of liquid phase sonochemical reactions and bubble dynamics is still not sufficiently known. To further clarify this issue, we image sonoluminescence and bubble oscillations, translations, and shapes in an acoustic cavitation setup at 23 kHz in sulfuric acid with dissolved sodium sulfate and xenon gas saturation. The colour of sonoluminescence varies in a way that emissions from excited non-volatile sodium atoms are prominently observed far from the acoustic horn emitter (“red region”), while such emissions are nearly absent close to the horn tip (“blue region”). High-speed images reveal the dynamics of distinct bubble populations that can partly be linked to the different emission regions. In particular, we see smaller strongly collapsing spherical bubbles within the blue region, while larger bubbles with a liquid jet during collapse dominate the red region. The jetting is induced by the fast bubble translation, which is a consequence of acoustic (Bjerknes) forces in the ultrasonic field. Numerical simulations with a spherical single bubble model reproduce quantitatively the volume oscillations and fast translation of the sodium emitting bubbles. Additionally, their intermittent stopping is explained by multistability in a hysteretic parameter range. The findings confirm the assumption that bubble deformations are responsible for pronounced sodium sonoluminescence. Notably the observed translation induced jetting appears to serve as efficient mixing mechanism of liquid into the heated gas phase of collapsing bubbles, thus potentially promoting liquid phase sonochemistry in general.     

Origin of the broad-band noise in acoustic cavitation

The broad-band noise has been experimentally used to monitor the cavitation activity in a sonochemical reactor, an ultrasonic cleaning bath, a biological tissue, etc. However, the origin of the broad-band noise is still under debate. In the present review, two models for the mechanism of the broad-band noise are discussed. One is acoustic emissions from chaotically (non-periodically) pulsating bubbles. The other is acoustic emissions from bubbles with temporal fluctuation in the number of bubbles. It is suggested that the latter mechanism is sometimes dominant. Further studies are required on the role for bubble cluster dynamics as well as the bubble–bubble interaction in the broad-band noise especially at relatively low ultrasonic frequencies.     

High-speed imaging of ultrasound driven cavitation bubbles in blind and through holes

The interest in application of ultrasonic cavitation for cleaning and surface treatment processes has increased greatly in the last decades. However, not much is known about the behavior of cavitation bubbles inside the microstructural features of the solid substrates. Here we report on an experimental study on dynamics of acoustically driven (38.5 kHz) cavitation bubbles inside the blind and through holes of PMMA plates by using high-speed imaging. Various diameters of blind (150, 200, 250 and 1000 µm) and through holes (200 and 1000 µm) were investigated. Gas bubbles are usually trapped in the holes during substrate immersion in the liquid thus preventing their complete wetting. We demonstrate that trapped gas can be successfully removed from the holes under ultrasound agitation. Besides the primary Bjerknes force and acoustic streaming, the shape oscillations of the trapped gas bubble seem to be a driving force for bubble removal out of the holes. We further discuss the bubble dynamics inside microholes for water and Cu²⁺ salt solution. It is found that the hole diameter and partly the type of liquid media influences the number, size and dynamics of the cavitation bubbles. The experiments also showed that a large amount of the liquid volume inside the holes can be displaced within one acoustic cycle by the expansion of the cavitation bubbles. This confirmed that ultrasound is a very effective tool to intensify liquid exchange processes, and it might significantly improve micro mixing in small structures. The investigation of the effect of ultrasound power on the bubble density distribution revealed the possibility to control the cavitation bubble distribution inside the microholes. At a high ultrasound power (31.5 W) we observed the highest bubble density at the hole entrances, while reducing the ultrasound power by a factor of ten shifted the bubble locations to the inner end of the blind holes or to the middle of the through holes.     

Fluid dynamics phenomena induced by power ultrasounds

Propagation of power ultrasound (from 20 to 800 kHz) through a liquid inside a cylindrical reactor initiates acoustic cavitation and also fluid dynamics phenomena such as free surface deformation, convection, acoustic streaming, etc. Mathematical modelling is performed as a new approach to predict where active bubbles are and how intense cavitation is. A calculation based on fluid dynamics equations is undertaken using computational fluid dynamics code; this is of great interest because such code provides not only the pressure field but also velocity and temperature fields. The link between the acoustic pressure and the cavitation field is clearly established. Moreover, the pressure profile near a free surface allows one to predict the shape of the acoustic fountain. The influence of the acoustic fountain on the wave propagation is shown to be important. The convective flow inside a reactor is numerically obtained and agrees well with particle image velocity measurements. Non-linearities arising from the dissipation of the acoustic wave are computed and lead to the calculation of the acoustic streaming. The superimposed velocity field (convective flow and acoustic streaming) succeeds in simulating the bubble behaviour at 500 kHz, for instance.     

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

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

Spiraling bubbles: how acoustic and hydrodynamic forces compete

Experiments to study the effect of acoustic forces on individual bubbles in shear flows have been carried out. In the system that we have used, the competition between acoustic and fluid dynamical forces results in a spiraling bubble trajectory. This dynamics is modeled by expressing the balance between Bjerknes and hydrodynamic forces in terms of an ordinary differential equation model, to which a separation of time scales is applied. The success of this model shows that the simple force-balance approach is still meaningful when bubbles are subjected to sound fields.     

Modeling bubble dynamics and radical kinetics in ultrasound induced microalgal cell disruption

Microalgal cell disruption induced by acoustic cavitation was simulated through solving the bubble dynamics in an acoustical field and their radial kinetics (chemical kinetics of radical species) occurring in the bubble during its oscillation, as well as calculating the bubble wall pressure at the collapse point. Modeling results indicated that increasing ultrasonic intensity led to a substantial increase in the number of bubbles formed during acoustic cavitation, however, the pressure generated when the bubbles collapsed decreased. Therefore, cumulative collapse pressure (CCP) of bubbles was used to quantify acoustic disruption of a freshwater alga, Scenedesmus dimorphus, and a marine alga, Nannochloropsis oculata and compare with experimental results. The strong correlations between CCP and the intracellular lipid fluorescence density, chlorophyll-a fluorescence density, and cell particle/debris concentration were found, which suggests that the developed models could accurately predict acoustic cell disruption, and can be utilized in the scale up and optimization of the process.     

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.     

Bubble interaction dynamics in Lagrangian and Hamiltonian mechanics

Two models of interacting bubble dynamics are presented, a coupled system of second-order differential equations based on Lagrangian mechanics, and a first-order system based on Hamiltonian mechanics. Both account for pulsation and translation of an arbitrary number of spherical bubbles. For large numbers of interacting bubbles, numerical solution of the Hamiltonian equations provides greater stability. The presence of external acoustic sources is taken into account explicitly in the derivation of both sets of equations. In addition to the acoustic pressure and its gradient, it is found that the particle velocity associated with external sources appears in the dynamical equations.     

A dual passive cavitation detector for localized detection of lithotripsy-induced cavitation in vitro

A passive cavitation detector (PCD) identifies cavitation events by sensing acoustic emissions generated by the collapse of bubbles. In this work, a dual passive cavitation detector (dual PCD), consisting of a pair of orthogonal confocal receivers, is described for use in shock wave lithotripsy. Cavitation events are detected by both receivers and can be localized to within 5 mm by the nature of the small intersecting volume of the focal areas of the two receivers in association with a coincidence detection algorithm. A calibration technique, based on the impulse response of the transducer, was employed to estimate radiated pressures at collapse near the bubble. Results are presented for the in vitro cavitation fields of both a clinical and a research electrohydraulic lithotripter. The measured lifetime of the primary growth-and-collapse of the cavitation bubbles increased from 180 to 420 microseconds as the power setting was increased from 12 to 24 kV. The measured lifetime compared well with calculations based on the Gilmore-Akulichev formulation for bubble dynamics. The radiated acoustic pressure 10 mm from the collapsing cavitation bubble was measured to vary from 4 to 16 MPa with increasing power setting; although the trends agreed with calculations, the predicted values were four times larger than measured values. The axial length of the cavitation field correlated well with the 6-dB region of the acoustic field. However, the width of the cavitation field (10 mm) was significantly narrower than the acoustic field (25 mm) as bubbles appeared to be drawn to the acoustic axis during the collapse. The dual PCD also detected signals from "rebounds," secondary and tertiary growth-and-collapse cycles. The measured rebound time did not agree with calculations from the single-bubble model. The rebounds could be fitted to a Rayleigh collapse model by considering the entire bubble cloud as an effective single bubble. The results from the dual PCD agreed well with images from high-speed photography. The results indicate that single-bubble theory is sufficient to model lithotripsy cavitation dynamics up to time of the main collapse, but that upon collapse bubble cloud dynamics becomes important.     

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

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

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.     

Experimental investigation of cavitational bubble dynamics under multi-frequency system

Acoustic emission spectra measurements have been carried out under mono and multi-frequency acoustic sources to understand the fundamental difference in bubble/cavity dynamics. The effect of introducing the dual and triple frequency acoustic waves of different frequency on the sono-chemical yield has also been investigated experimentally. The introduction of a second wave has increased the number of cavitating bubbles and as well as the collapsing intensity of cavities resulting into higher sono-chemical yield, and better effective utilization of reactor volume with a large number of resonating cavitating bubbles. To get the information about the intensity of each of the cavity oscillating events, decomposition of the pressure signal measured by the hydrophone in the frequency domain of the FFT power spectrum has been carried out. Inverse fourier reconstruction technique, has been used to elaborate the dynamics of the cavitating bubbles in the multi-frequency system.     

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.