## Shock formation in the collapse of a vapor nano-bubble

F. Magaletti, L. Marino, and C.M. Casciola ∗ Dipartimento di Ingegneria Meccanica e Aerospaziale, Università di Roma La Sapienza , Via Eudossiana 18, 00184 Roma Italy

(Dated: December 2, 2014)

In this Letter a diffuse-interface model featuring phase change, transition to supercritical conditions, thermal conduction, compressibility effects and shock wave propagation is exploited to deal with the dynamics of a cavitation bubble. At variance with previous descriptions, the model is uniformly valid for all phases (liquid, vapor and supercritical) and phase transitions involved, allowing to describe the non-equilibrium processes ongoing during the collapse. As consequence of this unitary description, rather unexpectedly for pure vapor bubbles, the numerical experiments show that the collapse is accompanied by the emission of a strong shock wave in the liquid and by the oscillation of the bubble that periodically disappears and reappears, due to transition to super/sub critical conditions. The mechanism of shock wave formation is strongly related to the transition of the vapor to supercritical state, with a progressive steepening of the compression wave to form the shock which is eventually reflected as an outward propagating wave in the liquid.

Vapor bubble collapse is a fascinating classical problem [1] involving vapor-liquid phase transition and extreme pressures and temperatures [2, 3]. Typical experiments concern ultra-fast imaging and the analysis of light and sound emitted after the collapse [4-6], see also the review [7]. Both free cavitation bubbles and nano bubbles at solid walls are increasingly investigated [8-10]. In ordinary conditions, gas bubble nucleation is associated with the metastability of the mixture of liquid and dissolved gas which, once the free energy barrier between the two states is overcome and the critical nucleus is formed, evolves toward a finite size bubble. The same mechanism is at work in forming pure vapor bubbles when the (ultra pure) liquid is kept in metastable conditions [11], i.e. its pressure is below the equilibrium vapor pressure at the given temperature. In these conditions, away from solid walls (see e.g. [12] for the role of asperities on solid surfaces as a catalyst of bubble nucleation), local density fluctuations can generate the critical vapor nucleus from which the eventual bubble is formed.

Intermingled phenomenologies, [13, 14], such as interface dynamics [15, 16], thermodynamics of phase change [17], and dissolved gas diffusion [18], are a challenge to theoretical modeling of cavitation. The available descriptions combine two distinct adjoining regions, liquid and vapor phase, respectively, with vapor pressure taken to be the saturation pressure [19] and phase transition accounted for through suitable kinetic equations and latent heat release [18].

Contrary to available models, the diffuse interface approach discussed in the present Letter encompasses all phases (liquid, vapor and supercritical) and phase transitions involved, embedding capillary forces, compressibility effects and shock wave propagation. The approach enables an unprecedented analysis of collapse, where the bubble interface speed may exceed the speed of sound. This leads to the formation of a shock wave focused towards the bubble that is successively reflected back in the liquid. Latent heat of condensation and rapid compression locally bring the vapor in supercritical conditions. This explains the observed rebounds usually considered a typical feature of incondensable gas bubbles. Indeed, the liquid-vapor interface may disappear and reappear again, according to the local thermodynamic conditions.

The physical model. We exploit an unsteady diffuse interface description [20] of the multiphase flow based on the van der Waals gradient approximation of the free energy functional [21, 22]

$$F [ \rho, \theta ] = \int _ { \Omega } \left ( \hat { f } _ { 0 } \left ( \rho, \theta \right ) + \frac { \lambda } { 2 } | \boldsymbol \nabla \rho | ^ { 2 } \right ) \, d V \,, \quad ( 1 )$$

where λ is a coefficient controlling surface tension and interface thickness. For a van der Waals fluid, the free energy per unit volume at temperature θ and density ρ is

$$\hat { f } _ { 0 } \left ( \rho, \theta \right ) = \bar { R } \rho \theta \left [ - 1 + \log \left ( \frac { \rho \, K \, \theta ^ { \delta } } { 1 - b \rho } \right ) \right ] - a \rho ^ { 2 } \,, \quad ( 2 )$$

