Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “magnetic 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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Dependence of ion-cyclotron range of frequencies instabilities on species mix and fast-ion distribution: I. Stability

Ion cyclotron emission (ICE) driven by fast ions has the potential to be a reactor relevant diagnostic of confined and escaping fast ions; moreover, controlled experiments in magnetic fusion devices can shed light on similar magnetospheric instabilities below and above the ion cyclotron frequency f ci . Dedicated DIII-D experiments studied instabilities below f ci and at cyclotron harmonics l f ci by injecting thirteen different neutral beam populations into plasmas with different thermal compositions of H, D, and 3He at five different values of toroidal field. High-frequency magnetic loops diagnose mode properties. Each beam type is individually injected for an approximate slowing-down time (~100 ms), resulting in a database of 2529 conditions that are analyzed in frequency bands corresponding to ICE harmonics and < f ci global Alfvén eigenmodes. The sub-cyclotron modes are more unstable at low field and in plasmas with large H concentration. For D injection, first harmonic ICE at f cD is more unstable in hydrogen plasmas but, for H injection, the converse is true (ICE at f cH is more unstable in D), suggesting that thermal same-species ions damp fundamental ICE. ICE stability also depends upon the pitch-angle anisotropy gradient and upon the ‘bump-on-tail’ gradient ∂f/∂v. Although the frequency of unstable ICE in these plasmas always occurs near the central cyclotron frequency, the precise value also depends on the Alfvén speed and beam direction, consistent with the idea that instability occurs for modes that simultaneously satisfy ω ≃ 2πlf ci , the magnetosonic dispersion relation, and the Doppler-shifted resonance condition.

electromagnetic ion cyclotron instabilities↗
DOE OSTI2022

Neutron time-of-flight detectors (nTOF) used at Sandia’s Z-Machine

Neutron time-of-flight (nTOF) detectors have been used on Sandia National Laboratories’ Z-Machine for inertial confinement fusion and magnetized liner fusion experiments to infer physics parameters including the apparent fuel-ion temperature, neutron yield, the magnetic-radius product (BR), and the liner rho-r. Single-paddle, dual-paddle, and co-axial scintillation nTOF detectors are used in axial lines-of-sight (LOS) and LOS that are 12° from the midplane. Detector fabrication, characterization, and calibration are discussed.

Chandler, G. A. (ORCID:0000000219640927)↗
DOE OSTI2024

Final CRADA Report – NFE-21-08693

TAE Technologies is developing a magnetic fusion energy concept known as the beam-driven field-reversed configuration (FRC) with the ultimate goal of developing a reactor for commercial electricity production capable of burning aneutronic pB11 fuel. To achieve the high plasma temperatures this requires, auxiliary radiofrequency (RF) heating will likely be needed. High Harmonic Fast Wave (HHFW) heating has been identified as a candidate RF heating scheme to overcome the unique challenges posed to RF heating by the FRC, including the large distance from the plasma edge to the last closed flux surface and a magnetic field profile with strength decreasing from edge to core and reversing sign at a null point inside the plasma. The purpose of this project was to develop the experimental capabilities to test HHFW on TAE’s C-2W device through the design of a phased array antenna and accompanying matching network. The design was performed by ORNL and informed by experiments with a prototype four-strap phased antenna-array that was manufactured and installed on the LArge Plasma Device (LAPD) at UCLA and simulations conducted with the Petra-M code under the purview of a previous INFUSE grant. The ORNL team completed the conceptual design of the antenna and matching network which was then handed off to the TAE Mechanical Design team. The design was then iterated on to ensure changes to the mechanical design did not interfere with the RF performance. This process is now complete, and, with mechanical design in hand, TAE is proceeding with plans for final integration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Validity condition for the local sheath impedance boundary condition and a nonlocal generalization

ICRF sheaths can cause unwanted interactions of high-power RF waves with material surfaces in magnetic fusion devices. In previous work, a local RF sheath impedance boundary condition (BC) was derived for use in ICRF codes together with a microscale (i.e., Debye or sheath width scale) model for obtaining the sheath impedance used in that BC. This local RF sheath BC matches the normal component of current and electrostatic potential across the sheath-plasma interface. Collapsing the matching conditions at the sheath-plasma interface to a BC depends on the assumption of scale separation, which can be violated when conditions along the local radius of curvature of the surface vary sufficiently rapidly. The validity condition is explored in this contribution, with special attention to the case where the magnetic field approaches being tangent to the surface. When the local sheath BC no longer applies, a non-local sheath BC is developed under the assumption of a more relaxed scale separation assumption. It is shown that the non-local sheath BC reduces to the previous local sheath BC under appropriate conditions. Furthermore, a surface-integrated sheath admittance parameter describes the 2D physics in the new BC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Avoiding fusion plasma tearing instability with deep reinforcement learning

