No Arabic abstract
A study of turbulent impurity transport by means of quasilinear and nonlinear gyrokinetic simulations is presented for Wendelstein 7-X (W7-X). The calculations have been carried out with the recently developed gyrokinetic code stella. Different impurity species are considered in the presence of various types of background instabilities: ITG, TEM and ETG modes for the quasilinear part of the work; ITG and TEM for the nonlinear results. While the quasilinear approach allows one to draw qualitative conclusions about the sign or relative importance of the various contributions to the flux, the nonlinear simulations quantitatively determine the size of the turbulent flux and check the extent to which the quasilinear conclusions hold. Although the bulk of the nonlinear simulations are performed at trace impurity concentration, nonlinear simulations are also carried out at realistic effective charge values, in order to know to what degree the conclusions based on the simulations performed for trace impurities can be extrapolated to realistic impurity concentrations. The presented results conclude that the turbulent radial impurity transport in W7-X is mainly dominated by ordinary diffusion, which is close to that measured during the recent W7-X experimental campaigns. It is also confirmed that thermo-diffusion adds a weak inward flux contribution and that, in the absence of impurity temperature and density gradients, ITG- and TEM-driven turbulence push the impurities inwards and outwards, respectively.
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 attemps 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 W7-X magnetic geometry and includes different flux tubes, linear ion-temperature-gradient (ITG) and trapped-electron-mode (TEM)} stability analyses, computation of linear zonal flow responses and calculation of ITG-driven heat fluxes.
We study the effect of turbulent transport in different magnetic configurations of the Weldenstein 7-X stellarator. In particular, we performed direct numerical simulations with the global gyrokinetic code GENE-3D, modeling the behavior of Ion Temperature Gradient turbulence in the Standard, High-Mirror, and Low-Mirror configurations of W7-X. We found that the Low-Mirror configuration produces more transport than both the High-Mirror and the Standard configurations. By comparison with radially local simulations, we have demonstrated the importance of performing global nonlinear simulations to predict the turbulent fluxes quantitatively.
Nonlinear gyrokinetic simulations have been conducted to investigate turbulent transport in tokamak plasmas with rotational shear. At sufficiently large flow shears, linear instabilities are suppressed, but transiently growing modes drive subcritical turbulence whose amplitude increases with flow shear. This leads to a local minimum in the heat flux, indicating an optimal E x B shear value for plasma confinement. Local maxima in the momentum fluxes are also observed, allowing for the possibility of bifurcations in the E x B shear. The sensitive dependence of heat flux on temperature gradient is relaxed for large flow shear values, with the critical temperature gradient increasing at lower flow shear values. The turbulent Prandtl number is found to be largely independent of temperature and flow gradients, with a value close to unity.
In the Wendelstein 7-X magnetic confinement experiment, a reduction of turbulent density fluctuations as well as anomalous impurity diffusion is associated with a peaking of the plasma density profile. These effects correlate with improved confinement and appear largely due to a reduction of anomalous transport as the change in neoclassical transport is small. The observed decrease of turbulent heat flux with increased density gradients is in agreement with nonlinear gyrokinetic simulations, and has been attributed to the unique geometry of W7-X that limits the severity of trapped electron modes.
This work addresses linear transport in turbulent media, with emphasis on neutral particle (atoms, molecules) transport in magnetized fusion plasmas. A stochastic model for turbulent plasmas, based upon a multivariate Gamma distribution, is presented. The geometry is a 2D slab and turbulence is assumed to be statistically homogeneous. The average neutral density and ionization source, which are the quantities relevant for integrated simulations and diagnostic applications, are calculated analytically in the scattering free case. The boundary conditions and the ratio of the turbulence correlation length to the neutral mean free path are identified as the main control parameters in the problem. The non trivial relationship between the average neutral density and the ionization source is investigated. Monte Carlo calculations including scattering are then presented, and the main trends obtained in the scattering free case are shown to be conserved.