Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plasma accelerators”

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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Physics of runaway electrons with shattered pellet injection at JET

Runaway electrons (REs) created during tokamak disruptions pose a threat to the reliable operation of future larger machines. Experiments using shattered pellet injection (SPI) have been carried out at the JET tokamak to investigate ways to prevent their generation or suppress them if avoidance is not sufficient. Avoidance is possible if the SPI contains a sufficiently low fraction of high-Z material, or if it is fired early in advance of a disruption prone to runaway generation. These results are consistent with previous similar findings obtained with Massive Gas Injection. Suppression of an already accelerated beam is not efficient using High-Z material, but deuterium leads to harmless terminations without heat loads. This effect is due to the combination of a large magnetohydrodynamic instability scattering REs on a large area and the absence of runaway regeneration during the subsequent current collapse thanks to the flushing of high-Z impurities from the runaway companion plasma. This effect also works in situations where the runaway beam moves upwards and undergoes scraping-off on the wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A Performance Portable, Fully Implicit Landau Collision Operator with Batched Linear Solvers

Modern accelerators use hierarchical parallel programming models that enable massive multithreading within a processing element (PE), with multiple PEs per device driven by traditional processes. Batching is a technique for exposing PE-level parallelism in algorithms that have traditionally run on MPI processes or multiple threads within a single process. Opportunities for batching arise in, for example, kinetic discretizations of magnetized plasmas where collisions are advanced in velocity space at each spatial point independently. This paper builds on previous work on a high-performance, fully nonlinear, Landau collision operator by batching the linear solver, as well as batching the spatial point problems and adding new support for multiple grids for multiscale, multispecies problems. An anisotropic relaxation verification test that agrees well with previously published results and analytical models is presented. The performance results from NVIDIA A100 and AMD MI250X nodes are presented with hardware utilization analysis for each architecture. Finally, the entire implicit Landau operator time advance is implemented in Kokkos for performance portability, running entirely on the device and is available in the PETSc numerical library.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Studying the Richtmyer–Meshkov instability in convergent geometry under high energy density conditions using the Decel platform

The “Decel” platform at Sandia National Laboratories investigates the Richtmyer–Meshkov instability (RMI) in converging geometry under high energy density conditions [Knapp et al., Phys. Plasmas 27, 092707 (2020)]. In Decel, the Z machine magnetically implodes a cylindrical beryllium liner filled with liquid deuterium, launching a converging shock toward an on-axis beryllium rod machined with sinusoidal perturbations. The passage of the shock deposits vorticity along the Be/D2 interface, causing the perturbations to grow. Here, we present platform improvements along with recent experimental results. To improve the stability of the imploding liner to the magneto Rayleigh–Taylor instability, we modified its acceleration history by shortening the Z electrical current pulse. Next, we introduce a “split rod” configuration that allows two axial modes to be fielded simultaneously in different axial locations along the rod, doubling our data per experiment. We then demonstrate that asymmetric slots in the return current structure modify the magnetic drive pressure on the surface of the liner, advancing the evolution on one side of the rod by multiple ns compared to its 180° counterpart. This effectively enables two snapshots of the instability at different stages of evolution per radiograph with small deviations of the cross-sectional profile of the rod from the circular. Using this platform, we acquired RMI data at 272 and 157 μm wavelengths during the single shock stage. Finally, we demonstrate the utility of these data for benchmarking simulations by comparing calculations using ALEGRA MHD and RageRunner.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Explainable tokamak-agnostic forecasting of fusion plasma instability via megahertz turbulent fluctuations

