Kelvin-Helmholtz Instability: Shear Rolling Into Vortices

Two fluids sliding past each other

Whenever two fluid layers move past each other at different speeds — wind over water, warm air over cold, fast plasma next to slow — the interface between them is not stable. Any tiny ripple in the boundary gets amplified: on the crest of the ripple, streamlines bunch together and the flow speeds up, which by Bernoulli's principle drops the local pressure and sucks the crest further up; in the trough, the opposite happens. The result is the Kelvin-Helmholtz instability, named after Lord Kelvin and Hermann von Helmholtz who worked out the underlying fluid dynamics in the 1860s-70s, and it turns a straight shear line into a train of rolling billows that eventually curl into vortices and mix the two layers together.

The velocity shear as the energy source

The instability draws its energy directly from the kinetic energy of the shear — the relative velocity difference across the interface. For two inviscid, incompressible fluid layers with velocities U₁ and U₂ and densities ρ₁ and ρ₂, linear stability analysis shows that a sinusoidal perturbation of wavenumber k grows exponentially with growth rate:

γ = k · √( ρ₁ρ₂ (U₁−U₂)² / (ρ₁+ρ₂)² − g k (ρ₂−ρ₁) / (ρ₁+ρ₂) )

(second term drops out for equal densities and no gravity — pure shear)

With no density difference and no gravity, every wavelength is unstable — the classic textbook result. Once gravity and a density difference are added (heavier fluid on the bottom), the second term inside the root becomes a stabilising restoring force, and the instability only grows once the shear is strong enough to overcome it.

The Richardson number: when stratification wins

That competition between destabilising shear and stabilising buoyancy is captured in a single dimensionless number, the gradient Richardson number: Ri = N² / (dU/dz)², where N is the buoyancy (Brunt-Väisälä) frequency of the density stratification and dU/dz is the local shear. The Miles-Howard theorem shows that a stratified shear flow is linearly stable if Ri > 1/4 everywhere in the flow — strong enough stratification can suppress the instability outright, while a shear layer with Ri < 1/4 anywhere is a candidate for billows to form there.

From ripple to roll-up to turbulence

Linear theory only describes the initial exponential growth of a small ripple. As the amplitude grows, the interface steepens, folds over on itself, and rolls into discrete cat-eye vortices spaced roughly one unstable wavelength apart — the textbook swirling-cloud pattern you can spot in satellite photos of cloud decks over mountains, or in the banded, wavy edges of Jupiter's cloud belts. Neighbouring vortices then interact: they can merge (pairing), which doubles the characteristic size of the structures and is one route by which a laminar shear layer transitions to fully developed, three-dimensional turbulence.

Where it shows up

The same instability appears at wildly different scales because the governing physics — shear plus a density or velocity jump — is scale-free. In the atmosphere it produces the braided, billowing Kelvin-Helmholtz clouds sometimes visible at the edge of temperature inversions; in the ocean it mixes layers of different salinity and temperature at the thermocline; in astrophysics it appears at the boundary of the solar wind and Earth's magnetosphere, at the edges of jets from active galactic nuclei, and inside supernova remnants where the shocked ejecta shears past the surrounding interstellar medium. In every case, the instability is one of the primary ways nature turns an ordered, layered flow into a turbulent, mixed one.

Frequently asked questions

Try it live

Everything above runs in your browser — open Kelvin-Helmholtz Instability and change the parameters while it is running. Nothing is installed, nothing is uploaded, the whole model lives in one tab.