with δ = ¯ R/c v , ¯ R the gas constant, c v the constant volume specific heat, a and b the van der Waals coefficients, with K a suitable constant [23]. The evolution is governed by the conservation equations for mass ∂ ρ t + ∇ · ( ρ u ) = 0 , momentum ∂ t ( ρ u ) + ∇ · ( ρ u ⊗ u ) = ∇ · τ , and total energy, ∂ E t + ∇ · ( u E ) = ∇ · [ τ · u -λρ ∇ · u ∇ ρ + k ∇ θ ] , where k is the thermal conductivity. The stress tensor, τ = -p 0 I + τ v + τ c , with p 0 = ρ ∂ f /ρ /∂ρ 2 ( ˆ 0 ) | θ the pressure and τ v the classical viscous component, features the diffused capillary stress (see Supplemental Material (SM) [24]),

$$\tau ^ { c } = \lambda \left [ \left ( \frac { 1 } { 2 } | \nabla \rho | ^ { 2 } + \rho \nabla ^ { 2 } \rho \right ) I - \nabla \rho \otimes \nabla \rho \right ] \,.$$

τ c accounts for surface tension σ , defined starting from the deviatoric component of the equilibrium contribution ( τ eq = -p 0 I + τ c ) to the stress tensor, τ dev = τ eq -

1 / 3 tr ( τ eq ) = -λ ∇ ρ ⊗ ∇ ρ , as the integral along a onedimensional cut normal to the interface

$$\sigma = - \int _ { - \infty } ^ { + \infty } \nu \cdot \left ( \tau ^ { d e v } \cdot \nu \right ) \, \mathrm d n = \int _ { - \infty } ^ { + \infty } \lambda \left ( \frac { \mathrm d \rho } { \mathrm d n } \right ) ^ { 2 } \mathrm d n \,, \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

with n the normal distance across the interface, and ν the unit normal to the cut (tangent to the interface). In fact, the integral is restricted to the region of thickness ∆ n ∝ √ λ where d ρ/ d n is substantially concentrated, resulting in σ = ∫ ρ L ρ V √ 2 λw ρ dρ ( ) ∝ √ λ [24, 25], where L and V denote liquid and vapor, respectively.

The thermodynamic condition ( ρ, θ ) identifies completely the bulk state, and Eq. (2) provides the latent heat released during the phase change, θ ∆ η , with η = -∂f /∂θ 0 | ρ the entropy per unit mass. As a consequence, the present model intrinsically features the phase transition whenever the local thermodynamic conditions are appropriate, as shown by the temperature equation

$$\rho c _ { v } \frac { D \theta } { D t } = - \rho \theta \frac { \partial \eta } { \partial \rho } \Big | _ { \theta } \frac { D \rho } { D t } + \tau ^ { v } \colon \nabla \mathbf u + \nabla \cdot ( k \nabla \theta ) \. \ \ ( 3 ) \quad \text{as th.}$$

Along the phase transition, the factor θ∂η/∂ρ | θ , basically consists of the amount of heat required to change the state by dρ . Assuming, for the sake of definiteness, an isothermal transformation in subcritical conditions, the integral of this term from vapor to liquid state is indeed the latent heat θ ∆ η . Consistently with physical expectation, the van der Waals model prescribes a latent heat release only in subcritical conditions. In supercritical conditions, the same term is related to the isochoric thermal pressure coefficient (1 /p )( ∂η/∂v ) | θ = (1 /p )( ∂p/∂θ ) | v , with v the specific volume. Clearly, during the actual evolution, the thermodynamic transformation followed by a material element is not isothermal, the actual transformation being dictated by mass, momentum and energy conservation.

The model can in principle be extended to a more sophisticated equation of state (EoS) than the van der Waals', when a direct comparison with experiments is sought for, see, e.g., [26] where an empirical EoS is used to include dissociation and ionization effects which become relevant in the late stage of the collapse.

Simulation details . The system of conservation laws to be solved constitutes a quite non-standard problem that calls for specialized numerical techniques. We list here a few issues to provide a flavor of the numerics used in the simulations, see Supplemental Material, Section B [24] for a more complete discussion of the numerical method: i) The extremely thin liquid-vapor interface calls for a high numerical resolution; ii) Liquid, vapor, and supercritical fluid compressibility gives rise to shock waves propagating in a non-uniform, multiphase environment; iii) The system manifests a compound nature, partly controlled by acoustics (hyperbolic behavior) and partly induced by viscosity and capillarity (diffusion and dispersion, respectively); iv) Although the sound speed, c 2 = ∂p /∂ρ 0 | η , is well defined in most of the phase space ( c 2 &gt; 0 ), a region exists below the spinodal where c 2 &lt; 0 . This hybrid behavior [27] is cause of failure for the standard hyperbolic solvers.

