Search NASA⌕ Search

SEARCH · Search NASA

Results for “Compressible flows”

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

Numerical analysis of air dehumidification through electrospray-enhanced vortical flow cyclone separator

The dehumidification of air is a highly energy-intensive process, typically requiring vapor compression systems to cool the air significantly below its dew point temperature. Here, this study numerically investigates air dehumidification achieved by electrostatic and dielectrophoresis forces in a converging duct and centrifugal principles in cyclone separators. An electrospray was employed to energize water droplets, which then captured water vapor via dielectrophoresis. A novelty of the present study is the inclusion of a cyclone separator positioned downstream of the electrospray. This combination increased the dehumidification efficiency and controlled the flow of the condensate. By integrating electrospray technology and applying an electric field, this study examines the impact of particle separation efficiency, pressure drops, and cyclone performance on air dehumidification. The control parameters for air dehumidification were investigated. Three cyclones were modelled for airflow rates of 5 cubic feet per minute [cfm] (0.0024 m 3 /s), 38 cfm (0.0179 m 3 /s), and 200 cfm (0.094 m 3 /s). The results indicated air dehumidification was up to 5.7 % under environmental conditions, similar to those used to validate the CFD model. The degree of air dehumidification was primarily influenced by the electric charge acquired by the water droplets, their flight trajectory, and the vorticity intensity within the cyclone. Air dehumidification could increase by up to 13 % by increasing the number of droplets injected into the air stream, the flight time, and the electrical charge. The breakdown electric field strength threshold near the electrospray was a limiting factor for further enhancing dehumidification performance. This study examines the optimization of an electro-assisted dehumidification system that includes cyclone-based droplet separation through numerical analysis. The results provide insights into the design of energy-efficient air treatment technologies.

42 ENGINEERING↗
DOE OSTI2024

Concurrent measurement of strain and chemical reaction rates in a calcite grain pack undergoing pressure solution: Evidence for surface-reaction controlled dissolution

Pressure solution is inferred to be a significant contributor to sediment compaction and lithification, especially in carbonate sediments. For a sediment deforming primarily by pressure solution, the compaction rate should be directly related to the rate of calcite dissolution, transport along grain contacts, and calcite reprecipitation. Previous experimental work has shown that there is evidence that deformation in wet calcite grain packs is consistent with control by pressure solution, but considerable ambiguity remains regarding the rate limiting mechanism. We present the results of laboratory compaction experiments designed to directly measure calcite dissolution and precipitation rates (recrystallization rates) concurrently with strain rate to test whether measured rates are consistent with predicted rates both in absolute magnitude and time evolution. Recrystallization rates are measured using trace element chemistry (Sr/Ca, Mg/Ca) and isotopes (87Sr/86Sr) of fluids flowing slowly through a compacting grain pack as it is being triaxially compressed. Imaging techniques are used to characterize the grain contacts and strain effects in the post-experiment grain pack. Our data show that calcite recrystallization rates calculated from all three geochemical parameters are in approximate agreement and that the rates closely track strain rate. The geochemically inferred rates are close to predicted rates in absolute magnitude. Uncertainty in grain contact dimensions makes distinguishing between surface reaction control and diffusion control difficult. Measured reaction rates decrease faster than predicted from standard pressure solution creep flow laws. This inconsistency may indicate that calcite dissolution rates at grain contacts are more complex, and more time-dependent, than suggested by geometric models designed to predict grain contact stresses.

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

Numerical Solution of the Steady-State Network Flow Equations for a Non-Ideal Gas

Herein we formulate a steady-state network flow problem for non-ideal gas that relates injection rates and nodal pressures in the network to flows in pipes. For this problem, we present and prove a theorem on uniqueness of generalized solution for a broad class of non-ideal pressure-density relations that satisfy a monotonicity property. Further, we develop a Newton-Raphson algorithm for numerical solution of the steady-state problem, which is made possible by a systematic non-dimensionalization of the equations. The developed algorithm has been extensively tested on benchmark instances and shown to converge robustly to a generalized solution. Previous results [1]-[4], indicate that the steady-state network flow equations for an ideal gas are difficult to solve by the Newton-Raphson method because of its extreme sensitivity to the initial guess. In contrast, we find that non-dimensionalization of the steady-state problem is key to robust convergence of the Newton-Raphson method. We identify criteria based on the uniqueness of solutions under which the existence of a non-physical generalized solution found by a non-linear solver implies non-existence of a physical solution, i.e., infeasibility of the problem. Finally, we compare pressure and flow solutions based on ideal and non-ideal equations of state to demonstrate the need to apply the latter in practice. The solver developed in this article is open-source and is made available for both the academic and research communities as well as the industry.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

