Numerical study of the effect of liquid compressibility on acoustic droplet vaporization

In acoustic droplet vaporization (ADV), a cavitated bubble grows and collapses depending on the pressure amplitude of the acoustic pulse. During the bubble collapse, the surrounding liquid is compressed to high pressure, and liquid compressibility can have a significant impact on bubble behavior and ADV threshold. In this work, a one-dimensional numerical model considering liquid compressibility is presented for ADV of a volatile microdroplet, extending our previous Rayleigh-Plesset based model [Ultrason. Chem. 71 (2021) 105361]. The numerical results for bubble motion and liquid energy change in ADV show that the liquid compressibility highly inhibits bubble growth during bubble collapse and rebound, especially under high acoustic frequency conditions. The liquid compressibility effect on the ADV threshold is quantified with varying acoustic frequencies and amplitudes.     

Erratum: transient pulsations of small gas bubbles in water [J. Acoust. Soc. Am. 84, 985-998 (1988)] [corrected and republished

Transient behavior of small gas bubbles in a liquid set into violent motion by ultrasonic pressure waves is of interest because of widespread use of microsecond pulses in diagnostic ultrasound. Such pulses contain only a few pressure cycles and the transient pulsations of bubbles set in motion by such pulses would determine the bubble-ultrasound interaction. A computer study has been made to obtain a global representation of the pulsation amplitudes R (t) of small gas bubbles (nuclei) in water during the first few cycles of a cw ultrasonic pressure. One objective was to obtain a better understanding of cavitation phenomena where many nuclei with initial radii Rn from 0.1-20 microns are set in motion at pressures ranging from 0.5-5 bars and at frequencies from 0.5-10 MHz. Results allowed construction of surfaces showing the relative bubble amplitude R/Rn as a function of Rn and of the time t/TA, where TA is the acoustic period. One finding is that, in the range of peak pressures found in diagnostic pulses, transient cavities would be generated during the first pressure cycle from nuclei with initial radii as small as a few microns (micron). Nuclei that grow into transient cavities in the first pressure cycle are here called "prompt" nuclei. At a specified pressure, the size range of radii Rn in which they occur decreases with increasing frequency. At 5 bars, the range of Rn for prompt nuclei is 0.166-11.35 microns at 0.5 MHz and vanishes at 10 MHz.     

Thermodynamic effect of single bubble near a rigid wall

The objective of this paper is to numerically investigate the thermodynamic effect during bubble collapse near a rigid boundary. A compressible fluid model is introduced to accurately capture the transient process of bubble shapes and temperature, as well as corresponding pressure, and velocity. The accuracy of the numerical model is verified by the experimental data of bubble shapes, and Keller-Kolodner equation as well as its thermodynamic equation. The results show that a bubble near the rigid boundary presents high-speed jet in collapse stage and counter jet in rebound stage, respectively. In the collapse stage, the bubble margin will shrink rapidly and do the positive work on the compressible vapor inside the bubble, then a significant amount of heat will be generated, and finally the generation of high-speed jet drives the low-temperature liquid outside the bubble to occupy the position of high-temperature vapor inside the bubble. In the rebound stage, the counter jet moving away from the rigid boundary takes part of heat away from the sub-bubble, which avoids the external work of the expansion of the sub-bubble and the temperature reduction caused by the dissipation effect of the vortex structure. In addition, the initial standoff has a significant effect on the thermodynamics of bubble oscillation. The temperature keeps increasing with the increase of the initial standoff in the collapse stage, while it shows a downward trend with the increase of the initial standoff in the rebound stage. That’s because the high-speed jet and counter jet of bubble gradually disappear when the initial standoff increases, which is the important reason for the opposite evolution trend of temperature in collapse and rebound stage.     

Occurrence of transient cavitation in pulsed sawtooth ultrasonic fields

Thus far, studies conducted to assess the safety of diagnostic ultrasound have employed sinusoidal sound fields. To evaluate the influence of nonlinearly distorted acoustic fields, this article compares the responses of microbubbles of variable size, exposed to (1) a sinusoidal pulse and (2) a sawtooth pulse. The nonlinear oscillations of a spherical bubble in a viscous compressible liquid stimulated into motion by an ultrasonic pulse are predicted, using a theoretical model for bubble dynamics. The maximum gas pressures inside the bubble when it collapses under the influence of a sinusoid or a sawtooth are deduced. Experimental work on Drosophila larvae exposed to sinusoidal and to sawtooth fields is consistent with the theoretical analysis.     

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

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

Experimental investigation on the effects of the standoff distance and the initial radius on the dynamics of a single bubble near a rigid wall in an ultrasonic field

Bubble behaviors near a boundary in an ultrasonic field are the fundamental forms of acoustic cavitation and of substantial importance in various applications, such as industry cleaning, chemical engineering and food processing. The effects of two important factors that strongly affect the dynamics of a single acoustic cavitation bubble, namely, the initial bubble radius and the standoff distance, were investigated in this work. The temporal evolution of the bubble was recorded using high speed microphotography. Meanwhile, the time of bubble collapse and the characteristics of the liquid jets were analyzed. The results demonstrate that the intensity of the acoustic cavitation, which is characterized by the time of bubble collapse and the liquid jet speed, reaches the optimum level under suitable values of the initial bubble radius and the normalized standoff distance. As the initial bubble radius and the normalized standoff distance increase or decrease from the optimal values, the time of the bubble collapse increases, and the first liquid jet’s speed decreases substantially, whereas the speeds of the second and third liquid jets exhibit no substantial changes. These results on bubble dynamics in an ultrasonic field are important for identifying or correcting the mechanisms of acoustic cavitation and for facilitating its optimization and application.     

Physical investigation of acoustic waves induced by the oscillation and collapse of the single bubble

The objective of this paper is to apply both experimental and numerical methods to investigate acoustic waves induced by the oscillation and collapse of a single bubble. In the experiments, the schlieren technique is used to capture the temporal evolution of the bubble shapes, and the corresponding acoustic waves. The results are presented for the single bubble generated by a low-voltage bubble generator in the free field of water. During the numerical simulations, a three-dimensional (3D) weakly compressible model is introduced to investigate the single bubble dynamics, including the generation and propagation of acoustic waves. The results show that (1) Compression wave, rarefaction wave and shock wave are generated during expansion stage, collapse stage and rebound stage of the bubble respectively. (2) Compression waves are induced by the rapid expansion of the bubble and eventually steepen into one shock wave propagating outward in the liquid, then another strong shock wave is emitted at the final collapse stage. The velocity and pressure of the liquid field increases after the shock wave. (3) Rarefaction waves are generated during the collapse stage due to the contraction of the bubble. The rarefaction wave reduces the liquid pressure and its spatial distribution is dispersive. The pressure of these acoustic waves and their effect on the liquid velocity attenuate with the increase of propagation distance.     

声驻波场中空化泡的动力学特性

以水为工作介质, 考虑了液体的可压缩性, 研究了驻波声场中空化泡的运动特性, 模拟了驻波场中各位置处空化泡的运动状态以及相关参数对各位置处空化泡在主Bjerknes力作用下运动方向的影响. 结果表明: 驻波声场中, 空化泡的运动状态分为三个区域, 即在声压波腹附近空化泡做稳态空化, 在偏离波腹处空化泡做瞬态空化, 在声压波节附近, 空化泡在主Bjerknes 力作用下, 一直向声压波节处移动, 显示不发生空化现象; 驻波场中声压幅值增加有利于空化的发生, 但声压幅值增加到一定上限时, 压力波腹区域将排斥空化泡, 并驱赶空化泡向压力波节移动, 不利于空化现象的发生; 当声频率小于初始空化泡的共振频率时, 声频率越高, 由于主Bjerknes 力的作用将有更多的空化泡向声压波节移动, 不利于空化的发生, 尤其是驻波场液面的高度不应是声波波长的1/4; 当声频率一定时, 空化泡初始半径越大越有利于空化现象的发生, 但当空化泡的初始半径超过声频率的共振半径时, 由于主Bjerknes力的作用将有更多的空化泡向声压波节移动, 不利于空化的发生.     

Collective bubble dynamics near a surface in a weak acoustic standing wave field

The transport of bubbles to a neighboring surface is very important in surface chemistry, bioengineering, and ultrasonic cleaning, etc. This paper proposes a multi-bubble transport method by using an acoustic standing wave field and establishes a model that explains the multi-bubble translation by expressing the balance between Bjerknes forces and hydrodynamic forces on a bubble in a liquid medium. Results indicated that the influence of primary Bjerknes force, secondary Bjerknes force, and buoyancy force on the bubble translation depends on the position of the target bubble in the acoustic field. Moreover, it was found that increasing the size of a bubble or pressure amplitude can accelerate the bubble motion and enhance the bubble-bubble interaction. The secondary Bjerknes force between two bubbles can switch from an attractive one when they oscillate in phase to a repulsive one when the bubble oscillations are out of phase. These findings provide an insight into the multi-bubble translation near a surface and can be applied to future bubble motion control studies, especially in drug delivery, sonoporation, and ultrasonic cleaning.     

超声振动珩磨作用下空化泡动力学及影响参数

为了合理利用超声振动珩磨作用下的空化效应,以磨削区单个空化泡为研究对象,考虑珩磨头合成扰动速度和珩磨压力的作用建立了磨削区空化泡的动力学模型。数值模拟了空化泡初始半径,珩磨压力,液体静压力和超声声压幅值对磨削区空化效应的影响。研究表明考虑超声振动珩磨作用时,空化泡膨胀的幅值会受到抑制,其溃灭时间也会缩短,而且较容易出现稳态空化。珩磨压力和液体静压力对磨削区空化主要起抑制作用,超声波声压幅值在一定范围内能够促进磨削区空化效果的提升。本文的研究为进一步理解超声振动珩磨的空化机理提供了理论支持。     

Implosion of an underwater spark-generated bubble and acoustic energy evaluation using the Rayleigh model

The growth, collapse, and rebound of a vapor bubble generated by an underwater spark is studied by means of high-speed cinematography, simultaneously acquiring the emitted acoustic signature. Video recordings show that the growth and collapse phases are nearly symmetrical during the first two or three cycles, the bubble shape being approximately spherical. After 2-3 cycles the bubble behavior changes from a collapsing/rebounding regime with sound-emitting implosions to a pulsating regime with no implosions. The motion of the bubble wall during the first collapses was found to be consistent with the Rayleigh model of a cavity in an incompressible liquid, with the inclusion of a vapor pressure term at constant temperature within each bubble cycle. An estimate of the pressure inside the bubble is obtained measuring the collapse time and maximum radius, and the amount of energy converted into acoustical energy upon each implosion is deduced. The resulting value of acoustic efficiency was found to be in agreement with measurements based on the emitted acoustic pulse.     

Optimum bubble temperature for the sonochemical production of oxidants

Numerical simulations of bubble oscillations in liquid water irradiated by an ultrasonic wave are performed for various acoustic amplitudes and various ambient pressures. In the numerical simulations, effect of non-equilibrium evaporation and condensation of water vapor at the bubble wall and that of chemical reactions of gases and vapor inside a bubble are taken into account. The oxidants such as OH radicals, O radicals, H(2)O(2) molecules, and O(3) molecules are created from water vapor inside a heated bubble when a bubble collapses strongly. They are dispersed into the liquid and solutes are oxidized by the oxidants, which is called sonochemical reactions. The computer simulations have revealed that there exists the optimum bubble temperature, which is about 5500 K, for the production of the oxidants inside an air bubble because at higher bubble temperature the oxidants are strongly consumed inside a bubble by oxidizing nitrogen. Correspondingly, there exists an optimum acoustic amplitude for the production of the oxidants, which is about 2.2 atm when the ultrasonic frequency is 140 kHz and the ambient pressure is 1 atm. For an oxygen bubble, on the other hand, the amount of the oxidants created inside a bubble becomes nearly independent of the bubble temperature at the collapse above about 6000 K because nitrogen is absent.     

Theoretical model of ice nucleation induced by inertial acoustic cavitation. Part 2: Number of ice nuclei generated by a single bubble

In the preceding paper (part 1), the pressure and temperature fields close to a bubble undergoing inertial acoustic cavitation were presented. It was shown that extremely high liquid water pressures but quite moderate temperatures were attained near the bubble wall just after the collapse providing the necessary conditions for ice nucleation. In this paper (part 2), the nucleation rate and the nuclei number generated by a single collapsing bubble were determined. The calculations were performed for different driving acoustic pressures, liquid ambient temperatures and bubble initial radius. An optimal acoustic pressure range and a nucleation temperature threshold as function of bubble radius were determined. The capability of moderate power ultrasound to trigger ice nucleation at low undercooling level and for a wide distribution of bubble sizes has thus been assessed on the theoretical ground.     

声场中双空化泡的运动特性

空化泡的运动特性是声场作用下的动力学行为,受空化泡初始半径,声压幅值,驱动声压频率,液体特性等众多因素的影响,是个复杂工程。本文从双空化泡运动方程出发,考虑到液体粘滞系数、空化泡辐射阻尼项的影响,研究了不同初始半径、驱动声压频率、驱动声压幅值、液体粘滞系数下空化泡泡壁的运动情况,研究结果表明不同初始半径、外界驱动声压频率、驱动声压幅值、液体粘滞系数均会对空化泡的膨胀比和空化泡的溃灭时间有一定影响。     

Acoustic cavitation mechanism: a nonlinear model

During acoustic cavitation process, bubbles appear when acoustic pressure reaches a threshold value in the liquid. The ultrasonic field is then submitted to the action of the bubbles. In this paper we develop a model to analyze the cavitation phenomenon in one-dimensional standing waves, based on the nonlinear code SNOW-BL. Bubbles are produced where the minimum rarefaction pressure peak exceeds the cavitation threshold. We show that cavitation bubbles appear at high amplitude and drastically affect (dissipation, dispersion, and nonlinearity) the ultrasonic field. This paper constitutes the first work that associates the nonlinear ultrasonic field to a bubble generation process.     

The Rayleigh-like collapse of a conical bubble

Key to the dynamics of the type of bubble collapse which is associated with such phenomena as sonoluminescence and the emission of strong rebound pressures into the liquid is the role of the liquid inertia. Following the initial formulation of the collapse of an empty spherical cavity, such collapses have been termed "Rayleigh-like." Today this type of cavitation is termed "inertial," reflecting the dominant role of the liquid inertia in the early stages of the collapse. While the inertia in models of spherical bubble collapses depends primarily on the liquid, experimental control of the liquid inertia has not readily been achievable without changing the liquid density and, consequently, changing other liquid properties. In this paper, novel experimental apparatus is described whereby the inertia at the early stages of the collapse of a conical bubble can easily be controlled. The collapse is capable of producing luminescence. The similarity between the collapses of spherical and conical bubbles is investigated analytically, and compared with experimental measurements of the gas pressures generated by the collapse, the bubble wall speeds, and the collapse times.     

Dependence of the characteristics of bubbles on types of sonochemical reactors

Computer simulations of bubble oscillations in liquid water irradiated by an ultrasonic wave have revealed that the characteristic of bubbles depends on types of sonochemical reactors: a horn-type reactor and a standing-wave type reactor. When the acoustic amplitude is large at 20 kHz, the bubble content is mostly water vapor even at the end of the bubble collapse and the temperature inside a bubble at the collapse is relatively low. On the other hand, when the acoustic amplitude is relatively low, the bubble content is mostly noncondensable gas at the end of the bubble collapse and the bubble temperature is relatively high. In a horn-type sonochemical reactor, the former type of bubbles are dominant because many bubbles exist near the horn-tip where the acoustic amplitude is large, while in a standing-wave type reactor the latter type of bubbles are dominant because the Bjerknes force gathers bubbles at a region where acoustic amplitude is relatively low.     

On the influence of stochastic pulsations of a bubble on its translational motion

This communication is devoted to theoretical analysis of the dynamics of a solitary cavitation bubble pulsating in a compressible viscous liquid under the action of a nonuniform acoustic field. The system of two nonlinear ordinary second-order differential equations is integrated numerically. In the range of acoustic field parameters corresponding to the principal resonance region, the bubble performs large-scale spatial oscillations. It is shown that in a very small range of initial radii, the bubble stops its oscillatory motion due to stochastic pulsations and is expelled into the region of the acoustic-pressure block. Therefore, stochastic pulsations of the bubble radically change the form of the solution to the system of the above-mentioned equations.     

声场中水力空化泡的动力学特性

以水为工作介质,考虑了液体黏性、表面张力、可压缩性及湍流作用等情况,对文丘里管反应器中空化泡在声场作用下的动力学行为特性进行了数值研究.分析了超声波频率、声压及喉径比对空化泡运动特性以及空化泡崩溃时所形成泡温以及压力脉冲的影响.结果表明,超声将水力空化泡运动调制成稳态空化,有利于增强空化效果. 关键词： 超声波 水力空化 湍流 气泡动力学     

Study on the bubble transport mechanism in an acoustic standing wave field

The use of bubbles in applications such as surface chemistry, drug delivery, and ultrasonic cleaning etc. has been enormously popular in the past two decades. It has been recognized that acoustically-driven bubbles can be used to disturb the flow field near a boundary in order to accelerate physical or chemical reactions on the surface. The interactions between bubbles and a surface have been studied experimentally and analytically. However, most of the investigations focused on violently oscillating bubbles (also known as cavitation bubble), less attention has been given to understand the interactions between moderately oscillating bubbles and a boundary. Moreover, cavitation bubbles were normally generated in situ by a high intensity laser beam, little experimental work has been carried out to study the translational trajectory of a moderately oscillating bubble in an acoustic field and subsequent interactions with the surface. This paper describes the design of an ultrasonic test cell and explores the mechanism of bubble manipulation within the test cell. The test cell consists of a transducer, a liquid medium and a glass backing plate. The acoustic field within the multi-layered stack was designed in such a way that it was effectively one dimensional. This was then successfully simulated by a one dimensional network model. The model can accurately predict the impedance of the test cell as well as the mode shape (distribution of particle velocity and stress/pressure field) within the whole assembly. The mode shape of the stack was designed so that bubbles can be pushed from their injection point onto a backing glass plate. Bubble radial oscillation was simulated by a modified Keller–Miksis equation and bubble translational motion was derived from an equation obtained by applying Newton’s second law to a bubble in a liquid medium. Results indicated that the bubble trajectory depends on the acoustic pressure amplitude and initial bubble size: an increase of pressure amplitude or a decrease of bubble size forces bubbles larger than their resonant size to arrive at the target plate at lower heights, while the trajectories of smaller bubbles are less influenced by these factors. The test cell is also suitable for testing the effects of drag force on the bubble motion and for studying the bubble behavior near a surface.