Scientific applications of artificial intelligence (AI) often remain limited by device-specific training and unexplained “black-box” approaches, creating fundamental barriers to cross-system generalization. This challenge is critical for nuclear fusion, where future reactors will have limited operational data for AI training. Here, we demonstrate that our neural network, trained solely on megahertz-scale turbulence measurements from one machine (DIII-D), forecasts Type-I edge localized mode (ELM) onsets in a different tokamak (KSTAR) through zero-shot weight transfer following physics-consistent preprocessing without device-specific retraining. Through an explainable AI framework combining gradient-weighted class activation mapping with physics validation, we reveal that our network can internalize physics relationships governing the ELM instabilities rather than memorizing device-specific patterns. The network perceives spatiotemporal features that correlate consistently with independently calculated instability growth rates, magnetohydrodynamic stability limits, and pedestal structure dynamics. Statistical analyses of dimensionally-reduced saliency features reveal the identical triangular features between the saliency representations, instability growth rates, and prediction probability across tokamaks, providing evidence that our forecasting system can show tokamak-agnostic generalization. This work contributes to a foundation for explainable scientific AI systems, where cross-system developments are essential for transcending traditional domain-specific constraints.

AI↗
DOE OSTI2022

Modeling ablator grain structure impacts in ICF implosions

High-density carbon is a leading ablator material for inertial confinement fusion (ICF). This and some other ablator materials have grain structure which is believed to introduce very small-scale (~nm) density inhomogeneity. In principle, such inhomogeneity can affect key ICF metrics like fuel compression and yield, by, for example, acting as a seed for instabilities and inducing mix between ablator and fuel. However, assessments of such effects are uncertain due to the difficulty of modeling this small-scale structure in ICF simulations, typically requiring reduced-resolution modeling that scales these features. We present a grain model and show both the impact of de-resolving grains and the complex mixing dynamics such structures can induce. We find that different methods for de-resolving grains can yield both different total deposition of kinetic energy perturbations and different fuel–ablator mixing. We then show a simple-to-implement approach for approximately conserving the deposition of perturbed kinetic energy and demonstrate that, for the present grain model and test cases, this approach yields a reasonably matched time history of mix width between less and more resolved grain models. The simulations here also demonstrate the complex interaction history between grain-induced mixing and instability around the fuel–ablator interface, showing, for example, that the grain-induced perturbations typically trigger instability of conduction-driven density gradients in the DT fuel, enhancing mix penetration early in the acceleration of the shell. Furthermore simulating both microscale and nanoscale grains, we find initial evidence for larger mixing in the microscale case of the present model, despite smaller deposited kinetic energy perturbation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A process optimization framework for laser direct energy deposition: Densification, microstructure, and mechanical properties of an Fe-Cr alloy

Laser Direct Energy Deposition (DED) is a metal additive manufacturing technique with the ability to fabricate large and complex parts through deposition of metal powders. However, achieving high-density parts and targeted build heights using DED can be challenging due to the large number of highly sensitive process variables. This work proposes a robust fabrication parameter optimization framework to generate process maps for primary parameters in DED, including laser power, scan speed, mass flow rate, hatch spacing, and layer height. Simple single-track experiments were utilized to map out the parameter space, and a combination of geometric criteria for hatch spacing and layer height were proposed to determine parameter sets that achieve both targeted build heights and mitigate porosity formation. Using this framework, specimens with >99 % density and consistent mechanical properties were successfully fabricated over a wide range of process parameters for an Fe-9wt.%Cr (Fe9Cr) alloy, a surrogate for radiation damage-resistant reduced activation ferritic/martensitic (RAFM) steels. Processing these materials using DED is of particular interest in the development of plasma facing components for nuclear fusion applications. The microstructure and mechanical properties of as-printed Fe9Cr were characterized using optical and electron microscopy, X-ray diffraction, and uniaxial tensile tests. As-printed Fe9Cr displayed ~25 % elongation and ultimate tensile strengths of up to 475 MPa which is comparable to similar wrought alloys. Finally, the proposed framework will allow for accelerated DED parameter optimization for novel alloy systems, as well as open the possibility for local microstructure control while simultaneously mitigating defect formation.

42 ENGINEERING↗
DOE OSTI2024

Focused Ion Beam Tomography of Alloy 617 Corroded in Molten Chloride Salt

