Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tokamaks”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Linear gyrokinetic simulations of toroidal Alfvén eigenmodes in the Mega-Amp Spherical Tokamak

Linear gyrokinetic (GK) simulations using the Gyrokinetic Toroidal Code (GTC) have been performed to investigate Toroidicity-driven Alfvén Eigenmodes (TAEs) driven by the neutral beam injection (NBI) induced fast ions in the Mega-Amp Spherical Tokamak (MAST) to identify the non-perturbative and kinetic effects of thermal plasma. A specific TAE in MAST discharge 26887, with an on-axis NBI power of approximately 1.5 MW and plasma current around 800 kA, exhibited frequency chirping, and the tangential soft x-ray camera array resolved the radial mode structure peaked near |q|=1.5. Various excitation methods were used in the GTC linear simulations, illustrating this code's capability to realistically represent the mechanisms and behaviors of fast ion-driven TAEs in spherical tokamaks. The radial structures from these GK simulations closely match measurements and calculations performed using the NOVA ideal MHD code, though with the frequencies approximately 10 kHz lower, likely due to various kinetic and non-perturbative effects. The simulations measured the damping rates due to continuum damping, radiative damping, and ion Landau damping, revealing that ion Landau damping has the most significant contribution to the total damping rate of the TAE. A comparison of growth rates of TAEs excited by fast ion Maxwellian and slowing-down distributions shows that the TAEs excited by a fast ion anisotropic pitch distribution (as part of the slowing-down distributions) are more unstable compared to those excited by a Maxwellian distribution with an equivalent fast ion beta. This shows that the use of fast ion anisotropy alters the number of fast ions to be in shear Alfvén resonance, and hence, it can greatly affect the stability of TAEs. These tests can be performed with the GTC but impossible with ideal MHD simulations, highlighting the necessity of kinetic simulations such as the GTC for a precise prediction of the TAE stability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Vertical instability forecasting and controllability assessment of multi-device tokamak plasmas in DECAF with data-driven optimization

Abstract Reliable vertical position control will be an essential element of any future tokamak-based fusion power plant in order to reduce disruptions and maximize performance. We investigate methods to improve vertical controllability boundary determination in plasma operational space and demonstrate a data-driven approach based on direct pseudoinversion of operational space data that is rigorously quantitative, applicable in real-time plasma control systems, and physically intuitive to interpret. Applied to historical shot data from entire run campaigns on the MAST-U, KSTAR, and NSTX tokamaks, this approach, implemented in DECAF, improves vertical displacement event identification accuracy to 98.9%–100%. Further, we explore the application of a physics-based vertical stability metric as an early warning forecaster for vertical displacement events. The development of a linear surrogate model for the plasma current density profile, with a coefficient of determination of 0.992 on the training dataset, enables potential employment of this forecaster in real-time. The application of this approach on historical data from the MAST-U MU02 campaign yields a forecaster with 62.6% accuracy, indicating promise for this method when further refined and potentially coupled with other stability metrics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A 3D helical filament surrogate model for 3D tokamak equilibria

A novel approach for efficient representation of three-dimensional (3D) tokamak equilibria is investigated, where a set of helical current filaments occupying the plasma region are employed to resolve deviations from the two-dimensional (2D) axi-symmetric state. A discrete set of 3D filaments, located at rational surfaces for a given toroidal mode number n and following the 2D equilibrium field lines (thus forming closed current loops), are found to provide a surrogate model of 3D equilibria with reasonable accuracy. Specifically, application of the filament model to 3D perturbed equilibria, due to the resonant magnetic perturbation (RMP) in DIII-D and MAST-U discharges, reveals that (1) a single helical filament per rational surface is sufficient; (2) 21 such helical filaments are capable of representing the n = 2 3D response field in MAST-U with less than 10% relative error as compared to that computed by a full magnetohydrodynamic code; (3) optimizing currents (both amplitude and phase) flowing in 3D filaments with fixed geometry, the highest accuracy fitting is found to depend on the characteristics of the 3D equilibria such as the coil current phasing of the RMP coils in our case studies. Here, whis filament approach is also applicable for generating surrogate models of other type of 3D tokamak equilibria, including those during the initial phase of the plasma disruption.

MARS-F↗
DOE OSTI2025

Runaway electron-induced plasma facing component damagein tokamaks

