Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tokamak edge”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Effect of edge plasma density on hot spot in LHCD plasma in EAST

Hot spots are a serious challenge limiting long pulse operation with lower hybrid current drive (LHCD) in tokamaks. Here, in order to mitigate hot spots in the guard limiter and improve LHCD capability, the effect of edge plasma density on inducing hot spots and current drive has been studied in EAST. The temperature in the guard limiter of the LH antenna, inducing hot spot directly, increases with edge density and LH power. Studies show that the hot spot is mainly ascribed to the heat flux in front of LH antenna. Further simulation indicates that such spots correspond to the peak position of edge density due to local LH electric field. In addition, due to the stronger parametric instability (PI) behavior in the case of higher edge density, the current drive capability decreases with edge density. Strike-point splitting behaviour appears as density increase, in agreement with current profiles in the edge region and the reduction of total driven current, suggesting that more power is deposited in the edge region, which then contributes more to hot spot formation. These studies offer one possible idea to optimize the edge density so as to satisfy the coupling, mitigate the heat flux in the guard limiter, and improve current drive capability in fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Nonlinear second order electromagnetic gyrokinetic theory for a tokamak plasma

The steep plasma pressure gradient that forms at the edge of the high confinement, H-mode regime of tokamak operation provides free energy to drive electromagnetic micro-instabilities that are widely believed to influence the transport processes in this so-called pedestal region. This high pressure gradient also provides a high current density (bootstrap current), known to influence ballooning mode stability and to be important for driving kink modes in the ideal magneto-hydrodynamic plasma model (so-called peeling-ballooning modes). Furthermore, efficient, steady state future tokamak power plants must operate with a large bootstrap current in the core and especially concerning spherical tokamaks, confinement will be influenced by electromagnetic turbulence. To accommodate these important situations, conventional electromagnetic gyrokinetic theory is extended to incorporate neoclassical effects in the equilibrium drives, allowing $B_{\vartheta} ~ B_0$ (B0 is the confining magnetic field, and $B_{\vartheta}$ is its poloidal component). This provides a global gyrokinetic model that self-consistently captures the consequences of large bootstrap current fractions on the equilibrium distribution functions.

kink modes↗
DOE OSTI2021

Confinement in electron heated plasmas in Wendelstein 7-X and ASDEX Upgrade; the necessity to control turbulent transport

In electron (cyclotron) heated plasmas, in both ASDEX Upgrade (L-mode) and Wendelstein 7-X, clamping of the ion temperature occurs at $T_\mathrm{i}$ ~ 1.5 keV independent of magnetic configuration. The ions in such plasmas are heated through the energy exchange power as ${n}_{\mathrm{e}}^{2}({T}_{\mathrm{e}}-{T}_{\mathrm{i}})/{T}_{\mathrm{e}}^{3/2}$, which offers a broad ion heating profile, similar to that offered by alpha heating in future thermonuclear fusion reactors. However, the predominant electron heating may put an additional constraint on the ion heat transport, as the ratio $T_\mathrm{e}T_\mathrm{i}$ > 1 can exacerbates ITG/TEM core turbulence. Therefore, in practical terms the strongly 'stiff' core transport translates into Ti-clamping in electron heated plasmas. Due to this clamping, electron heated L-mode scenarios, with standard gas fueling, in either tokamaks or stellarators may struggle to reach high normalized ion temperature gradients required in a compact fusion reactor. The comparison shows that core heat transport in neoclassically optimized stellarators is driven by the same mechanisms as in tokamaks. The absence of a strong H-mode temperature edge pedestal in stellarators, sofar (which, like in tokamaks, could lift the clamped temperature-gradients in the core), puts a strong requirement on reliable and sustainable core turbulence suppression techniques in stellarators.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Further modeling of q 95 windows for the suppression of edge localized modes by resonant magnetic perturbations in the DIII-D tokamak

