Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tokamak integrated modeling”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Flexible, integrated modeling of tokamak stability, transport, equilibrium, and pedestal physics

The STEP (Stability, Transport, Equilibrium, and Pedestal) integrated-modeling tool has been developed in OMFIT to predict stable, tokamak equilibria self-consistently with core-transport and pedestal calculations. STEP couples theory-based codes to integrate a variety of physics, including magnetohydrodynamic stability, transport, equilibrium, pedestal formation, and current-drive, heating, and fueling. The input/output of each code is interfaced with a centralized ITER-Integrated Modelling & Analysis Suite data structure, allowing codes to be run in any order and enabling open-loop, feedback, and optimization workflows. This paradigm simplifies the integration of new codes, making STEP highly extensible. STEP has been verified against a published benchmark of six different integrated models. Core-pedestal calculations with STEP have been successfully validated against individual DIII-D H-mode discharges and across more than 500 discharges of the H98,y2 database, with a mean error in confinement time from experiment less than 19%. STEP has also reproduced results in less conventional DIII-D scenarios, including negative-central-shear and negative-triangularity plasmas. Predictive STEP modeling has been used to assess performance in several tokamak reactors. Simulations of a high-field, large-aspect-ratio reactor show significantly lower fusion power than predicted by a zero-dimensional study, demonstrating the limitations of scaling-law extrapolations. STEP predictions have found promising scenarios for an EXhaust and Confinement Integration Tokamak Experiment, including a high-pressure, 80%-bootstrap-fraction plasma. ITER modeling with STEP has shown that pellet fueling enhances fusion gain in both the baseline and advanced-inductive scenarios. Finally, STEP predictions for the SPARC baseline scenario are in good agreement with published results from the physics basis.

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

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 OSTI2026

Towards self-consistent integrated modeling of the tokamak pedestal, scrape-off layer, and divertor using SOLPS-ITER and EPED

A predictive, self-consistent integrated model of the tokamak edge plasma and neutral system has been developed and validated against experimental data from DIII-D, spanning from the divertor target and scrape-off layer to the top of the pedestal. The model uses the SOLPS-ITER fluid boundary plasma and neutral code to constrain the EPED prediction of the pedestal height and width by taking the separatrix temperature and density calculated from SOLPS as inputs to EPED. An empirical proxy function based on data from a subset of these discharges is then used to infer the ratio of density at the pedestal top to that at the separatrix based on the electron temperature at the divertor target. This allows the density at the top of the pedestal to be provided to EPED as well, so that knowledge of this quantity is no longer required a priori, while SOLPS radial transport coefficients are needed instead. The gas puff rate provided to SOLPS is then a primary model input for this forward model coupling (along with radial transport coefficients, which are kept fixed through each scan in this study). Comparisons with measurements from DIII-D density scans through a range of target conditions find significant pedestal pressure degradation as the density is increased towards detachment in three different divertor configurations, consistent with the experimental data in these ballooning-limited pedestal regimes.

core-edge integration↗
DOE OSTI2025

Integrated modeling of advanced tokamaks for tearing mode avoidance on DIII-D

Understanding the physics of the onset and evolution of tearing modes (TMs) in tokamak plasmas is crucial for high-performance steady-state operations. The onset condition of (m,n) = (2,1) tearing stability, $Δ′$ > $Δ_c^′$> 0, has been studied in DIII-D steady-state hybrid discharges with accurate equilibrium reconstruction and well-measured plasma profiles. Here, m is the poloidal mode number and n is the toroidal mode number. The onset of n = 1 TMs is carefully identified using the Mirnov array to determine in-phase signals for the onset time and using electron cyclotron emission measurements to trace back from the saturated island width in the experiments for both initial island width and its growth rate. The tearing stability index $Δ′$ at TM onset is calculated in multiple ways, using the modified Rutherford equation and using the PEST-III or resistive DCON magnetohydrodynamic codes, and the results are in reasonable agreement with each other. The calculated $Δ′$ is compared with an analytical formula of the tearing stability threshold $Δ_c^′$ to determine the condition for the mode onset, $Δ′$ > $Δ_c^′$. Finally, the local gradient effect of the plasma current and the pressure profiles on the tearing stability is investigated in terms of the difference between $Δ′$ and $Δ_c^′$ through a series of modeled equilibria with the parametric variations. A database of experimental mode instability onsets is compatible with the simulated regions of instability obtained from the parametric variations.

