Three-dimensional numerical simulation of crown spike due to coupling effect between bubbles and free surface

The motion of gas bubbles beneath a free surface will lead to a spike of fluid on the free surface. The distance of the bubbles to the free surface is the key factor to different phenomena. When the inception distance varies in some range, crown phenomenon would happen after the impact of weak buoyancy bubbles, so this kind of spike is defined as crown spike in the present paper. Based on potential flow theory, a three-dimensional numerical model is established to simulate the motion of the free-surface spike generated by one bubble or a horizontal line of two in-phase bubbles. After the downward jet formed near the end of the collapse phase, the simulation of the free surface is performed to study the crown spike without regard to the toroidal bubble＇s effect. Calculations about the interaction between one bubble and free surface agree well with the experimental results conducted with a high-speed camera, and relative error is within 15%. Crown spike in both single- and two-bubble cases are simulated numerically. Different features and laws of the motion of crown spike, depending on the bubble-boundary distances and the inter-bubble distances, have been investigated.     

有倾角的竖直壁面附近气泡与自由面相互作用研究

本文采用边界元方法, 研究了具有倾角的竖直壁面附近气泡与自由面的相互耦合作用. 首先基于不可压缩势流理论, 建立了边界元气泡动力学模型, 并针对无限壁面附近的气泡和自由面作用问题, 采用镜像法模拟了倾斜壁面的作用. 然后, 基于本文所建立的数值模型, 分别计算了不同倾角的壁面对气泡和自由面水冢形态的影响, 发现倾斜角度会导致自由面的边界条件不连续, 从而对其运动产生明显的非线性影响. 最后, 分别分析了无浮力和有浮力情况下倾斜角度对其影响规律. 关键词： 水下爆炸 气泡动力学 倾斜壁面 自由面     

Effects of surface tension on the dynamics of a single micro bubble near a rigid wall in an ultrasonic field

Acoustic cavitation is a very important hydrodynamic phenomenon, and is often implicated in a myriad of industrial, medical, and daily living applications. In these applications, the effect mechanism of liquid surface tension on improving the efficiency of acoustic cavitation is a crucial concern for researchers. In this study, the effects of liquid surface tension on the dynamics of an ultrasonic driven bubble near a rigid wall, which could be the main mechanism of efficiency improvement in the applications of acoustic cavitation, were investigated at the microscale level. A synchronous high-speed microscopic imaging method was used to clearly record the temporary evolution of single acoustic cavitation bubble in the liquids with different surface tension. Meanwhile, the bubble dynamic characteristics, such as the position and time of bubble collapse, the size and stability of the bubbles, the speed of bubble boundaries and the micro-jets, were analyzed and compared. In the case of the single bubbles near a rigid wall, it was found that low surface tension reduces the stability of the bubbles in the liquid medium. Meanwhile, the bubbles collapse earlier and farther from the rigid wall in the liquids with lower surface tension. In addition, the surface tension has no significant influence on the speed of the first micro-jet, but it can substantially increase the speed of second and the third micro-jets after the first collapse of the bubble. These effects of liquid surface tension on the bubble dynamics can explain the mechanism of surfactants in numerous fields of acoustic cavitation for facilitating its optimization and application.     

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

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

Cavitation Bubble Collapse near a Curved Wall by the Multiple-Relaxation-Time Shan-Chen Lattice Boltzmann Model

The cavitation bubble collapse near a cell can cause damage to the cell wall. This effect has received inereasing attention in biomedical supersonics. Based on the lattice Boltzmann method, a multiple-relaxation-time Shan-Chen model is built to study the cavitation bubble collapse. Using this model, the cavitation phenomena induced by density perturbation are simulated to obtain the coexistence densities at certain temperature and to demonstrate the Young-Laplace equation. Then, the cavitation bubble collapse near a curved rigid wall and the consequent high-speed jet towards the wall are simulated. Moreover, the influences of initial pressure difference and bubble-wall distance on the cavitation bubble collapse are investigated.     

Experimental study of bubble dynamics in the neighbourhood of a vertical incomplete boundary

The bubbles have been widely used in biomedical field, military and chemical industry. The liquid jet generated by the bubble collapse through an orifice is utilized in needle-free injections and inkjet printing. In this paper we devised synchronized triggering equipment, experimentally investigated the mechanism in the interaction of an electric-spark generated a single bubble and a vertical wall with an air-back opening. Detailed observations were recorded and described for bubble oscillation, migration, jetting, as well as the high-speed water spike penetrating through the opening. The results revealed that there was a critical value of the bubble-wall distance, below which the bubble was directed away from the incomplete boundary, while the bubble may tear from the middle for larger distance. As the distance varied, we studied the volume of the water that rushed through the opening, the velocity at the tip of the water spike, and the center of the bubble as well as the migration of the bubble boundary. This work reveals that the high-speed water spike caused by the bubble may be a potential threat to the structures, specifically for cases with a small opening size and short bubble-boundary distance.     