Materials qualification of reactor structural materials is a critical step in rapid implementation of advanced nuclear reactor technologies, particularly to assess the corrosion performance in these designs. Accelerated qualification of reactor structural materials requires incorporating powerful computational toolsets, such as phase field modelling in the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework, to predict the evolution of structural materials due to corrosion. Accordingly, computational toolsets will require experimental data generated at appropriate length scales to validate accuracy. Focused ion beam (FIB) provides a high degree of control over manipulation of materials for analytical purposes, including capturing data on the evolution in the microstructure and elemental composition of materials at the mesoscale, an appropriate length scale for phase field modelling of intergranular diffusion phenomena using the MOOSE framework. For instance, the FEI Helios G4 UX dual beam plasma FIB microscope at the Irradiated Materials Characterization Laboratory (IMCL) is capable of backscatter diffraction (EBSD) and energy-dispersive x-ray spectroscopy (EDS) documenting the evolution in the microstructure and elemental composition, respectively. The Helios can perform EDS and EBSD three-dimensionally (3D) using tomography, which is then combined using different software packages to visualize 3D volumes correlating elemental composition to microstructural data. The purpose of this investigation was to develop a streamlined characterization and data processing workflow for 3D tomography studies on the FEI Helios G4 plasma FIB. The investigation is segmented into three parts: 1) Optimizing the data collection workflow, 2) identifying appropriate data processing and visualization software (i.e. DREAM.3D, MIPAR, and VGStudioMax), and 3) establishing an infrastructure for public release. The optimization of the data collection workflow is in collaboration with members of the U220 department to setup formal training on the tomography operation of the G4, through ThermoFisher Scientific, and exploring DREAM.3D, MIPAR, and VGStudioMax data processing/visualization software packages. VGStudioMax currently demonstrates the most promise for future use. Optimization of the data collection and processing workflow is still ongoing. A collaboration with INL High Performance Computing (HPC) established an open-source license for expediting the public release of FIB tomography datasets through HPC. FIB tomography data generated by the G4 will provide comprehensive data for validating 3D phase field mesoscale modelling tools within the MOOSE framework for accelerated qualification of reactor structural materials.

Copeland-Johnson, Trishelle↗
DOE OSTI2021

How Hall electric fields intrinsically chaotize and heat ions during collisionless magnetic reconnection

The origin of extreme and anomalous ion heating during magnetic reconnection has been a longstanding mystery. It is shown in this paper that, at the core of the process, the stochastic heating mechanism is responsible for this energization. The in-plane Hall electric fields that intrinsically develop during collisionless reconnection not only accelerate ions, but also thermalize them by inducing chaotic motion. The mechanism explains various observed features of ion heating, e.g., its relation to the guide field strength and the ion mass-to-charge ratio.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Plasma self-driven current in tokamaks with magnetic islands

Magnetic island perturbations may cause a reduction in plasma self-driven current that is needed for tokamak operation. A novel effect on tokamak self-driven current revealed by global gyrokinetic simulations is due to magnetic-island-induced 3D electric potential structures, which have the same dominant mode numbers as that of the magnetic island, whereas centered at both the inner and outer edge of the island. The non-resonant potential islands are shown to drive a current through an efficient nonlinear parallel acceleration of electrons. In large aspect ratio (large-A) tokamak devices, this new effect can result in a significant global reduction of the electron bootstrap current when the island size is sufficiently large, in addition to the local current loss across the island region due to the pressure profile flattening. It is shown that there exists a critical magnetic island width for large-A tokamaks beyond which the electron bootstrap current loss is global and increases rapidly with the island size. As such, this process may introduce a size limit for tolerable magnetic islands in large-A tokamak devices in the context of steady state operation. On the other hand, the current loss caused by magnetic islands in low-A tokamaks such as spherical tokamak (ST) NSTX/U is minor. The reduction of the axisymmetric current by magnetic islands scales with the square of island width. However, the loss of the current is mainly local to the island region, and the pace of current loss as the island size increases is substantially slower compared to large-A tokamaks. In particular, the bootstrap current reduction in STs is even smaller in the reactor-relevant high-β p regime where neoclassical tearing modes are more likely to develop.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