Linear stability analysis↗
DOE OSTI2020

Predictive Integrated Modeling of Heated Tokamak Plasmas (Final Technical Report)

The research carried out by the Lehigh University fusion group during this reporting period: (1) Self-consistent predictive core-pedestal ITER scenario modeling; (2) Microtearing instabilities and electron thermal transport in standard and spherical tokamak discharges; (3) Neural network model of the Multi-Mode anomalous transport module; (4) Anomalous transport in KSTAR, JET, EAST, Alcator C-Mod and DIII-D discharges; (5) Theory for transport in magnetized plasmas; (6) Control-oriented modeling of the q-profile dynamic response in preparation for advanced scenario control in KSTAR and in EAST; (7) Modeling of edge harmonic oscillation in DIII-D QH mode discharges; (8) Elucidating plasma dynamics in Hasegawa-Wakatani turbulence by information geometry; (9) Non-linear Langevin and Fractional Fokker–Planck equations for anomalous diffusion by Levy stable processes; and (10) Investigation of the plasma shaping effects on the H-mode pedestal structure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Results and lessons learned from accelerating radio frequency modeling using machine learning [slides]

The “advanced tokamak” reactor concept is a leading candidate for a steady state fusion pilot plant. An advanced tokamak (AT) sustains a majority of the required plasma current with effects resulting from maintenance of the peaked pressure at the device center. This current is augmented by auxiliary current drive sources. These auxiliary actuators may consist of neutral particle beams and/or radio frequency (RF) systems such as lower hybrid current drive (LHCD) and high harmonic fast wave (HHFW) current drive using radio and microwaves from antennas. The primary focus of this work is to develop models of RF current profile control suitable for use in integrated modeling frameworks and for real-time control in experiments. Direct physics models of RF current drive can be computationally intensive. In order to achieve predictive times appropriate for the thousands of calls needed in real-time control of experiments and for use in integrated models, we will apply modern machine learning (ML) techniques to accelerate these models and interpolate their results. To generate the fast and accurate models for use in control level algorithms and integrated modeling we need to replace present models with high dimensional interpolation of their results. We will perform additional simulations across a broader parameter range for EAST and other tokamaks in different physics regimes (Alcator C-Mod, DIII-D, WEST, CFETR, ARC, ITER) and combine them into a larger database for training and testing of the ML models. Further testing of the control level models with experimental current profile data from EAST and C-Mod tokamaks will provide additional confirmation of the control level model before integration in a tokamak control system or integrated modeling suite. ML will be used to optimize the selection of training data consisting of RF current driven at different values of density profile, temperature profile, plasma current, and wavenumber. ML will also be used to facilitate classification of current drive from these input data. The output of this effort will be a validated classifier capable of determining the current drive profiles for HHFW CD and LHCD on a mille-second timescale. This will provide a breakthrough capability enabling real-time control of RF driven current profiles in experiments including ITER ICRF and use integrated modeling frameworks requiring thousands of current profile calculations in discharge simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Doppler Backscattering Data Analysis and Integrated Modeling with OMFIT

One Modeling Framework for Integrated Tasks (OMFIT) is a widely used software tool in the magnetic fusion research community. OMFIT provides magnetic fusion energy researchers with a framework for the development of special-purpose physics modules. This paper describes an OMFIT physics module pertaining to the Doppler Backscattering (DBS) fusion plasma diagnostic. DBS measures density fluctuations and flow velocity through plasma scattering of electromagnetic waves. The OMFIT DBS module was developed to analyze experimental DBS data and facilitate modeling of DBS systems installed on multiple tokamak devices. The OMFIT DBS module is designed to support several analysis workflows: detailed analysis of experimental data, experimental planning, and theory-based synthetic diagnostic modeling. The DBS module uses integrated modeling by leveraging other OMFIT physics modules to perform tasks related to DBS, e.g. ray/beam–tracing simulations, edge-localized mode–synchronized data analysis, magnetic equilibrium reconstruction, and fitting kinetic profile data. Furthermore, this paper describes several supported workflows and serves a reference for the OMFIT DBS module.

