Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gas and solid 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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

An experimental investigation of gas jets in confined swirling air flow

The fluid dynamics of jets in confined swirling flows which is of importance to designers of turbine combustors and solid fuel ramjets used to power missiles fired from cannons were examined. The fluid dynamics of gas jets of different densities in confined swirling flows were investigated. Mean velocity and turbulence measurements are made with a one color, one component laser velocimeter operating in the forward scatter mode. It is shown that jets in confined flow with large area ratio are highly dissipative which results in both air and helium/air jet centerline velocity decays. For air jets, the jet like behavior in the tube center disappears at about 20 diameters downstream of the jet exit. This phenomenon is independent of the initial jet velocity. The turbulence field at this point also decays to that of the background swirling flow. A jet like behavior in the tube center is noticed even at 40 diameters for the helium/air jets. The subsequent flow and turbulence field depend highly on the initial jet velocity. The jets are fully turbulent, and the cause of this difference in behavior is attributed to the combined action swirl and density difference. This observation can have significant impact on the design of turbine combustors and solid fuel ramjets subject to spin.

Mongia, H.↗
DOE OSTI2022

Memory access optimization for particle operations in computational fluid dynamics-discrete element method simulations

Computational Fluid Dynamics - Discrete Element Method is used to model gas-solid systems in several applications in energy, pharmaceutical and petrochemical industries. Computational performance bottlenecks often limit the problem sizes that can be simulated at industrial scale. The data structures used to store several millions of particles in such large-scale simulations have a large memory footprint that does not fit into the processor cache hierarchies on current high-performance-computing platforms, leading to reduced computational performance. This paper specifically addresses this aspect of memory access bottlenecks in industrial scale simulations. The use of space-filling curves to improve memory access patterns is described and their impact on computational performance is quantified in both shared and distributed memory parallelization paradigms. The Morton space filling curve applied to uniform grids and k-dimensional tree partitions are used to reorder the particle data-structure thus improving spatial and temporal locality in memory. The performance impact of these techniques when applied to two benchmark problems, namely the homogeneous-cooling-system and a fluidized-bed, are presented. We report these optimization techniques lead to approximately two-fold performance improvement in particle focused operations such as neighbor-list creation and data-exchange, with ~ 1.5 times overall improvement in a fluidization simulation with 1.27 million particles.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2002

Assessment of Various Flow Solvers Used to Predict the Thermal Environment inside Space Shuttle Solid Rocket Motor Joints

It is very important to accurately predict the gas pressure, gas and solid temperature, as well as the amount of O-ring erosion inside the space shuttle Reusable Solid Rocket Motor (RSRM) joints in the event of a leak path. The scenarios considered are typically hot combustion gas rapid pressurization events of small volumes through narrow and restricted flow paths. The ideal method for this prediction is a transient three-dimensional computational fluid dynamics (CFD) simulation with a computational domain including both combustion gas and surrounding solid regions. However, this has not yet been demonstrated to be economical for this application due to the enormous amount of CPU time and memory resulting from the relatively long fill time as well as the large pressure and temperature rising rate. Consequently, all CFD applications in RSRM joints so far are steady-state simulations with solid regions being excluded from the computational domain by assuming either a constant wall temperature or no heat transfer between the hot combustion gas and cool solid walls.

Wang, Qun-Zhen↗
DOE OSTI2025

Atomic dynamics of gas-dependent oxide reducibility

Understanding oxide reduction is critical for advancing metal production, catalysis and energy technologies. Although carbon monoxide (CO) and hydrogen (H 2 ) are widely used reductants, the mechanisms by which they work are often presumed to be similar, both involving lattice oxygen removal. However, because of growing interest in replacing CO with H 2 to lower CO 2 emissions, distinguishing gas-specific reduction pathways is critical. Yet, capturing these atomic-scale processes under reactive gas and high-temperature conditions remains challenging. Here we use environmental transmission electron microscopy, which is capable of real-time, atomic-resolution imaging of gas–solid redox reactions to directly visualize the gas-dependent oxide reduction dynamics in NiO. We show that CO drives surface nucleation and the growth of metallic Ni islands, leading to self-limiting surface metallization. Conversely, H 2 activates a coupled surface-to-bulk transformation, where protons from dissociated H 2 infiltrate the oxide lattice to promote the inward migration of surface-generated oxygen vacancies and enabling bulk metallization. By contrast, oxygen vacancies formed by CO remain confined near the surface, where they rapidly form a metallic Ni layer that inhibits further reduction. Furthermore, these results reveal distinct atomistic pathways for CO and H 2 and provide insights that may guide metallurgical processes and catalyst design.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Correction to: Investigating the rheology of fluidized and non-fluidized gas-particle beds: implications for the dynamics of geophysical flows and substrate entrainment

