Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fusion plasma”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Containers for Massive Ensemble of I/O Bound Hierarchical Coupled Simulations

We present our experience using containers to scale up a massive ensemble of coupled I/O bound workloads on the NERSC Cori supercomputer. We describe the design of a hierarchical simulation structure using the Integrated Plasma Simulator (IPS) that enables the flexible execution of coupled simulations at the system, node, and core level using the same coupling abstraction and API. The hierarchical design allows for the node-level execution to be efficiently executed using containers while not impacting the structure of the simulation at the system level. We demonstrate the viability of the approach by presenting experimental results from applications in coupled fusion plasma simulations that illustrate the performance impact of using containers to deploy the node-level workloads, in conjunction with the user mountable XFS file systems to ameliorate the load on the Lustre parallel file system. We also present results from production runs showing the ability of the ensemble simulations to scale to hundreds of Cori Haswell nodes, with little or no overhead.

Elwasif, Wael↗
DOE OSTI2021

The five line-of-sight neutron time-of-flight (nToF) suite on the National Ignition Facility (NIF)

Measurement of the neutron spectrum from inertial confinement fusion implosions is one of the primary diagnostics of implosion performance. Analysis of the spectrum gives access to quantities such as neutron yield, hot-spot velocity, apparent ion temperature, and compressed fuel ρr through measurement of the down-scatter ratio. On the National Ignition Facility, the neutron time-of-flight suite has been upgraded to include five independent, collimated lines of sight, each comprising a high dynamic range bibenzyl/diphenylacetylene-stilbene scintillator [R. Hatarik et al., Plasma Fusion Res. 9, 4404104 (2014)] and high-speed fused silica Cherenkov detectors [A. S. Moore et al., Rev. Sci. Instrum. 89, 10I120 (2018)].

Moore, A. S. (ORCID:0000000221258882)↗
DOE OSTI2022

In situ calibration of charged particle spectrometers on the OMEGA Laser Facility using 241 Am and 226 Ra sources

Charged particle spectrometry is a critical diagnostic to study inertial-confinement-fusion plasmas and high energy density plasmas. The OMEGA Laser Facility has two fixed magnetic charged particle spectrometers (CPSs) to measure MeV-ions. In situ calibration of these spectrometers was carried out using 241 Am and 226 Ra alpha emitters. The alpha emission spectrum from the sources was measured independently using surface-barrier detectors (SBDs). The energy dispersion and broadening of the CPS systems were determined by comparing the CPS measured alpha spectrum to that of the SBD. The calibration method significantly constrains the energy dispersion, which was previously obtained through the measurement of charged particle fusion products. Overall, a small shift of 100 keV was observed between previous and the calibration done in this work.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1985

Nondispersive X-ray spectroscopy and imaging of plasmas using a charge-coupled device

A virtual-phase charge-coupled device (CCD) was used to obtain pinhole images and X-ray spectra of laser-produced, solid target plasmas. With the CCD used in the single-photon counting mode, the spectrum in the energy range 2-10 keV was obtained without a dispersive element. Typical spectra reveal two distinct temperatures: a cold component of approximately 200 eV and a hot component of approximately 5 keV. Also, multiline spectra comprising characteristic line emission (K alpha, K beta) from a multilayer target bombarded by beta-rays were recorded using a three-phase CCD. The results demonstrate the potential of CCDs as imaging spectrometers with application in space, laboratory, and fusion-plasma research.

Marsh, K.↗
DOE OSTI2024

First measurement in a magnetic confinement fusion experiment of the H 3 + H 3 → He 5 + n intermediate two-body resonant reaction

