Search NASA⌕ Search

SEARCH · Search NASA

Results for “TGLF”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Successful prediction of tokamak transport in the L-mode regime

A long standing shortfall in the predicted L-mode edge energy transport by reduced quasi-linear models of gyrokinetic turbulent transport has been resolved. The improved model TGLF-SAT2 has higher fidelity to gyrokinetic simulations of the electron-scale contribution to the electron energy transport and the ion-scale flux surface shape dependence of energy transport. The success of TGLF-SAT2 in predicting the L-mode and Ohmic edge profiles is critical to whole pulse simulation and opens the door to prediction of the H-mode power threshold.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Validation of energetic particle transport models for time-dependent integrated simulations of burning plasmas (Final ST Report)

The major goal of this project is to deliver a suite of validated energetic particle (EP) modeling tools for time-dependent integrated tokamak simulations (with the code TRANSP) that can be used for modeling plasma discharges featuring a substantial EP content, such as burning plasmas. During FY20, issues were resolved with an energetic particle transport model called TGLF-EP, with improvements made to the TGLF-EP workflow, generation of input files, and interpretation of runs on the NERSC supercomputing cluster. The recipient also gave input on simulations with the code TRANSP to reproduce the reference JET discharge, and assessed the accuracy and uncertainties of the simulated B field, current evolution, thermal plasma profiles, stored energy, etc. Input was also given on interpretation of fast-ion instability calculations and recommendations for modifications to JET afterglow discharge for the upcoming DT campaign.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Multi-scale transport in the DIII-D ITER baseline scenario with direct electron heating and projection to ITER

Multi-scale fluctuations measured by turbulence diagnostics spanning long and short wavelength spatial scales impact energy confinement and the scale-lengths of plasma kinetic profiles in the DIII-D ITER baseline scenario with direct electron heating. Contrasting discharge phases with ECH + neutral beam injection (NBI) and NBI only at similar rotation reveal higher energy confinement and lower fluctuations when only NBI heating is used. Modeling of the core transport with TGYRO using the TGLF turbulent transport model and NEO neoclassical transport reproduces the experimental profile changes upon application of direct electron heating and indicates that multi-scale transport mechanisms are responsible for changes in the temperature and density profiles. Intermediate and high-k fluctuations appear responsible for the enhanced electron thermal flux, and intermediate-k electron modes produce an inward particle pinch that increases the inverse density scale length. Projection to ITER is performed with TGLF and indicates a density profile that has a finite scale length due to intermediate-k electron modes at low collisionality and increases the fusion gain. For a range of E×B shear, the dominant mechanism that increases fusion performance is suppression of outward low-k particle flux and increased density peaking.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Predictions of core plasma performance for the SPARC tokamak

SPARC is designed to be a high-field, medium-size tokamak aimed at achieving net energy gain with Ion Cyclotron Range-of-Frequencies (ICRF) as its primary auxiliary heating mechanism. Empirical predictions with conservative physics indicate that SPARC baseline plasmas would reach Q~11, well above its mission objective of Q>2. To build confidence that SPARC will be successful, physics-based integrated modeling has also been performed. The TRANSP code coupled with the theory-based TGLF turbulence model and EPED predictions for pedestal stability find that Q~9 is attainable in standard H-mode operation and confirms Q>2 operation is feasible even with adverse assumptions. In this analysis, ion cyclotron waves are simulated with the full wave TORIC code and alpha heating is modeled with the Monte-Carlo fast ion NUBEAM module. Detailed analysis of expected turbulence regimes with linear and nonlinear CGYRO simulations is also presented, demonstrating that profile predictions with the TGLF reduced model are in reasonable agreement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Equilibrium reconstruction of DIII-D plasmas using predictive modeling of the pressure profile

