Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “space 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 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

CFD Simulation of Helium Flow Loop Test Section

A helium flow loop is being assembled at Oak Ridge National Laboratory to analyze heat transfer enhancement for systems such as blanket and divertor components. To efficiently identify optimum geometries for heat transfer enhancement in these applications, simulation work is performed to optimize test section designs that are built and tested in the helium flow loop that operates at 4 MPa and a mass flow rate of 100 g/s. Different ribbed geometries that examine rib shape, rib height, rib orientation, rib spacing, and three dimensional orientation are modeled and simulated in STAR-CCM+ to compare their ability to remove heat and mitigate pressure drop. Following the simulations, models are selected and manufactured for the helium flow loop tests. Simulations initially focus on a hydrodynamic study to determine the appropriate mesh and physics models and then add a heat flux to analyze the heat transfer abilities of the models. The simulations are run in steady state and use a Reynolds-averaged Navier-Stokes k-ε turbulence model. The helium is modeled as an ideal gas. The simulation explores models of geometries that enhance the heat transfer and decrease pressure drop with an overall goal of increasing fluid collision with the wall. Enhanced geometries are simulated to select appropriate designs for manufacturing, and preliminary experimental results are used to validate the simulations. Furthermore, the factors that are being analyzed in the comparison between the experimental and the simulated results include matching thermocouple temperatures, pressure drop, roughness, and fluid velocity.

42 ENGINEERING↗
DOE OSTI2024

Hot Hydrogen Testing of W-Coated UN Kernels in a Mo30W Matrix

Ceramic uranium mononitride (UN) is being considered as a reactor fuel for nuclear thermal propulsion. To avoid or reduce the dissociation of UN at the high temperatures needed, embedding it in a metallic matrix (cermet) has been proposed. To assess the viability of this concept, hot hydrogen testing of tungsten-coated UN kernels embedded in a Mo-30 wt% W (Mo30W) alloy matrix has been performed at temperatures from 1800°C to 2300°C. Both the isolated kernels and kernels consolidated by spark plasma sintering in the Mo30W matrix were tested. In addition to direct observations and mass loss measurements, the samples were analyzed by X-ray diffraction (XRD) and scanning electron microscopy (SEM)/energy dispersive X-ray spectroscopy (EDS) after each run. The decomposition of UN started at 1800°C despite the coating and matrix, and increased at 2000°C. Uranium seeped through the tungsten grain boundaries of the coating at all temperatures. The consolidated sample expanded irregularly at 2000°C through the formation of voids, and SEM/EDS analysis showed uranium-containing veins in the matrix consisting of U 2 Mo according to the XRD data. The observed pore generation at 2000°C was explained by the formation of water vapor from residual oxides and diffused hydrogen. At 2200°C and above, both the kernels and the consolidated samples melted through the formation of uranium or low–melting point uranium-molybdenum alloys.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

In-medium bottomonium properties from lattice NRQCD calculations with extended meson operators

We calculate the temperature dependence of bottomonium correlators in (2+1)-flavor lattice QCD with the aim to constrain in-medium properties of bottomonia at high temperature. The lattice calculations are performed using HISQ action with physical strange quark mass and light quark masses twenty times smaller than the strange quark mass at two lattice spacings a = 0.0493 fm and 0.0602 fm, and temporal extents N τ = 16 − 30, corresponding to the temperatures T = 133 − 250 MeV. We use a tadpole-improved NRQCD action including spin-dependent v 6 corrections for the heavy quarks and extended meson operators in order to be sensitive to in-medium properties of the bottomonium states of interest. We find that within estimated errors the bottomonium masses do not change compared to their vacuum values for all temperatures under our consideration; however, we find different nonzero widths for the various bottomonium states.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Evaluation of SPARC divertor conditions in H-mode operation using SOLPS-ITER

