Professor Felix Parra Diaz

Abstract:

The first experimental campaigns have proven that, due to the optimization of the magnetic configuration with respect to neoclassical transport, the contribution of turbulence is essential to understand and predict the total particle and energy transport in Wendelstein 7-X (W7-X). This has spurred much work on gyrokinetic modelling for the interpretation of the available experimental results and for the preparation of the next campaigns. At the same time, new stellarator gyrokinetic codes have just been or are being developed. It is therefore desirable to have a sufficiently complete, documented and verified set of gyrokinetic simulations in W7-X geometry against which new codes or upgrades of existing codes can be tested and benchmarked. This paper attempts to provide such a set of simulations in the form of a comprehensive benchmark between the recently developed code stella and the well-established code GENE. The benchmark consists of electrostatic gyrokinetic simulations in the W7-X magnetic geometry and includes different flux tubes, linear ion-temperature-gradient (ITG) and trapped-electron-mode stability analyses, computation of linear zonal-flow responses and calculation of ITG-driven heat fluxes.

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.

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.