In the edge of an L-mode tokamak plasma, particle transport and ion energy transport are shown to follow a strong microturbulence (SMT) scaling, whereas in the plasma core the transport is shown to ...follow quasilinear turbulence scaling. The dependence of diffusivity on potential fluctuation amplitude is linear in the SMT regime, and quadratic in the quasilinear regime. The transition to strong microturbulence results from larger E×B drift velocities in the edge compared to the plasma core. At these larger velocities, ions traverse the spatially correlated range faster than the stochastic evolution of the electric potential. Hence, these particles do not experience a time-stochastic field as required by the quasilinear approximation. Instead, scattering of particles in the SMT regime is caused by spatial stochasticity. In contrast, electron energy transport remains quasilinear due to decorrelations caused by collisions and fast parallel motion. Improved understanding of transport beyond quasilinear theory opens the path to more accurate modeling of transport in the tokamak plasma edge.
We describe a new approach to solve the electromagnetic gyrokinetic equations which is optimized for accurate treatment of multispecies Fokker–Planck collisions including both pitch-angle and energy ...diffusion. The new algorithm is spectral/pseudospectral in four of the five phase space dimensions, and in the fieldline direction a novel 5th-order conservative upwind scheme is used to permit high-accuracy electromagnetic simulation even in the limit of very high plasma β and vanishingly small perpendicular wavenumber, k⊥→0. To our knowledge, this is the first pseudospectral implementation of the collision operator in a gyrokinetic code. We show that the new solver agrees closely with GYRO in the limit of weak Lorentz collisions, but gives a significantly more realistic description of collisions at high collision frequency. The numerical methods are also designed to be efficient and scalable for multiscale simulations that treat ion-scale and electron–scale turbulence simultaneously.
The findings of an investigation into the properties of the three dimensional (3D) saturated fluctuation intensity of the electric potential in gyrokinetic turbulence simulations is presented. Scans ...in flux surface elongation and Shafranov shift are used to isolate the tokamak geometric dependencies. The potential intensity required in order to compute exact fluxes by a quasilinear method is determined using linear eigenmodes computed with the gyrokinetic code. A model of this non-linear intensity is constructed using the linear eigenmode properties and the geometry shape functions obtained from the 3D intensity spectrum. The model computes the poloidal wavenumber spectrum of the electron and ion energy fluxes with unprecedented accuracy. New insights are gained into the way zonal flow mixing saturates ion-scale turbulence by controlling the radial wavenumber width of the turbulence spectrum.
In this work, we explore both the potential improvements and pitfalls that arise when using advanced collision models in gyrokinetic simulations of plasma microinstabilities. Comparisons are made ...between the simple-but-standard electron Lorentz operator and specific variations of the advanced Sugama operator. The Sugama operator describes multi-species collisions including energy diffusion, momentum and energy conservation terms, and is valid for arbitrary wavelength. We report scans over collision frequency for both low and high k θ s modes, with relevance for multiscale simulations that couple ion and electron scale physics. The influence of the ion-ion collision terms-not retained in the electron Lorentz model-on the damping of zonal flows is also explored. Collision frequency scans for linear and nonlinear simulations of ion-temperature-gradient instabilities including impurity ions are presented. Finally, implications for modeling turbulence in the highly collisional edge are discussed.
The transport of heat in laboratory and astrophysical plasmas is dominated by the complex nonlinear dynamics of plasma turbulence. In magnetically confined plasmas used for fusion energy research, ...turbulence is responsible for cross-field transport that limits the performance of tokamak reactors. We report a set of novel gyrokinetic simulations that capture ion and electron-scale turbulence simultaneously, revealing the dynamics of cross-scale energy transfer and zonal flow modification that give rise to heat losses. Multi-scale simulations are required to match experimental ion and electron heat fluxes and electron profile stiffness, establishing the applicability of the newly discovered physics to experiment. Importantly, these results provide a likely explanation for the loss of electron heat from tokamak plasmas, the 'great unsolved problem' (Bachelor et al (2007 Plasma Sci. Technol. 9 312-87)) in plasma turbulence and the projected dominant loss channel in ITER.