The predicted divertor conditions for the SPARC tokamak are calculated using SOLPS-ITER for a range of scrape-off-layer (SOL) heat flux widths λ q , input powers, and particle fueling locations. Under H-mode scenario conditions with an upstream separatrix density of 1 x 10 20 m -3 , the most conservative range of λ q extrapolations ( 0.15 mm) results in extremely high unmitigated particle and energy fluxes to the divertor, both under full field (12.2 T) and power (P SOL = 29 MW) conditions, and 2/3 field with P SOL = 10 MW. Increasing the cross-field SOL diffusivities by 2–10× reduces the magnitude of the mitigation challenge, however strategies such as impurity seeding or strike-point-sweeping will likely still be required. A combination of steady-state and time-dependent SOLPS-ITER simulations are used to map out phase space diagrams of upstream and divertor conditions. The simulations include parallel currents but neglect cross-field drifts. At low upstream density the inner and outer divertor conditions are highly asymmetric, with a large temperature difference and significant heat fluxes driven by parallel currents. The solution has sharp bifurcations with a region of hysteresis, depending on whether the initial state is at a low or high density. This behavior is observed even when the fueling location, cross-field diffusivity, and impurity level is changed, although the density window with asymmetry is reduced with increasing diffusivity. The addition of neon impurity seeding reduces the divertor heat fluxes, but also causes a drop in the upstream electron density with fixed particle throughput. This drop can be counteracted by increased main ion throughput, however too much neon results in a back transition into the asymmetric divertor regimes suggesting a need for control of both main ion and impurity seeding levels to achieve a desired divertor state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Automated and highly parallelized Bayesian optimization scheme for direct drive fusion experiments on OMEGA

Finding the optimal implosion design on existing experimental facilities for inertial confinement fusion requires an exhaustive search of the vast design parameter space. This is infeasible both with experiments and with simulations. Consequently, a large fraction of the experimentally realizable design space remains unexplored, and new design schemes are challenging to optimize in a reasonable time frame. On the OMEGA laser facility, predictive machine learning models have been developed to accurately forecast the result of an experiment using only inexpensive simulations and the large dataset of prior experimental data. However, the full design space remains vast enough to be unassailable with simple optimization techniques. Here we develop an automated and optimally parallel Bayesian optimization algorithm that can entirely optimize the target and pulse shape of a direct-drive ICF implosion under a given design paradigm. We use this algorithm to find a markedly improved design for the performance implosions on OMEGA that is predicted to hydroequivalently scale to ignition at 2.15 MJ.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Simultaneous Observation of Ion-scale Wave Packets with Opposite Polarizations and Their Implications on the Generation Region in the Inner Heliosphere

This paper reports a dispersion analysis of two wave packets simultaneously observed near the local proton gyrofrequency by the Parker Solar Probe. The observed wave event exhibits clear two-banded wave packets both propagating along the magnetic field, characterized by left-handed (L-mode) and right-handed (R-mode) polarizations simultaneously. By incorporating the Doppler shift effect into a linear dispersion analysis, we find two possible scenarios that explain these simultaneous opposite polarizations: (1) Two inherently L-mode waves in the plasma frame, propagate parallel and antiparallel to the solar wind velocity, with similar wave frequencies and wave numbers. The polarization of the antiparallel propagating wave reverses as it moves sunward in the plasma frame while still comoving with the solar wind in the stationary frame. This reversal manifests the polarization of the wave as an R-mode in the spacecraft frame. (2) Simultaneous L-mode and R-mode waves propagate parallel to the solar wind velocity, with different wave frequencies and wave numbers. Concurrent proton observations during the wave event reveal a dominant anisotropic ($T$⟂/$T$ ∥ > 1) core distribution with a drifting beam population. Estimation of the linear growth rate for both L-mode and R-mode waves suggests that both scenarios are plausible, indicating that the observation is near the wave-generation region. We explore the potential impact of these simultaneous waves on solar wind heating and scattering effects, hypothesizing that such waves might enhance efficiency compared to waves with a single wave packet, contingent upon the statistical significance of such waves.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Inertial Confinement Fusion Design Search Using Bayesian Optimization

