Engineering Papers⌕ Search

SEARCH · Engineering Papers

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.

13 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

TRANSP-TGLF core predictive modeling of the JET DT baseline scenario

In recent years, an intense modeling activity has been focused on preparing and analyzing the second JET Deuterium–Tritium (D–T) experimental campaign DTE2. Among the numerous scientific outcomes of this campaign was the unique opportunity to test and validate the state of the art of modeling tools with fusion-relevant DT plasmas using the full metallic ITER-like wall in different scenarios. This work reports on the core predictive modeling of plasma density, electron and ion temperatures (n e , T e , T i ) performed using TRANSP (Pankin et al 2025 Comput. Phys. Commun. 312 109611) coupled with Trapped Gyro-Landau Fluid (TGLF)-SAT2 (Kinsey et al 2008 Phys. Plasmas 15 055908), (Staebler et al 2021 Nucl. Fusion 61 116007) for the JET D–T baseline scenario (I p = 3.5 MA, q 95 = 3, β N < 2, with pellet pacing) (Garzotti et al 2025 Plasma Phys. Control. Fusion 67 075011). The sensitivity to different input parameters, as the $\vec{E}$ x $\vec{B}$ shear parameterization and the values of the kinetic quantities at the boundary of the prediction domain (ρ = 0.85) has been assessed, identifying the confidence interval of the prediction results. In particular, the dependence of the electron density profile on the particle source parameters has been studied, identifying the ionization source as the main cause for the density gradient under-prediction obtained by TRANSP-TGLF, reported in (Hyun-Tae et al 2023 Nucl. Fusion 63 112004).

JET↗
DOE OSTI2026

Benchmarking core turbulence and transport predictions for an inductive compact tokamak reactor plasma

Motivated by the need for accurate, timely, and efficient calculations of plasma transport, predictions of plasma turbulence properties made using different TGLF saturation rules are benchmarked against corresponding predictions from linear and nonlinear gyrokinetic CGYRO simulations. This benchmarking is carried out using parameters taken from an inductive burning plasma scenario in a hypothetical compact high-field (R maj = 4 m, B T = 8 T) tokamak, lying in a much different regime of parameter space than either the TGLF calibration regime or current-day experiments. The core turbulent transport in this scenario is predicted to be dominated by ion temperature gradient (ITG) turbulence. In general, the ITG critical gradients predicted by various TGLF saturation rules are quite close to the CGYRO predictions. Both codes predict similar linear ITG growth rates and frequency spectra, as well as their scaling with R/L T i = −Rd ln(T i )/dr. However, TGLF systematically predicts unstable trapped-electron modes (TEMs) above k y ρ s ≃ 0.5 not seen by CGYRO for the same parameters, due to TGLF predicting a lower threshold in R/L T e than CGYRO for TEM onset. It is shown that for this scenario, nonlinear CGYRO simulations predict stiffer ITG turbulence than the TGLF SAT0 and SAT1 saturation rules, with energy fluxes close in magnitude and scaling with R/L T i to what is predicted by the SAT2 saturation rule. Self-consistent core profiles calculated using nonlinear CGYRO flux predictions and the PORTALS transport solver are shown to agree fairly well with corresponding predictions made using the TGLF SAT2 model, including a similar level of density peaking.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Assessing time-dependent temperature profile predictions using reduced transport models for high performing NSTX plasmas

