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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows

High-order Godunov methods for gas dynamics have become a standard tool for simulating different classes of astrophysical flows. Their accuracy is mostly determined by the spatial interpolant used to reconstruct the pair of Riemann states at cell interfaces and by the Riemann solver that computes the interface fluxes. In most Godunov-type methods, these two steps can be treated independently, so that many different schemes can in principle be built from the same numerical framework. Because astrophysical simulations often test out the limits of what is feasible with the computational resources available, it is essential to find the scheme that produces the numerical solution with the desired accuracy at the lowest computational cost. However, establishing the best combination of numerical options in a Godunov-type method to be used for simulating a complex hydrodynamic problem is a nontrivial task. In fact, formally more accurate schemes do not always outperform simpler and more diffusive methods, especially if sharp gradients are present in the flow. For this work, we used our fully compressible Seven-League Hydro (SLH) code to test the accuracy of six reconstruction methods and three approximate Riemann solvers on two- and three-dimensional (2D and 3D) problems involving subsonic flows only. We considered Mach numbers in the range from 10 −3 to 10 −1 , which are characteristic of many stellar and geophysical flows. In particular, we considered a well-posed, 2D, Kelvin–Helmholtz instability problem and a 3D turbulent convection zone that excites internal gravity waves in an overlying stable layer. Although the different combinations of numerical methods converge to the same solution with increasing grid resolution for most of the quantities analyzed here, we find that (i) there is a spread of almost four orders of magnitude in computational cost per fixed accuracy between the methods tested in this study, with the most performant method being a combination of a low-dissipation Riemann solver and a sextic reconstruction scheme; (ii) the low-dissipation solver always outperforms conventional Riemann solvers on a fixed grid when the reconstruction scheme is kept the same; (iii) in simulations of turbulent flows, increasing the order of spatial reconstruction reduces the characteristic dissipation length scale achieved on a given grid even if the overall scheme is only second order accurate; (iv) reconstruction methods based on slope-limiting techniques tend to generate artificial, high-frequency acoustic waves during the evolution of the flow; and (v) unlimited reconstruction methods introduce oscillations in the thermal stratification near the convective boundary, where the entropy gradient is steep.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Advanced manufacturing provides tailor-made solutions for crystallography with x-ray free-electron lasers

Serial crystallography at large facilities, such as x-ray free-electron lasers and synchrotrons, evolved as a powerful method for the high-resolution structural investigation of proteins that are critical for human health, thus advancing drug discovery and novel therapies. However, a critical barrier to successful serial crystallography experiments lies in the efficient handling of the protein microcrystals and solutions at microscales. Microfluidics are the obvious approach for any high-throughput, nano-to-microliter sample handling, that also requires design flexibility and rapid prototyping to deal with the variable shapes, sizes, and density of crystals. Here, we discuss recent advances in polymer 3D printing for microfluidics-based serial crystallography research and present a demonstration of emerging, large-scale, nano-3D printing approaches leading into the future of 3D sample environment and delivery device fabrication from liquid jet gas-dynamic virtual nozzles devices to fixed-target sample environment technology.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Collisionless adiabatic afterglow

In this work, we study, by numerical and analytical means, the evolution of a collisionless plasma initiated between absorbing walls. The ensuing flow is described by rarefaction waves that travel inward from the boundaries, interact, and eventually vanish after crossing through, leading up to the asymptotic stage of the decay. Particle simulations indicate that the kinetic evolution strongly resembles one found in isentropic gas dynamics. Namely, a very gradual density profile forms in the expanding central region where the rarefaction waves interact, with an accompanying linear velocity profile. Asymptotically, the density falls off as 1/t. The density and the flux at the boundary show little variation over the period when rarefaction waves still exist. Plasma potential, on the other hand, drops quite rapidly (on the underlying ion-acoustic timescale) to less than initial electron temperature Te when over 70% of the particles still remain in the system. This is due to electron kinetics being governed by conservation of adiabatic invariant in a slowly varying potential well. Analytical model of the velocity distribution is presented to explain the simulations. The results have implications for afterglow plasmas used in material processing and also for ion-extraction devices. One property of potential interest is good uniformity of the decaying plasma that occurs after approximately one ion-acoustic time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A drift kinetic model for the expander region of a magnetic mirror