Inertial confinement fusion (ICF) experiments rely on complex multi-physics simulation codes such as the Lawrence Livermore National Laboratory-developed HYDRA to guide design work. However, these simulations have several dozen tunable parameters and can be computationally expensive. This makes searching the parameter space challenging and time-consuming. Recently developed automated tools utilize Bayesian optimization to search these high-dimensional parameter spaces for optimal designs. The optimization tools run 2D integrated simulations in HYDRA to converge on a design that produces specified scalar or vector outputs. In this paper, we apply the Bayesian optimization tools to two common tuning scenarios. First, we tune simulation inputs to match measurements of a well-characterized experiment at the National Ignition Facility. This type of tuning is commonly performed to compensate for the use of simplified simulation settings (e.g. reduced resolution) or to account for missing physics in the simulations. Second, we search for an ICF simulation design that has a particular radiation drive profile. These optimizations replicate the kinds of tuning researchers routinely perform, but do so with significantly reduced manual effort. This approach demonstrates a powerful and efficient pathway toward autonomous, high-fidelity design optimization for future ICF experiments.

Bayesian optimization↗
DOE OSTI2025

Laboratory Measurements of n >= 3 K-shell Transition Energies of Sulfur Ions from F-like S viii to Li-like S xiv

Inner-shell transitions are ubiquitous in nonequilibrium collisionally ionized plasmas, such as supernova remnants, and in photoionized plasmas, such as outflows from active galactic nuclei and X-ray binaries. Inner-shell X-ray emission can help determine key parameters of these systems, such as ionization time, τ, and ionization parameter, ξ. Despite their importance, only theoretical inner-shell transition energies are available for many ions. To provide experimental benchmarks, we have measured the dominant n → 1 K-shell transitions of sulfur ions where n ≥ 3 from Li-like S xiv to F-like S viii using LLNL’s SuperEBIT electron beam ion trap and the NASA/GSFC EBIT Calorimeter Spectrometer (ECS). We identify over 30 spectral features and measure their energies with uncertainties in the ∼0.1–1 eV range. We compare these results to Flexible Atomic Code (FAC) and multireference Møller–Plesset (MR-MP) calculations and find differences between theory and experiment of ∼1 eV for FAC and <0.5 eV for most MR-MP calculations. We also compare these results to two widely used atomic databases, AtomDB and CHIANTI, and find discrepancies as high as 7 eV. Furthermore, many transitions are missing from these databases despite being prominent in our data.

Atomic data benchmarking↗
DOE OSTI2020

Achieving a robust grassy-ELM operation regime in CFETR

We have identified a robust grassy-edge localized mode (ELM) operation regime for future tokamak reactors. The regime exists within a pedestal top electron collisionality (ν*) window at high global poloidal beta (β p ). The existence of an upper ν* limit for grassy-ELMs is consistent with results previously reported in experiments (Oyama et al 2010 Nucl. Fusion 50 064014), while the existence of a lower ν* limit has not been reported previously. Using EPED and BOUT++, a theoretical model that quantitatively explains the physics of the grassy-ELMs within the window, which distinguishes them from the small mixed-ELMs at lower ν*, is presented for the first time. A peeling-ballooning stability boundary is obtained by scanning the operating density space. The change in density corresponds to a change in ν* that affects the pedestal bootstrap current. High β p leads to a strong Shafranov shift, which affects the flux surface averaged pressure drive. The two effects combine to create a peeling-dominated window in intermediate ν* buffered by ballooning-dominated regimes. Only the peeling-dominated regime shows a cyclic behavior in the perturbed pressure during the nonlinear simulation of an ELM crash, reminiscent of grassy-ELM dynamics. Similarly, the energy released across the separatrix is demonstrated to be significantly smaller. The quick recovery of the ELM crash is explainable by the rapid rise of a low n kink-peeling instability when the pedestal current I ped exceeds a threshold at high β p . It minimizes the excursion beyond marginal stability and is absent in the ballooning-dominated regime. Comparison with recent experiments over a range of β p and ν* strongly supports the physical picture proposed by the modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Time-domain modelling and thermometry of the CH 4 v 1 Q-branch using hybrid femtosecond/picosecond coherent anti-Stokes Raman scattering