Bubble dynamics and shock evolution. The model described above is used to perform numerical simulations of spherically symmetric, pressure-induced bubble collapse. A vapor nano-bubble with radius R eq = 100nm is initially in equilibrium and the collapse is initiated by an overpressure enforced on the liquid, ( p ∞ -p e ) /p e = ∆ p/p e &gt; 0 , where p e is the liquid equilibrium pressure. Cases differing for overpressure and for thermal conductivity, measured by Pr = 3 µR/ ¯ (8 k ) , are considered to analyze their effects on the collapse dynamics. In all cases the initial temperature is θ /θ e c = 0 5 . and the surface tension is σ/ p R ( c eq ) = 0 045 . , corresponding to a realistic value σ = 0 09 . N/m for water. Here and throughout the subscript c denotes critical values.

Along the evolution, the bubble radius is measured as the distance of the liquid region from the center of the bubble, see Fig. 1 where data are shown for different ∆ p/p e and Pr . The dynamics consists of a sequence of rebounds and collapses associated with shock formation. The collapse time of a macroscopic bubble is estimated as t c = 0 915 . R eq √ ρ ∞ / p ( ∞ -p v ) , with p v the bubble equilibrium pressure, where capillary, viscous and compressibility effects are neglected [13, 14]. For nano-bubbles, however, surface tension is crucial and the numerical results suggest the scaling

$$\theta ) | _ { v }, \quad t _ { c } = 0. 9 1 5 R _ { e q } \sqrt { \frac { \rho _ { \infty } } { p _ { \infty } - p _ { v } + 2 \sigma / R _ { e q } } } = 0. 9 1 5 R _ { e q } \sqrt { \frac { \rho _ { \infty } } { \Delta p } } \.$$

Before the first collapse, the radius evolution is independent of thermal conductivity with a slight sensitivity to the overpressure. Although predicted by models of incondensable gas bubbles [14], rebounds are missed by simplified models which neglect the inner vapor dynamics. The rebounds are affected by thermal conductivity and overpressure, Fig. 1. The radius where the first collapse phase ends, and the successive rebound starts, increases with the overpressure, suggesting the presence of an incondensable gaseous phase inside the bubble. The increase of the overpressure leads to faster dynamics, see the expression of t c , resulting in increased pressure inside the bubble. An enhanced thermal conductivity, solid lines in Fig. 1, reduces the subsequent oscillations of the bubble by diffusing thermal energy from the hotter bubble to the colder liquid, recovering the isothermal Keller model with no-rebounds in the limit Pr → 0 .

The shock position s = ( ) s t is provided in the inset of Fig. 1. The seemingly different velocity is an artifact of the ∆ p -dependent time scale t c . In fact w = ˙ s is fairly constant. Away from the bubble, the shock propagates in the still liquid that, after the expansion wave, relaxed back to p ∞ , θ ∞ = θ e . The shock speed w is determined

Figure 1. Top: Time evolution of the bubble radius for various cases with the prediction of the Keller model [15] shown by the dotted line. The collapse of the bubble is induced by an overpressure in the liquid, ∆ = p p ∞ -p e . The other control parameters are Pr defined in the text, a measure of thermal conductivity, and Re = √ p ρ R c c eq /µ and C = λρ / 2 c ( p R c 2 eq ) , playing the role of a Reynolds and a Weber number, see Supplemental Material, Section A for their detailed definition. In correspondence with the first bubble rebound a shock wave is radiated in the liquid. The shock position s t ( ) is shown in the inset. Bottom: Successive snapshots of the system configuration. Density field (upper row) and pressure gradient intensity (lower row) in arbitrary units for p ∞ -p v = 0 01 . p c and Pr = 0 2 . . The plots in the lower row highlight the position of bubble interface and radiated shock.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0643v2-figure-6e824968-3b55-48f4-ae3b-40180ad60da8.jpg)