StOKeDMD: Streaming Occupation kernel dynamic mode decomposition

Dynamic mode decomposition (DMD) has become a common technique for constructing surrogate models for dynamical systems from observed system states. The Occupation Kernel DMD (OKDMD) method proposed in (Rosenfeld et al., 2022) and (Rosenfeld et al., 2024) is a Liouville operator based method that builds surrogate models from system state trajectories. Here, this paper proposes an extension of OKDMD to the case when the system states are observed in a streaming fashion, i.e., only a small fraction of the state trajectory is available at a given time. The developed method, Streaming Occupation Kernel DMD (StOKeDMD), accommodates the streaming data input by leveraging properties of specific choices of kernel functions and occupation kernels. We apply the StoKeDMD method as a compression method for streaming data, analyze the memory complexity, and demonstrate the performance of StoKeDMD in the compression of streaming data generated from a Lorenz system and a fluid flow simulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Characterization of Biomass-Based Fuels and Fuel Blends for Low-Emissions, Advanced Compression Ignition Engines (Final Report)

The goal of the Co-Optima program was to accelerate deployment of co-optimized fuels and engines that will reduce fuel consumption, and criterion pollutants and greenhouse gas (GHG) emissions for both spark and compression-ignited engines. In this project, combustion measurements from a constant pressure flow rig (CPFR) with optical access were used to assess fundamental combustion behavior of number of different candidate fuels and fuel blends. This chamber and its peripheral systems allowed stable control of ambient air pressure and temperature and fuel injection pressure and timing for hundreds of repeat injections allowing statistical interpretation of the results. The optical diagnostics enabled capturing of liquid length, vapor length, ignition timing, lift-off length as well as localized (temporal and spatially resolved) mixture refractive index, flame structure, and soot mass. The operating conditions allowed the study of combustion conditions promoting lean lifted flame combustion (LLFC), a mixing-controlled low-temperature combustion (LTC) strategy to reduce soot production in Advanced Compression Ignition (ACI) engines.

09 BIOMASS FUELS↗
DOE OSTI2021

Constraint and size effects in confined layer plasticity

Plane strain finite element analyses are used to model the experiments for a thin metal layer confined between elastic solids. Here, the thin metal layer undergoes elastic-plastic deformations, here modelled by strain gradient plasticity, while the rest of the solid deforms only elastically. Plane strain is assumed and finite strains are accounted for. The microscopic boundary condition at the band interface is initially taken to be zero plastic strain increment, but plastic straining at the interface is allowed when the plastic strain gradient at the interface has exceeded a critical value. When the narrow metal band is inclined 45° relative to the direction of compression of the solid simple shear develops in the band and predictions are rather similar to those of a one-dimensional analysis for an infinite band. For a band at 90° the elastic parts give significant constraint on plastic flow, the stress fields are non-uniform along the band and the average nominal compressive stress shows strong sensitivity to the layer thickness and to the width of the specimen. Remarkably good qualitative, and in some cases quantitative, agreement is found with the experimental observations. Considering also other band angles, both the average nominal shear and compressive stress show variations with the angle.

36 MATERIALS SCIENCE↗
DOE OSTI2023

MOOSE Navier–Stokes module

The MOOSE Navier–Stokes module solves mass, momentum, energy, and passive scalar conservation equations in the context of fluid flow. The module supports solution of these equations in both free flow and porous medium contexts and for a range of fluid compressibility. The conservation equations can be discretized in space using continuous Galerkin finite elements or with cell centered finite volumes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Volumetric efficiency and pressure pulsations of a triplex reciprocating supercritical CO 2 pump