We report on the first experimental measurements made at a magnetic confinement fusion device of the tritium(T)-tritium(T) reaction T + T → He 4 + 2 n indicating the presence of the intermediate two-body resonant reaction T + T → He 5 + n . During the second deuterium-tritium campaign (DTE2) at the Joint European Torus, measurements of fusion plasmas with high tritium concentrations, n T / ( n T + n D ) ≈ 0.99 , heated with tritium neutral beam injection, were performed using the neutron time-of-flight (TOF) spectrometer TOFOR. We detect a peak in the neutron emission TOF spectrum consistent with the two-body resonant reaction. The TT neutron emission energy spectrum is modeled using an R -matrix framework where the distributions of the most likely model parameters given our experimental TOF data are determined utilizing a Markov chain Monte Carlo approach. We compare our best estimate of the T + T neutron emission energy spectrum with results obtained at inertial confinement fusion experiments at the OMEGA facility and find a spectral shape that is consistent with the energy dependency in the neutron spectrum observed at OMEGA. Published by the American Physical Society 2024

Physics↗
DOE OSTI2023

Stationary Bragg reflection of laser light in inhomogeneous absorbing plasmas inside inertial confinement fusion Hohlraums

Laser-produced plasma in inertial confinement fusion (ICF) Hohlraums are marked with density non-uniformity whose length scale can go down to micrometers. This scale is of the order of the laser wavelength. The WKB approximation, which is classically used in radiation-hydrodynamic codes to compute the laser trajectory, cannot correctly take into account such small-scale inhomogeneity of the plasma. Going beyond this approximation, we predict a novel mechanism for the laser reflection. Here we show that an electromagnetic plane wave with wave number k resonates with the k B =2 k Fourier component of a multimode perturbation of the background density and generates a reflected wave. It is the first time that this reflection is considered for stationary inhomogeneous ICF plasmas, and the energy absorption is taken into account. This mechanism, which is a form of Bragg reflection, can occur away from the critical surface and generate a drift of the location of the laser absorption. Furthermore, this absorption will be periodically modulated with a k B wave number. The stationary Bragg reflection can explain ongoing discrepancies between experimental and numerical data about laser trajectory and absorption in ICF Hohlraums.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A Fast Algebraic Multigrid Solver and Accurate Discretization for Highly Anisotropic Heat Flux I: Open Field Lines

We present a novel solver technique for the anisotropic heat flux equation, aimed at the high level of anisotropy seen in magnetic confinement fusion plasmas. Such problems pose two major challenges: (i) discretization accuracy and (ii) efficient implicit linear solvers. We simultaneously address each of these challenges by constructing a new finite element discretization with excellent accuracy properties, tailored to a novel solver approach based on algebraic multigrid (AMG) methods designed for advective operators. We pose the problem in a mixed formulation, introducing the directional temperature gradient as an auxiliary variable. The temperature and auxiliary fields are discretized in a scalar discontinuous Galerkin space with upwinding principles used for discretizations of advection. We demonstrate the proposed discretization’s superior accuracy over other discretizations of anisotropic heat flux, achieving error 1000x smaller for anisotropy ratio of 10 9 , for closed field lines. The block matrix system is reordered and solved in an approach where the two advection operators are inverted using AMG solvers based on approximate ideal restriction, which is particularly efficient for upwind discontinuous Galerkin discretizations of advection. To ensure that the advection operators are nonsingular, in this paper we restrict ourselves to considering open (acyclic) magnetic field lines for the linear solvers. We demonstrate fast convergence of the proposed iterative solver in highly anisotropic regimes where other diffusion-based AMG methods fail.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Materials Genomics Search for Possible Helium‐Absorbing Nano‐Phases in Fusion Structural Materials

Abstract Civilian fusion demands structural materials that can withstand the harsh environments imposed inside fusion plasma reactors. The structural materials often transmute under 14.1 MeV fast neutrons, producing helium (He), which embrittles the grain boundary (GB) network. Here, it is shown that neutron‐friendly and mechanically strong nano‐phases with atomic‐scale free volume can have low He‐embedding energy and >10 at.% He‐absorbing capacity, and can be especially advantageous for soaking up He on top of resisting radiation damage and creep, provided they have thermodynamic compatibility with the matrix phase, satisfactory equilibrium wetting angle, as well as a high enough melting point. The preliminary experimental demonstration proves that is a good ab initio predictor of He shielding potency in nano‐heterophase materials, and thus, is used as a key feature for computational screening. In this context, a list of viable compounds expected to be good He‐absorbing nano‐phases is presented, taking into account , the neutron absorption and activation cross‐sections, the elastic moduli, melting temperature, the thermodynamic compatibility, and the equilbrium wetting angle of the nano‐phases with the Fe matrix as an example.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Energies and transition parameters of fusion interest in Cr-like ions between Hf XLIX and Au LVI

