Search NASA⌕ Search

SEARCH · Search NASA

Results for “PERFECT GAS”

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 19 records

Thermodynamic properties of “near-perfect” gas

When applied thermodynamics requires simple representations of thermodynamic quantities, one often finds quantities such as pressure and the thermal energy density expressed as a monomial linear product of density and temperature power laws. This is a simple generalization of the product of density and temperature used for the perfect gas that admits a broader range of thermodynamic behavior into analytic fluid calculations or provides an analytic form that is readily fit to tabulated equation-of-state data for an arbitrary material. This paper reviews the thermodynamic properties of this generalized perfect-gas model, treating it as a class of “near-perfect-gas” models that are defined and unified here, for the first time, through a shared Helmholtz free energy. This generalization from perfect to near-perfect-gas models preserves some perfect-gas properties (e.g., constant adiabatic exponents) but not others (e.g., constant specific heats) and is constrained by thermodynamic consistency. It is important to be aware of the properties of this class of near-perfect-gas models when using them as surrogates for real nonperfect materials that may have properties that cannot be captured by near-perfect-gas models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Exact Solution for a Magnetogasdynamical Cylindrical Shock Wave in a Self-Gravitating Rotating Perfect Gas with Radiation Heat Flux and Variable Density

An exact similarity solution for a magnetoradiative cylindrical shock wave in a self-gravitating rotating perfect gas is obtained. The density, azimuthal velocity, and magnetic field strength are assumed to vary in an undisturbed medium. It is shown that the flow variables, namely, the radial velocity, pressure, magnetic field strength, azimuthal velocity, mass, and the radiation flux, decrease from the highest values at the shock front to zero; however, the density tends to infinity as the symmetry axis is approached. The effects of variation in the magnetic field strength, gravitational parameter, rotational parameter, and in the adiabatic exponent on the flow variables and shock strength are discussed. The solutions obtained for self-gravitating and nongravitating media are compared. The total energy of the shock wave is shown to be not constant.

42 ENGINEERING↗

Outflow Boundary Conditions for Turbine-Integrated Rotating Detonation Combustors

This study examines outflow boundary conditions (BCs) in computational fluid dynamics (CFD) simulations of a transition duct with and without guide vanes that converts supersonic flow exiting a rotating detonation combustor (RDC) to subsonic flow to drive a turbine. Since the flow exiting the transition duct has swirling shock waves with significant spatial and temporal variations in pressure, temperature, and Mach number, imposing proper BCs poses a challenge. To ensure all swirling shock waves exit the transition duct without creating non-physical reflected waves at its outlet, this study examined three outflow BCs: (1) the average pressure imposed at the duct’s outlet, (2) a nonreflecting BC (NRBC) with a specified average pressure imposed at the duct’s outlet, (3) the average pressure imposed at the outlet of an extension duct made up of a buffer layer and a sponge layer. This study is based on the three-dimensional, unsteady density-weighted-ensemble-averaged continuity, Navier–Stokes, and energy equations for a thermally perfect gas closed by the realizable k–ε model and “enhanced” wall functions. The results obtained show that imposing an average pressure at the transition duct’s outlet produces spurious waves that degrade the physical meaningfulness of the solution. When the NRBC was applied, swirling shock waves exited the duct’s outlet without creating spurious waves. However, its usage requires the gas to be thermally, as well as calorically, perfect, which this study shows could be a concern. By imposing the average pressure at the outlet of an extension duct, the gas does not need to be calorically perfect. The results obtained show the effects of the sponge layer’s length and coarsening ratio on damping nonuniformities in non-physical reflected waves to ensure the flow exiting the transition duct’s outlet can do so as if there are no boundaries present and has the desired average pressure—even though the BC is applied at the extension duct’s outlet.

gas turbines↗

Dalton’s and Amagat’s laws fail in gas mixtures with shock propagation

A shock propagating through a gas mixture leads to pressure, temperature, and density increases across the shock front. Rankine-Hugoniot relations correlating pre- and post-shock quantities describe a calorically perfect gas but deliver a good approximation for real gases, provided the pre-shock conditions are well characterized with a thermodynamic mixing model. Two classic thermodynamic models of gas mixtures are Dalton’s law of partial pressures and Amagat’s law of partial volumes. We measure post-shock temperature and pressure in experiments with nonreacting binary mixtures of sulfur hexafluoride and helium (two dramatically disparate gases) and show that neither model can accurately predict the observed values, on time scales much longer than that of the shock front passage, due to the models’ implicit assumptions about mixture behavior on the molecular level. However, kinetic molecular theory can help account for the discrepancy. Our results provide starting points for future theoretical work, experiments, and code validation.