This Roadmap article addresses the critical and multifaceted challenge of plasma facing component (PFC) damage caused by runaway electrons (REs) in tokamaks, a phenomenon that poses a significant threat to the viability and longevity of future fusion reactors such as ITER and DEMO. The dramatically increased RE production expected in future high-current tokamaks makes it very difficult to avoid or mitigate REs in such devices when a plasma discharge terminates abnormally. Preventing damage from the intense localised heat loads they can cause requires a holistic approach that considers plasma, REs and PFC damage. Despite decades of progress in understanding the physics of REs and the thermomechanical response of PFCs separately, their complex interplay remains poorly understood. This document aims to initiate a coordinated, interdisciplinary approach to bridge this gap by reviewing experimental evidence, advancing diagnostic capabilities, and improving modelling tools across different scales, dimensionalities, and fidelities. Key topics include RE beam formation and transport, damage mechanisms in both brittle and metallic PFCs, and observed effects in major facilities such as JET, DIII-D, WEST and EAST. The Roadmap emphasises the urgency of predictive, high-fidelity modelling validated against well-diagnosed controlled experiments, particularly in the light of recent changes in ITER’s wall material strategy and the growing importance of private sector fusion initiatives. Each section of the Roadmap article is written to provide a concise overview of one area of this multidisciplinary subject, with an assessment of the status, a look at current and future challenges, and a brief summary. The ultimate goal of this initiative is to guide future mitigation strategies and design resilient components that can withstand the intense localised loads imposed by REs, thus ensuring the safe and sustainable operation of the next generation of fusion power plants.

