Extended electron tails in electrostatic microinstabilities and the nonadiabatic response of passing electrons
Plasma Physics and Controlled Fusion IOP Publishing 64:5 (2022) 055004
Abstract:
Ion-gyroradius-scale microinstabilities typically have a frequency comparable to the ion transit frequency. Due to the small electron-to-ion mass ratio and the large electron transit frequency, it is conventionally assumed that passing electrons respond adiabatically in ion-gyroradius-scale modes. However, in gyrokinetic simulations of ion-gyroradius-scale modes in axisymmetric toroidal magnetic fields, the nonadiabatic response of passing electrons can drive the mode, and generate fluctuations in narrow radial layers, which may have consequences for turbulent transport in a variety of circumstances. In flux tube simulations, in the ballooning representation, these instabilities reveal themselves as modes with extended tails. The small electron-to-ion mass ratio limit of linear gyrokinetics for electrostatic instabilities is presented, in axisymmetric toroidal magnetic geometry, including the nonadiabatic response of passing electrons and associated narrow radial layers. This theory reveals the existence of ion-gyroradius-scale modes driven solely by the nonadiabatic passing electron response, and recovers the usual ion-gyroradius-scale modes driven by the response of ions and trapped electrons, where the nonadiabatic response of passing electrons is small. The collisionless and collisional limits of the theory are considered, demonstrating parallels in structure and physical processes to neoclassical transport theory. By examining initial-value simulations of the fastest-growing eigenmodes, the predictions for mass-ratio scaling are tested and verified numerically for a range of collision frequencies. Insight from the small electron-to-ion mass ratio theory may lead to a computationally efficient treatment of extended modes.Input files for "A novel approach to radially global gyrokinetic simulation using the flux-tube code stella": arXiv 2201.01506
University of Oxford (2022)
Abstract:
A collection of input files, and instructions on how to use them, for the manuscript entitled 'A novel approach to radially global gyrokinetic simulation using the flux-tube code stella'. These input files, which are designed for the gyrokinetic flux-tube code stella, are specifically for the simulations resulting in figures 2-10.Supplementary data for "extended electron tails in electrostatic microinstabilities and the nonadiabatic response of passing electrons": arXiv 2108.02822
University of Oxford (2022)
Abstract:
Supplementary data for the article "Extended electron tails in electrostatic microinstabilities and the nonadiabatic response of passing electrons": arXiv 2108.02822. The dataset includes a readme, GS2 FORTRAN namelist input files necessary to reproduce the simulations presented in the article, as well as scripts (using a mixture of Mathematica, MATLAB, and Python) for the calculation of collisional transport coefficients that appear in the collisional theory of the studied microinstabilities.Gyrokinetic simulations in stellarators using different computational domains
Nuclear Fusion IOP Publishing 61:11 (2021) 116074
Abstract:
In this work, we compare gyrokinetic simulations in stellarators using different computational domains, namely, flux tube (FT), full-flux-surface (FFS), and radially global (RG) domains. Two problems are studied: the linear relaxation of zonal flows (ZFs) and the linear stability of ion temperature gradient (ITG) modes. Simulations are carried out with the codes EUTERPE, GENE, GENE-3D, and stella in magnetic configurations of LHD and W7-X using adiabatic electrons. The ZF relaxation properties obtained in different FTs are found to differ with each other and with the RG result, except for sufficiently long FTs, in general. The FT length required for convergence is configuration-dependent. Similarly, for ITG instabilities, different FTs provide different results, but the discrepancy between them diminishes with increasing FT length. FFS and FT simulations show good agreement in the calculation of the growth rate and frequency of the most unstable modes in LHD, while for W7-X differences in the growth rates are found between the FT and the FFS domains. RG simulations provide results close to the FFS ones. The radial scale of unstable ITG modes is studied in global and FT simulations finding that in W7-X, the radial scale of the most unstable modes depends on the binormal wavenumber, while in LHD no clear dependency is found.Turbulent transport of impurities in 3D devices
Nuclear Fusion IOP Publishing 61:11 (2021) 116019