42 ENGINEERING↗

A scalable compressible volume of fluid solver using a stratified flow model

Here, in this paper, a Volume of Fluid method which uses a Stratified Flow model for flux calculation is proposed to solve a compressible multiphase flow in a way that has high parallel efficiency. The solver is part of an open-source computing tool to study fuel entrainment and combustion in the application of paraffin-based hybrid rocket motors. This paper focuses on the compressible Volume of Fluid solver for the liquid and gaseous phases in this environment. This solver uses perfect gas and stiffened gas models for equation of state, with the ability to easily add other models as desired. Additionally a comparison of a Riemann solver versus an AUSM+up scheme is conducted, showing improvement with the later in the solution of multiphase flows. This solver successfully reproduces shock tube cases with the expected accuracy, as well as simulating two-dimensional shear and gravity-driven flows. Results also show that the code is able to utilize a Stratified Flow model to evaluate shear flow without the need for interface reconstruction or gradient calculations. Finally, it is shown that the solver has near-ideal scaling under strong scaling tests as well as good performance in static scaling, giving improved performance over current options and implying future advancement of high-performance computing efficiency for multiphase flow solvers.

97 MATHEMATICS AND COMPUTING↗

A finite scale model for shock structure

In this paper we explore the use of a finite scale model of fluid dynamics to predict the finite structure of a shock wave in a perfect gas. We begin by documenting the history and issues that have arisen when Navier–Stokes theory is applied to the shock structure problem, and continue by motivating the improvement that finite scale theory might provide from its representation of inviscid (anomalous) dissipation. Furthermore, our primary results include the formulation of a traveling wave equation from finite scale theory, an analysis of the solutions of that equation as regards shock width and monotonicity properties, and an estimation of the critical parameter of the theory from experimental data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Finite scale theory: Predicting nature’s shocks

In this paper, we describe a continuum model that accurately reproduces the experimentally measured structure of physical shocks in a perfect gas. In this work, we begin by presenting a history of shock structure research, theoretical, experimental and numerical, to quantify the significant discrepancies between Navier–Stokes predictions and laboratory measurements. In our first main result, we discuss modifications that generalize the Chapman–Enskog approximation and lead to a new continuum model termed the Finite Scale Equations. In our second main result, we use this continuum model to calculate shock structure with results comparable in accuracy to numerical simulations based on the Boltzmann equation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Direct numerical simulation of supercritical oxy-methane mixing layers with CO 2 substituted counterparts

Direct numerical simulations (DNS) of temporally developing, three dimensional, CH 4 /CO 2 , CH 4 /O 2 and CO 2 /O 2 mixing layers, are conducted at a supercritical pressure of 300 atm. To effectively model the supercritical regime, the employed formulation includes the compressible form of the governing equations, the cubic Peng-Robinson equation of state and a generalized formulation for heat and mass flux vectors derived from non-equilibrium thermodynamics and fluctuation theory. A linear inviscid stability analysis is also performed for each case, to determine its most unstable wavelength. Flow visualizations reveal the presence of high density gradient magnitude regions for all three mixing layers, with conditional averages indicating increased presence of heavier fluid species within these regions. No significant departures are observed from perfect gas behavior, with compressibility factors very close to unity for all three mixing cases. Applicability of presumed probability density function methods (PDF) is examined for the three supercritical mixing layers. An a priori analysis is also conducted to investigate various simplifying assumptions employed in modeling various subgrid scale (SGS) flux models. Two additional terms are identified in the large eddy simulations (LES) equations, the gradient of SGS contribution of pressure in the momentum equation and the gradient of SGS contribution of heat flux in energy equation, whose magnitudes are similar and comparable with their respective resolved terms. The performance of the scale similarity model to represent these additional terms is investigated. Lastly, the performance of Smagorinsky, gradient and scale similarity models are also investigated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Efficient GPU Implementation of Automatic Differentiation for Computational Fluid Dynamics