Using fully relativistic multiconfiguration Dirac–Fock (MCDF) wavefunctions in an active space approximation, we perform systematic calculations of the energies, wavelengths, and lifetimes for the lowest 143 levels belonging to the 3p{sup 6}3d{sup 6} and 3p{sup 5}3d{sup 7} configurations of Cr-like ions between Hf XLIX and Au LVI. Additionally, line strengths, oscillator strengths, and radiative rates for the electric dipole (E1), electric quadrupole (E2), magnetic dipole (M1), and magnetic quadrupole (M2) transitions among these levels are reported. In the calculations, the contributions from Breit interaction (BI) and quantum electrodynamics (QED) corrections to the wave functions are taken into account (the addition of the BI and QED in CI affects both the wave functions and the level energies). For some excited states, there are no other theoretical and/or experimental data to compare with to illustrate the quality of our results/assess the accuracy. A similar calculation has been performed with the Flexible Atomic Code (FAC), using W LI as an example. The results obtained in this work are compared with the data compiled in the Atomic Spectra Database (ASD) of the National Institute of Standards and Technology (NIST) and other available theoretical/experimental information. Satisfactory agreement is found between them. The energies agree within 0.8% with the experimental results, and the M1 transition rates agree within 5% with other theoretical results. The present results are more extensive than those from previous calculations and considerably expand the existing data sets for Cr-like ions. This study should be beneficial for the analysis of experimental spectra and modeling of fusion plasmas.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Electron capture, ionization and excitation cross sections for keV collisions between fully stripped ions and atomic hydrogen in ground and excited states

Using a semiclassical close-coupling approach, we have calculated electron capture, excitation and ionization cross sections for collisions of fully stripped hydrogen, helium and lithium ions with atomic hydrogen in the ground state and in all excited states up to n=3. The cross sections for collision energies between 1 and 100 keV/u are given in table form. Furthermore, we provide estimates of the accuracy of the cross sections. The set of data presented in this work represent the first complete and consistent quantum study of these collision systems and will find use in the modeling and diagnosis of thermonuclear fusion plasma reactors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Multi-species collisions for delta-f gyrokinetic simulations: Implementation and verification with GENE

we report a multi-species linearized collision operator based on the model developed by Sugama et al. has been implemented in the nonlinear gyrokinetic code, GENE. Such a model conserves particles, momentum, and energy to machine precision, and is shown to have negative definite free energy dissipation characteristics, satisfying Boltzmann’s H-theorem, including for realistic mass ratio. Finite Larmor Radius (FLR) effects have also been implemented into the local version of the code. For the global version of the code, the collision operator has been developed to allow for block-structured velocity space grids, allowing for computationally tractable collisional global simulations. The validity of the collision operator has been demonstrated by relaxation and conservation tests, as well as appropriate benchmarks. The newly implemented operator shall be used in future simulations to study magnetically confined fusion plasma turbulence and transport in more extreme regions with higher collisionality.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

A hybrid nodal-staggered pseudo-spectral electromagnetic particle-in-cell method with finite-order centering

Electromagnetic particle-in-cell (PIC) codes are widely used to perform computer simulations of a variety of physical systems, including fusion plasmas, astrophysical plasmas, plasma wakefield particle accelerators, and secondary photon sources driven by ultra-intense lasers. In a PIC code, Maxwell's equations are solved on a grid with a numerical method of choice. This article focuses on pseudo-spectral analytical time-domain (PSATD) algorithms and presents a novel hybrid PSATD PIC scheme that combines the respective advantages of standard nodal and staggered methods. The novelty of the hybrid scheme consists in using finite-order centering of grid quantities between nodal and staggered grids, in order to combine the solution of Maxwell's equations on a staggered grid with the deposition of charges and currents and the gathering of electromagnetic forces on a nodal grid. The correctness and performance of the novel hybrid scheme are assessed by means of numerical tests that employ different classes of PSATD equations in a variety of physical scenarios, ranging from the modeling of electron-positron pair creation in vacuum to the simulation of laser-driven and particle beam-driven plasma wakefield acceleration. It is shown that the novel hybrid scheme offers significant numerical and computational advantages, compared to purely nodal or staggered methods, for all the test cases presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Generative AI models for learning flow maps of stochastic dynamical systems in bounded domains