A reciprocating triplex pump has been integrated into a closed-loop supercritical carbon dioxide test facility. The pump is designed to compress liquid or supercritical carbon dioxide and has a maximum volumetric flow rate of 51 L/min. A sweep across the pump’s speed range of 100–1000 rev/min was performed with inlet conditions between 60 and 160 bar. The volumetric efficiency was found to increase as pump speed increased from 10 to 40 % before it plateaued at higher pump speeds. The mechanical efficiency increased linearly with pump speed and was not affected by system pressure. Pressure pulsations caused by the plungers and valves were also measured at each speed. Additionally, a 35-liter buffer volume, connected near the pump’s discharge was isolated at each speed to determine its effectiveness as a pulsation dampener. Overall the reciprocating pump was shown to perform exceptionally with liquid or supercritical inlet conditions applicable to many supercritical fluid applications.

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

Dark Energy Survey Year 3 results: Simulation-based 𝑤CDM inference from weak lensing and galaxy clustering maps with deep learning: Analysis design

Data-driven approaches using deep learning are emerging as powerful techniques to extract non-Gaussian information from cosmological large-scale structure. Here, this work presents the first simulation-based inference (SBI) pipeline that combines weak lensing and galaxy clustering maps in a realistic Dark Energy Survey Year 3 (DES Y3) configuration and serves as preparation for a forthcoming analysis of the survey data. We develop a scalable forward model based on the CosmoGridV1 suite of N-body simulations to generate over one million self-consistent mock realizations of DES Y3 at the map level. Leveraging this large dataset, we train deep graph convolutional neural networks on the full survey footprint in spherical geometry to learn low-dimensional features that approximately maximize mutual information with target parameters. These learned compressions enable neural density estimation of the implicit likelihood via normalizing flows in a ten-dimensional parameter space spanning cosmological 𝑤CDM, intrinsic alignment, and linear galaxy bias parameters, while marginalizing over baryonic, photometric redshift, and shear bias nuisances. To ensure robustness, we extensively validate our inference pipeline using synthetic observations derived from both systematic contaminations in our forward model and independent Buzzard galaxy catalogs. Our forecasts yield significant improvements in cosmological parameter constraints, achieving 2−3× higher figures of merit in the 𝛺 𝑚 − 𝑆 8 plane relative to our implementation of baseline two-point statistics and effectively breaking parameter degeneracies through probe combination. These results demonstrate the potential of SBI analyses powered by deep learning for upcoming Stage-IV wide-field imaging surveys.

Thomsen, A. [Zurich, ETH] (ORCID:0000000203099021)↗
DOE OSTI2021

Strongly Anisotropic Thermomechanical Response to Shock Wave Loading in Oriented Samples of the Triclinic Molecular Crystal 1,3,5-Triamino-2,4,6-trinitrobenzene

All-atom molecular dynamics (MD) simulations were used to study shock wave loading in oriented single crystals of the highly anisotropic triclinic molecular crystal 1,3,5-triamino-2,4,6-trinitrobenzene (TATB). The crystal structure consists of planar hydrogen-bonded sheets of individually planar TATB molecules that stack into graphitic-like layers. Shocks were studied for seven systematically prepared crystal orientations with limiting cases that correspond to shock propagation exactly perpendicular and exactly parallel to the graphitic-like layers. The simulations were performed for initially defect-free crystals using a reverse-ballistic configuration that generates explicit, supported shocks. Final longitudinal stress components are between ≈8.5 and ≈10.5 GPa for the 1.0 km s –1 impact speed studied. Orientation-dependent properties are reported including shock speeds, stresses, temperatures, compression ratios, and local material strain rates. Spatiotemporal maps of the temperature, stress tensor, material flow, and molecular orientations reveal complicated processes that arise for specific shock directions. Furthermore, the results indicate that TATB shock response is highly sensitive to crystal orientation, with significant qualitative differences for the time evolution of the stress tensor and temperature, elastic/inelastic compression response, defect formation and growth, critical von Mises stress, and strain rates during shock rise that span nearly an order of magnitude. A variety of inelastic deformation mechanisms are identified, ranging from crumpling of graphitic-like layers to dislocation-mediated plasticity to intense shear strain localization. To our knowledge, these are the first systematic MD simulations and analysis of explicit shock wave propagation along nontrivial crystal directions in a triclinic molecular crystal.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Probing turbulence with lasers and coherent light sources

