Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plasma fusion”

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 1,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

The Radiative Magneto-Hydrodynamics of Colliding Flow: Instabilities, Reconnection and Exoplanet Atmosphere Connections

Hypersonic flows of magnetized plasmas naturally produce strong shocks when obstructions or self‐interactions (flow collisions) occur. Such complex flows occur in a wide variety of HEDP settings such as Z‐pinches and laser driven implosion experiments that include magnetic fields via explosive Helmholtz coils (MIFEDS) When thermal energy behind the shocks can be lost to optically thin radiation, the flow dynamics becomes even more complex, admitting a range of new instabilities such as the Nonlinear Thin Shell Instability (NTSI). Because of their inherent non‐ linearity, multi‐dimensionality and time‐dependence (i.e. turbulence), the ability to articulate the underlying physics of colliding hypersonic magnetized flows represents a significant frontier in HEDP studies. The addition of radiative cooling to these flows propels the study even further into the domain of extreme plasma physics. Exploration of these forms of HED plasmas is of interest to basic plasma physics as well as of interest for their application to Inertial Confinement Fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Plasma pressure profiles in a sheared-flow-stabilized Z-pinch

We report the plasma pressure reached inside the central plasma column of a sheared-flow-stabilized Z-pinch using Thomson scattering measurements. Building on previously reported experimental results and the analysis methods established for the high temperature and moderate density plasmas generated on the FuZE device, we show evidence of a central plasma region with higher electron temperature and density, which is consistent with a pinch behavior. Elevated electron temperatures up to 2.25 ± 0.8 keV and densities up to (4.9±0.2)×1017 cm−3 are observed to temporally coincide with the fusion neutron production from the plasma. Reconstructed plasma pressure profiles highlight the presence of a several millimeter-wide column with elevated pressure whose location varies shot-to-shot. The plasma pressure rises as neutron production from the deuterium plasma increases, reaching a peak value of 2.6 kBar. This peak value is consistent with a radially force-balanced pinch equilibrium model based on the measured ∼320 kA pinch current. Complete datasets were obtained at two axial locations, 10 and 20 cm axial position from the tip of the central electrode, which corroborate the estimated neutron source axial lengths.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

CRADA Final Report: CRADA Number NFE-22-09330 with General Fusion

General Fusion is developing a magnetized target fusion (MTF) approach that involves compressing an initial magnetically confined plasma inside a cavity formed in liquid metal. This approach builds from concepts initially developed under the Linus program at the U.S. Naval Research Laboratory and combines it with advances from compact toroid experiment (CTX) and sustained spheromak plasma experiment (SSPX) in compact toroid plasmas and coaxial Marshall gun systems. Modeling the tokamak during compression is central to designing a successful MTF device. The plasma is formed by coaxial helicity injection in the General Fusion device. Immediately after formation, the plasma has a diverted tokamak configuration with a single null. As the wall moves inwards, the plasma is repelled from the conducting surface and driven inwards by currents induced by its magnetic field in the liquid metal wall. As the liquid metal closes (or bridges) the opening of the coaxial plasma injector, the magnetic field topology alters to remove the null. Due to this, the plasma moves from a diverted to a wall-limited configuration. The liquid metal liner continues to close in and change shape, reducing in radius by a factor of ten at the peak of plasma compression. A model of the MTF plasma must be able to handle this continually varying geometry, and to be predictive, it must faithfully include the real imperfections arising in the process. In this project, we pursued a Monte Carlo approach to closures for MHD by computing kinetic electron trajectories in an MHD plasma background from simulations of GF devices. This requires enhancing the capabilities of the KORC-T code for running large ensembles of kinetic trajectories by porting it to GPU architectures and enabling workflows for large ensembles on OLCF machines. With these capabilities, it is possible to produce a large library of kinetic calculations of electron orbits evolving in plasma configurations spanning the magnetic configurations and plasma density profiles, including non-axisymmetry, arising in the General Fusion’s existing PI3 spherical tokamak device. Using ensembles will capture particles passing a single point in space in a given magnetic configuration, and the entire dataset will cover a range of global magnetic field geometries. By sampling around many starting points, this dataset will capture the spatial dependence of the plasma parameters. From this large dataset, it is possible to produce a reduced model for the kinetic effects not captured in MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Effect of outer divertor leg detachment on the high field side scrape-off layer in DIII-D and ASDEX Upgrade