Time-dependent, predictive simulations were performed with the 1.5D tokamak integrated modeling code TRANSP on a large set of well-analyzed, high performing discharges from the National Spherical Torus Experiment (NSTX) in order to evaluate how well modern reduced transport models can reproduce experimentally observed temperature profiles in spherical tokamaks. Overall, it is found that simulations using the Multi-Mode Model (MMM) more consistently agree with the NSTX observations than those using the Trapped Gyro-Landau Fluid (TGLF) model, despite TGLF requiring orders of magnitude greater computational cost. When considering all examined discharges, MMM has median overpredictions of electron temperature (T e ) and ion temperature (T i ) profiles of 28% and 27%, respectively, relative to the experiment. TGLF overpredicts T e by 46%, with much larger variance than MMM, and underpredicts T i by 25%. As the ratio of kinetic to magnetic field pressure (β) is increased across NSTX discharges, TGLF predicts lower T e and significant flattening of the T i profile, conflicting with NSTX observations. When using an electrostatic version of TGLF, both T e and T i are substantially overpredicted, underscoring the importance of electromagnetic turbulence in the high β spherical tokamak regime. Additionally, calculations with neural net surrogate models for TGLF were performed outside of TRANSP with a time slice flux matching transport solver, finding better agreement with experiment than the TRANSP simulations, highlighting the impact of different transport solvers and simulation techniques. Furthermore, the reasonable agreement with experiment of temperature profiles predicted by MMM motivates a more detailed examination of the sensitivities of the TRANSP simulations with MMM to different NSTX plasma regimes in a companion paper (Lestz et al Plasma Phys. Control. Fusion 67 105030), in preparation for self-consistent, time-dependent predictive modeling of NSTX-U scenarios.