Experimental studies of bubble dynamics inside a corner

The dynamics of a bubble near a corner formed by two flat rigid boundaries (walls), is studied experimentally using a spark-generated bubble. The expansion, collapse, rebound, re-collapse and migration of the bubble, along with jetting and protrusion, are captured using a high-speed camera. Our experimental observations reveal the behaviour of the bubble in terms of the corner angle and the dimensionless standoff distances to the near and far walls in terms of the maximum bubble radius. The bubble remains approximately spherical during expansion except for its surface becoming flattened when in close proximity to a wall. When a bubble is initiated at the bisector of the two walls, the bubble becomes oblate along the bisector during the late stages of collapse. A jet forms towards the end of collapse, pointing to the corner. The closer the bubble to the two walls, the more oblate along the bisector the bubble becomes, and the wider the jet. A bubble initiated near one of the two walls is mainly influenced by the nearer wall. The jet formed is pointing to the near wall but inclined towards the corner. After the jet penetrates through the bubble surface, the bubble becomes a bubble ring, and a bubble protrusion forms following the jet. The bubble ring collapses and subsequently disappears, while the protrusion firstly expands, and then collapses and migrates to the corner.     

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.     

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.     

Dynamics of a cavitation bubble near a solid wall

The cavitation bubble dynamics, the variation of pressure and velocity fields of the surrounding liquid in the process of the bubble axisymmetric compression near a planar solid wall are considered. It is assumed that the liquid is at rest at the initial moment of time, and the bubble has a spheroidal shape. The liquid is assumed inviscid and incompressible, its motion being potential. The bubble surface deformation and the liquid velocity on the surface are computed by the Euler scheme using the boundary element method until the moment of the collision of some parts of the bubble surface with one another. The influence of the distance of the bubble from the wall and its initial nonsphericity on the liquid pressure and velocity fields, the bubble shape, and the pressure inside the bubble at the end of the time interval under consideration are studied. The maximum pressure in liquid is shown to realize at the bottom of the cumulative jet arising at the bubble collapse with direction to the wall. In the upper part of this jet, the velocity and pressure are practically constant, and the pressure in the jet is approximately equal to the pressure in the bubble.     

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.     

曲面固壁附近空化泡溃灭动力学的 流体体积数值研究*

流体体积法(VOF)可以便捷、高效地实现对多相流界面的捕捉和追踪。本文基于VOF方法,对单个空化泡在曲面固壁附近的运动进行了数值模拟,从实验对比、压力场、速度场、温度场演化、溃灭时间、射流速度、固壁温度等方面分析了空化泡溃灭过程的热动力学影响。结果表明,数值模拟得到的空化泡形态演化与实验观测到的现象一致,随着位置参数、泡内外压差及曲面固壁尺寸的改变,空化泡热动力学行为也将发生变化,受到流体运动及射流冲击的影响,溃灭瞬间产生的高温高压使得曲面固壁温度升高。本文研究的曲面固壁附近空化泡溃灭效应,揭示了空化泡与曲面固壁间的相互作用规律,对学术研究及工程应用都具有重要意义。     

Nonlinear dynamics of a cavitation bubble pair near a rigid boundary in a standing ultrasonic wave field

The dynamics of a micrometer-sized bubble pair in water near a rigid boundary under standing ultrasonic wave excitation is investigated in this study. The viscous effect in the boundary layer at the air-water interface is considered following the viscous correction model. The evolution of the bubble surface at the collapsing stage of the bubble pair is presented for different parameter sets. The field pressure near the rigid boundary, which is induced by the oscillating bubble pair, and the high-speed water jet at the collapse stage, form the main focus of the analysis. This reveals that a horizontal configuration of the bubble pair retards the strength of the bubble jet towards the boundary, whilst a vertical configuration, especially with differently-sized bubbles, can enhance the bubble collapse. This study may help to understand the interaction of multiple bubbles in an acoustic field and its application to surface cleaning.     

上浮气泡在壁面处的弹跳特性研究