New workflows have been developed for predictive modeling of magnetohydrodynamic (MHD) equilibrium in tokamak plasmas. The goal of this work is to predict the MHD equilibrium in tokamak discharges without having measurements of the kinetic profiles. The workflows include a cold start tool, which constructs all the profiles and power flows needed by transport codes; a Grad–Shafranov equilibrium solver; and various codes for the sources and sinks. For validation purposes, a database of DIII-D tokamak discharges has been constructed that is comprised of scans in the plasma current, toroidal magnetic field, and triangularity. Initial efforts focused on developing a workflow utilizing an empirically derived pressure model tuned to DIII-D discharges with monotonic safety factor profiles. This workflow shows good agreement with experimental kinetic equilibrium calculations, but is limited in that it is a single fluid (equal ion and electron temperatures) model and lacks H-mode pedestal predictions. The best agreement with the H-mode database is obtained using a theory-based workflow utilizing pressure profile predictions from a coupled TGLF turbulent transport and EPED pedestal models together with external magnetics and Motional Stark Effect (MSE) data to construct the equilibrium. Here, we obtain an average root mean square error of 5.1% in the safety factor profile when comparing the predicted and experimental kinetic equilibrium. We also find good agreement with the plasma stored energy, internal inductance, and pressure profiles. Including MSE data in the theory-based workflow results in noticeably improved agreement with the q-profiles in high triangularity discharges in comparison with the results obtained with magnetic data only. The predictive equilibrium workflow is expected to have wide applications in experimental planning, between-shot analysis, and reactor studies.

Kinsey, J. E. (ORCID:0000000193347473)↗
DOE OSTI2024

Core performance predictions in projected SPARC first-campaign plasmas with nonlinear CGYRO

This work characterizes the core transport physics of SPARC early-campaign plasmas using the PORTALS-CGYRO framework. Empirical modeling of SPARC plasmas with L-mode confinement indicates an ample window of breakeven (Q > 1) without the need of H-mode operation. Extensive modeling of multi-channel (electron energy, ion energy, and electron particle) flux-matched conditions with the nonlinear CGYRO code for turbulent transport coupled to the macroscopic plasma evolution using PORTALS reveals that the maximum fusion performance to be attained will be highly dependent on the near-edge pressure. Stiff core transport conditions are found, particularly when fusion gain approaches unity, and predicted density peaking is found to be in line with empirical databases of particle source-free H-modes. Impurity optimization is identified as a potential avenue to increase fusion performance while enabling core-edge integration. Extensive validation of the quasilinear TGLF model builds confidence in reduced-model predictions. The implications of projecting L-mode performance to high-performance and burning-plasma devices is discussed, together with the importance of predicting edge conditions.

Rodriguez-Fernandez, P. (ORCID:0000000273611131)↗
DOE OSTI2025

Molybdenum density profiles in Alcator C-Mod plasmas

Molybdenum (Z = 42) density spatial profiles have been obtained for Alcator C-Mod L-, I-, and H-mode plasmas from x-ray observations of the neon-like Mo 32+ 4D transition. Profile shapes range from flat to modestly peaked in the core. The observed peaking factors for molybdenum are in good agreement with those calculated by TGLF across all confinement regimes, and there is little dependence seen with core density gradient, sawtooth amplitude, or rotation velocity. Predictions for tungsten peaking in SPARC are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Electron ITB formation in EAST high poloidal beta plasmas under dominant electron heating

Abstract Plasma confinement and transport in tokamaks play a crucial role in the development of high poloidal beta steady-state operation scenarios. Therefore, it is very important to study the relevant mechanism of internal transport barriers (ITBs), which can help plasmas to obtain better confinement in order to achieve higher fusion gain. This paper mainly introduces the analysis of the characteristics of electron heat transport of discharges with ITB in high β P operation regime in Experimental Advanced Superconducting Tokamak (EAST). Based on the statistical analysis of stable discharges with β P > 1.5, it is found that there is an obvious bifurcation of the normalised electron temperature gradient (ETG) ( R / L T e ) in the range of β P = 2–2.2. Then the discharges of lower β P ( β P < 2, where the value of β P is below the bifurcation threshold) and of higher β P ( β P > 2.2, where the value of β P is above the threshold) were selected for analysis. The diagnostic data provided by Thomson scattering, x-ray crystal spectrometry and charge exchange recombination spectroscopy are used to provide reliable parameter profiles and then combined with the data of external magnetic probe measurements and the polarisation interferometer diagnosis system to fully reconstruct the balance. A relevant plasma current calculation model is used to calculate and analyse the current density profiles and power deposition, and then the transport analysis is carried out. Interestingly, in the higher β P discharges, it is found that the turbulence intensity provided by the CO 2 laser collective scattering system gradually decreases and the normalised ETG R / L T e gradually increases. It is also found that the electron heat transport coefficient decreases in the discharge with higher β P and the growth rate of the electron-scale turbulence calculated by transport gyro-Landau fluid (TGLF) is significantly reduced. Meanwhile, similar conclusions are also obtained in the discharges when the β P is further increased.

