Gyrokinetic Solvers
Gyrokinetic solvers evolve a reduced kinetic equation that averages over fast gyration, cutting phase-space dimension while retaining turbulence physics.
Averaging away the gyration
A magnetized plasma particle spirals rapidly around field lines while drifting slowly across them. Following the full gyration in a six-dimensional kinetic simulation is prohibitively expensive for the slow turbulent transport of interest. Gyrokinetic theory averages over the fast gyromotion analytically, reducing the problem to a five-dimensional distribution function on guiding centers with the magnetic moment as a conserved parameter. Gyrokinetic solvers discretize and evolve this reduced equation.
The reduction removes the fastest time scale (the gyrofrequency) and one velocity dimension, making first-principles turbulence simulation of realistic devices feasible while keeping the kinetic effects, wave-particle resonance and finite-Larmor-radius physics, that fluid models miss.
Numerical approaches
Two broad strategies exist. Particle-in-cell gyrokinetics samples the distribution with marker particles, usually in delta-f form to cut noise, coupled to a field solve for the gyro-averaged potential. Continuum (Eulerian) gyrokinetics discretizes the distribution on a fixed phase-space grid using finite-volume, finite-difference, spectral, or semi-Lagrangian methods. Both must handle the gyro-averaging operator and the field equation that closes the system.
Coordinates and challenges
Solvers typically use field-aligned coordinates so the grid follows the strong anisotropy of turbulence, which varies slowly along field lines and rapidly across them; this alignment saves enormous resolution. Challenges include the gyro-averaging (a Bessel-function operator), maintaining conservation of energy and particles, and controlling numerical noise or dissipation over long turbulent runs.
- Fast gyration averaged out, phase space cut to five dimensions
- Particle (delta-f) versus continuum discretizations
- Field-aligned coordinates exploit turbulence anisotropy
- Gyro-averaging and conservation are key numerical hurdles
Gyrokinetic simulation is the primary first-principles tool for predicting turbulent transport in magnetic-confinement devices, informing how heat and particles leak across confining fields.