Lestz, J. B. [General Atomics, San Diego, CA (Unit↗
DOE OSTI2025

Analysis of the impact of parallel magnetic fluctuations on linear gyrokinetic stability in NSTX-U and verification of gyro-fluid models

In this work, we use the CGYRO gyrokinetic code to analyze two L- and one H-mode discharges from the National Spherical Torus Experiment (NSTX) and NSTX-Upgrade (NSTX-U) selected due to their different mix of ion-scale driftwaves, ion temperature gradient (ITG) mode and trapped electron mode (TEM), and electromagnetic instabilities, kinetic ballooning mode (KBM), and micro-tearing mode (MTM) in the plasma core. It is found that the effect of parallel magnetic fluctuations is strongly destabilizing to the unstable KBMs compared to calculations with only perpendicular magnetic fluctuations. Two discharges have a mix of ITG/TEM and MTMs that are predicted to be dominant instability across the plasma radius. The parallel magnetic fluctuations are found to have little effect on the MTM stability but are destabilizing to ITG/TEM modes. To test the validity of the gyro-fluid linear stability codes TGLF and GFS at low aspect ratio, a database of linear growth rates has been created using the CGYRO gyrokinetic code. The database is comprised of various parameter scans around a standardized set of NSTX-U core parameters. It contains a group of electrostatic cases and an electromagnetic group that includes the effects of perpendicular and parallel magnetic fluctuations. Comparing the results from the GFS and TGLF models, we find that GFS exhibits the best agreement with the database of CGYRO linear growth rates. Comparing the model results for the electromagnetic scans shows that GFS captures the effects of parallel magnetic fluctuations accurately, while the TGLF model does not, as it lacks sufficient perpendicular energy resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Enhanced shear stabilization of turbulence in NSTX

In studying a particular non-stationary NSTX L-mode plasma, we observed unexpectedly high levels of flux—first with quasilinear (TGLF) modeling and subsequently with nonlinear gyrokinetic simulations. Upon more detailed analysis, a novel confinement regime was discovered in which a modest increase in E x B shear (beyond baseline experimental estimates) rapidly reduced turbulent transport to levels consistent with power balance. This modest increase is plausible given the errors inherent to the estimation of shearing rates, and the added complexity of the non-stationary (time-dependent) power balance. Remarkably, an additional small increase in shear yields the familiar ion-neoclassical transport level with what appears to be the onset of high-k electron transport only. Although analyses using the TGLF-SAT2 model successfully capture numerous parametric dependencies of this plasma, TGLF does not reproduce the rapid E x B stabilization seen in CGYRO. We believe the results presented should help to better characterize the nonlinear physics of spherical tokamak confinement regimes, provide useful ST datasets for reduced model development, and motivate more accurate experimental diagnosis of E x B shearing rates.

Atomic and molecular collisions↗
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 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 OSTI2025

WEST L-mode record long pulses guided by predictions using Integrated Modeling

A new record was set on the WEST Tokamak, designed to operate long duration plasmas in a tungsten (W) environment, with an injected energy of 1.15 GJ and a plasma duration 364s. Scenario development was supported by integrated modeling using the High Fidelity Plasma Simulator (HFPS), the European IMAS-coupled version of JETTO/JINTRAC, which integrates physics-driven modules into a unified framework. In particular, a reduced model for Lower-Hybrid heating and Current-Drive (LHCD) and the quasi-linear turbulent transport model TGLF are crucial for long pulses predictions up to the Last Closed Flux Surface (LCFS). Using this workflow, a 100 s reference discharge was modeled and plasma kinetic profiles and loop voltage were quantitatively well matched. In preparation for the recent long duration experiments, non-inductive current-drive actuators (I P , n e , P LHCD ) were varied to determine the operational domain going towards fully non-inductive discharges. In particular, decreasing the plasma current is shown to ease the access to such conditions, with a careful monitoring of (n e , P LHCD ) to avoid machine limitations. In addition, post-prediction experiments conducted within the investigated parameter range validated the predicted dependencies and were shown to be in quantitative agreement. Exploratory work on the use of ECCD for MHD stability purpose is also introduced.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Shafranov-shift destabilization of ballooning-type micro-instabilities

Contrary to common wisdom, we argue that the Shafranov shift is destabilizing for the ambient ballooning-type instabilities, which account for most of the transport flux in tokamak core plasmas. Higher Shafranov shift indeed reduces the magnetic drift frequency ω d , i.e. improves the bad curvature mildly, around the outboard midplane. However, this improvement is limited to a very narrow spatial region while ω d is increased over the remaining poloidal space. The eigenfunction averaged $\langle$ω d $\rangle$ is effectively enhanced due to the finite mode width of the eigenfunction and is thus destabilizing to the ballooning-type mode, as demonstrated by gyrokinetic simulations with the CGYRO code using local Miller equilibrium geometry. The predicted nonlinear flux also increases with Shafranov shift, consistent with linear simulations. The reduced transport model TGLF can capture the physics reasonably well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

On FIRE mode in KSTAR

We report on the status of Fast Ion Regulated Enhancement (FIRE) mode experiments in the Korea Superconducting Tokamak Advanced Research. This regime is being developed for high-performance, steady-state operation which features a stationary ion internal transport barrier, enabling a central ion temperature approaching 10 keV to be sustained for up to 50 s, without the need for delicate profile control and with no significant impurity accumulation. As its key novelty lies in the significant contribution of fast ions that stabilize core turbulence, the regime has been named FIRE mode. To achieve this regime, neutral beam injection is applied at moderate power levels near the L–H power threshold, while maintaining low plasma density to avoid the L–H transition. The scenario is typically established in diverted magnetic configurations. The core features of FIRE mode were investigated through power balance analysis and fluctuation measurements, revealing a clear transport bifurcation in the ion channel. At the plasma edge, FIRE mode occasionally exhibits I-mode characteristics, particularly in unfavorable magnetic null configurations with q 95 ∼ 4, including the presence of weakly coherent modes. In terms of MHD activity, sawtooth oscillations are observed but appear to be stabilized during the high-performance phase. Fast ion-driven Alfvénic eigenmodes (AEs), indicated by strong frequency chirping near 200 kHz, are also observed. Additionally, lower-frequency MHD activities, distinct from the fast ion-driven AEs, are present and have some influence on plasma performance. The enhancement of core confinement is primarily attributed to fast ion effects, which were evaluated with respect to dilution, alpha stabilization, and resonant interactions with turbulence using gyrokinetic analyses. Among these, the dilution effect was found to be the most dominant. The characteristics of FIRE mode were compared with those of other hot ion plasma scenarios, such as supershot and hot ion mode. While they share many similarities, FIRE mode is distinguished by the accessibility to conditions with T i ≈ T e , presence of I-mode edge features, and its long-duration sustainment. Predictive simulations of FIRE mode were performed using integrated transport modeling with TRIASSIC, employing the TGLF anomalous transport model. These simulations confirmed the critical role of fast ions in achieving this regime. The future prospect of FIRE mode for application in fusion reactors is discussed, with an emphasis on possibly extending the regime to higher density operation.

FIRE mode↗
DOE OSTI2025

FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment (Final Report)

This report documents the FY25 Theory and Simulation Performance Target (TSPT) of developing an integrated modeling framework for fusion reactor design and assessment (FREDA). Over Q1-Q4, new capabilities were developed across both plasma and engineering domains and demonstrated on an example representation of a Compact Advanced Tokamak with a Dual Cooled Lead Lithium blanket. This represents a first-of-a-kind demonstration of coupled core-to-wall-to-engineering for a reactor. Self-consistent CESOL workflows were applied to provide core, pedestal, and SOL prediction; new modules were developed for energetic particle stability (FAR3D) and transport (TGLF-EP) analysis; and boundary plasma modeling (SOLPS-ITER, BOUT++/Hermes-3) was expanded to evaluate wall and divertor heat fluxes and interface with engineering thermal analysis. A parameterized CAD tool, TRACER, was expanded to generate medium-fidelity divertor, blanket, and coil geometries; OpenFOAM and Diablo workflows were applied for first-wall and divertor thermal analyses with helium cooling; and reduced-order models were created for high-mass-flux divertor cooling. Magnet multiphysics capabilities were verified between Elmer, Diablo, and a new MFEM-based solver, and workflows enable stress, thermal, and neutron-fluence analysis of TF coils with neutronics-driven heating. Nuclear and blanket analysis workflows were demonstrated, including tritium breeding, transport, and CFD-informed thermo-mechanical assessment. Preliminary multi-fidelity uncertainty quantification workflows were applied to boundary modeling codes and shown to achieve variance reductions with fewer high-fidelity boundary simulations. Key findings highlight the challenges of resolving the ITEP gap to find suitable balance between wall and divertor loads, neutron heating, and practical limits of PFC cooling. Next step priorities are to develop automated workflows to check boundary code convergence and detachment, implement tighter physics-engineering CAD provenance tracking, and inclusion of plasma-material interface models for SLAG and tungsten cracking behavior. Collectively, these developments establish sophisticated capabilities for predictive, multi-fidelity, whole-device modeling that integrates plasma physics, materials, magnets, and nuclear engineering to guide pathways to viable Fusion Pilot Plant design points.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Combining physics-based and data-driven models for quantitatively accurate plasma profile prediction that extrapolates well; with application to DIII-D, AUG, and ITER tokamaks

For design, scenario planning, and control, ITER and all other envisioned tokamaks rely on a variety of statistical and physics-based models to extrapolate to unseen regimes; most notably from low plasma current to high. A 'meta-learning' methodology for combining the accuracy of data-driven models with the generalizability of physics-based models is described and tested, yielding a 5–10 percent improvement in performance beyond either alone for the task of extrapolating time-dependent plasma profile prediction from low- to high- plasma current DIII-D tokamak discharges. Meanwhile, it is shown that both machine learning models extrapolated far-distribution and state-of-the-art 'physics-based' profile predictors fare worse than merely assuming plasma profiles do not change from their initial values. Finally, a variety of other mechanisms for helping data-driven models generalize—transfer learning, adding contextual information from physics simulators, and adding data from the ASDEX Upgrade tokamak—are attempted for similar extrapolation tasks but, in the methodology used in this paper, yield no significant improvement beyond simple data-driven models. Results are summarized in figures 15 and 16.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