Temperature Gradient Scale Research Articles

Impurity transport driven by kinetic ballooning mode in the strong gradient pedestal of tokamak plasmas

Abstract The impurity transport driven by kinetic ballooning mode (KBM) is theoretically studied in the DIII-D H-mode strong gradient pedestal plasmas. From the electromagnetic gyrokinetic equation, including the correction of the strong radial electric field, the dispersion relationship of KBM instability with non-trace impurity is firstly derived. Then, the turbulent impurity flux and ion heat flux, as well as the associated transport coefficients, are further calculated. Through the parametric dependence analysis of analytical results, it is found that dilution effects of light fully ionized impurities can reduce the drive of KBM by affecting the kinetic pressure gradient parameter α and diamagnetic effects, thus leading to a decrease in both the absolute value of the real frequency | ω r ′ | and the growth rate γ k ′ of KBM instability. Stronger dilution effects by increasing the impurity charge number Z or steepening the impurity density profile correspond to stronger effects. Moreover, the removal efficiency of light fully ionized impurities, quantified by the ratio between the impurity diffusivity and effective ion heat conductivity D z χ i eff ≈ 1 − b z + 4 ω Dz / ω r ′ 3 ( 1 + 1 / η i ) / 2 , increases with an increase of Z mainly due to the smaller impurity finite Larmor radius (FLR) effects reflected by b z ∝ 1 / Z . Besides, the increase of the impurity density gradient can significantly enhance D z / χ i eff , and this is because stronger impurity dilution effects make a larger magnetic drift term | ω Dz | / | ω r ′ | ( ω Dz is the magnitude of impurity magnetic drift frequency) and η i (the ratio of ion density gradient scale length to ion temperature gradient scale length). For heavy metal impurities with a concentration of 10 − 4 , the peaking factor (PF) is positive, which means that its density profile is inwardly peaked, and the PF decreases with the enhancement of impurity FLR effects. These results may provide some theoretical reference on understanding the physical mechanism of impurity transport in the pedestal of H-mode plasmas.

Nonlocal suppression of Biermann battery magnetic-field generation for arbitrary atomic numbers and magnetization

The Biermann battery term of magnetohydrodynamics (MHD) generates a magnetic field where electron density gradients and electron temperature gradients are perpendicular to one another. Kinetic simulations and experiments have shown that the rate of magnetic-field generation is lower than Biermann when the electron mean free path becomes comparable to or greater than the temperature gradient scale length, known as the nonlocal regime. We investigate the nonlocal suppression of the Biermann term using simplified Fokker–Planck simulations covering a wide range of parameters. We provide the first fit for nonlocal Biermann suppression that has physically accurate behavior for small and large values of a suitable nonlocality parameter, valid for an arbitrary atomic number, and that includes the effect of magnetization on nonlocality. The fit is intended to provide an approximate method to account for reduced magnetic-field generation in MHD codes and theory.

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

Experiments on the dynamics and scaling of spontaneous-magnetic-field saturation in laser-produced plasmas.

In laser-produced high-energy-density plasmas, large-scale strong magnetic fields are spontaneously generated by the Biermann battery effects when temperature and density gradients are misaligned. Saturation of the magnetic field takes place when convection and dissipation balance field generation. While theoretical and numerical modeling provide useful insight into the saturation mechanisms, experimental demonstration remains elusive. In this letter, we report an experiment on the saturation dynamics and scaling of Biermann battery magnetic field in the regime where plasma convection dominates. With time-gated charged-particle radiography and time-resolved Thomson scattering, the field structure and evolution as well as corresponding plasma conditions are measured. In these conditions, the spatially resolved magnetic fields are reconstructed, leading to a picture of field saturation with a scaling of B∼1/L_{T} for a convectively dominated plasma, a regime where the temperature gradient scale (L_{T}) exceeds the ion skin depth.

The role of ETG modes in JET–ILW pedestals with varying levels of power and fuelling

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 and normalised density gradients around their experimental values, we demonstrate that it is the ratio of these two quantities 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 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 , it follows a linear trend in correspondence with previous work by different authors.

Reduced models for ETG transport in the tokamak pedestal

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.

Explaining the lack of power degradation of energy confinement in wide pedestal quiescent H-modes via transport modeling

