Numerical investigation on acoustic cavitation characteristics of an air-vapor bubble: Effect of equation of state for interior gases

The cavitation dynamics of an air-vapor mixture bubble with ultrasonic excitation can be greatly affected by the equation of state (EOS) for the interior gases. To simulate the cavitation dynamics, the Gilmore-Akulichev equation was coupled with the Peng–Robinson (PR) EOS or the Van der Waals (vdW) EOS. In this study, the thermodynamic properties of air and water vapor predicted by the PR and vdW EOS were first compared, and the results showed that the PR EOS gives a more accurate estimation of the gases within the bubble due to the less deviation from the experimental values. Moreover, the acoustic cavitation characteristics predicted by the Gilmore-PR model were compared to the Gilmore-vdW model, including the bubble collapse strength, the temperature, pressure and number of water molecules within the bubble. The results indicated that a stronger bubble collapse was predicted by the Gilmore-PR model rather than the Gilmore-vdW model, with higher temperature and pressure, as well as more water molecules within the collapsing bubble. More importantly, it was found that the differences between both models increase at higher ultrasound amplitudes or lower ultrasound frequencies while decreasing as the initial bubble radius and the liquid parameters (e.g., surface tension, viscosity and temperature of the surrounding liquid) increase. This study might offer important insights into the effects of the EOS for interior gases on the cavitation bubble dynamics and the resultant acoustic cavitation-associated effects, contributing to further optimization of its applications in sonochemistry and biomedicine.     

Sonoluminescence

Abstract

Sonoluminescence is the light emission phenomenon from collapsing bubbles in liquid irradiated by an ultrasonic wave. In the present review, theoretical and experimental studies of the two types of sonoluminescence [single‐bubble sonoluminescence (SBSL) and multibubble sonoluminescence (MBSL)] are described. SBSL is a sonoluminescence from a single stably pulsating bubble trapped at the pressure antinode of a standing ultrasonic wave. MBSL is a sonoluminescence occurring from many bubbles in liquid irradiated by an ultrasonic wave. The theoretical and experimental studies suggest that SBSL originates in emissions from plasma inside the heated bubble at the bubble collapse, whereas MBSL originates both in emissions from plasma and in chemiluminescence inside heated bubbles at the bubble collapse. Unsolved problems of sonoluminescence have also been explained in detail.     

Theoretical prediction of the yield of strong oxides under acoustic cavitation

Considering liquid viscosity, surface tension, and liquid compressibility, the effects of dynamical behaviors of cavitation bubbles on temperature and the amount of oxides inside the bubble are numerically investigated by acoustic field,regarding water as a work medium. The effects of acoustic frequency, acoustic pressure amplitude, and driving waveforms on bubble temperature and the number of oxides inside the bubbles by rapid collapse of cavitation bubbles are analysed.The results show that the changes of acoustic frequency, acoustic pressure amplitude, and driving waveforms not only have an effect on temperature and the number of oxides inside the bubble, but also influence the degradation species of pollution,which provides guidance for improving the degradation of water pollution.     

Interpreting the influence of liquid temperature on cavitation collapse intensity through bubble dynamic analysis

The violent collapse of inertial bubbles generates high temperature inside and emits strong impulsive pressure. Previous tests on sonoluminescence and cavitation erosion showed that the influence of liquid temperature on these two parameters is different. In this paper, we conducted a bubble dynamic analysis to explore the mechanism of the temperature effect and account for the above difference. The results show that the increase of vapor at higher liquid temperatures changes both the external compression pressure and the internal cushion and is responsible for the variation of bubble collapse intensity. The different trends of the collapsing temperature and emitted sound pressure are caused by the energy distribution during the bubble collapse. Moreover, a series of simulations are conducted to establish the distribution map of the optimum liquid temperature where the collapse intensity is maximized. The relationship between the collapse intensity and the radial dynamics of the bubble is discussed and the reliable indicator is identified. This study provides a clear picture of how the thermodynamic process changes cavitation aggressiveness and enriches the understanding of this complex thermal-hydrodynamic phenomenon.     

Analysis of the effect of impact of near-wall acoustic bubble collapse micro-jet on Al 1060