We report natural geophysical mass flows are among the most complex granular systems and their dynamics are often modified by the presence of an interstitial fluid. Prediction of their runout requires the development of models estimating the solid stresses in these hazardous currents wherein excess pore-fluid pressure can develop. We use discrete element modelling (DEM-CFD) with a Coarse-Graining post-processing technique (CG) to investigate the rheology of unsteady gas-particle fluidized to non-fluidized granular beds placed on horizontal and inclined planes. Similar to fluidized beds immersed in viscous fluids, the effective friction coefficient of air-fluidized beds can be defined as a function of the classic μ(I)-rheology and the non-dimensional fluid or solid pressure to explain the failure and dynamics of granular flows with excess pore pressure on inclines. However, dilation imposed by fluid drag and particle collisions in gas-particle fluidized beds can drastically change its effective frictional properties. In contrast with the common assumption in water-particle flows that granular temperature is negligible, in our gas-particle simulations, the contribution of the velocity fluctuations to the stress tensor is significant. Hence, the shear stress is found to be non-zero even when the flow is fully fluidized in the inertial regime. These results suggest the need to better understand velocity fluctuations to predict the effective viscosity of sheared fluidized mixtures and are relevant for many applications. Notably, a unified approach is useful for many geophysical flows that encompass a range of fluidization conditions in a single flow such as pyroclastic density currents and snow avalanches.

58 GEOSCIENCES↗
DOE OSTI2020

Multiresolution GPC-Structured Control of a Single-Loop Cold-Flow Chemical Looping Testbed

Chemical looping is a near-zero emission process for generating power from coal. It is based on a multi-phase gas-solid flow and has extremely challenging nonlinear, multi-scale dynamics with jumps, producing large dynamic model uncertainty, which renders traditional robust control techniques, such as linear parameter varying H ∞ design, largely inapplicable. This process complexity is addressed in the present work through the temporal and the spatiotemporal multiresolution modeling along with the corresponding model-based control laws. Namely, the nonlinear autoregressive with exogenous input model structure, nonlinear in the wavelet basis, but linear in parameters, is used to identify the dominant temporal chemical looping process dynamics. The control inputs and the wavelet model parameters are calculated by optimizing a quadratic cost function using a gradient descent method. The respective identification and tracking error convergence of the proposed self-tuning identification and control schemes, the latter using the unconstrained generalized predictive control structure, is separately ascertained through the Lyapunov stability theorem. The rate constraint on the control signal in the temporal control law is then imposed and the control topology is augmented by an additional control loop with self-tuning deadbeat controller which uses the spatiotemporal wavelet riser dynamics representation. The novelty of this work is three-fold: (1) developing the self-tuning controller design methodology that consists in embedding the real-time tunable temporal highly nonlinear, but linearly parametrizable, multiresolution system representations into the classical rate-constrained generalized predictive quadratic optimal control structure, (2) augmenting the temporal multiresolution loop by a more complex spatiotemporal multiresolution self-tuning deadbeat control loop, and (3) demonstrating the effectiveness of the proposed methodology in producing fast recursive real-time algorithms for controlling highly uncertain nonlinear multiscale processes. The latter is shown through the data from the implemented temporal and augmented spatiotemporal solutions of a difficult chemical looping cold flow tracking control problem.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1982

Wind tunnel tests of space shuttle solid rocket booster insulation material in the aerothermal tunnel c

