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.     

Energy partitioning in laser-induced millimeter-sized spherical cavitation up to the fourth oscillation

To investigate the energy partitioning up to the fourth oscillation of a millimeter-scale spherical cavitation bubble induced by laser, we used nanosecond laser pulses to generate highly spherical cavitation bubbles and shadowgraphs to measure the radius-time curve. Using the extended Gilmore model and considering the continuous condensation of the vapor in the bubble, the time evolution of the bubble radius, bubble wall velocity, and pressure in the bubble is calculated till the 4th oscillation. Using Kirkwood-Bethe hypothesis, the evolution of velocity and pressure of shock wave at the optical breakdown, the first and second collapses are calculated. The shock wave energy at the breakdown and bubble collapse is directly calculated by numerical method. We found the simulated radius-time curve fits well with experimental data for the first four oscillations. The energy partition at the breakdown is the same as that in previous studies, the ratio of shock wave energy to bubble energy is about 2:1. In the first collapse and the second collapse, the ratio of shock wave energy to bubble energy is 14.54:1 and 2.81:1 respectively. In the third and fourth collapses, the ratio is less, namely than 1.5:1 and 0.42:1 respectively. The formation mechanism of the shock wave at the collapse is analyzed. The breakdown shock wave is mainly driven by the expansion of the supercritical liquid resulting from the thermalization of the energy of the free electrons in the plasma, and the collapse shock wave is mainly driven by the compressed liquid around the bubble.     

Shock-wave model of acoustic cavitation

Shock-wave model of liquid cavitation due to an acoustic wave was developed, showing how the primary energy of an acoustic radiator is absorbed in the cavitation region owing to the formation of spherical shock-waves inside each gas bubble. The model is based on the concept of a hypothetical spatial wave moving through the cavitation region. It permits using the classical system of Rankine-Hugoniot equations to calculate the total energy absorbed in the cavitation region. Additionally, the model makes it possible to explain some newly discovered properties of acoustic cavitation that occur at extremely high oscillatory velocities of the radiators, at which the mode of bubble oscillation changes and the bubble behavior approaches that of an empty Rayleigh cavity. Experimental verification of the proposed model was conducted using an acoustic calorimeter with a set of barbell horns. The maximum amplitude of the oscillatory velocity of the horns' radiating surfaces was 17 m/s. Static pressure in the calorimeter was varied in the range from 1 to 5 bars. The experimental data and the results of the calculations according to the proposed model were in good agreement. Simple algebraic expressions that follow from the model can be used for engineering calculations of the energy parameters of the ultrasonic radiators used in sonochemical reactors.     

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.     

Phase transitions of perfluorocarbon nanoemulsion induced with ultrasound: A mathematical model

While ultrasound has been used in many medical and industrial applications, only recently has research been done on phase transformations induced by ultrasound. This paper presents a numerical model and the predicted results of the phase transformation of a spherical nanosized droplet of perfluorocarbon in water. Such a model has applications in acoustic droplet vaporization, the generation of gas bubbles for medical imaging, therapeutic delivery and other biomedical applications.The formation of a gas phase and the subsequent bubble dynamics were studied as a function of acoustic parameters, such as frequency and amplitude, and of the physical aspects of the perfluorocarbon nanodroplets, such as chemical species, temperature, droplet size and interfacial energy. The model involves simultaneous applications of mass, energy and momentum balances to describe bubble formation and collapse, and was developed and solved numerically. It was found that, all other parameters being constant, the maximum bubble size and collapse velocity increases with increasing ultrasound amplitude, droplet size, vapor pressure and temperature. The bubble size and collapse velocity decreased with increasing surface tension and frequency. These results correlate with experimental observations of acoustic droplet vaporization.     

Defining the unknowns of sonoluminescence