Here, we present time-domain measurements and modelling of the CH 4 v 1 vibrational Q-branch, as well as simultaneous temperature acquisition from the CH 4 and N 2 vibrational Q-branches using broadband hybrid femtosecond/picosecond coherent anti-Stokes Raman scattering (fs/ps CARS). Accurate measurements with fs/ps CARS require knowledge of the time-domain response of the probed molecule at different temperatures and pressures. In this work, a fs/ps CARS model was developed for theCH 4 v 1 Q-branch and validated using time-resolved probe delay scans at pressures between 70 and 600 Torr and temperatures between 295 and 1000 K. A simplified modified exponential energy gap model was used to determine the linewidths of the Raman transitions within the Q-branch. Total Q-branch linewidths for the entire explored parameter space were also determined from the probe delay scans. This enabled quantitative measurements of CH 4 with fs/ps CARS as well as the temperature through the hot bands of the CH 4 v 1 Q-branch. With a broadband fs laser pulse in a two-beam phase matching scheme, the N 2 Q-branch and CH 4 v 1 Q-branch were simultaneously measured in a heated cell, and the temperatures from each molecule's spectra was evaluated. Here, the fitted temperatures corresponded well with differences ranging from 4 to 17%. It was also shown that exchanging 90% of CH 4 for N 2 or Ar at 295 K and 500 Torr did not significantly impact the normalized time-domain signal prior to a 100 picosecond probe delay. Spatially resolved one-dimensional (1-D) thermometry using the CH 4 Q-branch was also demonstrated with an imaged spatial resolution of 40 µm. With the developed model, 1-D quantitative measurements of CH 4 concentration and temperature can be measured with time resolution of tens of picoseconds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Cognitive simulation models for inertial confinement fusion: Combining simulation and experimental data

The design space for inertial confinement fusion (ICF) experiments is vast, and experiments are extremely expensive. Researchers rely heavily on computer simulations to explore the design space in search of high-performing implosions. However, ICF multiphysics codes must make simplifying assumptions, and thus deviate from experimental measurements for complex implosions. For more effective design and investigation, simulations require input from past experimental data to better predict future performance. In this work, we describe a cognitive simulation method for combining simulation and experimental data into a common, predictive model. This method leverages a machine learning technique called “transfer learning,” the process of taking a model trained to solve one task, and partially retraining it on a sparse dataset to solve a different, but related task. In the context of ICF design, neural network models are trained on large simulation databases and partially retrained on experimental data, producing models that are far more accurate than simulations alone. Here, we demonstrate improved model performance for a range of ICF experiments at the National Ignition Facility and predict the outcome of recent experiments with less than 10% error for several key observables. We discuss how the methods might be used to carry out a data-driven experimental campaign to optimize performance, illustrating the key product—models that become increasingly accurate as data are acquired.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Long range frequency chirping of Alfvén eigenmodes