The bubble collapse near a wall will generate strong micro-jet in a liquid environment under ultrasonic field. To explore the effect of the impact of near-wall acoustic bubble collapse micro-jet on an aluminum 1060 sheet, the cavitation threshold formula and micro-jet velocity formula were first proposed. Then the Johnson-Cook rate correlation material constitutive model was considered, and a three-dimensional fluid-solid coupling model of micro-jet impact on a wall was established and analyzed. Finally, to validate the model, ultrasonic cavitation test and inversion analysis based on the theory of spherical indentation test were conducted. The results show that cavitation occurs significantly in the liquid under ultrasonic field, as the applied ultrasonic pressure amplitude is much larger than liquid cavitation threshold. Micro pits appear on the material surface under the impact of micro-jet. Pit depth is determined by both micro-jet velocity and micro-jet diameter, and increases with their increase. Pit diameter is mainly related to the micro-jet diameter and d_p/d_j ≈ 0.95–1.2, while pit’s diameter-to-depth ratio is mainly negatively correlated with the micro-jet velocity. Wall pressure distribution is mostly symmetric and its maximum appears on the edge of micro-jet impingement. Obviously, the greater the micro-jet velocity is, the greater the wall pressure is. Micro pits formed after the impact of micro-jet on aluminum 1060 surface were assessed by ultrasonic cavitation test. Inversion analysis results indicate that equivalent stress, equivalent strain of the pit and impact strength, and velocity of the micro-jet are closely related with pit’s diameter-to-depth ratio. For the pit’s diameter-to-depth ratio of 16–68, the corresponding micro-jet velocity calculated is 310–370 m/s.     

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.     

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.     

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 estimation of the temperature and pressure within collapsing acoustical bubbles

Formation of highly reactive species such as OH, H, HO₂ and H₂O₂ due to transient collapse of cavitation bubbles is the primary mechanism of sonochemical reaction. The crucial parameters influencing the formation of radicals are the temperature and pressure achieved in the bubble during the strong collapse. Experimental determinations estimated a temperature of about 5000 K and pressure of several hundreds of MPa within the collapsing bubble. In this theoretical investigation, computer simulations of chemical reactions occurring in an O₂-bubble oscillating in water irradiated by an ultrasonic wave have been performed for diverse combinations of various parameters such as ultrasound frequency (20–1000 kHz), acoustic amplitude (up to 0.3 MPa), static pressure (0.03–0.3 MPa) and liquid temperature (283–333 K). The aim of this series of computations is to correlate the production of OH radicals to the temperature and pressure achieved in the bubble during the strong collapse. The employed model combines the dynamic of bubble collapse in acoustical field with the chemical kinetics of single bubble. The results of the numerical simulations revealed that the main oxidant created in an O₂ bubble is OH radical. The computer simulations clearly showed the existence of an optimum bubble temperature of about 5200 ± 200 K and pressure of about 250 ± 20 MPa. The predicted value of the bubble temperature for the production of OH radicals is in excellent agreement with that furnished by the experiments. The existence of an optimum bubble temperature and pressure in collapsing bubbles results from the competitions between the reactions of production and those of consumption of OH radicals at high temperatures.     

Proximity effect in Nb/Zr multilayers with variable Nb/Zr ratio

We have investigated the effect of varying the individual layer thickness on the superconducting transition temperature (T_C) of Nb/Zr multilayers. These thin film multilayer structures were deposited using UHV magnetron sputtering with layer thickness ranging from 0.5 to 8 nm. In conformity with the predictions of the de Gennes equations in the Cooper limit (layer thickness small compared to coherence length), we find that the T_C increases with increasing thickness of the Nb layer (when the Zr layer thickness is constant), and decreases with increasing thickness of the Zr layer (when the Nb layer thickness is constant). The possible effect of the existence of an interfacial Nb-Zr layer is discussed. We also point out the marked influence of the in-plane grain dimension on the T_C in these multilayers.     

A comprehensive numerical analysis of heat and mass transfer phenomenons during cavitation sono-process

