Trapped Electron Mode Turbulence Research Articles

Simultaneous reproduction of experimental profiles, fluxes, transport coefficients, and turbulence characteristics via nonlinear gyrokinetic profile predictions in a DIII-D ITER similar shape plasma

Experimental conditions obtained on the DIII-D tokamak in the ITER Similar Shape (ISS) have been compared extensively with nonlinear gyrokinetic simulation using the CGYRO code [Candy et al., J. Comput. Phys. 324, 73–93 (2016)] with comparisons spanning ion and electron heat fluxes, electron and impurity particle transport, and turbulent fluctuation levels and characteristics. Bayesian optimization techniques [Rodriguez-Fernandez et al., Nucl. Fusion 62(7), 076036 (2022)], combined with nonlinear gyrokinetics, have been used to obtain simultaneously Qi, Qe, and Γe flux-matched profiles that are found to be in good agreement with experimental profile measurements. Synthetic diagnostics were used to compare measured beam emission spectroscopy and correlation electron cyclotron emission turbulent fluctuations with nonlinear simulation. Although some disagreements exist, nonlinear simulations are found to be in generally good agreement with measured fluctuation levels, spectral shapes, and measured radial trends in low-k δne/ne and δTe/Te. Low (Li and C) and mid-Z (Ca) impurity transport was also compared with these flux-matched simulations. Fully stripped, low-Z impurities are well reproduced by the gyrokinetic modeling while clear disagreement exists in comparisons with mid-Z impurities. Nonlinear gyrokinetic investigation into the Z dependence of impurity transport in the ISS conditions is also performed, demonstrating clear trends of impurity diffusion with impurity Z (both D∝Z and D∝1/Z) that vary with the radial location studied. These trends are shown to result from the local dominance of ion temperature gradient or ∇n driven trapped electron mode turbulence and may contribute to the disagreement between simulation and experiment in mid-Z impurity transport. The results of this work represent one of the most complete validation studies of the gyrokinetic model performed to date and provide an example of new capabilities for predicting performance in future fusion devices.

Isotope effects on energy transport in the core of ASDEX-Upgrade tokamak plasmas: Turbulence measurements and model validation

Design and operation of future tokamak fusion reactors using a deuterium–tritium 50:50 mix requires a solid understanding of how energy confinement properties change with ion mass. This study looks at how turbulence and energy transport change in L-mode plasmas in the ASDEX Upgrade tokamak when changing ion species between hydrogen and deuterium. For this purpose, both experimental turbulence measurements and modeling are employed. Local measurements of ion-scale (with wavevector of fluctuations perpendicular to the B-field k⊥< 2 cm−1, k⊥ρs< 0.2, where ρs is the ion sound Larmor radius using the deuterium ion mass) electron temperature fluctuations have been performed in the outer core (normalized toroidal flux ρTor=0.65−0.8) using a multi-channel correlation electron cyclotron emission diagnostic. Lower root mean square perpendicular fluctuation amplitudes and radial correlation lengths have been measured in hydrogen vs deuterium. Measurements of the cross-phase angle between a normal-incidence reflectometer and an ECE signal were made to infer the cross-phase angle between density and temperature fluctuations. The magnitude of the cross-phase angle was found larger (more out-of-phase) in hydrogen than in deuterium. TRANSP power balance simulations show a larger ion heat flux in hydrogen where the electron-ion heat exchange term is found to play an important role. These experimental observations were used as the basis of a validation study of both quasilinear gyrofluid trapped gyro-Landau fluid-SAT2 and nonlinear gyrokinetic GENE codes. Linear solvers indicate that, at long wavelengths (k⊥ρs<1), energy transport in the deuterium discharge is dominated by a mixed ion-temperature-gradient (ITG) and trapped-electron mode turbulence while in hydrogen transport is exclusively and more strongly driven by ITG turbulence. The Ricci validation metric has been used to quantify the agreement between experiments and simulations taking into account both experimental and simulation uncertainties as well as four different observables across different levels of the primacy hierarchy.

Full-flux-surface effects on electrostatic turbulence in Wendelstein 7-X-like plasmas