Fluid and plasma turbulence consists of multi-scale interactions where the energy of the large scales is dissipated through a cascade to the smallest scales set by the microphysics. Turbulence and mixing can be detrimental to system performance, such as in inertial confinement fusion (ICF). Our current understanding of turbulence is mostly derived from studies of incompressible flows, where the intermediate-scaled eddies in the “inertial range” bridge the transfer of energy from the large eddies to the dissipative range, and that there is no communication between eddies of disparate scales. Whether a similar cascade process exists for compressible turbulence, however, remains to be seen; experiments have been limited by diagnostics and supersonic flow facilities. This project has focused on the development of techniques to measure the evolution of flows and microphysical transport properties using particle tracking velocimetry. This project has investigated three main areas: (1) particle-tracking experiments with x-ray radiography, (2) shock-particle interactions, (3) quantifying viscous effects on particles in shocked flows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

SPH Modeling of Biomass Granular Flow: Engineering Application in Hoppers and Augers

Numerical modeling of granular biomass material flow in handling operations is indispensable to decipher flow upsets, commonly manifested as clogging and jamming in hoppers and augers. With a computational tool developed based on smoothed particle hydrodynamics (SPH), we simulated the hopper flow and auger feeding of six granular biomass materials. The good agreement between the experimental and numerical flow rates demonstrated the capability of the developed SPH solver in modeling the complicated flow of granular biomass materials. Further, the impact of physical and numerical parameters is investigated, and the major results show that the hopper flow pattern is controlled by shear band evolution; hopper clogging is collectively influenced by the opening size, wall friction, and material packing, and the impact of material compressibility on the auger feed rate is minimal. These parametric studies validate the solver’s robustness in simulating biomass flow in handling equipment and demonstrate that the SPH computational tool can provide insights about granular flow mechanics to facilitate handling equipment design and optimize handling operations.

09 BIOMASS FUELS↗
DOE OSTI2022

Modeling atomically mixed graded density impactors

Graded density impactors (GDIs) are multi-material composite impactors used in gas gun experiments to tailor the drive conditions imparted to a sample test material. Previous graded density impactors generally rely on thin, but discrete, layers of different materials. The thinner and the greater number of layers will result in smoother compression. Taken to the limit of very thin layers would be pure material 1 at one surface, such as the front surface of an impactor, smoothly transitioning at the atomic scale to pure material 2 on the back surface. Such an impactor can initially shock, then smoothly compresses a material during a dynamic experiment. This type of experiment can serve to explore a larger region of thermodynamic space than a single or even multi-shock experiments. An overview of how graded density impactors are made is reviewed and sample results are given. Additionally, a strategy for modeling these kinds of impactors is presented. The length scales of constituent mixing are given from the experimental build through electrochemical-deposition. Equation of state models for pure constituents and their subsequent mixtures are presented. It is demonstrated that the time scales for pressure and temperature equilibration, for atomically mixed GDIs, are short enough to be a justifiable closure for the resulting multiphase flow. Furthermore, we present simulation results of dynamic shock followed by a ramp compression, utilizing a silver/gold graded density impactor, onto a tantalum sample.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Elastic Flow Modeling for Hydropower Digital Twins

This report details the elastic, unsteady, one-dimensional flow equations that are used to model flow through a penstock in a hydropower facility. The elastic flow model is accurate even in cases of fast transients, long penstock length, and high gravitation head. The compressibility of the water and elasticity of the pipe walls are explicitly accounted for in the developed models so that the water hammer phenomena can be accurately captured. The elastic flow model is coupled to a mechanistic turbine model to resolve the dynamic feedbacks between the elastic water column and the turbine rotation rate. A finite volume method is employed to solve the governing equations, and the method is shown to have low numerical diffusion in sharp gradient phenomena, such as those encountered when simulating water hammer. The method is applied to single dimensional linear advection benchmark problem and numerical results are compared with analytic results. This is followed on by an application to a full hydropower system with a coupled turbine. The elastic flow model results are compared against an inelastic model.

