About the Simulation
The Kelvin-Helmholtz (KH) instability in plasma arises whenever two magnetised plasma layers flow past each other at different speeds. The velocity discontinuity at their boundary acts as a vortex sheet: any tiny perturbation is amplified by the Bernoulli pressure imbalance between the fast and slow sides, causing the interface to roll up into characteristic cat-eye vortices. The linear growth rate is γ = k·ΔV·√(ρ₁ρ₂)/(ρ₁+ρ₂), so higher wavenumber modes (shorter wavelengths) grow fastest — until nonlinear saturation and vortex merging transfer energy to larger scales.
This simulation models the plasma interface using a regularised vortex sheet (Birkhoff-Rott equation with blob desingularisation). N marker points distributed along the interface carry circulation proportional to the local velocity jump ΔV. Their positions evolve via the discretised Biot-Savart kernel, advanced in time with a Runge-Kutta integrator. The dispersion panel below the main canvas shows the theoretical growth rate γ(k) curve, with a marker at the currently selected mode. You can control the velocity shear, the density ratio between the two plasma layers, and the initial perturbation wavenumber k using the sliders in the HUD panel.
Frequently Asked Questions
What is the Kelvin-Helmholtz instability in a plasma?
The Kelvin-Helmholtz (KH) instability in a plasma occurs at the interface between two plasma layers moving at different velocities. The velocity shear at the boundary amplifies small perturbations: pressure drops where flow is faster (Bernoulli effect), pulling the interface further in, generating vortices. In magnetised plasma, the magnetic field can stabilise or destabilise depending on field orientation relative to the shear flow.
What is the growth rate formula γ = k·ΔV/2?
For two incompressible fluid layers of equal density with a velocity difference ΔV = V₁ − V₂, the KH growth rate is γ = k·ΔV/2, where k = 2π/λ is the wavenumber of the surface perturbation. This means shorter wavelengths (larger k) grow faster in the equal-density limit. When the densities differ, the growth rate becomes γ = k·ΔV·√(ρ₁ρ₂)/(ρ₁+ρ₂), reducing the growth for large density contrasts.
Where does plasma KH instability occur in nature?
Plasma KH instability is ubiquitous in space plasmas. At the Earth's magnetopause, the solar wind shearing past the magnetosphere drives KH vortices that transport solar plasma inward. Jupiter's cloud belts show KH billows at jet-stream boundaries. In the solar corona, KH vortices have been directly imaged in plumes and coronal rain. Pulsar wind nebulae and astrophysical jets show characteristic KH patterns at their boundaries.
How does density ratio affect KH instability?
The density ratio ρ₂/ρ₁ modifies the growth rate through the factor √(ρ₁ρ₂)/(ρ₁+ρ₂). At equal density (ratio = 1) this factor equals 1/2, giving the maximum growth rate k·ΔV/2. As the density contrast increases (large ρ₂/ρ₁), the effective growth rate decreases and the heavy layer resists deflection. A very large density ratio suppresses the instability.
What is vortex rollup and how does it develop?
Vortex rollup is the nonlinear saturation of the KH instability. Initially, sinusoidal perturbations grow exponentially (linear phase). As amplitudes become comparable to wavelength, vorticity concentrates into discrete cores: the interface winds up into mushroom-shaped cat-eye vortices. Adjacent vortices eventually merge (vortex pairing) to form larger structures — an inverse energy cascade that drives turbulent mixing.
How does the perturbation mode k affect the simulation?
The integer perturbation mode k sets how many complete sinusoidal wavelengths fit across the simulation domain. Mode k=1 shows a single large vortex rollup; k=4 shows four simultaneous vortices. Higher modes grow faster (γ ∝ k) in the linear phase but also saturate at smaller amplitude relative to domain size. The most visually dramatic rollup typically occurs at intermediate modes k=2–4.
What stabilises the Kelvin-Helmholtz instability in magnetised plasma?
A magnetic field component parallel to the flow provides magnetic tension that resists the bending of field lines during vortex formation. The instability is stabilised when the Alfvén speed V_A = B/√(μ₀ρ) exceeds ΔV/2. A perpendicular magnetic field has no stabilising effect because the field lines can slip past each other.
What is the mixing layer growth rate in the nonlinear phase?
After initial exponential growth, the mixing layer width δ grows linearly with time as δ = α·ΔV·t, where α ≈ 0.06–0.1 is the turbulent mixing coefficient (Brown-Roshko constant). This linear growth phase persists until the shear layer fills the domain. The mixing layer momentum thickness θ grows as dθ/dt ≈ 0.017·ΔV.
How is plasma KH instability relevant to fusion energy?
In tokamak edge plasmas, the KH instability drives turbulent transport across the separatrix, mixing hot core plasma with cooler edge plasma. In inertial confinement fusion (ICF), KH instabilities at the ablator-fuel interface degrade implosion symmetry and mix cold ablator material into the hot fuel, quenching ignition. Controlling KH instability is therefore central to achieving fusion energy.
What numerical method does this simulation use?
This simulation uses a 2D vortex sheet method on a periodic domain. The interface is represented by N marker points, each carrying vorticity proportional to the local velocity jump. Their positions evolve via the Birkhoff-Rott integral, discretised using the Biot-Savart law with a blob desingularisation parameter (δ²) for regularisation. A Runge-Kutta time integrator advances the sheet positions, capturing rollup through to the saturation phase.