Part 4, Chapter 7

Drift Waves

Density gradient-driven micro-instabilities in inhomogeneous plasmas

7.1 Physical Mechanism

Drift waves are low-frequency oscillations that arise whenever a plasma has a density gradient perpendicular to a magnetic field. They are ubiquitous in magnetically confined plasmas and are often called "universal instabilities" because the free energy source -- the density gradient -- is always present in confined plasmas.

Consider a slab geometry with a uniform magnetic field $\mathbf{B} = B\hat{z}$ and a density gradient in the $x$-direction:$dn_0/dx < 0$. We define the density gradient scale length:

$$L_n = -\left(\frac{1}{n_0}\frac{dn_0}{dx}\right)^{-1}$$

The mechanism involves a feedback loop: a density perturbation creates a potential perturbation (via the Boltzmann response of electrons along field lines), which drives$\mathbf{E}\times\mathbf{B}$ drifts that advect the background density gradient, reinforcing the original perturbation.

The characteristic frequency is the diamagnetic drift frequency:

$$\boxed{\omega_* = \frac{k_y \rho_s c_s}{L_n} = k_y v_{de}}$$

where $\rho_s = c_s/\Omega_i$ is the ion sound Larmor radius, $c_s = \sqrt{T_e/m_i}$ is the sound speed, and $v_{de} = \rho_s c_s / L_n$ is the electron diamagnetic drift velocity.

7.2 Drift Wave Dispersion Relation

Starting from the electron Boltzmann response and the ion vorticity equation in a slab, the linearized drift wave dispersion relation is:

$$\omega = \frac{\omega_*}{1 + k_\perp^2 \rho_s^2}$$

For long wavelengths ($k_\perp \rho_s \ll 1$), the frequency approaches $\omega_*$. For short wavelengths, the frequency decreases -- this is because finite Larmor radius effects average out the perturbation.

The phase velocity in the $y$-direction is:

$$v_{ph,y} = \frac{\omega}{k_y} = \frac{v_{de}}{1 + k_\perp^2 \rho_s^2}$$

Drift waves propagate in the electron diamagnetic drift direction, which is perpendicular to both the magnetic field and the density gradient.

7.3 Hasegawa-Mima Equation

The nonlinear evolution of drift waves is described by the Hasegawa-Mima equation, which governs the electrostatic potential $\phi$ in normalized form:

$$\boxed{\frac{\partial}{\partial t}\left(\nabla_\perp^2 \phi - \phi\right) + v_{de}\frac{\partial \phi}{\partial y} + \left[\phi, \nabla_\perp^2 \phi\right] = 0}$$

where the Poisson bracket is $[\phi, \psi] = \hat{z}\cdot(\nabla\phi \times \nabla\psi)$, representing the $\mathbf{E}\times\mathbf{B}$ nonlinearity.

This equation is the plasma analogue of the quasi-geostrophic equation in geophysical fluid dynamics, with the Coriolis parameter replaced by the density gradient. Key properties:

Conserves two invariants: energy $E = \frac{1}{2}\int(|\nabla\phi|^2 + \phi^2)\,d^2x$ and enstrophy
Supports dual cascade: energy to large scales, enstrophy to small scales
Forms zonal flows (structures with $k_y = 0$) that regulate transport
Is the simplest model capturing drift-wave turbulence physics

7.4 Computational Exploration

The code below visualizes the drift wave structure in slab geometry, showing the density perturbation pattern and the relationship between potential and density.

Drift Waves Simulation

Python

script.py60 lines

import numpy as np
import matplotlib.pyplot as plt

# Slab geometry drift wave visualization
Ly = 20.0    # box size in y (drift direction)
Lx = 10.0    # box size in x (gradient direction)
Nx, Ny = 200, 400

x = np.linspace(-Lx/2, Lx/2, Nx)
y = np.linspace(0, Ly, Ny)
X, Y = np.meshgrid(x, y, indexing='ij')

