Albert Mollén

Comparison of particle trajectories and collision operators for collisional transport in nonaxisymmetric plasmas

In this work, we examine the validity of several common simplifying assumptions used in numerical neoclassical calculations for nonaxisymmetric plasmas, both by using a new continuum drift-kinetic code and by considering analytic properties of the kinetic equation. First, neoclassical phenomena are computed for the LHD and W7-X stellarators using several versions of the drift-kinetic equation, including the commonly used incompressible-E × B-drift approximation and two other variants, corresponding to different effective particle trajectories. It is found that for electric fields below roughly one third of the resonant value, the different formulations give nearly identical results, demonstrating the incompressible E × B-drift approximation is quite accurate in this regime. However, near the electric field resonance, the models yield substantially different results. We also compare results for various collision operators, including the full linearized Fokker-Planck operator. At low collisionality, the radial transport driven by radial gradients is nearly identical for the different operators; while in other cases, it is found to be important that collisions conserve momentum.

Effect of poloidal asymmetries on impurity peaking in tokamaks

Poloidal impurity asymmetries are frequently observed in tokamaks. In this paper, the effect of poloidal asymmetry on electrostatic turbulent transport is studied, including the effect of the E×B drift. Collisions are modeled by a Lorentz operator, and the gyrokinetic equation is solved with a variational approach. The impurity transport is shown to be sensitive to the magnetic shear and changes sign for s≳0.5 in the presence of inboard accumulation. The zero-flux impurity density gradient (peaking factor) is shown to be rather insensitive to collisions in both ion temperature gradient and trapped electron mode driven cases. Our results suggest that the asymmetry (both the location of its maximum and its strength) and the magnetic shear are the two most important parameters that affect the impurity peaking.

A possible mechanism responsible for generating impurity outward flow under radio frequency heating

The effect of poloidal asymmetry of impurities on impurity transport driven by electrostatic turbulence in tokamak plasmas is analyzed. It is found that in the presence of in–out asymmetric impurity populations the zero-flux impurity density gradient (the so-called peaking factor) is significantly reduced. A sign change in the impurity flux may occur if the asymmetry is sufficiently large. This may be a contributing reason for the observed outward convection of impurities in the presence of radio frequency heating. This paper extends a previous work (Fulop and Moradi 2011 Phys. Plasmas 18 030703), by including the effect of ion parallel compressibility on the peaking factor, which is found to have a significant contribution in the presence of poloidal asymmetry. It is shown here that in the ion temperature gradient mode dominated plasmas the presence of an in–out poloidal asymmetry can lead to a negative impurity peaking factor, and it becomes more negative in regions with larger ion temperature gradients. In the trapped electron mode dominated plasmas an in–out poloidal asymmetry results in a strong reduction of the peaking factor; however, it remains positive for typical experimental parameters. Furthermore, it is shown that an up–down asymmetry reduces the peaking factor while an out–in asymmetry increases it.

Microtearing modes in spherical and conventional tokamaks

The onset and characteristics of microtearing modes (MTM) in the core of spherical (NSTX) and conventional tokamaks (ASDEX Upgrade and JET) are studied through local linear gyrokinetic simulations with GYRO (Candy and Belli 2011 General Atomics Report GA-A26818). For experimentally relevant core plasma parameters in the NSTX and ASDEX Upgrade tokamaks, in agreement with previous works, we find MTMs as the dominant linear instability. Also, for JET-like core parameters considered in our study an MTM is found as the most unstable mode. In all of these plasmas, finite collisionality is needed for MTMs to become unstable and the electron temperature gradient is found to be the fundamental drive. However, a significant difference is observed in the dependence of the linear growth rate of MTMs on electron temperature gradient. While it varies weakly and non-monotonically in JET and ASDEX Upgrade plasmas, in NSTX it increases with the electron temperature gradient.

Turbulent transport of impurities and their effect on energy confinement

By presenting linear and nonlinear gyrokinetic studies, based on a balanced neutral beam injection deuterium discharge from the DIII-D tokamak, we demonstrate that impurities alter the scaling of the transport on the charge and mass of the main species, and even more importantly, they can dramatically change the energy transport even in relatively small quantities. A poloidally varying equilibrium electrostatic potential can lead to a strong reduction or sign change of the impurity peaking factor due to the combined effect of the in–out impurity density asymmetry and the E × B drift of impurities. We present an approximate expression for the impurity peaking factor and demonstrate that impurity peaking is not significantly affected by impurity self-collisions.