by the state ahead of the shock ( p ∞ , ρ ∞ , u ∞ = 0 ) and by an additional parameter, the density ρ b behind the shock, say, see [28] and [23] for details concerning a van der Waals fluid. The small compressibility of the liquid ( ρ b -ρ ∞ ) /ρ ∞ glyph[lessmuch] 1 , allows the linearization w = c ∞ + α ∞ ( ρ b -ρ ∞ ) /ρ ∞ , where c ∞ is the unperturbed sound speed in the liquid and α ∞ ( p ∞ , ρ ∞ ) = ρ ∞ dw/dρ | ∞ . It turns out that w glyph[similarequal] c ∞ for the cases explicitly reported here (more precisely, ( w -c ∞ ) /c ∞ ≤ 3% ).

The main plot in the top panel of Fig. 2 shows the pressure profile through the bubble center for the initial phase of the process up to the bubble collapse. The corresponding density profiles are provided in the inset. Initially (solid, dashed and dash-dotted lines), the location of the interface is identified by an extremely sharp density drop. The bubble shrinks while the vapor den-

Figure 2. Radial profiles before the first bubble rebound. The different line styles correspond to successive time instants (solid: t/t c = 0 , dashed: t/t c = 0 64 . , dash-dotted: t/t c = 0 85 . , dotted: t/t c = 0 96 . , long-dashed: t/t c = 0 97 . , dash-dotdotted: t/t c = 0 98 . ). The top panel shows the pressure with density in the inset. The bottom panel shows the temperature and the two contributions to the heat release in the insets: a) the non dimensional latent heat q ∗ l = H q ρ/ p l ˙ ( 3 / 2 c ρ -1 / 2 c R -1 eq ) ; b) α ∗ = αρ/ p ˙ ( 3 / 2 c ρ -1 / 2 c R -1 eq ) .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0643v2-figure-091af75a-bd59-4e84-b0cf-00c6969372bf.jpg)

Figure 3. Radial profiles after the first bubble rebound. The different line styles correspond to successive time instants (dash-dotted: t/t c = 1 005 . , dotted: t/t c = 1 01 . , long-dashed: t/t c = 1 02 . , dash-dot-dotted: t/t c = 1 03 . ). The top panel shows the temperature with density in the inset. The bottom panel shows the pressure. In the inset the Hickling-Plesset power law for the shock peak attenuation is compared with the simulations.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0643v2-figure-5a145714-d846-499c-b635-14353f8190fc.jpg)

sity, pressure, and temperature (main plot in the bottom panel) increase. Successively, the shrinkage accelerates, the gaseous phase is compressed, its temperature raises and the fluid transitions to supercritical conditions ( p/p c &gt; 1 , θ/θ c &gt; 1 ). During this phase the inner vapor core is surrounded by a shell of supercritical fluid whose density increases through a strong density gradient to eventually adjoin the external liquid. Later, the vapor disappears altogether, transformed into a low density supercritical fluid. Subsequently pressure and temperature peaks in the profile keep increasing. The late collapse phase is dominated by a strong pressure wave propagat-

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0643v2-figure-2df62d2e-ea0e-453b-951d-7167d539acd9.jpg)

c

Figure 4. Comparison between the bubble radius and the position of the liquid-vapor interface. The bubble radius is defined irrespective of the nature of the gaseous phase inside the bubble, either the vapor or the supercritical fluid. The interface is defined as the boundary between vapor and liquid. The bubble radius changes continuously in time, up to complete condensation and disappearance of the bubble (blue region). The interface disappears (orange regions) when the fluid inside the bubble transitions to supercritical state.

ing in a homogeneous supercritical fluid which focuses into a converging shock (see Supplemental Material, Section C [24] for the comparison with a supercritical converging shock). Pressure and temperature extrema are reached at the first collapse time t c , when the inner low density core disappears.

