Electron Temperature Gradient Turbulence Research Articles

Modeling electron temperature profiles in the pedestal with simple formulas for ETG transport

This paper reports on the refinement (building on (Hatch D.R. et al 2022 Phys. Plasmas 29 062501)) and application of simple formulas for electron heat transport from electron temperature gradient (ETG) driven turbulence in the pedestal. The formulas are improved by (1) improving the parameterization for certain key parameters and (2) carefully accounting for the impact of geometry and shaping in the underlying gyrokinetic simulation database. Comparisons with nonlinear gyrokinetic simulations of ETG transport in the MAST pedestal demonstrate the model’s applicability to spherical tokamaks in addition to standard aspect ratio tokamaks. We identify bounds for model applicability: the model is accurate in the steep gradient region, where the ETG turbulence is largely slab-like, but accuracy decreases as the temperature gradient becomes weaker in the pedestal top. We use the formula to model the electron temperature profile in the pedestal for four experimental scenarios while extensively varying input parameters to represent uncertainties. In all cases, the predicted electron temperature profile exhibits extreme sensitivity to separatrix temperature and density, which has implications for core-edge integration. The model reproduces the electron temperature profile for high ηe=Lne/LTe scenarios but not for low η e scenarios in which microtearing modes (MTMs) have been identified. We develop a proof-of-concept model for MTM transport and explore the relative roles of ETG and MTM in setting the electron temperature profile.

Multi-scale gyrokinetic simulations of JT-60U L-mode plasma: reduction of the ion scale energy loss due to the nonlinear coupling with the electron scale turbulence

We investigate the effect of the electron temperature gradient (ETG) driven turbulence on the energy transport in JT-60U L-mode plasma by means of the multi-scale gyrokinetic simulation. In the core region at r/a=0.5 , the instability in the ion scale is driven by the ion temperature gradient (ITG), meanwhile strong unstable ETG modes are found in the electron scale. The nonlinear multi-scale gyrokinetic simulation shows that ETG modes are stabilized in the nonlinear phase and the energy transport in the multi-scale simulation is similar to that obtained in the ion scale ITG simulation. In an outer region at r/a=0.6 , the ion scale instability changes to be the trapped electron mode (TEM). The multi-scale simulation shows that both the ion and the electron energy flux are reduced by ∼30% compared to those obtained in the single scale TEM simulation. Interestingly, the electron energy flux is close to the experimental value after this reduction. From the data analyses, we find that ETG turbulence damps the energy of TEM modes through the ion scale/electron scale coupling and the electron scale/electron scale coupling, and can be modeled as a turbulent diffusion of TEM modes. These results suggest that the single ion scale simulation seems to be still valid in the inner region with r/a<0.5 . However, in the outer region it is necessary to include the ETG modes in the gyrokinetic simulations to explain the energy transport in this L-mode plasma. This is the first result showing that ETG turbulence can reduce the electron energy loss via the cross-scale interaction in a real tokamak equilibrium profile.

ETG turbulent transport in the Mega Ampere Spherical Tokamak (MAST) pedestal

The results of a gyrokinetic analysis of turbulent transport driven by the electron temperature gradient (ETG) in the MAST pedestal are presented. Local nonlinear simulations from the gyrokinetic code, GENE, show that heat flux produced by ETG turbulence is 10%–30% of the total applied heating power in the upper pedestal and pedestal top during both the pre-edge-localized mode (ELM) (80%–99% inter-ELM period) and post-ELM (0%–20%) periods. Increasing strongly with the ETG, the ETG transport appears to be stiff. Considering radiation losses, ion thermal transport, and the strong sensitivity of the transport to the ETG, we propose that ETG transport is a plausible mechanism mediating the inter-ELM temperature profile on MAST. Cognizant of the possibility that sharp variations in the pedestal profiles may violate the local approximation, we conducted global nonlinear simulations; the results are in good agreement with local simulations except near the pedestal top, where extended radial structures and high transport levels (far beyond experimental) develop in the local simulations. We quantify and parameterize the discrepancy between local and global simulations by calculating the ratio of the radial correlation length to a length scale representative of the profile curvature. When this ratio is sufficiently small, local and global simulations agree as expected.

