Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reactive Flow”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Validation of Reactive Flow Capabilities in PAGOSA

Several reactive flow modeling capabilities are studied and applied to physics problems of interest in the PAGOSA hydrocode with an emphasis on validation. The Davis equation of state and the Arrhenius shock temperature dependent WSD (AWSD) reactive flow model have been implemented by the author in PAGOSA over the last year. These high-explosive modeling options along with the previously existing HE-JWL equation of state and the SURF/SURFplus reactive flow models are discussed, exercised and compared. Validation of the models analyze numerical results against experimental data from cylinder tests, two-stage gas gun experiments and proton radiography (pRad) experiments. Applications of the modeling options are explored involving both HMX based PBX 9501 and TATB based PBX 9502. The effect of mesh size is considered for prediction of phenomena such as corner turning and dead zones, and the relative code performance of the modeling options is discussed.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2023

Evolution of Two-Phase High Explosive Reactive Flow [Slides]

Reactive hydrodynamics involve rapid conversion of reactants to products along a detonation wave and needs a closure rule: P-T equilibrium. To calculate P-T equilibrium, we use a root finding method, which can be a computational bottleneck. We evolve the guess for the root finder rather than using the previous value to reduce iterations. Evolving products volume fraction Φ p for an initial root finding guess does reduce root finding iterations and seems to give the best results immediately following a detonation wave (reducing root finding iterations from 5 to 2).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Coupling Surface Flow with High-performance Subsurface Reactive Flow and Transport Code PFLOTRAN

Water exchange between the surface and subsurface is important for both water resource management and environmental protection. In this paper, we develop coupled surface and subsurface flow simulation capability in a parallel subsurface flow and reactive transport code PFLOTRAN. We sequentially couple the diffusion wave-based surface flow with the subsurface flow governedby the Richards equation in PFLOTRAN. These two flow domains are linked with a boundary condition switching method that ensures continuity of pressure and flux at the surface-subsurface interface. We verify the coupled code against other existing hydrologic models and observation data using a number of numerical experiments. The coupled hydrological model exhibits good performance in strong parallel scaling tests. The new coupled surface and subsurface simulator significantly advance community simulation capability towards improving integrated hydrologic and biogeochemical understanding of complex systems such as watersheds and river corridors. Keywords: Surface flow, Integrated hydrological modeling, Boundary condition switching, Parallel computing

Wu, Runjian↗
DOE OSTI2021

Reactive Flows in Porous Media: Challenges in Theoretical and Numerical Methods

We review theoretical and computational research, primarily from the past 10 years, addressing the flow of reactive fluids in porous media. The focus is on systems where chemical reactions at the solid–fluid interface cause dissolution of the surrounding porous matrix, creating nonlinear feedback mechanisms that can often lead to greatly enhanced permeability. Here, we discuss insights into the evolution of geological forms that can be inferred from these feedback mechanisms, as well as some geotechnical applications such as enhanced oil recovery, hydraulic fracturing, and carbon sequestration. Until recently, most practical applications of reactive transport have been based on Darcy-scale modeling, where averaged equations for the flow and reactant transport are solved. We summarize the successes and limitations of volume averaging, which leads to Darcy-scale equations, as an introduction to pore-scale modeling. Pore-scale modeling is computationally intensive but offers new insights as well as tests of averaging theories and pore-network models. We include recent research devoted to validation of pore-scale simulations, particularly the use of visual observations from microfluidic experiments.

wormhole formation↗
DOE OSTI2025

Temperature-based reactive flow model for triaminotrinitrobenzene (TATB) plastic bonded explosives

A new reactive flow model is presented for triaminotrinitrobenzene (TATB)-based plastic bonded explosives, applicable to shock initiation and steady detonation problems of differing initial temperature. Temperature disequilibrium is assumed between unreacted explosive, material in the vicinity of compressed defects (called hot spots), and reaction products. The model incorporates temperature-dependent decomposition reaction rates. Particularly, Arrhenius model parameters were derived from quantum-based molecular dynamics simulations of TATB decomposition. Further, a model of detonation carbon aggregation is incorporated, describing the slow release of energy inherent to detonation in TATB-based materials. Model parameters were calibrated against gas gun shock initiation experiments and steady detonation rate stick tests. The predictive ability of the model in the shock initiation regime is tested against recent thin pulse experiments. The model is found to perform equally well in predicting the size-effect curve of ambient, cold, and hot rate sticks. The present work demonstrates the viability of incorporating results from subscale simulations into a continuum-scale reactive flow model.

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