We present the first nonlinear, gyrokinetic, surface-global simulations of a Wendelstein 7-X-like stellarator with kinetic electrons. As a first application, we investigate the interplay between Ion Temperature Gradient (ITG) and Trapped Electron Mode (TEM) driven turbulence in a Full-Flux-Surface (FFS) approach, as well as the effect of a neoclassical radial electric field, something that escapes the capabilities of flux-tube simulations. We find that even in this more complex setup, ITG turbulence is stabilised through a finite density gradient while TEM turbulence remains relatively weak. Furthermore, we show that the effect of the radial electric field itself is small in comparison with the variation of the gradients. Nevertheless, we observe that for some of the cases shown here, there is not only quantitative but also qualitative disagreement between flux-tube and FFS simulations, in contrast to earlier studies with an adiabatic electron model. These results emphasise the potential importance of retaining geometrical variations on the flux-surface when describing stellarator turbulence under realistic conditions.

Particle pump-out induced by trapped electron mode turbulence in electron cyclotron heated plasmas on XuanLong-50 spherical torus

Particle pump-out effects induced by low-frequency (<200 kHz) density fluctuations were observed in solely electron cyclotron wave (ECW)-heated plasmas on the spherical torus XuanLong-50 (EXL-50) without a central solenoid. The intensity of the relative density fluctuations increases with increasing ECW power and decays when the ECW is turned off while sustaining the plasma current. The electron densities are maintained relatively high and steady when the density fluctuations are completely absent, indicating that the outward transport of electrons is dominated by the particle pump-out effect of the ECW. The density fluctuations are modulated by a supersonic molecular beam injection pulse and the modulation amplitude decreases with increasing electron density at the same ECW injection power and decreasing ECW power at the same electron density, respectively. Analysis revealed that a critical value of electron temperature gradient (ETG) triggers the density fluctuations, and the intensity of the relative density fluctuations is positively correlated with the ETG and approximately inversely proportional to the effective collision frequency. With plasma parameters similar to those of EXL-50 experiments, the HD7 code simulations demonstrate that trapped electron mode (TEM) turbulence can be excited by ETG higher than the critical value observed in the experiment. In addition, the dependence of the mode growth rate (supposed to be proportional to the saturation level of fluctuations in quasi-linear theory) and the measured intensity of the density fluctuations is comparable. The simulated outward particle flux integrated over the poloidal wave number spectrum is significant and proportional to ETG. These observations demonstrate that the density fluctuation is TEM turbulence, which is driven by ETG and induces particle pump-out when the electron density/effective electron collision frequency is low. The potential relevance of this work with the controls of plasma profiles, impurities, helium ash, and heat transport in future reactors of similar low effective collision frequency is also discussed.

Peaking and hollowness of low-Z impurity profiles: an interplay between ITG and TEM induced turbulent transport

The transport of low-Z impurities in fusion plasmas is dominated, in many regimes of interest, by ion temperature gradient (ITG) and trapped electron mode (TEM) turbulence. We use a statistical test-particle approach to analyze the effects of these drift instabilities on the transport coefficients. It is found that the convection of motion by the polarization drift drives an outward radial pinch via ITG and an inward pinch via TEM turbulence. Opposite radial pinches are driven by turbulent motion along magnetic field-lines. The interplay between these effects can explain the co-existence of peaked and hollow impurity density profiles observed in tokamak plasmas. The theoretical framework is validated through comparison with experimental data on the transport of Boron impurities.

ITB formation in gyrokinetic flux-driven ITG/TEM turbulence

The formation mechanism of internal transport barriers (ITBs) in flux-driven turbulence is studied by means of the full-f gyrokinetic code GKNET. In the adiabatic electron case with a weak magnetic shear configuration, toroidal momentum injection can change the radial mean electric field through radial force balance, leading to a kind of driven ITB formation in which the ion thermal diffusivity by ion temperature gradient (ITG) turbulence decreases to the neoclassical transport level. Only cocurrent toroidal rotation in the outer core region can benefit the ITB formation, and this mechanism is identified to originate from a positive feedback loop between the radial shear and resultant momentum flux. On the other hand, in the kinetic electron case with a reversed magnetic shear configuration, robust co-intrinsic rotation is driven near the surface in ITG turbulence and sustains the shear through the radial force balance, leading to the spontaneous reduction of ion turbulent thermal diffusivity, while this is not observed in the adiabatic electron case. In the presence of electron heating, counter-intrinsic rotation by trapped electron mode turbulence is selectively driven in the negative magnetic shear region, which provides steeper shear formation and a resultant larger reduction of ion turbulent thermal diffusivity. This indicates that the co-existence of different modes can trigger the ‘discontinuity’ of mode structure, intrinsic rotation, and resultant mean near , leading to spontaneous ITB formation.