In this work, evidence is presented that detachment of the outer divertor leg leads to a reduction of electron density and neutral pressure in the high-field side (HFS) scrape-off layer (SOL) of the ASDEX Upgrade (AUG) and DIII-D tokamaks with ion $B$ x $∇B$ drift directed toward the X-point (favorable configuration). These results are observed across multiple diagnostics and without the use of impurity seeding to reach detachment. In AUG, outer divertor leg detachment correlates with a decrease in electron density near the separatrix at the inner midplane, measured with HFS reflectometry. A concurrent reduction in inner divertor density and neutral pressure at the inner target is observed using divertor Thomson scattering and neutral pressure gauges. These effects are present in both L- and H-mode plasmas. In DIII-D, a similar reduction is detected through analysis of line-integrated hydrogenic emission measured by multiple diagnostics in the HFS SOL close to the separatrix. The consistent trends in both devices indicate that high electron density and strong hydrogenic emissivity in the HFS SOL are common features of H-mode plasmas in the favorable configuration, independent of wall material and divertor geometry. In L-mode plasmas, the reduction in electron density and neutral pressure is not observed in DIII-D, possibly due to differences in wall material. These results emphasize the importance of the divertor state in determining the two-dimensional neutral distribution and edge fueling.

detachment↗
DOE OSTI2021

Lithium, a path to make fusion energy affordable

In this tutorial article, we review the technological, physics, and economic basis for a magnetic fusion device utilizing a flowing liquid lithium divertor (molten metal velocity in the range of cm/s) and operating in a low-recycling plasma regime. When extrapolated to magnetic fusion reactor scale, the observed effects of a liquid lithium boundary on recycling reduction, confinement increase, and anomalous heat transport mitigation may offer a fundamentally distinct and promising alternative route to fusion energy production. In addition, this lithium-driven low recycling regime could accelerate fusion’s commercial viability since such a device would be smaller, dramatically decreasing plant and electricity costs if all technological complexities are solved. First, the theoretical basis of the energy confinement and fusion performance as well as the related possibilities of low recycling regimes driven by flowing lithium plasma-facing components are reviewed. Then the paper emphasizes the technological obstacles that need to be overcome for developing the necessary systems for such a flowing liquid lithium solution at reactor scale and details how many of these have been overcome at laboratory and/or proof-of-concept scale. Finally, the current and planned scientific and engineering endeavors being performed at the University of Illinois at Urbana-Champaign regarding this alternative reactor option are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Autoregressive long-horizon prediction of plasma edge dynamics *

Accurate modeling of scrape-off layer (SOL) and divertor-edge dynamics is vital for designing plasma-facing components in fusion devices. High-fidelity edge fluid/neutral codes such as SOLPS-ITER capture SOL physics with high accuracy, but their computational cost limits broad parameter scans and long transient studies. We present transformer-based, autoregressive surrogates for efficient prediction of 2D, time-dependent plasma edge state fields. Trained on SOLPS-ITER spatiotemporal data for the KSTAR tokamak, the surrogates forecast electron temperature, electron density, and radiated power over extended horizons. We evaluate model variants trained with increasing autoregressive horizons (1–100 steps) on short- and long-horizon prediction tasks. Longer-horizon training systematically improves rollout stability and mitigates error accumulation, enabling stable predictions over hundreds to thousands of steps and reproducing key dynamical features such as the motion of high-radiation regions. Measured end-to-end wall-clock times show the surrogate is orders of magnitude faster than SOLPS-ITER, enabling rapid parameter exploration. Prediction accuracy degrades when the surrogate enters physical regimes not represented in the training dataset, motivating future work on data enrichment and physics-informed constraints. Overall, this approach provides a fast, accurate surrogate for computationally intensive plasma edge simulations, supporting rapid scenario exploration, control-oriented studies, and progress toward real-time applications in fusion devices.

autoregressive deep learning↗
DOE OSTI2023

Implementation of extreme ultraviolet spectroscopy on a sheared-flow-stabilized Z pinch

A diagnostic for extreme ultraviolet spectroscopy was fielded on the sheared-flow-stabilized (SFS) fusion Z-pinch experiment (FuZE-Q) for the first time. The spectrometer collected time-gated plasma emission spectra in the 5–40 nm wavelength (30–250 eV) range for impurity identification, radiative power studies, and for plasma temperature and density measurements. The unique implementation of the diagnostic included fast (10 ns risetime) pulsed high voltage electronics and a multi-stage differential pumping system that allowed the vacuum-coupled spectrometer to collect three independently timed spectra per FuZE-Q shot while also protecting sensitive internal components. Analysis of line emission identifies oxygen (N-, C-, B-, Be-, Li-, and He-like O), peaking in intensity shortly after maximum current (>500 kA). This work provides a foundation for future high energy spectroscopy experiments on SFS Z-pinch devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Direct Measurement of ICRF-Enhanced Plasma Potentials on WEST Using Reciprocating Emissive Probes