Doppler backscattering↗
DOE OSTI2020

Neural-network accelerated coupled core-pedestal simulations with self-consistent transport of impurities and compatible with ITER IMAS

pedestal structure, current profile, and plasma equilibrium physics has been developed and tested against a DIII-D discharge. Here, key features of the achieved core-pedestal coupled workflow are its ability to account for the transport of impurities in the plasma self-consistently, as well as its use of machine learning accelerated models for the pedestal structure and for the turbulent transport physics. Notably, the coupled workflow is implemented within the OMFIT framework, and makes use of the ITER integrated modeling and analysis suite (IMAS) data structure for exchanging data among the physics codes that are involved in the simulations. Such technical advance has been facilitated by the development of a new numerical library named OMAS.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

SICAS, a new integrated divertor-SOL-pedestal-core framework for self-consistent modeling of ion and impurity transport in fusion devices

We present a new framework for core-edge integration studies named SICAS (SOLPS-ITER coupled to ASTRA-STRAHL) which enables high fidelity simulations of the core, edge, and divertor regions encompassing the transport of ions as well as the impurities through the entire plasma domain. SICAS handles the exchanging of the particle and power fluxes as well as transport coefficients to ensure consistency through the codes. An overlapping region is defined from the inside of the separatrix of the plasma to the SOLPS-ITER core boundary allowing for matching profiles and fluxes between the two codes for a self-consistent approach. The results presented here demonstrated the flexibility of SICAS to simulate different configurations, scenarios, divertor geometries, and plasma species with good agreement with DIII-D experimental data. This tool opens new possibilities in integrated modeling of fusion devices integrating all relevant phenomena in the core and the divertor plasmas. These capabilities are required for the interpretation of current experiments as well as the design of new devices.

core-edge integration↗
DOE OSTI2026

Validation of an integrated modeling framework for investigating 3D plasma responses in tokamak plasmas

As contemporary experimental tokamaks are pushed toward reactor-relevant operation, they provide essential testbeds for demonstrating and exploring ELM control strategies for deployment in future fusion pilot plants. Accurate predictions of full plasma responses are essential to guide and optimize these demonstrations. This paper introduces and validates an integrated modeling framework over a historical range of DIII-D operational space. The integrated modeling framework uses only scalar plasma parameters and optional reference boundary to self-consistently and flexibly scan through tokamak operational space and estimate the corresponding ELM-suppression relevant plasma response over 3D coil phase space. The framework generates tightly converged equilibria that satisfy a target set of plasma parameters (I p , β N , l i ), with kinetic profiles constrained by an EPED(NN)-computed pedestal and empirical core model. The plasma response of these modeled equilibria is calculated with GPEC. Validated against 55 distinct DIII-D equilibria, the framework consistently reproduces experimentally constrained equilibrium pressure, q, and other representative profiles using the time-varying information from only evolving scalar plasma parameters; plasma response validation is performed for the n=3 perturbation. In particular, the plasma response of a DIII-D discharge scanning q 95 in search of ELM suppression windows is quantitatively reproduced. By accurately reproducing experimental equilibria and plasma responses across wide parameter variation, the framework supports potential for synthetic parameter scans in key operational, stability, and plasma response dimensions to investigate RMP ELM-suppression experiments and inform predictive RMP scenario optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Self-consistent modeling of tokamak edge plasma transport with lithium sources

