Hostname: page-component-74656f9884-npj9w Total loading time: 0 Render date: 2026-11-07T10:55:36.734Z Has data issue: false hasContentIssue false

A gyrokinetic model for the plasma periphery of tokamak devices

Published online by Cambridge University Press:  16 March 2020

B. J. Frei*
Affiliation:
École Polytechnique Fédérale de Lausanne (EPFL), Swiss Plasma Center (SPC), CH-1015Lausanne, Switzerland
R. Jorge
Affiliation:
École Polytechnique Fédérale de Lausanne (EPFL), Swiss Plasma Center (SPC), CH-1015Lausanne, Switzerland Instituto de Plasmas e Fusão Nuclear, Instituto Superior Técnico, Universidade de Lisboa, 1049-001Lisboa, Portugal
P. Ricci
Affiliation:
École Polytechnique Fédérale de Lausanne (EPFL), Swiss Plasma Center (SPC), CH-1015Lausanne, Switzerland
*
†Email address for correspondence: baptiste.frei@epfl.ch
Rights & Permissions [Opens in a new window]

Abstract

A gyrokinetic model is presented that can properly describe large and small amplitude electromagnetic fluctuations occurring on scale lengths ranging from the electron Larmor radius to the equilibrium perpendicular pressure gradient scale length, and the arbitrarily large deviations from thermal equilibrium that are present in the plasma periphery of tokamak devices. The formulation of the gyrokinetic model is based on a second-order accurate description of the single charged particle dynamics, derived from Lie perturbation theory, where the fast particle gyromotion is decoupled from the slow drifts assuming that the ratio of the ion sound Larmor radius to the perpendicular equilibrium pressure scale length is small. The collective behaviour of the plasma is obtained by a gyrokinetic Boltzmann equation that describes the evolution of the gyroaveraged distribution function. The collisional effects are included by a nonlinear gyrokinetic Dougherty collision operator. The gyrokinetic model is then developed into a set of coupled fluid equations referred to as the gyrokinetic moment hierarchy. To obtain this hierarchy, the gyroaveraged distribution function is expanded onto a Hermite–Laguerre velocity-space polynomial basis. Then, the gyrokinetic equation is projected onto the same basis yielding the spatial and temporal evolution of the Hermite–Laguerre expansion coefficients. A closed set of fluid equations for the lowest-order coefficients is presented. The Hermite–Laguerre projection is performed accurately at arbitrary perpendicular wavenumber values. Finally, the self-consistent evolution of the electromagnetic fields is described by a set of gyrokinetic Maxwell equations derived from a variational principle where the velocity integrals are explicitly evaluated.

Information

Type
Research Article
Copyright
© Cambridge University Press 2020
Figure 0

Table 1. Comparison between the gyrokinetic model presented herein and previous gyrokinetic theories. The present model is fully electromagnetic (EM), is in both drift-kinetic (DK) and gyrokinetic (GK) regimes, second-order accurate in the collisionless part (CLO), in both $\unicode[STIX]{x1D716}$ and $\unicode[STIX]{x1D716}_{\unicode[STIX]{x1D6FF}}$, and includes the effects of strong flows (SF). The exact structure the phase-space volume element (PSV) is preserved. A closed and analytical form of the gyroaverage operator, $\langle \cdot \rangle$, in terms of gyro-moments (G-M) is provided, while a FLR model, $\unicode[STIX]{x1D6E4}$, is often used in other models, or a closed expression is not given. A gyrokinetic Dougherty (D) collision operator (CO) is used, while a Coulomb (C) operator will be developed in future work. The polarization and magnetization corrections in the gyrokinetic Poisson equations (GKP) and Ampere laws (GKA) are analytically evaluated at arbitrary wavelengths and as functions of G-M. In other theories, these terms appear as velocity integrals ($\int$), or are $O(\unicode[STIX]{x1D716}_{\bot }^{2})$ accurate.