Here, we present a drift kinetic model for the free expansion of a thermal plasma out of a magnetic nozzle. This problem relates to plasma space propulsion systems, natural environments such as the solar wind, and end losses from the expander region of mirror magnetically confined fusion concepts such as the gas dynamic trap. The model incorporates trapped and passing orbit types encountered in the mirror expander geometry and maps to an upstream thermal distribution. This boundary condition and quasineutrality require the generation of an ambipolar potential drop of ~5T e /e, forming a thermal barrier for the electrons. The model for the electron and ion velocity distributions and fluid moments is confirmed with data from a fully kinetic simulation. Finally, the model is extended to account for a population of fast sloshing ions arising from neutral beam heating within a magnetic mirror, again resulting in good agreement with a corresponding kinetic simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Characterization of a new ultra-high pressure shock tube facility for combustion and propulsion studies

A new shock tube facility has been designed, constructed, and characterized at the University of Central Florida. This facility is capable of withstanding pressures of up to 1000 atm, allowing for combustion diagnostics of extreme conditions, such as in rocket combustion chambers or in novel power conversion cycles. For studies with toxic gas impurities, the high initial pressures required the development of a gas delivery system to ensure the longevity of the facility and the safety of the personnel. Data acquisition and experimental propagation were implemented with remote access to ensure safety, paired with a LabVIEW- and Python-based user interface. Thus far, test pressures of 270 atm, blast pressures of 730 atm, and temperatures approaching 10 000 K have been achieved. The extreme limitations of this facility allow for emission spectroscopy to be performed during the oxidation of fuel mixtures, e.g., alkanes diluted in argon and carbon dioxide. Ignition delay times were determined and compared to simulations using chemical kinetic mechanisms. The design, experimental procedures, processes of analysis, and uncertainty determination are outlined, and typical pressure profiles are compared with a new gas dynamics solver and empirical correlations developed across multiple shock tube facilities. Preliminary reactive mixture analyses are included with further investigation of the mixtures outlined.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Thermal-fluctuation effects on small-scale statistics in turbulent gas flow

Kolmogorov's theory of turbulence assumes that the small-scale turbulent structures in the energy cascade are universal and are determined by the energy dissipation rate and the kinematic viscosity alone. However, thermal fluctuations, absent from the continuum description, terminate the energy cascade near the Kolmogorov length scale. Here, we propose a simple superposition model to account for the effects of thermal fluctuations on small-scale turbulence statistics. Furthermore, for compressible Taylor–Green vortex flow, we demonstrate that the superposition model in conjunction with data from direct numerical simulation of the Navier–Stokes equations yields spectra and structure functions that agree with the corresponding quantities computed from the direct simulation Monte Carlo method of molecular gas dynamics, verifying the importance of thermal fluctuations in the dissipation range.

42 ENGINEERING↗
DOE OSTI2021

First evidence of dominant influence of E × B drifts on plasma cooling in an advanced slot divertor for tokamak power exhaust

Addressing power exhaust in tokamaks is presently recognized as one of the major remaining open issues for the development of fusion reactors. At the forefront of this endeavor is the effort to develop an advanced divertor by maximizing dissipation of plasma power and momentum inside the divertor. Here, we demonstrate, for the first time, that the electromagnetic (E × B) drifts exert a key influence on plasma and gas dynamics in a new advanced slot divertor in the DIII-D tokamak, named the small angle slot (SAS). SAS leverages the effect of drifts to achieve a highly dissipative divertor with electron temperature T e ≲ 10 eV over a wide range of plasma densities, for ion B × ∇B away from the divertor, as used for advanced tokamak operation on DIII-D. Modeling with the SOLPS-ITER code shows that for this drift direction, the E × B flow carries particles toward the outer common flux region (CFR) via the private flux region (PFR), reinforcing neutral recycling and enhancing divertor dissipation. In contrast, for the opposite field direction, the E × B flow carries particles away from the outer CFR into the PFR, offsetting the anticipated SAS geometric effects. Furthermore, this finding is an important step in the understanding of the behavior of advanced divertors for achieving a power exhaust solution for fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Numerical investigation of liquid wall ablation in inertial fusion energy chambers