For stable and efficient fusion energy production using a tokamak reactor, it is essential to maintain a high-pressure hydrogenic plasma without plasma disruption. Therefore, it is necessary to actively control the tokamak based on the observed plasma state, to manoeuvre high-pressure plasma while avoiding tearing instability, the leading cause of disruptions. This presents an obstacle-avoidance problem for which artificial intelligence based on reinforcement learning has recently shown remarkable performance. However, the obstacle here, the tearing instability, is difficult to forecast and is highly prone to terminating plasma operations, especially in the ITER baseline scenario. Previously, we developed a multimodal dynamic model that estimates the likelihood of future tearing instability based on signals from multiple diagnostics and actuators. Here we harness this dynamic model as a training environment for reinforcement-learning artificial intelligence, facilitating automated instability prevention. We demonstrate artificial intelligence control to lower the possibility of disruptive tearing instabilities in DIII-D, the largest magnetic fusion facility in the United States. The controller maintained the tearing likelihood under a given threshold, even under relatively unfavourable conditions of low safety factor and low torque. In particular, it allowed the plasma to actively track the stable path within the time-varying operational space while maintaining H-mode performance, which was challenging with traditional preprogrammed control. This controller paves the path to developing stable high-performance operational scenarios for future use in ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Beam-ion Studies in NSTX and NSTX-U (Final Technical Report)

The confinement of high energy "fast ions" is crucial for the success of magnetic fusion as a practical energy source. The research performed by UC Irvine on NSTX and NSTX-U provided new information about this important topic in the configuration known as a "spherical tokamak." Injected neutral beams provided the fast ions (also known as "beam ions"). There were five overarching goals of the research. One objective was to measure the confinement of beam ions in order to ascertain if they were accomplishing their desired purpose of transferring their energy to the bulk plasma. A second goal was to better understand instabilities that are driven unstable by the fast ions. A related goal was to use the similarities and differences between NSTX and the DIII-D conventional tokamak to better understand the fast-ion driven instabilities. A fourth goal was to develop new instruments to measure fast ions and techniques that facilitate analysis of the data. The fifth objective was to measure and understand acceleration of beam ions by RF waves. Much was accomplished in all of these five areas. In the first, it was shown that the large magnetic moment of the beam ions did not harm confinement but instabilities known as the "sawtooth" and "long-lived mode" do. In the second area, much attention was devoted to the fast-ion instabilities known as "Alfv&eacute;n eigenmodes." In particular, Alfv&eacute;n eigenmode "avalanches" can cause ~ 30% of the fast ions to be lost in a single explosive burst. The comparative studies with DIII-D showed that an instability discovered on NSTX is also important in conventional tokamaks. In the fourth category, several instruments were developed, a number of effects that complicate interpretation of the data were understood, and progress toward a new method to infer the fast-ion distribution function from the data was made. In the fifth category, the measured profile of accelerated fast ions was initially much broader than theoretical predictions but subsequent improvements in the theoretical modeling achieved good agreement with the data.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Hall interchange instability as a seed for helical magneto-Rayleigh–Taylor instabilities in magnetized liner inertial fusion Z-Pinches scaled from Z-Machine parameters to a next generation pulsed power facility

Magnetized liner inertial fusion (MagLIF) is a magneto-inertial-fusion concept that is studied on the 20-MA, 100-ns rise time Z Pulsed Power Facility at Sandia National Laboratories. Given the relative success of the platform, there is a wide interest in studying the scaled performance of this concept at a next-generation pulsed-power facility that may produce peak currents upward of 60 MA. An important aspect that requires more research is the instability dynamics of the imploding MagLIF liner, specifically how instabilities are initially seeded. It has been shown in magnetized 1-MA thin-foil liner Z-pinch implosion simulations that a Hall interchange instability (HII) effect can provide an independent seeding mechanism for helical magneto-Rayleigh–Taylor instabilities. Here in this paper, we explore this instability at higher peak currents for MagLIF using 2D discontinuous Galerkin PERSEUS simulations, an extended magneto-hydrodynamics code, which includes Hall physics. Our simulations of scaled MagLIF loads show that the growth rate of the HII is invariant to the peak current, suggesting that studies at 20-MA are directly relevant to 60-MA class machines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Experimental evaluation of thermal-fluids performance of Helium-Cooled Flat Plate (HCFP) divertor

