缓变主流中三维气泡的非线性振动

空化现象和水下噪声机制与液体中气泡的动力学行为密切相关．在无粘势流的假定下，采用多参数摄动分析，研究了缓变主流中三维气泡的非线性体积模态振动．推导了关于缓变泡形展开的各阶扰动方程，获得了一阶振动的演化方程和一些特殊情况下的解析解；并采用高阶有限元离散的边界积分方程方法，对平面固壁和自由面附近三维气泡的固有频率进行了数值计算     

On the modeling and simulation of a laser‐induced cavitation bubble

Single cavitation bubbles exhibit severe modeling and simulation difficulties. This is due to the small scales of time and space as well as due to the involvement of different phenomena in the dynamics of the bubble. For example, the compressibility, phase transition, and the existence of a noncondensable gas inside the bubble have strong effects on the dynamics of the bubble. Moreover, the collapse of the bubble involves the occurrence of critical conditions for the pressure and temperature. This adds extra difficulties to the choice of equations of state. Even though several models and simulations have been used to study the dynamics of the cavitation bubbles, many details are still not clearly accounted for. Here, we present a numerical investigation for the collapse and rebound of a laser‐induced cavitation bubble in liquid water. The compressibility of the liquid and vapor are involved. In addition, great focus is devoted to study the effects of phase transition and the existence of a noncondensable gas on the dynamics of the collapsing bubble. If the bubble contains vapor only, we use the six‐equation model for two‐phase flows that was modified in our previous work [A. Zein, M. Hantke, and G. Warnecke, J. Comput. Phys., 229(8):2964‐2998, 2010]. This model is an extension to the six‐equation model with a single velocity of Kapila et al. (Phys. Fluid, 13:3002‐3024, 2001) taking into account the heat and mass transfer. To study the effect of a noncondensable gas inside the bubble, we add a third phase to the original model. In this case, the phase transition is considered only at interfaces that separate the liquid and its vapor. The stiffened gas equations of state are used as closure relations. We use our own method to determine the parameters to obtain reasonable equations of state for a wide range of temperatures and make them suitable for the phase transition effects. We compare our results with experimental ones. Also our results confirm some expected physical phenomena. Copyright © 2013 John Wiley & Sons, Ltd.     

水平均流中细管排放气泡的三维数值模拟

在液体为无粘不可压，流动有势和气体遵循完全气体绝热关系的假定下，本文应用边界积分方程方法数值模拟了水平均流中垂直细管排放气泡的三维动力学问题，边界采用高阶有限元表达。文中介绍了有关泡面法向矢量、切向速度、曲率和接触线等的计算技术。与已知解的比较，表明了这一数值方法的高精度和优越性。算例显示了水平均流对于气泡形状和体积的影响     

Dynamics of large-scale eddies in turbulent flows

A model equation based on the equipartition of the turbulent dissipation is proposed for describing the dynamics of large-scale eddies in turbulent flows. The equation is reducible to the equation of motion of an inviscid fluid, so that the motion of the large-scale eddies can be described in terms of inviscid fluid dynamics. It is found that the large-scale eddies are always weakened by the background turbulence and their evolution is slowed down compared with the corresponding inviscid motion. In the case of turbulent mixing layer, its linear growth in downstream direction is accounted for by the exponential growth in time of the perturbation in an inviscid plane vortex sheet.     

水下爆炸气泡运动的理论研究

在适当深度的无黏、无旋的流体中对水下爆炸气泡运动特性进行理论研究。综合运用势流理论、能量方程以及拉格朗日方程建立气泡在不可压缩流体中的运动方程。并以此为基础,考虑重力、浮力以及阻力等多种因素对气泡运动特性的影响,通过引入新的边界积分方程,结合分析力学中完整非保守系统的Hamilton原理建立气泡在可压缩流体中的运动微分方程,并对微分方程进行求解。将方程的数值解与MSC.DYTRAN非线性有限元软件的计算结果以及经验公式进行对比,方程数值解与二者都具有较好的一致性。结果表明,基于非保守系统可压缩流体建立的气泡运动方程正确、可行,相关的理论研究和计算具有一定参考价值。     

Dynamics of a “collective” bubble in a magma melt flow behind the decompression wave front