The present study treats the effects of mass transport, heat transfer and chemical reactions heat on the bubble dynamics by spanning a range of ambient bubble radii. The thermodynamic behavior of the acoustic bubble was shown for three wave frequencies, 355, 515 and 1000 kHz. The used acoustic amplitude ranges from 1 to 3 atm. It has been demonstrated that the ambient bubble radius, R₀, of the maximal response (i.e., maximal bubble temperature and pressure, T_max and P_max) is shifted toward lower values if the acoustic amplitude (at fixed frequency) or the ultrasonic frequency (at fixed amplitude) are increased. The range of the ambient bubble radius narrows as the ultrasonic frequency increases. Heat exchange at the bubble interface was found to be the most important mechanism within the bubble internal energy balance for acoustic amplitudes lower than 2.5 and 3 atm for ultrasonic frequencies of 355 and 515 kHz, respectively. For acoustic amplitudes greater or equal to 2.5 and 3 atm, corresponding to 355 and 515 kHz, respectively, mass transport mechanism (i.e., evaporation and condensation of water vapor) becomes dominant compared to the other mechanisms. At 1000 kHz, the mechanism of heat transfer persists to be dominant for all the used acoustic amplitudes (from 1 to 3 atm). Practically, all the above observations were maintained for bubbles at and around the optimum bubble radius, whereas no significant impact of the three energetic mechanisms was observed for bubbles of too lower and too higher values of R₀ (limits of the investigated ranges of R₀).     

Monte Carlo study of surface roughening in the three-dimensional Ising model

Above a roughening temperature of about 0.56T _c the thickness of the two-dimensional interface separating domain on a simple cubic lattice is found to increase roughly logarithmically with system size. Also, the interface tension is determined for temperatures belowT _c.     

The influence of pressure on the acoustic cavitation in saturated CO2-expanded N,N-dimethylformamide

CO₂-expanded organic solvent is a kind of important fluid medium and has broad applications in chemical industry, environmental protection and other fields. Ultrasonic cavitation in gas expanded liquids (GXLs) is conducive to enhancing mass transfer and producing many exciting phenomena. In this paper, the ultrasonic cavitations and streaming in the saturated CO₂-expanded liquid N, N-dimethylformamide (DMF) at 4.2 MPa and 5.2 MPa are observed by a high-speed camera. The cavitation intensity and time trace of pressure pulses are recorded using a PZT hydrophone. The influences of gas–liquid equilibrium pressure and ultrasonic power on the cluster dynamics of transient and stable cavitation are examined. The excess molar enthalpies required for CO₂ dissociation from DMF are calculated by Peng-Robinson equations of state and the change of surface free energy of CO₂-expanded DMF is predicted. The results show that the excess enthalpy of the mixture is one of the key factors to control ultrasonic cavitation at high pressurized conditions, while the surface tension is the key factor for low pressure. As the increase of applied ultrasonic power, the formation and collapsing frequency of bubble clusters increases, and the amplitude and cyclic frequency of pressure pulse are enhanced. The transient cavitation intensity increases as it reaches a maximum value at a certain ultrasonic power and then decreases. The change trends of stable cavitation intensity under different pressures are basically same. It can be concluded from the evidence that ultrasonic cavitation in CO₂-expanded DMF is affected by the combined effect of compression and substitution: compression depresses the nucleation and growth of bubbles, while the high solubility of CO₂ in DMF is conducive to the generation of bubbles in cavitation.     

Surface and interface effects on structural transformation of vapor-deposited ethylbenzene films

We have investigated how the structures of vapor-deposited glassy films change with increasing temperature by using time-of-flight secondary ion mass spectrometry and ion scattering spectroscopy. It is found that intermixing of the topmost layer of an ethylbenzene film occur at temperature (~ 80 K) considerably lower than the glass transition temperature (T_g = 118 K) when the film is deposited at 20 K. This phenomenon can be interpreted as the occurrence of a two-dimensional liquid that diffuses into pores of the film, which is evidenced from comparison with surface diffusivity measurements using a porous silicon layer. For nonporous films deposited at higher temperatures, the molecules intermix gradually prior to the abrupt film morphology change at T_g. This phenomenon can be interpreted as decoupling between translational diffusivity and viscosity in the bulk. The film thickness has no significant effects on the evolution of supercooled liquid at T_g except for the monolayer film, whereas crystallization is quenched for the films thinner than 8 monolayers. The roles of the 2D liquid on the surface and an immobilized layer formed at the interface are discussed in finite-size effects on the glass-liquid transition and crystallization.     

HYDROELASTIC COUPLED VIBRATIONS IN A CYLINDRICAL CONTAINER WITH A MEMBRANE BOTTOM, CONTAINING LIQUID WITH SURFACE TENSION

