Evolution of deviations from the spherical shape of a vapor bubble in supercompression

This paper considers the evolution of small deviations of a cavitation bubble from a spherical shape during its single compression under conditions of experiments on acoustic cavitation of deuterated acetone. Vapor motion in the bubble and the surrounding liquid is defined as a superposition of the spherical component and its non-spherical perturbation. The spherical component is described taking into account the nonstationary heat conductivity of the liquid and vapor and the nonequilibrium nature of the vaporization and condensation on the interface. At the beginning of the compression process, the vapor in the bubble is considered an ideal gas with a nearly uniform pressure. In the simulation of the high-rate compression stage, realistic equations of state are used. The non-spherical component of motion is described taking into account the effect of liquid viscosity, surface tension, vapor density in the bubble, and nonuniformity of its pressure. Estimates are obtained for the amplitude of small perturbations (in the form of harmonics of degree n = 2, 3, ... with the wavelength λ = 2πR/n, where R is the bubble radius) of the spherical shape of the bubble during its compression until reaching extreme values of pressure, density, and temperature. These results are of interest in the study of bubble fusion since the non-sphericity of the bubble prevents its strong compression.     

Evolution of a Small Distortion of the Spherical Shape of a Gas Bubble under Strong Expansion-Compression

The evolution of a small distortion of the spherical shape of a gas bubble which undergoes strong radial expansion-compression upon a single oscillation of the ambient liquid pressure under a harmonic law are analyzed by numerical experiments. It is assumed that the distortions of the spherical bubble shape are axisymmetric and have the form of individual spherical surface harmonics with numbers of 2–5. Bubble-shape oscillations prior to the beginning of expansion are taken into account. Generally, the distortion value during bubble expansion-compression depends on the phase of bubble-shape oscillation at the beginning of the expansion (initial phase). Emphasis is placed on the dependence of the maximum distortions in the initial phase at certain characteristic times of bubble expansion-compression on the amplitude of the external excitation, liquid viscosity, and distortion mode (harmonic number). The parameters of the problem are typical of the stable periodic sonolumiescence of an individual air bubble in water at room temperature. An exception is the liquid pressure oscillation amplitude, which is varied up to values that are five times the static pressure. That large excitation amplitudes are beyond the stability threshold of periodic oscillations of spherical bubbles. Their consideration is of interest from the point of view of increasing the compression ratio of the bubble gas, i. e., increasing the maximum temperature and density achievable in the final compression stage.__________Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 4, pp. 17–28, July–August, 2005.     

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.     

Bubbles in liquids with phase transition

In the forthcoming second part of this paper a system of balance laws for a multi-phase mixture with many dispersed bubbles in liquid is derived where phase transition is taken into account. The exchange terms for mass, momentum and energy explicitly depend on evolution laws for total mass, radius and temperature of single bubbles. Therefore in the current paper we consider a single bubble of vapor and inert gas surrounded by the corresponding liquid phase. The creation of bubbles, e.g. by nucleation is not taken into account. We study the behavior of this bubble due to condensation and evaporation at the interface. The aim is to find evolution laws for total mass, radius and temperature of the bubble, which should be as simple as possible but consider all relevant physical effects. Special attention is given to the effects of surface tension and heat production on the bubble dynamics as well as the propagation of acoustic elastic waves by including slight compressibility of the liquid phase. Separately we study the influence of the three phenomena heat conduction, elastic waves and phase transition on the evolution of the bubble. We find ordinary differential equations that describe the bubble dynamics. It turns out that the elastic waves in the liquid are of greatest importance to the dynamics of the bubble radius. The phase transition has a strong influence on the evolution of the temperature, in particular at the interface. Furthermore the phase transition leads to a drastic change of the water content in the bubble. It is shown that a rebounding bubble is only possible, if it contains in addition an inert gas. In Part 2 of the current paper the equations derived are sought in order to close the system of equations for multi-phase mixture balance laws for dispersed bubbles in liquids involving phase change.     

Oscillations of a gas bubble in a non-Newtonian liquid under the action of an acoustic field