Specific features of the dynamics of the wave field structure and growth of a “collective” bubble behind the decompression wave front in the “Lagrangian” section of the formed cavitation zone are numerically analyzed. Two cases are considered: with no diffusion of the dissolved gas from the melt to cavitation nuclei and with the diffusion flux providing an increase in the gas mass in the bubbles. In the first case, it is shown that an almost smooth decompression wave front approximately 100 m wide is formed, with minor perturbations that appear when the front of saturation of the cavitation zone with nuclei is passed. In the case of the diffusion process, the melt state behind the saturation front is principally different: jumps in mass velocity and viscosity are observed in the vicinity of the free surface, and the pressure in the “collective” cavitation bubble remains unchanged for a sufficiently long time interval, despite the bubble growth and intense diffusion of the gas from the melt. It is assumed that the diffusion process (and, therefore, viscosity) actually become factors determining the dynamics of growth of cavitation bubbles beginning from this time interval. A pressure jump is demonstrated to form near the free surface.     

nfluence of exciting parameters on the vibration of single cavitation bubble driven by intensive sound

根据考虑了液体可压缩性的改进的微气泡动力学方程,采用改进的初始半径对单泡超声空化现象进行了数值计算研究. 结果表明,微气泡振动对一些参量很敏感：微气泡振动半径与初始半径的比值随振动频率的增大而减小;提高声场声压会加剧气泡崩塌程度,但过高的声压又不能使微气泡崩塌;微气泡崩塌速率随气泡初始半径的增加而增大,在一定范围内能保证空化泡稳定振动,在初始半径为1.6\,$\mu$m 处空化程度最强,如果继续增大初始半径则空化程度减弱、甚至消失;微气泡崩塌程度随黏滞系数和表面张力的增大而减弱,过大的黏滞系数和表面张力会使微气泡崩塌难以发生. 计算结果与他人的实验数据相比,发现液体的可压缩性使单泡空化强度增强, 对最佳空化区域范围的确定有较大的影响.     

单空泡与自由液面相互作用规律研究进展

空化与空泡溃灭现象普遍存在于自然界、标识码械和生物医学等领域.空泡与自由面相互作用会产生瞬态强烈耦合,涉及到空泡非球形溃灭、自由面非线性变形及失标识码象,是流体力学领域重要的前沿与基础问题. 本文围绕这一热点,从空泡非球形演化和自由面变形规律角度出发,概标识码纳近年该领域的研究进展与成果. 对于近自由面空泡的非球形演化,基于表征开尔文冲量的无量纲参数,重点关注了体积振荡、射流生成、水锤效应及溃灭标识码生成等关键过程,介绍了关键参数的理论建模方法,获得了空泡溃灭过程中能量分配机制. 针对自由液面变形演化,根据细射流和粗射流生成和发展,归纳了 4 种典型现象及特点:透明水层及水柱生成、不稳定与稳标识码水裙结构. 进一步总结了开尔文冲量理论、界面凹陷奇点概念和泰勒不稳定性等理论模型的建立和应用,讨论了气泡溃灭过程、液面标识码界面稳定性等主要机制. 此外,本文也概述了空泡脉动对球状、圆柱状等非平面液面变形行为的影响,归纳了曲率对于液面变形的影响机制. 最后,针对目前研究状况提出该领域研究中尚未解决的问题,期望对将来的空泡及空泡群与自由液面相互作用深入研究提供借鉴.      

The effect of viscoelasticity on a rising gas bubble

This paper investigates the role of viscoelasticity on the dynamics of rising gas bubbles. The dynamics of bubbles rising in a viscoelastic liquid are characterised by three phenomena: the trailing edge cusp, negative wake, and the rise velocity jump discontinuity. There is much debate in the literature over the cause of the jump discontinuity, which is observed once the bubble exceeds a certain critical volume. In this paper, the employment of some choice modelling assumptions allows insights into the mechanisms of the jump discontinuity which cannot be ascertained experimentally. The ambient fluid is assumed incompressible and the flow irrotational, with viscoelastic effects included through the stress balance on the bubble surface. The governing equations are solved using the boundary element method. Some Newtonian predictions are discussed before investigating the role of viscoelasticity. The model predicts the trademark cusp at the trailing end of a rising bubble to a high resolution. However, the irrotational assumption precludes the prediction of the negative wake. The corresponding absence of the jump discontinuity supports the hypothesis that the negative wake is primarily responsible for the jump discontinuity, as mooted in previous studies.     

A novel numerical scheme for the investigation of surface tension effects on growth and collapse stages of cavitation bubbles

