91̽»¨

Abstract:

Multiple space and time scales arise in plasma turbulence in magnetic confinement fusion devices because of the smallness of the square root of the electron-to-ion mass ratio ${\left({m}_{{\rm{e}}}/{m}_{{\rm{i}}}\right)}^{1/2}$ and the consequent disparity of the ion and electron thermal gyroradii and thermal speeds. Direct simulations of this turbulence that include both ion and electron space–time scales indicate that there can be significant interactions between the two scales. The extreme computational expense and complexity of these direct simulations motivates the desire for reduced treatment. By exploiting the scale-separation between ion scales (IS) and electron scales (ES), and expanding the gyrokinetic equations for the turbulence in ${\left({m}_{{\rm{e}}}/{m}_{{\rm{i}}}\right)}^{1/2}$, we derive such a reduced system of gyrokinetic equations that describes cross-scale interactions. The coupled gyrokinetic equations contain novel terms which provide candidate mechanisms for the observed cross-scale interaction. The ES turbulence experiences a modified drive due to gradients in the IS distribution function, and is advected by the IS $E\times B$ drift, which varies in the direction parallel to the magnetic field line. The largest possible cross-scale term in the IS equations is sub-dominant in our ${\left({m}_{{\rm{e}}}/{m}_{{\rm{i}}}\right)}^{1/2}$ expansion. Hence, in our model the IS turbulence evolves independently of the ES turbulence. To complete the scale-separated approach, we provide and justify a parallel boundary condition for the coupled gyrokinetic equations in axisymmetric equilibria based on the standard 'twist-and-shift' boundary condition. This approach allows one to simulate multi-scale turbulence using ES flux tubes nested within an IS flux tube.

Abstract:

Here we present details of an operator-split, implicit-explicit numerical scheme for the solution of the gyrokinetic-Poisson system of equations in the local limit. This scheme has been implemented in a new code called stella, which is capable of evolving electrostatic fluctuations with full kinetic electron effects and an arbitrary number of ion species in general magnetic geometry. We demonstrate the advantages of this mixed approach over a fully explicit treatment and provide linear and nonlinear benchmark comparisons for both axisymmetric and non-axisymmetric magnetic equilibria.

Abstract:

© 2018 IOP Publishing Ltd. Differential rotation is induced in tokamak plasmas when an underlying symmetry of the governing gyrokinetic-Maxwell system of equations is broken. One such symmetry-breaking mechanism is considered here: the turbulent acceleration of particles along the mean magnetic field. This effect, often referred to as the 'parallel nonlinearity', has been implemented in the δf gyrokinetic code stella and used to study the dependence of turbulent momentum transport on the plasma size and on the strength of the turbulence drive. For JET-like parameters with a wide range of driving temperature gradients, the momentum transport induced by the inclusion of turbulent acceleration is similar to or smaller than the ratio of the ion Larmor radius to the plasma minor radius. This low level of momentum transport is explained by demonstrating an additional symmetry that prohibits momentum transport when the turbulence is driven far above marginal stability.

Abstract:

The control of impurity accumulation is one of the main challenges for future stellarator fusion reactors. The standard argument to explain this accumulation relies on the, in principle, large inward pinch in the neoclassical impurity flux caused by the typically negative radial electric field in stellarators. This simplified interpretation was proven to be flawed by Helander et al. [Phys. Rev. Lett. 118, 155002 (2017)], who showed that in a relevant regime (low-collisionality main ions and collisional impurities) the radial electric field does not drive impurity transport. In that reference, the effect of the component of the electric field that is tangent to the magnetic surface was not included. In this Letter, an analytical calculation of the neoclassical radial impurity flux incorporating such effect is given, showing that it can be very strong for highly charged impurities and that, once it is taken into account, the dependence of the impurity flux on the radial electric field reappears. Realistic examples are provided in which the inclusion of the tangential electric field leads to impurity expulsion.

Abstract:

Using a kinetic model for the ions and adiabatic electrons, we solve a steady state, electron-repelling magnetic presheath in which a uniform magnetic field makes a small angle $\alpha \ll 1$ (in radians) with the wall. The presheath characteristic thickness is the typical ion gyroradius $\rho_{\text{i}}$. The Debye length $\lambda_{\text{D}}$ and the collisional mean free path of an ion $\lambda_{\text{mfp}}$ satisfy the ordering $\lambda_{\text{D}} \ll \rho_{\text{i}} \ll \alpha \lambda_{\text{mfp}}$, so a quasineutral and collisionless model is used. We assume that the electrostatic potential is a function only of distance from the wall, and it varies over the scale $\rho_{\text{i}}$. Using the expansion in $\alpha \ll 1$, we derive an analytical expression for the ion density that only depends on the ion distribution function at the entrance of the magnetic presheath and the electrostatic potential profile. Importantly, we have added the crucial contribution of the orbits in the region near the wall. By imposing the quasineutrality equation, we derive a condition that the ion distribution function must satisfy at the magnetic presheath entrance --- the kinetic equivalent of the Chodura condition. Using an ion distribution function at the entrance of the magnetic presheath that satisfies the kinetic Chodura condition, we find numerical solutions for the self-consistent electrostatic potential, ion density and flow across the magnetic presheath for several values of $\alpha$. Our numerical results also include the distribution of ion velocities at the Debye sheath entrance. We find that at small values of $\alpha$ there are substantially fewer ions travelling with a large normal component of the velocity into the wall.