11.3 Bubble dynamics

Bubble formation and growth

Bubble dynamics describes how gas pockets nucleate, grow, deform, and collapse within liquids. This topic sits at the intersection of surface tension, pressure fields, heat transfer, and fluid inertia, and it matters for everything from boiling heat transfer to cavitation damage in turbomachinery.

The process of bubble formation and growth is governed by a complex interplay of surface tension, pressure, and heat transfer. Getting a handle on these fundamentals is essential before tackling more advanced phenomena like cavitation or acoustic bubble behavior.

Nucleation sites

Nucleation sites are preferential locations where bubbles tend to form. Nucleation can be homogeneous (occurring spontaneously within the bulk liquid) or heterogeneous (occurring at solid surfaces, impurities, or pre-existing gas pockets). In practice, heterogeneous nucleation dominates because it requires far less energy.

• Surface roughness, crevices, and hydrophobic patches act as favorable nucleation sites by trapping small gas pockets that lower the energy barrier for bubble formation.
• The presence and spatial distribution of nucleation sites strongly influence the resulting bubble size distribution and the overall rate of bubble production.

Bubble growth rate

Once a bubble nucleates, its growth rate depends on the surrounding liquid properties, temperature, and pressure conditions.

• Heat-transfer-driven growth: The bubble expands as liquid evaporates at the bubble interface. This is the dominant mechanism in boiling.
• Diffusion-driven growth: Dissolved gas molecules diffuse across the interface into the bubble, contributing to expansion. This is what happens when you open a carbonated drink.
• In the early stages of growth, the bubble radius typically follows a square-root time dependence: $R \propto t^{1/2}$. This classic scaling emerges from the diffusion-limited transport of heat or mass to the interface.

Factors affecting bubble growth

Several physical properties of the liquid control how quickly a bubble can grow:

• Viscosity: Higher viscosity slows bubble growth by exerting greater drag on the expanding interface.
• Surface tension: Acts to minimize bubble surface area and opposes growth, particularly for small bubbles where the Laplace pressure ($\Delta P = 2\sigma / R$) is large.
• Thermal conductivity: Determines how fast heat reaches the bubble interface, directly controlling the evaporation rate.

External factors also matter. Pressure fluctuations, bulk flow velocity, and the concentration of dissolved gases all influence growth dynamics.

Bubble shapes and oscillations

Bubbles don't always stay perfectly round. Their shape and oscillation behavior depend on the balance between surface tension (which favors a sphere), inertial and viscous forces, and any external disturbances. These characteristics affect mass transfer rates, acoustic signatures, and whether a bubble will break apart.

Spherical vs. non-spherical bubbles

Surface tension drives bubbles toward a spherical shape because a sphere minimizes surface area for a given volume. However, bubbles deviate from this ideal under several conditions:

• High rise velocity distorts the bubble into an oblate ellipsoid.
• Strong shear flows or bubble-bubble interactions can further deform the interface.
• Surfactants alter interfacial mobility and can change the effective shape.

Common non-spherical shapes include ellipsoidal, cap-shaped, and skirted bubbles, each associated with different size ranges and Reynolds number regimes.

Natural frequency of oscillations

Every bubble has a natural oscillation frequency that depends on its size, the surrounding liquid properties, and the ambient pressure. The Minnaert frequency gives this for a spherical bubble:

$f_0 = \frac{1}{2\pi R_0} \sqrt{\frac{3\gamma P_0}{\rho}}$

where $R_0$ is the equilibrium radius, $\gamma$ is the specific heat ratio of the gas inside the bubble, $P_0$ is the ambient pressure, and $\rho$ is the liquid density.

Notice that smaller bubbles oscillate at higher frequencies. When an external acoustic field matches $f_0$, the bubble resonates and its oscillation amplitude grows dramatically. This resonance is central to both cavitation physics and medical ultrasound applications.

Damping effects on oscillations

Bubble oscillations don't persist forever. Three main damping mechanisms dissipate energy:

• Viscous damping: Shear stresses at the bubble-liquid interface resist the oscillatory motion. This effect is strongest for small bubbles and high-viscosity liquids.
• Thermal damping: During compression the gas heats up, and during expansion it cools. Heat exchange with the surrounding liquid during each cycle removes energy from the oscillation.
• Acoustic radiation damping: The oscillating bubble acts like a tiny loudspeaker, emitting sound waves that carry energy away.

Surfactants or impurities adsorbed on the bubble surface can introduce additional damping by stiffening the interface and modifying its rheological properties.

Bubble rise and motion

A bubble in a liquid rises because the buoyancy force (from the density difference between gas and liquid) exceeds the drag force resisting its motion. The details of how fast it rises, what shape it takes, and whether its path stays straight depend on bubble size and fluid properties.

Terminal velocity of rising bubbles

A rising bubble accelerates until drag balances buoyancy, at which point it reaches terminal velocity. For small spherical bubbles at low Reynolds numbers, Stokes' law applies:

$v_t = \frac{2}{9} \frac{(\rho_l - \rho_g) g R^2}{\mu}$

where $\rho_l$ and $\rho_g$ are the liquid and gas densities, $g$ is gravitational acceleration, $R$ is the bubble radius, and $\mu$ is the liquid dynamic viscosity.

This equation tells you that terminal velocity scales with $R^2$, so doubling the bubble radius quadruples the rise speed (in the Stokes regime). For larger bubbles at higher Reynolds numbers, the bubble deforms away from a sphere and the simple Stokes result no longer holds. Empirical drag correlations become necessary.

Drag coefficient and bubble shape

The drag coefficient $C_D$ quantifies the resistance a rising bubble experiences from the surrounding fluid.

• For spherical bubbles, $C_D$ is a function of the Reynolds number. Common empirical correlations include the Schiller-Naumann and Clift-Grace-Weber models.
• For non-spherical bubbles (ellipsoidal, cap-shaped), the drag coefficient differs significantly from the spherical case and depends on the specific shape regime.

Getting the drag coefficient right is critical for accurate predictions of bubble rise velocity in computational models of multiphase flows.

Bubble trajectory and path instability

Bubbles don't always rise in a straight line.

• At low Reynolds numbers, buoyancy and drag are well-balanced, and the bubble follows a rectilinear (straight vertical) path.
• Above a critical size (roughly 1–2 mm in water), the wake behind the bubble becomes unstable. Vortex shedding causes the bubble to follow a zigzag or helical trajectory.
• In shear flows, bubbles experience lateral forces (lift) that push them across streamlines.

These path instabilities enhance mixing and dispersion in bubble columns and other multiphase systems, which directly affects mass transfer performance.

Bubble coalescence and breakup

Coalescence and breakup control the bubble size distribution in any multiphase system. Coalescence merges smaller bubbles into larger ones; breakup fragments large bubbles into smaller ones. The competition between these two processes determines the equilibrium size distribution, which in turn governs interfacial area and transfer rates.

Bubble collision and coalescence

Coalescence requires two steps: collision, then merging.

1. Collision: Two bubbles are brought together by turbulence, wake entrainment, or differential buoyancy (larger bubbles rise faster and overtake smaller ones).
2. Film drainage: A thin liquid film is trapped between the colliding bubbles. This film must thin and eventually rupture before the bubbles can merge.
3. Merging: Once the film ruptures, surface tension rapidly pulls the two bubbles into a single larger bubble.

Coalescence is promoted by high bubble concentration, low liquid viscosity (faster film drainage), and the absence of surfactants. Surfactants stabilize the thin film by creating surface tension gradients (Marangoni effects) that resist drainage.

Mechanisms of bubble breakup

Breakup occurs when disruptive forces overcome the cohesive force of surface tension. The main mechanisms are:

• Turbulent breakup: Energetic eddies comparable in size to the bubble deform and fragment it. This is the dominant breakup mechanism in stirred tanks and turbulent pipe flows.
• Shear-induced breakup: When shear stresses from the surrounding flow exceed the restoring surface tension force, the bubble elongates into a filament and pinches off.
• Interfacial instabilities: Rayleigh-Taylor instability (when a heavy fluid accelerates into a light one) or Kelvin-Helmholtz instability (from velocity differences along the interface) can distort the bubble surface enough to cause fragmentation.