In the present study the effects of surface tension on the growth and collapse stages of cavitation bubbles are studied individually for both spherical and nonspherical bubbles. The Gilmore equation is used to simulate the spherical bubble dynamics by considering mass diffusion and heat transfer. For the collapse stage near a rigid boundary, the Navier–Stokes and energy equations are used to simulate the flow domain, and the VOF method is adopted to track the interface between the gas and the liquid phases. Simulations are divided into two cases. In the first case, the collapse stage alone is considered in both spherical and nonspherical situations with different conditions of bubble radius and surface tension. According to the results, surface tension has no significant effects on the flow pattern and collapse rate. In the second case, both the growth and collapse stages of bubbles with different initial radii and surface tensions are considered. In this case surface tension affects the growth stage considerably and, as a result, the jet velocity and collapse time decrease with increasing surface tension coefficient. This effect is more significant for bubbles with smaller radii.     

瞬态空化泡演变过程的数值模拟

采用边界积分方程方法，对无粘流体中三个空化泡以及自由面附近二个空化泡相互作用的演变过程进行了数值模拟。计算中边界用二阶有限元离散，影响系数矩阵非对角线元素用六点高斯数值积分方法计算，对第一类、第二类完全椭圆积分用高次多项式近似，对计算系数矩阵对角线元素中遇到的奇异积分进行了特殊处理。结果表明，在不同的给定参数下，空化泡的溃灭形态各异，柱状射流和环形射流都有可能发生，使空化泡演变成双泡或环形泡。     

Bubble collapse in compressible fluids using a spectral element marker particle method. Part 2. Viscoelastic fluids

This paper is concerned with the development of a high‐order numerical scheme for two‐phase viscoelastic flows. In the companion paper, herein referred to as Part 1, the scheme is applied to the modelling of two‐phase Newtonian flows. The particular problem of the collapse of a 2D bubble in the vicinity of a rigid boundary is considered. Attention is given to the construction of the most general form of the compressible Oldroyd B model that is consistent with the compressible Newtonian and upper‐convected Maxwell models in the appropriate limits. The governing equations are discretized using the spectral element method, and the two phases are modelled using a marker particle method. A comprehensive set of results is presented for the problem of bubble collapse near a rigid wall, and qualitative agreement is obtained with other numerical studies and experimental observations. Viscoelastic effects that are predicted include increased bubble oscillation with increasing Weissenberg number and considerable bubble deformation and cusping near the wall. Most importantly, it has been shown that viscoelasticity has the ability to prevent jet formation and therefore is likely to have a mitigating effect on cavitation damage. Copyright © 2012 John Wiley & Sons, Ltd.     

Vibrations of a gas-filled spherical cavity in pseudoplastic and dilatant liquids

Questions of the dynamics of bubbles in a liquid are connected with problems of cavitation [1]. In connection with cavitation phenomena in non-Newtonian media, in particular in polymeric liquids [2, 3], a study is made of the pulsations of a bubble in a polymeric liquid with an exponential rheological law. The equation of the motion of the boundary of the gas cavity is integrated numerically; here, the cases of pseudo-plastic and dilatant liquids are discussed separately. The results obtained can be used in the analysis of acoustical cavitation in aqueous solutions of polymers.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 146–148, January–February, 1975.     

Modeling and computation of cavitation in vortical flow

The paper presents an investigation of Euler–Lagrangian methods for cavitating two-phase flows. The Euler–Euler methods, widely used for simulations of cavitating flows in ship technology, perform well in regions of moderate flow changes but fail in zones of strong, vortical flow. Reasons are the strong approximations of cavitation models in the Euler concept. Alternatively, Euler–Lagrangian concepts enable more detailed formulations for transport, dynamics and acoustic of discrete vapor bubbles. Test calculations are performed to study the influence of different parameters in the equations of motion and in the Rayleigh–Plesset equation for bubble dynamics. Results confirm that only Lagrangian models are able to describe correctly the bubble behavior in vortices, while Eulerian results deviate strongly. Lagrangian formulations enable additionally the determination of acoustic pressure of cavitation noise. Two-way coupling between the phases is required for large regions of the vapor phase. A new coupling concept between continuous fluid flow and discrete bubble phase is developed and demonstrated for flow through a nozzle. However, the iterative coupling between the phases via volume fractions is computationally expensive and should therefore be applied only in regions where Eulerian treatment fails. A corresponding local concept for combination with an Euler–Euler method is outlined and is in progress.     

Cavitation bubble dynamics — new tools for an intricate problem