Many scientific and engineering applications require repeated calculations of derivatives of output functions with respect to input parameters. Automatic Differentiation (AD) is a method that automates derivative calculations and can significantly speed up code development. In Computational Fluid Dynamics (CFD), derivatives of flux functions with respect to state variables (Jacobian) are needed for efficient solutions of the nonlinear governing equations. AD of flux functions on graphics processing units (GPUs) is challenging as flux computations involve many intermediate variables that create high register pressure and require significant memory traffic because of the need to store the derivatives. This paper presents a forward-mode AD method based on multivariate dual numbers that addresses these challenges and simultaneously reduces the floating-point operation count. The dimension of the multivariate dual numbers is optimized for performance. The flux computations are restructured to minimize the number of temporary variables and reduce register pressure. For effective utilization of memory bandwidth, shared memory is used to store the local flux Jacobian. This AD implementation is compared with several other Jacobian implementations on an NVIDIA V100 GPU (V100). For three-dimensional perfect-gas compressible-flow equations implemented in a practical CFD code, the AD implementation of a flux Jacobian based on multivariate dual numbers of dimension 5 outperforms all other GPU AD implementations on V100. Its performance is comparable with the optimized hand-differentiated version. Finally, the implementation achieves 75% of the peak floating-point throughput and 61 % of the peak global device memory bandwidth usage.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A Polytropic Approximation of Compressible Flow in Pipes With Friction

This paper demonstrates the usefulness of treating subsonic Fanno flow (adiabatic flow, with friction, of a perfect gas in a constant-area pipe) as a polytropic process. It is shown that the polytropic model allows an explicit equation for mass flow rate to be developed. The concept of the energy transfer ratio is used to develop a close approximation to the polytropic index. Explicit equations for mass flow rate and net expansion factor in terms of upstream properties and pressure ratio are developed for Fanno and isothermal flows. An approximation for choked flow is also presented. Here, the deviation of the results of this polytropic approximation from the values obtained from a traditional gas dynamics analysis of subsonic Fanno flow is quantified and discussed, and a typical design engineering problem is analyzed using the new method.

42 ENGINEERING↗

Modern infinitesimals and the entropy jump across an inviscid shock wave

Here, this article applies nonstandard analysis to study the generalized solutions of entropy and energy across one-dimensional shock waves in a compressible, inviscid, perfect gas. Nonstandard analysis is an area of modern mathematics that studies number systems that contain both infinitely large and infinitely small numbers. For an inviscid shock wave, it is assumed that the shock thickness occurs on an infinitesimal interval and that the jump functions for the field variables are smoothly defined on this interval. A weak converse to the existence of the entropy peak is derived and discussed. Generalized solutions of the Euler equations for entropy and energy are then derived for both theoretical and realistic normalized velocity profiles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Symmetry Determining Equations of the Rankine-Hugoniot Equations for Variable Velocity Shock Waves

The “constant velocity piston” problem (Fig. 1), also known as the “piston problem,” is a standard model for a one dimensional, in our case linear, symmetric shock wave moving through an inviscid, perfect gas. The model can be divided into two regions - a perturbed section on the left and an unperturbed section on the right - by a moving shock wave moving left to right. Both the perturbed and unperturbed sections, i.e. the shocked and unshocked regions, respectively, obey the Eulerian conservation equations; however, at the exact location of the shock, there is a mathematical discontinuity not satisfied by the Euler equations. To ensure continuity and conservation of certain quantities when crossing between the unshocked and shocked regions, we evoke a series of equations derived from the Eulerian conservation equations, called the Rankine-Hugoniot equations, or “jump” equations as it is often referred to in the literature on the topic. The classical constant-velocity piston problem assumes the piston features a constant driving velocity (among many other willing suspensions of belief required in the pursuit of a first principles equation model); consequent to this assumption is a constant-velocity shock and a constant-velocity shocked flow state. However, using Lie Group Theory (LGT), also known as symmetry analysis, we can attempt to reinterpret the model with a shock wave of variable velocity in time and space. An extension of the model in this way opens up the possibility for obtaining new analytical solutions to the piston problem for certain shock velocity models. In this report, we use LGT to derive the symmetry determining equations (SDEs), whose solutions are Lie groups, which permit analytical solutions. In the future, we can then use the SDEs to define constraint equations on the shock velocity model and what the successive solutions to the Euler equations might be based off such constraints. This report is structured as follows: Section 2 provides a brief derivation of the Rankine-Hugoniot (“jump”) equations; Section 3 gives an overview of Lie group theory; Section 4 derives the SDEs of the jump equations; Section 5 derives the Euler conservation equations for fluids; and Section 6 presents concluding remarks and opportunities for future studies.