Wide pedestal quiescent H (WPQH)-mode is an attractive scenario for future burning plasmas as they operate without ELMs. WPQH is characterized by formation of a wider and higher pedestal (than quiescent H-mode), and broadband fluctuations in the pedestal. Unlike conventional H-modes, where the energy confinement time reduces with increasing heating power, the WPQH plasmas reported in this paper do not show power degradation of the energy confinement. As the injected neutral beam power was increased, reduced core (ρ ⩽ 0.45) transport calculated by transp, as well as increased core temperatures, pressure gradient and diamagnetic E × B shear rate were observed. The reduction in the heat transport and rapid decrease in the ion temperature gradient scale length suggest the formation of an ion internal transport barrier (ITB) that was accompanied by increased stored energy in the core. Quasilinear turbulent transport modeling using the trapped gyro Landau fluid (tglf) code was used to predict the ITB and its turbulence stability properties. By using profiles and equilibria produced by matching the transp transport fluxes with the tglf transport model within the tgyro transport solver, the energy confinement time captures the experimentally observed insensitivity to the increased P NBI. Linear stability analysis reveals that drift-wave instabilities in the core are stabilized by E × B shear, T i/T e ratio and Shafranov shift; the latter was found to have the strongest effect on the turbulence suppression at the highest heating level.

Fast modulating electron cyclotron emission (FMECE) diagnostic for tokamaks.

Utilizing variable-frequency channels, e.g., yttrium iron garnet (YIG) bandpass filters, in the intermediate frequency (IF) section of an electron cyclotron emission (ECE) radiometer facilitates flexibility in the volume viewed by the ECE channels as well as high resolution electron temperature and temperature fluctuation measurements in tokamaks. Fast modulating electron cyclotron emission (FMECE), a stand-alone IF section with eight channels, is a novel application of YIG filters for real-time electron temperature gradient and gradient scale length measurements. Key to FMECE is a simultaneous input/output data acquisition unit, as well as a modified type of YIG filters, which is capable of fast switching of their center (set) frequencies with a frequency slew rate of 600 µs/GHz. A new FMECE has been implemented and tested on the DIII-D tokamak, demonstrating its capability in real-time gradient measurements. The data presented here shows that FMECE can identify flattening in the electron temperature profile; the latter can be used as a sensor for real time monitoring and control of plasma instabilities. Implementation and application are planned for the EAST tokamak.

Temperature distributions and gradients in laser-heated plasmas relevant to magnetized liner inertial fusion.

We present two-dimensional temperature measurements of magnetized and unmagnetized plasma experiments performed at Z relevant to the preheat stage in magnetized liner inertial fusion. The deuterium gas fill was doped with a trace amount of argon for spectroscopy purposes, and time-integrated spatially resolved spectra and narrow-band images were collected in both experiments. The spectrum and image data were included in two separate multiobjective analysis methods to extract the electron temperature spatial distribution T_{e}(r,z). The results indicate that the magnetic field increases T_{e}, the axial extent of the laser heating, and the magnitude of the radial temperature gradients. Comparisons with simulations reveal that the simulations overpredict the extent of the laser heating and underpredict the temperature. Temperature gradient scale lengths extracted from the measurements also permit an assessment of the importance of nonlocal heat transport.

Propagation of input parameter uncertainties in transport models

The many sources of uncertainty in validation studies of plasma turbulence in magnetically confined fusion devices are well-known. In this paper, we investigate how to efficiently transform uncertainties in experimentally derived transport model inputs into model prediction uncertainties, using the quasilinear trapped-gyro-Landau-fluid (TGLF) turbulent transport model [Staebler et al., Phys. Plasmas 14, 055909 (2007)]. We use the rapidly converging and computationally inexpensive non-intrusive probabilistic collocation method (PCM) to propagate input parameter uncertainty probability distribution functions (PDFs) through TGLF, yielding PDFs of predicted transport fluxes. We observe in many cases that the flux PDFs exhibit significant non-normal features such as strong skewness, even when the input distributions were normal. To illustrate the utility of the PCM approach, we apply this methodology to transport predictions for a DIII-D ITER baseline plasma [Grierson et al., Phys. Plasmas 25, 022509 (2018)] in which the mix of neutral beam injection (NBI) and electron cyclotron heating (ECH) was varied. The model predictions show clear changes in the parametric dependencies and sensitivities of the turbulence between the two heating mixes. Specifically, when only NBI heating was used, the transport fluxes responded significantly only to the ion temperature gradient scale length. However, when both NBI and ECH were applied, the electron transport channels demonstrate a strong sensitivity to the electron temperature and density gradients not observed in the NBI-only case. Additional context for the PCM approach is provided by comparing its predictions with those obtained via a local flux-matching approach. A new set of validation metrics based on the Wasserstein distance is proposed for PDF-based comparisons.