Various helium (He)-cooled solid-tungsten (W) divertor concepts have been proposed for long-pulse magnetic fusion energy reactors. Among these concepts, the He-cooled flat plate divertor (HCFP) modules have the largest plasma-facing surface area of ~0.2 m 2 . Simulations have shown that the most recent version of the design can withstand heat fluxes as great as 8 MW/m 2 . Earlier experimental studies of a single shortened HCFP cooling unit with a slot length of 7.6 cm used air at ambient temperature and pressures below 0.6 MPa. Here, we present initial experimental studies of a copper alloy and steel test section modeling a shortened HCFP cooling unit using He at prototypical pressure of 10 MPa, inlet temperatures T i ≤ 200 °C and steady-state incident heat fluxes q" ≤ 1.2 MW/m 2 . Results for Nusselt number Nu as a function of Reynolds number Re were obtained for Re = 1.2×10 4 –3.4×10 4 , and used to develop a Nu(Re) correlation and validate numerical models of the test section using commercial computational fluid dynamics (CFD) software. Simulations with this model are performed to evaluate the effect of the shortened slot. Furthermore, these analyses are used to estimate the thermal-fluids performance of the HCFP under prototypical conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Estimation of the Thermal-Fluids and Thermal-Structural Performance of Helium-Cooled Modular Finger-Type Divertors

Over the past decade, our group has investigated the thermal-fluid performance of the helium-cooled modular divertor with multiple jets (HEMJ) and a simplified “flat” design of the HEMJ for long-pulse magnetic fusion energy (MFE) reactors. Experimental studies were performed in a helium (He) loop at the prototypical pressure of 10 MPa, nearly prototypical He temperatures and incident heat fluxes using test sections made from stainless steel and tungsten alloys. Correlations for average Nusselt numbers and pressure loss coefficients were developed from the data and are used to validate computational fluid dynamics (CFD) models. Here, this work presents updated thermal-fluids performance curves based on these correlations that estimate the maximum heat flux that can be accommodated by the plasma-facing surface and coolant pumping power requirements at prototypical operating conditions. Thermal-structural performance curves developed from ITER structural design criteria are introduced, which include protection against ductile and non-ductile failure, ratcheting fatigue, and creep fatigue. The performance design curves for these finger-type divertors demonstrate that the “flat” design, with a significantly less complicated geometry than the HEMJ, has thermal-fluid and thermal-structural performance comparable to the original HEMJ concept, and may be superior in terms of non-ductile failure criteria.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Quasi-optical beam tracing module development for millimeter-wave high-wavenumber collective scattering on the NSTX-U and EAST tokamaks

A Python3-based beam tracing code utilizing Quasi-Optics has been developed to track both incident and receiving beams in high-k collective millimeter wave scattering systems within magnetic fusion plasmas. In contrast to existing ray tracing codes that solely consider refraction, this beam tracing code incorporates diffraction phenomena, providing a more comprehensive calculation. Here, this enhanced capability allows for a more accurate calculation of the scattering volume and spatial resolution in high-k collective scattering systems, crucial for evaluating system performance and facilitating data analysis. Unlike Geometrical Optics, Quasi-Optics employs the complex eikonal method, representing a Gaussian beam as a collection of coupled rays to accurately preserve diffraction characteristics. The developed code is intended for application in NSTX-Upgrade and EAST high-k beam tracing analyses, targeting frequencies of 693 GHz and 270 GHz, respectively. The high-k system's primary objective is the observation of electron-scale instabilities. Employing a symplectic integrator, the code ensures numerical accuracy, assessed through the conservation of the Hamiltonian. With its precision and efficiency, the code facilitates rapid inter-shot analyses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A scoping study of far-SOL main-wall protection limiters for steady-state operation of compact pilot plant tokamaks