Simulating stochastic differential equations (SDEs) in bounded domains, presents significant computational challenges due to particle exit phenomena, which requires accurate modeling of interior stochastic dynamics and boundary interactions. Despite the success of machine learning-based methods in learning SDEs, existing learning methods are not applicable to SDEs in bounded domains because they cannot accurately capture the particle exit dynamics. We present a unified hybrid data-driven approach that combines a conditional diffusion model with an exit prediction neural network to capture both interior stochastic dynamics and boundary exit phenomena. Our ML model consists of two major components: a neural network that learns exit probabilities using binary cross-entropy loss with rigorous convergence guarantees, and a training-free diffusion model that generates state transitions for non-exiting particles using closed-form score functions. The two components are integrated through a probabilistic sampling algorithm that determines particle exit at each time step and generates appropriate state transitions. Here, the performance of the proposed approach is demonstrated via three test cases: a one-dimensional simplified problem for theoretical verification, a two-dimensional advection-diffusion problem in a bounded domain, and a three-dimensional problem of interest to magnetically confined fusion plasmas.

Bounded domains↗
DOE OSTI2026

Achieving high tensile strength and ductility in refractory alloys by tuning electronic structure

The energy efficiency of heat engines (gas and steam turbines) for electricity production and propulsion is determined by the Carnot cycle and scales with operating temperature. Commercial nickel- and cobalt-based superalloys melt near 1,500 °C and rapidly lose mechanical strength beyond 1,000 °C. Refractory metals melt well above 2,000 °C but have inherent manufacturability challenges that are barriers to adoption, such as high ductile-to-brittle transition temperatures. Using density functional theory-guided design, we demonstrate tailored local lattice distortions that promote phase-stable, non-equiatomic refractory concentrated solid solutions with both high ductility and strength. Here, we exemplify this for single-phase, body-centred cubic Nb 4 Ta 4 V 3 Ti that exhibits castability, excellent room-temperature tensile yield strength (∼1 GPa) and ductility (approaching 20% uniform strain), and exceptional high-temperature tensile strength (500 MPa at 1,000 °C). These findings illustrate a path for designing materials that hold great potential for advancing next-generation technologies such as Generation IV fission reactors, first-generation fusion-plasma reactors, and more efficient gas turbines for electricity generation and propulsion.

DFT↗
DOE OSTI2021

A novel photomultiplier tube neutron time-of-flight detector

A traditional neutron time-of-flight (nTOF) detector used in inertial confinement fusion consists of a scintillator coupled with a photomultiplier tube (PMT). The instrument response function (IRF) of such a detector is dominated by the scintillator-light decay. At DT neutron yields larger than 1013 a novel detector consisting of a microchannel-plate photomultiplier tube in a housing without a scintillator (PMT nTOF) can be used to measure DT yield, ion temperature, and neutron velocity. Most of the neutron signal in PMT nTOF is produced from neutron interaction with a PMT window. The direct interaction of neutrons with MCP has negligible contribution. The elimination of the scintillator removes the scintillator decay from the instrument response function and makes the IRF of the PMT nTOF faster, which makes the ion temperature and neutron velocity measurements more accurate. Three PMT nTOFs were deployed on the OMEGA Laser System for the first time to diagnose inertial confinement fusion plasma. Here, the design details, characteristics, and calibration results of these detectors in DT implosions on OMEGA are presented. Recommendations on the use of different PMTs for specific applications are provided.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