42 ENGINEERING↗

Contributions to the mechanistic understanding of the microstructural evolution in irradiated U-Mo dispersion fuel

Here, advanced microstructural characterization techniques, such as scanning electron microscopy (SEM) and scanning transmission electron microscopy - energy dispersive x-ray spectroscopy (STEM-EDS), were used to interpret the fuel microstructure evolution and fission products behavior in U-Mo dispersion fuel irradiated in the Advanced Test Reactor (ATR) as part of the European Mini-Plate Irradiation Experiment (EMPIrE) test. The larger as-fabricated fuel grain size achieved by heat-treating the U-Mo powder resulted in slower high burnup structure (HBS) development and reduced fission gas porosity. Slower HBS kinetics was observed at the fuel kernels’ periphery, which contained smaller and less fission gas bubbles at all fission densities (FDs) investigated and was attributed to a locally reduced damage density and fission products concentration, as corroborated with Monte Carlo simulations. The non-refined grains at the fuel kernel periphery hosted a perfectly ordered fission Gas Bubble Superlattice (GBS) up to 6.3 × 10 21 fissions/cm 3 . Nano-scale STEM-EDS analysis presented in this study provided useful information on the GBS characteristic morphology and evolution in U-Mo fuel. The concentration of fission gas in the GBS progressively increased with FD, pointing to an evolution of the nanobubble pressure status with irradiation. A possible connection between the GBS collapse and HBS onset is proposed for which there exists a threshold in the misorientation of the refined sub-grains above which the GBS stability during irradiation is no longer preserved, resulting in the GBS collapse.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Uncovering dark matter density profiles in dwarf galaxies with graph neural networks

Dwarf galaxies are small, dark-matter-dominated galaxies, some of which are embedded within the Milky Way. Their lack of baryonic matter (e.g., stars and gas) makes them perfect test beds for probing the properties of dark matter—understanding the spatial dark matter distribution in these systems can be used to constrain microphysical dark matter interactions that influence the formation and evolution of structures in our Universe. Here, we introduce a new method that leverages simulation-based inference and graph-based machine learning in order to infer the dark matter density profiles of dwarf galaxies from observable kinematics of stars gravitationally bound to these systems. Our approach aims to address some of the limitations of established methods based on dynamical Jeans modeling. We show that this novel method can place stronger constraints on dark matter profiles and, consequently, has the potential to weigh in on some of the ongoing puzzles associated with the small-scale structure of dark matter halos, such as the core-cusp discrepancy.

79 ASTRONOMY AND ASTROPHYSICS↗

Investigation of high-pressure turbulent jets using direct numerical simulation

Direct numerical simulations of free round jets at a Reynolds number (Re D ) of 5000, based on jet diameter (D) and jet-exit bulk velocity (U e ), are performed to study jet turbulence characteristics at supercritical pressures. The jet consists of nitrogen (N 2 ) that is injected into N 2 at the same temperature. To understand turbulent mixing, a passive scalar is transported with the flow at unity Schmidt number. Two sets of inflow conditions that model jets issuing from either a smooth contraction nozzle (laminar inflow) or a long pipe nozzle (turbulent inflow) are considered. By changing one parameter at a time, the simulations examine the jet-flow sensitivity to the thermodynamic condition (characterized in terms of the compressibility factor ($Z$) and the normalized isothermal compressibility), inflow condition and ambient pressure (p ∞ ) spanning perfect- to real-gas conditions. The inflow affects flow statistics in the near field (containing the potential core closure and the transition region) as well as further downstream (containing fully developed flow with self-similar statistics) at both atmospheric and supercritical p ∞ . The sensitivity to inflow is larger in the transition region, where the laminar-inflow jets exhibit dominant coherent structures that produce higher mean strain rates and higher turbulent kinetic energy than in turbulent-inflow jets. Decreasing $Z$ at a fixed supercritical p ∞ enhances pressure and density fluctuations (non-dimensionalized by local mean pressure and density, respectively), but the effect on velocity fluctuations depends also on the local flow dynamics. When $Z$ is reduced, large mean strain rates in the transition region of laminar-inflow jets significantly enhance velocity fluctuations (non-dimensionalized by local mean velocity) and scalar mixing, whereas the effects are minimal in jets from turbulent inflow.