This paper presents a novel approach for modeling liquid wall ablation in liquid wall-protected inertial fusion energy (IFE) chambers. These systems are promising candidates for the implementation of fusion technology, yet significant gaps remain in understanding the underlying physical processes and their implications for design. Following target ignition, a portion of the fusion energy is released as x-rays, which deposit their energy into an array of liquid jets, leading to partial vaporization. Accurately modeling this heat deposition and vaporization process remains challenging due to the complex geometries typical of (pre-conceptual) IFE chamber designs. Furthermore, the subsequent expansion of vaporized material into the chamber’s vacuum environment poses difficulties for conventional CFD methods based on continuum assumptions, which can lead to significant inaccuracies. To address some aspects of these challenges, this work introduces a ray-tracing-based methodology to map the spatial distribution of ablated material in liquid wall-protected systems. In addition, a vacuum-tracking scheme is developed to extend the applicability of an OpenFOAM-based solver to gas dynamics in rarefied environments. The proposed approach has been verified through numerical benchmarks and applied to a practical case involving the HYLIFE-II (High Yield Lithium Injection Fusion Energy) chamber. The methodology advances the modeling capabilities for liquid wall-protected IFE systems and provides valuable tools to support their design and optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

flash -light on the ring : hydrodynamic simulations of expanding supernova shells near supermassive black holes

ABSTRACT The way supermassive black holes (SMBHs) in Galactic Centres (GCs) accumulate their mass is not completely determined. At large scales, it is governed by galactic encounters, mass inflows connected to spirals arms and bars, or due to expanding shells from supernova (SN) explosions in the central parts of galaxies. The investigation of the latter process requires an extensive set of gas dynamical simulations to explore the multidimensional parameter space needed to frame the phenomenon. The aims of this paper are to extend our investigation of the importance of SNe for inducing accretion on to an SMBH and carry out a comparison between the fully hydrodynamic code flash and the much less computationally intensive code ring, which uses the thin shell approximation. We simulate 3D expanding shells in a gravitational potential similar to that of the GC with a variety of homogeneous and turbulent environments. In homogeneous media, we find convincing agreement between flash and ring in the shapes of shells and their equivalent radii throughout their whole evolution until they become subsonic. In highly inhomogeneous, turbulent media, there is also a good agreement of shapes and sizes of shells, and of the times of their first contact with the central 1-pc sphere, where we assume that they join the accretion flow. The comparison supports the proposition that an SN occurring at a galactocentric distance of 5 pc typically drives 1–3 M⊙ into the central 1 pc around the GC.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Effect of slip on vortex shedding from a circular cylinder in a gas flow

Most studies of vortex shedding from a circular cylinder in a gas flow have explicitly or implicitly assumed that the no-slip condition applies on the cylinder surface. To investigate the effect of slip, vortex shedding is simulated using molecular gas dynamics (the direct simulation Monte Carlo method) and computational fluid dynamics (the incompressible Navier-Stokes equations with a slip boundary condition). A Reynolds number of 100, a Mach number of 0.3, and a corresponding Knudsen number of 0.0048 are examined. For these conditions, compressibility effects are small, and periodic laminar vortex shedding is obtained. Slip on the cylinder is varied using combinations of diffuse and specular molecular reflections with accommodation coefficients from zero (maximum slip) to unity (minimum slip). Although unrealistic, bounce-back molecular reflections are also examined because they approximate the no-slip boundary condition (zero slip). The results from both methods are in reasonable agreement. The shedding frequency increases slightly as the accommodation coefficient is decreased, and shedding ceases at low accommodation coefficients (large slip). The streamwise and transverse forces decrease as the accommodation coefficient is decreased. Based on the good agreement between the two methods, computational fluid dynamics is used to determine the critical accommodation coefficient below which vortex shedding ceases for Reynolds numbers of 60–100 at a Mach number of 0.3. Finally, conditions to observe the effect of slip on vortex shedding appear to be experimentally realizable, although challenging.

42 ENGINEERING↗
DOE OSTI2023

