Search NASA⌕ Search

SEARCH · Search NASA

Results for “Gas dynamics”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Numerical simulation of compressible fluid-dynamics in the chamber of inertial fusion energy systems

Here, this paper aims to establish new and innovative modeling capabilities for analyzing chambers in Inertial Fusion Energy (IFE) systems. IFE is emerging as a promising method to achieve fusion power production, but several challenges must be overcome to develop an IFE pilot plant or deploy commercial IFE systems. These challenges are both theoretical and technical, encompassing a deeper understanding of the underlying physical phenomena and the development of new technologies and materials. One of the needs is to develop mathematical models to describe IFE systems and numerical tools to simulate them. This paper contributes to this endeavor by presenting a new OpenFOAM solver for IFE systems, focusing on gas dynamics in their chambers. The analysis and development of chamber designs will play a significant role in the transition from single-shot experiments to high-repetition rates, as there is a need to protect the chamber walls from the intense radiation fields produced by fusion reactions. A promising design option, normally referred to as thick wall chamber design, consists in using lithium or molten salt jet arrays within the chamber. A critical phenomenon is the venting of high-pressure gases from the center to the external part of the chamber, passing through the blanket array. This process involves the propagation and attenuation of strong pressure waves, requiring suitable modeling approaches for compressible fluid-dynamics. The solver proposed in this work implements a multi-material hydrodynamics model tailored to accurately describe the non-linear propagation of pressure waves while avoiding numerical oscillation issues typical of high-velocity compressible simulation. This solver is verified against numerical test cases, validated against experimental data, and applied to the analysis of the High-Yield Lithium-Injection Fusion-Energy (HYLIFE-I) concept. The relevance of this paper is threefold. Firstly, it contributes to developing and testing modeling approaches for compressible fluid-dynamics phenomena, with specific focus on the new and unexplored topic of IFE thick-liquid-wall blanket modeling. Secondly, it marks one of the first applications of the OpenFOAM library in the research field of IFE systems. Finally, the investigated problem is of practical interest for IFE developers, as it provides useful indications about relevant phenomena in pressure wave propagation in the chamber of these systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Shock physics and shadowgraphic measurements of laser-produced cerium plasmas

Shadowgraphic measurements are combined with theory on gas-dynamics to investigate the shock physics associated with nanosecond laser ablation of cerium metal targets. Time-resolved shadowgraphic imaging is performed to measure the propagation and attenuation of the laser-induced shockwave through air and argon atmospheres at various background pressures, where stronger shockwaves characterized by higher propagation velocities are observed for higher ablation laser irradiances and lower pressures. The Rankine-Hugoniot relations are also employed to estimate the pressure, temperature, density, and flow velocity of the shock-heated gas located immediately behind the shock front, predicting larger pressure ratios and higher temperatures for stronger laser-induced shockwaves.

42 ENGINEERING↗
DOE OSTI2025

Design and Engineering of LUPIN: A Test-Bed Radio-Frequency Ion Source for Enhanced Neutral Beam Injection on DIII-D

The Large, Uniform Plasma for Ionizing Neutrals (LUPIN) is a radio-frequency (RF) inductively coupled plasma (ICP) chamber for demonstrating plasma performance of an RF ICP positive ion source upgrade for the DIII-D neutral beam injection (NBI) system. LUPIN will be used to investigate ion source physics, including neutral gas dynamics, plasma density uniformity, interactions with Faraday shields, and power coupling to novel RF antenna designs. LUPIN has an RF generator capable of delivering 20 kW of power at 2 MHz, which is coupled into a cylindrical quartz vessel measuring 20 cm in length and 10 cm in radius. This configuration matches the power density requirements for a full-scale ion source. Target hydrogen and deuterium plasma densities exceeding 10 18 m -3 would relate to extracted ion current densities of 2100 A/m 2 for 10s. Vacuum conductance and gas flow calculations predict a maximum achievable neutral gas flow rate of 1675 Pa ⋅ L/s at 5 Pa of He, which mimics the gas flow of the DIII-D NBI system. Designs have been developed for an internal Faraday shield to mitigate heat flux and ion sputtering on the dielectric vessel. Thermomechanical finite element simulations demonstrated the Faraday shield design to be capable of withstanding anticipated heat loads from worst-case operation scenarios. Finally, results of upcoming experimental investigations on LUPIN will guide the design of a full-scale prototype for DIII-D integration.