42 ENGINEERING↗

Non-cratonic Diamonds from UHP Metamorphic Terranes, Ophiolites and Volcanic Sources

The discovery of diamonds in metamorphic rocks of continental affinities occurred shortly after the discovery of coesite in similar rocks. These important discoveries led to a revolution in our understanding of the subduction and exhumation of continental materials and the establishment of a new discipline, ultra-high-pressure metamorphism (UHPM). After these discoveries more ultra-high pressure (UHP) minerals were found and dozens of UHPM terranes were established. Diamond is a remarkable material which is chemically inert and stable over geological timescales making it the perfect “geological container” where gas, fluid, and solid inclusions can be preserved. Moreover, its presence is indicative of specific pressure and temperature conditions which implie subduction to a minimum depth of ~120 km. The inclusions trapped inside these diamonds can shed light on the composition, redox state, and evolution of the fluids related to UHPM diamond formation while the carbon isotope ratios of the diamond itself can inform us about the source of carbon. The purpose of this review is to provide a brief history on the discovery of microdiamonds starting with work that began in the mid-1960s in Kokchetav massif, Kazakhstan, through to the present day. Particular attention will be on more recent micro-diamond discoveries in the last decade. The current state of understating of the mechanisms of UPHM diamond formation, the misidentification of microdiamonds due to contamination from sample extraction/preparation, and the future of UHPM diamond research are discussed. The paper also considers the controversial topic of the occurrence of lonsdaleite within the population of microdiamonds from UHPM terranes as well as recently reported ophiolite- and volcanic-hosted diamonds.

58 GEOSCIENCES↗

An Uncertainty Management Framework for Integrated Gas-Electric Energy Systems

In many parts of the world, electric power systems have seen a significant shift toward generation from renewable energy and natural gas. Because of their ability to flexibly adjust power generation in real time, gas-fired power plants are frequently seen as the perfect partner for variable renewable generation. However, this reliance on gas generation increases interdependence and propagates uncertainty between power grids and gas pipelines and brings coordination and uncertainty management challenges. To address these issues, we propose an uncertainty management framework for uncertain, but bounded gas consumption by gas-fired power plants. The admissible ranges are computed based on a joint optimization problem for the combined gas and electricity networks, which involves chance-constrained scheduling for the electric grid and a novel robust optimization formulation for the natural-gas network. This formulation ensures feasibility of the integrated system with a high probability, while providing a tractable numerical formulation. A key advance with respect to existing methods is that our method is based on a physically accurate, validated model for transient gas pipeline flows. Our case study benchmarks our proposed formulation against methods that ignore how reserve activation impacts the fuel use of gas power plants and only consider predetermined gas consumption. Here, the results demonstrate the importance of considering uncertainty to avoid operating constraint violations and curtailment of gas to the generators.

24 POWER TRANSMISSION AND DISTRIBUTION↗

SWIFT - Shaped-Waveform Interrogation of Fluids Technique

Researchers at Los Alamos have developed the Shaped-Waveform Interrogation of Fluids Technique (SWIFT), a technology enabling unprecedented process control capability. This technology provides important composition, flow rate, and pipe thickness measurements in real-time, at an extremely low cost and maintenance cycle. By utilizing a specialized shaped-waveform, as well as state-of-the-art signal processing, this information can be gathered without any frequency-domain signal processing, allowing for simplified electronics compared to nearly any other meter on the market. Since this technology is in expensive and non-invasive, with the capability of measuring liquids in both flow pipes and sealed containers, it is perfect for several industrial applications in the oil and gas, beverage, chemical, and pharmaceutical spaces. We are seeking a commercialization partner to license the technology or fund us through a Cooperative Research and Development Agreement (CRADA) to a further develop a usable product for commercial purposes.

47 OTHER INSTRUMENTATION↗