An extensive documentation of ICRF-enhanced plasma potentials has been conducted over two experimental campaigns on the WEST tokamak using reciprocating emissive probes magnetically connected to two ICRF antennas. The collected data spans a wide range of antenna electrical settings (coupled power, toroidal phasing, left–right power balance) and plasma parameters (density at the antenna limiter above and below the lower hybrid resonance, plasma current, minority fraction). By scanning the edge safety factor across multiple probe plunges, the magnetic connection between the probe and the antenna varied, enabling the construction of a 2D map of the plasma and floating potentials around an active ICRF antenna. This dataset will be used to validate RF simulation tools equipped with the sheath boundary condition and used to predict RF rectified potentials and ICRF-induced impurity sputtering in future machines. This paper presents the diagnostic and some initial measurements, while the rest will be reported elsewhere.

Diab, Raymond [Massachusetts Inst. of Technology (↗
DOE OSTI2025

Enhanced pedestal transport driven by edge collisionality on Alcator C-Mod and its role in regulating H-mode pedestal gradients

Experimental measurements of plasma and neutral profiles across the pedestal are used in conjunction with 2D edge modeling to examine pedestal stiffness in Alcator C-Mod H-mode plasmas. Enhanced D α experiments on Alcator C-Mod observed pedestal degradation and loss in confinement below a critical value of net power crossing the separatrix, P net = $P^{crit}_{net}$ ≈ 2.3 MW, in the absence of any external fueling. New analysis of ionization and particle flux profiles reveal saturation of the pedestal electron density, $n^{ped}_{e}$, despite continuous increases in ionization throughout the pedestal, inversely related to P net . A limi to the pedestal $\nabla$n e emerges as the particle flux, Γ D , continues to grow, implying increases in the effective particle diffusivity, D eff . This is well-correlated with the separatrix collisionality, $v^*_{sep}$ and a turbulence control parameter, α t , implying a possible transition in type of turbulence. The transition is well correlated with the experimentally observed value of $P^{crit}_{net}$. SOLPS-ITER modeling is performed for select discharges from the power scan, constrained with experimental electron and neutral densities, measured at the outer midpane. The modeling confirms general growth in D eff , consistent with experimental findings, and additionally suggests even larger growth in Χ e at the same $P^{crit}_{net}$.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Fast ion relaxation in ITER mediated by Alfvén instabilities

Abstract We address the critical issue for future burning plasmas of whether high-energy fusion products or auxiliary heating-beam ions will be confined for a sufficiently long time to compensate for thermal plasma energy losses. This issue can be mitigated by one of the most deleterious collective phenomena—the instability of low, sub-cyclotron frequency Alfvén eigenmodes (AEs), such as toroidicity-induced AEs and reversed-shear AEs in the ITER steady-state scenario. Using a revised quasi-linear (QL) theory applied to energetic particle (EP) relaxation in the presence of AEs, we find that the AE instabilities can affect both neutral beam ions and alpha particles, although the resulting fast ion transport is expected to be modest if classical particle slowing down is assumed. On the other hand, the QL theory predicts that the AE amplitudes will be enhanced by the background microturbulence, although this topic remains outside our scope due to the significant numerical effort required to evaluate these effects. We report our results for EP relaxation dynamics obtained utilizing several tools: (i) a comprehensive linear stability study of the sub-cyclotron Alfvénic spectrum as computed by ideal magnetohydrodynamic NOVA simulations for the AE eigenproblem, (ii) drift kinetic NOVA-C calculations for wave–particle interaction and AE growth/damping rates, and (iii) predictive QL modeling coupled with the global transport code TRANSP to assess the EP relaxation on the equilibrium timescale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

The Development of Kinetic and Radiation Hydrodynamics Modeling of Thermonuclear Burn Propagation in Isochoric p - 11 B Through the Support of the INFUSE Program

The report summarizes DOE INFUSE-supported work between HB11 Energy and the University of Rochester’s TriForce Institute to improve computational modeling of advanced fusion fuels, especially proton–boron-11 (p- 11 B). The project extended the TriForce particle-in-cell/Monte Carlo collision code to include physics needed for dense, high-temperature p- 11 B burn studies, including p- 11 B fusion reactions, three-alpha-particle reaction products, relativistic Coulomb collisions, large-angle nuclear scattering, bremsstrahlung radiation, inverse bremsstrahlung absorption, and photon transport. The upgraded models were verified against focused physics tests and against known deuterium–tritium burn behavior. The study then used one-dimensional spherical simulations to estimate the conditions required for thermonuclear burn propagation in isochoric p- 11 B fuel. The calculations found that burn propagation is possible in the model, but only under very extreme hot-spot conditions, such as about 7000 g/cm 3 at 500 keV or 9000 g/cm 3 at 300 keV for a 20-micron hot spot. These conditions are much more demanding than current demonstrated inertial confinement fusion hot spots. The report concludes that the INFUSE collaboration successfully advanced kinetic and radiation modeling capabilities for p- 11 B fusion and provided useful estimates of ignition requirements. However, the simulated fuel gains remain below what would be needed for practical inertial fusion energy, and further work is needed to reconcile differences among kinetic, radiation-hydrodynamic, and analytic models and to identify more achievable target designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Pedestal origin and extrapolation of high-density small edge-localised-modes peak parallel energy fluence in ITER and SPARC

Experimental analysis and simulations with the BOUT++ code show that small edge-localised modes (ELMs) in reactor-relevant high-density regimes originate in a region close to the separatrix and only marginally perturb the pedestal structure. The measured divertor peak parallel energy fluence (ε ∥,peak ) for a database of small ELM scenarios in DIII-D and ASDEX Upgrade can be reproduced, within 40 % accuracy on average, if an ad hoc modification of the Eich peak parallel ELM energy fluence model is applied to account for the small ELM pedestal birth location. This allows for first-order extrapolation of small-ELM divertor ε ∥,peak to ITER and SPARC, resulting in values that satisfy the nominal melting threshold of tungsten monoblocks of 12 MJ m −2 . The findings reported in this study, both via modelling and direct measurements, constitute a step forward in assessing small ELMs in high edge-collisionality scenarios as a viable plasma regime for the operation of next-generation fusion machines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Propagation of radio frequency waves through turbulent plasmas

The practical and economic viability of tokamak fusion reactors depends, in a significant way, on the efficiency of radio frequency (RF) waves to deliver energy and momentum to the plasma in the core of the reactor. The RF electromagnetic waves, excited by antenna structures placed near the wall of a tokamak, have to propagate through the turbulent edge plasma along their path to the core of the fusion device. In present day experiments, the radial width of the edge region and scrape-off layer is of the order of a few centimeters. In ITER, and in future fusion reactors, this width will be of the order of tens of centimeters. Any effects on RF waves due to plasma turbulence have to be properly understood in order to optimize the delivery of RF energy and momentum into the core. This paper is on a multi-pronged, theoretical and computational, approach that is being pursued to quantify the effect of edge plasma turbulence on the propagation of RF waves. The theoretical and analytical models are based on solutions of the Faraday-Ampere equation in a magnetized plasma and on the Kirchhoff tangent plane approximation. An effective medium approach has been developed so as to approximate the permittivity of a turbulent plasma by analytical expressions. The computations are being carried out with a newly developed code ScaRF that is based on the finite difference finite domain technique for solving Maxwell's equations. The plasma permittivity can be assigned as desired. The code is being used to validate the analytical and theoretical models and to evaluate their limitations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Reduced mixing in inertial confinement fusion with early-time interface acceleration

In inertial confinement fusion (ICF) implosions, the interface between the cryogenic DT fuel and the ablator is unstable to shock acceleration (the Richtmyer-Meshkov instability, RM) and constant acceleration (Rayleigh-Taylor instability, RT). Instability growth at this interface can reduce the final compression, limiting fusion burnup. If the constant acceleration is in the direction of the lighter material (negative Atwood number), the RT instability produces oscillatory motion that can stabilize against RM growth. Theory and simulations suggest this scenario occurred at early times in some ICF experiments on the National Ignition Facility, possibly explaining their favorable performance compared to one-dimensional simulations. In conclusion, this characteristic is being included in newer, lower adiabat designs, seeking to improve compression while minimizing ablator mixing into the fuel.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Implementation of a High-Speed Multichannel Data Acquisition System for Magnetic Diagnostics and Plasma Centroid Position Control in ISTTOK

In tokamak and other fusion devices, magnetic control is the main tool that allows to regulate the plasma current, position and shape; it is in charge of actuating the desired plasma current waveform, steering the plasma position to a given set point and maintain the plasma shape close to a prescribed plasma equilibrium. This work describes the application of several physics concepts and computational tools in order to obtain a novel optimal controller for the plasma centroid position, which has been implemented and tested in the real-time plasma control system at the ISTTOK tokamak. A key point for the development of the new control system was the installation of a recently upgraded hardware, that numerically integrates in real-time the magnetic probes signals.

poloidal field coils↗
DOE OSTI2021

The role of edge plasma parameters in H-mode density limit on the JET-ILW

A study of a dataset of JET H-mode plasma with the Be/W ITER-like wall (JET-ILW) shows that reaching the edge MHD ballooning limit leads to confinement degradation. However, unlike JET plasmas with a carbon wall (JET-C), the JET-ILW plasmas stay in a marginal dithering phase for a relatively long period, associated with a higher (≈20%) H-mode density limit (HDL) than JET-C equivalents. This suggests that ITER could be operated in H-mode with higher density than the scaling based on carbon wall devices, but likely with a dithering phase plasma with lower confinement. A new, reliable estimator for JET E r, min has been derived by combining HRTS measurements of pedestal gradient and edge-SOL decay lengths. JET radial E r ETB wells are observed in the range of -15 to -60 kV m -1 in high performance H-modes, consistent with previous CXRS results in ASDEX Upgrade. The results imply that a higher positive E × B shear in the near SOL plays a role in sustaining a marginal phase in JET-ILW which leads to a higher HDL than that in JET-C. The results of the JET-ILW dataset show agreement with the Goldston finite collisionality HD model for SOL broadening at high collisionality. A hypothesis for the dithering H-mode phase is proposed: as n e,SOL increases, ν *,SOL increases, SOL broadens, E r shear decreases, triggers L-mode; n e drops, ν *,SOL decreases, SOL becomes narrower, and E r shear increases, triggering H-mode, resulting in a cycle of H–L–H- oscillations. For burning plasma devices, such as ITER, operating just below the MHD limit for the dithering phase could be a promising regime for maximising core density, and fusion performance while minimising plasma-material interaction. The oscillatory signal during the dithering phase could be used as a precursor of undesirable plasma performance for control purposes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The Verification and Validation of a Magnetic Plasma Fluid Model Utilizing the MOOSE (Multiphysics Object Oriented Simulation Environment) Framework

As the goal of achieving fusion power on the grid comes closer to fruition, fully coupled multiphysics models of fusion devices will be crucial. Currently, there are two main approaches to developing these platforms: (1) loosely coupled, where one couples existing codes and solvers together through input and output parameters and data, and (2) tightly coupled, where one develops the necessary models within a singular, integrated framework. This work focuses on the latter approach for magnetically confined fusion devices by developing a fluid-based plasma-edge model within the Multiphysics Object Oriented Simulation Environment (MOOSE) Framework. This effort is coordinated with other efforts to develop, test, demonstrate, and deploy fusion relevant multiphysics capabilities including electromagnetics, particle-in-cell plasma, tritium transport, and fusion blanket design. This new model is an expansion of the MOOSE-based plasma application, Zapdos, which was originally formulated to model low-temperature, non-magnetized plasma processes. Verification, benchmarking, and validation studies have been conducted. Verification studies involved utilizing the method of manufactured solutions and comparing the convergence slope of a known solution to the theoretical slope. Benchmarking consists of comparisons to existing edge codes, namely BOUT++ and SOLEDGE3X. Validation efforts focused on comparisons against open-source data from the TCV tokamak.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Understanding the formation of a low-pressure pedestal in the presence of a strong internal transport barrier in DIII-D high β p plasmas

As a promising scenario for fusion reactors, the high poloidal-beta (β P ) scenario is characterized by a strong large radius internal transport barrier (ITB), which significantly enhances the overall confinement quality and the bootstrap current fraction for fully non-inductive operation. It is frequently observed that in the presence of a strong ITB, the pedestal height is lower and is accompanied by small edge localized modes (ELMs), which further improves the compatibility of a high performance core with an edge solution. A mechanism for the formation of the low pedestal is proposed in this paper. It is found that the strong ITB creates an off-axis bootstrap current to clamp the local safety factor q, and thus the magnetic shear in the outer core/pedestal region is increased. Gyrokinetic simulations with the CGYRO code show that the higher magnetic shear brings the experimental profiles into the range where the growth rate of drift-wave instabilities and thus transport is higher, and therefore a lower pedestal gradient is expected. Here, the combination of low pedestal and high magnetic shear further enhances the turbulent transport across the whole pedestal, consistent with power balance analysis. Such a positive feedback mechanism ultimately results in a lower pressure pedestal as observed in experiments. Under such a low pedestal, linear simulations with BOUT++ predict the growth rates of peeling–ballooning modes to be lower across the whole toroidal mode number spectra, and the nonlinear BOUT++ simulation exhibits lower saturated fluctuation intensity as well, consistent with the experimentally observed lower ELM size.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