No Arabic abstract
The quasilinear particle flux arising from gyrokinetic instabilities is calculated in the electrostatic and collisionless approximation, keeping the geometry of the magnetic field arbitrary. In particular, the flux of electrons and heavy impurity ions is studied in the limit where the former move quickly, and the latter slowly, along the field compared with the mode frequency. Conclusions are drawn about how the particle fluxes of these species depend on the magnetic-field geometry, mode structure and frequency of the instability. Under some conditions, such as everywhere favourable or unfavourable magnetic curvature and modest temperature gradients, it is possible to make general statements about the fluxes independently of the details of the instability. In quasi-isodynamic stellarators with favourable bounce-averaged curvature for most particles, the particle flux is always outward if the temperature gradient is not too large, suggesting that it might be difficult to fuel such devices with gas puffing from the wall. In devices with predominantly unfavourable magnetic curvature, the particle flux can be inward, resulting in spontaneous density peaking in the centre of the plasma. In the limit of highly charged impurities, ordinary diffusion (proportional to the density gradient) dominates over other transport channels and the diffusion coefficient becomes independent of mass and charge. An estimate for the level of transport caused by magnetic-field fluctuations arising from ion-temperature-gradient instabilities is also given and is shown to be small compared with the electrostatic component.
In order to predict and analyze turbulent transport in tokamaks, it is important to model transport that arises from microinstabilities. For this task, quasilinear codes have been developed that seek to calculate particle, angular momentum, and heat fluxes both quickly and accurately. In this tutorial, we present a derivation of one such code known as QuaLiKiz, a quasilinear gyrokinetic transport code. The goal of this derivation is to provide a self-contained and complete description of the underlying physics and mathematics of QuaLiKiz from first principles. This work serves both as a comprehensive overview of QuaLiKiz specifically as well as an illustration for deriving quasilinear models in general.
A family of rigorous upper bounds on the growth rate of local gyrokinetic instabilities in magnetized plasmas is derived from the evolution equation for the Helmholtz free energy. These bounds hold for both electrostatic and electromagnetic instabilities, regardless of the number of particle species, their collision frequency, and the geometry of the magnetic field. A large number of results that have earlier been derived in special cases and observed in numerical simulations are thus brought into a unifying framework. These bounds apply not only to linear instabilities but also imply an upper limit to the nonlinear growth of the free energy.
Self-consistent simulations of neoclassical and electrostatic turbulent transport in a DIII-D H-mode edge plasma under resonant magnetic perturbations (RMPs) have been performed using the global total-f gyrokinetic particle-in-cell code XGC, in order to study density-pump out and electron heat confinement. The RMP field is imported from the extended magneto-hydrodynamics (MHD) code M3D-C1, taking into account the linear two-fluid plasma response. With both neoclassical and turbulence physics considered together, the XGC simulation reproduces two key features of experimentally observed edge transport under RMPs: increased radial particle transport in the pedestal region that is sufficient to account for the experimental pump-out rate, and suppression of the electron heat flux in the steepest part of the edge pedestal. In the simulation, the density fluctuation amplitude of modes moving in the electron diamagnetic direction increases due to interaction with RMPs in the pedestal shoulder and outward, while the electron temperature fluctuation amplitude decreases.
Global electromagnetic gyrokinetic simulations show the existence of near threshold conditions for both a high-$n$ kinetic ballooning mode (KBM) and an intermediate-$n$ kinetic version of peeling-ballooning mode (KPBM) in the edge pedestal of two DIII-D H-mode discharges. When the magnetic shear is reduced in a narrow region of steep pressure gradient, the KPBM is significantly stabilized, while the KBM is weakly destabilized and hence becomes the most-unstable mode. Collisions decrease the KBMs critical $beta$ and increase the growth rate.
The nonlinear gyrokinetic code GS2 has been extended to treat non-axisymmetric stellarator geometry. Electromagnetic perturbations and multiple trapped particle regions are allowed. Here, linear, collisionless, electrostatic simulations of the quasi-axisymmetric, three-field period National Compact Stellarator Experiment (NCSX) design QAS3-C82 have been successfully benchmarked against the eigenvalue code FULL. Quantitatively, the linear stability calculations of GS2 and FULL agree to within ~10%.