Search

This paper reports on the development of reduced models for electron temperature gradient (ETG) driven transport in the pedestal. Model development is enabled by a set of 61 nonlinear gyrokinetic simulations with input parameters taken from pedestals in a broad range of experimental scenarios. The simulation data have been consolidated in a new database for gyrokinetic simulation data, the multiscale gyrokinetic database (MGKDB), facilitating the analysis. The modeling approach may be considered a generalization of the standard quasilinear mixing length procedure. The parameter η, the ratio of the density to temperature gradient scale length, emerges as the key parameter for formulating an effective saturation rule. With a single order-unity fitting coefficient, the model achieves an error of 15%. A similar model for ETG particle flux is also described. We also present simple algebraic expressions for the transport informed by an algorithm for symbolic regression.

Nonlinear multiscale gyrokinetic simulations of a Joint European Torus edge pedestal are used to show that electron-temperature-gradient (ETG) turbulence has a rich three-dimensional structure, varying strongly according to the local magnetic-field configuration. In the plane normal to the magnetic field, the steep pedestal electron temperature gradient gives rise to anisotropic turbulence with a radial (normal) wavelength much shorter than in the binormal direction. In the parallel direction, the location and parallel extent of the turbulence are determined by the variation in the magnetic drifts and finite-Larmor-radius (FLR) effects. The magnetic drift and FLR topographies have a perpendicular-wavelength dependence, which permits turbulence intensity maxima near the flux-surface top and bottom at longer binormal scales, but constrains turbulence to the outboard midplane at shorter electron-gyroradius binormal scales. Our simulations show that long-wavelength ETG turbulence does not transport heat efficiently, and significantly decreases overall ETG transport -- in our case by $\sim$40 \% -- through multiscale interactions., 17 pages, 14 figures

We present a comparative transport analysis of the isotope mass dependence in the pedestal of two pairs of deuterium/hydrogen type I ELMy H-mode discharges in JET with ITER-like wall, one characterized by the same input power, the other one by the same stored energy. The investigation, carried out using the gyrokinetic code GENE, focuses on the steep profile region of the pedestal. While large wavenumber modes mainly contribute to the electron heat flux and are scarcely influenced by the main gas isotope, an effect of the ion mass in agreement with the experimental (so called anti-gyro-Bohm) scaling is revealed in the low wavenumber range. In this context, the major role played by the E × B shear in regulating the ion-temperature-gradient turbulence is analyzed in some detail. The competing level of turbulent and neoclassical transport is quantified to shed light on the experimental features of the pedestal profiles at different ion mass, with the particle transport found to be consistent with a higher pedestal top density for increasing isotope masses, and the heat transport shown to match the roughly unaltered observed temperature profiles.

We present the results of GENE gyrokinetic calculations based on a series of JET–ITER-like-wall (ILW) type I ELMy H-mode discharges operating with similar experimental inputs but at different levels of power and gas fuelling. We show that turbulence due to electron-temperature-gradient (ETGs) modes produces a significant amount of heat flux in four JET–ILW discharges, and, when combined with neoclassical simulations, is able to reproduce the experimental heat flux for the two low gas pulses. The simulations plausibly reproduce the high-gas heat fluxes as well, although power balance analysis is complicated by short ELM cycles. By independently varying the normalised temperature gradients ( ω T e ) and normalised density gradients ( ω n e ) around their experimental values, we demonstrate that it is the ratio of these two quantities η e = ω T e / ω n e that determines the location of the peak in the ETG growth rate and heat flux spectra. The heat flux increases rapidly as η e increases above the experimental point, suggesting that ETGs limit the temperature gradient in these pulses. When quantities are normalised using the minor radius, only increases in ω T e produce appreciable increases in the ETG growth rates, as well as the largest increases in turbulent heat flux which follow scalings similar to that of critical balance theory. However, when the heat flux is normalised to the electron gyro-Bohm heat flux using the temperature gradient scale length L T e , it follows a linear trend in correspondence with previous work by different authors.

Mega Ampere Spherical Tokamak (MAST) pedestal data has been analysed, where a pedestal database of 892 shots was obtained, using the upgraded MAST Thomson scattering (TS) diagnostic. Various edge localised mode (ELM) types are discussed, where characteristics and trends of MAST pedestals are shown. The data from the upgraded TS diagnostic confirms pedestal characteristics found in earlier analysis, using previous TS systems. Using the database, calculations of the bootstrap current are obtained using the different analytical formulae (Sauter and Redl), and comparisons are performed. The upgraded MAST TS system now spans the full plasma mid-plane, such that direct comparisons between inboard and outboard pedestals can be obtained, and asymmetries in the density pedestal width were found. To increase confidence in spherical tokamak pedestal predictions, the assumptions of Europed have been validated from the MAST pedestal data, and a value for the kinetic ballooning mode constraint ( C ∼ 0.145 ) has been obtained. The first spherical tokamak pedestal predictions were performed in Europed and compared to experimental values. Using C = 0.145 the temperature pedestal height was predicted to within 10% of the experimental value. In addition type II ELMs on MAST are analysed, and stability analysis and parameter scans have been performed using ESSIVE. Similar magneto-hydrodynamic (MHD) stability properties are observed for type I and type II ELMs, originating from the mixed ELM regimes; it is therefore difficult to distinguish these ELMs using the ideal MHD codes.