Magnetic confinement fusion devices require effective heat and particle exhaust solutions on the divertor plates to operate sustainably, especially under reactor-relevant conditions. Liquid lithium divertors have been proposed to address two major challenges: control of excessive heat flux to plasma-facing components through vapor shielding and minimization of core plasma contamination from impurities. The National Spherical Torus Experiment-Upgrade (NSTX-U) will explore lithium as a divertor material due to its potential to meet both objectives. We present a self-consistent coupling framework between the plasma boundary transport code UEDGE and the lithium wall transport code Wall–Li to evaluate the feasibility and operational limits of lithium-based divertors. The model aims to optimize lithium sourcing levels to prevent core plasma contamination via fuel dilution while ensuring divertor protection through vapor shielding. This integrated framework, applicable to any tokamak with lithium sources, dynamically adjusts lithium sourcing based on local plasma conditions and surface temperature. The coupled model is tested using NSTX-like geometry and plasma conditions to assess its performance and reliability. Wall–Li calculates lithium fluxes from plasma-facing components, incorporating physical sputtering, thermally enhanced sputtering, and evaporation driven by surface temperature and ion flux. These fluxes are reintroduced into UEDGE as neutral lithium atoms, enabling simulation of their transport and distribution within the plasma. UEDGE computes plasma and neutral transport, surface heat flux, and iteratively feeds this information back to Wall–Li. A small time step is employed to ensure numerical stability and convergence, enabling accurate simulations over typical tokamak discharge durations. This integrated modeling approach provides a robust tool for identifying operational regimes that balance effective lithium sourcing with minimal core plasma contamination, offering critical insights for optimizing lithium-based divertor systems in current and future fusion devices.

Magnetic confinement fusion↗
DOE OSTI2025

TRANSP integrated modeling code for interpretive and predictive analysis of tokamak plasmas

This paper provides a comprehensive review of the TRANSP code, a sophisticated tool for interpretive and predictive analysis of tokamak plasmas, detailing its major capabilities and features. It describes the equations for particle, power, and momentum balance analysis, as well as the poloidal field diffusion equations. The paper outlines the spatial and time grids used in TRANSP and details the equilibrium assumptions and solvers. Various models for heating and current drive and radiation, including updates to the NUBEAM model, are discussed. The handling of large-scale events such as sawtooth crashes and pellet injections is examined, along with the predictive capabilities for advancing plasma profiles. The integration of TRANSP with the ITER Integrated Modeling and Analysis Suite (IMAS) is highlighted, demonstrating enhanced data access and analysis capabilities. Additionally, the paper discusses best practices and continuous integration techniques to enhance TRANSP's robustness. The suite of TRANSP tools, designed for efficient data analysis and simulation, further supports the optimization of tokamak operations and coupling with other tokamak codes. Continuous development and support ensure that TRANSP remains a major code for the analysis of experimental data for controlled thermonuclear fusion, maintaining its critical role in supporting the optimization of tokamak operations and advancing fusion research.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Integrated operation scenarios: Chapter 6 of the special issue: on the path to tokamak burning plasma operation

Here we report the progress of the development and optimization of operational scenarios for ITER and beyond, focusing upon baseline, hybrid, and steady-state scenarios since 2007. This includes advancements made by the integrated operation scenarios (IOS) topical group of the international tokamak physical activity as well as contributions from the broader tokamak community. The key area of research involves developing IOSs that encompass tokamak physics, operation, and technology by utilizing integrated modeling and control strategies. This requires leveraging available actuators to simultaneously control plasma position and shape, MHD activities that could lead to disruptions, transport, plasma-wall interaction and power exhaust, fuel cycle, fusion burn, and tritium breeding. The control extends from the plasma initiation phase, through the current ramp-up, flattop, start and end of the fusion burn, and current ramp-down, to the plasma termination phase. A review of the currently developed scenarios and modeling is provided in terms of (i) optimizing plasma initiation in ITER, (ii) preparing for the low activation phase to fully commission all tokamak systems and establish and validate physics and scenario conditions in preparation for deuterim-tritium (DT) operation, (iii) developing and preparing baseline and hybrid scenarios to demonstrate the feasibility of achieving these regimes within device constraints, (iv) exploring steady-state scenarios to meet ITER’s steady-state goals, (v) evaluating and preparing actuators for ITER, (vi) developing integrated control solutions using shared actuators. The most notable achievements include; (i) the development of ITER demonstration discharges by matching various dimensionless parameters, (ii) the development of scenarios in an ITER-like tungsten environment and DT operation, and (iii) the development of scenarios in superconducting tokamaks, enabling long-pulse operations with similar coil constraints to ITER. Along with these significant achievements, outstanding issues and recommendations for further research and development are provided. Importantly, this study goes beyond simply updating the ITER Physics Basis; it carries profound implications for the broader field of burning plasma research, offering valuable insights and guidance for the next generation of fusion experiments and devices.

ITER↗
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 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↗