An improved resonant plasma response model that more accurately captures the physics of the interaction between a tokamak plasma and a resonant magnetic perturbation (RMP) is developed. The model interpolates between the linear and the nonlinear response regimes and takes into account the fact that the slip-frequency is non-zero in the nonlinear regime. The improved model is incorporated into the extended perturbed equilibrium code (EPEC) toroidal asymptotic matching code. The modified EPEC code is used to investigate RMP-induced edge-localized-mode (ELM) suppression in DIII-D H-mode discharge #145380. Somewhat surprisingly, allowing for a finite slip-frequency (i.e., relaxing the so-called no-slip constraint) is found to only slightly facilitate the locking of driven magnetic island chains to the RMP, and, hence, to only slightly facilitate RMP-induced ELM suppression. This is true despite the fact that the nature of non-locked island solutions is radically different when the no-slip constraint is imposed compared to when it is relaxed (in the first case, the widths of the island chains driven at the rational surfaces pulsate, and in the second case, they remain steady). The previously obtained conclusion that the response of a typical H-mode tokamak plasma to an RMP cannot be accurately modeled by linear theory is confirmed. Here, the previously obtained conclusion that the best agreement between theory and observations is achieved by assuming that the natural frequencies of tearing modes, in the absence of the RMP, are determined by the local equilibrium E × B velocity is also confirmed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Suppression of edge localized modes with real-time boron injection using the tungsten divertor in EAST

We report an observation of robust suppression of edge-localized modes (ELMs) in the Experimental Advanced Superconducting Tokamak (EAST), enabled by continuous boron (B) powder injection. Edge harmonic oscillations appear during B powder injection, providing sufficient particle transport to maintain constant density and avoid impurity accumulation in ELM-stable plasmas. Here, quasi-steady ELM suppression discharges are demonstrated with modest energy confinement improvement and over a wide range of conditions: heating power and technique variation, electron density range over a factor ~3.5, deuterium or helium ion species, and with either direction of the toroidal magnetic field. ELM suppression is observed above a threshold edge B intensity and ceases within 0.5 s of termination of the B injection. In contrast to ELM suppression accompanied by recycling reduction during Li powder injection in NSTX and EAST (Maingi et al 2018 Nucl. Fusion 58 024003), reduced recycling due to hydrogenic species retention is unnecessary for the ELM suppression with B powder injection, paving the way for its consideration as an ELM control tool for future fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Real-time observation of toroidal current redistributions induced by three-dimensional MHD phenomena triggering vertical displacement events in tokamak plasmas

Three-dimensional MHD instabilities, including edge-localized modes (ELMs) and internal reconnection events (IREs), have been observed to precipitate loss of vertical stability in tokamak plasmas, resulting in vertical displacement events (VDEs). This vertical destabilization can occur due to toroidal current redistributions and/or shape changes resulting from these phenomena. Using a recently introduced method for rapidly reconstructing the two-dimensional toroidal plasma current density profile in real-time, results are presented that demonstrate the specific current distribution changes that occur during ELMs (on KSTAR) and IREs (on MAST-U) that lead to loss of vertical control. The method most efficiently reconstructs the toroidal current density profile by doing so on a basis of principal components of historical profiles. These principal components isolate dominant current profile dynamics, improving interpretability, increasing speed, and reducing dimensionality of the profile computation. On KSTAR, this computation is executed in the real-time plasma control system at a rate of 10 kHz (limited by available CPU cycle times), allowing the current profile evolution to be assessed at several times over the course of each ELM event. Further, by incorporating the reconstructions into a novel vertical stability metric, the contribution of specific current profile dynamics to the loss of vertical stability can be assessed in real-time for VDE avoidance and improved understanding of the causal relationship between three-dimensional MHD phenomena and VDEs. The success of this method in approximating toroidal current density profiles from kinetic equilibrium reconstructions is also presented ($R^2=0.990$), along with its capability to produce other equilibrium quantities of interest in real-time at high time resolution.

edge-localized modes↗
DOE OSTI2022

Neoclassical transport due to resonant magnetic perturbations in DIII-D

