Homogeneous Bubble Growth
From Thermal-FluidsPedia
Bubble growth within a superheated liquid drop will be considered first. Once a growing vapor bubble’s radius reaches that unstable equilibrium, it grows spontaneously. During the early stage, when the bubble radius is small, the Laplace-Young equation indicates that the pressure differential across the interface is at its maximum value. The resulting high interfacial velocity leads to significant inertia terms in the momentum equation. Meanwhile, the temperature of the interface is close to the superheat temperature of the surrounding liquid, so heat transfer into the vapor bubble experiences the highest driving temperature differential that occurs at any time during the process. As a result, early bubble growth tends to be limited by inertia, or the exchange of momentum between the vapor and liquid phases. In the physical model shown in the figure on the right, the instantaneous radius and the liquid-vapor interface velocity at time t are R(t) and dR / dt, respectively. For an incompressible, radially symmetric inviscid flow, the continuity equation in the spherical coordinate system is

where u is the radial velocity in the liquid phase. Assuming , this velocity at the liquid-vapor interface is equal to the growth rate, i.e.,

Integrating eq. (1) over the interval of (R,r), where r is an arbitrary location in the liquid phase, one obtains

The momentum equation for the surrounding liquid is

Substituting eq. (3) into eq. (4) yields
![\frac{{{\rho _\ell }}}{{{r^2}}}\left[ {2R{{\left( {\frac{{dR}}{{dt}}} \right)}^2} + {R^2}\frac{{{d^2}R}}{{d{t^2}}}} \right] - \frac{{2{\rho _\ell }{R^4}}}{{{r^5}}}{\left( {\frac{{dR}}{{dt}}} \right)^2} = - \frac{{\partial p}}{{\partial r}} \qquad \qquad(5)](/encyclopedia/images/math/6/0/0/6003a48735d60ec5bcf3da4adfa290e1.png)
Integrating eq. (5) from the bubble surface (r = R) to infinity () yields
![R\frac{{{d^2}R}}{{d{t^2}}} + \frac{3}{2}{\left( {\frac{{dR}}{{dt}}} \right)^2} = \frac{1}{{{\rho _\ell }}}\left[ {{p_\ell }(R) - {p_\ell }(\infty )} \right] \qquad \qquad(6)](/encyclopedia/images/math/4/e/a/4eaef8877f154c012fa8af1e4c526c0f.png)
Substituting Nucleation and Inception Eq (1)

with pg = 0 into eq. (6) gives
![R\frac{{{d^2}R}}{{d{t^2}}} + \frac{3}{2}{\left( {\frac{{dR}}{{dt}}} \right)^2} = \frac{1}{{{\rho _\ell }}}\left[ {{p_v} - {p_\ell }(\infty ) - \frac{{2\sigma }}{R}} \right] \qquad \qquad(7)](/encyclopedia/images/math/a/a/1/aa18e591e725217f5f16497dc2d2bbbf.png)
which is referred to as Rayleigh’s equation. During the inertia-controlled stage of bubble growth, the surface tension term 2σ / R is negligible compared to The pressure difference
can be evaluated by using the Clausius-Clapeyron equation, i.e.,
![{p_v} - {p_\ell }(\infty ) = \frac{{{\rho _v}{h_{\ell v}}\left[ {{T_\infty } - {T_{sat}}[{p_\ell }(\infty )]} \right]}}{{{T_{sat}}[{p_\ell }(\infty )]}}\qquad\qquad(8)](/encyclopedia/images/math/5/d/f/5dfeb77153b942609896afd09c5553be.png)
Substituting eq. (8) into eq. (7) and neglecting the surface tension term 2σ / R in eq. (7), one can obtain an ordinary differential equation for the bubble radius:
![R\frac{{{d^2}R}}{{d{t^2}}} + \frac{3}{2}{\left( {\frac{{dR}}{{dt}}} \right)^2} = \frac{{{\rho _v}{h_{\ell v}}}}{{{\rho _\ell }}}\frac{{{T_\infty } - {T_{sat}}[{p_\ell }(\infty )]}}{{{T_{sat}}[{p_\ell }(\infty )]}} \qquad \qquad(9)](/encyclopedia/images/math/3/0/f/30f359ed5ca1ef20dc5bced4e77d81af.png)
which is subject to the following initial condition:

The instantaneous bubble radius can be obtained by solving eq. (9):
![R(t) = {\left\{ {\frac{2}{3}\frac{{{\rho _v}{h_{\ell v}}}}{{{\rho _\ell }}}\frac{{{T_\infty } - {T_{sat}}[{p_\ell }(\infty )]}}{{{T_{sat}}[{p_\ell }(\infty )]}}} \right\}^{\frac{1}{2}}}t \qquad \qquad(11)](/encyclopedia/images/math/d/8/a/d8aaa857a5df1b2e21c23afe34515c12.png)
which indicates that the radius of the bubble is a linear function of time during the early stage of the bubble’s growth.
During the later stage of bubble growth, the inertial force becomes insignificant as the pressure difference (remember Laplace-Young’s equation) diminishes and the interfacial motion slows. At the same time, the temperature difference between the vapor inside the bubble and the interface diminishes. With a lower level of forcing fraction, heat transfer becomes the process that limits bubble growth.

where the radial velocity, u, is obtained by eq. (3). The initial and the boundary conditions for eq. (12) are

The energy balance at the liquid-vapor interface is

The instantaneous bubble radius is

where CR is a constant. For a large Jakob number (corresponding to low pressure), CR can be obtained by (Plesset and Zwick, 1954)

and the Jakob number is defined as
![Ja = \frac{{{\rho _\ell }{c_{p\ell }}\left[ {{T_\ell }(\infty ) - {T_{sat}}({p_v})} \right]}}{{{\rho _v}{h_{\ell v}}}} \qquad \qquad(17)](/encyclopedia/images/math/2/a/1/2a1d4fdf0faabb06dd1a69a4ec3ea107.png)
For a small Jakob number (corresponding to high pressure)

Equation (15) indicates that the bubble radius is proportional to the square root of time. Bubble growth in the two limiting cases described above has been analyzed. The inertia-controlled bubble growth model ignored heat transfer, which is only valid in the early stage of bubble growth. In the heat transfer-controlled bubble growth model, on the other hand, the effect of inertia on bubble growth was neglected; this assumption is appropriate only during the late stage of bubble growth. A comprehensive bubble growth model must include both inertia and heat transfer and provide a smooth transition between the two regimes. For bubble growth in an initially uniformly superheated liquid, Mikic et al. (1970) obtained the following expression of the relationship between bubble radius and time:
![{R^ + } = \frac{2}{3}\left[ {{{({t^ + } + 1)}^{\frac{3}{2}}} - {{({t^ + })}^{\frac{3}{2}}} - 1} \right] \qquad \qquad(19)](/encyclopedia/images/math/0/a/0/0a0b3e009f3ff650119606dd9ce30857.png)
where R + and t + are nondimensional radius and time, respectively, defined as

And
![A = {\left\{ {\frac{{b\left[ {{T_\ell }(\infty ) - {T_{sat}}({p_\ell }(\infty ))} \right]{h_{\ell v}}{\rho _v}}}{{{\rho _\ell }{T_{sat}}({p_\ell }(\infty ))}}} \right\}^{\frac{1}{2}}} \qquad \qquad(21)](/encyclopedia/images/math/2/a/8/2a8d54d8ebdf463efb4125005cf2ce24.png)

![{\rm{Ja}} = \frac{{{\rho _\ell }{c_{p\ell }}[{T_\ell }(\infty ) - {T_{sat}}({p_\ell }(\infty ))]}}{{{\rho _v}{h_{\ell v}}}} \qquad \qquad(23)](/encyclopedia/images/math/3/9/8/3986ce59d6accf9357515df985fd8919.png)
where b in eq. (21) is a geometric parameter with a value of 2/3 for a perfect spherical bubble. Equation (19) provides a generalized solution for the bubble radius because it can be reduced to eq. (11) for small t + and to eq. (15) for larger t + . Prosperetti and Plesset (1978) developed a more advanced model for vapor bubble growth in a superheated liquid.
References
Faghri, A., and Zhang, Y., 2006, Transport Phenomena in Multiphase Systems, Elsevier, Burlington, MA
Faghri, A., Zhang, Y., and Howell, J. R., 2010, Advanced Heat and Mass Transfer, Global Digital Press, Columbia, MO.
Mikic, B.B., Rohsenow, W.M., and Griffith, P., 1970, “On Bubble Growth Rate,” International Journal of Heat and Mass Transfer, Vol. 13, pp. 657-666.
Plesset, M.S., and Zwick, S.A., 1954, “The Growth of Vapor Bubble in Superheated Liquids,” Journal of Applied Physics, Vol. 25, pp. 493-500.
Prosperetti, A., and Plesset, M.S., 1978, “Vapor Bubble Growth in a Superheated Liquid,” Journal of Fluid Mechanics, Vol. 85, pp. 349-368.