Faraday shield↗
DOE OSTI2024

Design and Engineering of LUPIN: A Test-Bed Radiofrequency Ion Source for Enhanced Neutral Beam Injection on DIII-D

The Large, Uniform Plasma for Ionizing Neutrals (LUPIN) is an RF inductively coupled plasma (ICP) chamber for demonstrating plasma performance of an RF ICP positive ion source upgrade for the DIII-D Neutral Beam Injection (NBI) system. LUPIN will be used to investigate ion source physics, including neutral gas dynamics, plasma density uniformity, interactions with Faraday shields, and power coupling to novel RF antenna designs. LUPIN has an RF generator capable of delivering 20 kW of power at 2 MHz, which is coupled into a cylindrical quartz vessel measuring 20 cm in length and 10 cm in radius. This configuration matches the power density requirements for a full-scale ion source. Target hydrogen and deuterium plasma densities exceeding 10 18 m -3 would relate to extracted ion current densities of 2100 A/m 2 for 10 s. Vacuum conductance and gas flow calculations predict a maximum achievable neutral gas flow rate of 1675 Pa·L/s at 5 Pa of He, which mimics the gas flow of the DIII-D NBI system. Designs have been developed for an internal Faraday shield to mitigate heat flux and ion sputtering on the dielectric vessel. Thermomechanical finite element simulations demonstrated the Faraday shield design to be capable of withstanding anticipated heat loads from worst-case operation scenarios. Results of upcoming experimental investigations on LUPIN will guide the design of a full-scale prototype for DIII-D integration.

Faraday shield↗
DOE OSTI2020

A Novel use of Direct Simulation Monte-Carlo to Model Dynamics of COVID-19 Pandemic Spread

In this report, we evaluate a novel method for modeling the spread of COVID-19 pandemic. In this new approach we leverage methods and algorithms developed for fully-kinetic plasma physics simulations using Particle-In-Cell (PIC) Direct Simulation Monte-Carlo (DSMC) models. This approach then leverages Sandia-unique simulation capabilities, and High-Performance Computer (HPC) resources and expertise in particle-particle interactions using stochastic processes. Our hypothesis is that this approach would provide a more efficient platform with assumptions based on physical data that would then enable the user to assess the impact of mitigation strategies and forecast different phases of infection. This work addresses key scientific questions related to the assumptions this new approach must make to model the interactions of people using algorithms typically used for modeling particle interactions in physics codes (kinetic plasma, gas dynamics). The model developed uses rational/physical inputs while also providing critical insight; the results could serve as inputs to, or alternatives for, existing models. The model work presented was developed over a four-week time frame, thus far showing promising results and many ways in which this model/approach could be improved. This work is aimed at providing a proof-of-concept for this new pandemic modeling approach, which could have an immediate impact on the COVID-19 pandemic modeling, while laying a basis to model future pandemic scenarios in a manner that is timely and efficient. Additionally, this new approach provides new visualization tools to help epidemiologists comprehend and articulate the spread of this and other pandemics as well as a more general tool to determine key parameters needed in order to better predict pandemic modeling in the future. In the report we describe our model for pandemic modeling, apply this model to COVID-19 data for New York City (NYC), assess model sensitivities to different inputs and parameters and , finally, propagate the model forward under different conditions to assess the effects of mitigation and associated timing. Finally, our approach will help understand the role of asymptomatic cases, and could be extended to elucidate the role of recovered individuals in the second round of the infection, which is currently being ignored.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Gauging the Impact of Cosmic-ray Feedback on the Stellar Initial Mass Function