In this work, a theoretical framework has been developed for an NBI scenario to model the hard non-linear evolution of global Alfvén eigenmodes (GAEs) where the adiabatic motion of phase-space structures (holes and clumps), associated with frequency chirping, occurs in generalized phase-space of slowing down energetic particles. The radial profile of the GAE is expanded using finite elements, which allows an update of the mode structure as the mode frequency chirps. Constants of motion are introduced to track the dynamics of energetic particles during frequency chirping by implementing proper action-angle variables and canonical transformations which reduce the dynamics essentially to 1D. Consequently, we specify whether the particles are drifting inward/outward as the frequency deviates from the initial MHD eigenfrequency. Using the principle of least action, we have derived the non-linear equation describing the evolution of the radial profile by varying the total Lagrangian of the system with respect to the weights of the finite elements. For the choice of parameters in this work, it is shown that the peak of the radial profile is shifted and also broadens due to frequency chirping. The time rate of frequency change is also calculated using the energy balance and we show that the adiabatic condition remains valid once it is satisfied. This model clearly illustrates the theoretical treatment to study the long-range adiabatic frequency sweeping events observed for Alfvén gap modes in real experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The Illumination of Thunderclouds by Lightning: 1. The Extent and Altitude of Optical Lightning Sources

Abstract Optical space‐based lightning sensors including NOAA's Geostationary Lightning Mapper (GLM) detect lightning though its transient illumination of the surrounding clouds. What space‐based optical lightning sensors measure is influenced by the physical attributes of the light source, the location of the source within the cloud scene, and the spatial variations in cloud composition. We focus on the lightning channels that serve as optical sources for GLM groups and flashes in this first part of our thundercloud illumination study. We match Lightning Mapping Array (LMA) sources with GLM groups and flashes during two thunderstorms to examine channel segments that are active during optical emission. We find that in each storm, the LMA sources matched with GLM groups are small (median: 3 km) compared to GLM pixels (nominal: 8 km), and preferentially come from high altitudes in the cloud (>8–10 km). The detection advantage for high‐altitude sources permits GLM to resolve faint optical pulses near the cloud top that might be missed from lower altitudes. However, the most energetic groups can be detected from all altitudes, and the largest groups largely originate at low altitudes. The relationship between group energy and illuminated area depends on flash development within the cloud medium, and flash development into different cloud regions can be identified by tracking GLM metrics of cloud illumination over time.

58 GEOSCIENCES↗
DOE OSTI2026

Suppressing proximity effects during rapid serial two-photon lithography through tuning of reaction–diffusion kinetics

The ability of two-photon lithography (TPL) to print deterministic nanoporous 3D structures is highly valuable for many applications. However, it is challenging to print such structures rapidly due to the proximity effects that cause closely spaced features to enlarge and merge. A key challenge is the limited understanding of the origins and spatiotemporal dynamics of the long-range proximity effects that extend beyond the optical focal spot. Here, we empirically investigate these proximity effects in serial TPL using custom-made acrylate-based photoresists. We demonstrate that the complex spatiotemporal dynamics of the proximity effects can be explained through the kinetics of the reaction–diffusion photopolymerization mechanisms that underlie the curing process. Specifically, we show that the proximity effects arise due to comparable timescales of reaction and diffusion of oxygen. Furthermore, we demonstrate that long-range proximity effects can be suppressed by introducing phenolic inhibitors, which reduce oxygen consumption through non-inhibiting side reactions and promote faster termination of curing reactions. These insights enabled improving the linewidths from >400 to 260 nm during printing of nanoporous 3D woodpiles at scanning speeds of 50 mm s −1 . Thus, the knowledge generated here can be applied to deterministically tune the proximity effects and enable rapid printing of high-resolution nanoporous 3D structures.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Conservative projection-based data-driven model order reduction of a fluid-kinetic spectral solver

Kinetic simulations are computationally intensive due to six-dimensional phase space discretization. Many kinetic spectral solvers use the asymmetrically weighted Hermite expansion due to its conservation and fluid-kinetic coupling properties, i.e., the lower-order Hermite moments capture and describe the macroscopic fluid dynamics, and higher-order Hermite moments describe the microscopic kinetic dynamics. We leverage this structure by developing a parametric data-driven reduced-order model based on the proper orthogonal decomposition, which projects the higher-order kinetic moments while retaining the fluid moments intact. We demonstrate analytically and numerically that the method ensures local and global mass, momentum, and energy conservation. The numerical results show that the proposed method effectively replicates the high-dimensional spectral simulations at a fraction of the computational cost and memory, as validated on the weak Landau damping and two-stream instability benchmark problems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Theoretical analysis of key factors achieving reversed magnetic shear q-profiles sustained with lower hybrid waves on EAST

