Critically balanced ion temperature gradient turbulence in fusion plasmas

Comparing pedestal structure in JET-ILW H-mode plasmas with a model for stiff ETG turbulent heat transport

A predictive model for the electron temperature profile of the H-mode pedestal is described, and its results are compared with the pedestal structure of JET-ILW plasmas. The model is based on a scaling for the gyro-Bohm normalized, turbulent electron heat flux [Formula: see text] resulting from electron temperature gradient (ETG) turbulence, derived from results of nonlinear gyrokinetic (GK) calculations for the steep gradient region. By using the local temperature gradient scale length [Formula: see text] in the normalization, the dependence of [Formula: see text] on the normalized gradients [Formula: see text] and [Formula: see text] can be represented by a unified scaling with the parameter [Formula: see text], to which the linear stability of ETG turbulence is sensitive when the density gradient is sufficiently steep. For a prescribed density profile, the value of [Formula: see text] determined from this scaling, required to maintain a constant electron heat flux [Formula: see text] across the pedestal, is used to calculate the temperature profile. Reasonable agreement with measurements is found for different cases, the model providing an explanation of the relative widths and shifts of the [Formula: see text] and [Formula: see text] profiles, as well as highlighting the importance of the separatrix boundary conditions. Other cases showing disagreement indicate conditions where other branches of turbulence might dominate. This article is part of a discussion meeting issue 'H-mode transition and pedestal studies in fusion plasmas'.

Gyrokinetic electrostatic turbulence close to marginality in the Wendelstein 7-X stellarator

The transition from strong (fluidlike) to nearly marginal (Floquet-type) regimes of ion-temperature-gradient (ITG) driven turbulence is studied in the stellarator Wendelstein 7-X by means of numerical simulations. Close to marginality, extended (along magnetic field lines) linearly unstable modes are dominant, even in the presence of kinetic electrons, and provide a drive that results in finite turbulent transport. A total suppression of turbulence above the linear stability threshold of the ITG modes, commonly present in tokamaks and known as the "Dimits shift," is not observed. We show that this is mostly due to the peculiar radial structure of marginal turbulence, which is more localized than in the fluid case and therefore less likely to be stabilized by shearing flows.

Stellarators Resist Turbulent Transport on the Electron Larmor Scale

Electron temperature gradient (ETG)-driven turbulence, despite its ultrafine scale, is thought to drive significant thermal losses in magnetic fusion devices-but what role does it play in stellarators? The first numerical simulations of ETG turbulence for the Wendelstein 7-X stellarator, together with power balance analysis from its initial experimental operation phase, suggest that the associated transport should be negligible compared to other channels. The effect, we argue, originates essentially from the geometric constraint of multiple field periods, a generic feature of stellarators.

Distinct Turbulence Saturation Regimes in Stellarators

In the complex 3D magnetic fields of stellarators, ion-temperature-gradient turbulence is shown to have two distinct saturation regimes, as revealed by petascale numerical simulations and explained by a simple turbulence theory. The first regime is marked by strong zonal flows and matches previous observations in tokamaks. The newly observed second regime, in contrast, exhibits small-scale quasi-two-dimensional turbulence, negligible zonal flows, and, surprisingly, a weaker heat flux scaling. Our findings suggest that key details of the magnetic geometry control turbulence in stellarators.

Linear technique to understand non-normal turbulence applied to a magnetized plasma

In nonlinear dynamical systems with highly nonorthogonal linear eigenvectors, linear nonmodal analysis is more useful than normal mode analysis in predicting turbulent properties. However, the nontrivial time evolution of nonmodal structures makes quantitative understanding and prediction difficult. We present a technique to overcome this difficulty by modeling the effect that the advective nonlinearities have on spatial turbulent structures. The nonlinearities are taken as a periodic randomizing force with time scale consistent with critical balance arguments. We apply this technique to a model of drift wave turbulence in the Large Plasma Device [W. Gekelman et al., Rev. Sci. Instrum. 62, 2875 (1991)], where nonmodal effects dominate the turbulence. We compare the resulting growth rate spectra to the spectra obtained from a nonlinear simulation, showing good qualitative agreement, especially in comparison to the eigenmode growth rate spectra.