In this work, the role of neoclassical physics in the particle and energy transport during the application of resonant magnetic perturbations (RMPs) to suppress the edge localised modes in a tokamak is analysed. The neoclassical fluxes in non-axisymmetric DIII-D equilibria with applied RMPs are calculated using the NEO code. The magnetic field provided to NEO as an input is calculated using M3D-C1 and includes the nonlinear one-fluid plasma response. Neoclassical fluxes obtained in this study are found to dramatically increase in the presence of applied resonant magnetic perturbations, and are in same range as the total radial particle fluxes calculated in comparable RMP discharges in DIII-D. This suggests that neoclassical transport plays a significant role in edge transport when RMPs are present. An increase in neoclassical fluxes during the ELM suppressed phase in DIII-D plasmas is calculated and is strongly correlated with the observation of density pump-out in the experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Using convolutional neural networks to detect edge localized modes in DIII-D from Doppler backscattering measurements

In H-mode tokamak plasmas, the plasma is sometimes ejected beyond the edge transport barrier. These events are known as edge localized modes (ELMs). ELMs cause a loss of energy and damage the vessel walls. Understanding the physics of ELMs, and by extension, how to detect and mitigate them, is an important challenge. In this paper, we focus on two diagnostic methods—deuterium-alpha (D α ) spectroscopy and Doppler backscattering (DBS). The former detects ELMs by measuring Balmer alpha emission, while the latter uses microwave radiation to probe the plasma. DBS has the advantages of having a higher temporal resolution and robustness to damage. These advantages of DBS diagnostic may be beneficial for future operational tokamaks, and thus, data processing techniques for DBS should be developed in preparation. In sight of this, we explore the training of neural networks to detect ELMs from DBS data, using D α data as the ground truth. With shots found in the DIII-D database, the model is trained to classify each time step based on the occurrence of an ELM event. The results are promising. When tested on shots similar to those used for training, the model is capable of consistently achieving a high f1-score of 0.93. Furthermore, this score is a performance metric for imbalanced datasets that ranges between 0 and 1. We evaluate the performance of our neural network on a variety of ELMs in different high confinement regimes (grassy ELM, RMP mitigated, and wide-pedestal), finding broad applicability. Beyond ELMs, our work demonstrates the wider feasibility of applying neural networks to data from DBS diagnostic.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A new discovery of edge localized modes suppression using ICRH

Here, the high-confinement mode (H-mode) is very important for the fusion performance of ITER and future fusion reactors. The H-mode pedestal is a region featuring strongly reduced turbulence and transport just inside the limiting flux surface with strong plasma gradients, which drive edge localized modes (ELMs) in tokamaks. The ELMs, however, are quasiperiodic relaxations of the pedestal, resulting in a series of hot plasma eruptions similar to magnetospheric substorms, solar and stellar flares, which could potentially damage the ITER divertor plates and first walls. To mitigate and suppress the ELMs in H-mode plasmas, the fusion community has tried to develop effective actuators over the past two decades. Two separate excellent inventions have been awarded to recognize their achievements. The 2018 APS excellence award in Plasma Physics Research was given to Drs. Todd E. Evans, Max E. Fenstermacher, and Richard Alan Moyer for their first experimental demonstration of the stabilization of ELMs in high-confinement diverted discharges by application of very small edge-resonant magnetic perturbations, which led to the adoption of suppression coils in the ITER design. The 2018 Fusion Technology Award from the IEEE’s Nuclear and Plasma Science Society (NPSS) was given to Dr. Larry Baylor for his work designing the fueling, pumping and disruption mitigation system for the US ITER Project and for the implementation of pellet ELM mitigation for ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Edge turbulence velocity preceding the L-H transition in NSTX

