Ion-scale turbulence in MAST: anomalous transport, subcritical transitions, and comparison to BES measurements

Abstract

We investigate the effect of varying the ion temperature gradient (ITG) and toroidal equilibrium scale sheared flow on ion-scale turbulence in the outer core of MAST by means of local gyrokinetic simulations. We show that nonlinear simulations reproduce the experimental ion heat flux and that the experimentally measured values of the ITG and the flow shear lie close to the turbulence threshold. We demonstrate that the system is subcritical in the presence of flow shear, i.e., the system is formally stable to small perturbations, but transitions to a turbulent state given a large enough initial perturbation. We propose that the transition to subcritical turbulence occurs via an intermediate state dominated by low number of coherent long-lived structures, close to threshold, which increase in number as the system is taken away from the threshold into the more strongly turbulent regime, until they fill the domain and a more conventional turbulence emerges. We show that the properties of turbulence are effectively functions of the distance to threshold, as quantified by the ion heat flux. We make quantitative comparisons of correlation lengths, times, and amplitudes between our simulations and experimental measurements using the MAST BES diagnostic. We find reasonable agreement of the correlation properties, most notably of the correlation time, for which significant discrepancies were found in previous numerical studies of MAST turbulence. * ferdinand.vanwyk@physics.ox.ac.uk † highcock@chalmers.se ‡ alex.schekochihin@physics.ox.ac.uk arXiv:1704.02830v2 [physics.plasm-ph] 1 Aug 2017 2 V 0 dV B · ∇φ is the toroidal magnetic flux, V is the volume enclosed by the flux surface, B is the magnetic field, φ is the toroidal angle, and ψ tor,LCFS is the toroidal flux enclosed by the last closed flux surface [see figure 1(b)], ψ pol = (1/2π) 2 V 0 dV B · ∇θ is the poloidal magnetic flux, θ is the poloidal angle, and ψ pol,LCFS is the poloidal flux enclosed by the LCFS. 4 http://w3.pppl.gov/transp/ c , including two cases that match the experimental level of heat flux. This is a considerable improvement over previous nonlinear gyrokinetic simulations of this MAST discharge [51], which overpredicted τ SYNTH c by two orders of magnitude. Examining figure 24(d), we see that (δn i /n i ) SYNTH rms increases with increasing κ T or de-

Overview of new MAST physics in anticipation of first results from MAST Upgrade

“…This can also result in a skewed distribution of the density fluctuations. This is not inconsistent with gyro-kinetic simulations of MAST equilibria, close to threshold, where the turbulence is dominated by a few, isolated, long-lived structures [7]. The simulations suggest that further from the threshold, the symmetry is effectively restored, reducing both the tilt of the correlation function and skewness of the distribution.…”

Simple advecting structures and the edge of chaos in subcritical tokamak plasmas

“…Those two mechanisms would be compatible with a traveling wave in either direction [18,24], but when the nonlinearity in the simulation is turned off, the k y = 0 mode amplitude continues to propagate (not shown) in the same direction for 10a/c s due to the group velocity, which depends on the mean k x value and thus the sign of S [18]. Note that narrow simulations do not permit a vortex pair [13] advection mechanism, where a spatially localised 'blob' self-advects across the domain, as they are dominated by one k y mode (unlike y-localised features seen elsewhere [39]). Time snapshots (fig 3b) of the mid-plane potential for narrow and standard simulations show similar tilting and localization in x for the two simulations, but despite close similarities in y-averaged diagnostics, the standard simulation does not decay towards the narrow edge state.…”

stella: An operator-split, implicit–explicit δf-gyrokinetic code for general magnetic field configurations

“…This local, or flux tube, approximation is routinely used to model microinstabilities and turbulence in axisymmetric magnetic field configurations and has been successfully validated across a range of experiments (cf. [24,25,26,27,28,29]). For non-axisymmetric field configurations, a single flux tube would in principle need to extend infinitely far along the magnetic field line as it ergodically samples a flux surface.…”

“…The dynamics of parallel streaming described by Eqs. (29) and (33) are treated implicitly in stella following a similar approach to that taken by the local, δ f gyrokinetic code GS2 [4]. We discretize Eq.…”