CO Laser Absorption Measurements During Syngas Combustion at High Pressure

Abstract Syngas is a desirable fuel for combustion in the Allam-Fetvedt cycle, which involves combustion under supercritical-CO2 conditions. While some work has been conducted in collecting ignition delay times (IDT) at the extreme pressures required by these systems, significant model deficiencies remain. Additionally, considerable barriers in terms of nonideal gas dynamic effects have been shown for these experiments in shock tubes. Further investigation into the fundamental combustion kinetics of H2/CO/CO2 mixtures is required. Time-resolved speciation measurements for target species have been shown to better aid in improving the understanding of underlying chemical kinetics than global ignition delay time measurements. Therefore, laser absorption measurements of CO were measured behind reflected shock waves during combustion of syngas at 5 and 10 bar and temperatures between 1080 and 2100 K. The mixtures investigated utilized H2-to-CO ratios of 1:1 and 1:4, respectively, each at stoichiometric conditions, allowing for discussions of the effect of initial fuel composition. A ratio of fuel to CO2 of 1:2 was also utilized to represent commercially available syngas. The mixtures were diluted in helium and argon (20% He, 76.5% Ar) to minimize thermal effects and to expedite CO thermal relaxation during the experiment. The resulting CO time histories were then compared to modern chemical kinetics mechanisms, and disagreement is seen for this system, which is assumed to be fairly well known. This study elucidates particular chemistry that needs improvement in moving toward a better understanding of syngas combustion at elevated pressures.

Engineering↗
DOE OSTI2025

Parallel-in-Time Solution of Hyperbolic PDE Systems via Characteristic-Variable Block Preconditioning

We consider the parallel-in-time solution of both linear and nonlinear hyperbolic partial differential equation (PDE) systems in one spatial dimension. In the nonlinear setting, the discretized equations are solved with a preconditioned residual iteration based on a global linearization. The linear(ized) equation systems are approximately solved parallel-in-time using a block preconditioner applied in the characteristic variables of the underlying linear(ized) hyperbolic PDE. This change of variables is motivated by the observation that intervariable coupling between characteristic variables is weak, at least locally where spatio-temporal variations in the eigenvectors of the associated flux Jacobian are sufficiently small, while that between the original variables is not. For an ℓ-dimensional system of PDEs, applying the preconditioner consists of solving a sequence of ℓ scalar linear(ized)-advection-like problems, each associated with a different characteristic wave-speed in the underlying linear(ized) PDE. Furthermore, we approximately solve these linear advection problems using multigrid reduction-in-time (MGRIT); however, any other suitable parallel-in-time method could be used. Numerical examples are shown for the (linear) acoustics equations in heterogeneous media and for the (nonlinear) shallow water equations and Euler equations of gas dynamics with shocks and rarefactions. For many test problems, the solver converges in just a handful of iterations and with mesh-independent convergence rates.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

On High-Order/Low-Order and Micro-Macro Methods for Implicit Time-Stepping of the BGK Model

In this paper, a high-order/low-order (HOLO) method is combined with a micro-macro (MM) decomposition to accelerate iterative solvers in fully implicit time-stepping of the Bhatnagar–Gross–Krook (BGK) equation for gas dynamics. The MM formulation represents a kinetic distribution as the sum of a local Maxwellian and a perturbation. In highly collisional regimes, the perturbation away from initial and boundary layers is small and can be compressed to reduce the overall storage cost of the distribution. The convergence behavior of the MM methods, the usual HOLO method, and the standard source iteration method is analyzed on a linear BGK model. Both the HOLO and MM methods are implemented using a discontinuous Galerkin (DG) discretization in phase space, which naturally preserves the consistency between high- and low-order models required by the HOLO approach. Furthermore, the accuracy and performance of these methods are compared on the Sod shock tube problem and a sudden wall heating boundary layer problem. Overall, the results demonstrate the robustness of the MM and HOLO approaches and illustrate the compression benefits enabled by the MM formulation when the kinetic distribution is near equilibrium.

BGK model↗
DOE OSTI2025

Cosmological Hydrodynamics at Exascale: A Trillion-Particle Leap in Capability