Existence of finite and anisotropic heavy ion parallel compressibility pinch in the gyrokinetic turbulence

Turbulence-driven heavy ion transport in hot magnetized plasma is investigated by means of the gyrokinetic theory and simulations. A finite heavy ion parallel compressibility pinch (${\mathrm{\ensuremath{\Gamma}}}_{s,\ensuremath{\parallel}}$) is found in the gyrokinetic framework, in contrast to the conventional understanding that ${\mathrm{\ensuremath{\Gamma}}}_{s,\ensuremath{\parallel}}$ is negligible. A perturbation theory clarifies the turbulence frequency dependence of ${\mathrm{\ensuremath{\Gamma}}}_{s,\ensuremath{\parallel}}$, resolving the discrepancy with experimental observations. It is also predicted by a nonlocal approach of the parallel advection term that ${\mathrm{\ensuremath{\Gamma}}}_{s,\ensuremath{\parallel}}$ is strongly anisotropic on a magnetic surface. The parameter dependence shows that decreasing the heavy ion mass ${m}_{s}$ strongly enhances ${\mathrm{\ensuremath{\Gamma}}}_{s,\ensuremath{\parallel}}$ through kinetic effects, leading to a deviation from the $1/{m}_{s}$ scaling. Moreover, the pinch direction can be reversed in nonlinear trapped electron mode turbulence through the inverse cascade.

Effects of alpha particles on the transport of helium ash driven by collisionless trapped electron mode turbulence

The effects of alpha (α) particles on the transport of helium ash driven by collisionless trapped electron mode (CTEM) turbulence are analytically studied using quasi-linear theory in tokamak deuterium (D) and tritium (T) plasmas. Under the parameters used in this work, the transport of helium ash is mainly determined by the diffusion due to very weak convection. It is found that the ratio between helium ash diffusivity and effective electron thermal conductivity (D He/χ eff) driven by CTEM turbulence, which is a proper normalized parameter for quantifying the efficiency of helium ash removal, is smaller than unity. This indicates the less efficient removal of helium ash through CTEM turbulence as compared with ion temperature gradient (ITG) turbulence in [Angioni et al 2009 Nucl. Fusion 49 055013]. However, the efficiency of helium ash removal is increased 55% by the presence of 3% α particles with their density gradient being equivalent to that of electrons, and this enhancement can be further strengthened by steeper profile of α particles. This is mainly because the enhancement of helium ash diffusivity by α particles is stronger than that of the effective electron thermal conductivity. Moreover, the higher fraction of T ions, higher temperature ratio between electrons and thermal ions as well as flatter electron density profile, the stronger enhancement of D He/χ eff, and α particles further strengthen the favorable effects of these parameters on the removal of helium ash.

Global gyrokinetic simulations of electrostatic microturbulent transport using kinetic electrons in LHD stellarator

Global gyrokinetic simulations of ion temperature gradient (ITG) and trapped electron mode (TEM) in the LHD stellarator are carried out using the gyrokinetic toroidal code (GTC) with kinetic electrons. ITG simulations show that kinetic electron effects increase the growth rate by more than 50% and more than double the turbulent transport levels compared with simulations using adiabatic electrons. Zonal flow dominates the saturation mechanism in the ITG turbulence. Nonlinear simulations of the TEM turbulence show that the main saturation mechanism is not the zonal flow but the inverse cascade of high to low toroidal harmonics. Further nonlinear simulations with various pressure profiles indicate that the ITG turbulence is more effective in driving heat conductivity whereas the TEM turbulence is more effective for particle diffusivity.