Ratynskaia, Svetlana [KTH Royal Inst. of Technolog↗
DOE OSTI2026

Impact of impurities on peeling–ballooning modes and turbulence in tokamak plasmas

This study investigates the impact of various impurity species on peeling–ballooning (PB) modes and microturbulence in tokamak plasmas through the extension of traditional two-fluid and gyro-landau-fluid (GLF) models. By incorporating finite Larmor radius (FLR) effects, the analysis provides a comprehensive understanding of impurity-driven impact and its interaction with plasma turbulence. Depending on charge state and local plasma conditions, heavy impurities may exhibit gyro-radii larger than those of main ions, which are captured in the extended GLF model presented. Following the presentation of modified two-fluid equations incorporating impurity effects, we systematically analyze the distinctions between impurity and main ion dynamics and their resultant feedback mechanisms on plasma behavior. Derivation of the linear dispersion relation enables quantification of impurity-mediated modifications to: plasma vorticity, diamagnetic drift and gyroviscous effects, electron Hall physics, and FLR effects. BOUT++ – based linear simulations corroborate this formalism, demonstrating systematic stabilization of PB modes upon impurity seeding. And then operational implications for practical impurity control strategies in tokamak devices are proposed. The results underscore the necessity of impurity management to maintain stability and optimize plasma confinement, with specific focus on how FLR effects contribute to transport dynamics. This work paves the way for enhanced modeling and simulation efforts, supporting the development of strategies to control impurity-induced turbulence and improve overall reactor performance.

BOUT++ simulation↗
DOE OSTI2024

Overview of physics results from the ADITYA-U tokamak and future experiments

The ADITYA upgrade (ADITYA-U), a medium-sized (R 0 = 75 cm, a = 25 cm) conventional tokamak facility in India, has been consistently producing experiments findings by usingcircular and shaped-plasmas. Recognizing the plasma parameters aligning closely with the design parameters of circular limited plasmas, ADITYA-U shifted its focus toward exploring the operational regime for experimentation on saw-tooth and MHD phenomena. Moreover, ADITYA-U has made consistent advancements toward conducting preliminary plasma shaping experiments through the activation of top and bottom divertor coils utilizing hydrogen as well as deuterium fuels. Confinement is improved by a factor of ~1.5 in D 2 plasmas when compared to H 2 plasmas of ADITYA-U. Further, ADITYA-U operations emphasize preventing disruptions and runaway electrons (REs) to ensure safe operations for future fusion devices. Significant suppression of REs has been achieved in ADITYA-U with the application of pulsed localized vertical magnetic field (LVF) perturbation, thereby establishing the technique's independence from the tokamak device. The successful RE mitigation requires a critical threshold of LVF pulse magnitude, which is approximately 1% of the toroidal magnetic field, and a minimum duration of ~5 ms. Apart from this, several novel findings have been achieved in the ADITYA-U experiments, including the modification of sawtooth duration through gas-puff, the emergence of MHD-induced geodesic acoustic mode-like oscillations, the propagation of fast heat pulses induced by MHD activity, the control of RE dynamics through Gas-puffs, the propagation of pinch-driven cold-pulses, the transport and core accumulations of argon impurities, the mass dependency of plasma toroidal rotation and the detection of 'RICE' scaling, as well as the characterization of edge plasma using wall conditioning methods, such as glow discharge cleaning using a combination of Ar-H 2 mixture, localized wall cleaning by electron cyclotron resonant plasma, and the development of machine learning-based disruption predictions, will be discussed in this paper.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

MHD stability analysis against pressure and current-driven modes in the SMall Aspect Ratio Tokamak

Linear magnetohydrodynamic (MHD) simulations for the SMall Aspect Ratio Tokamak (SMART) have been carried out for the first time, for both positive (PT) and negative triangularity (NT) shaped plasmas using the MARS-F code. The MHD stability of projected SMART plasmas against internal kinks, infernal modes and edge peeling-ballooning modes have been analyzed for a wide range of realistic equilibria. A stabilization of internal kinks and infernal modes is observed when increasing the safety factor profile and reducing plasma beta. PT shaped plasmas are more stable against both internal kinks and infernal modes than their counterpart NT shaped plasmas. Toroidal flows have little impact on the MHD stability of the internal kinks, but they have a strong stabilizing effect on infernal modes, which can be further mitigated in NT shaped plasmas. The MHD stability of peeling-ballooning modes is reduced in NT shaped plasmas, as observed in conventional tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Radiofrequency sheath rectification on WEST: application of the sheath-equivalent dielectric layer technique in tokamak geometry *

Radiofrequency sheath rectification is a phenomenon relevant to the operation of Ion Cyclotron Range of Frequencies (ICRFs) actuators in tokamaks. Techniques to model the sheath rectification on 3D ICRF antenna geometries have only recently become available (Shiraiw et al 2023 Nucl. Fusion 63 026024; Beers et al 2021 Phys. Plasmas 28 093503). In this work, we apply the 'sheath-equivalent dielectric layer' technique, used previously only on linear devices (Beers et al 2021 Phys. Plasmas 28 103508), in tokamak geometry, computing rectified sheath potentials on the WEST ICRF antenna. Advancing the state of the art in sheath rectification modeling, we compute the sheath potentials not just on the limiters, but also on the Faraday Screen bars. The calculations show a peak rectified DC potential of 300 V on the limiters and 500 V on the Faraday screen. Assuming a typical sputtering yield curve, the RF sheath rectification increases the sputtering yield from the limiters by a factor of 2.6 w.r.t. the sputtering due to the non-rectified thermal sheath.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Design workflow of a symmetric traveling wave antenna for fast ion production on DD tokamaks

Initial computational plasma physics scoping and a finite element method antenna modeling design workflow for a symmetric center-fed high-field side high harmonic fast wave traveling wave array (TWA) antenna are reported here. The TWA is designed to generate a test population of fast deuterium ions in an existing D–D tokamak by heating neutral beam deuterium ions, accelerating them from 80 keV to several hundred keV. The resulting fast particles are tailored to mimic key reactor energetic particle parameters with regards to exciting Alfven eigenmode instabilities, allowing for a D–D tokamak like DIII-D or ASDEX-U to replicate reactor-relevant conditions experimentally. Initial scenario scoping for high single-pass absorption as well as good preferential fast ion damping relative to electron damping was completed using the ray-tracing/Fokker–Planck codes GENRAY and CQL3D. Python RF network analysis packages were used to create a custom TWA optimization tool to inform a COMSOL flat antenna design, and Petra-M was used to study cold plasma effects. The TWA produced by this workflow has several novel features when compared to previous TWA studies, including symmetric center feeding, and passive end straps for image current cancellation for reduced impurity production. We show here that the antenna design workflow can readily produce TWA antennas optimized for reflection coefficient, image current cancellation, and launched power spectrum shape; and that a population of fast ions can be generated in the correct region of parameter space, warranting future more detailed studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

On the role of ion temperature gradient turbulence in driving ion thermal transport in neutral beam injection-heated L-mode plasmas in a superconducting tokamak

In this paper, we report a detailed experimental study of the role of ion-temperature-gradient driven turbulence in driving ion thermal transport in neutral beam injection (NBI)-heated L-mode plasmas in Experimental Advanced Superconducting Tokamak (EAST) (Wan et al 2000 Nucl. Fusion 40 1057) coupled with linear and nonlinear gyrokinetic simulations. Significant ion-scale turbulence k ⟂ < cm −1 (k ⟂ < 1.5, where k ⟂ where is the perpendicular wavenumber and ρ s is the ion gyroradius calculated using local electron temperature T e ), measured with a microwave reflectometer, is observed in the plasma core at for r/a ≈ 0.265 and 0.5 (where r is half the diameter of the closed flux surface at a given radial position, and a is half the diameter of the last closed flux surface). Local linear stability analysis with experimental equilibrium quantities at these two radial locations using the GS2 gyrokinetic code shows that the most unstable ion-scale micro-instability is the ion temperature gradient (ITG) mode. Since the computed maximum ion-scale linear growth rates are larger than the local Waltz-Miller E × B shearing rate (Waltz and Miller 1999 Phys. Plasmas 6 4265), the E × B shear is unable to suppress ion-scale turbulence, consistent with the experimental observation of ion-scale turbulence. Experimental ion and electron thermal transports, calculated with power balance analysis, are both anomalous with ion thermal transport being the dominant channel. Nonlinear gyrokinetic simulations using the electrostatic global particle-in-cell gyrokinetic tokamak simulation (GTS) code (Weixing Wang et al 2010 Phys. Plasmas 17 072511) show good quantitative agreement between predicted ion thermal transport and that from the power balance analysis. These results demonstrate for the first time the important role of ITG turbulence in driving thermal transport in NBI-heated L-mode plasmas in EAST, showing the effectiveness of GTS code in simulating electrostatic turbulence in EAST and the need of suppressing ITG turbulence to improve energy confinement in EAST.

Sun, P. J. [Chinese Academy of Sciences (CAS), Hef↗
DOE OSTI2025

Stability evaluation and mitigation strategies in advanced tokamaks using 3D MHD spectroscopy

Multi-modal, active 3D MHD spectroscopy is applied in high-performance advanced tokamak scenarios to study their stability time evolution, revealing an intriguing dependence on both and . A tailored applied 3D field provides a 3D plasma response to extract the growth rate of the least stable mode. The estimated growth rate finds a decrease in stability when the minimum in the safety factor (q) passes through 2.0 and reveals inherent risks of crossing an additional rational surface at integer q min , even above the usual q = 1 sawtooth condition. Based on this result, the potential scenario in which q min ~ 2 can be safely crossed during a more stable lower β N phase was investigated, and the improved stability of this scenario is confirmed by the estimated growth rate. This shows that 3D MHD spectroscopy can offer insights into strategies for improving stability by identifying the vulnerable aspects of such scenarios. In addition, the method highlights its potential for instability avoidance by enabling early detection of multiple modes, even before magnetic coils can measure them. The measured growth rate by the 3D MHD spectroscopy shows its reliability by exhibiting a correlation with the programmed rises in plasma beta across various high β N and high q min discharges. In addition, this method is successfully applied during rapidly evolving I p ramp-up phases, a key part of the scenario development. By achieving reasonable growth rate measurements at high-performance scenario developments, this technique contributes to the development of advanced diagnostic tools for tokamak scenario stability, which will help identify an effective pathway to stable, high-performance scenarios.

Yang, S. M. [Princeton Plasma Physics Laboratory (↗
DOE OSTI2025

Conceptual study on using Doppler backscattering to measure magnetic pitch angle in tokamak plasmas

We introduce a new approach to measure the magnetic pitch angle profile in tokamak plasmas with Doppler backscattering (DBS), a technique traditionally used for measuring flows and density fluctuations. The DBS signal is maximised when its probe beam's wavevector is perpendicular to the magnetic field at the cutoff location, independent of the density fluctuations [Hillesheim \emph{et al} 2022 \emph{Nucl. Fusion} \textbf{55} 073024]. Hence, if one could isolate this effect, DBS would then yield information about the magnetic pitch angle. By varying the toroidal launch angle, the DBS beam reaches cutoff with different angles with respect to the magnetic field, but with other properties remaining similar. Hence, the toroidal launch angle which gives maximum backscattered power is thus that which is matched to the pitch angle at the cutoff location, enabling inference of the magnetic pitch angle. We performed systematic scans of the DBS toroidal launch angle for repeated DIII-D tokamak discharges. Experimental DBS data from this scan were analysed and combined with Gaussian beam-tracing simulations using the Scotty code [Hall-Chen \emph{et al} 2022 \emph{Plasma Phys. Control. Fusion} \textbf{64} 095002]. The pitch-angle inferred from DBS is consistent with that from magnetics-only and motional-Stark-effect-constrained (MSE) equilibrium reconstruction in the edge. In the core, the pitch angles from DBS and magnetics-only reconstructions differ by one to two degrees, while simultaneous MSE measurements were not available. The uncertainty in these measurements was under a degree; we show that this uncertainty is primarily due to the error in toroidal steering, the number of toroidally separated measurements, and shot-to-shot repeatability. We find that the error of pitch-angle measurements can be reduced by optimising the poloidal launch angle and initial beam properties. Since DBS has high spatial and temporal resolutions, is non-perturbative, does not require neutral beams, and is likely robust to neutron damage of and debris on the first mirrors, using DBS to measure the pitch angle in future fusion energy systems is especially appealing.

beam tracing↗
DOE OSTI2026

Empirical scaling of the L–H threshold power for metal wall tokamaks using a multi-device database

The empirical scaling for the H-mode power threshold in tokamaks has been revisited using a database with threshold data from machines with a metallic first wall as part of International Tokamak Physics Activity (ITPA) task TC-26. The database contains discharges from ASDEX Upgrade (AUG) (W), JET (Be/W) and Alcator C-Mod (Mo). This was motivated by reports that in like-for-like discharges the power threshold was reduced by approximately 30% after the change from carbon based to metallic first wall materials on AUG (Ryter et al 2013 Nucl. Fusion 53 113003) and JET (Maggi et al 2014 Nucl. Fusion 54 023007). The database contains L–H transition data for all hydrogen isotopes and mixtures, including T and DT from the recent JET campaigns. Compared to the ITPA 2008 scaling (Martin et al 2008 J. Phys.: Conf. Ser. 123 012033), the metal wall scaling has a smaller magnetic field exponent but a larger density exponent. We present an additional parameter to capture the strong dependence of the L–H power threshold (approx. factor 2) on the magnetic configuration in the divertor on JET. The scaling recovers the approximate inverse isotope mass scaling of the threshold power. Alternative scalings involving the plasma current and poloidal magnetic field are explored. Despite the reduction in threshold observed earlier, the scalings based on the metal wall database do not necessarily extrapolate to a lower threshold for ITER compared to the ITPA 2008 scaling, especially at high density. The divertor configuration effect induces the largest uncertainty in the extrapolation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Macroscopic trends of neoclassical tearing stability in high-field H-mode tokamak pilot plants

The neoclassical tearing mode (NTM) stability metric—minimum marginally stable island width $w$$^{*}_{m}$—was compared across 14651 inductive high-field tokamak pilot plant equilibria. Larger devices with reduced elongation and/or increased minor radius demonstrated an order-of-magnitude increase in $w$$^{*}_{m}$, primarily due to a reduction in bootstrap drive. This work is part of an ongoing effort to ensure passive NTM-stability in the ARC tokamak, in which the technology to achieve active tearing-suppression with localised electron cyclotron current drive does not yet exist. The equilibrium scenarios in the database were Monte Carlo generated and normalised to the same >400MW fusion power, minimum pressure scenario at a range of plasma currents, before tearing analysis using the modified Rutherford equation was applied for all resonant poloidal and toroidal m, n modes up to n = 4. Single-helicity toroidal Δ' calculations in resistive DCON set the minimum marginally stable island width, and a simple modal scaling proportional to –m 2 n –1 was identified for high-m Δ' values. The dominant correlates of $w$$^{*}_{m}$ and Δ' across the database were analysed using interpretable machine learning techniques.

NTM seeding↗
DOE OSTI2026

Variability of MHD instabilities in benign termination of high-current runaway electron beams in the JET and DIII-D tokamaks

Benign termination, in which magnetohydrodynamic (MHD) instabilities deconfine runaway electrons (REs) following hydrogenic injections, is a promising strategy for mitigating dangerous RE loads after disruptions. Recent experiments on the Joint European Torus (JET) have explored this scenario at higher pre-disruptive plasma currents than are achievable on other devices, revealing challenges in obtaining benign terminations at I p ≥2.5 MA. This work analyzes the evolution of these high-current RE beams and their terminating MHD events using fast magnetic sensor measurements and EFIT equilibrium reconstructions for approximately 40 JET and 20 DIII-D tokamak discharges. On JET, unsuccessful non-benign terminations occur at low edge safety factor (q edge ≈ 2), and are preceded by intermittent, non-terminating MHD events at higher rational qedge. Trends in the internal inductance I i indicate more peaked RE current profiles in the high-I p non-benign population, which may hinder successful recombination through re-ionization of the companion plasma. In contrast, benign terminations on JET typically occur at higher q edge ≥3 and exhibit less peaked RE current profiles. DIII-D displays a broader range of terminating edge safety factors, again correlated with the measured values. Across both tokamaks, the RE current peaking is therefore found to determine which MHD instability boundary is encountered, a result confirmed by linear resistive MHD modeling with the CASTOR3D code. Measured growth rates are similar for benign and non-benign cases, indicating that ideal MHD timescales at low density after hydrogenic injection do not alone explain efficient RE deconfinement. Instead, non-benign cases are most readily characterized by their comparably lower overall MHD perturbation amplitudes δB. These observations suggest that the interplay between ideal and resistive dynamics governs the termination process, with implications for extrapolating benign RE termination to high-I p reactor scenarios.

MHD instabilities↗
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 OSTI2026

Results from the last DD and DT JET campaigns in the framework of the EUROfusion Tokamak Exploitation Work Package activity

JET, the only tokamak capable of operating with deuterium–tritium (D–T) fuel (since TFTR was shutdown in 1999), has provided essential experimental data to support ITER and DEMO design and operation. Within the EUROfusion Tokamak Exploitation Work Package, JET completed its final campaigns (2022–2023), culminating in the third D–T campaign (DTE3). These experiments addressed key challenges in plasma scenarios, exhaust control, and tritium management under reactor-relevant conditions. Significant progress was achieved in demonstrating ITER-like integrated scenarios with impurity seeding, achieving partial divertor detachment and high confinement ($H_{98}(y,2)$ ≈ 0.85) at 3 MA in D–T plasmas. Advanced exhaust regimes such as quasi-continuous exhaust (QCE) and X-point radiator (XPR) were successfully achieved first in D–D and then extended to D–T operation, confirming their relevance for mixed isotope operation. Operational milestones included a new world record of 69 MJ fusion energy in tritium-rich hybrid plasmas and long-pulse H-mode operation up to 60 s, contributing with unique data to the CICLOP database. Physics studies focused on peeling-limited pedestals in support of ITER and improved understanding of edge stability and impurity screening in metallic environments. Extensive usage of the shattered pellet injector (SPI) on JET provided critical information for the design of the ITER disruption mitigation system (DMS). Real-time control systems for D/T ratio control and plasma exhaust were deployed and demonstrated in D–D and D–T, while energetic particle physics investigations unfolded the role of fast ions in turbulence suppression mechanisms. Comprehensive tritium retention studies using gas balance method, post-mortem analysis, and ITER-relevant laser induced desorption spectroscopy (LIDS) diagnostics provided essential input for tritium accountancy strategies. These results are validating the ITER operational concepts, inform DEMO design, and deliver critical experience in nuclear operation and scenario integration.

disruptions↗
DOE OSTI2026

The impact of kinetic and global effects on ideal ballooning 2nd stable pedestals of conventional and low aspect-ratio tokamaks

The EPED model (Snyder et al 2011 Nucl. Fusion 51 103016) has had success in describing the pedestal structure of Type-I ELM and QH-mode plasmas in tokamaks, by combining kinetic ballooning mode (KBM) and peeling–ballooning (PB) constraints. Within EPED, the KBM constraint is often approximated by a technique using calculated ideal ballooning mode (IBM) thresholds, fit to a functional form designed to capture non-local effects. It has been noted that quantitative differences between local ideal MHD and gyro-kinetic (GK) ballooning stability can be larger at low aspect ratio. KBM critical pedestals are consistent with observations in initial studies on conventional and spherical tokamaks. In this work, the application of a reduced model for the calculation of the kinetic ballooning stability boundary is presented based on a novel and newly developed gyro-fluid system (GFS) code (Staebler et al 2023 Phys. Plasmas 30 102501). GFS is observed to capture KBMs in DIII-D as well as the NSTX pedestals, opening a route to integrating this model into EPED. Finally, high but finite $n$ global ballooning modes are observed to limit the access to the local 2nd stability and thus provide a transport mechanism that constrains the pedestal evolution with $β_{p,\textrm{ped}}$. The high $n$ global ballooning stability is approximated by its ideal MHD analog using ELITE. It is shown that nearly-local high $n$ modes with $k_yρ_s$ ~ $0.25 - 0.5$ can provide a proxy for the critical $β_{p,\textrm{ped}}$ when a 2nd stable access exists on DIII-D plasmas. The use of GFS and ELITE scaling in EPED provides improved agreement to EPED1 in an initial comparison with DIII-D pedestal data.

Anastopoulos Tzanis, Michail Savvas [Oak Ridge Nat↗