Wind tunnel tests of the space shuttle Solid Rocket Booster Insulation were conducted in the von Karman Gas Dynamics Facility Tunnel C. For these tests, Tunnel C was run at Mach 4 with a total temperature of 1100-1440 and a total pressure of 100 psia. Cold wall heating rates were changed by varying the test article support wedge angle. Selected results are presented to illustrate the test techniques and typical data obtained.

Hartman, A. S.↗
DOE OSTI2024

Solid, structured composite neutron detectors with high dynamic range capability

Neutron detectors are essential across disciplines such as fundamental science, nuclear security, safeguards, and civilian applications. While 3 He-filled gas proportional counters have long been revered for their efficacy in detecting thermal neutrons and praised for their efficiency, neutron/gamma discrimination, and stability, the scarcity of 3He has spurred a search for alternatives. Here, we explore a solid structured scintillating particle composite (SPC) consisting of 6 Li-containing scintillating glass particles within an acrylic matrix as a neutron detector for high dynamic range applications. We show for the first time that an SPC neutron detector can boast an intrinsic detection efficiency of 0.261% for pure 252 Cf fission neutrons and an overall neutron detection efficiency of (0.546 ± 0.003)% at the Neutron Free-in-Air facility while being able to function in an intense gamma-ray environment. We also show that the SPC neutron detector supports fast neutron capture times and enables a dual-readout scheme that extends the detector dynamic range to high incident neutron fluxes. A scalable fabrication process allows for tailoring the SPC detector properties to the requirements of specific applications. Good agreement is found between the experimental results taken with a National Institute of Standards and Technology traceable 252 Cf source and the coupled MCNP6 and optical-ray-tracing simulations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Halide Perovskites Breathe Too: The Iodide–Iodine Equilibrium and Self-Doping in Cs 2 SnI 6

The response of an oxide crystal to the atmosphere can be personified as breathing-a dynamic equilibrium between O 2 gas and O 2– anions in the solid. We characterize the analogous defect reaction in an iodide double-perovskite semiconductor, Cs 2 SnI 6 . Here, I 2 gas is released from the crystal at room temperature, forming iodine vacancies. The iodine vacancy defect is a shallow electron donor and is therefore ionized at room temperature; thus, the loss of I2 is accompanied by spontaneous n-type self-doping. Conversely, at high I 2 pressures, I 2 gas is resorbed by the perovskite, consuming excess electrons as I 2 is converted to 2I – . Halide mobility and irreversible halide loss or exchange reactions have been studied extensively in halide perovskites. However, the reversible exchange equilibrium between iodide and iodine [2I – (s) ↔ I 2(g) + 2e – ] described here has often been overlooked in prior studies, though it is likely general to halide perovskites and operative near room temperature, even in the dark. An analysis of the 2I – (s) /I 2(g) equilibrium thermodynamics and related transport kinetics in single crystals of Cs 2 SnI 6 therefore provides insight toward achieving stable composition and electronic properties in the large family of iodide perovskite semiconductors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1988

Planetary environments and the conditions of life

Geophysical models of the first 600 Ma ofthe earth's history following accretion and core formation point to a period of great environmental disequilibrium. In such an environment, the passage of energy from the earth's interior and from the sun through gas-liquid-solid domains and their boundaries with each other generated a dynamically interacting, complex hierarchy of self-organized structures ranging from bubbles at the sea-air interface to tectonic plates. The ability of a planet to produce such a hierarchy is speculated to be a prerequisite to the origin and sustenance of life. The application of this criterion to Mars argues against the origin of Martian life.

Chang, S.↗
NASA NTRS1999

Dynamic Nucleation of Supercooled Melts and Measurement of the Surface Tension and Viscosity