Conceptual design and performance predictions for 2D beam emission spectroscopy turbulence measurements at Wendelstein 7-X

A conceptual design for a 2D beam emission spectroscopy diagnostic system to measure ion gyro-scale plasma turbulence at Wendeslstein 7-X is described. The conceptual design identifies field-aligned viewing geometries and ports for cross-field turbulence measurements in the neutral beam volume. A 2D sightline grid covers the outer plasma region, and the grid configuration provides sufficient k-space coverage in radial and poloidal directions for ion temperature gradient and trapped-electron mode turbulence measurements. Emission intensity estimates, optical transmission losses, and detector noise levels indicate that the measurements will be sensitive to plasma density fluctuations as small as δn/n ≈ 0.5% with a bandwidth of 1MHz. Implementation challenges include a small beam emission Doppler shift due to nearly radial heating beams and reduced optical throughput due to collection aperture limitations.

Limit cycle oscillations, response time, and the time-dependent solution to the Lotka–Volterra predator–prey model

In this work, the time-dependent solution for the Lotka–Volterra predator–prey model is derived with the help of the Lambert W function. This allows an exact analytical expression for the period of the associated limit cycle oscillations and also for the response time between predator and prey population. These results are applied to the predator–prey interaction of zonal density corrugations and turbulent particle flux in gyrokinetic simulations of the collisionless trapped-electron mode turbulence. In the turbulence simulations, the response time is shown to increase when approaching the linear threshold, and the same trend is observed in the Lotka–Volterra model.

Simulation prediction of micro-instability transition and associated particle transport in tokamak plasmas

Two reduced simulation approaches are exploited to predict the parametric boundary of dominant instability regime with global effects and the characteristics of corresponding turbulent particle fluxes in tokamak plasmas. One is usual numerical simulation of coexisting ion temperature gradient (ITG) mode and trapped electron mode (TEM) turbulence employing an extended fluid code (ExFC) based on the so-called Landau–Fluid model including the trapped electron dynamics. Here the density gradient (i.e. R/L n ) driven TEM (∇n-TEM) is emphasized. The other one is a surrogate turbulence transport model, taking a neural network (NN) based approach with speeding calculation. It is shown that the turbulent particle flux, particularly their directions depend on the type of micro-instability as ITG and/or TEM. On the other hand, the density gradient may govern the direction of the turbulent particle fluxes in general circumstances. Specifically, in the parameter regime explored here, the ITG and the electron temperature gradient driven TEM (∇T e-TEM) are destabilized for flat density profile, generally causing an inward particle flux, i.e. particle pinch. Contrarily, for steep density profile, the ∇n-TEM or coexisting ITG and TEM turbulence are dominant so that the particle always diffuses outwards. An empirical criterion is obtained to predict the dominant instability and the direction of particle flux for medium density gradients, involving the gradients of both ion and electron temperature as well as the density. These two transport models are applied to analyze the spontaneous excitation of a quasi-coherent mode in the turbulence modulation discharge by MHD magnetic island observed on tokamak HL-2A, clearly showing a dynamic transition from ITG to TEM. Furthermore, the ExFC-NN model can predict and speed up the analysis of the turbulence transport in tokamak experiments.

Machine learning of turbulent transport in fusion plasmas with neural network

Turbulent transport resulting from drift waves, typically, the ion temperature gradient (ITG) mode and trapped electron mode (TEM), is of great significance in magnetic confinement fusion. It is also well known that turbulence simulation is a challenging issue in both the complex physical model and huge CPU cost as well as long computation time. In this work, a credible turbulence transport prediction model, extended fluid code (ExFC-NN), based on a neural network (NN) approach is established using simulation data by performing an ExFC, in which multi-scale multi-mode fluctuations, such as ITG and TEM turbulence are involved. Results show that the characteristics of turbulent transport can be successfully predicted including the type of dominant turbulence and the radial averaged fluxes under any set of local gradient parameters. Furthermore, a global NN model can well reproduce the radial profiles of turbulence perturbation intensities and fluxes much faster than existing codes. A large number of comparative predictions show that the newly constructed NN model can realize rapid experimental analysis and provide reference data for experimental parameter design in the future.