Upgrade of the ECE diagnostic on EAST.

The electron cyclotron emission (ECE) diagnostic on the experimental advanced superconducting tokamak (EAST) was upgraded recently to provide electron temperature profile measurement with wider radial coverage and better precision. The lower limit of the ECE detection frequency band was extended from 104 GHz to 97 GHz by adding a new 8-channel heterodyne radiometer, which ensures capability for the measurement of the second harmonic ECE with toroidal magnetic field down to 1.75 T. Also, the existing 32-channel heterodyne radiometer has been upgraded, with the frequency interval for the lower frequency range up to 120 GHz reduced from 2 GHz to 1 GHz by introducing eight channels in the intermediate frequency part. In addition, a plan is presented to incorporate tunable yttrium iron garnet filters into the existing heterodyne radiometer to obtain detailed measurements of the electron temperature gradient scale length as well as finer spatial pinpointing of magnetohydrodynamic modes. Examples from DIII-D are provided where similar ECE diagnostic allowed precise measurement of the center and width of neoclassical tearing modes.

Electron critical gradient scale length measurements of ICRF heated L-mode plasmas at Alcator C-Mod tokamak

A profile for the critical gradient scale length (Lc) has been measured in L-mode discharges at the Alcator C-Mod tokamak, where electrons were heated by an ion cyclotron range of frequency through minority heating with the intention of simultaneously varying the heat flux and changing the local gradient. The electron temperature gradient scale length (LTe−1 = |∇Te|/Te) profile was measured via the BT-jog technique [Houshmandyar et al., Rev. Sci. Instrum. 87, 11E101 (2016)] and it was compared with electron heat flux from power balance (TRANSP) analysis. The Te profiles were found to be very stiff and already above the critical values, however, the stiffness was found to be reduced near the q = 3/2 surface. The measured Lc profile is in agreement with electron temperature gradient (ETG) models which predict the dependence of Lc−1 on local Zeff, Te/Ti, and the ratio of the magnetic shear to the safety factor. The results from linear Gene gyrokinetic simulations suggest ETG to be the dominant mode of turbulence in the electron scale (k⊥ρs > 1), and ion temperature gradient/trapped electron mode modes in the ion scale (k⊥ρs < 1). The measured Lc profile is in agreement with the profile of ETG critical gradients deduced from Gene simulations.

Properties of ion temperature gradient and trapped electron modes in tokamak plasmas with inverted density profiles

We perform a numerical study of linear stability of the ion temperature gradient (ITG) mode and the trapped electron mode (TEM) in tokamak plasmas with inverted density profiles. A local gyrokinetic integral equation is applied for this study. From comprehensive parametric scans, we obtain stability diagrams for ITG modes and TEMs in terms of density and temperature gradient scale lengths. The results show that, for the inverted density profile, there exists a normalized threshold temperature gradient above which the ITG mode and the TEM are either separately or simultaneously unstable. The instability threshold of the TEM for the inverted density profile is substantially different from that for normal and flat density profiles. In addition, deviations are found on the ITG threshold from an early analytic theory in sheared slab geometry with the adiabatic electron response [T. S. Hahm and W. M. Tang, Phys. Fluids B 1, 1185 (1989)]. A possible implication of this work on particle transport in pellet fueled tokamak plasmas is discussed.

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

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.

Temperature gradient scale length measurement: A high accuracy application of electron cyclotron emission without calibration.