Cross phases of temperature-gradient-driven turbulence as a model basis for I-mode particle transport

As a model for understanding the type of transport behavior characteristic of the tokamak I mode, cross-phase physics for particle-transport is studied analytically for turbulence dominated by either ion-temperature-gradient (ITG) or electron-temperature-gradient (ETG) instability. I mode is a transport-barrier regime of reduced thermal transport but essentially unaffected particle transport. It is assumed that ITG turbulence applies to the baseline L mode, ETG to I mode, and that E × B flow shear is stronger in I mode, lowering all fluxes. In ITG turbulence, particle transport is governed by trapped electrons. Sensitivity to collisions produces the well-known temperature-gradient-driven pinch that offsets density-gradient-driven outward diffusion, weakening particle transport in L mode. In ETG turbulence, nonadiabatic ions are collisionless. Nonzero transport requires an ion spectrum feature whose magnetic-drift resonance supplies the necessary cross phase. If frequencies of order the ion diamagnetic drift frequency dominate the ion part of the spectrum, as would occur with weakly unstable ITG turbulence, all components of the particle transport are outward and can offset flow-shear-induced flux reductions to produce a flux that is similar to the ITG L-mode particle flux. Nonlinear frequencies are potentially relevant and discussed in relation to I mode.

ETG turbulence in a tokamak pedestal

This paper explores the fundamental characteristics of electron-temperature-gradient (ETG)-driven turbulence in the tokamak pedestal. The extreme gradients in the pedestal produce linear instabilities and nonlinear turbulence that are distinct from the corresponding ETG phenomenology in the core plasma. The linear system exhibits multiple (greater than ten) unstable eigenmodes at each perpendicular wave vector, representing different toroidal and slab branches of the ETG instability. Proper orthogonal decomposition of the nonlinear fluctuations reveals no clear one-to-one correspondence between the linear and nonlinear modes for most wave vectors. Moreover, nonlinear frequencies deviate strongly from those of the linear instabilities, with spectra peaking at positive frequencies, which is opposite the sign of the ETG instability. The picture that emerges is one in which the linear properties are preserved only in a narrow range of k-space. Outside this range, nonlinear processes produce strong deviations from both the linear frequencies and eigenmode structures. This is interpreted in the context of critical balance, which enforces alignment between the parallel scales and fluctuation frequencies. We also investigate the nonlinear saturation processes. We observe a direct energy cascade from the injection scale to smaller scales in both perpendicular directions. However, in the bi-normal direction, there is also nonlocal inverse energy transfer to larger scales. Neither streamers nor zonal flows dominate the saturation.

Stabilization of trapped electron mode through effective diffusion in electron temperature gradient turbulence

Multi-scale gyrokinetic theory and simulations of a toroidal magnetized plasma have revealed the existence of cross-scale interactions of the trapped electron mode (TEM) and the electron temperature gradient (ETG) turbulence. Reduction of the TEM instability growth rate in the ETG turbulence is clearly identified, and is well represented in the form of effective diffusion. A theoretical model based on the stochastic forcing by the ETG turbulence well describes the turbulent diffusion coefficient observed in multi-scale turbulence simulations.

Zonal flow excitation in electron-scale tokamak turbulence