Reactive Flow Characteristic Equations

The characteristic equations are derived for the reactive flow PDEs in 1-dimension with 1 irreversible reaction. They differ from the non-reacting case only by source terms proportional to the reaction rate.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Transforming Energy Through Computational Excellence: Bringing Low Mach Number Reactive Flow Simulations at the Exascale

PeleLMeX's unique capabilities are allowing for reactive flow modeling at unprecedented scales and a reasonable time and cost. The code is currently being extended to tackle more practical, design-oriented simulations by implementing Large Eddy Simulation and data-driven chemical models, providing a fast but accurate tool for engineers considering the emergence of GPU-accelerated platforms. These extensions are critical for enabling the physical insight required to design the next generation of combustion devices as a key component of a renewable energy future.

MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Temperature-based reactive flow model for TATB plastic bonded explosives

A new reactive flow model is presented for TATB-based plastic bonded explosives, applicable to shock initiation and detonation problems of differing initial temperature. Temperature disequilibrium is assumed between unreacted explosive, material in the vicinity of compressed defects (called hot spots), and detonation products. The model incorporates temperature-dependent decomposition reaction rates. Arrhenius model parameters were derived from subscale quantum-based molecular dynamics simulations of shock-induced TATB decomposition. Further, a model of detonation carbon aggregation is incorporated, describing the late-time slow release of energy inherent to detonation in TATB-based materials. The predictive ability of the model in the shock initiation regime is tested against recent thin pulse experiments on PBX 9502. Here, the model is found to perform equally well in predicting the size effect curve of ambient, cold, and hot PBX 9502 rate sticks.

Chemistry - Chemical explosives↗
DOE OSTI2025

AWSD reactive flow model for PBX 9404

An Arrhenius–Wescott–Stewart–Davis (AWSD) reactive flow model for high explosive PBX 9404 is developed. We specifically calibrate an AWSD model for PBX 9404 by fitting equations of state for reactants and detonation products to the results of thermochemical calculations and to experimental data from multiple sources. The calibrated equations of state are then coupled with an Arrhenius rate law based on shock temperature that describes the reaction progress during PBX 9404 detonation. The parameter values in the rate law are calibrated to experimental gas-gun data and diameter effect results. The results of the calibrated AWSD model are in strong agreement with available experimental data for PBX 9404. A similar level of agreement between predicted and experimental results is observed when the calibrated model is validated on data that were not used in the model parameterization procedure. Our results illustrate that the AWSD model is capable of accurately describing the many important properties and observables in the reactive burn of PBX 9404. Because of the historical significance of PBX 9404 in high explosives research and its current use in aging studies, this work provides an important model of a legacy material, which can be used to make comparisons to new high explosive formulations.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2021

A reactive flow model for the 3,3'-diamino-4,4'-azoxyfurazan based plastic bonded explosive (PBX 9701)

Here, this paper describes the calibration of the Arrhenius Wescott–Stewart–Davis (AWSD) reactive flow model for the recently developed high explosive PBX 9701, which consists of 97% 3,3'-diamino-4,4'-azoxyfurazan (DAAF) and 3% FK-800 binder by weight. DAAF-based explosive formulations have several desirable qualities as they are relatively insensitive to non-shock insults but have higher performance than triaminotrinitrobenzene based formulations. Equations of state for the explosive reactants and products are calibrated using a combination of existing and new experimental data and theoretical calculations. The AWSD rate law calibration utilizes both one-dimensional shock-to-detonation and multidimensional rate stick experiments to capture the shock initiation and propagating detonation regimes. Validation of the calibrated model is demonstrated through comparison with recent gas-gun experiments.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2022

Pore-Scale Modeling of Mineral Growth and Nucleation in Reactive Flow