13 HYDRO ENERGY↗
DOE OSTI2021

Improvement of Numerical Methods in Pronghorn

During the fiscal year of 2021, the finite-volume method (FVM) was deployed in Pronghorn to provide improved efficiency, stability, and accuracy for coarse-mesh, thermal-hydraulics problems. While the main goals of fiscal year 2021 were met, several issues emerged from the early deployment of the finite-volume method in Pronghorn. These issues were: The compressible and incompressible formulations that were implemented are inadequate for many nuclear reactor flow problems; Omission of terms accounting for the porosity and Darcy-Forchheimer body force discontinuities in the Rhie-Chow interpolation lead to oscillations in pressure and velocity at these discontinuities; The FVM lacks a correction for non-orthogonal grids for computing accurate pressure gradient leading to loss of accuracy in regions with skewed elements; The FVM currently uses a monolithic solver. Monolithic solvers have issues dealing with the saddle-point nature of the discretized fluid equations. The result are bad convergence if direct factorization is not used and large memory consumption when direct factorization is used. These four issues are addressed in this report. In particular, we report the completion of the following task: Implementation of a weakly compressible formulation in the MOOSE Navier-Stokes module; Implementation of Moukalled’s method for including body forces in the RCI. Additionally, we identified the need to smooth the porosity using Moukalled’s face-cell smoothing operator; Implementation of a non-orthogonal correction for the Green-Gauss gradient computation; Preliminary implementation of a SIMPLE segregated solver.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Effect of PBX prill size on the compressive properties of PBX 9502

Plastic-Bonded Explosive (PBX) prills are agglomerates a few millimeters in diameter composed of High Explosive crystals non-uniformly distributed in and surrounded by a polymer binder. When pressed at elevated temperatures at or above the binder melt temperature the polymer flows to uniformly coat the HE crystals to create a well-consolidated compaction. In reality it does so imperfectly, such that x-ray tomographic scans often look like a collection of prills mashed together. It is speculated that the larger the prill, the more non-uniform the binder distribution, but this has never been proven. Even if that were not so, the larger the prill the farther binder must flow toward the middle to homogenize during pressing. Thus, the degree of homogenization depends in part on prill size. The degree to which binder flows during pressing in turn affects the void distribution within pressed charges, which one suspects will affect shock sensitivity and material strength. In this paper we explore, for ~30 PBX 9502 formulation batches using the same TATB powder lot, how prill size affects the compressive mechanical properties.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Understanding latent timescales in neural ordinary differential equation models of advection-dominated dynamical systems

The neural ordinary differential equation (ODE) framework has shown considerable promise in recent years in developing highly accelerated surrogate models for complex physical systems characterized by partial differential equations (PDEs). For PDE-based systems, state-of-the-art neural ODE strategies leverage a two-step procedure to achieve this acceleration: a nonlinear dimensionality reduction step provided by an autoencoder, and a time integration step provided by a neural-network based model for the resultant latent space dynamics (the neural ODE). This work explores the applicability of such autoencoder-based neural ODE strategies for PDEs in which advection terms play a critical role. More specifically, alongside predictive demonstrations, physical insight into the sources of model acceleration (i.e., how the neural ODE achieves its acceleration) is the scope of the current study. Such investigations are performed by quantifying the effects of both autoencoder and neural ODE components on latent system time-scales using eigenvalue analysis of dynamical system Jacobians. To this end, the sensitivity of various critical training parameters – de-coupled versus end-to-end training, latent space dimensionality, and the role of training trajectory length, for example – to both model accuracy and the discovered latent system timescales is quantified. Furthermore, this work specifically uncovers the key role played by the training trajectory length (the number of rollout steps in the loss function during training) on the latent system timescales: larger trajectory lengths correlate with an increase in limiting neural ODE time-scales, and optimal neural ODEs are found to recover the largest time-scales of the full-order (ground-truth) system. Demonstrations are performed across fundamentally different unsteady fluid dynamics configurations influenced by advection: (1) the Kuramoto–Sivashinsky equations (2) Hydrogen-Air channel detonations (the compressible reacting Navier–Stokes equations with detailed chemistry), and (3) 2D Atmospheric flow.

Advection-dominated dynamical systems↗