Wang, Z. H. (ORCID:0009000203431036)↗
DOE OSTI2025

Role of perturbed parallel magnetic field effects in predicting turbulent transport in NSTX

This study presents analysis of gyrokinetic simulations on the National Spherical Torus Experiment (NSTX) to investigate the effects of electromagnetic fields on plasma turbulence and transport. The simulations, performed with varying levels of fidelity using the gyrokinetic CGYRO code, include electrostatic (ES), single-field electromagnetic (EM1), and two-field electromagnetic (EM2) models. A detailed comparison across the simulation database reveals that electromagnetic effects increase both predicted growth rates and quasilinear fluxes, with EM2 simulations producing stronger turbulence than ES and EM1 cases. Quasilinear modeling using QLGYRO demonstrates that while the perturbed parallel magnetic field (δB ∥ ) does not drastically affect the total flux at experimental gradients, it leads to a shift in the dominant instability, altering mode structures from microtearing to kinetic ballooning modes (KBMs). The proximity of the plasma profiles to the KBM threshold is explored, with the experimental conditions being near the onset of KBM-driven transport. The KBM, with its large growth rates, is identified as a potential driver of electron temperature flattening, as it can rapidly transport heat across flux surfaces. Performing stability analysis shows core-localized unstable a low- mode that could contribute to the flattening at the early times of the discharge. TGYRO predictive modeling, incorporating both TGLF and QLGYRO, indicates that the inclusion of δB ∥ significantly improves the accuracy of temperature profile predictions in NSTX high-beta plasmas, although challenges remain in modeling the sharp flux discontinuities caused by KBM-driven instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Sensitivities of time-dependent temperature profile predictions for NSTX with the multi-mode model

The Multi-Mode Model (MMM) for turbulent transport was applied to a large set of well-analyzed discharges from the National Spherical Torus Experiment (NSTX) in order to evaluate its sensitivities to a wide range of plasma conditions. MMM calculations were performed for hundreds of milliseconds in each discharge by performing time-dependent predictive simulations with the 1.5D tokamak integrated modeling code TRANSP. A closely related study (Lestz et al 2025 Plasma Phys. Control. Fusion 67 105029) concluded that MMM predicted electron and ion temperature profiles that were in reasonable agreement with NSTX observations, generally outperforming a different reduced transport model, TGLF. This finding motivates the more thorough investigation of the characteristics of the MMM predictions conducted in this work. The simulations with MMM have electron energy transport dominated by electron temperature gradient modes in the examined discharges with relatively low plasma β (ratio of kinetic plasma pressure to magnetic field pressure) and high collisionality, transitioning to a mixture of different modes for higher β and lower collisionality. The thermal ion diffusivity predicted by MMM is much smaller than the neoclassical contribution, in line with previous experimental analysis of NSTX. Nonetheless, the electron and ion temperature profiles are coupled via collisional energy exchange and thus sensitive to which transport channels are predicted. The time-dependent simulations with MMM are robust to the simulation start time, converging to remarkably similar temperature profiles later during the discharge. MMM typically overpredicts confinement relative to NSTX observations, leading to the prediction of overly steep temperature profiles. Plasmas with spatially broader temperature profiles, higher plasma β, and longer energy confinement times tend to be predicted by MMM with better agreement with the experiment. As a result, these findings provide useful context for understanding the regime-dependent tendencies of MMM in anticipation of self-consistent, time-dependent predictive simulations of NSTX-U discharges with these same modeling tools.

MMM↗
DOE OSTI2020

Predicting the rotation profile in ITER

Determining the toroidal rotation for future tokamaks like ITER is a challenging and important problem. By combining empirical scalings for the intrinsic rotation at the top of the pedestal with the expected neutral beam torque and modeling of momentum transport, the toroidal rotation profile for ITER is predicted with TGYRO using TGLF (SAT0 and SAT1). On axis rotation exceeds 20 krad s –1 and the E x B shear is significant enough to reduce turbulent transport and significantly increase confinement and fusion power when comparing to cases that ignore the effect of rotation. Here, the prediction of the rotation at the top of the pedestal is made with increased confidence due to experiments and modeling in DIII-D that have determined the importance of fast-ion and neutral particle transport effects on intrinsic rotation at this location. In particular, the effect of neutral particles on momentum transport in the pedestal region is found to be insignificant.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Collisionality driven turbulent particle transport changes in DIII-D H-mode plasmas

