---
title: "How a Tiny Ripple Rolls into a Shear Vortex Interactive Lab | Sandbox Physics"
description: "When does a ripple between two fluid layers sliding past each other start to grow?"
canonical_url: "https://sandboxphysics.com/labs/kelvin-helmholtz-instability/"
language: "en"
page_type: "lab"
alternate_languages:
  zh-CN: "https://sandboxphysics.com/zh/labs/kelvin-helmholtz-instability/"
last_modified: "2026-10-27"
---

Fluid dynamics · interface stability

# Kelvin–Helmholtz Shear-Layer Instability

**When does a ripple between two fluid layers sliding past each other start to grow?** When, for a ripple of a given wavelength, the shear feeds it faster than gravity and surface tension can pull the interface back flat, and viscosity cannot damp it out. Gravity resists through the density difference between the layers. The growth is exponential only at first; once the ripple rolls up into vortices, the motion turns nonlinear and the simple growth rate no longer applies.

Watch two counterflowing layers, the interface between them, the vorticity and the velocity vectors in two dimensions, next to a Fourier spectrum and a growth-rate plot. Five preset cases isolate the roles of density, gravity, surface tension, viscosity and the seed spectrum.

> Interactive simulation: open this page in a web browser to change the parameters and run the model.

Early growth follows $\eta_k(t)=\eta_k(0)e^{\gamma_k t}$ and the gravity–capillary Kelvin–Helmholtz dispersion relation. Once $|\eta_k|/\lambda_k\ge 0.075$, the display switches visibly to a regularized Birkhoff–Rott sheet. This late-stage model does not resolve mixing, viscous diffusion or interface breakup.

Physics tutorial

## Kelvin–Helmholtz Instability: from perturbation growth to sheet roll-up

**Background** Two incompressible layers flow in opposite directions along a horizontal interface. A small ripple with wavenumber $k$ is amplified by the velocity jump, while stable density stratification, gravity and surface tension resist the deformation. The Lab predicts the growth of each Fourier mode and visibly switches to a different model once the small-amplitude assumption fails.

**Why it matters** Cloud decks, ocean internal waves, jet edges and astrophysical plasmas can all show shear-driven billows. Looking alike does not mean one model fits them all. What you can test is the growth rate and the range where the model holds.

Start with the essentials

- **Focus question**: When does amplification from the velocity difference overcome gravity, surface tension, and viscous damping?
- **One-sentence intuition**: In the linear phase, each wavenumber evolves on its own. A positive rate $\gamma_k>0$ only tells you the perturbation grows exponentially. Once the interface rolls up, modes couple and the geometry turns nonlinear, so the linear prediction has to stop.

### Core mathematical model

#### Two-layer gravity–capillary dispersion relation

$$
\omega=k\frac{\rho_u U_u+\rho_\ell U_\ell}{\rho_u+\rho_\ell}\pm\sqrt{\frac{(\rho_\ell-\rho_u)gk+\sigma k^3}{\rho_u+\rho_\ell}-\frac{\rho_u\rho_\ell}{(\rho_u+\rho_\ell)^2}(\Delta U)^2k^2}
$$

Under the square root, the first term is the restoring effect of stable stratification and surface tension; the second is the shear drive. When the quantity under the root turns negative, the frequency gains an imaginary part and the disturbance grows exponentially.

#### Teaching corrections for thickness and damping

$$
\gamma_k=\max\!\left[0,\sqrt{\frac{\rho_u\rho_\ell}{(\rho_u+\rho_\ell)^2}(\Delta U)^2k^2\operatorname{sech}(k\delta)-\frac{(\rho_\ell-\rho_u)gk+\sigma k^3}{\rho_u+\rho_\ell}}-2\nu k^2\right]
$$

The thickness factor and the viscous term are simple low-order teaching approximations, not a solution of the viscous Orr–Sommerfeld eigenvalue problem. They show how sensitive short waves are to thickness and dissipation.

#### Amplitude growth in the linear phase