Importance of gyrokinetic exact Fokker-Planck collisions in fusion plasma turbulence.

Gyrokinetic simulations of turbulence are fundamental to understanding and predicting particle and energy loss in magnetic fusion devices. Previous works have used model collision operators with approximate field-particle terms of unknown accuracy and/or have neglected collisional finite Larmor radius effects. This work moves beyond models to demonstrate important corrections using a gyrokinetic Fokker-Planck collision operator with the exact field-particle terms, in realistic simulations of turbulence in magnetically confined fusion plasmas. The exact operator shows significant corrections for temperature-gradient-driven trapped electron mode turbulence and zonal flow damping, and for microtearing modes in a Joint European Torus pedestal under ITER-like wall conditions. Analysis of the corrections using parameter scans motivates an accurate model which closely reproduces the exact results while reducing computational demands.

Heat pulse propagation and anomalous electron heat transport measurements on the optimized stellarator W7-X

The optimized stellarator Wendelstein 7-X (W7-X) is designed to have an approximately quasi-isodynamic magnetic configuration with reduced neoclassical transport in comparison to a classical stellarator, and turbulent transport is expected to be a significant source of anomalous heat transport across the plasma minor radius. The ion temperature gradient driven mode and the trapped electron mode (TEM) are thought to be responsible for the ion-scale turbulence in W7-X plasmas with volume averaged pressure below 1%. In this work, the electron temperature gradient driven turbulence is shown to be a good candidate for the explanation of the observed electron heat flux, in the inner plasma region where the density gradient is weak (in the outer region, a relatively stronger density gradient would drive additional TEM turbulence). The experimental electron heat transport measured during electron cyclotron resonant heating power and plasma density scans is compared to neoclassical predictions, and the stiffness in the electron heat transport measured during transient transport experiments is presented in three common magnetic configurations of W7-X. In low-⟨β⟩ plasma discharges, the stiffness in the electron heat flux, quantified by the ratio of the heat pulse to power balance diffusivity, , is measured to be less than 2, and trend downwards with increasing collisionality.

The role of collisionless trapped electron mode turbulence on removal of helium ash and transport of deuterium-tritium ions

Removal of helium ash and the anomalous transport of deuterium (D) and tritium (T) ions driven by collisionless trapped electron mode (CTEM) turbulence in tokamak plasmas with weak magnetic shear are studied. We derive the eigenvalue of CTEM with helium ash, and calculate the quasi-linear turbulent fluxes of helium ash, D and T ions simultaneously. Based on the analytical results, the parametric dependence of CTEM instability as well as the anomalous transport of helium ash and D-T ions is investigated, in order to explore the parameter region that is favorable for expelling more helium ash than D and T ions. It is found that helium ash with higher temperature and steeper density profile plays a role of destabilizing CTEM instability, and has higher transport level than that of T ions. We also find that increasing electron temperature and flattening electron density profile are favorable for exhausting helium ash. Isotopic effects (i.e. increasing the fraction of T ions) enhance the transport of both helium ash and D-T ions. Moreover, the trend of stronger transport level of helium ash than that of D-T ions is enhanced by raising electron temperature and flattening electron density profile as well as isotopic effects. Besides, the diffusivity is much larger than the convection. This indicates that the CTEM turbulence driven helium ash transport is favorable for removing helium ash under the parameter region used in the present paper. The possible relevance of our theoretical results to experimental observations is also discussed.

Study of the origin of quasi-coherent modes in low-density KSTAR ECH plasmas

Gradual or abrupt transitions from a coherent mode to quasi-coherent mode (QCM, a type of trapped electron mode turbulence) are observed in the core region of low-density plasmas in Korea Superconducting Tokamak Advanced Research. The coherent mode is initiated in the vicinity of the q = 1 surface by off-axis electron cyclotron resonant heating (ECH), and gradually or abruptly evolves into the QCM as the collisionality is slowly or rapidly increased through electron density control. When the collisionality further increases above a critical level, the QCM is changed into broadband turbulence, which has been preferentially observed in tokamak plasmas. The observed evolution of the QCMs suggests that (1) QCMs can originate from a coherent mode in ECH applied low-collisionality plasmas and (2) quasi-coherent spectra of QCMs are closely related with the collisionality in a finite range (‘collisional broadening of the wavenumber and frequency’). The coherent mode appears to be a metastable electrostatic mode at the rational surface excited by off-axis ECH.