Influence of fluid properties

• Viscosity: Higher liquid viscosity slows film drainage, which actually promotes coalescence. At the same time, it increases shear stresses on bubbles, which can promote breakup. The net effect depends on the flow regime.
• Surface tension: Acts as the cohesive force resisting both deformation and breakup. Lower surface tension makes bubbles easier to break apart.
• Viscosity ratio (gas-to-liquid): Influences how the bubble deforms internally during breakup and affects the daughter bubble size distribution.
• Surfactants: Modify interfacial properties in complex ways. They generally inhibit coalescence (by stabilizing films) and can either promote or inhibit breakup depending on how they alter interfacial rigidity.

Heat and mass transfer

Bubbles dramatically enhance heat and mass transfer by increasing interfacial area and promoting mixing. This is why bubbling is used in so many industrial processes, from boilers to bioreactors.

Heat transfer across bubble interface

Heat moves between a bubble and the surrounding liquid through three mechanisms:

• Conduction: Heat flows through the thin thermal boundary layer surrounding the bubble, driven by the temperature gradient.
• Convection: The bubble's motion and the associated liquid circulation thin the boundary layer and enhance transport. Faster-moving bubbles transfer heat more effectively.
• Latent heat transfer: During bubble growth, liquid evaporates at the interface, absorbing latent heat from the surrounding liquid. During collapse or condensation, latent heat is released. This mechanism can dominate in boiling systems.

Mass transfer and gas diffusion

Mass transfer across the bubble interface is driven by concentration differences of dissolved species between the bubble interior and the surrounding liquid.

• The transfer rate depends on the concentration gradient, the interfacial area, and the diffusion coefficient of the gas in the liquid.
• Bubble growth can be sustained by dissolved gas diffusing into the bubble (supersaturated conditions). Conversely, a bubble will dissolve if the surrounding liquid is undersaturated.
• Bubble motion enhances mass transfer by continuously refreshing the liquid at the interface, reducing the thickness of the concentration boundary layer.

These principles are central to the design of gas-liquid reactors, bubble columns, and wastewater aeration systems.

Bubble collapse and micromixing

Bubble collapse is a rapid, violent event that occurs when the surrounding pressure suddenly increases or the bubble enters a high-pressure region.

• During collapse, the gas inside the bubble is compressed to extreme temperatures (thousands of Kelvin in some cases) and pressures.
• The collapse can generate high-velocity liquid microjets and shock waves that intensely mix the surrounding fluid at very small scales.
• Near solid surfaces, asymmetric collapse produces a liquid jet directed toward the surface. Repeated impacts cause cavitation erosion, which is a major concern for pump impellers, ship propellers, and hydraulic valves.
• The extreme local conditions during collapse are exploited in sonochemistry, where they drive chemical reactions that would not occur under normal conditions.

Cavitation and bubble dynamics

Cavitation is the formation and subsequent collapse of vapor bubbles in a liquid due to local pressure dropping below the liquid's vapor pressure. It's distinct from boiling: boiling is driven by heating, while cavitation is driven by pressure reduction at roughly constant temperature. Cavitation can be destructive (eroding turbine blades) or useful (ultrasonic cleaning), so understanding it is essential.

Cavitation inception and bubble formation

Cavitation inception is the point at which vapor bubbles first appear in a flow.

• Hydrodynamic cavitation: Pressure drops in accelerating flows (e.g., around a propeller blade tip or through a venturi).
• Acoustic cavitation: High-intensity sound waves create alternating pressure cycles; bubbles form during the low-pressure (rarefaction) phase.
• Optic cavitation: High-energy laser pulses locally vaporize liquid, creating a bubble.

Inception depends on liquid purity, dissolved gas content, and the availability of nucleation sites. Very pure, degassed water can withstand significant tension before cavitating, while "dirty" water with abundant nuclei cavitates much more readily.

Once formed, cavitation bubbles grow rapidly by evaporation and can cluster into cavitation clouds.

Rayleigh-Plesset equation

The Rayleigh-Plesset equation is the foundational model for spherical bubble dynamics in an infinite liquid:

$\rho \left( R \ddot{R} + \frac{3}{2} \dot{R}^2 \right) = P_b - P_\infty - \frac{2\sigma}{R} - \frac{4\mu \dot{R}}{R}$

where:

• $R$ = bubble radius (dots denote time derivatives)
• $\rho$ = liquid density
• $P_b$ = pressure inside the bubble
• $P_\infty$ = far-field liquid pressure
• $\sigma$ = surface tension
• $\mu$ = liquid dynamic viscosity

The left side represents liquid inertia. On the right side, the four terms are: the driving pressure difference, the far-field pressure, the Laplace pressure from surface tension, and viscous damping. This equation captures both the growth and collapse phases of a cavitation bubble and forms the basis for nearly all computational cavitation models.

Cavitation damage and erosion

Cavitation damage results from the repeated collapse of bubbles near solid surfaces.

1. A bubble collapses asymmetrically near a wall, producing a high-speed liquid microjet directed at the surface.
2. The jet impacts the surface at velocities that can exceed 100 m/s, along with pressure pulses from the associated shock wave.
3. Repeated impacts cause fatigue, pitting, and progressive material removal.

The severity of erosion depends on bubble size, collapse intensity, material hardness, and the frequency of cavitation events. Mitigation strategies include:

• Designing flow passages to avoid low-pressure regions where cavitation initiates
• Using cavitation-resistant materials (e.g., stainless steel, stellite coatings)
• Injecting air or polymer additives to cushion bubble collapse

Acoustic and ultrasonic cavitation

Acoustic cavitation occurs when high-intensity sound waves (typically ultrasound, above 20 kHz) drive bubble formation and collapse in a liquid. The extreme local conditions generated during collapse make ultrasonic cavitation a powerful tool for cleaning, chemical processing, and biomedical applications.

Acoustic cavitation threshold

The acoustic cavitation threshold is the minimum acoustic pressure amplitude needed to nucleate bubbles in a given liquid.

• Liquids with higher vapor pressures and lower surface tensions cavitate more easily (lower threshold).
• Higher ultrasonic frequencies generally require higher pressure amplitudes to induce cavitation, because the rarefaction phase is shorter and bubbles have less time to grow.
• Pre-existing nuclei (dissolved gas, microparticles) significantly lower the threshold. Degassed, filtered liquids are much harder to cavitate.

Knowing the threshold is essential for designing ultrasonic systems: you need enough acoustic power to exceed it, but excessive power wastes energy and can cause unwanted erosion.

Ultrasonic bubble dynamics

Under ultrasonic irradiation, bubbles experience rapid expansion and compression each acoustic cycle. Two distinct regimes emerge:

• Stable cavitation: Bubbles oscillate around an equilibrium size over many cycles without collapsing. They generate steady microstreaming currents in the surrounding liquid.
• Transient (inertial) cavitation: Bubbles grow rapidly during rarefaction and then collapse violently during compression. The collapse can produce localized temperatures above 5000 K and pressures exceeding 500 atm, along with shock waves and liquid microjets.

The bubble behavior depends on the acoustic pressure amplitude, frequency, bubble size relative to the resonance size, and liquid properties. The Rayleigh-Plesset equation (with a time-varying $P_\infty$ representing the acoustic field) is the standard modeling tool.

Applications in cleaning and processing

Ultrasonic cavitation is used across many industries:

• Ultrasonic cleaning: Collapsing bubbles near contaminated surfaces generate microjets and shock waves that dislodge particles, grease, and biofilms. The mechanical scrubbing action reaches into crevices and geometries that conventional cleaning cannot access.
• Sonochemistry: The extreme temperatures and pressures inside collapsing bubbles drive chemical reactions, including the generation of free radicals. This is used for degradation of pollutants, synthesis of nanoparticles, and acceleration of organic reactions.
• Materials processing: Ultrasonic cavitation is applied in emulsification, dispersion of nanoparticles, degassing of melts, and cell disruption in biotechnology.

In each case, the key is controlling the cavitation intensity to achieve the desired effect without causing damage to the equipment or the product.