A fundamental understanding of mineral precipitation kinetics relies largely on microscopic observations of the dynamics of mineral surfaces exposed to supersaturated solutions. Deconvolution of tightly bound transport, surface reaction, and crystal nucleation phenomena still remains one of the main challenges. Particularly, the influence of these processes on texture and morphology of mineral precipitate remains unclear. This study presents a coupling of pore-scale reactive transport modeling with the Arbitrary Lagrangian-Eulerian approach for tracking evolution of explicit solid interface during mineral precipitation. It incorporates a heterogeneous nucleation mechanism according to Classical Nucleation Theory which can be turned “on” or “off.” This approach allows us to demonstrate the role of nucleation on precipitate texture with a focus at micrometer scale. In this work precipitate formation is modeled on a 10 micrometer radius particle in reactive flow. The evolution of explicit interface accounts for the surface curvature which is crucial at this scale in the regime of emerging instabilities. The results illustrate how the surface reaction and reactive fluid flow affect the shape of precipitate on a solid particle. It is shown that nucleation promotes the formation of irregularly shaped precipitate and diminishes the effect of the flow on the asymmetry of precipitation around the particle. The observed differences in precipitate structure are expected to be an important benchmark for reaction-driven precipitation in natural environments.

58 GEOSCIENCES↗
DOE OSTI2022

Predictive Large-Eddy Simulation of Supercritical-Pressure Reactive Flows in the Cold Ignition Regime

This report describes a one-year study performed under DOE sponsorship, continuing the investigation of high-pressure turbulent reactive flows. The interest is in the effect of the chemical species distribution in high-pressure turbulent flows in the presence of strong temperature gradients as would occur during reactions in realistic flows. The prime example of such flows are boundary layers in which the wall is at a lower temperature than that of the fluid, as would be the case in Diesel engines. Previous DOE BES work in the program further highlighted the importance of the boundary layer as a configuration for fundamental studies: soot formation in boundary layers is still a problem poorly understood and depends on the availability of particular chemical species at that location, and when boundary layer Large Eddy Simulations results were compared with measurements, the agreement was unfavorable, showing that this important ‘unit’ problem is not well understood. To understand this unit problem, modeling and Direct Numerical Simulations of this unit problem were conducted for the simplest possible multispecies system, that is a binary-species system. The results discovered a new phenomenon, that is, Soret effect induced uphill diffusion. The far-reaching implication is that through an imposed wall/fluid temperature difference it is possible to control the distribution of the species in the boundary layer. Because the results have been documented in one paper published in the refereed literature, and also in conference papers, this final report is in the form of an Executive Summary succinctly describing the results and putting them in perspective with respect to existing information. The refereed and conference papers published are individually listed as Appendices and attached to this report. One manuscript is still in progress and is thus not listed.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Integration and validation of some modules for modelling of high-speed chemically reactive flows in two-phase gas-droplet mixtures

Three modules are integrated into the built-in OpenFOAM rhoCentralFoam solver towards accurate and efficient modelling of high-speed chemically reactive flows in two-phase gas-droplet mixtures within the OpenFOAM 10.0 framework. The first module is the mixture-averaged diffusion model. The second module is the built-in OpenFOAM Lagrangian solver coupled with optimised droplet drag coefficient and convective heat transfer coefficient sub-models. The last module is a sparse stiff chemistry solver based on dynamic adaptive hybrid integration (AHI-S). The optimised droplet sub-models are first verified in correct implementation for subsequent simulations in this work. Further, they show good accuracy against experimental and analytical data in the modelling of ammonia droplet acceleration and cooling in the flowing and/or low-temperature air. The accuracy and efficiency gains related to the mixture-averaged diffusion model and the AHI-S chemistry solver are examined by simulating 1-D detonation propagation in ammonia droplet-free/laden ammoniaoxygen mixtures. Numerical results of detonation propagation speed, gaseous temperature, density, and species distributions around the induction zone show good agreement with experimental data and analytical solutions. Compared to the built-in OpenFOAM diffusion model, the mixture-averaged diffusion model provides different numerical predictions of pulsating instabilities in detonation propagation. It shows better accuracy in depicting the detonation structure within the droplet-free section attributed to improved multi-component diffusion modelling. Compared to the built-in OpenFOAM solver EulerImplicit (backward Euler), the AHI-S chemistry solver reduces the computational cost by around 50%. It achieves satisfactory accuracy in calculating detonation propagation speed within the droplet-free section with the optimal efficiency when the safety factor, β, equals 0.5.