We present a novel method for handling steady-state heat fluxes incident on the main wall of pilot plant-scale magnetic fusion devices, based on the utilization of protection limiters in the far scrape-off layer (SOL). This method helps avoid large plasma-wall gaps, without excessively compromising blanket performance. We present an optimization algorithm for determining the appropriate size and scale of these protection limiters given (1) probability distributions of SOL plasma parameters and (2) assumed risk tolerance. As part of this optimization, we have developed an analytic description of parallel heat fluxes across limiter shadows, and an objective cost function (the ‘Far-SOL Marginal Cost’) to quantify the impact that different main-wall thermal management design choices have on reactor capital cost. Applying the model to a midscale fusion pilot plant concept shows that making use of far-SOL protection limiters can reduce capital costs on the order of $500 M, relative to naively increasing the plasma-wall gap. Our analysis demonstrates that the far-SOL power decay length is the highest-leverage plasma assumption for thermal loading of the first wall, and the primary cost driver for main wall thermal management. The relative cost efficiency of protection limiters increases as assumptions on the far-SOL heat flux become more pessimistic. The concepts described in this paper motivate the further development of far-SOL protection limiters as part of larger efforts to design economical core-edge-wall compatible solutions for a fusion pilot plant.

Design under uncertainty↗
DOE OSTI2022

A conservative implicit-PIC scheme for the hybrid kinetic-ion fluid-electron plasma model on curvilinear meshes

We report that the hybrid kinetic-ion fluid-electron plasma model is widely used to study challenging multi-scale problems in space and laboratory plasma physics. Here, a novel conservative scheme for this model employing implicit particle-in-cell techniques is extended to arbitrary coordinate systems via curvilinear maps from logical to physical space. The scheme features a fully non-linear electromagnetic formulation with a multi-rate time advance - including sub-cycling and orbit-averaging for the kinetic ions. By careful choice of compatible particle-based kinetic-ion and mesh-based fluid-electron discretizations in curvilinear coordinates, as well as particle-mesh interpolations and implicit midpoint time advance, the scheme is proven to conserve total energy for arbitrary curvilinear meshes. In the electrostatic limit, the method is also proven to conserve total momentum for arbitrary curvilinear meshes. Although momentum is not conserved for arbitrary curvilinear meshes in the electromagnetic case, it is for an important subset of Cartesian tensor-packed meshes. The scheme and its novel conservation properties are demonstrated for several challenging numerical problems using different curvilinear meshes, including a merging flux-rope simulation for a space weather application, and a helical m = 1 mode simulation for magnetic fusion energy application.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Effect of alloying on intrinsic ductility in WTaCrV high entropy alloys

Tungsten (W) exhibits desirable properties for extreme applications, such as the divertor in magnetic fusion reactors, but its practicality remains limited due to poor formability and insufficient irradiation resistance. In this work, we study the intrinsic ductility of body-centered cubic WTaCrV based high entropy alloys (HEAs), which are known to exhibit excellent irradiation resistance. The ductility evaluations are carried out using a criterion based on the competition between the critical stress intensity factors for emission (K Ie ) and cleavage (K Ic ) in the {110} slip planes and {110} crack planes, which are evaluated within the linear elastic fracture mechanics framework and computed using density functional theory. The results suggest that increasing the alloying concentrations of V and reducing the concentrations of W can significantly improve the ductility in these HEAs. The elastic anisotropy for these HEAs is analyzed using the Zener anisotropy ratio and its correlation with the concentration of W in the alloys is studied. Results indicate that these alloys tend to be fairly isotropic independently from the concentration of W in them. The computed data for the elastic constants of these HEAs is also compared against available experimental data. The results are in good agreement, validating the robustness and accuracy of the computational methods. Multiple phenomenological ductility metrics were also computed and analyzed against the analytical model. Some metrics, mainly the surrogate D parameter, show a good correlation with the Rice model. The potential of these empirical metrics to serve as surrogate screening models for optimizing the compositional space is also discussed.

Anisotropy↗
DOE OSTI2021

Implementation of higher-order velocity mapping between marker particles and grid in the particle-in-cell code XGC

The global total-f gyrokinetic particle-in-cell code XGC, used to study transport in magnetic fusion plasmas or to couple with a core gyrokinetic code while functioning as an edge gyrokinetic code, implements a five-dimensional continuum grid to perform the dissipative operations, such as plasma collisions, or to exchange the particle distribution function information with a core code. To transfer the distribution function between marker particles and a rectangular two-dimensional velocity-space grid, XGC employs a bilinear mapping. The conservation of particle density and momentum is accurate enough in this bilinear operation, but the error in the particle energy conservation can become undesirably large and cause non-negligible numerical heating in a steep edge pedestal. In the present work we update XGC to use a novel mapping technique, based on the calculation of a pseudo-inverse, to exactly preserve moments up to the order of the discretization space. Here we describe the details of the implementation and we demonstrate the reduced interpolation error for a tokamak test plasma using first- and second-order elements with the pseudo-inverse method and comparing with the bilinear mapping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Collisional broadening of nonlinear resonant wave–particle interactions