The evolution of the radius of a spherical cavitation bubble in an incompressible non-Newtonian liquid under the action of an external acoustic field is investigated. Non-Newtonian liquids having relaxation properties and also pseudoplastic and dilatant liquids with powerlaw equation of state are studied. The equations for the oscillation of the gas bubble are derived, the stability of its radial oscillation and its spherical form are investigated, and formulas are given for the characteristic frequency of oscillations of the cavitation hollow in a relaxing liquid. The equations are integrated numerically. It is shown that in a relaxing non-Newtonian liquid the viscosity may lead to the instability of the radial oscillations and the spherical form of the bubble. The results obtained here are compared with the behavior of a gas bubble in a Newtonian liquid.     

Effect of homogeneous condensation on the dynamics of a hot vapor bubble in a cold liquid jet

The problem of initiating cavitation bubbles in a cold liquid jet by injecting hot steam into high-pressure zone specially organized at the nozzle outlet is considered. Previously, in [1], a plane flowfield in which vapor bubbles were formed at the cusp of the cavity (high-pressure zone) and propagated together with the liquid along the axis of symmetry was considered. In certain cases, in the bubble expansion process the vapor temperature drops below the saturation temperature. In the present paper, vapor condensation in the bubble volume (homogeneous condensation) is also taken into account.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 6, pp. 56–61, November–December, 1996.     

自由场中液氮单空泡动力学特性的实验研究

本文专门设计搭建了低温介质空泡演化实验测试平台, 对液氮单空泡非定常演化过程和动力学特性开展了实验研究. 实验中利用电火花瞬态放电激发液氮汽化形成单空泡, 通过高速摄影系统对单空泡的瞬态特征进行了精细化捕捉. 为了进一步揭示低温介质独特的物理性质以及强热力学效应对单空泡演化过程的影响机制, 对比分析了在相同环境压力下, 77.41 K液氮和298.36 K水单空泡的演化过程和动力学特性. 基于实验得到空泡半径与界面速度等定量数据, 阐明了液氮单空泡球形与非球形演化阶段的非定常特性. 研究结果表明: (1) 在相同输入电压下, 液氮单空泡的整体尺寸比常温水更小, 当输入电压为400 V时, 液氮空泡的最大半径约为常温水空泡的0.69倍; 同时, 液氮单空泡经历了膨胀阶段?收缩阶段?振荡阶段以及上升阶段的演化过程. (2)液氮空泡的收缩过程主要由相界面的热传导主导, 没有明显的塌陷现象, 收缩阶段液氮空泡的最小收缩半径约为常温水的5.5倍. (3)在液氮空泡振荡初期, 空泡相界面传热增强, Rayleigh-Taylor不稳定与热力学效应共同引起了空泡界面的表面粗化效应; 在整个振荡阶段, 空泡界面附近存在破碎的小泡. 当输入电压较高时, 空泡底部的小泡数量显著增多. (4)由于液氮空泡浮力系数较大, 液氮空泡在演化后期空泡整体向上迁移显著, 液氮空泡底部收缩更快产生凹陷, 促使空泡变为环状.      

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.     

Free nonlinear oscillations of a thermally relaxing spherical gas bubble in an incompressible liquid

The quasi-adiabatic regime of free oscillation of a bubble in the presence of irreversible interphase heat transfer between the bubble and the ambient liquid is studied. On the basis of simplified model equations of a rarefield bubble mixture, a nonlinear-oscillation equation of the relaxation type is obtained. In constructing an exact particular solution of this equation, the heat transfer law associated with bubble compression is established. For studying the harmonic oscillations, the Krylov-Bogolyubov-Mitropol’skii asymptotic method is used. It is shown that, for a small bubble, the viscosity and heat transfer effects are of the same order. For a small bubble, the influence of these effects on the formation of the natural-oscillation frequency, which is small in the linear approximation, may be significant in the nonlinear formulation. For a large bubble, the influence of these effects is negligible in both approximations. For the approximate solution of the nonlinear equation, a uniformly valid second-order expansion is constructed.     

Forced oscillations of a gas bubble in a spherical volume of a compressible liquid

A spherically symmetric problem of oscillations of a single gas bubble at the center of a spherical flask filled with a compressible liquid under the action of pressure oscillations on the flask wall is considered. A system of differential-difference equations is obtained that extends the Rayleigh-Plesset equation to the case of a compressible liquid and takes into account the pressure-wave reflection from the bubble and the flask wall. A linear analysis of solutions of this system of equations is performed for the case of harmonic oscillations of the bubble. Nonlinear resonance oscillations and nearly resonance nonharmonic oscillations of the bubble caused by harmonic pressure oscillations on the flask wall are analyzed. Ufa Scientific Center, Russian Academy of Sciences, Ufa 450000. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 40, No. 2, pp. 111–118, March–April, 1999.     

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.     