Cosmic rays (CRs) drive ionization and influence gas dynamics in molecular clouds (MCs), potentially impacting the resulting star formation outcomes. Although previous simulations of individual star formation have included methods for CR transport (CRT), none have been large enough to resolve the stellar initial mass function (IMF). We conduct numerical simulations following the collapse of a 20,000 M ⊙ MC and the subsequent star formation including CRT, both with and without CRs accelerated by winds from the young massive stars, and compare against a non-CRT simulation. We show that after the first massive stars form, the cavity produced by feedback is more pronounced in the CRT simulations because the external CRs are able to propagate inward and compress the gas into higher-density structures. This increases the subsequent star formation in the cloud; by the end of the simulation, the star formation efficiency (SFE) in the CRT simulation including stellar wind CRs is 43% higher than the non-CRT simulation. The IMF is also top-heavy in comparison, with a slope above 1 M ⊙ that is shallower by ∼20%. These effects are also present in the simulation without wind-accelerated CRs, but they are not as pronounced; the SFE is only 16% higher than the non-CRT simulation, and the IMF high-mass slope is shallower by ∼10%. These results may explain some of the observed top-heavy IMFs, which typically occur in high-CR environments such as the Galactic center.

Fitz Axen, Margot [Oak Ridge National Laboratory (↗
DOE OSTI2019

Detonation propagation for shock-driven, subsonic and supersonic confiner flow

In this work, we study the compressible flow dynamics of two-dimensional, steady detonation wave propagation in a high explosive (HE) confined by aluminium (Al) or stainless steel (SS), outside of which is an air layer. We examine how the thickness of the confinement affects the subsonic detonation driving zone structure (DDZ) and the detonation speed (D 0 ), demonstrating a strong dependence on whether the oblique shock-driven flow in the confiner is supersonic, as for SS, or subsonic, as for Al. A characteristic path analysis is used to examine the information flow from the material boundaries through the supersonic flow regions in both the HE and confiner that can impact the sonic surfaces bounding the subsonic flow regions. It is shown that the nature of gas-dynamic wave reflection off the SS–air or Al–air boundary can significantly influence the DDZ and D 0 .

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Fractal Scaling of Explosively Driven Product Gases

ABSTRACT Characterization of the interface between explosive product gases and ambient air in an explosion is a complicated task due to the turbulent mixing and inherently three‐dimensional expansion of the interface. This study aims to quantify the evolution of the interface as a temporally varying Hausdorff dimension. Two test series were conducted with Composition C‐4 charges with masses of 105 and 880 g. Imaging data were collected from the time of detonation until shock wave detachment using ultra‐high‐speed cameras. Gas cloud profiles were extracted using automated image processing algorithms, and the Hausdorff dimension of these two‐dimensional slices of the gas cloud was then estimated using boxcounting algorithms. When scaled with standard gas dynamic nondimensional scalings, the Hausdorff dimension of all explosive events appears to collapse towards a single curve. The fireball was initially nonfractal and began to develop fractal properties as the shock wave separated from the detonation products. Artificial perturbation of the charged surface had no detectable impact on the evolution of the Hausdorff dimension in the early development of the fireball outside of error, despite visible phenomenological differences in the early development of mixing on the fireball surface.

42 ENGINEERING↗
DOE OSTI2023

Combustion of 10–100 μm aluminum droplets in detonation products gases

Here we describe a two-phase model of combustion effects in aluminized high explosive (HE) charges. It is based on: (i) a Gas Dynamic Model of the expansion of the detonation product gases and their turbulent combustion with air; and (ii) a Heterogeneous Continuum Model of aluminum (Al) droplets and their combustion with the detonation product gases. Initial conditions are based on an analytical similarity solution for a cylindrical Chapman-Jouguet (CJ) detonation propagating at the CJ detonation velocity. The CJ jump conditions are computed at the thermodynamic equilibrium state by the Cheetah code, assuming the Al droplets are inert. We assume that the Al is 10 % of the charge mass and occurs as droplets at the CJ state. Different initial droplet diameters, ranging from 10 to 100 microns, are studied. A hydrodynamic combustion model based on large Damköhler numbers is employed in this study, -3σ(1 + 0.276 $\sqrt{Re}$)/($Kd^2_w$). It has a square-root dependence on the Reynolds number (Re) and inversed-squared dependence on the droplet diameter (d w ). The burnout time (t B ) of the Al droplets has a three-halves dependence on the droplet diameter, t B ~ ($d^0_w$) 3/2 . After burnout, the detonation products act as detonation products of the HE charge with active Al. They turbulently mix with air and form a combustion layer on the outer edge of the fireball. Details of the two-phase model, initial conditions and evolution of the flow field will be described.

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

Neural-network based collision operators for the Boltzmann equation

Kinetic gas dynamics in rarefied and moderate-density regimes have complex behavior associated with collisional processes. These processes are generally defined by convolution integrals over a high-dimensional space (as in the Boltzmann operator), or require evaluating complex auxiliary variables (as in Rosenbluth potentials in Fokker-Planck operators) that are challenging to implement and computationally expensive to evaluate. In this work, we develop a data-driven neural network model that augments a simple and inexpensive BGK collision operator with a machine-learned correction term, which improves the fidelity of the simple operator with a small overhead to overall runtime. The composite collision operator has a tunable fidelity and, in this work, is trained using and tested against a direct-simulation Monte-Carlo (DSMC) collision operator.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Continuously bounds-preserving discontinuous Galerkin methods for hyperbolic conservation laws

For finite element approximations of transport phenomena, it is often necessary to apply a form of limiting to ensure that the discrete solution remains well-behaved and satisfies physical constraints. However, these limiting procedures are typically performed at discrete nodal locations, which is not sufficient to ensure the robustness of the scheme when the solution must be evaluated at arbitrary locations (e.g., for adaptive mesh refinement, remapping in arbitrary Lagrangian–Eulerian solvers, overset meshes, etc.). In this work, a novel limiting approach for discontinuous Galerkin methods is presented which ensures that the solution is continuously bounds-preserving (i.e., across the entire solution polynomial) for any arbitrary choice of basis, approximation order, and mesh element type. Through a modified formulation for the constraint functionals, the proposed approach requires only the solution of a single spatial scalar minimization problem per element for which a highly efficient numerical optimization procedure is presented. Here, the efficacy of this approach is shown in numerical experiments by enforcing continuous constraints in high-order unstructured discontinuous Galerkin discretizations of hyperbolic conservation laws, ranging from scalar transport with maximum principle preserving constraints to compressible gas dynamics with positivity-preserving constraints.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A note on higher-order and nonlinear limiting approaches for continuously bounds-preserving discontinuous Galerkin methods

In Dzanic (2024), a limiting approach for high-order discontinuous Galerkin schemes was introduced which allowed for imposing constraints on the solution continuously (i.e., everywhere within the element). While exact for linear constraint functionals, this approach only imposed a sufficient (but not the minimum necessary) amount of limiting for nonlinear constraint functionals. This short note shows how this limiting approach can be extended to allow exactness for general nonlinear quasiconcave constraint functionals through a nonlinear limiting procedure, reducing unnecessary numerical dissipation. Finally, some examples are shown for nonlinear pressure and entropy constraints in the compressible gas dynamics equations, where both analytic and iterative approaches are used.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

High-order limiting methods using maximum principle bounds derived from the Boltzmann equation I: Euler equations

The use of limiting methods for high-order numerical approximations of hyperbolic conservation laws generally requires defining an admissible region/bounds for the solution. In this work, we present a novel approach for computing solution bounds and limiting for the Euler equations through the kinetic representation provided by the Boltzmann equation, which allows for extending limiters designed for linear advection directly to the Euler equations. Given an arbitrary set of solution values to compute bounds over (e.g., numerical stencil) and a desired linear advection limiter, the proposed approach yields an analytic expression for the admissible region of particle distribution function values, which may be numerically integrated to yield a set of bounds for the density, momentum, and total energy. Further, these solution bounds are shown to preserve positivity of density/pressure/internal energy and, when paired with a limiting technique, can robustly resolve strong discontinuities while recovering high-order accuracy in smooth regions without any ad hoc corrections (e.g., relaxing the bounds). This approach is demonstrated in the context of an explicit unstructured high-order discontinuous Galerkin/flux reconstruction scheme for a variety of difficult problems in gas dynamics, including cases with extreme shocks and shock-vortex interactions. Furthermore, this work presents a foundation for limiting techniques for more complex macroscopic governing equations that can be derived from an underlying kinetic representation for which admissible solution bounds are not well-understood.

42 ENGINEERING↗
DOE OSTI2025

Drift-cyclotron loss-cone instability in 3-D simulations of a sloshing-ion simple mirror

The kinetic stability of collisionless, sloshing beam-ion (45° pitch angle) plasma is studied in a three-dimensional (3-D) simple magnetic mirror, mimicking the Wisconsin high-temperature superconductor axisymmetric mirror experiment. The collisional Fokker–Planck code CQL3D-m provides a slowing-down beam-ion distribution to initialize the kinetic-ion/fluid-electron code Hybrid-VPIC, which then simulates free plasma decay without external heating or fuelling. Over 1 – 10 μs, drift-cyclotron loss-cone (DCLC) modes grow and saturate in amplitude. The DCLC scatters ions to a marginally stable distribution with gas-dynamic rather than classical-mirror confinement. Sloshing ions can trap cool (low-energy) ions in an electrostatic potential well to stabilize DCLC, but DCLC itself does not scatter sloshing beam-ions into the said well. Instead, cool ions must come from external sources such as charge-exchange collisions with a low-density neutral population. Manually adding cool ∼1keV ions improves beam-ion confinement several-fold in Hybrid-VPIC simulations, which qualitatively corroborates prior measurements from real mirror devices with sloshing ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Quantification of Hydrogen Isotopes Utilizing Raman Spectroscopy Paired with Chemometric Analysis for Application across Multiple Systems

On-line and real-time analysis of a chemical process is a major analytical challenge that can drastically change the way the chemical industry or chemical research operates. With in situ analyses, new and powerful understanding of chemistry can be gained; however, building robust tools for long-term monitoring faces many challenges that include compensating for instrument drift, instrument replacement, and sensor or probe replacement. Accounting for these changes by recollecting calibration data and rebuilding quantification models can be costly and time consuming. Here, in this study, methods to overcome these challenges are demonstrated with an application of Raman spectroscopy to monitoring hydrogen isotopes with varied speciation within dynamic gas streams. Specifically, chemical data science tools such as chemometric modeling are leveraged along with several examples of calibration transfer approaches. Furthermore, the optimization of instrument and sensors cell parameters for targeted gas phase analyses is discussed. While the particular focus on hydrogen is highly beneficial within the nuclear energy sector, mechanisms built and demonstrated here are widely applicable to optical spectroscopy monitoring in numerous other chemical systems that can be leveraged in other hazardous processes.

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

Plasma Properties in the Earth's Magnetosheath Near the Subsolar Magnetopause: Implications for Geocoronal Density Estimates

Combined in situ ion measurements and remote sensing of energetic neutral atoms are used to determine the geocoronal Hydrogen density at large (∼10 R E ) distances from the Earth. This method for determining the geocoronal density requires global magnetospheric modeling. Observations in the Earth's subsolar magnetosheath from the Magnetospheric Multiscale mission are used to determine the accuracy of using global models to predict the geocoronal density. On average, gas dynamic and magnetohydrodynamic (MHD) models and observations are in reasonable agreement, with differences <25%. In addition, the MHD model subsolar magnetopause is about 0.5 R E sunward of the observed location. However, variations around averages are large (up to a factor of 2), indicating that global models introduce relatively large uncertainties in geocoronal density estimates. Finally, the critical ion flux in the Interstellar Boundary Explorer IBEX‐Hi energy range is often minimally affected by fluctuations of a factor of 2 in the density.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Feld-induced modulation of two-dimensional electron gas at LaAlO 3 /SrTiO 3 interface by polar distortion of LaAlO 3

Since the discovery of two-dimensional electron gas at the LaAlO 3 /SrTiO 3 interface, its intriguing physical properties have garnered significant interests for device applications. Yet, understanding its response to electrical stimuli remains incomplete. Our in-situ transmission electron microscopy analysis of a LaAlO 3 /SrTiO 3 two-dimensional electron gas device under electrical bias reveals key insights. Inline electron holography visualized the field-induced modulation of two-dimensional electron gas at the interface, while electron energy loss spectroscopy showed negligible electromigration of oxygen vacancies. Instead, atom-resolved imaging indicated that electric fields trigger polar distortion in the LaAlO 3 layer, affecting two-dimensional electron gas modulation. This study refutes the previously hypothesized role of oxygen vacancies, underscoring the lattice flexibility of LaAlO 3 and its varied polar distortions under electric fields as central to two-dimensional electron gas dynamics. These findings open pathways for advanced oxide nanoelectronics, exploiting the interplay of polar and nonpolar distortions in LaAlO 3 .

36 MATERIALS SCIENCE↗