New, improved analysis of correlation ECE data to accurately determine turbulent electron temperature spectra and magnitudes (invited)

Turbulent electron temperature fluctuation measurement using a correlation electron cyclotron emission (CECE) radiometer has become an important diagnostic for studying energy transport in fusion plasmas, and its use is widespread in tokamaks (DIII-D, ASDEX Upgrade, Alcator C-Mod, Tore Supra, EAST, TCV, HL-2A, etc.). The CECE diagnostic typically performs correlation analysis between two closely spaced (within the turbulent correlation length) ECE channels that are dominated by uncorrelated thermal noise emission. This allows electron temperature fluctuations embedded in the thermal noise to be revealed and fluctuation level and spectra determined. We have demonstrated a new, improved CECE coherency-based analysis for calculating the temperature fluctuation frequency spectrum and level, which has been verified both numerically through the simulation of synthetic ECE radiometer data and through analysis of experimental data from the CECE system on DIII-D. The new formulation places coherency-based analysis on a firm foundational footing and corrects some currently published methodologies. This new method accurately accounts for bias error in the coherence function and correctly calculates noise levels for a fixed data record length. It provides excellent accuracy in determining temperature fluctuation level (e.g., <10% error) even for a small realization number in the ensemble average. The method also has a smaller uncertainty (i.e., error bar) in the power spectrum when compared to the more standard cross-power method when evaluated at low coherency. Direct calculation of system noise level using correlation between randomized intermediate frequency signals is recommended.

Wang, G. (ORCID:0000000225739827)↗
DOE OSTI2025

Axisymmetric gyrokinetic simulation of ASDEX-Upgrade scrape-off layer using a conservative implicit BGK collision operator

Collisions play an important role in turbulence and transport of fusion plasmas. For kinetic simulations, as the collisionality increases in the domain of interest, the size of the time step to resolve the collisional physics can become overly restrictive in an explicit time integration scheme, leading to high computational cost. With the aim of overcoming such restriction, we have implemented an implicit Bhatnagar–Gross–Krook (BGK) collision operator for use in the discontinuous Galerkin full-f gyrokinetic solver within the Gkeyll framework, which, when combined with Gkeyll's traditional explicit time integrator for collisionless advection, can significantly increase the time step in gyrokinetic simulations of highly collisional regimes. To ensure conservation of density, momentum, and energy, we utilize an iterative scheme to correct the discretized approximation to the equilibrium Maxwellian distribution to which the BGK collision operator relaxes. We have further generalized the BGK infrastructure, both the implicit scheme and the correction routine, to handle cross-species collisions. This improved implicit and conservative BGK operator is benchmarked against the more accurate but more computationally expensive Lenard–Bernstein–Dougherty (LBD) operator, which has been utilized in prior studies with Gkeyll. The implicit BGK operator enables 2D axisymmetric simulations of the ASDEX-Upgrade scrape-off layer to run 56 times faster to completion than the simulations with the LBD operator, because the BGK operator is more robust and converges at a lower resolution than is required by the LBD operator. Additionally, in this more collisional limit, we demonstrate that the results of our simulations utilizing the implicit BGK operator agreed well with simulations utilizing the more computationally expensive LBD operator.

Gyrokinetic simulations↗
DOE OSTI2026

Inferring performance metrics for laser direct drive experiments on OMEGA

Quantifying performance improvements on the OMEGA laser facility requires robust inference of established no-alpha performance metrics, which requires, at minimum, a model to infer the shocked fuel mass and pressure of the confined fusion plasma. In this work, we describe the methodology used to infer performance metrics on OMEGA and present the current state-of-the art model used to infer these metrics from OMEGA experiments. In particular, since neutron images of cryogenic implosions are not available on OMEGA at present, we present how x-ray sizes are determined on OMEGA using a Gaussian Process regression model and how the neutron production region's size is inferred from them. As a result, we end by benchmarking the model using synthetic data and 1-D LILAC simulations and test its experimental self-consistency across available x-ray diagnostic channels.

Gopalaswamy, V. [Laboratory for Laser Energetics, ↗