The local nature of the plasma response to cold pulses with electron and ion heating at ASDEX Upgrade

In low density Ohmic heated tokamak plasmas a fast increase of the central electron temperature is regularly observed in response to a drop at the edge produced by the cold pulse of an impurity ablation. This observation has challenged the physical understanding and the theory-based local transport modelling for several decades. Recent experimental results in C-Mod and related modelling (Rodriguez-Fernandez et al 2018 Phys. Rev. Lett. 120 075001) have for the first time provided an explanation for this phenomenology, which results from the stabilization of trapped electron mode turbulence caused by the sudden flattening of the electron density profile as a consequence of the impurity ablation. Starting from this breakthrough result, recent cold pulse experiments in ASDEX Upgrade produced by laser ablation of carbon in deuterium plasmas demonstrate that the fast increase of the central electron temperature is only observed in the presence of dominant electron heating, whereas no fast increase is observed in dominantly ion heated plasmas, even at the same low density. The measured fast flattening of the density profile is shown to be large enough to produce a strong stabilization of the trapped electron mode turbulence in the case where the electron temperature largely exceeds the ion temperature, allowing the fast increase of the central electron temperature, consistent with the C-Mod modelling results. In contrast, it produces a destabilization of the electron temperature gradient modes when the electron and ion temperatures are closer, preventing the increase of the central electron temperature. Finally, modelling which also includes the dynamical evolution of the density profile with a local theory-based transport model shows that the impurity penetration transiently modifies the local turbulence causing a fast penetration of the particles, which explains the observation of a very fast flattening of the electron density profile. By this, the complex density and temperature behaviours in response to a cold pulse are demonstrated to be fully explained by multi-channel interactions within a local model.

Effects of zonal flows on transport crossphase in dissipative trapped-electron mode turbulence in edge plasmas

Confinement regimes with edge transport barriers occur through the suppression of turbulent (convective) fluxes in the particle and/or thermal channels, i.e. , and , respectively, for drift-wave turbulence. The quantity is the turbulence intensity, while and are the crossphases. For H-mode, standard decorrelation theory predicts that it is the turbulence intensity that is mainly affected via flow-induced shearing of turbulent eddies. However, for other regimes (e.g. I-mode, characterized by high energy confinement but low particle confinement), this decrease of turbulence amplitude cannot explain the decoupling of particle versus thermal flux, since a suppression of turbulence intensity would necessarily affect both fluxes the same way. Here, we explore a possible new stabilizing mechanism: zonal flows may directly affect the transport crossphase. We show the effect of this novel mechanism on the turbulent particle flux, by using a simple fluid model (Baver et al 2002 Phys. Plasmas 9 3318) for dissipative trapped-electron mode (DTEM), including zonal flows. We first derive the evolution equation for the transport crossphase δk between density and potential fluctuations, including contributions from the E × B nonlinearity. By using a parametric interaction analysis including the back-reaction on the pump, we obtain a predator–prey like system of equations for the pump amplitude ϕp, the pump crossphase δp, the zonal amplitude ϕz and the triad phase-mismatch Δδ. The system displays limit-cycle oscillations where the instantaneous DTEM growth rate—proportional to the crossphase—shows quasi-periodic relaxations where it departs from that predicted by linear theory.

On the cascading of collisionless trapped-electron mode turbulence in tokamak plasmas

Cascading of collisionless trapped-electron mode turbulence via particle-induced scatterings is analytically investigated for axisymmetric tokamak plasmas. Our results demonstrate that the nonlinear mode-coupling processes depend crucially on both the parallel and radial mode structures associated with the realistic toroidal geometry. As a consequence, the scattering coefficients are found to be significantly suppressed with respect to those derived for slab-like models. For typical tokamak parameters, the wave energy cascades inversely towards lower toroidal mode numbers.