42 ENGINEERING↗
DOE OSTI2022

Simulations of Low Mach Number Reactive Flows Coupled with Electric Fields

Applying electric fields to flames can provide an effective, non-invasive control of the combustion process. In this work, an algorithm recently introduced to handle the interactions between reactive flows and an electric field is extended to handle multiple dimensions and adaptive mesh refinement (AMR), and applied to a series of laminar flames. Experimental data allows to validate the ability of the algorithm to capture the mechanisms driving the flame response to electric fields.

electric fields↗
DOE OSTI2022

Simulations of Low Mach Number Reactive Flows Coupled with Electric Fields: Preprint

Applying electric fields to flames can provide an effective, non-invasive control of the combustion process. In this work, an algorithm recently introduced to handle the interactions between reactive flows and an electric field is extended to handle multiple dimensions and adaptive mesh refinement (AMR), and applied to a series of laminar flames. Experimental data allows to validate the ability of the algorithm to capture the mechanisms driving the flame response to electric fields.

electric field↗
DOE OSTI2023

PeleLMeX: an AMR Low Mach Number Reactive Flow Simulation Code without level sub-cycling

PeleLMeX simulates chemically reacting low Mach number flows with block-structured adaptive mesh refinement (AMR). The code is built upon the AMReX library, which provides the underlying data structures and tools to manage and operate on them across massively parallel computing architectures. PeleLMeX algorithmic features are inherited from its predecessor PeleLM but key improvements allow representation of more complex physical processes. Together with its compressible flow counterpart PeleC, the thermo-chemistry library PelePhysics and the multi-physics library PeleMP, it forms the Pele suite of open-source reactive flow simulation codes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A Fully Explicit Integrator for Modeling Astrophysical Reactive Flows

Simulating complex astrophysical reacting flows is computationally expensive—reactions are stiff and typically require implicit integration methods. The reaction update is often the most expensive part of a simulation, which motivates the exploration of more economical methods. In this research note, we investigate how the explicit Runge–Kutta–Chebyshev (RKC) method performs compared to an implicit method when applied to astrophysical reactive flows. These integrators are applied to simulations of X-ray bursts arising from unstable thermonuclear burning of accreted fuel on the surface of neutron stars. We show that the RKC method performs with similar accuracy to our traditional implicit integrator, but is more computationally efficient when run on CPUs.

GPU computing↗
DOE OSTI2020

Transported PDF Modeling of Compressible Turbulent Reactive Flows by using the Eulerian Monte Carlo Fields Method

Although the transported probability density function (PDF) method has been developed for decades, its application has been mainly focused on the low-Mach number flow problems. This work extends the transported PDF method to compressible flow problems. The Eulerian Monte Carlo fields (EMCF) solution method is employed to solve the transported PDF equation for compressible flow problems. A pseudo stagnation enthalpy is introduced and its stochastic partial differential equation is derived to ensure total energy conservation numerically. A new mixing model called interaction by partial exchange with mean (IPEM) is introduced to expand the available choices of mixing models for the EMCF method. The consistency of the EMCF method is examined for solving the transported PDF equation. Numerical implementation details are discussed, such as the density coupling between the compressible flow solver and the EMCF solver, discretization schemes for the mixing terms and the stochastic terms. The implemented compressible flow solver coupled with the EMCF solver is verified and validated in a series of test cases with increasing level of complexity, ranging from a statistically one-dimensional turbulent mixing layer to a self-excited resonance model rocket combustor. It is observed that in general with the increase of compressibility, there is an increase in the sensitivity of the modeling results to the different models and algorithms. This makes it necessary to develop a thorough understanding of the model sensitivity in order to develop a robust and accurate simulation solver for highly compressible turbulent reactive flows. The thermo-acoustic instability inside the model rocket combustor case is captured reasonably, which demonstrates the overall capability of the developed compressible turbulent combustion solver based on the transported PDF method.

42 ENGINEERING↗