A general procedure for understanding plasma behaviour when resonant wave–particle interactions are the sole destabilizing and transport mechanism or only heating and/or current drive source is highlighted without recourse to involved numerical or analytical treatments. These phenomena are characterized by transport that appears to be collisionless even though collisions play a central role in narrow collisional boundary layers. The order of magnitude estimates, which include nonlinear effects, are shown to provide expressions in agreement with the principal results of recent toroidal Alfvén eigenmode (TAE), toroidal magnetic field ripple, and heating and current drive treatments. More importantly, the retention of nonlinearities leads to new estimates of the alpha particle energy diffusivity at saturation for TAE modes, and the ripple threshold at which superbanana plateau evaluations of alpha particle transport are modified by nonlinear radial drift effects. In addition, the estimates indicate when quasilinear descriptions for heating and current drive will begin to fail. The phenomenological procedure demonstrates that in magnetic fusion relevant plasmas, narrow collisional boundary layers must be retained for resonant wave–particle interactions as they enhance the role of collisions, and make stochastic particle motion unlikely to be more important than other nonlinear processes.

Physics↗
DOE OSTI2022

Data-driven model for divertor plasma detachment prediction

We present a fast and accurate data-driven surrogate model for divertor plasma detachment prediction leveraging the latent feature space concept in machine learning research. Our approach involves constructing and training two neural networks: an autoencoder that finds a proper latent space representation (LSR) of plasma state by compressing the multi-modal diagnostic measurements and a forward model using multi-layer perception (MLP) that projects a set of plasma control parameters to its corresponding LSR. By combining the forward model and the decoder network from autoencoder, this new data-driven surrogate model is able to predict a consistent set of diagnostic measurements based on a few plasma control parameters. In order to ensure that the crucial detachment physics is correctly captured, highly efficient 1D UEDGE model is used to generate training and validation data in this study. The benchmark between the data-driven surrogate model and UEDGE simulations shows that our surrogate model is capable of providing accurate detachment prediction (usually within a few per cent relative error margin) but with at least four orders of magnitude speed-up, indicating that performance-wise, it has the potential to facilitate integrated tokamak design and plasma control. Comparing with the widely used two-point model and/or two-point model formatting, the new data-driven model features additional detachment front prediction and can be easily extended to incorporate richer physics. This study demonstrates that the complicated divertor and scrape-off-layer plasma state has a low-dimensional representation in latent space. Understanding plasma dynamics in latent space and utilising this knowledge could open a new path for plasma control in magnetic fusion energy research.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

A physics-constrained deep learning surrogate model of the runaway electron avalanche growth rate

A surrogate model of the runaway electron avalanche growth rate in a magnetic fusion plasma is developed. This is accomplished by employing a physics-informed neural network (PINN) to learn the parametric solution of the adjoint to the relativistic Fokker–Planck equation. The resulting PINN is able to evaluate the runaway probability function across a broad range of parameters in the absence of any synthetic or experimental data. This surrogate of the adjoint relativistic Fokker–Planck equation is then used to infer the avalanche growth rate as a function of the electric field, synchrotron radiation and effective charge. Predictions of the avalanche PINN are compared against first principle calculations of the avalanche growth rate with excellent agreement observed across a broad range of parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A closer look at turbulence spreading: How bistability admits intermittent, propagating turbulence fronts

In magnetic fusion plasmas, mounting evidence indicates the possibility of sustained turbulence below the linear stability threshold or more generally global turbulence bistability. The usual reduced models for turbulence spreading are unistable/supercritical and incompatible with this result. The older models further cannot realistically support fronts connecting laminar and turbulent domains. In this work, a minimal model for “subcritical” turbulence spreading is introduced and analyzed. The model may be viewed as phenomenological or derived directly by considering the effect of profile corrugations in an E×B staircase. The model, which is related to the FitzHugh–Nagumo system, supports the robust coexistence of multiple turbulence levels via bistability. We show that this model predicts stronger penetration of turbulence into a linearly stable region as well as the formation of intermittent turbulence fronts that resemble avalanches. We derive the critical size that a localized slug of turbulence must exceed in order to spread. Lastly, we make a prediction of global hysteretic behavior associated with the bistability, which should be testable via experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