3D ion gyro-orbit heat load predictions for NSTX-U

Abstract High power tokamaks operate with divertor heat loads capable of destroying the plasma facing components (PFCs). High fidelity heat load predictions are necessary to ascertain the PFC state for design and during operation. Typical heat flux calculations are 2D, time invariant, and assume that power flows directly along the magnetic field lines (the optical approximation). These assumptions neglect the complex 3D geometries employed to protect the PFCs, the time varying nature of the plasma and PFC thermal state, and the helical trajectories of ions with finite Larmor radii (the gyro-orbit approximation). An integrated software framework, the heat flux engineering analysis toolkit (HEAT), was developed to generate time varying optical heat loads applied to real engineering computer aided design (CAD) (Looby et al 2022 Fusion Sci. Technol. 78 10–27). Recently, an ion-gyro orbit module has been added to HEAT. This module calculates the helical trajectories of ions as they gyrate about the magnetic field lines using kinetic theory macro-particles to accelerate the calculation. First, the new gyro-orbit module will be presented. Next, a comparison to existing research is performed. Finally, an analysis of the gyro-orbit heat loads for NSTX-U is presented for diverted discharges using the engineering CAD models utilized for PFC fabrication. Including these gyro-orbit effects can enhance the PFC performance by ‘smearing’ out the magnetic shadows associated with the castellated fish-scaled geometry. Simultaneously, the helical trajectories can degrade performance when they load narrow regions on edges and corners with high heat fluxes. Analysis of the trade-offs between these competing effects is included, and regions for further investigation are identified.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Time-resolved ion energy distribution functions during a HiPIMS discharge with cathode voltage reversal

Abstract The effect on the ion energy distribution function (IEDF) of plasma produced during a high-power impulse magnetron sputtering (HiPIMS) discharge as the pulse conditions are varied is reported. Pressure was varied from 0.67–2.00 Pa (5–15 mTorr), positive kick pulses up to 200 V tested with a constant 4 μ s delay between negative and positive cycles. The results demonstrate that the resulting plasma during the positive cathode voltage reversal is the result of expansion through the largely neutral gas species between the end of the magnetic trap of the target and the workpiece. The plasma potential rises on similar time scale with the evolution of a narrow peak in the IEDF close to the applied bias. The peak of the distribution function remains narrow close to the applied bias irrespective of pulse length, and with only slight pressure dependence. One exception discovered is that the IEDF contains a broad high energy tail early in the kick pulse due to acceleration of ions present beyond the trap from the main pulse separate from the ionization front that follows.

Physics↗
DOE OSTI2020

Corrosion of AISI 310 and 316 in Molten Hydroxide Under Steam Electrolysis Conditions

Steam electrolyzers, utilizing molten hydroxide electrolyte, offers potential for improvement in electrochemical efficiency, cost reduction, use of conventional materials of construction and process scale up for large scale hydrogen production. Stainless steels, used for the fabrication of cell components (current collector, gas separator, wet seals, and manifolds) experience accelerated corrosion in the presence of molten hydroxide electrolyte in both oxidizing (anodic) and reducing (cathodic) atmospheres. In this study, the corrosion behavior of AISI 310 and 316 in hydroxide melt has been studied at 600°C for 50 h. Melt immersion tests revealed formation of the porous lithium iron oxide at the melt-oxide interface for both AISI 310 and 316. Inductively coupled plasma optical emission spectrometry (ICP-OES) analysis showed the presence of Cr in the electrolyte melt obtained from samples exposed to oxidizing atmosphere. Thermochemical analysis was performed to validate the experimental results.

corrosion↗
DOE OSTI2021

Evidence of Alfvenic Poynting Flux as the Primary Driver of Auroral Motion During a Geomagnetic Substorm