Impurity transport in trapped electron mode driven turbulence

Trapped electron mode turbulence is studied by gyrokinetic simulations with the GYRO code and an analytical model including the effect of a poloidally varying electrostatic potential. Its impact on radial transport of high-Z trace impurities close to the core is thoroughly investigated, and the dependence of the zero-flux impurity density gradient (peaking factor) on local plasma parameters is presented. Parameters such as ion-to-electron temperature ratio, electron temperature gradient, and main species density gradient mainly affect the impurity peaking through their impact on mode characteristics. The poloidal asymmetry, the safety factor, and magnetic shear have the strongest effect on impurity peaking, and it is shown that under certain scenarios where trapped electron modes are dominant, core accumulation of high-Z impurities can be avoided. We demonstrate that accounting for the momentum conservation property of the impurity-impurity collision operator can be important for an accurate evaluation of the impurity peaking factor.

Impurity transport due to electromagnetic drift wave turbulence

Finite β effects on impurity transport are studied through local linear gyrokinetic simulations with GYRO [J. Candy and E. Belli, General Atomics Report No. GA-A26818, 2011]; in particular, we investigate the parametric dependences of the impurity peaking factor (zero-flux density gradient) and the onset of the kinetic ballooning modes (KBMs). We find that electromagnetic effects even at low β can have significant impact on the impurity transport. The KBM instability threshold depends on the plasma parameters, particularly strongly on plasma shape. We have shown that magnetic geometry significantly influences the results, and the commonly used s-α model overestimates the KBM growth rates and ITG stabilization at high β. In the β range, where the KBM is the dominant instability the impurity peaking factor is strongly reduced, with very little dependence on β and the impurity charge.

Impurities in a non-axisymmetric plasma: Transport and effect on bootstrap current

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.

On collisional impurity transport in nonaxisymmetric plasmas

The presence of impurity species in magnetic confinement fusion devices leads to radiation losses and plasma dilution. Thus it is important to analyze impurity dynamics, and search for means to control them. In stellarator plasmas the neoclassical ambipolar radial electric field often points radially inwards (referred to as the ion root regime), causing impurities to accumulate in the core. This can limit the performance of nonaxisymmetric devices. In the present work we analyze neoclassical impurity transport in stellarator plasmas using a recently developed continuum drift-kinetic solver, the SFINCS code (the Stellarator Fokker- Planck Iterative Neoclassical Conservative Solver). The study is performed for a case close to the edge of W7-X using the standard configuration magnetic geometry. We investigate the sensitivity of impurity transport to impurity charge, main species density and temperature gradients, as well as ion temperature. At the studied radial location we find that the neoclassical impurity peaking factor can be very large, particularly for high-Z impurities. The ambipolar radial electric field is in the ion root regime, and impurity accumulation can thus be expected. The accumulation is strengthened by the large main species density and temperature gradients. Moreover we find that the size of the bootstrap current is affected by the value of the plasma effective charge, suggesting that employing a realistic ion composition can be important when calculating the bootstrap current.

Impurity transport and plasma flow in a mixed collisionality stellarator plasma

The accumulation of impurities in the core of magnetically confined plasmas, resulting from standard collisional transport mechanisms, is a known threat to their performance as fusion energy sources. Whilst the axisymmetric tokamak systems have been shown to benefit from the effect of temperature screening, that is an outward flux of impurities driven by the temperature gradient, impurity accumulation in stellarators was thought to be inevitable, driven robustly by the inward pointing electric field characteristic of hot fusion plasmas. We have shown in Helander et al. (Phys. Rev. Lett, vol. 118, 2017a, 155002) that such screening can in principle also appear in stellarators, in the experimentally relevant mixed collisionality regime, where a highly collisional impurity species is present in a low collisionality bulk plasma. Details of the analytic calculation are presented here, along with the effect of the impurity on the bulk ion flow, which will ultimately affect the bulk contribution to the bootstrap current.