With the help of laser produced bubbles in water and high speed photography and holography sophisticated experiments on cavitation bubble dynamics can be conducted. The observation of a bubble vortex ring after jet formation upon collapse of a spherical bubble in front of a plane solid boundary is reported. The vortex ring may expand and contract several times until it disintegrates into a ring of bubbles by some instability finally taking over. A critical discussion of our qualitative understanding of jet formation is included. In a second part the problem of the acoustic cavitation noise spectrum is discussed. Numerically obtained ‘visible cavitation noise’ plots from a single bubble already resemble those obtained experimentally from acoustic cavitation. A discussion shows that the theory should be extended to self-consistency.     

Bubble dynamics: a review and some recent results

Several aspects of small-amplitude oscillations of bubbles containing gas, vapor, or a gas-vapor mixture are discussed. An application to pressure-wave propagation in a bubbly liquid is described. Nonlinear forced oscillations are considered in the light of recent research on forced oscillations of nonlinear systems. The growth of vapor bubbles, an extension of the Rayleigh-Plesset equation to non-Newtonian liquids and appreciable mass transfer at the interface, and a boundary integral numerical method for nonspherical cavitation bubble dynamics are also briefly discussed.     

水下爆炸近壁面流场局部空化形成机理

为深入认识水下爆炸近壁面流场局部空化的形成机理,采用自行研制的转镜式分幅相机,获得了炸药水下爆炸近壁面流场局部空化效应的光学图像,结合数值模拟和Taylor平面波理论、空泡动力学理论,分析了近壁面空化效应的形成过程。结果表明：界面反射的稀疏波作用和水中空化核的膨胀发展是水下爆炸近壁面流场空化效应形成的原因;外界流场压力对空泡初期膨胀运动影响较小,对空泡后期运动行为影响较大;低压环境下不同尺度空泡的运动行为存在较大差异,小尺度空泡（半径小于10μm）在低压环境下处于快速膨胀、溃灭状态,对流场空化影响较小;大尺度空泡（半径大于10μm）可失去稳定性,半径持续增大,对流场空化区的形成影响较大;水中不同尺寸空泡空间分布的随机性可导致空化区成长过程中呈现非规则形状。     

Particle tracking velocimetry of the flow field around a collapsing cavitation bubble

The velocity field in the vicinity of a laser-generated cavitation bubble in water is investigated by means of particle tracking velocimetry (PTV). Two situations are explored: a bubble collapsing spherically and a bubble collapsing aspherically near a rigid wall. In the first case, the accuracy of the PTV method is assessed by comparing the experimental data with the flow field around the bubble as obtained from numerical simulations of the radial bubble dynamics. The numerical results are matched to the experimental radius–time curve extracted from high-speed photographs by tuning the model parameters. Trajectories of tracer particles are calculated and used to model the experimental process of the PTV measurement. For the second case of a bubble collapsing near a rigid wall, both the bubble shape and the velocity distribution in the fluid around the bubble are measured for different standoff parameters γ at several instants in time. The results for γ > 1 are compared with the corresponding results of a boundary-integral simulation. For both cases, good agreement between simulation and experiment is found.     

The Peculiar Dynamics of Cavitation Bubbles

Details from cavitation bubble dynamics are reported: jet formation, counterjet formation, shock wave radiation and light emission. Multiple shock wave radiation from single bubble collapse with jet formation could be time resolved by high speed photography with 20 million frames per second. An explanation of counterjet formation is given. Pictures of the light emission (sonoluminescence) in acoustic cavitation are presented.     

Level set method for numerical simulation of a cavitation bubble,its growth,collapse and rebound near a rigid wall

A level set method of non-uniform grids is used to simulate the whole evolution of a cavitation bubble, including its growth, collapse and rebound near a rigid wall. Single-phase Navier–Stokes equation in the liquid region is solved by MAC projection algorithm combined with second-order ENO scheme for the advection terms. The moving interface is captured by the level set function, and the interface velocity is resolved by “one-side” velocity extension from the liquid region to the bubble region, complementing the second-order weighted least squares method across the interface and projection inside bubble. The use of non-uniform grid overcomes the difficulty caused by the large computational domain and very small bubble size. The computation is very stable without suffering from large flow-field gradients, and the results are in good agreements with other studies. The bubble interface kinematics, dynamics and its effect on the wall are highlighted, which shows that the code can effectively capture the “shock wave”-like pressure and velocity at jet impact, toroidal bubble, and complicated pressure structure with peak, plateau and valley in the later stage of bubble oscillating. The project supported by the National Natural Science Foundation of China (10272032 and 10672043). The English text was polished by Keren Wang.