# Parameters
Ln = 3.0      # density gradient scale length
rho_s = 0.5   # ion sound Larmor radius
ky = 2 * np.pi / Ly * 3   # wavenumber in y
kx = np.pi / Lx            # wavenumber in x

# Background density profile
n0 = np.exp(-X / Ln)

# Drift wave perturbation (propagating in +y direction)
# Potential perturbation
phi_amp = 0.15
phi1 = phi_amp * np.sin(ky * Y) * np.cos(kx * X)

# Density perturbation (Boltzmann response: n1/n0 = e*phi/Te)
n1 = n0 * phi1

# Total density
n_total = n0 * (1 + phi1)

fig, axes = plt.subplots(1, 3, figsize=(15, 5))

# Background density
im0 = axes[0].pcolormesh(Y, X, n0, cmap='inferno', shading='auto')
axes[0].set_xlabel('y (drift direction)', fontsize=11)
axes[0].set_ylabel('x (gradient direction)', fontsize=11)
axes[0].set_title('Background Density $n_0(x)$', fontsize=12)
plt.colorbar(im0, ax=axes[0], shrink=0.8)

# Potential perturbation
im1 = axes[1].pcolormesh(Y, X, phi1, cmap='RdBu_r', shading='auto', vmin=-phi_amp, vmax=phi_amp)
axes[1].set_xlabel('y (drift direction)', fontsize=11)
axes[1].set_ylabel('x (gradient direction)', fontsize=11)
axes[1].set_title(r'Potential $\phi_1$', fontsize=12)
plt.colorbar(im1, ax=axes[1], shrink=0.8)

# Total perturbed density
im2 = axes[2].pcolormesh(Y, X, n_total, cmap='inferno', shading='auto')
axes[2].set_xlabel('y (drift direction)', fontsize=11)
axes[2].set_ylabel('x (gradient direction)', fontsize=11)
axes[2].set_title('Total Density $n_0 + n_1$', fontsize=12)
plt.colorbar(im2, ax=axes[2], shrink=0.8)

plt.tight_layout()
plt.savefig('output.png', dpi=100, bbox_inches='tight', facecolor='#0a0a1a')
plt.close()
print("Drift wave slab geometry visualization complete.")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

7.5 Drift Wave Instability

The ideal drift wave described by the simple dispersion relation is marginally stable. Instability requires a mechanism that introduces a phase shift between the density and potential perturbations. Common destabilizing mechanisms include:

Resistive drift waves: Finite parallel resistivity (or collisions) prevents electrons from responding adiabatically, introducing a phase lag between $n_1$ and $\phi_1$. The growth rate scales as $\gamma \propto \nu_{ei}/k_\parallel^2$.
Kinetic drift waves: Wave-particle resonance ($\omega = k_\parallel v_\parallel$) with electrons provides Landau-type drive, which can destabilize the mode even in collisionless plasmas.
Magnetic curvature: In toroidal geometry, the magnetic field curvature and $\nabla B$ drifts provide an additional coupling that drives the interchange-like component of drift modes.

7.6 Relevance to Fusion Transport

Drift waves and their nonlinear cousins (ion temperature gradient modes, trapped electron modes, electron temperature gradient modes) are responsible for the anomalous transport that limits confinement in tokamaks and stellarators:

Transport levels: Drift-wave turbulence produces heat and particle transport 10-100 times larger than classical collisional predictions.
Zonal flow regulation: Self-generated zonal flows ($\mathbf{E}\times\mathbf{B}$ shear layers) regulate turbulence and can trigger transport barriers (H-mode, ITBs).
Gyrokinetic simulations: Modern codes like GENE, GS2, and GYRO solve the gyrokinetic equation to predict drift-wave turbulent transport from first principles, guiding the design of next-generation fusion reactors.
Critical gradient: Drift-wave instabilities exhibit a threshold -- below a critical density or temperature gradient, the modes are stable, leading to the concept of "stiff" transport profiles.

← Previous Next →