The low-to-high mode or L-H transition in tokamaks involves a sudden reduction in the edge turbulence level and decrease in the edge plasma transport. The mechanism for the L-H transition is widely believed to be associated changes in the the poloidal velocity of the edge turbulence. Using a gas puff imaging (GPI) diagnostic, the poloidally averaged poloidal turbulence velocity < V pol > pol was measured for a set of 16 shots with L-H transitions in NSTX. These measurements were focused -2 cm inside the separatrix where the relative fluctuation level as seen in GPI decreases significantly and consistently at the transition. Herein, the magnitude of this poloidal velocity preceding the transition varied widely, ranging from |V pol | ~1-4 km/sec during the last 1 msec before the transition, and the magnitude of the radial gradient also had a wide shot-to-shot range of |grad V pol | ~ 0.34 to 1.9 km/sec/cm during this same time. The frequency spectrum of < V pol > pol during the 30 msec before the transition had a near-coherent peak at ~2-5 kHz, but only in about half of these shots. In general, there were no clear and consistent changes in the poloidal velocity of the turbulence preceding the L-H transition in this database.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A model investigation of the impact of lower hybrid wave scattering angle on current drive profile in EAST and Alcator C-Mod

Lower hybrid current drive (LHCD) is beneficial for developing a steady-state operation scenario in a tokamak. This paper conducts a modeling investigation to identify an optimum rotation angle of the initial lower hybrid perpendicular (to the background magnetic field $\overrightarrow{B}$) wavevector for best matching the experimental RF current profile. It is hypothesized that central RF power deposition widely observed in the present-day LHCD experiments arises from wave scattering by turbulence. In a standard model without considering such interactions, the predicted power deposition profile is generally broad with off-axis peaking, not in agreement with experimental observations. A heuristic approach is adopted by introducing a spectral broadening mechanism by modifying the initial orientation of the perpendicular wavevector. The ray-tracing/Fokker–Planck solver GENRAY/CQL3D is utilized within the python-based π-scope framework. A focus is given to identify the perpendicular wavenumber orientation angle with respect to the magnetic surface normal vector at the initial ray location. Our modeling study shows that rotating the perpendicular wavevector in such a way as to increase the initial poloidal component is effective in reproducing the centrally peaked current profile observed in normal shear plasmas on both EAST and C-Mod. These waves can readily be absorbed to the central plasma, which reduces the sensitivity of the power deposition profile to a slight change of the plasma condition. The same approach is also found to help broaden the off-axis power deposition profile in a reverse-shear EAST plasma, leading to a better agreement with the experiment. The results presented here suggest that spectral modification arising from edge density fluctuations in a tokamak may need to be considered in understanding wave propagation and absorption. A further experimental and theoretical/modeling study is vital as a reverse approach is adopted in this study. Finally, our work suggests that mitigation or control measures are critical for parasitic effects occurring on the first pass in a reactor regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Power and particle exhaust in the ST-E1 fusion power plant

Power exhaust challenges and potential solutions for a 5 m major radius, low-aspect ratio burning tokamak have been explored. 1D edge plasma models have been used to screen for access to detachment using short and long outer divertor legs in double and single null configurations, using Ar as the primary impurity and assuming tungsten plasma-facing components (PFCs). These show that detachment access can be accessed for all but the most conservative assumptions on scrape-off layer (SOL) width and power, but that trade-offs will be required between magnet engineering and the size of the acceptable window of as-yet uncertain plasma parameters. SOLPS-ITER was used to further model selected plasma scenarios, confirming that Ar seeding can be used to achieve dissipative divertor scenarios with peak deposited heat fluxes below 15 MWm -2 . Initial scoping of first wall loads and positioning of limiters has been carried out, showing the feasibility of protecting the breeding blanket wall during steady state without impeding tritium breeding. Initial PFC technology selection is also presented, identifying this as a critical area where further work is needed to find an attractive solution for helium-cooled PFCs that can handle high heat fluxes without excessive power requirements. Key questions and trade-offs for concept development have been identified, including: how to achieve high radiation for reduction of SOL power without core performance degradation; whether power exhaust can be well-controlled in a double null plasma; mechanical design and materials challenges of high-heat flux PFCs; and control of material erosion, redeposition and tritium retention.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Analysis of RF Sheath-Driven Tungsten Erosion at RF Antenna in the WEST Tokamak

