Rayleigh-Plesset Equation: The Life and Violent Death of a Bubble

What It Is and Why It Matters

The Rayleigh-Plesset equation is the fundamental equation of motion for a single spherical bubble in a liquid. It reduces the full Navier-Stokes problem — in principle a partial-differential problem for the whole fluid — to a single ordinary differential equation for one variable: the bubble radius R(t). This collapse in dimensionality is possible because of one powerful assumption: spherical symmetry in an incompressible liquid of infinite extent. Once the wall moves, incompressibility fixes the radial velocity everywhere as u(r,t) = R²Ṙ/r².

Why does it matter? Because bubbles do not merely fizz — they can implode with astonishing violence. When a low-pressure region grows a bubble that then collapses, the concentrated inertia of the inrushing water focuses energy into a shrinking volume, generating microjets, shock waves, temperatures of tens of thousands of kelvin, and pressures of thousands of atmospheres. This is cavitation, and it pits ship propellers and pump impellers, but is also harnessed in medical lithotripsy, ultrasound imaging, and sonochemistry.

The Mechanism, Step by Step

Start from mass conservation. Incompressibility forces radial flow u(r) = R²Ṙ/r², so all the fluid's kinetic energy is slaved to the wall motion. Integrate the radial momentum equation (the Navier-Stokes r-component) from the bubble wall out to infinity. The pressure boundary condition at the wall is not simply the gas pressure: surface tension adds a Laplace jump 2σ/R, and viscosity contributes a normal stress 4μṘ/R.

The kinetic energy of the entire liquid is (KE) = 2πρR³Ṙ². Its time derivative must equal the work done by the pressure difference across the wall — this is the energy route to the equation. The two inertial terms, R R̈ and (3/2)Ṙ², emerge directly from differentiating that energy. The crucial physics of collapse lives in the (3/2)Ṙ² term: as R shrinks, mass conservation demands Ṙ grow, and this term drives a runaway. Meanwhile the gas inside resists as p_B ∝ R⁻³ᵞ (near-adiabatic, γ ≈ 5/3 for a monatomic gas), providing the stiff spring that eventually reverses the implosion in a violent rebound.

The Governing Equation and Characteristic Scales

In its standard form the equation reads:

R R̈ + (3/2) Ṙ² = (1/ρ) [ p_B(R,t) − p_∞(t) − 2σ/R − 4μṘ/R ]

Here R is the radius, ρ the liquid density, p_B the pressure inside the bubble (gas + vapor), p_∞ the driving pressure far away, σ the surface tension, and μ the dynamic viscosity. Rayleigh's original 1917 result is the special case of an empty cavity (p_B = 0, σ = μ = 0) collapsing under a fixed overpressure Δp. He derived the celebrated Rayleigh collapse time:

τ = 0.915 R_max √(ρ/Δp)

For a 45-μm bubble in water (ρ = 1000 kg/m³) under Δp ≈ 1 atm ≈ 10⁵ Pa, τ ≈ 4 μs. The wall velocity Ṙ formally diverges as R → 0 in the empty-cavity limit — real bubbles cut this off when the compressed gas pressure and, ultimately, liquid compressibility (accounted for in the Keller-Miksis extension) intervene at wall speeds comparable to the ~1500 m/s sound speed in water.

How It Is Realized and Measured

The cleanest laboratory realization is single-bubble sonoluminescence (SBSL), discovered by Gaitan and Crum around 1990. A single air bubble is trapped at the pressure antinode of a standing acoustic wave in a flask of degassed water, driven at 20-30 kHz. The bubble breathes in sync with the sound: an ambient radius R₀ ≈ 4.5 μm expands to R_max ≈ 45 μm during the rarefaction half-cycle, then collapses inertially to R_min ≈ 0.5 μm.

At each collapse it emits a flash of broadband ultraviolet-to-visible light lasting less than 50 picoseconds, clock-locked to the acoustic period with femtosecond regularity. Spectra imply a gas temperature of at least 20,000 K — some models push toward 10⁴-10⁵ K. Experimentally, R(t) is tracked by Mie scattering of a laser off the bubble wall, and the light pulse by fast photomultipliers and streak cameras. The measured R(t) curve — slow growth, sharp catastrophic collapse, damped afterbounces — is quantitatively reproduced by integrating the Rayleigh-Plesset equation (with the Keller-Miksis compressible correction), one of the great successes of the theory.

Where It Operates and How It Differs from Related Effects

The Rayleigh-Plesset equation governs any near-spherical bubble whose surroundings can be treated as incompressible on the relevant timescale: cavitation bubbles behind propellers and hydrofoils, laser- and spark-induced bubbles, ultrasound contrast microbubbles, and lithotripsy cavities. Cavitation inception is set by the Blake threshold, the critical tension below which surface tension can no longer hold a nucleus in equilibrium and explosive growth ensues.

It is important to distinguish the regimes. Non-inertial (stable) cavitation means small oscillations about equilibrium — the linearized RP equation gives the Minnaert resonance frequency ω₀ = √(3γ p₀/ρ)/R₀. Inertial (transient) cavitation is the violent growth-and-collapse the full nonlinear equation predicts. The RP equation itself breaks down when the wall goes supersonic (needing the compressible Keller-Miksis or Gilmore models), when the bubble deviates from sphericity (Rayleigh-Taylor and parametric shape instabilities produce microjets near walls), or when internal chemistry and shock structure inside the gas matter.

Applications, Open Questions, and Significance

Cavitation is a double-edged sword. Its violence erodes turbine blades, pump impellers, and propellers — Rayleigh's 1917 study was itself motivated by Royal Navy propeller damage. Yet the same energy focusing is deliberately exploited: shock-wave lithotripsy shatters kidney stones, targeted microbubbles enhance ultrasound imaging and can transiently open the blood-brain barrier for drug delivery, and sonochemistry uses collapse hot spots to drive radical reactions and clean surfaces.

Open questions remain genuinely deep. The precise light-emission mechanism of sonoluminescence — bremsstrahlung, recombination, and blackbody-like emission from a partially ionized plasma — is still debated, as is whether shock waves converge inside the gas. The most provocative claim, sonofusion (bubble-collapse-driven deuterium fusion, reported by Taleyarkhan in 2002), was never independently reproduced and is not accepted. Modern work couples the RP framework to detailed gas thermodynamics, mass diffusion (the rectified-diffusion stability that selects the ambient radius), and shape-stability analysis — extending a 1917 equation into an active research frontier a century later.

Phases of a cavitation bubble and the dominant physics in each regime, as captured by the Rayleigh-Plesset equation

Phase	Dominant balance	Characteristic scale	What the RP equation predicts
Nucleation / inception	Tension vs. surface tension 2σ/R	Blake threshold p ≈ p_v − 4σ/(3√3 R₀)	Unstable equilibrium; bubble grows once tension exceeds threshold
Explosive growth	Driving underpressure vs. inertia	R_max ≈ 10× R₀ over ~half acoustic cycle	Slow, near-isothermal expansion to maximum radius
Inertial collapse	Liquid inertia (R R̈, 3/2 Ṙ²) vs. gas pressure	τ ≈ 0.915 R_max√(ρ/Δp) ≈ μs; Ṙ → supersonic	Runaway acceleration; wall speed exceeds sound speed in liquid
Rebound / minimum	Compressed gas p_B ~ R⁻³ᵞ vs. inertia	R_min ≈ 0.5 μm; T ≳ 20,000 K; flash < 50 ps	Sharp turnaround, shock launch, damped ringing afterbounces