Abstract:This work reports on the development and numerical implementation of the linear electromagnetic gyrokinetic (GK) model in a tokamak flux-tube geometry using a moment approach based on the expansion of the perturbed distribution function on a velocity-space Hermite-Laguerre polynomials basis. A hierarchy of equations of the expansion coefficients, referred to as the gyro-moments (GM), is derived. We verify the numerical implementation of the GM hierarchy in the collisionless limit by performing a comparison with the continuum GK code GENE, recovering the linear properties of the ion-temperature gradient, trapped electron, kinetic ballooning, and microtearing modes, as well as the collisionless damping of zonal flows. The present investigation reveals the ability of the GM approach to describe fine velocity-space scale structures appearing near the trapped and passing boundary and kinetic effects associated with parallel and perpendicular particle drifts. In addition, the effects of collisions are studied using advanced collision operators, including the GK Coulomb collision operator. The main findings are that the number of GMs necessary for convergence decreases with plasma collisionality and is lower for pressure gradient-driven modes, such as in H-mode pedestal regions, compared to instabilities driven by trapped particles and magnetic gradient drifts often found in the core. The accuracy of approximations often used to model collisions (relative to the GK Coulomb operator) is studied, showing differences between collision operator models that increase with collisionality and electron temperature gradient in the case of TEM. The present linear analysis demonstrates that the GM approach efficiently describes the plasma dynamics for typical parameters of the tokamak boundary, ranging from the low-collisionality banana H-mode to the high-collisionality Pfirsch-Schlüter conditions.

Computing MHD equilibria of stellarators with a flexible coordinate frame

Constructing stellarators with quasisymmetry to high order

High-order Discretization of a Gyrokinetic Vlasov Model in Edge Plasma Geometry

Asymptotic quasisymmetric high-beta 3D MHD equilibria near axisymmetry

Energetic Particle Tracing in Optimized Quasisymmetric Stellarator Equilibria

A family of quasi-axisymmetric stellarators with varied rotational transform

A framework for discrete optimization of stellarator coils

Moment-Based Approach to the Flux-Tube linear Gyrokinetic Model

Asymptotic quasisymmetric high-beta three-dimensional magnetohydrodynamic equilibria near axisymmetry

Optimization of Nonlinear Turbulence in Stellarators

Modeling Non-planar Coils in a Full-Scale Stellarator

Going Beyond the MHD Approximation: Physics-Based Numerical Solution of the CGL Equations

An asymptotic Grad-Shafranov equation for quasisymmetric stellarators

Energetic particle tracing in optimized quasi-symmetric stellarator equilibria

A well-balanced and exactly divergence-free staggered semi-implicit hybrid finite volume/finite element scheme for the incompressible MHD equations

Stellarator equilibrium axis-expansion to all orders in distance from the axis for arbitrary plasma beta

An adaptive Newton-based free-boundary Grad-Shafranov solver

Technologies for supporting high-order geodesic mesh frameworks for computational astrophysics and space sciences

Optimal shape of stellarators for magnetic confinement fusion