Geomagnetic substorms are major energy transfer events where energy stored in the Earth's magnetotail is released into the ionosphere. Substorm phenomena, including auroral activities, earthward Poynting flux, magnetic field dipolarization, etc, have been extensively studied. However, the complex interplay among them is not fully understood. In a fortuitous event on June 7, 2013, the twin Van Allen Probes (separated by 0.4 h in local time) observed bursts of earthward Alfvenic Poynting flux in the vicinity of the plasma sheet boundary layer (PSBL). The Poynting flux bursts correlate with enhancements of auroral brightness around the footpoints of both spacecraft. This indicates a temporal and spatial correlation between the auroral brightening and Poynting flux bursts, and that the auroral motion is directly linked to the perpendicular expansion of the Alfven wave. These observations suggest that the Alfvenic Poynting flux is a primary driver for the auroral electron acceleration. Furthermore, around the time of auroral brightening, a dipolarization was seen to propagate more than 4 h in local time during a 20 min period. The azimuthal phase speed of this dipolarization (2 deg/min) is too small to explain the azimuthal motion of the aurora (13.6 deg/min), but the dipolarization could be related to the generation of the Alfvenic Poynting flux through phase mixing at strong density gradients like those in the PSBL.

58 GEOSCIENCES↗
DOE OSTI2023

Density Measurement of NaCl-MgCl 2 -PuCl 3 and NaCl-UCl 3 -PuCl 3 Molten Salt Systems by Neutron Radiographic Dilatometry

A novel neutron radiographic technique was developed and implemented to measure the change in volume as a function of temperature of select plutonium-containing liquid chloride mixtures from melt to 1250 K. These measurements were performed at Los Alamos National Laboratory by members of the molten salt research group under support from Gateway for Accelerated Innovation in Nuclear (GAIN) voucher NE-21-25117: “Density Measurements of Plutonium Bearing Salts via Neutron Beam Dilatometry”. This technical report contains a description of the experimental method and a summary of the results. Initial characterization of the {NaCl, MgCl 2 , UCl 3 , PuCl 3 } salts used in this experiment includes a combination of x-ray diffraction, differential scanning calorimetry, gamma ray spectroscopy, and inductively coupled plasma mass spectrometry. Four samples of NaCl + MgCl 2 + x mol% PuCl 3 and four samples of NaCl + UCl 3 + x mol% PuCl 3 were prepared at the Plutonium Facility, shipped to the Los Alamos Neutron Science Center (LANSCE), and radiographed with the density via neutron radiography (DvNR) apparatus. Density data in the form of meniscus heights were recorded as a function of temperature. The results show a linear dependence of density with temperature and a third order polynomial dependence of density with composition as a function of PuCl 3 addition in the corresponding halide salt mixture. Additions of PuCl 3 result in a significant increase in density of the salt mixture and a slight reduction in the coefficient of volumetric expansion. These results represent a first measurement of the density of plutonium chloride salt mixtures and, in combination with additional measurements, may lead to a generalized description of the density of such salt systems for arbitrary compositions and temperature.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Durability of PGM catalyst MEAs of polymer electrolyte membrane fuel cells for heavy-duty vehicles

Polymer electrolyte membrane fuel cells (PEMFCs) are promising power sources for heavy-duty vehicles (HDVs) owing to cleanliness and efficiency. However, the degradation of membrane electrode assemblies (MEAs) under HDV conditions remains a huge challenge. Here, this work investigated MEA durability under HDV conditions using a US Department of Energy standard accelerated stress test for 180,000 cycles (equivalent to 1 million miles of HDV operation). Effects of catalyst Pt content on MEA durability were examined using homemade 30% Pt/C (H-Pt/C) and commercial 46% Pt/C (C-Pt/C) catalysts. Both MEAs experienced H 2 /air and H 2 /O 2 performance loss over cycles. Analysis with scanning transmission electron microscopy, X-ray diffraction, inductively coupled plasma mass spectrometry, and mercury intrusion porosimetry revealed severe degradation of Pt nanoparticles (NPs), support structures, and the catalyst layer. Two degradation stages for NPs were proposed: Ostwald ripening dominated the initial 60,000 cycles, followed by combined Ostwald ripening and particle migration. Measurements with ion chromatography, high-frequency resistance, and oxygen-diffusion resistance revealed degradation of membrane and ionomer, respectively.