Unsteady waves in a liquid containing vapor bubbles

Unsteady wave processes in vapor-liquid media containing bubbles are investigated taking into account the unsteady interphase heat and mass transfer. A single velocity model of the medium with two pressures is used for this, which takes into account the radial inertia of the liquid with a change in volume of the medium and the temperature distribution in it [1]. The system of original differential equations of the model is converted into a form suitable for carrying out numerical integration. The basic principles governing the evolution of unsteady waves are studied. The determining influence of the interphase heat and mass transfer on the wave behavior is demonstrated. It is found that the time and distance at which the waves reach a steady configuration in a vapor-liquid bubble medium are considerably less than the correponding characteristics in a gas-liquid medium. The results of the calculation are compared with experimental data. The propagation of acoustic disturbances in a liquid with vapor bubbles was studied theoretically in [2]. The evolution of waves of small but finite amplitude propagating in one direction in a bubbling vapor-liquid medium is investigated in [3, 4] on the basis of the generalization of the Burgers-Korteweg-de Vries equation obtained by the authors. An experimental investigation of shock waves in such a medium is reported in [5, 6], and the structure of steady shock waves is discussed [7].Translated from Izvestiya Akademii Nauk SSSR, Hekhanika Zhidkosti i Gaza, No. 5, pp. 117–125, September–October, 1984.     

Numerical analysis of an asymptotic model for the collapse of a vapor bubble

In this work, we analyze the thermal collapse of a vapor bubble immersed in a unbounded and subcooled liquid. In this thermal regime, controlled basically by Jakob number (Ja), we present an asymptotic limit of the governing equations by identifying the appropriate temporary and spatial scales to solve numerically the mathematical model. In the limit of Ja ≫ 1, the governing equations describe the spatial and temporal evolution of the adjacent thermal boundary layer to the radius of the bubble. In particular, we prove that the influence of curvature effects due to conductive and convective heat terms of the energy equation for the liquid are responsible to characterize the thermal collapse regime. The numerical results for the evolution of the nondimensional radius of the bubble, a, and the corresponding nondimensional temperature profiles, θ, for different values of the Ja, show that the ending collapse state has a singular behavior, which we have denoted as a “thermal runaway”.     

Direct 2D simulation of small gas bubble clusters: from the expansion step to the equilibrium state

A numerical strategy, based on an adaptive finite element method, is proposed for the direct two‐dimensional simulation of the expansion of small clusters of gas bubbles within a Newtonian liquid matrix. The velocity and pressure fields in the liquid are first defined through the Stokes equations and are subsequently extended to the gas bubbles. The liquid–gas coupling is imposed through the stress exerted on the liquid by gas pressure (ruled by an ideal gas law) and by surface tension. A level set method, combined with a mesh adaptation technique, is used to track liquid–gas interfaces. Many numerical simulations are presented. The single bubble case allows to compare the simulations to an analytical model. Simulations of the expansion of small clusters are then presented showing the interaction and evolution of the gas bubbles to an equilibrium state, involving topological rearrangements induced by Plateau's rule. Copyright © 2006 John Wiley & Sons, Ltd.     

Method for calculating the dynamics of vapor bubbles when a liquid boils in a centrifugal force field

The specific feature of the study of the dynamics of vapor bubbles during boiling of a liquid in a centrifugal force field is the fact that the velocity of a bubble is much greater than the rate of change of its radius, and its movement occurs in fields of variable pressure and underheating that have to be determined in the solution of the problem. In addition, when investigating processes occurring when liquid helium boils in a centrifugal force field, its thermodynamic parameters may be close to the critical values, and the dependences of the thermophysical properties of the liquid and vapor on the temperature and pressure must be taken into consideration. The equation of state of a substance close to its critical thermodynamic point cannot be an approximation to the equation of state of an ideal gas, as has been suggested in a series of articles. The nonequilibrium nature of the phase transition must be taken into consideration in the case of substances existing at near-critical parameters and substances with a low coefficient of accommodation. A marked deformation of the bubbles, which also has to be taken into account, will be observed in strong centrifugal force fields. Such studies have not appeared in the specialist journals. Equations of the two-temperature and two-velocity hydrodynamics of two-phase media in a one-dimensional form for substances obeying the equation of state for an ideal gas were discussed in [1, 2] with allowance for the dependence of the thermophysical properties on the temperature and pressure. In strong centrifugal force fields the one-dimensional approach is unacceptable and the flow of liquid around a buoyant bubble must be taken into account. A joint examination of the change in the temperature field with time in the vicinity of a vapor bubble with changes in its dimensions and position was made for the first time in [3–8]. The present article is an extension of the latter work and takes the aforementioned factors into account.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 66–71, July–August, 1984.     