We investigate the phenomenon of acoustic pressure-induced nucleation by using a novel approach involving the large amplitude resonant radial oscillations and collapse of a single bubble intentionally injected into a supercooled liquid. Using a combination of previously developed and proven techniques, the bubble is suspended in a fluid host by an ultrasonic field which supplies both the levitation capability as well as the forcing of the radial oscillations. We observe the effects of an increase in pressure (due to bubble collapse) in a region no larger than 100 microns within the supercooled melt to rigorously probe the hypothesis of pressure-induced nucleation of the solid phase. The use of single bubbles operating in narrow temporal and spatial scales will allow the direct and unambiguous correlation between the origin and location of the generation of the disturbance and the location and timing of the nucleation event. In a companion research effort, we are developing novel techniques for the non-contact measurements of the surface tension and viscosity of highly viscous supercooled liquids. Currently used non-invasive methods of surface tension measurement for the case of undercooled liquids generally rely of the quantitative determination of the resonance frequencies of drop shape oscillations, of the dynamics of surface capillary waves, or of the velocity of streaming flows. These methods become quickly ineffective when the liquid viscosity rises to a significant value. An alternate and accurate method which would be applicable to liquids of significant viscosity is therefore needed. We plan to develop such a capability by measuring the equilibrium shape of levitated undercooled melt droplets as they undergo solid-body rotation. The experimental measurement of the characteristic point of transition (bifurcation point) between axisymmetric and two-lobed shapes will be used to calculate the surface tension of the liquid. Such an approach has already been validated through the experimental verification of numerical modeling results. The experimental approach involves levitation, melting, and solidification of undercooled droplets using a hybrid ultrasonic-electrostatic technique in both a gaseous as well as a vacuum environment. A shape relaxation method will be investigated in order to derive a reliable method to measure the viscosity of undercooled melts. The analysis of the monotonic relaxation to equilibrium shape of a drastically deformed and super-critically damped free drop has been used to derive interfacial tension of immiscible liquid combinations where one of the component has high viscosity. A standard approach uses the initial elongation of a droplet through shear flows, but an equivalent method could involve the initial deformation of a drop levitated in a gas by ultrasonic radiation pressure, electric stresses, or even solid body rotation. The dynamic behavior of the free drop relaxing back to equilibrium shape will be modeled, and its characteristic time dependence should provide a quantitative means to evaluate the liquid viscosity.

Trinh, E. H.↗
DOE OSTI2025

Coarse Graining Discrete Element Method Information in Particle-in-Cell Length Scales Using a Machine Learning Approach

This report details the development of a machine learning (ML)-driven framework to coarse-grain inter-particle collision dynamics from high-fidelity Discrete Element Method (DEM) simulations to Particle-in-Cell (PIC) scales for gas-solid systems. Traditional PIC models, while computationally efficient, rely on empirical granular stress formulations that fail to capture the full complexity of collision physics, particularly the heterogeneity in particle dynamics. This study adopts a bottom-up approach, integrating insights from DEM simulations to improve the physical fidelity and interpretability of PIC-scale models.

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

First-principles thermodynamic assessment of Sr-containing secondary phase formation in strontium-substituted lanthanum manganites for solid oxide cell applications

Sr-secondary phase formation is a potentially significant degradation mode with direct impact upon solid-oxide cell (SOC) commercial viability. A first-principles based thermodynamic study was performed for La 1−x Sr x MnO 3±δ (LSM) perovskites to assess their stability against formation of different Sr-secondary phases, including SrO, SrCrO 4 , SrSO 4 , SrCO 3 , and Sr(OH) 2 , for SOC applications. The Sr-secondary phase formation reaction free energies were determined via a thermodynamic model by combining ab initio lattice dynamics calculations for the solid phases and ab initio thermodynamic data for the gas phases. The current approach expands the previously reported thermodynamic modeling studies by integrating first-principles based point defect equilibria into the thermodynamic analysis. The modeling results obtained using this new approach indicate an increased tendency to form the SrO oxide upon decreasing the oxygen partial pressure. Additionally, the enhancing factors to form the Sr-related secondary phase from the associated SrO activity in LSM are further quantified by considering the equilibrium of SrO reacting with the contaminant gas species as a function of temperature and pressure.

defect thermodynamics modeling↗
DOE OSTI2020

Three Jahn-Teller States of Matter in Spin-Crossover System Mn(taa)

