research36 min studyLesson 308 in Physics

Classical and neoclassical plasma transport

Compare collision-induced steps along and across a magnetic field, then derive the geometric enhancement and ordering of banana-regime transport.

On this page

What you’ll understand

  • Estimate classical parallel and perpendicular diffusivities with their validity conditions.
  • Explain how toroidal trapped orbits modify collisional transport.
  • Distinguish neoclassical geometry effects from turbulent transport and nonlocal closure failure.

Strong magnetization makes transport anisotropic because collisions interrupt very different orbits along and across the field. In a straight field, the parallel collisional step can be a mean free path, while the perpendicular step is associated with a gyroradius. Toroidal geometry adds trapped guiding-center orbits whose radial width exceeds that radius. Neoclassical transport combines these orbit effects with collisions. It is distinct from transport driven by fluctuating turbulent electric and magnetic fields.

A diffusion estimate needs an orbit and a decorrelation time

A random-walk estimate multiplies squared step length by decorrelation rate, up to coefficients determined by the detailed kinetic problem. Along a field, a thermal particle travels approximately one mean free path between scattering events, producing an inverse-collision-frequency diffusivity. Across a strong uniform field, collisions shift the guiding center by an orbit-scale distance and give a diffusivity proportional to collision frequency. More collisions can therefore reduce parallel transport while increasing perpendicular transport within the same ordered regime.

These estimates do not establish a valid local transport closure by themselves. A classical parallel Fourier-like closure needs the mean free path small compared with the parallel gradient length, with sufficiently slow evolution and a near-Maxwellian state. Strong magnetization requires cyclotron frequency much greater than collision frequency, and perpendicular gradients should be slow across an orbit. If the mean free path exceeds the system, extending the inverse-frequency formula indefinitely predicts an unphysical local diffusion process where streaming and boundaries actually dominate.

λmfp∼vTν,ρT=vT∣Ω∣,D∥∼vT2ν,D⊥,cl∼νρT2\lambda_{\mathrm{mfp}}\sim\frac{v_T}{\nu},\quad \rho_T=\frac{v_T}{|\Omega|},\quad D_\parallel\sim\frac{v_T^2}{\nu},\quad D_{\perp,\mathrm{cl}}\sim\nu\rho_T^2
vT uses the one-component thermal convention, ν is a relevant scattering rate in s⁻¹, lengths are in metres, and diffusivities in m²/s. These are scalings, not universal numerical coefficients. Require |Ω| ≫ ν and appropriate local-gradient orderings.↔ Scroll wide equations horizontally.

Toroidal trapping creates wider radial steps

In a large-aspect-ratio axisymmetric tokamak, field strength varies along a field line and traps a fraction of particles between mirror points. Their guiding centers form banana-shaped orbits in a poloidal cross-section. The characteristic radial width contains the safety factor and inverse aspect ratio, and can greatly exceed a gyroradius while remaining small compared with equilibrium profile scales. Collisions scatter particles between trapped and passing regions, decorrelating those wider orbits.

A useful banana-regime scaling combines trapped fraction, squared banana width, and the enhanced rate at which pitch-angle scattering crosses the narrow trapped region. Multiplying these factors gives the familiar geometric enhancement of perpendicular diffusion. It is an order estimate for a specified tokamak regime, not an exact coefficient for every species, field geometry, or radial electric field. Stellarators and other nonaxisymmetric systems can have qualitatively different orbit losses and ambipolarity constraints.

ft∼ε,Δb∼qρTε,νeff∼νε,Db∼q2ε3/2νρT2f_t\sim\sqrt\varepsilon,\quad \Delta_b\sim\frac{q\rho_T}{\sqrt\varepsilon},\quad \nu_{\mathrm{eff}}\sim\frac\nu\varepsilon,\quad D_b\sim\frac{q^2}{\varepsilon^{3/2}}\nu\rho_T^2
ε = r/R ≪ 1 is local inverse aspect ratio, q is the dimensionless safety factor rather than particle charge, Δb is banana width in metres, and Db is a banana-regime diffusion scaling. Order-unity coefficients and coupled transport coefficients are omitted.↔ Scroll wide equations horizontally.

Collisionality and spatial locality select a regime

The banana regime requires scattering out of the trapped region to be slow compared with bounce motion. A dimensionless collisionality compares the effective scattering rate with the bounce frequency. Increasing it moves through intermediate plateau behavior and toward a more collisional regime; the low-collisionality scaling should not be extrapolated across all of them. Even at low collisionality, a radially local neoclassical calculation needs orbit width small relative to the profile scale.

A collisionless closed orbit does not by itself produce a secular random walk across flux surfaces. Collisions, fluctuating fields, asymmetries, or boundaries supply the mechanism that breaks repeated orbit return. Turbulent fluxes can exceed a classical or neoclassical estimate and depend on correlations not present in an equilibrium orbit calculation. Comparing measured transport with a baseline is useful only after checking geometry, collisionality, orbit width, gradients, and the definition of the measured particle or heat flux.