On the deformation of two droplets in a quasisteady stokes flow

A theoretical method is given for the determination of the shape of two drops (bubbles) moving with constant velocities parallel to their line of centres in a quiescent viscous fluid. The Reynolds numbers for the motions within the fluids are assumed to be sufficiently small that the equations governing these motions are quasisteady Stokes' equations. It is also assumed that the maximum deviation of the interfaces from spherical form is small when compared with the radius of the “equivalent” spherical drop. The paper deduces the first-order pressure distribution exterior and interior to the droplets. Effects of fluid viscosities, capillary numbers and distances between the droplets are taken into account. Special attention is paid to the influence of a solid plane or solid sphere on the shape of a drop (bubble) approaching or receding away from the solid boundary. The obtained solutions may serve as a first iteration of an iterative procedure for determining more accurate flow fields, taking into account the deviation from sphericity of the deformed particles.     

Formation and evolution of the liquid-phase particle size spectrum in the presence of nonequilibrium homogeneous condensation

The evolution of the liquid droplet size distribution during nucleation, growth and re-evaporation is considered. A closed system of equations, which consistently takes into account the dependence of the growth rate on particle size and ensures a fairly accurate solution of the basic kinetic equation, is obtained for the four lowest moments of the distribution function. Comparative calculations of the condensation in an expanding volume of vapor for constant and periodic expansion rates are carried out on the basis of the system proposed and, moreover, by directly solving the kinetic equation.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 108–116, January–February, 1990.     

Dynamics of a bubble in a liquid under laser pulse action

Simulation was performed of the behavior of a vapor bubble in a liquid under laser irradiation in laboratory experiments. A mathematical model was developed to analyze the effect of heat conduction, diffusion, and mass transfer on the bubble dynamics under compression and expansion. It is found that at the stage of collapse, intense condensation occurs on the bubble wall, which results in a significant (more than 15fold) decrease in bubble mass and an increase in pressure (to 10⁵ atm) and temperature (to 10⁴ K(. Results of numerical calculations of the radius of the first rebound and the amplitude of the divergent shock wave in water are compared with experimental data. It is shown that small (:about 1%) additives of an incondensable gas lead to a considerable decrease in mass transfer on the bubble wall.     

Flow regime transition criteria for two-phase flow in a vertical annulus

In this work, a new flow regime transition model is proposed for two-phase flows in a vertical annulus. Following previous works, the flow regimes considered are bubbly (B), slug (S) or cap-slug (CS), churn (C) and annular (A). The B to CS transition is modeled using the maximum bubble package criteria of small bubbles. The S to C transition takes place for small annulus perimeter flow channels and it is assumed to occur when the mean void fraction over the entire region exceeds that over the slug–bubble section. If the annulus perimeter is larger that the distorted bubble limit the cap-slug flow regime will be considered since in these conditions it is not possible to distinguish between cap and partial-slug bubbles. The CS to C transition is modeled using the maximum bubble package criteria. However, this transition considers the coalescence of cap and spherical bubbles in order to take into account the flow channel geometry. Finally, the C to A transition is modeled assuming two different mechanisms, (a) flow reversal in the liquid film section along large bubbles; (b) destruction on liquid slugs or large waves by entrainment or deformation. In the S to C and C to A flow regime transitions the annulus flow channel is considered as a rectangular flow channel with no side walls. In all the modeled transitions the drift-flux model is used to obtain the final correlations. The final equations for every flow regime transition are easy to be implemented in computational codes and not experimental input is needed. The prediction accuracy of the newly developed model has been checked against air–water as well as boiling flow regime maps. In all the cases, the new developed model shows better predicting capabilities than the existing correlations most used in literature.