Three high-spin phases recently discovered in the spin-crossover system Mn(taa) are identified through analysis by a combination of first-principles calculations and Monte Carlo simulation as a low-temperature Jahn-Teller ordered (solid) phase, an intermediate-temperature dynamically correlated (liquid) phase, and an uncorrelated (gas) phase. In particular, the Jahn-Teller liquid phase arises from competition between mixing with low-spin impurities, which drive the disorder, and intermolecular strain interactions. Additionally, the latter are a key factor in both the spin-crossover phase transition and the magnetoelectric coupling. Jahn-Teller liquids may exist in other spin-crossover materials and materials that have multiple equivalent Jahn-Teller axes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Operation and Dynamic Control of Hybrid Power Systems

Poster describing the planned startup process of a solid oxide gas turbine hybrid power system to be constructed at a 20 kW scale. The work is part of the ARPA-E project award to NETL under the INTEGRATE program.

Harun, Nor Farida↗
NASA NTRS1973

The starting transient of solid propellant rocket motors with high internal gas velocities

A comprehensive analytical model which considers time and space development of the flow field in solid propellant rocket motors with high volumetric loading density is described. The gas dynamics in the motor chamber is governed by a set of hyperbolic partial differential equations, that are coupled with the ignition and flame spreading events, and with the axial variation of mass addition. The flame spreading rate is calculated by successive heating-to-ignition along the propellant surface. Experimental diagnostic studies have been performed with a rectangular window motor (50 cm grain length, 5 cm burning perimeter and 1 cm hydraulic port diameter), using a controllable head-end gaseous igniter. Tests were conducted with AP composite propellant at port-to-throat area ratios of 2.0, 1.5, 1.2, and 1.06, and head-end pressures from 35 to 70 atm. Calculated pressure transients and flame spreading rates are in very good agreement with those measured in the experimental system.

Peretz, A.↗
DOE OSTI2020

A Method for Obtaining Liquid–Solid Adsorption Rates from Molecular Dynamics Simulations: Applied to Methanol on Pt(111) in H 2 O

Adsorption is an important step in heterogeneous catalysis as it predetermines how many reactant molecules can participate in a surface reaction per unit time. While the rate of adsorption processes is well studied in gas–solid adsorption in both theory and experiment, such rates are still not well studied for liquid–solid adsorption. This is partly because the ever-changing configurations of liquid-phase solvent molecules impede the ability to study a molecule approaching a surface from a liquid phase by either experiment or theory. In this work, we develop a method using molecular dynamics (MD) simulations to study the rate of adsorption in liquid–solid adsorption processes. Specifically, we use MD to model the diffusion of a methanol molecule in aqueous solvent and its adsorption to a Pt(111) surface. We find that by approximating the solute motion as following the same displacement rates as a random walk model, the adsorbed and non-adsorbed states of the methanol molecule near the Pt(111) surface can be discerned and quantified. In particular, this methodology enables extracting a sticking coefficient and a macroscopically relatable adsorption rate. This method can be applied to arbitrary types of reactants and surfaces, as well as different liquid environments, thus providing a general tool for predicting quantitative adsorption rates of liquid–solid adsorption systems.

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

First-Principles Thermodynamic Assessments of Sr-Containing Secondary Phase Formation in La1-xSrxMnO3±δ Perovskites for Solid Oxide Cell Applications

Sr-secondary phase formation is a potentially significant degradation mode threatening solid-oxide cell (SOC) commercial viability. A first-principles thermodynamic study was performed for rhombohedral perovskite (La1-xSrx) MnO3±δ (LSM) to assess its stability against Sr secondary phase formation in SOC applications. In this work, the Sr secondary phase formation reaction free energies were determined by combining ab initio lattice dynamics calculations for the solid phases and an ab initio thermodynamics approach for the gas phases. Furthermore, this approach goes beyond previous thermodynamic modeling studies by integrating first-principles based point-defect equilibria into the analyses. The modeling results indicate an increased tendency to form SrO oxide from LSM upon decreasing the oxygen partial pressure. Additionally, enhancing factors to form the Sr-related secondary phase from the associated SrO activity in LSM are further quantified by considering the equilibrium of SrO reacting with contaminant gas species as a function of temperature and gas pressure.

Defect and phase stability↗