Calibration is a crucial procedure in electron temperature (Te) inference from a typical electron cyclotron emission (ECE) diagnostic on tokamaks. Although the calibration provides an important multiplying factor for an individual ECE channel, the parameter ΔTe/Te is independent of any calibration. Since an ECE channel measures the cyclotron emission for a particular flux surface, a non-perturbing change in toroidal magnetic field changes the view of that channel. Hence the calibration-free parameter is a measure of Te gradient. BT-jog technique is presented here which employs the parameter and the raw ECE signals for direct measurement of electron temperature gradient scale length.

Electron temperature critical gradient and transport stiffness in DIII-D

In a continuing effort to validate turbulent transport models, the electron energy flux has been probed as a function of electron temperature gradient on the DIII-D tokamak. In the scan of gradient, a critical electron temperature gradient has been found in the electron heat fluxes and stiffness at various radii in L-mode plasmas. The TGLF reduced turbulent transport model (Staebler et al 2007 Phys. Plasmas 14 055909) and full gyrokinetic GYRO model (Candy and Waltz 2003 J. Comput. Phys. 186 545) recover the general trend of increasing electron energy flux with increasing electron temperature gradient scale length, but they do not predict the absolute level of transport at all radii and gradients. Comparing the experimental observations of incremental (heat pulse) diffusivity and stiffness to the models’ reveals that TGLF reproduces the trends in increasing diffusivity and stiffness with increasing electron temperature gradient scale length with a critical gradient behavior. The critical gradient of TGLF is found to have a dependence on q95, contrary to the independence of the experimental critical gradient from q95.

How to apply a turbulent transport model based on a gyrokinetic simulation for the ion temperature gradient mode in helical plasmas

How to apply a reduced model for the turbulent ion heat diffusivity [Nunami M et al. 2013 Phys. Plasmas, 20 092307] derived from a gyrokinetic code to a transport simulation is proposed. The reduced model is given by the function of the linear growth rate of the ion temperature gradient (ITG) mode and the decay time of zonal flows. The ion temperature gradient scale length is chosen for the additional modeling to include the linear growth rate of the ITG mode from a gyrokinetic code. The formula for the zonal flow decay time is derived at the given magnetic field configuration. The calculation of an extremely low computational cost in this article reproduces the results of the reduced model for the turbulent ion heat diffusivity within allowable errors.

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.

Hydrogen isotope effects on ITG scale length, pedestal and confinement in JT-60 H-mode plasmas

The dependence of heat transport, edge pedestal and confinement on isotopic composition was investigated in conventional H-mode plasmas. Identical profiles for the electron density, electron temperature and ion temperature were obtained for hydrogen and deuterium plasmas, whereas the required power clearly increased for hydrogen, which resulted in a reduction in heat diffusivity for deuterium. The determination of identical temperature profiles, despite the different heating power, suggested that the characteristics of heat conduction essentially differ for hydrogen and deuterium, even at the same scale length of temperature gradient. The self-regulating physics mechanism determining the overall H-mode confinement was also addressed. The inverse of the ion temperature gradient (ITG) scale length, or ∇Ti/Ti, which is required for a given ion heat diffusivity, increased by a factor of approximately 1.2 for deuterium compared with that for hydrogen. The relationship between edge pedestal pressure and global βp holds true consistently regardless of the difference in the isotopic composition. A higher value of βp was obtained for deuterium because of its smaller ITG scale length and because of the additional stored energy in the thermal and fast ion components, the latter due to an increase in the slowing down time with an increase in isotopic mass.

Gyrokinetic Theory of Turbulent Acceleration of Parallel Rotation in Tokamak Plasmas

A mechanism for turbulent acceleration of parallel rotation is discovered using gyrokinetic theory. This new turbulent acceleration term cannot be written as a divergence of parallel Reynolds stress. Therefore, turbulent acceleration acts as a local source or sink of parallel rotation. The physics of turbulent acceleration is intrinsically different from the Reynolds stress. For symmetry breaking by positive intensity gradient, a positive turbulent acceleration, i.e., cocurrent rotation, is predicted. The turbulent acceleration is independent of mean rotation and mean rotation gradient, and so constitutes a new candidate for the origin of spontaneous rotation. A quasilinear estimate for ion temperature gradient turbulence shows that the turbulent acceleration of parallel rotation is explicitly linked to the ion temperature gradient scale length and temperature ratio Ti0/Te0. Methods for testing the effects of turbulent parallel acceleration by gyrokinetic simulation and experiment are proposed.