The phase space analysis technique is applied to provide new insights into the fully non-inductive lower hybrid current driven (LHCD) discharges on EAST (Garofalo et al 2017 Nucl. Fusion 57 076037). The analysis shows that there are bounded and unbounded topologies of lower hybrid (LH) waves in phase space. For typical parameters on EAST, the propagation domain for 4.6 GHz LH waves is bounded, while unbounded for 2.45 GHz LH waves and one of the conditions is recognized to be good for achieving an off-axis current profile driven by 4.6 GHz LH waves on EAST. The parametric analysis on the potential power deposition (PPD) region for those experiments demonstrates that the reversed magnetic shear dominates the wave behavior, and confines the LH power absorption to the far off-axis region as long as it occurs in the 4.6 GHz dominated LHCD discharge. GENRAY/CQL3D simulations also confirm this effect of reversed shear on the power absorption of 4.6 GHz waves. A bounded propagation domain and a positive feedback loop between magnetic shear reversal and off-axis LHCD profile could explain the LH-only-sustained strong reversed magnetic shear observed in experiments on multiple machines. In contrast, there is no such effect for 2.45 GHz waves on EAST since the PPD region is less sensitive to the change of reversed shear due to the unbounded propagation domain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Synthetic modeling of soft x-ray emissivity for magnetic island analysis in LTX-β

We present a synthetic soft x-ray (SXR) forward-modeling framework that characterizes the emissivity structure of rotating magnetic islands in the Lithium Tokamak eXperiment-β across space and time. Magnetic islands associated with tearing modes produce modulations in line-integrated SXR brightness. Traditional tomographic methods struggle to resolve these structures in devices with limited sightlines, resulting in an under-determined inversion problem, or otherwise necessitate equilibrium reconstruction data, preventing their use in active control. Our approach removes this challenge by combining a fully three-dimensional ray-tracing model with a time-dependent emissivity prescription of a $m/n$ = 2/1 magnetic island. The model is validated through helical island geometry, rotation measured by magnetic diagnostics, and equilibrium constraints from PSI-Tri reconstructions. Synthetic brightness signals are generated for each photodiode sightline from a modeled emissivity profile and directly compared with experimental data from a tangential midplane SXR array. By fitting the synthetic diagnostic output to the observed brightness across time, we simultaneously infer key geometric parameters-such as island radial location, island width, and rotation frequency-without requiring full tomographic reconstruction, all with an average absolute deviation of less than 3%. This work demonstrates that forward modeling with a single tangential array can extract key island parameters in a spherical tokamak, providing a computationally efficient alternative to conventional SXR tomography and providing a pathway toward real-time magnetic-island characterization in future devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Physical mechanisms for the transition from type-III to large ELMs induced by impurity injection on EAST

Transition from type-III to large-amplitude ELMs induced by neon injection has been observed in the EAST tokamak at overlapping q 95 space between large and small ELMs. With neon injection, pedestal density gradient shows a remarkable increase accompanied by some decrease of pedestal electron temperature, and consequently the pressure gradient increases moderately and edge bootstrap current has minimal change. Further experiment demonstrates that the occurrence of large ELMs after neon injection is highly correlated with the change in edge density. Linear peeling-ballooning stability analysis indicates that the large ELM case is more unstable than the type-III ELM case during the ELM transition. A scan of pedestal density gradient in linear stability analysis shows that the direct destabilizing effect of steep pedestal density gradient on peeling-ballooning instabilities via two-fluid effects could also facilitate the transition to large ELMs. In conclusion, these results could provide more insight into the role of pedestal density gradient on pedestal stability and ELM behavior.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