25 ENERGY STORAGE↗
DOE OSTI2025

Quantum physics of stars

Stars are slowly developing objects; the lifetimes of the different burning phases are determined by the strength of nuclear reactions, which in turn are defined by the quantum structure of the associated nuclei at the threshold and the respective reaction mechanisms. Stars, from the nuclear physics perspective, are cold environments where only a few of the key nuclear reactions have been measured at the actual stellar plasma temperatures. This is also the case for more dynamic astrophysical phenomena from the big bang to stellar explosions. Most of the nuclear reaction rates are therefore based on theoretical extrapolations. A number of discrepancies between these predictions and the associated stellar signatures have been observed, and many may be due to low-energy or near-threshold quantum effects. These effects need to be understood in order to reliably model nuclear reaction processes, not only for stars but also for low-temperature plasma environments such as controlled magnetic or inertial confinement fusion systems, which operate in similar temperature regimes. This review summarizes the various theoretical techniques presently used for deriving reaction rates and discusses possible quantum effects that may impact the reaction cross section near the reaction threshold. These resemble enhanced single-particle and cluster structures near threshold and associated interference effects. New experimental techniques such as deep-underground accelerators or the study of transfer reactions to mimic the quantum-mechanical transition strength, the so-called Trojan horse method, provide ways to directly or indirectly probe the reaction features that determine the reaction rates at stellar energies. Furthermore, this is demonstrated on a number of key nuclear reactions for different nucleosynthesis environments. Finally, current inconsistencies between experimental predictions and observations are discussed.

Models & methods for nuclear reactions↗
DOE OSTI2023

Axion-mediated Transport of Fast Radio Bursts Originating in Inner Magnetospheres of Magnetars

Abstract Among magnetar models of fast radio bursts (FRBs), there is ongoing debate about whether the site of coherent radio emission lies within or beyond the light cylinder. We propose a mechanism by which FRBs produced near the magnetar surface are transported out of the magnetosphere by axions, which are hypothetical particles that couple to photons. If the emission site hosts strong accelerating electric fields, a considerable fraction of the FRB energy budget is converted to an axion burst. Once produced, the axion burst free streams out of the magnetosphere due to the rapidly decreasing magnetic field. The burst may escape through either the open or closed magnetosphere while retaining the temporal signature of the original FRB. In the wind region, axions resonantly excite ordinary (O) modes that escape as the plasma density decreases. The radio efficiency of this mechanism satisfies energetics constraints from FRB 121102 for axion−photon coupling strengths that have not been excluded by other astrophysical probes.

Prabhu, Anirudh (ORCID:0000000191157844)↗
DOE OSTI2022

High-Throughput Ion Irradiation and Microstructural Characterization of Multi-Principal-Element Alloys

Multi-principal-element alloys (MPEAs) have garnered interest in the nuclear community for their apparent resistance to microstructural damage under irradiation. While both experimentation and simulation have demonstrated reduced hardening and void swelling in MPEAs, the relationship between composition and irradiation response is still poorly understood since a rich experimental database for MPEAs does not yet exist. In this work, to accelerate the generation of experimental irradiation-response data for MPEAs, high-throughput synthesis, irradiation, and characterization techniques have been employed. Specifically, additive manufacturing has been used to produce arrays of various MPEAs (1-cm2 coupons), which have been subsequently irradiated using 4-MeV Ni2+ ions to a peak damage of 200 dpa at 500 °C. Hardening and void swelling have been assessed using nanoindentation and trenching via plasma focused ion beam (PFIB), respectively. Trends in the experimental data relating hardening and void swelling behavior to composition have been analyzed using different machine learning approaches.

36 MATERIALS SCIENCE↗