No Arabic abstract
Impurities cause radiation losses and plasma dilution, and in stellarator plasmas the neoclassical ambipolar radial electric field is often unfavorable for avoiding strong impurity peaking. In this work we use a new continuum drift-kinetic solver, the SFINCS code (the Stellarator Fokker-Planck Iterative Neoclassical Conservative Solver) [M. Landreman et al., Phys. Plasmas 21 (2014) 042503] which employs the full linearized Fokker-Planck-Landau operator, to calculate neoclassical impurity transport coefficients for a Wendelstein 7-X (W7-X) magnetic configuration. We compare SFINCS calculations with theoretical asymptotes in the high collisionality limit. We observe and explain a 1/nu-scaling of the inter-species radial transport coefficient at low collisionality, arising due to the field term in the inter-species collision operator, and which is not found with simplified collision models even when momentum correction is applied. However, this type of scaling disappears if a radial electric field is present. We also use SFINCS to analyze how the impurity content affects the neoclassical impurity dynamics and the bootstrap current. We show that a change in plasma effective charge Zeff of order unity can affect the bootstrap current enough to cause a deviation in the divertor strike point locations.
The bootstrap current and flow velocity of a low-collisionality stellarator plasma are calculated. As far as possible, the analysis is carried out in a uniform way across all low-collisionality regimes in general stellarator geometry, assuming only that the confinement is good enough that the plasma is approximately in local thermodynamic equilibrium. It is found that conventional expressions for the ion flow speed and bootstrap current in the low-collisionality limit are accurate only in the $1/ u$-collisionality regime and need to be modified in the $sqrt{ u}$-regime. The correction due to finite collisionality is also discussed and is found to scale as $ u^{2/5}$.
An asymptotic expansion is performed to obtain quasi-axisymmetric magnetic configurations that are weakly non-axisymmetric. A large space of solutions is identified, which satisfy the condition of quasi-axisymmetry on a single magnetic flux surface, while (non-axisymmetric) globally quasi-axisymmetric solutions are shown to not exist, agreeing with the conclusions of previous theoretical work. The solutions found are shown to be geometrically constrained at low aspect ratio or high toroidal period number. Solutions satisfying the more general condition of omnigeneity (generalized quasi-axisymmetry) are also shown to exist, and it is found that quasi-axisymmetric deformations can be superposed with an omnigenous solution, while preserving the property of omnigeneity, effectively extending the space of good configurations. A numerical solution of the first order quasi-axisymmetry problem is demonstrated and compared with solutions found with a widely used MHD equilibrium solver, independently verifying that quasi-axisymmetry is satisfied at the appropriate order. It is thereby demonstrated that approximately quasi-axisymmetric solutions can be directly constructed, i.e. without using numerical search algorithms.
We present a numerical study of collisional transport in a tokamak pedestal in the presence of non-trace impurities, using the radially global $delta f$ neoclassical solver PERFECT [M. Landreman et al. 2014 Plasma Phys. Control. Fusion 56 045005]. It is known that in a tokamak core with non-trace impurities present the radial impurity flux opposes the bulk ion flux to provide an ambipolar particle transport, with the electron transport being negligibly small. However, in a sharp density pedestal with sub-sonic ion flows the electron transport can be comparable to the ion and impurity flows. Furthermore, the neoclassical particle transport is not intrinsically ambipolar, and the non-ambipolarity of the fluxes extends outside the pedestal region by the radial coupling of the perturbations. The neoclassical momentum transport, which is finite in the presence of ion orbit-width scale profile variations, is significantly enhanced when impurities are present in non-trace quantities, even if the total parallel mass flow is dominated by the bulk ions.
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. In the present work the diffusive contribution is also calculated in other stellarator plasmas. In particular, the diffusion (D) and convection (V) coefficients of carbon and iron impurities produced by ion-temperature-gradient (ITG) turbulence are obtained for W7-X, LHD, TJ-II and NCSX. 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 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 V and a stronger decrease of D 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.
Carbon and nitrogen impurity transport coefficients are determined from gas puff experiments carried out during repeat L-mode discharges on the Mega-Amp Spherical Tokamak (MAST) and compared against a previous analysis of helium impurity transport on MAST. The impurity density profiles are measured on the low-field side of the plasma, therefore this paper focuses on light impurities where the impact of poloidal asymmetries on impurity transport is predicted to be negligible. A weak screening of carbon and nitrogen is found in the plasma core, whereas the helium density profile is peaked over the entire plasma radius.