This study applies the newly developed STRIPE (Simulated Transport of RF Impurity Production and Emission) framework to analyze tungsten (W) erosion at RF antenna structures in the WEST tokamak. STRIPE integrates SolEdge3x for edge plasma backgrounds, COMSOL for 3D RF sheath potentials, RustBCA for sputtering yields, and GITR for impurity transport and ion energy–angle distributions. Building on prior work by Kumar et al. (2025) Nuclear Fusion, 65, 076039, which validated STRIPE for WEST ICRH discharge #57877, the present study provides a spatially resolved assessment of gross W erosion at both Q2 antenna limiters under ohmic and ICRH conditions. Simulations using 2D SolEdge3x profiles in COMSOL capture rectified sheath potentials exceeding 300 V, leading to strong upper-limiter localization. Both poloidal and toroidal asymmetries are observed and attributed to RF sheath effects, with modeled erosion patterns deviating from experiment—highlighting sensitivity to sheath geometry and plasma resolution. Erosion is driven primarily by high-charge-state oxygen ions (O6+–O8+), while D+ plays a negligible role. Assuming a plasma composition of 1% oxygen and 98% deuterium, STRIPE predicts a 30-fold increase in gross W erosion from ohmic to ICRH phases, consistent with a >25-fold rise in W-I (400.9 nm) brightness. Quantitative agreement is within 5% in the ohmic phase and 30% under ICRH, demonstrating predictive capability. Importantly, the study shows that the magnitude of ICRH-driven W erosion depends strongly on the concentration of light impurities (O, B, N, C), which drive sputtering through high charge states. Cleaner plasma conditions with reduced impurity content are therefore expected to substantially mitigate antenna W sources in WEST and other toroidal fusion devices. These findings establish STRIPE as a predictive framework for RF-induced plasma–material interactions and support its application to reactor-scale antenna design.

Kumar, Atul [ORNL] (ORCID:0000000277475210)↗
DOE OSTI2021

Two-color scattering for the measurement of neutrals at the edge of fusion devices

Laser two-color scattering (TCS) is proposed to detect the neutral species in the edge of fusion devices, namely, tokamaks. Here, TCS uses two wavelengths to probe both the laser Rayleigh scattering and Thomson scattering of the neutral-electron bath, with emphasis on neutral density measurements such as that of hydrogen and deuterium. Modeling of the Rayleigh scattering of tokamak neutral species under various plasma conditions (electron density and temperature) shows that, with an appropriate filtering of the Thomson signal and by going to ultraviolet-region wavelengths, identification of the Rayleigh signal can be achieved. Photon count and signal fractions are calculated in two test cases, one in the midplane region of the National Spherical Torus Experiment and one in the divertor region of DIII-D. An uncertainty analysis and discussion of the feasibility of the TCS diagnostic is also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Interpreting AI for fusion: An application to plasma profile analysis for tearing mode stability

Artificial intelligence models have demonstrated strong predictive capabilities for various instabilities in fusion devices such as Tokamaks, including tearing modes (TM), edge localized modes, and disruptive events, but their opaque nature raises concerns about safety and trustworthiness when applied to fusion power plants. Here, we present a physics-based interpretation framework using a TM prediction model as a demonstration that is validated through a dedicated DIII-D TM avoidance experiment. By applying Shapley analysis, we identify how profiles such as rotation, temperature, and density contribute to the model's prediction of TM stability. Our analysis shows that in our experimental scenario, core electron temperature and rotation peaking play the primary role in TM stability, while density changes have smaller effects on stability. We show that off-axis ion temperature stabilizes TMs, suggesting that off-axis neutral beam heating can further stabilize this scenario. This work presents a generalizable ML-based event prediction methodology, from training to physics-driven interpretation, bridging the gap between physics understanding and opaque ML models.