As the intensity of a standing sound wave is increased the pulsations of a bubble of gas trapped at a velocity node attain sufficient amplitude so as to emit picosecond flashes of light with a broadband spectrum that increases into the ultraviolet. The acoustic resonator can be tuned so that the flashes of light occur with a clocklike regularity: one flash for each cycle of sound with a jitter in the time between flashes that is also measured in picoseconds. This phenomenon (sonoluminescence or “SL”) is remarkable because it is the only means of generating picosecond flashes of light that does not use a laser and the input acoustic energy density must be concentrated by twelve orders of magnitude in order to produce light. Light scattering measurements indicate that the bubble wall is collapsing at more than 4 times the ambient speed of sound in the gas just prior to the light emitting moment when the gas has been compressed to a density determined by its van der Waals hard core. Experiments indicate that the collapse is remarkably spherical, water is the best fluid for SL, some noble gas is essential for stable SL, and that the light intensity increases as the ambient temperature is lowered. In the extremely stable experimental configuration consisting of an air bubble in water, measurements indicate that the bubble chooses an ambient radius that is not explained by mass diffusion. Experiments have not yet been able to map out the complete spectrum because above 6 eV it is obscured by the cutoff imposed by water, and furthermore experiments have only determined an upper bound on the flash widths. In addition to the above puzzles, the theory for the light emitting mechanism is still open. The scenario of a supersonic bubble collapse launching an imploding shock wave which ionizes the bubble contents so as to cause it to emit Bremsstrahlung radiation is the best candidate theory but it has not been shown how to extract from it the richness of this phenomenon. Most exciting is the issue of whether SL is a classical effect or whether Planck's constant should be invoked to explain how energy which enters a medium at the macroscopic scale holds together and focuses so as to be emitted at the microscopic scale.     

Cavitation bubble collapse studied at 20 million frames per second

The collapse of laser-induced bubbles in water is investigated by high speed photography at framing rates as high as 20 million frames per second. The case of a spherical bubble in an unbounded liquid is compared with the Gilmore model. Bubbles collapsing in front of a solid wall show a rich dynamics depending on their normalized distance. Unprecedented details are given of the generic sequence of events leading to multiple shock waves and bubble shape metamorphosis upon collapse.     

Cavitation bubble behavior and bubble–shock wave interaction near a gelatin surface as a study of in vivo bubble dynamics

The collapse of a single cavitation bubble near a gelatin surface, and the interaction of an air bubble attached to a gelatin surface with a shock wave, were investigated. These events permitted the study of the behavior of in vivo cavitation bubbles and the subsequent tissue damage mechanism during intraocular surgery, intracorporeal and extracorporeal shock wave lithotripsy. Results were obtained with high-speed framing photography. The cavitation bubbles near the gelatin surface did not produce significant liquid jets directed at the surface, and tended to migrate away from it. The period of the motion of a cavitation bubble near the gelatin surface was longer than that of twice the Rayleigh's collapse time for a wide range of relative distance, L/R_max, excepting for very small L/R_max values (L was the stand-off distance between the gelatin surface and the laser focus position, and R_max was the maximum bubble radius). The interaction of an air bubble with a shock wave yielded a liquid jet inside the bubble, penetrating into the gelatin surface. The liquid jet had the potential to damage the gelatin. The results predicted that cavitation-bubble-induced tissue damage was closely related to the oscillatory bubble motion, the subsequent mechanical tissue displacement, and the liquid jet penetration generated by the interaction of the remaining gas bubbles with subsequent shock waves. The characteristic bubble motion and liquid jet formation depended on the tissue's mechanical properties, resulting in different damage mechanisms from those observed on hard materials.     

Laser-induced cavitation bubbles and shock waves in water near a concave surface

The interplay among the cavitation structures and the shock waves following a nanosecond laser breakdown in water in the vicinity of a concave surface was visualized with high-speed shadowgraphy and schlieren cinematography. Unlike the generation of the main cavitation bubble near a flat or a convex surface, the concave surface refocuses the emitted shock waves and causes secondary cavitation near the acoustic focus which is most pronounced when triggered by the shock wave released during the first main bubble collapse. The shock wave propagation, reflection from the concave surface and its scattering on the dominant cavity is clearly resolvable on the shadowgraphs. The schlieren approach revealed the pressure build up in the last stage of the collapse and the first stage of the rebound. A persistent low-density watermark is left behind the first collapse. The observed effects are important wherever cavities collapse near indented surfaces, such as in cavitation peening, cavitation erosion and ophthalmology.     