The derivation of an intermediate-scale gyrokinetic-electron theory in nonuniform tokamak plasmas (Chen et al 2021 Nucl. Fusion 61 066017) has shown that a Navier–Stokes type nonlinearity couples electron-temperature-gradient (ETG) modes and zonal flow (ZF) modes with wavelengths much shorter than the ion gyroradius but much longer than the electron gyroradius. This intermediate-scale ETG-ZF coupling is typically stronger than the Hasegawa–Mima type nonlinearity characteristic of the fluid approximation and is predicted to lead to relevant ZF generation and ETG mode regulation. Electron-scale, continuum, gyrokinetic simulation results are presented here which include both single-mode ETG and full-spectrum ETG turbulence. The ZF generation due to single ETG modes is investigated and the single-mode intermediate-scale results are found to be in agreement with theory. The full-spectrum results are then presented and explained qualitatively in terms of the single-mode results. It is found that the ETG-driven ZFs regulate intermediate-scale electron heat flux transport to levels in the predicted range.

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'.

Seeking turbulence optimized configurations for the Wendelstein 7-X stellarator: ion temperature gradient and electron temperature gradient turbulence

We examine the turbulence driven by the ion and electron temperature gradients in selected magnetic configurations of the Wendelstein 7-X (W7-X) stellarator. The inherent flexibility in the configuration space of W7-X enables us to find candidate configurations manifesting low turbulent transport. We follow insights gained by stellarator optimization techniques, in order to identify key geometric features, which are directly related to the ion and electron heat fluxes produced by plasma turbulence. One such a feature is the flux expansion at locations where the curvature is particularly unfavourable. Starting from a configuration routinely used in the W7-X experiment, we end up with an optimized configuration. Based on this equilibrium, we select a configuration from W7-X configuration database with a similar feature as the optimized one. With the help of nonlinear gyrokinetic simulations, we show that the heat flux in this configuration is less stiff than in the initial configuration, both for ion temperature gradient and electron temperature gradient turbulence.

Electromagnetic instabilities and plasma turbulence driven by electron-temperature gradient

Electromagnetic (EM) instabilities and turbulence driven by the electron-temperature gradient (ETG) are considered in a local slab model of a tokamak-like plasma. Derived in a low-beta asymptotic limit of gyrokinetics, the model describes perturbations at scales both larger and smaller than the electron inertial length$d_e$, but below the ion Larmor scale$\rho _i$, capturing both electrostatic and EM regimes of turbulence. The well-known electrostatic instabilities – slab and curvature-mediated ETG – are recovered, and a new instability is found in the EM regime, called the thermo-Alfvénic instability (TAI). It exists in both a slab version (sTAI, destabilising kinetic Alfvén waves) and a curvature-mediated version (cTAI), which is a cousin of the (electron-scale) kinetic ballooning mode. The cTAI turns out to be dominant at the largest scales covered by the model (greater than$d_e$but smaller than$\rho _i$), its physical mechanism hinging on the fast equalisation of the total temperature along perturbed magnetic field lines (in contrast to kinetic ballooning mode, which is pressure balanced). A turbulent cascade theory is then constructed, with two energy-injection scales:$d_e$, where the drivers are slab ETG and sTAI, and a larger (parallel system size dependent) scale, where the driver is cTAI. The latter dominates the turbulent transport if the temperature gradient is greater than a certain critical value, which scales inversely with the electron beta. The resulting heat flux scales more steeply with the temperature gradient than that due to electrostatic ETG turbulence, giving rise to stiffer transport. This can be viewed as a physical argument in favour of near-marginal steady-state in electron-transport-controlled plasmas (e.g. the pedestal). While the model is simplistic, the new physics that is revealed by it should be of interest to those attempting to model the effect of EM turbulence in tokamak-relevant configurations with high beta and large ETGs.

Three-dimensional inhomogeneity of electron-temperature-gradient turbulence in the edge of tokamak plasmas

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 ∼40%—through multiscale interactions.

Integrated modelling and multiscale gyrokinetic validation study of ETG turbulence in a JET hybrid H-mode scenario