A linear free hydroelastic vibration analysis of a frictionless liquid with a free surface contained in a cylindrical tank with a flexible bottom has been performed. The side-wall has been treated as rigid and the effect of surface tension taken into consideration. The container bottom was treated as a membrane, while for the free liquid surface the effect of two contact line conditions has been investigated. One edge condition was that of a slipping contact line, while the other one treated the contact line as fixed, ie., an anchored edge. The vibration characteristics of a membrane-liquid coupled system have been investigated for various system parameters, i.e., membrane tension parameter T, liquid surface tension parameter σ, material density parameter ρ, liquid height ratio ?₀ and vibration mode numbers m and n. The degree of coupling between a membrane and a liquid was represented with vibration mode diagrams as well as with frequency diagrams. For axisymmetric coupled vibrations with anchored edge condition, vibration mode exchanging of both a membrane and liquid free surface with membrane bottom tension parameter T has been investigated. An interesting phenomenon which is only observed for a flexible bottom container and an anchored edge free surface condition is presented.     

Shock wave emission upon spherical bubble collapse during cavitation-induced megasonic surface cleaning

When a gas bubble in a liquid interacts with an acoustic wave near a solid surface, the bubble first expands and then collapses. In this paper, a mathematical framework combining the Gilmore model and the method of characteristics is presented to model the shock wave emitted at the end of the bubble collapse. It allows to describe the liquid velocity at the shock front as a function of the radial distance to the bubble center in the case of spherical bubble collapse. Numerical calculations of the liquid velocity at the shock front have shown that this velocity increases with the acoustic amplitude and goes through a maximum as a function of the initial bubble radius. Calculations for different gas state equations inside the bubble show that the Van der Waals law predicts a slightly higher liquid velocity at the shock front than when considering a perfect gas law. Finally, decreasing the value of the surface tension at the bubble/liquid interface results in an increase of the liquid velocity at the shock front. Our calculations indicate that the strength of the shock waves emitted upon spherical bubble collapse can cause delamination of typical device structures used in microelectronics.     

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

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

On the thermodynamic behaviors and interactions between bubble pairs: A numerical approach

Thermodynamic behaviors and interactions between bubble pairs are important to better understand the cavitation phenomena. In this study, a compressible two-phase model, accounting for thermal effects to investigate the thermodynamic behaviors and interactions between bubble pairs, is developed in OpenFOAM. The volume of fluid (VOF) method is adopted to capture the interface. Validations are performed by comparing the simulation results of a single bubble and bubble pairs with corresponding experimental data. The dynamical behaviors of bubble pairs and their thermodynamic effect at different relative distances γ are investigated and discussed, which help reveal the bubble cloud dynamics. The quantitative analysis of γ effects on the maximum temperature during bubble collapse is performed with three distinct stages identified. For a single bubble collapsing near the rigid surface, the thermodynamic characteristics at different relative distances are similar to that of the bubble pairs, but the maximum temperature is higher since the single bubble can collapse to a smaller volume.     

Acoustic cavitation in 1-butyl-3-methylimidazolium bis(triflluoromethyl-sulfonyl)imide based ionic liquid

In this work, a comparison between the temperatures/pressures within acoustic cavitation bubble in an imidazolium-based room-temperature ionic liquid (RTIL), 1-butyl-3-methylimidazolium bis(triflluoromethyl-sulfonyl)imide ([BMIM][NTf₂]), and in water has been made for a wide range of cavitation parameters including frequency (140–1000 kHz), acoustic intensity (0.5–1 W cm⁻²), liquid temperature (20–50 °C) and external static pressure (0.7–1.5 atm). The used cavitation model takes into account the liquid compressibility as well as the surface tension and the viscosity of the medium. It was found that the bubble temperatures and pressures were always much higher in the ionic liquid compared to those predicted in water. The valuable effect of [BMIM][NTf₂] on the bubble temperature was more pronounced at higher acoustic intensity and liquid temperature and lower frequency and external static pressure. However, confrontation between the predicted and the experimental estimated temperatures in ionic liquids showed an opposite trend as the temperatures measured in some pure ionic liquids are of the same order as those observed in water. The injection of liquid droplets into cavitation bubbles, the pyrolysis of ionic liquids at the bubble-solution interface as well as the lower number of collapsing bubbles in the ionic liquid may be the responsible for the lower measured bubble temperatures in ionic liquids, as compared with water.