本文针对毫米量级的上浮气泡在壁面处的弹跳现象进行数值研究.基于势流方法求解气泡的运动,同时考虑气泡的表面张力作用.在伯努利方程中,对气泡与壁面之间水膜中因黏性引起的压力梯度进行修正,开发相应的计算程序,计算值与实验值符合良好.从气泡弹跳的基本现象入手,研究了特征参数对气泡弹跳过程的动态特性以及最终平衡形态的影响.发现随着泡在撞击壁面之前上浮距离增大,气泡回弹距离和弹跳周期增加,但是当上浮距离增加到一定程度后将不会影响气泡的弹跳特性;表面张力是影响气泡弹跳特性的重要因素,气泡的弹跳周期随其增大逐渐减小,但回弹距离却呈现先增后减的规律;最后,影响气泡最终平衡形态的主要因素是气泡的浮力参数与韦伯数.     

Numerical investigation of the deformation mechanism of a bubble or a drop rising or falling in another fluid

A numerical method for simulating the motion and deformation of an axisymmetric bubble or drop rising or falling in another infinite and initially stationary fluid is developed based on the volume of fluid （VOF） method in the frame of two incompressible and immiscible viscous fluids under the action of gravity, taking into consideration of surface tension effects. A comparison of the numerical results by this method with those by other works indicates the validity of the method. In the frame of inviseid and incompressible fluids without taking into consideration of surface tension effects, the mechanisms of the generation of the liquid jet and the transition from spherical shape to toroidal shape during the bubble or drop deformation, the increase of the ring diameter of the toroidal bubble or drop and the decrease of its cross-section area during its motion, and the effects of the density ratio of the two fluids on the deformation of the bubble or drop are analysed both theoretically and numerically.     

Investigation of noninertial cavitation produced by an ultrasonic horn

This paper reports on noninertial cavitation that occurs beyond the zone close to the horn tip to which the inertial cavitation is confined. The noninertial cavitation is characterized by collating the data from a range of measurements of bubbles trapped on a solid surface in this noninertial zone. Specifically, the electrochemical measurement of mass transfer to an electrode is compared with high-speed video of the bubble oscillation. This gas bubble is shown to be a "noninertial" event by electrochemical surface erosion measurements and "ring-down" experiments showing the activity and motion of the bubble as the sound excitation was terminated. These measurements enable characterization of the complex environment produced below an operating ultrasonic horn outside of the region where inertial collapse can be detected. The extent to which solid boundaries in the liquid cause the frequencies and shapes of oscillatory modes on the bubble wall to differ from their free field values is discussed.     

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.     

Modeling of surface cleaning by cavitation bubble dynamics and collapse

Surface cleaning using cavitation bubble dynamics is investigated numerically through modeling of bubble dynamics, dirt particle motion, and fluid material interaction. Three fluid dynamics models; a potential flow model, a viscous model, and a compressible model, are used to describe the flow field generated by the bubble all showing the strong effects bubble explosive growth and collapse have on a dirt particle and on a layer of material to remove. Bubble deformation and reentrant jet formation are seen to be responsible for generating concentrated pressures, shear, and lift forces on the dirt particle and high impulsive loads on a layer of material to remove. Bubble explosive growth is also an important mechanism for removal of dirt particles, since strong suction forces in addition to shear are generated around the explosively growing bubble and can exert strong forces lifting the particles from the surface to clean and sucking them toward the bubble. To model material failure and removal, a finite element structure code is used and enables simulation of full fluid–structure interaction and investigation of the effects of various parameters. High impulsive pressures are generated during bubble collapse due to the impact of the bubble reentrant jet on the material surface and the subsequent collapse of the resulting toroidal bubble. Pits and material removal develop on the material surface when the impulsive pressure is large enough to result in high equivalent stresses exceeding the material yield stress or its ultimate strain. Cleaning depends on parameters such as the relative size between the bubble at its maximum volume and the particle size, the bubble standoff distance from the particle and from the material wall, and the excitation pressure field driving the bubble dynamics. These effects are discussed in this contribution.     

Effect of surface tension on a liquid-jet produced by the collapse of a laser-induced bubble against a rigid boundary

The effect of surface tension on the behavior of a liquid-jet is investigated experimentally by means of a fiber-coupled optical beam deflection (OBD) technique. It is found that a target under water is impacted in turn by a laser-plasma ablation force and by a high-speed liquid-jet impulse induced by bubble collapse in the vicinity of a rigid boundary. The liquid-jet impact is found to be the main damage mechanism in cavitation erosion. Furthermore, the liquid-jet increases monotonously with surface tension, so cavitation erosion rises sharply with increasing surface tension. Surface tension also reduces bubble collapse duration. From the experimental results and the modified Rayleigh theory, the maximum bubble radius is obtained and it is found to reduce with increasing surface tension.     

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.