Previous studies with first-principle-based integrated modelling suggested that electron temperature gradient (ETG) turbulence may lead to an anti-gyroBohm isotope scaling in JET high-performance hybrid H-mode scenarios. A dedicated comparison study against higher-fidelity turbulence modelling invalidates this claim. Ion-scale turbulence with magnetic field perturbations included, can match the power balance fluxes within temperature gradient error margins. Multiscale gyrokinetic simulations from two distinct codes produce no significant ETG heat flux, demonstrating that simple rules-of-thumb are insufficient criteria for its onset.

Simulation of transport in the FT-2 tokamak up to the electron scale with GENE

Prior experimental work on the FT-2 tokamak has observed electron density fluctuations at electron Larmor radius scales using the enhanced scattering (ES) diagnostic (Gusakov et al 2006 Plasma Phys. Control. Fusion 48 A371–6, Gurchenko and Gusakov 2010 Plasma Phys. Control. Fusion 52 124035). Gyrokinetic GENE simulations of conditions at the upper hybrid resonance layer probed by the ES diagnostic show the presence of the anticipated turbulence from the electron temperature gradient (ETG) driven instability in linear and nonlinear simulations. Ion-scale turbulence is responsible for majority of the transport via trapped electron modes, while impurities act to merge the spectrum of the ion and the electron scale instabilities into a continuum. The linear spectrum at electron scales is characterized by maximal growth rate at a significant ballooning angle θ 0, and at ion scales the turbulence is broad in the ballooning angle distribution. The neoclassical shearing rate obtained from GENE breaks symmetry in nonlinear simulations of ETG turbulence, which manifests itself as an asymmetric turbulence spectrum. The electron density fluctuation spectrum obtained with GENE corresponds well to the ES measurement at electron scales, as do the fluxes obtained from the ion-scale simulations.

Gyrokinetic benchmark of the electron temperature-gradient instability in the pedestal region

Transport from turbulence driven by the electron temperature-gradient (ETG) instability is likely a major source of electron heat losses through the pedestal. Due to extreme gradients and strong shaping, ETG instabilities in the pedestal are distinct from those in the core, having, for example, multiple branches (toroidal and slab) in different wavenumber ranges. Due to its importance for pedestal transport, and its rather exotic character, a rigorous multi-code benchmarking exercise is imperative. Here, we describe such an exercise, wherein we have carried out a detailed comparison of local linear pedestal ETG simulations using three gyrokinetic codes, CGYRO, GEM, and GENE and testing different geometric parameters (such as circular, Miller, and equilibrium EFIT geometry). The resulting linear frequencies, growth rates, and eigenfunctions show very good agreement between the codes in the three types of employed geometries. A nonlinear benchmark between CGYRO and GENE is also described, exhibiting good agreement (a maximum of 20% difference in the heat fluxes computed) at two locations in the pedestal. This lays the foundation for confidently modeling ETG turbulence in the pedestal.

Excitation of zonal flow by intermediate-scale toroidal electron temperature gradient turbulence

On the basis of gyrokinetic theory, we derive nonlinear equations for the zonal flow (ZF) generation in intermediate-scale electron temperature gradient (ETG) turbulence (with wavelength much shorter than the ion Larmor radius but much longer than the electron Larmor radius) in nonuniform tokamak plasmas. Both the spontaneous and forced generation of ZFs are kept on the same footing. The resultant Schrödinger equation for the ETG amplitude is characterized by a Navier–Stokes type nonlinearity, which is typically stronger than the Hasegawa–Mima type nonlinearity resulting from the fluid approximation. The physics underlying the three stages of ZF generation process is clarified, and the role of parallel mode structure decoupling is discussed. It is found that ZFs can be more easily excited in the intermediate-scale ETG turbulence than in the short wavelength regime.

Testing predictions of electron scale turbulent pedestal transport in two DIII-D ELMy H-modes

