This paper presents a complete theoretical framework for studying turbulence and transport in rapidly rotating tokamak plasmas. The fundamental scale separations present in plasma turbulence are codified as an asymptotic expansion in the ratio ε = ρi/α of the gyroradius to the equilibrium scale length. Proceeding order by order in this expansion, a set of coupled multiscale equations is developed. They describe an instantaneous equilibrium, the fluctuations driven by gradients in the equilibrium quantities, and the transport-timescale evolution of mean profiles of these quantities driven by the interplay between the equilibrium and the fluctuations. The equilibrium distribution functions are local Maxwellians with each flux surface rotating toroidally as a rigid body. The magnetic equilibrium is obtained from the generalized Grad-Shafranov equation for a rotating plasma, determining the magnetic flux function from the mean pressure and velocity profiles of the plasma. The slow (resistive-timescale) evolution of the magnetic field is given by an evolution equation for the safety factor q. Large-scale deviations of the distribution function from a Maxwellian are given by neoclassical theory. The fluctuations are determined by the 'high-flow' gyrokinetic equation, from which we derive the governing principle for gyrokinetic turbulence in tokamaks: the conservation and local (in space) cascade of the free energy of the fluctuations (i.e. there is no turbulence spreading). Transport equations for the evolution of the mean density, temperature and flow velocity profiles are derived. These transport equations show how the neoclassical and fluctuating corrections to the equilibrium Maxwellian act back upon the mean profiles through fluxes and heating. The energy and entropy conservation laws for the mean profiles are derived from the transport equations. Total energy, thermal, kinetic and magnetic, is conserved and there is no net turbulent heating. Entropy is produced by the action of fluxes flattening gradients, Ohmic heating and the equilibration of interspecies temperature differences. This equilibration is found to include both turbulent and collisional contributions. Finally, this framework is condensed, in the low-Mach-number limit, to a more concise set of equations suitable for numerical implementation.

This paper describes a conceptual framework for understanding kinetic plasma turbulence as a generalized form of energy cascade in phase space. It is emphasized that conversion of turbulent energy into thermodynamic heat is only achievable in the presence of some (however small) degree of collisionality. The smallness of the collision rate is compensated for by the emergence of a small-scale structure in the velocity space. For gyrokinetic turbulence, a nonlinear perpendicular phase-mixing mechanism is identified and described as a turbulent cascade of entropy fluctuations simultaneously occurring at spatial scales smaller than the ion gyroscale and in velocity space. Scaling relations for the resulting fluctuation spectra are derived. An estimate for the collisional cutoff is provided. The importance of adequately modelling and resolving collisions in gyrokinetic simulations is briefly discussed, as well as the relevance of these results to understanding the dissipation-range turbulence in the solar wind and the electrostatic microturbulence in fusion plasmas.

Demonstrating improved confinement of energetic ions is one of the key goals of the Wendelstein 7-X (W7-X) stellarator. In the past campaigns, measuring confined fast ions has proven to be challenging. Future deuterium campaigns would open up the option of using fusion-produced neutrons to indirectly observe confined fast ions. There are two neutron populations: 2.45 MeV neutrons from thermonuclear and beam-target fusion, and 14.1 MeV neutrons from DT reactions between tritium fusion products and bulk deuterium. The 14.1 MeV neutron signal can be measured using a scintillating fiber neutron detector, whereas the overall neutron rate is monitored by common radiation safety detectors, for instance fission chambers. The fusion rates are dependent on the slowing-down distribution of the deuterium and tritium ions, which in turn depend on the magnetic configuration via fast ion orbits. In this work, we investigate the effect of magnetic configuration on neutron production rates in W7-X. The neutral beam injection, beam and triton slowing-down distributions, and the fusion reactivity are simulated with the ASCOT suite of codes. The results indicate that the magnetic configuration has only a small effect on the production of 2.45 MeV neutrons from DD fusion and, particularly, on the 14.1 MeV neutron production rates. Despite triton losses of up to 50 %, the amount of 14.1 MeV neutrons produced might be sufficient for a time-resolved detection using a scintillating fiber detector, although only in high-performance discharges.

Electrostatic turbulence in weakly collisional, magnetized plasma can be interpreted as a cascade of entropy in phase space, which is proposed as a universal mechanism for dissipation of energy in magnetized plasma turbulence. When the nonlinear decorrelation time at the scale of the thermal Larmor radius is shorter than the collision time, a broad spectrum of fluctuations at sub-Larmor scales is numerically found in velocity and position space, with theoretically predicted scalings. The results are important because they identify what is probably a universal Kolmogorov-like regime for kinetic turbulence; and because any physical process that produces fluctuations of the gyrophase-independent part of the distribution function may, via the entropy cascade, result in turbulent heating at a rate that increases with the fluctuation amplitude, but is independent of the collision frequency. Introduction. Turbulence is inherently nonlinear and dynamically complicated. In the general case, a broad spectrum of fluctuations is excited, in both wave number and frequency. For turbulent, magnetized plasma, the equations of magnetohydrodynamics provide a pedagogically rich description of the dynamics. However, for those turbulent eddies whose parallel wavelengths (relative to the magnetic field) are comparable to or smaller than the collisional mean free path and whose perpendicular wavelengths are comparable to or smaller than the Larmor radius of one of the constituent species of the plasma, magnetohydrodynamic theory breaks down. In such cases, the gyrokinetic (GK) theory [1, 2] represents a rigorous limit of plasma kinetics for anisotropic (k ≪ k ⊥ ), lowfrequency (ω ≪ Ω, the ion cyclotron frequency) fluctuations. In this Letter, we present a GK description of turbulence in a simplified situation, chosen to isolate a novel phenomenon which is a generic component of all GK turbulence: the simultaneous cascade of entropy to smaller scales in both real space and velocity space. This phase-space cascade is the mechanism by which turbulent energy associated with fluctuating fields is brought to small scales in velocity space, where even very infrequent collisions are sufficient to provide irreversibility and thus heating. Below, we present the theory and first-principles simulations of the phase-space cascade in a homogeneous, electrostatic, magnetized plasma.

Two-dimensional gyrokinetics is a simple paradigm for the study of kinetic magnetised plasma turbulence. In this paper, we present a comprehensive theoretical framework for this turbulence. We study both the inverse and direct cascades (the `dual cascade'), driven by a homogeneous and isotropic random forcing. The key characteristic length of gyrokinetics, the Larmor radius, divides scales into two physically distinct ranges. For scales larger than the Larmor radius, we derive the familiar Charney--Hasegawa--Mima (CHM) equation from the gyrokinetic system, and explain its relationship to gyrokinetics. At scales smaller than the Larmor radius, a dual cascade occurs in phase space (two dimensions in position space plus one dimension in velocity space) via a nonlinear phase-mixing process. We show that at these sub-Larmor scales, the turbulence is self-similar and exhibits power law spectra in position and velocity space. We propose a Hankel-transform formalism to characterise velocity-space spectra. We derive the exact relations for third-order structure functions, analogous to Kolmogorov's four-fifths and Yaglom's four-thirds laws and valid at both long and short wavelengths. We show how the general gyrokinetic invariants are related to the particular invariants that control the dual cascade in the long- and short-wavelength limits. We describe the full range of cascades from the fluid to the fully kinetic range.Comment: Being considered for publication in the Journal of Fluid Mechanics. Version 3 contains more cross referencing including forward referencing, to make the paper easier to read. Final pre-production versio