In the inset of the bottom row of Fig. 2 we analyze the thermal aspects, Eq. (3), in particular the heat release rate -ρθ ∂η/∂ρ ( ) | θ Dρ/Dt . Let us consider the splitting -ρθ∂η/∂ρ ∣ ∣ θ = H q glyph[lscript] + α , where H is the characteristic function of the coexistence region ( H = 1 for states below the coexistence curve and 0 otherwise). The contribution due to the phase change is q glyph[lscript] = ρθ η [ V ( θ ) -η L ( θ )] / ρ [ L ( θ ) -ρ V ( θ )] and the corresponding heat release rate H q ρ glyph[lscript] ˙ is displayed in inset a) of Fig. 2 where αρ ˙ is too small to be appreciated in the plot. At later times, inset b), the fluid in the bubble becomes supercritical, i.e. H q glyph[lscript] = 0 , and the term αρ ˙ becomes substantial (for ideal gases α = ¯ Rθ , while in general α = θβ / κ ρ p ( θ ) , where thermal compressibility and thermal expansion coefficient are κ θ = -(1 /v ∂v/∂p ) | θ and β p = (1 /v ) ∂v/∂θ | p , respectively).

After reaching the bubble center, the shock is reflected and propagates toward the liquid, Fig. 3. The pressure peaks ahead of the shock follows the scaling law p s ∝ 1 /s , inset of Fig. 3, as predicted by the compressible HicklingPlesset model [29] for incondensable gas bubbles. After reflection, the fluid near the bubble center expands back and reduces quickly its temperature and pressure coming back to vapor with reappearance of the bubble. The expansion continues up to a maximum radius when the motion reverts and the process repeats itself.

As already mentioned, during the dynamics the liquid- vapor interface may disappear and reappear several times. In Fig. 4 the interface position, defined only when vapor contiguous to the liquid exists, is compared with the bubble radius, which is well defined also when the inner fluid becomes supercritical. The time intervals corresponding to supercritical conditions are highlighted in orange. Soon after the shock emission, the inner core transitions back into vapor and the liquid-vapor interface appears again. The phenomenon repeats until a full condensation is achieved.

Discussion &amp; conclusions. Phase change and transition to supercritical conditions play a crucial role in the collapse of a vapor nanobubble. Indeed, independently of the intensity of the initial overpressure, a strong pressure and temperature increase is experienced that induces the transition to incondensable gaseous state. As a consequence, a vapor bubble substantially resembles an incondensable gas bubble, making the boundary between the two kinds less sharply defined than usually assumed. The pressure and temperature peaks increase with ( p ∞ -p e ) /p e as the strength of the emitted shock wave does, see Table I. At fixed thermal conductivity, a limiting overpressure exists below which the bubble condenses altogether. Above the critical overpressure oscillations set in, with the bubble periodically reforming and emitting a shock upon collapse. From the above considerations it should be expected that a collapsing bubble could trigger a synchronized collapse of its neighbors. Indeed, the pressure at the distance r = R eq from the bubble center is substantially larger than the initial overpressure, Table I. Accounting for the 1 /r decay of the pressure peak, the pressure of the wave exceeds the initial overpressure in a region extending for, typically, 20 R eq .

The present model for vapor-bubbles collapse can be extended under several respects: a) Incondensable gas dissolved in the liquid can be taken into account by extending the basic free energy functional (1) to include the composition of the mixture; b) More realistic transport coefficients can be assumed, e.g. dependence of viscosity and thermal conductivity on thermodynamic conditions can be included; c) A more general EoS can be adopted to take into account dissociation and ionization effects which are expected to quantitatively modify the dynam-

Table I. Pressure and temperature peaks in the collapsing bubble as a function of the initial overpressure. The intensity of the pressure wave at one radius from the bubble center is shown in the last column.

|   ( p ∞ - p e ) / | p e |   p max /p c |   θ max /θ c |   ( p r = Req - p e ) / | p e | |
|-------------------------|--------------|--------------|---------------------------------|
|                   1.43  |       1764   |         8.58 |                            31.7 |
|                   0.95  |        339   |         4.53 |                            22   |
|                   0.63  |        190   |         3.53 |                            15.9 |
|                   0.073 |         52.8 |         2.01 |                             6.6 |

ics of the collapse, reducing peak temperature and pressure, while maintaining the overall phenomenology basically unchanged.

Acknowledgments. The research leading to these results has received funding from the European Research Council under the European Union's Seventh Framework Programme (FP7/2007-2013) / ERC Grant agreement n. [339446].

