Kelvin-Helmholtz Instability — Two-Phase Shear Flow

Two fluid layers sliding past each other destabilize at their interface into the characteristic KH billows. Tune density ratio, shear, surface tension, and toggle magnetic (MHD) stabilization.

Interface Evolution & Velocity Field

The heavy lower fluid (blue) slides against the lighter upper fluid (orange). When the shear overcomes gravity, tension and magnetic stiffening, the interface curls into KH billows. Arrows show the perturbed velocity field.

Growth Rate γ(k)

Imaginary part of ω. Peaks at the most-unstable wavenumber k*; the interface preferentially amplifies that wavelength.

Dispersion ωᵣ(k)

Real part of the frequency. Where γ=0 this is a stable gravity-capillary-Alfvén wave; where γ>0 the wave is overwritten by exponential growth.

The Kelvin-Helmholtz Mechanism

When two fluid layers slide past each other, the velocity jump at the interface is a sheet of vorticity. Any small corrugation of the interface lowers the pressure on the crest (Bernoulli) and lifts the trough, so the corrugation grows — the classic KH instability. For equal densities with no gravity or tension, every wavelength is unstable and the growth rate grows with wavenumber k, so the shortest waves win until surface tension (or numerical viscosity) cuts them off. Linear theory gives the dispersion relation (ρ₁+ρ₂)ω² − 2(ρ₁U₁+ρ₂U₂)kω + (ρ₁U₁²+ρ₂U₂² − B²/μ₀)k² − (ρ₁−ρ₂)gk − σk³ = 0, whose complex roots ω = ωᵣ − iγ give a positive growth rate γ wherever the shear term dominates.

What Stabilizes the Interface

Three effects fight the instability. (1) Stable density stratification — when the heavier fluid sits below (ρ₁>ρ₂), gravity (ρ₁−ρ₂)gk acts as a restoring force on long waves; this is why the sea surface only breaks into waves above a critical wind speed. (2) Surface tension σk³ kills the shortest wavelengths, setting a high-k cutoff; this is why tiny ripples are smooth. (3) A horizontal magnetic field adds magnetic tension B²k²/μ₀ (the Alfvén effect), which can suppress KH entirely — this is why the solar-wind magnetopause can remain sharp. Increasing ΔU or the density ratio pushes back toward instability.

Where KH Billows Appear

KH instability is ubiquitous. Wind over water seeds ocean waves and sea-spray. Jupiter's banded atmosphere is a gallery of KH rolled-up vortices, including the Great Red Spot's shear boundary. The Earth's magnetopause — the boundary between the solar wind and the magnetosphere — is a magnetized shear layer whose KH waves let solar-wind plasma leak in. KH also drives mixing in ocean overflows, cloud shear layers (the classic "billow cloud"), jet-engine exhaust, and the fuel-air interface in scramjet combustion. In astrophysics it shapes supernova-remnant boundaries and accretion-disk turbulence.

How to Use

Start with Classic Shear: equal densities, pure shear — every wavelength is unstable and the billows roll up fast. Now add Gravity (ρ₂<ρ₁) to mimic a stable stratification and watch long waves get suppressed. Raise Surface tension to clip the shortest waves and see the growth curve gain a high-k cutoff. Toggle MHD Mode and increase the magnetic tension: the growth rate drops across all k, and a strong enough field fully stabilizes the interface (the magnetopause case). The Growth Rate panel highlights the most-unstable mode k*; the Dispersion panel shows the real wave frequency that survives when γ=0.