In this paper, we present linear and nonlinear gyrokinetic analyses in the pedestal region of two DIII-D ELMy H-mode discharges using the CGYRO code. The otherwise matched discharges employ different divertor configurations to investigate the impact of varying recycling and particle source on pedestal profiles. Linear gyrokinetic simulations find electrostatic ion-scale instabilities (ion temperature gradient and trapped electron modes, ITG–TEM) are present just inside the top of the pedestal with growth rates that are enhanced significantly by parallel velocity shear. In the sharp gradient region, E × B shearing rates are comparable or larger than ion scale growth rates, suggesting the suppression of ITG–TEM modes in this region. Instead, the electron temperature profiles are found to be correlated with and just above the electron temperature gradient (ETG) instability thresholds. Using gradients varied within experimental uncertainties, nonlinear electron-scale gyrokinetic simulations predict electron heat fluxes from ETG turbulence, that when added to neoclassical (NC) ion thermal transport simulated by NEO, account for 30%–60% of the total experimental heat flux. In addition, the NC electron particle flux is found to contribute significantly to the experimental fluxes inferred from SOLPS-ITER analysis. Additional nonlinear gyrokinetic simulations are run varying input gradients to develop a threshold-based reduced model for ETG transport, finding a relatively simple dependence on η e = L ne/L Te. Predictive transport simulations are used to validate this pedestal-specific ETG model, in conjunction with a model for NC particle transport. In both discharges, the predicted electron temperatures are always overpredicted, indicative of the insufficient stiffness in the ETG pedestal model to account for all of the experimental electron thermal transport. In the case of the closed divertor discharge with lower particle source, the predicted electron density is close to the experiment, consistent with the magnitude of NC particle transport in that discharge. However, the density profiles are overpredicted in the open divertor discharge (larger particle source), due to insufficient model transport. The implications for other mechanisms accounting for the remainder of transport in the sharp gradient region in the two discharges are discussed.

Quantitative comparisons of electron-scale turbulence measurements in NSTX via synthetic diagnostics for high-k scattering

Two synthetic diagnostics are implemented for the high-k scattering system in NSTX (Smith et al 2008 Rev. Sci. Instrum. 79 123501) allowing direct comparisons between the synthetic and experimentally detected frequency and wavenumber spectra of electron-scale turbulence fluctuations. Synthetic diagnostics are formulated in real-space and in wavenumber space, and are deployed in realistic electron-scale simulations carried out with the GYRO code (Candy and Waltz 2003 J. Comput. Phys. 186 545). A highly unstable electron temperature gradient (ETG) mode regime in a modest-β NSTX NBI-heated H-mode discharge is chosen for the analysis. Mapping the measured wavenumbers to field aligned coordinates shows that the high-k system is sensitive to fluctuations that are closer to the spectral peak in the density fluctuation wavenumber spectrum (streamers) than originally predicted. The analyses of synthetic spectra show that the frequency response of the detected fluctuations is dominated by Doppler shift and is insensitive to the turbulence drive. The shape of the high-k density fluctuation wavenumber spectrum is sensitive to the ETG turbulence drive conditions, and can be reproduced in a sensitivity scan of the most pertinent turbulent drive terms in the simulation.

Electron temperature gradient turbulence induced energy flux in the large volume plasma device

The Large Volume Plasma Device (LVPD) has successfully demonstrated excitation of the Electron Temperature Gradient (ETG) driven turbulence in the finite plasma beta (β∼0.06−0.4) condition, where the threshold condition for ETG turbulence is, ηETG=Ln/LT&amp;gt;2/3 satisfied, where, Ln=1ndndx−1 is the density scale length and LTe=1TedTedx−1 is the temperature scale lengths [Mattoo et al., Phys. Rev. Lett. 108, 255007 (2012)]. The observed mode follows wave vector scaling and frequency ordering as k⊥ρe≤1 ≪ k⊥ρi, Ωi&amp;lt;ω ≪ Ωe, where k⊥ is the perpendicular wave vector, ρe, ρi are Larmor radii of the electron and ion, respectively, and Ωi, Ωe, ω are the ion, electron gyro frequencies and the mode frequency, respectively. Simultaneous measurement of fluctuations in electron temperature, δTe ∼ (10−30) %, plasma density, δne ∼ (5−12) %, and potential δVf ∼ (1−10) % are obtained. A strong negative correlation with correlation coefficients Cδn−δφ ∼−0.8 and CδT−δφ ∼−0.9 is observed between the density and potential and temperature and potential fluctuations, respectively. These correlated density, temperature, and potential fluctuations lead to the generation of turbulent heat flux. The measured heat flux is compared with the theoretically estimated heat flux from ETG model equations. The experimental result shows that the net heat flux is directed radially outward.