- ∗ carlomassimo.casciola@uniroma1.it
- [1] L. Rayleigh, The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science 34 , 94 (1917).
- [2] R. P. Taleyarkhan, C. West, J. Cho, R. Lahey, R. Nigmatulin, and R. Block, Science 295 , 1868 (2002).
- [3] C.-D. Ohl, Physics 3 , 65 (2010).
- [4] R. Pecha and B. Gompf, Physical review letters 84 , 1328 (2000).
- [5] G. Sankin, W. Simmons, S. Zhu, and P. Zhong, Physical review letters 95 , 034501 (2005).
- [6] K. Weninger, C. Camara, and S. Putterman, Physical Review E 63 , 016310 (2000).
- [7] M. P. Brenner, S. Hilgenfeldt, and D. Lohse, Reviews of modern physics 74 , 425 (2002).
- [8] M. P. Brenner and D. Lohse, Physical review letters 101 , 214505 (2008).
- [9] J. H. Weijs and D. Lohse, Physical review letters 110 , 054501 (2013).
- [10] X. Zhang, H. Lhuissier, O. R. Enríquez, C. Sun, and D. Lohse, Langmuir 29 , 9979 (2013).
- [11] K. Davitt, A. Arvengas, and F. Caupin, EPL (Europhysics Letters) 90 , 16002 (2010).
- [12] A. Giacomello, M. Chinappi, S. Meloni, and C. M. Cas-
- ciola, Langmuir 29 , 14873 (2013).
- [13] C. E. Brennen, Cavitation and bubble dynamics (Cambridge University Press, 2013).
- [14] M. S. Plesset and A. Prosperetti, Annual Review of Fluid Mechanics 9 , 145 (1977).
- [15] J. B. Keller and I. I. Kolodner, Journal of Applied physics 27 , 1152 (1956).
- [16] M. S. Plesset and R. B. Chapman, J. Fluid Mech 47 , 283 (1971).
- [17] S. Fujikawa and T. Akamatsu, Journal of Fluid Mechanics 97 , 481 (1980).
- [18] I. Akhatov, O. Lindau, A. Topolnikov, R. Mettin, N. Vakhitova, and W. Lauterborn, Physics of Fluids 13 , 2805 (2001).
- [19] Y. Hao and A. Prosperetti, Physics of Fluids (1994present) 11 , 2008 (1999).
- [20] D. Anderson, G. McFadden, and A. Wheeler, Annual review of fluid mechanics 30 , 139 (1998).
- [21] F. Dell'Isola, H. Gouin, P. Seppecher, et al. , Comptes Rendus de l Académie des Sciences-Series IIB-Mechanics 320 (1995).
- [22] D. Jamet, O. Lebaigue, N. Coutris, and J. Delhaye, Journal of Computational Physics 169 , 624 (2001).
- [23] N. Zhao, A. Mentrelli, T. Ruggeri, and M. Sugiyama, Physics of fluids 23 , 086101 (2011).
- [24] see Supplemental Material at http://correct.site. url for details, .
- [25] J. Cahn and J. Hilliard, J. Chem. Phys 28 , 258 (1958).
- [26] W. C. Moss, D. B. Clarke, J. W. White, and D. A. Young, Physics of Fluids (1994-present) 6 , 2979 (1994).
- [27] M. Affouf and R. E. Caflisch, SIAM Journal on Applied Mathematics 51 , 605 (1991).
- [28] R. Courant and K. O. Friedrichs, Supersonic flow and shock waves , Vol. 21 (Springer, 1976).
- [29] R. Hickling and M. S. Plesset, Physics of Fluids (19581988) 7 , 7 (2004).
Francesco Magaletti Luca Marino Carlo Massimo Casciola
Shock formation in the collapse of a vapor nano-bubble
In this paper a diffuse-interface model featuring phase change, transition to
supercritical conditions, thermal conduction, compressibility effects and shock
wave propagation is exploited to deal with the dynamics of a cavitation bubble.
At variance with previous descriptions, the model is uniformly valid for all
phases (liquid, vapor and supercritical) and phase transitions involved,
allowing to describe the non-equilibrium processes ongoing during the collapse.
As consequence of this unitary description, rather unexpectedly for pure vapor
bubbles, the numerical experiments show that the collapse is accompanied by the
emission of a strong shock wave in the liquid and by the oscillation of the
bubble that periodically disappears and reappears, due to transition to
super/sub critical conditions. The mechanism of shock wave formation is
strongly related to the transition of the vapor to supercritical state, with a
progressive steepening of the compression wave to form the shock which is
eventually reflected as an outward propagating wave in the liquid.