The results of the experimental dimensionless scan in this paper confirm that there is an increase in density peaking towards lower collisionality and that this can be partly linked to a shift in the turbulence regime from ITG towards TEM. However at the lowest collisionality, the changes in turbulence and transport are much more pronounced than expected from direct collisionality effect on the turbulence. In this paper, the collisionality, ν * is varied by a factor 5, while keeping ρ * , q , β , M , fixed. Additionally, a 3 Hz gas puff modulation is applied to modulate the electron density profile and extract the perturbed transport coefficients using two diagnostics. The transport analysis shows that the increase in density peaking at low ν * is linked to an increase in the inward particle pinch and not an increase in core fueling. These observations are not only in agreement with prior modeling scans of how turbulence changes as a function of collisionality and its impact upon the particle fluxes, but also with the multi-machine database (Fable E. et al 2010 Plasma Phys. Control. Fusion 52 015007) (Angioni C. et al 2003 Phys. Rev. Lett. 90 205003). The changes in turbulence across the collisionality scan were captured at large scale by the BES and at smaller scale by the DBS. A comparison with gradient-driven GENE simulations showed similar trends at both scales. Moreover, the changes observed in overall transport are in agreement with gradient-driven TGLF particle flux simulations. This indicates that TGLF/GENE when given the gradients as input, are able to reproduce the experimentally observed turbulence changes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Examination of stiff ion temperature gradient mode physics in simulations of DIII-D H-mode transport

This paper presents comparisons of predictions from two models of turbulent transport in tokamaks. The focus is on assessing how close a reduced but significantly less computationally demanding model (TGLF) does in matching the more detailed model (CGYRO) in predicting scalings with driving gradient ("stiffness"). The impacts of other physics mechanisms such as flow shear, fast ions, and electromagnetic fluctuations on these predictions is also studied. Overall, while there are clear differences in the scaling, these do not appear to significantly impact the fidelity of the reduced model in prediction the experimental conditions studied relative to the more advanced model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Verification of a quasi-linear model for gyrokinetic turbulent transport

The verification and calibration of a new quasi-linear transport model with a large database of gyrokinetic turbulence simulations is presented in this paper. In a previous paper (Staebler et al 2020 Plasma Phys. Control. Fusion 63 015013), a model for the saturated spectrum of electric potential fluctuations was developed based on the properties of the non-linear 3D spectrum. In this paper, a modification to the overall multiplicative factor in this model is found to be necessary to improve the fit to scans of the temperature and density gradients. The error in the fit of the quasi-linear fluxes of electron and ion energy fluxes is significantly better than for previous saturation models. The spectral shift model for the impact of equilibrium E × B velocity shear (Staebler et al 2013 Phys. Rev. Lett. 110 055003) and the zonal flow mixing model for electron-scale turbulence (Staebler et al 2016 Phys. Plasmas 23 062518) are both revised to be compatible with this new model. Furthermore, the models for the loss of bounce averaging and electron collisions in the TGLF reduced linear equations (Staebler et al 2005 Phys. Plasmas 12 102508) are also changed to improve the linear eigenmodes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Experimental inference of neutral and impurity transport in Alcator C-Mod using high-resolution x-ray and ultra-violet spectra

We present experimental inferences of cross-field impurity transport coefficients for AlcatorC-Mod plasmas without edge-localized modes, leveraging a novel forward model for the entire Ca Kα spectrum, including satellite lines within the spectral range, to compare to high-resolution x-ray imaging crystal spectroscopy (XICS). These measurements are complemented by extreme ultra-violet (EUV) spectroscopy that constrains transport closer to the edge. Using new atomic data sets for both XICS and EUV analysis has enabled consideration of line ratios across both spectral ranges and has increased the accuracy of inferred transport coefficients. Inclusion of charge exchange between edge thermal neutrals and impurities is shown to be extremely important in C-Mod pedestals. Atomic D neutral densities from experimental D Lyα measurements at the midplane are compared to SOLPS-ITER simulations, finding good agreement. Bayesian inferences of impurity transport coefficients are presented for L-, EDA H-, and I-mode discharges, making use of the Aurora package for forward modeling and combining our spectroscopic constraints. Experimentally inferred diffusion profiles are found to match turbulent transport models at midradius within uncertainties, using both quasilinear gyro-fluid TGLF SAT-1 and nonlinear ion-scale gyrokinetic CGYRO simulations. Furthermore, significant discrepancies in convection are observed in some cases, suggesting difficulties in predictions of flat or hollow impurity profiles.

