Question 1

What is the Kelvin-Helmholtz instability in a plasma?

Accepted Answer

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.

Question 2

What is the growth rate formula γ = k·ΔV/2?

Accepted Answer

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.

Question 3

Where does plasma KH instability occur in nature?

Accepted Answer

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.

Question 4

How does density ratio affect KH instability?

Accepted Answer

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 (ρ₂/ρ₁ ≫ 1) suppresses the instability: the heavier layer barely moves while the light layer deforms around it.

Question 5

What is vortex rollup and how does it develop?

Accepted Answer

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 across the shear layer.

Question 6

How does the perturbation mode k affect the simulation?

Accepted Answer

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.

Question 7

What stabilises the Kelvin-Helmholtz instability in a magnetised plasma?

Accepted Answer

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, i.e., when the magnetic energy exceeds the kinetic energy of the shear. A perpendicular magnetic field has no stabilising effect because the field lines can slip past each other, but it adds compressibility effects.

Question 8

What is the mixing layer growth rate in the nonlinear phase?

Accepted Answer

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, a result confirmed by both experiments and direct numerical simulations.

Question 9

How is plasma KH instability relevant to fusion energy?

Accepted Answer

In tokamak edge plasmas, the KH instability drives turbulent transport across the separatrix, mixing hot core plasma with cooler edge plasma. Edge-localised modes (ELMs) involve KH-like instabilities at the pedestal boundary. 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.

Question 10

What numerical method does this simulation use?

Accepted Answer

This simulation uses a 2D vortex sheet method on a periodic domain. The interface between the two plasma layers is represented by N marker points, each carrying vorticity proportional to the local velocity jump. The positions evolve via the Birkhoff-Rott integral, discretised using the Biot-Savart law with a desingularisation parameter (vortex blob regularisation). A fourth-order Runge-Kutta time integrator advances the sheet, capturing the rollup to the saturation phase.

⚡ Plasma Instability — Kelvin-Helmholtz

About the Simulation

What did you find?