Farre-Kaga, Hiro J. [Princeton Univ., NJ (United S↗
DOE OSTI2023

Effect of resonant magnetic perturbations including toroidal sidebands on magnetic footprints and fast ion losses in HL-2M

Externally applied resonant magnetic perturbations (RMPs), generated by magnetic coils located outside the plasma (referred to as RMP coils), provide an effective way to control the edge localized mode (ELM) in tokamak devices. Due to the discrete nature of the toroidal distribution of these window-frame coils, toroidal sidebands always exist together with the fundamental harmonics designed for ELM control. In this work, the MARS-F code (Liu et al 2000 Phys. Plasmas7 3681) is applied to investigate the detailed features of the RMP spectra considering both the dominant harmonic (n = 2) and the associated sideband (n = 6), and the impact of the combined fields on magnetic footprints as well as on the fast ion losses for a reference double-null scenario in the HL-2M device. It is found that the sum of the n = 2 and n = 6 RMP fields splits the footprint and widens the footprint area, as compared to the single-n (n = 2) harmonic case. The resistive plasma response breaks the up–down symmetry of the footprint pattern on the outer divertor plates, which is otherwise symmetric assuming vacuum RMP fields. Considering fast ion losses, a threshold value exists for the initially launched radial position of test particles, as well as for the RMP coil current, before the loss occurs. When the threshold criterion is satisfied, the combined n = 2 and n = 6 RMP fields enhance the fast ion loss rate by , as compared to that of the n = 2 component alone. These results illustrate the important role of the sideband of RMP fields on the magnetic footprints and fast ion losses in tokamak plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Comparison of plasma response models for RMP effects on the divertor and scrape-off layer in KSTAR

Resonant magnetic perturbations (RMPs) are beneficial for control of edge localized modes (ELMs) in tokamaks. Nevertheless, a side effect of RMPs is the appearance of helical striations in the particle and heat loads onto divertor targets. The extent and field line connection of these striations is significantly altered by the plasma response to external perturbations. For an ELM suppressed high confinement (H-mode) plasma at KSTAR, magnetic footprints are computed by FLARE based on plasma response from GPEC, MARS-F, M3D-C1 and JOREK with substantial differences in the resulting footprints (from 2 cm to 14 cm). This is reflected in EMC3-EIRENE simulations of the resulting heat loads: it is found that either the peak value or the extent of the striations appears to be overestimated compared to infrared camera measurements. Reasonable agreement can only be achieved for the smallest footprint for lower input power and lower cross-field transport, or for higher upstream density and radiative power losses.

Frerichs, H. [Univ. of Wisconsin, Madison, WI (Uni↗
DOE OSTI2025

Effects beyond ideal MHD on stability of wide and enhanced pedestal regimes in NSTX

Stability of edge-localized modes (ELMs) in spherical tokamaks is explored using the extended MHD model. Linear NIMROD simulations have been performed for three NSTX discharges 132543, 132588, and 141133, to investigate the role of resistivity, diffusivity, and shear flows on the onset of ELMs. The first discharge represents the wide pedestal regime and the later two discharges represent ELM-free enhanced pedestal H-mode. We first present the effect of toroidal rotation shear and find a flow shear destabilizing effect in these NSTX discharges. Simulations are also extended to include the two-fluid and ion gyroviscosity effects. Simulations show that the flow shear can shift the mode spectrum and alter the critical condition of ELM onset. We also uncover that ELM onset prediction in spherical tokamaks requires effects beyond MHD, in particular gyroviscosity and diamagnetic terms could stabilize Peeling-Ballooning modes consistent with the experimental observation of ELM-free regimes in NSTX. The findings give new insight into the nature of the interplay between resistivity, flows, and diamagnetic stabilization in ELM suppression and have potential applications to ELM control schemes in NSTX-U and next-generation spherical tokamaks. This study identifies the essential physical effects that must be included in future predictive and validation simulations.

NSTX↗