ωb∼εvTqR,ν∗≡νqRε3/2vT≪1,ΔbLprofile≪1\omega_b\sim\frac{\sqrt\varepsilon v_T}{qR},\qquad \nu_*\equiv\frac{\nu qR}{\varepsilon^{3/2}v_T}\ll1,\qquad \frac{\Delta_b}{L_{\mathrm{profile}}}\ll1
ωb is a bounce-frequency scale in s⁻¹, R and profile length are in metres, and ν* is dimensionless. These are banana and radial-locality conditions in the large-aspect-ratio tokamak illustration. They do not imply a collisional local parallel closure.↔ Scroll wide equations horizontally.

Classical anisotropy with explicit length checks

Take vT = 10⁵ m/s, ν = 10³ s⁻¹, and |Ω| = 10⁷ rad/s. Estimate gyroradius, mean free path, and the two classical diffusion scalings. A parallel profile varies over 10³ m.

  1. The gyroradius is ρT=105/107=0.010 m\rho_T=10^5/10^7=0.010\,\mathrm m, while mean free path is λmfp=105/103=100 m\lambda_{\mathrm{mfp}}=10^5/10^3=100\,\mathrm m.

  2. The estimates are D∥∼1010/103=107 m2/sD_\parallel\sim10^{10}/10^3=10^7\,\mathrm{m^2/s} and D⊥,cl∼103(0.010)2=0.10 m2/sD_{\perp,\mathrm{cl}}\sim10^3(0.010)^2=0.10\,\mathrm{m^2/s}. Their ratio is 10⁸.

  3. Magnetization is strong since |Ω|/ν = 10⁴. The parallel mean-free-path ratio is 0.10, small enough for an indicative local estimate but not an asymptotically tiny error parameter; a precise closure needs additional care.

The result

The orbit radius is 1 cm, mean free path 100 m, and scaling estimates are 10⁷ m²/s parallel and 0.10 m²/s perpendicular.

Explore the model

Interactive lab

Interactive lab · heat

A diffusivity changes a mode’s relaxation time

Reset and set BOTH Spatial mode n to 2 and Diffusivity κ to 0.25; predict exponent t. Then Reset and set BOTH n to 2 and κ to 1; predict exponent 4t. Play each fresh run from zero. This exact periodic scalar heat mode illustrates how a diffusion coefficient affects relaxation; it has no magnetic field, orbit geometry, trapped fraction, tensor transport, or neoclassical closure.

A diffusivity changes a mode’s relaxation timeThe exact periodic heat-equation mode is u(x,t) = exp(−κn²t) sin(nx), on [0,2π]. Playback advances physical time. The initial sine is shown separately, and the unnormalized L² norm is √π exp(−κn²t). This analytic instrument does not run a finite-difference solver or test discretization stability. Current readings: Current amplitude: 1; Spatial L² norm: 1.772; Decay rate · κn²: 3.6 /s. All parameters are adjustable with the labeled controls below.Periodic heat equation mode 3 at time 0: the amplitude is 1.0π/2π3π/22π−1−0.500.51x · periodic intervalTemperature perturbation uExact mode · n = 3
Current temperature modeInitial sineCurrent amplitude
Current amplitude
1
Spatial L² norm
1.772
Decay rate · κn²
3.6 /s
t = 0 s

The exact periodic heat-equation mode is u(x,t) = exp(−κn²t) sin(nx), on [0,2π]. Playback advances physical time. The initial sine is shown separately, and the unnormalized L² norm is √π exp(−κn²t). This analytic instrument does not run a finite-difference solver or test discretization stability. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forA diffusion law predicts relaxation only after its coefficient and regime have been justified.

Check your understanding

Within the strong-field classical scaling regime, how does halving ν at fixed vT and magnetic field affect D∥ and D⊥?

Try it yourself

For a banana-regime order estimate, use q = 2, ε = 0.04, and classical perpendicular diffusivity 0.002 m²/s. Given ν* = 0.01 and orbit width much smaller than the profile scale, estimate the geometric enhancement and Db.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Evaluate ε^(3/2) as ε√ε and remember that this is a scaling comparison rather than an exact transport coefficient.

Compare with the solution

Here ε3/2=0.04(0.2)=0.008\varepsilon^{3/2}=0.04(0.2)=0.008. The enhancement factor is q2/ε3/2=4/0.008=500q^2/\varepsilon^{3/2}=4/0.008=500. Thus Db∼500(0.002)=1 m2/sD_b\sim500(0.002)=1\,\mathrm{m^2/s}. The supplied ν* and orbit-width assumptions support the stated banana and radial-locality orderings. The numerical value is only an order estimate; species-dependent coefficients and coupled particle, heat, and momentum transport are not included.

Take it with you

  • Parallel and perpendicular classical diffusion have opposite collision-frequency scalings.
  • Trapped toroidal orbits enhance collisional transport through a larger effective radial step.
  • Neoclassical, turbulent, and nonlocal transport require different physical closures.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.