Interaction of shock waves with bubble-liquid drops

The regular reflection of an air shock wave from a spherical drop of a bubble liquid is studied. In the framework of an extended equilibrium model, the effect of the shock waves on single drops of various shapes and on drop ensembles (drop screens) is numerically investigated. It is shown that, when subjected to shock waves, bubble-liquid drops and drops of a bubble-free liquid collapse in a radically different way.     

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

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

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.     

Calculation of shock wave propagation in water containing reactive gas bubbles

The entry of a shock wave from air into water containing reactive gas (stoichiometric acetylene–oxygen mixture) bubbles uniformly distributed over the volume of the liquid has been numerically investigated using equations describing two-phase compressible viscous reactive flow. It has been demonstrated that a steady-state supersonic self-sustaining reaction front with rapid and complete fuel burnout in the leading shock wave can propagate in this bubbly medium. This reaction front can be treated as a detonation-like front or “bubble detonation.” The calculated and measured velocities of the bubble detonation wave have been compared at initial gas volume fraction of 2 to 6%. The observed and calculated data are in satisfactory qualitative and quantitative agreement. The structure of the bubble detonation wave has been numerically studied. In this wave, the gas volume fraction behind the leading front is approximately 3–4 times higher than in the pressure wave that propagates in water with air bubbles when the other initial conditions are the same. The bubble detonation wave can form after the penetration of the shock wave to a small depth (~300 mm) into the column of the bubbly medium. The model suggested here can be used to find optimum conditions for maximizing the efficiency of momentum transfer from the pressure wave to the bubbly medium in promising hydrojet pulse detonation engines.     

On the Formation of Centered Shock Waves in Gas Bubbles under the Conditions of Acoustic Cavitation

The possibility of the formation of centered shock waves in collapsing gas bubbles under the conditions of acoustic cavitation is considered. In this case, the overturning of the front of compression waves occurs at the instant the waves reach the center of the cavitation bubble, resulting in the highest possible temperatures and pressures inside the bubble. Examination of the magnetohydrodynamic equations has shown that the law of the motion of the wall of a bubble at the final stage of compression, described by the Rayleigh–Plesset equation, has a universal form and coincides with the condition of the formation of a spherically symmetric centered shock wave with the adiabatic constant = 5/3. For < 5/3, the collapse of a bubble occurs within a shorter time than it takes for a spherically symmetric centered shock wave to form. In this case, the overturning of the front of compression waves occurs earlier than they reach the center of the bubble, and shock waves are formed inside the bubble at different points. The most appropriate condition for the detection of centered shock waves is the cavitation in cryogenic fluids, such as helium, for which 5/3.     

A theoretical study of cavitation generated by an extracorporeal shock wave lithotripter

The intense acoustic wave generated at the focus of an extracorporeal shock wave lithotripter is modeled as the impulse response of a parallel RLC circuit. The shock wave consists of a zero rise time positive spike that falls to 0 at 1 microsecond followed by a negative pressure component 6 microseconds long with amplitudes scaled to +1000 and -160 bars, P+ and P-, respectively. This pressure wave drives the Gilmore-Akulichev formulation for bubble dynamics; the zero-order effect of gas diffusion on bubble response is included. The negative pressure component of a 1000-bar shock wave will cause a preexisting bubble in the 1- to 10-microns range to expand to over 100 times its initial size, R0, for 250 microseconds, with a peak radius of approximately 1400 microns, then collapse very violently, emitting far UV or soft x-ray photons (black body). Gas diffusion does not appreciably mitigate the amplitude of the pressure wave radiated at the primary collapse, but does significantly reduce the collapse temperature. Diffusion also increases the bubble radius from R0 up to 40 microns and extends the duration of ringing following the primary collapse, assuming that the bubble does not break up or shed microbubbles. Results are sensitive to P+/P- and to the duration of the negative pressure cycle but not to rise time.     

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.     

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.