Transition between saturation regimes of gyrokinetic turbulence

A gyrokinetic model of ion temperature gradient driven turbulence in magnetized plasmas is used to study the injection, nonlinear redistribution, and collisional dissipation of free energy in the saturated turbulent state over a broad range of driving gradients and collision frequencies. The dimensionless parameter L(T)/L(C), where L(T) is the ion temperature gradient scale length and L(C) is the collisional mean free path, is shown to parametrize a transition between a saturation regime dominated by nonlinear transfer of free energy to small perpendicular (to the magnetic field) scales and a regime dominated by dissipation at large scales in all phase space dimensions.

Experimental signatures of critically balanced turbulence in MAST

Beam emission spectroscopy (BES) measurements of ion-scale density fluctuations in the MAST tokamak are used to show that the turbulence correlation time, the drift time associated with ion temperature or density gradients, the particle (ion) streaming time along the magnetic field, and the magnetic drift time are consistently comparable, suggesting a "critically balanced" turbulence determined by the local equilibrium. The resulting scalings of the poloidal and radial correlation lengths are derived and tested. The nonlinear time inferred from the density fluctuations is longer than the other times; its ratio to the correlation time scales as ν(*i)(-0.8 ± 0.1), where ν(*i) = ion  collision rate/streaming rate. This is consistent with turbulent decorrelation being controlled by a zonal component, invisible to the BES, with an amplitude exceeding those of the drift waves by ∼ ν(*i)(-0.8).

Zero-turbulence manifold in a toroidal plasma

Sheared toroidal flows can cause bifurcations to zero-turbulent-transport states in tokamak plasmas. The maximum temperature gradients that can be reached are limited by subcritical turbulence driven by the parallel velocity gradient. Here it is shown that q/ϵ (magnetic field pitch/inverse aspect ratio) is a critical control parameter for sheared tokamak turbulence. By reducing q/ϵ, far higher temperature gradients can be achieved without triggering turbulence, in some instances comparable to those found experimentally in transport barriers. The zero-turbulence manifold is mapped out, in the zero-magnetic-shear limit, over the parameter space (γ(E), q/ϵ, R/L(T)), where γ(E) is the perpendicular flow shear and R/L(T) is the normalized inverse temperature gradient scale. The extent to which it can be constructed from linear theory is discussed.

Locality and universality in gyrokinetic turbulence

The nature of nonlinear interactions in gyrokinetic turbulence, driven by the ion-temperature gradient instability, is investigated using direct numerical simulations in toroidal flux tube geometry. To account for the level of separation existing between scales involved in an energetic interaction, the degree of locality of the free energy scale flux is analyzed employing Kraichnan's infrared (IR) and ultraviolet locality functions. Because of the nontrivial dissipative nature of gyrokinetic turbulence, an asymptotic level for the locality exponents, indicative of a universal dynamical regime for gyrokinetics, is not recovered and an accentuated nonlocal behavior of the IR interactions is found instead, in spite of the local energy cascade observed.

Turbulent transport and heating of trace heavy ions in hot magnetized plasmas

Scaling laws for the transport and heating of trace heavy ions in low-frequency magnetized plasma turbulence are derived and compared with direct numerical simulations. The predicted dependences of turbulent fluxes and heating on ion charge and mass number are found to agree with numerical results for both stationary and differentially rotating plasmas. Heavy ion momentum transport is found to increase with mass, and heavy ions are found to be preferentially heated, implying a mass-dependent ion temperature for very weakly collisional plasmas and for partially ionized heavy ions in strongly rotating plasmas.

Scaling of spontaneous rotation with temperature and plasma current in tokamaks

Using theoretical arguments, a simple scaling law for the size of the intrinsic rotation observed in tokamaks in the absence of a momentum injection is found: The velocity generated in the core of a tokamak must be proportional to the ion temperature difference in the core divided by the plasma current, independent of the size of the device. The constant of proportionality is of the order of 10 km·s(-1)·MA·keV(-1). When the intrinsic rotation profile is hollow, i.e., it is countercurrent in the core of the tokamak and cocurrent in the edge, the scaling law presented in this Letter fits the data remarkably well for several tokamaks of vastly different size and heated by different mechanisms.