Abstract:Magnetic confinement fusion reactors suffer severely from heat and particle losses through turbulent transport, which has inspired the construction of ever larger and more expensive reactors. Numerical simulations are vital to their design and operation, but particle collisions are too infrequent for fluid descriptions to be valid. Instead, strongly magnetised fusion plasmas are described by the gyrokinetic equations, a nonlinear integro-differential system for evolving the particle distribution functions in a five-dimensional position and velocity space, and the consequent electromagnetic field. Due to the high dimensionality, simulations of small reactor sections require hundreds of thousands of CPU hours on High Performance Computing platforms. We develop a Hankel-Hermite spectral representation for velocity space that exploits structural features of the gyrokinetic system. The representation exactly conserves discrete free energy in the absence of explicit dissipation, while our Hermite hypercollision operator captures Landau damping with few variables. Calculation of the electromagnetic fields becomes purely local, eliminating inter-processor communication in, and vastly accelerating, searches for linear instabilities. We implement these ideas in SpectroGK, an efficient parallel code. Turbulent fusion plasmas may dissipate free energy through linear phase mixing to fine scales in velocity space, as in Landau damping, or through a nonlinear cascade to fine scales in physical space, as in hydrodynamic turbulence. Using SpectroGK to study saturated electrostatic drift-kinetic turbulence, we find that the nonlinear cascade suppresses linear phase mixing at energetically-dominant scales, so the turbulence is fluid-like. We use this observation to derive Fourier-Hermite spectra for the electrostatic potential and distribution function, and confirm these spectra with simulations.

The Gaussian radial basis function method for plasma kinetic theory

Higher-order moment convergent method in weakly anisotropic plasma and the NLVFP code for solution of the 0D-2V Vlasov-Fokker-Planck equation

Theory of spatially non-symmetric kinetic equilibria for collisionless plasmas

A Parallel Low-Rank Solver for the Six-Dimensional Vlasov-Maxwell Equations

Non-linear gyrokinetic theory of magnetoplasmas

A Fourier transformation based UGKS for Vlasov–Poisson equations in cylindrical coordinates (r,θ)

Covariant formulation of spatially non-symmetric kinetic equilibria in magnetized astrophysical plasmas

Phase space eigenfunctions with applications to continuum kinetic simulations

Generalized fluid theory including non-Maxwellian kinetic effects

A reduced kinetic method for investigating non-local ion heat transport in ideal multi-species plasmas

Kinetic equilibria of relativistic collisionless plasmas in the presence of non-stationary electromagnetic fields

A mass-conserving sparse grid combination technique with biorthogonal hierarchical basis functions for kinetic simulations

Radial drift of plasma blobs in a toroidal magnetic field with fully kinetic and reduced fluid models

Vlasov methods in space physics and astrophysics

Fokker-Planck-Boltzmann Model for the Global Electron Velocity Distribution Function Combining Ohmic and Stochastic Heating

Gyrocenter Gauge Theory and Algorithm for Nonlinear Particle Simulations of Radio-Frequency Waves in Plasmas

A new covariant formalism for kinetic plasma simulations in curved spacetimes

Gyrokinetic simulations of fusion plasmas using a spectral velocity space representation

Vlasov-Fokker-Planck-Maxwell Simulations for Plasmas in Inertial Confinement Fusion

Multiphysics simulations of collisionless plasmas