Initial transport and turbulence analysis and gyrokinetic simulation validation in NSTX-U L-mode plasmas

Transport analysis, ion-scale turbulence measurements, and initial linear and nonlinear gyrokinetic simulations are reported for a transport validation study based on low aspect ratio NSTX-U L-mode discharges. The relatively long, stationary L-modes enabled by the upgraded centerstack provide a more ideal target for transport validation studies that were not available during NSTX operation. Transport analysis shows that anomalous electron transport dominates energy loss while ion thermal transport is well described by neoclassical theory. Linear gyrokinetic GYRO analysis predicts that ion temperature gradient (ITG) modes are unstable around normalized radii ρ = 0.6–0.8, although E × B shearing rates are larger than the linear growth rates over much of that region. Deeper in the core (ρ = 0.4–0.6), electromagnetic microtearing modes (MTM) are unstable as a consequence of the relatively high beta and collisionality in these particular discharges. Consistent with the linear analysis, local, nonlinear ion-scale GYRO simulations predict strong ITG transport at ρ = 0.76, whereas electromagnetic MTM transport is important at ρ = 0.47. The prediction of ion-scale turbulence is consistent with 2D beam emission spectroscopy (BES) that measures the presence of broadband ion-scale fluctuations. Interestingly, the BES measurements also indicate the presence of bi-modal poloidal phase velocity propagation that could be indicative of two different turbulence types. However, in the region between (ρ = 0.56, 0.66), ion-scale simulations are strongly suppressed by the locally large E × B shear. Instead, electron temperature gradient (ETG) turbulence simulations predict substantial transport, illustrating electron-scale contributions can be important in low aspect ratio L-modes, similar to recent analysis at conventional aspect ratio. However, agreement within experimental uncertainties has not been demonstrated, which requires additional simulations to test parametric sensitivities. The potential need to include profile-variation effects (due to the relatively large value of ρ* = ρi/a at low aspect ratio), including electromagnetic and possibly multi-scale effects, is also discussed.

Observation of electromagnetic fluctuation induced particle transport in ETG dominated large laboratory plasma

The large volume plasma device (LVPD), a cylindrically shaped, linear plasma device of dimension (Φ = 2 m, L = 3 m) have successfully demonstrated the excitation of electron temperature gradient (ETG) turbulence (Mattoo et al 2012 Phys. Rev. Lett. 108 1–5). The observed ETG turbulence shows significant power between corresponding to the wave number, where, and ion cyclotron frequency and electron Larmor radius, respectively. The observed frequency and wave number matches well with theoretical estimates corresponding to Whistler-ETG mode. We investigated electromagnetic (EM) fluctuations induced plasma transport in high beta ETG mode suitable plasma in LVPD. The radial EM electron (ion) flux results primarily from the correlation between fluctuations of parallel electron current and radial magnetic field The EM particle flux is observed to be much smaller than the electrostatic particle flux, The EM flux is small but finite contrary to the conventional slab ETG model. The EM flux is non-ambipolar. A theoretical model is obtained for the EM particle flux in straight homogeneous magnetic field geometry. The theoretical estimates are seen to compare well with the experimental observations. Sluggish parallel ion response is identified as the key mechanism for generation of small but finite EM flux.