Prof Michael Barnes

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.

Abstract:

The evidence of a large diffusive turbulent contribution to the radial impurity transport in Wendelstein 7-X (W7-X) plasmas has been experimentally inferred during the first campaigns and numerically confirmed by means of gyrokinetic simulations with the code stella. In general, the absence of strong impurity accumulation during the initial W7-X campaigns is attributed to this diffusive term. Given the large variety of possible stellarator configurations, in the present work the diffusive contribution is also calculated in other stellarator plasmas. In particular, a numerical cross-device comparison is presented, where the diffusion (D) and convection (V) coefficients of carbon and iron impurities produced by ion-temperature-gradient (ITG) turbulence are obtained. The simulations have been performed for the helias W7-X, the heliotron LHD, the heliac TJ-II and the quasi-axisymmetric stellarator NCSX at the radial position r/a = 0.75. The results show that, although the size of D and V can differ across the four devices, inward convection is found for all of them. For W7-X, TJ-II and NCSX the two coefficients are comparable and the turbulent peaking factor is surprisingly similar. In LHD, appreciably weaker diffusive and convective impurity transport and significantly larger turbulent peaking factor, in comparison with the other three stellarators, are predicted. All this suggests that ITG turbulence, although not strongly, would lead to negative impurity density gradients in stellarators. Then, considering mixed ITG/trapped electron mode (TEM) turbulence for the specific case of W7-X, it has been quantitatively assessed to what degree pellet fueled reduced turbulence scenarios feature reduced turbulent transport of impurities as well. The results for trace iron impurities show that, although their turbulent transport is not entirely suppressed, a significant reduction of the convection and a stronger decrease of the diffusion term are found. Although the diffusion is still above neoclassical levels, the neoclassical convection would gain under such conditions a greater specific weight on the dynamics of impurities in comparison with standard ECRH scenarios without pellet fueling.