$$
\eta_k(t)=\eta_k(0)e^{\gamma_k t},\qquad \gamma_{\mathrm{fit}}=\frac{\mathrm{d}}{\mathrm{d}t}\ln|\eta_k|
$$

The growth plot fits the slope of the log amplitude and compares it with the prediction. The theory line ends at the mark where the model switches.

#### Regularized vortex-sheet phase

$$
\frac{\partial \overline{z}}{\partial t}=\frac{1}{2\pi i}\operatorname{PV}\!\int\frac{\gamma(\alpha')}{z(\alpha)-z(\alpha')}\,\mathrm{d}\alpha'
$$

The Birkhoff–Rott velocity feeds the flow induced by the sheet back into the interface. The code smooths the singularity with a periodic vortex-blob kernel of finite core radius; that radius is a numerical choice, not a resolved physical scale.

### Common difficulties

#### One growth rate does not describe the whole run

**Typical misconception** The early growth is exponential, so the same rate should predict the later vortex size, mixing and breakup.

**Better mental model** Linear theory needs the amplitude to be much smaller than the wavelength. When the amplitude-to-wavelength ratio reaches 0.075, the Lab ends the linear prediction and labels the later stage as a regularized vortex sheet.

#### A finite layer is not a sharp interface

**Typical misconception** Increasing the thickness only blurs the picture.

**Better mental model** A layer of finite thickness weakens the velocity jump that short waves feel. This Lab uses $\operatorname{sech}(k\delta)$ to show the trend; accurate values need an eigenvalue calculation for the actual velocity profile.

### Run the experiment

1. 01

   #### Measure single-mode exponential growth

   Choose Single-mode seed, pause, and press Step 30 several times. Compare the solid and dashed traces in the growth audit plot.

   **What to observe:** While the amplitude is small, the log amplitude rises almost along a straight line and the fitted rate approaches the theory. After the model switch, the two traces are no longer required to agree.

2. 02

   #### Find the stability threshold

   Choose Stable stratification, then raise the velocity jump or lower gravity while watching the fastest mode and the active phase.

   **What to observe:** When the restoring forces win, the disturbance only oscillates and dies away. Past the threshold, a growing mode appears and the phase label changes to linear growth.

3. 03

   #### Let the spectrum pick a scale

   Choose Noise spectrum, then increase the shear-layer thickness and the viscous damping one at a time while comparing the Fourier bars.

   **What to observe:** Short waves get a stronger shear drive, but thickness, surface tension and viscosity also hold them back harder. The scale that wins comes out of this competition.

Physical conclusion

Kelvin–Helmholtz instability is a competition you can measure. Shear drives growth, stratification and surface tension pull the interface back, and dissipation filters out short waves; once the interface rolls up, a separate nonlinear model takes over.

Conceptual inspiration: [FluidMotion](http://www.bugman123.com/FluidMotion/index.html). Original simulations and graphics; the model limits are described above.

Fields & matter · Starter path

## From a field to an electron microscope

Map an electric field and steer a charged particle. Then find the electron's wavelength and use electron beams to see a crystal and a surface.

[Start here → **Make an Electric Field Visible** Move the probe around a dipole and check that the field always points downhill on the potential landscape. Observe · 8 min](https://sandboxphysics.com/labs/electric-field/)

1. [1 Make an Electric Field Visible](https://sandboxphysics.com/labs/electric-field/)
2. [2 A Charge Bending through Fields](https://sandboxphysics.com/labs/charged-particle/)
3. [3 Electron Wavelength & Voltage](https://sandboxphysics.com/labs/electron-wavelength/)
4. [4 Electron Diffraction & Reciprocal Space](https://sandboxphysics.com/labs/electron-diffraction/)
5. [5 TEM Bright Field, Dark Field & SAED](https://sandboxphysics.com/labs/tem-mode-selector/)
6. [6 SEM Secondary-Electron Topography](https://sandboxphysics.com/labs/sem-topography/)

[View the full path →](https://sandboxphysics.com/topics/fields-and-matter/#starter-path)