Resolving the most fundamental questions in cosmology requires simulations that match the scale, fidelity, and physical complexity demanded by next-generation sky surveys. To achieve the realism needed for this critical scientific partnership, detailed gas dynamics must be treated self-consistently with gravity for end-to-end modeling of structure formation. Exascale computing enables simulations that span survey-scale volumes while incorporating key astrophysical processes that shape complex cosmic structures. We present results from CRK-HACC, a cosmological hydrodynamics code built for extreme scalability. Using separation-of-scale techniques, GPU-resident tree solvers, in situ analysis pipelines, and multi-tiered I/O, CRK-HACCexecuted Frontier-E: a four trillion particle full-sky simulation, over an order of magnitude larger than previous efforts. The run achieved 513.1 PFLOPs peak performance, processing 46.6 billion particles per second and writing more than 100 PB of data in just over one week of runtime. Frontier-E marks a significant advance in predictive modeling for next-generation cosmological science.

Frontiere, Nicholas [Argonne National Laboratory (↗
DOE OSTI2021

yuwangcn/C4_dynamic_model

A dynamic systems model of C4 photosynthesis was developed based on the previous NADP-ME metabolic model for maize (Wang et al., 2014 ab). The NADP-ME metabolic model is an ordinary differential equation model including all individual steps in C4 photosynthetic carbon metabolism. Here, the model is extended to include posttranslational regulation and temperature response of enzyme activities, dynamic stomata conductance, and leaf energy balance. This model is written in Matlab (R2019a) Steady-state and dynamic gas exchange data for maize (B73), sugarcane (CP88-1762) and sorghum (Tx430) were measured in a greenhouse in Urbana, IL from July 25, 2019 through August 8, 2019. Data include: CO2 response curves Light response curves Photosynthetic induction curves measured in the transition from darkness to high light (1800 μmol m-2 s-1), data logged every 1 min. Photosynthetic induction in the transition from darkness to high light (1800 μmol m-2 s-1) to determine the kinetics of rubisco activation in these C4 crops (τ_Rubisco), data logged every 10 s. Gas exchange under fluctuating light. After dark adaptation, the leaves undergo three light change steps, light intensity was set as 1800 µmol m-2 s-1, 200 µmol m-2 s-1 and 1800 µmol m-2 s-1 for each 1800 s step.

Wang, Yu↗
DOE OSTI2021

AmeriFlux CA-DB2 Delta Burns Bog 2

This is the AmeriFlux version of the carbon flux data for the site CA-DB2 Delta Burns Bog 2. Site Description - The Burns Bog Ecological Conservancy Area (BBECA) in the metro Vancouver area presents a unique opportunity for insight into greenhouse gas dynamics post restoration. The bog was significantly disturbed for peat mining, development and agriculture over the last century, however since 2005 there has been both active restoration and natural regeneration occurring within the bog. Hence areas of bog are at varying stages of ecological recovery. This site is a Pine-Sphagnum low shrub site characteristic of a later successional stage. Conversely, CA-DBB, another AmeriFlux site in BBECA, is a Beakrush-Sphagnum site characteristic of an earlier successional stage.

Knox, Sara↗
DOE OSTI2022

AmeriFlux FLUXNET-1F CA-DB2 Delta Burns Bog 2

This is the AmeriFlux Management Project (AMP) created FLUXNET-1F version of the carbon flux data for the site CA-DB2 Delta Burns Bog 2. This is the FLUXNET version of the carbon flux data for the site CA-DB2 Delta Burns Bog 2 produced by applying the standard ONEFlux (1F) software. Site Description - The Burns Bog Ecological Conservancy Area (BBECA) in the metro Vancouver area presents a unique opportunity for insight into greenhouse gas dynamics post restoration. The bog was significantly disturbed for peat mining, development and agriculture over the last century, however since 2005 there has been both active restoration and natural regeneration occurring within the bog. Hence areas of bog are at varying stages of ecological recovery. This site is a Pine-Sphagnum low shrub site characteristic of a later successional stage. Conversely, CA-DBB, another AmeriFlux site in BBECA, is a Beakrush-Sphagnum site characteristic of an earlier successional stage.

Knox, Sara↗