Heavy fluid over light fluid → mushroom plumes
What is the Rayleigh-Taylor instability? It is a hydrodynamic instability that occurs whenever a denser fluid is supported by, or accelerated into, a lighter fluid. Gravity pulls the heavy fluid down through the light fluid, so any small ripple on the interface grows over time into descending spikes and rising bubbles.
Why do mushroom-shaped plumes form? As a heavy spike sinks, the velocity difference between the falling spike and the surrounding light fluid creates a secondary Kelvin-Helmholtz shear instability along its flanks. This rolls the tip of the spike outward into the familiar mushroom cap with two curling vortices.
What is the Atwood number? The Atwood number A = (ρ_heavy − ρ_light)/(ρ_heavy + ρ_light) is the dimensionless density contrast. It ranges from 0 (equal densities, neutrally stable) to 1 (heavy fluid in vacuum). A larger Atwood number gives faster, more violent plume growth.
In the early linear stage a perturbation of wavenumber k grows exponentially as exp(n·t), with growth rate n = sqrt(A·g·k). Short wavelengths grow fastest in the inviscid case, while viscosity and surface tension suppress the smallest scales and select a dominant wavelength.
It shapes supernova remnants such as the Crab Nebula, drives mixing in inertial-confinement-fusion targets, governs salt domes rising through rock, sculpts mushroom clouds, and appears any time you invert a glass of water sealed by a card.
The simulation solves a 2D incompressible two-fluid flow on a grid. A density field is advected with the velocity using a semi-Lagrangian scheme, buoyancy forces are added from gravity acting on the density, and a Jacobi pressure projection makes the velocity field divergence-free each step.
Incompressible fluid cannot be compressed, so its velocity field must have zero divergence. The projection step solves a Poisson equation for pressure and subtracts the pressure gradient from the velocity, removing any compression introduced by advection and buoyancy.
Viscosity diffuses momentum and damps velocity gradients, especially at small scales. This slows the roll-up of mushroom caps, smooths fine filaments, and lowers the growth rate of short-wavelength perturbations, producing larger, smoother plumes.
Rayleigh-Taylor is driven by buoyancy when heavy fluid sits over light fluid, while Kelvin-Helmholtz is driven by velocity shear between two streams. The two are linked: the shear along the flanks of a Rayleigh-Taylor spike triggers a secondary Kelvin-Helmholtz roll-up that forms the mushroom cap.
If the light fluid sits on top of the heavy fluid the arrangement is stable and ripples merely oscillate. The classic instability requires heavy-over-light (or equivalently acceleration directed from the light fluid into the heavy fluid). Surface tension can also stabilise sufficiently short wavelengths.
A dense fluid (purple) rests on top of a lighter fluid (dark). Because heavy-over-light is gravitationally unstable, the slightly wavy interface cannot stay flat: heavy fingers plunge downward while light bubbles rise, rolling up into mushroom-shaped plumes that mix the two fluids together.
∂u/∂t + (u·∇)u = −∇p/ρ + ν∇²u + g·ρ̂
∇·u = 0 ∂ρ/∂t + (u·∇)ρ = 0
A = (ρ_h − ρ_l)/(ρ_h + ρ_l) n = √(A·g·k)
The same instability that makes a mushroom cloud also sculpts the filaments of the Crab Nebula — and is one of the main obstacles to igniting controlled nuclear fusion, because it mixes cold fuel into the hot core just as it is meant to compress.