No Arabic abstract
Boundary plasma physics plays an important role in tokamak confinement, but is difficult to simulate in a gyrokinetic code due to the scale-inseparable nonlocal multi-physics in magnetic separatrix and open magnetic field geometry. Neutral particles are also an important part of the boundary plasma physics. In the present paper, noble electrostatic gyrokinetic techniques to simulate the flux-driven, low-beta electrostatic boundary plasma is reported. Gyrokinetic ions and drift-kinetic electrons are utilized without scale-separation between the neoclassical and turbulence dynamics. It is found that the nonlinear intermittent turbulence is a natural gyrokinetic phenomenon in the boundary plasma in the vicinity of the magnetic separatrix surface and in the scrape-off layer.
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.
A gyrokinetic model is presented that can properly describe strong flows, 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 large deviations from thermal equilibrium. The formulation of the gyrokinetic model is based on a second order 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 behavior of the plasma is obtained by a gyrokinetic Boltzmann equation that describes the evolution of the gyroaveraged distribution function and includes a non-linear 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 velocity-space Hermite-Laguerre polynomial basis and the gyrokinetic equation is projected onto the same basis, obtaining the spatial and temporal evolution of the Hermite-Laguerre expansion coefficients. 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 Maxwells equations derived from a variational principle, with the velocity integrals of the gyroaveraged distribution function explicitly evaluated.
Magnetic islands (MIs), resulting from a magnetic field reconnection, are ubiquitous structures in magnetized plasmas. In tokamak plasmas, recent researches suggested that the interaction between the MI and ambient turbulence can be important for the nonlinear MI evolution, but a lack of detailed experimental observations and analyses has prevented further understanding. Here, we provide comprehensive two-dimensional observations that indicate various effects of the ambient turbulence on the nonlinear MI evolution. It is shown that the modified plasma turbulence around the MI can lead to either destabilization or stabilization of the MI instability in tokamak plasmas. In particular, significantly enhanced turbulence at the X-point of the MI results in a violent disruption through the fast magnetic reconnection and magnetic field stochastization.
Electrostatic turbulence in weakly collisional, magnetized plasma can be interpreted as a cascade of entropy in phase space, which is proposed as a universal mechanism for dissipation of energy in magnetized plasma turbulence. When the nonlinear decorrelation time at the scale of the thermal Larmor radius is shorter than the collision time, a broad spectrum of fluctuations at sub-Larmor scales is numerically found in velocity and position space, with theoretically predicted scalings. The results are important because they identify what is probably a universal Kolmogorov-like regime for kinetic turbulence; and because any physical process that produces fluctuations of the gyrophase-independent part of the distribution function may, via the entropy cascade, result in turbulent heating at a rate that increases with the fluctuation amplitude, but is independent of the collision frequency.
3D2V continuum gyrokinetic simulations of electrostatic plasma turbulence in a straight, open-field-line geometry have been performed using the full-$f$ discontinuous-Galerkin code Gkeyll. These simulations include the basic elements of a fusion-device scrape-off layer: localized sources to model plasma outflow from the core, cross-field turbulent transport, parallel flow along magnetic field lines, and parallel losses at the limiter or divertor with sheath model boundary conditions. The set of sheath boundary conditions used in the model allows currents to flow through the walls. In addition to details of the numerical approach, results from numerical simulations of turbulence in the Large Plasma Device (LAPD), a linear device featuring straight magnetic field lines, are presented.