neutrals↗
DOE OSTI2021

The dependence of confinement on the isotope mass in the core and the edge of AUG and JET-ILW H-mode plasmas

Experiments in ASDEX Upgrade (AUG) and JET with the ITER-like wall (JET-ILW) are performed to separate the pedestal and core contributions to confinement in H-modes with different main ion masses. A strong isotope mass dependence in the pedestal is found which is enhanced at high gas puffing. This is because the ELM type changes when going from D to H for matched engineering parameters, which is likely due to differences in the inter ELM transport with isotope mass. The pedestal can be matched in H and D plasmas by varying only the triangularity and keeping the engineering parameters relevant for core transport the same. With matched pedestals Astra/TGLF (Sat1geo) core transport simulations predict the experimental profiles equally well for H and D. These core transport simulations show a negligible mass dependence and no gyro-Bohm scaling is observed. However, to match the experimental observations at medium β it is required to take the fast-ion dilution and rotation into account. This is not enough for high β plasmas where for the first time a profile match between H and D plasmas was achieved experimentally. Under these conditions quasilinear modelling with TGLF over predicts the transport in the core of H and D plasmas alike.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Global gyrokinetic simulations of ASDEX Upgrade up to the transport timescale with GENE–Tango

An accurate description of turbulence up to the transport timescale is essential for predicting core plasma profiles and enabling reliable calculations for designing advanced scenarios and future devices. Here, we exploit the gap separation between turbulence and transport timescales and couple the global gyrokinetic code GENE to the transport-solver Tango, including kinetic electrons, collisions, realistic geometries, toroidal rotation and electromagnetic effects for the first time. This approach overcomes gyrokinetic codes' limitations and enables high-fidelity profile calculations in experimentally relevant plasma conditions, significantly reducing the computational cost. We present numerical results of GENE–Tango for two ASDEX Upgrade discharges, one of which exhibits a pronounced peaking of the ion temperature profile not reproduced by TGLF–ASTRA. We show that GENE–Tango can correctly capture the ion temperature peaking observed in the experiment. By retaining different physical effects in the GENE simulations, e.g., collisions, toroidal rotation and electromagnetic effects, we show that the ion temperature profile's peaking can be linked to electromagnetic effects of submarginal (stable) KBM modes. Furthermore, based on these results, the expected GENE–Tango speedup for the ITER standard scenario is larger than two orders of magnitude compared to a single gyrokinetic simulation up to the transport timescale, possibly making first-principles ITER simulations feasible on current computing resources.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Isotope mass scaling and transport comparison between JET Deuterium and Tritium L-mode plasmas

The dimensionless isotope mass scaling experiment between pure Deuterium and pure Tritium plasmas with matched $\rho^*$, $\nu^*$ , β n , q and $T_e/T_i$ has been achieved in JET L-mode with dominant electron heating (NBI+ohmic) conditions. 28% higher scaled energy confinement time $B_t\tau_{E,\mathrm{th}}/A$ is found in favour of the Tritium plasma. This can be cast in the form of the dimensionless energy confinement scaling law as $\Omega_i \tau_{E,\mathrm{th}} \sim A^{0.48 \pm 0.16}$. This significant isotope mass scaling is consequently seen in the scaled one-fluid heat diffusion coefficient $A \chi_{\mathrm{eff}}/B_t$ which is around 50% lower in the Tritium plasma throughout the whole plasma radius. The isotope mass dependence in the particle transport channel is negligible, supported also by the perturbative particle transport analysis with gas puff modulation. The comparison of the edge particle fuelling or ionisation profiles from the EDGE2D-EIRENE simulations show that the absolute density differences that are necessary for the dimensionless match in the confined plasma dominate over any isotope mass dependencies of particle fuelling and ionization profiles at the plasma edge. Local GENE simulation results indicate a mild anti-gyroBohm effect at $\rho_\mathrm{tor}$ = 0.6 and thereby a small isotope mass dependence in favour of Tritium on heat transport and a negligible effect on particle transport. A significant fraction of the isotope scaling and reduced heat transport observed in the Tritium plasma is not captured in the GENE and ASTRA-TGLF-SAT2 simulations by simply changing the isotope mass for the same input profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