Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

Incorporating geological structure into sensitivity analysis of subsurface contaminant transport

Simulating subsurface contaminant transport at the kilometer-scale often entails modeling reactive flow and transport within and through complex geologic structures. These structures are typically meshed by hand and as a result geologic structure is usually represented by one or a few deterministically generated geological models for uncertainty studies of flow and transport in the subsurface. Uncertainty in geologic structure can have a significant impact on contaminant transport. In this study, the impact of geologic structure on contaminant tracer transport in a shale formation is investigated for a simplified generic deep geologic repository for permanent disposal of spent nuclear fuel. An open-source modeling framework is used to perform a sensitivity analysis study on transport of two tracers from a generic spent nuclear fuel repository with uncertain location of the interfaces between the stratum of the geologic structure. The automated workflow uses sampled realizations of the geological structural model in addition to uncertain flow parameters in a nested sensitivity analysis. Concentration of the tracers at observation points within, in line with, and downstream of the repository are used as the quantities of interest for determining model sensitivity to input parameters and geological realization. Finally, the results of the study indicate that the location of strata interfaces in the geological structure has a first-order impact on tracer transport in the example shale formation, and that this impact may be greater than that of the uncertain flow parameters.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Balanced Fast-SpaciMS capillary configurations provide practically noninvasive channel-average measurements in catalytic monoliths

Spatially resolved capillary-inlet mass spectrometry (SpaciMS) provides a detailed picture of the spatiotemporal evolution of reaction network in catalytic monoliths. In the present work, we combine the SpaciMS experiments with a newly developed non-isothermal 3D CFD model for heterogeneously catalyzed reactive flows, including diffusion and permeation through the coated catalyst and channel wall. We explore how the capillary size and sampling rate can be balanced to minimize the impact on probed-channel conversions and provide species concentration measurements representative of the free channel average. In all studied configurations, the balanced sampling-rate fraction is noticeably higher than the corresponding capillary occlusion fraction. For a typical 350-micron capillary that occupies 11% of the channel cross- section, the balanced sampling rate represents 43% of the channel flow. Furthermore, the present work shows that together, a balanced SpaciMS configuration and numerical simulation provide accurate channel-averaged information of a monolithic catalyst under realistic operation conditions.

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

Operator inference for non-intrusive model reduction of systems with non-polynomial nonlinear terms

Here in this work we present a non-intrusive model reduction method to learn low-dimensional models of dynamical systems with non-polynomial nonlinear terms that are spatially local and that are given in analytic form. In contrast to state-of-the-art model reduction methods that are intrusive and thus require full knowledge of the governing equations and the operators of a full model of the discretized dynamical system, the proposed approach requires only the non-polynomial terms in analytic form and learns the rest of the dynamics from snapshots computed with a potentially black-box full-model solver. The proposed method learns operators for the linear and polynomially nonlinear dynamics via a least-squares problem, where the given non-polynomial terms are incorporated on the right-hand side. The least-squares problem is linear and thus can be solved efficiently in practice. The proposed method is demonstrated on three problems governed by partial differential equations, namely the diffusion–reaction Chafee–Infante model, a tubular reactor model for reactive flows, and a batch-chromatography model that describes a chemical separation process. The numerical results provide evidence that the proposed approach learns reduced models that achieve comparable accuracy as models constructed with state-of-the-art intrusive model reduction methods that require full knowledge of the governing equations.

42 ENGINEERING↗
DOE OSTI2022

Modeling of the effects of non-equilibrium excitation and electrode geometry on H 2 /air ignition in a nanosecond plasma discharge

In this work, we present the results of two-dimensional modeling of the effects of non-equilibrium excitation and electrode geometry on H 2 /air ignition in a nanosecond plasma discharge. A multiscale adaptive reduced chemistry solver for plasma assisted combustion (MARCS-PAC) based on PASSKEy discharge modeling package and compressible multi-component reactive flow solver ASURF+ is developed and validated. This model is applied to simulate the impact of non-equilibrium plasma excitation and electrode geometry and heat loss on the dynamics of the discharge from streamer to spark and ignition kernel development in a H 2 /air mixture with a pair of cylindrical electrodes. The results show that the plasmagenerated species (N 2 (A), N 2 (B), N 2 (a'), N 2 (C), O( 1 D), O and H) in the spark and afterglow significantly accelerate the ignition kernel development. The increase of discharge voltage at the same total discharge energy promotes the non-equilibrium active species production. It is found that the production of electronically excited species at higher reduced electric field strength is more efficient in enhancing ignition in comparison to the vibrational excitation and heating. Moreover, the 2D simulation clearly reveals that the electric field and active species distribution are highly non-uniform. The streamers are initiated at the sharp outer edges of the negative and positive electrodes by a strong electric field while the electric field is much weaker at the centerline of the electrodes. Furthermore, the simulations reveal that the ignition enhancement is sensitive to the variation of electrode shape, diameter, and gap size due to the changes of electric field distribution and location of streamer formation. A cylindrical electrode produces a larger discharge volume and ignition kernel than the parabolic and spherical electrodes, when the discharge is localized near the axis of the gap. It is found that there is a non-monotonic dependence of ignition kernel size on the electrode diameter and inter-electrode distance. The increase of electrode diameter and gap size above the optimal conditions leads to the reduction of ignition kernel volume, due to the decrease of active species concentration and gas temperature. At a larger electrode surface area and electrode diameter as well as smaller electrode gap size, the heat loss to electrode plays a greater role in reducing the ignition kernel size and slowing ignition kernel development. This work provides insights and guidance to understand the kinetic enhancement of non-equilibrium plasma and the effects of electrode geometries on ignition for the optimization ignitors in advanced engines.

42 ENGINEERING↗
DOE OSTI2026

Enabling probabilistic learning on manifolds through double diffusion maps

Here, we present a generative learning framework for probabilistic sampling that extends Probabilistic Learning on Manifolds (PLoM), which is designed to generate statistically consistent realizations of a random vector in a finite-dimensional Euclidean space, informed by a (representative) set of observations. In its original form, PLoM constructs a reduced-order probabilistic model by combining three main components: (a) kernel density estimation to approximate the underlying probability measure, (b) Diffusion Maps to characterize the manifold of the data, and (c) a reduced-order Itô Stochastic Differential Equation (ISDE) to sample from the learned distribution. However, its sampling dynamics are posed in the ambient space and the retained number of reduced coordinates is chosen by projection-reconstruction error. In practice, this often (i) requires more coordinates than the data’s intrinsic dimension to achieve stable sampling and (ii) lacks a smooth, basis-independent lifting back to the data domain; moreover, standard Diffusion Maps emphasize harmonic eigenfunctions and can miss non-harmonic latent structure. We address these limitations by decoupling geometry learning from sampling: a first Diffusion Maps pass identifies non-harmonic coordinates on which we formulate a full-order ISDE directly in the latent space, while Double Diffusion Maps captures multiscale geometric features and Geometric Harmonics (GH) learns a smooth lifting map to the ambient variables that is independent of the particular diffusion basis. This hybrid design preserves the system’s dynamical richness with a compact geometric representation and enables principled out-of-sample inference. The effectiveness and robustness of the proposed method are illustrated through two numerical studies: one based on data generated from two-dimensional Hermite polynomial functions and another based on high-fidelity simulations of a detonation wave in a reactive flow.

Double diffusion maps↗
DOE OSTI2022

Flash-X: A multiphysics simulation software instrument

Flash-X is a highly composable multiphysics software system that can be used to simulate physical phenomena in several scientific domains. It derives some of its solvers from FLASH, which was first released in 2000. Flash-X has a new framework that relies on abstractions and asynchronous communications for performance portability across a range of increasingly heterogeneous hardware platforms. Flash-X is meant primarily for solving Eulerian formulations of applications with compressible and/or incompressible reactive flows. It also has a built-in, versatile Lagrangian framework that can be used in many different ways, including implementing tracers, particle-in-cell simulations, and immersed boundary methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Molecular fluctuations inhibit intermittency in compressible turbulence

In the standard picture of fully developed turbulence, highly intermittent hydrodynamic fields are nonlinearly coupled across scales, where local energy cascades from large scales into dissipative vortices and large density gradients. Microscopically, however, constituent fluid molecules are in constant thermal (Brownian) motion, but the role of molecular fluctuations in large-scale turbulence is largely unknown, and with rare exceptions, it has historically been considered irrelevant at scales larger than the molecular mean free path. Recent theoretical and computational investigations have shown that molecular fluctuations can impact energy cascade at Kolmogorov length scales. Here, we show that molecular fluctuations not only modify energy spectrum at wavelengths larger than the Kolmogorov length in compressible turbulence, but also significantly inhibit spatio-temporal intermittency across the entire dissipation range. Using large-scale direct numerical simulations of computational fluctuating hydrodynamics, we demonstrate that the extreme intermittency characteristic of turbulence models is replaced by nearly Gaussian statistics in the dissipation range. These results demonstrate that the compressible Navier–Stokes equations should be augmented with molecular fluctuations to accurately predict turbulence statistics across the dissipation range. Our findings have significant consequences for turbulence modelling in applications such as astrophysics, reactive flows and hypersonic aerodynamics, where dissipation-range turbulence is approximated by closure models.

compressible turbulence↗
DOE OSTI2025

Cross-Scale Catalyst Modeling Applied to H 2 Storage and Release via Formic Acid

Here, we propose the Systems-to-Atoms (S2A) modeling framework that integrates the kinetics of reaction chemistry and structural configurations across various length scales with the aim of establishing a versatile template for multiscale modeling of reactive flow problems and to predict the operando activity of catalyst materials. The approach encompasses a microkinetic model to analyze surface reactions on individual facets of catalyst nanoparticles coupled with the computation of average surface reaction rates for catalyst nanoparticles of specific size distributions. Macro-homogeneous surface reaction kinetics are derived as a function of catalyst loading and used as input parameters for the continuum-scale reactor model. The cross-scale framework enables the optimization of catalyst utilization through reactor design and operating strategy. To demonstrate the framework, we studied the storage and release of hydrogen from formic acid, a promising liquid organic hydrogen carrier (LOHC), over Pd, Pt, and Cu catalysts. The framework predicts observed trends in formic acid dehydrogenation activity for catalysts with comparable weight loadings and metal particle diameters, demonstrating satisfactory quantitative alignment. Finally, the seamless transmission of parameter uncertainties between scales is also discussed.

08 HYDROGEN↗
DOE OSTI2021

Fully consistent Eulerian Monte Carlo fields method for solving probability density function transport equations in turbulence modeling

The probability density function (PDF) transport equation method is a sophisticated model for the closure of turbulent mixing and turbulent reactive flows. An efficient solution approach for solving the PDF transport equation has been vital for the method to be widely used in applications. The Eulerian Monte Carlo fields (EMCF) method has been developed to solve the PDF transport equation efficiently for decades. A recent work by Wang et al. revealed a serious issue of the EMCF method for not being fully consistent with the PDF transport equation for which the method is designed to solve. This work advances the state of the art by introducing fully consistent EMCF methods for solving the PDF transport equation. The fully consistent EMCF formulations are derived for two different PDF equation forms. The consistency of the EMCF formulations is mathematically confirmed by examining the derived moment transport equations from the EMCF formulations and from the PDF transport equation. The method of manufactured solutions is employed to further verify the consistency and convergence of the different EMCF formulations numerically. The newly introduced EMCF formulations bring the EMCF method to full consistency with the PDF transport equations for the first time.

42 ENGINEERING↗
DOE OSTI2020

Time and space resolved diagnostics for plasma thermal-chemical instability of fuel oxidation in nanosecond plasma discharges

Here, an instability in a nanosecond pulsed dielectric barrier discharge plasma occurring in methane–oxygen–argon mixtures is experimentally observed and measured by 1D time-resolved in situ electric field measurements. This instability, which seems to be created by the positive feedback between plasma kinetics and plasma-assisted low temperature fuel oxidation, is studied using electric field induced second harmonic generation and direct ICCD imaging. The rapid formation of streamers from an originally uniform discharge appears to be caused by the chemical kinetics of plasma-assisted low temperature methane oxidation, which may be resulting in a new type of plasma instability: a thermal-chemical instability. The results also revealed that the occurrence of this possible thermal-chemical instability in a reactive flow drastically changes the plasma properties by forming multiple secondary discharges and possibly leads to micron-sized non-uniform electric distributions. Single shot uncalibrated measurements of the electric field of the micron sized streamers appears to show much greater strengths than the average electric field. Furthermore, one-dimensional data analysis shows the positive feedback loop between the streamers and the low temperature plasma assisted oxidation chemistry in the plasma thermal-chemical instability. The present finding advances the understanding plasma instability growth and provides a new way to control plasma uniformity in plasma-assisted combustion and plasma fuel reforming.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

PeleLMeX [SWR-22-48]

PeleLMeX is a solver for high fidelity reactive flow simulations, namely direct numerical simulation (DNS) and large eddy simulation (LES). The solver combines a low Mach number approach, adaptive mesh refinement (AMR), embedded boundary (EB) geometry treatment and high performance computing (HPC) to provide a flexible tool to address research questions on platforms ranging from small workstations to the world's largest GPU-accelerated supercomputers. PeleLMeX has been used to study complex flame/turbulence interactions in RCCI engines and hydrogen combustion or the effect of sustainable aviation fuel on gas turbine combustion. PeleLMeX is part of the Pele combustion Suite (https://amrex-combustion.github.io/)

Day, Marcus↗
DOE OSTI2024

Non-resonant picosecond three-wave mixing in the gas phase

We report on the experimental observation of non-resonant, second-order optical sum-frequency generation (SFG) in five different atomic and molecular gases. The measured signal is attributed to a SFG process by characterizing its intensity scaling and its polarization behavior. We show that the electric quadrupole mechanism cannot explain the observed trends and suggest a mechanism based on symmetry breaking along the incident beam path arising from laser-induced species ground state number density gradients. Our results demonstrate that the SFG is about four orders of magnitude stronger than the third-harmonic generation (THG) and independent from any externally applied electric fields. In conclusion, these features make this method suitable for gas number density measurements at the picosecond time scale in reactive flows and plasmas.

LaCombe, Grayson↗
DOE OSTI2021

Simulating PBX 9501 Gap-stick Experiment

It has previously been shown that the rate parameters of the SURF reactive burn model can be adjusted to fit the gap-stick data for PBX 9501 [Johnson et al., 2018, fig 4]. A recent recalibration of the SURF model for PBX 9501 (lot 730-010 at ρ = 1.837 g/cc 3 ) used different reactants and products EOS and an updated fitting form for the burn rate [Menikoff, 2021]. The new calibration fit gap-stick data [Hill et al., 2021] together with data from SDT experiments and curvature effect data. Here the results of gap-stick simulations using the xRage code are shown. They provide insight into the reactive flow in the gap-stick experiment.

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

Stochastic Modeling in a Multimaterial Continuum Mixture Shock Physics Code

Stochastic modelling approaches are presented to capture random effects at multiple time and length scales. Random processes that occur at the microscale produce nondeterministic effects at the macroscale. Here we present three stochastic modeling approaches that describe random processes at microscopic length scales and map these processes to the macroscopic length scale. The first stochastic modeling approach is based upon a particle based numerical technique to solve a Stochastic Differential Equation (SDE) using an arbitrary diffusion process to capture random processes at the microstructural level. The second approach prescribes a Probability Density Function (PDF) for the drift and diffusion of the random variable derived using the forward and backward Kolmogorov equations. This method requires mean and drift evolution PDF transport equations. The third approach is the coupling of multiple random variables which are dependent on each other. The relationship of the PDFs and a coupling function, known as a copula, produces a Joint Probability Density Function (JPDF). These stochastic modeling approaches are implemented into a Multiple Component (MC) shock physics computational code and used to model statistical fracture and reactive flow applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Diffraction Measurement of Reaction Products in Shock Compressed TATB on the NIF

The shock-induced reaction of high explosives to gaseous and solid reaction products is a rapid, complex process. Understanding solid reaction product structure and formation kinetics is essential in determining the high-pressure equation of state of these multicomponent systems. We use the National Ignition Facility (NIF) to shock compress ~500-µm thick pressed powder high-explosive TATB samples to 70-135 GPa and collect in situ structural data on detonation byproducts using the TARDIS X-ray diffraction diagnostic. Velocimetry is used to record the transmitted compression wave profile, which is well described by Cheetah hydrocode simulations coupled with a reactive flow model, providing strong evidence of reaction in the TATB sample. While an unambiguous determination of the product phases was not possible owing to the low signal-to-noise quality of the diffraction signal, X-ray diffraction data of the product phases formed within the first 50 ns of the reaction is most consistent with a mixture of amorphous products and crystalline hexagonal diamond.

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

Understanding and modeling current and future coastal wetland methane dynamics (Final Report)

The coastal terrestrial-aquatic interface (TAI) is a highly dynamic component of the Earth system that plays a critical role in biogeochemical cycling. Due to its dynamic nature, the processes that regulate decomposition and methane (CH 4 ) emissions are of greater significance at the TAI than in upland systems. Despite this, we have limited mechanistic understanding of how climate stressors interact to regulate the electron acceptors and donors that determine decomposition pathways within TAIs, including the generation of hot spots and hot moments. Accurately modeling these processes is critical for incorporating the coastal TAI into Earth systems models, such as DOE’s E3SM. With previous DOE support, we adapted an aerobic terrestrial representation of decomposition using PFLOTRAN, a reactive flow and transport model, and added anerobic decomposition pathways, salinity, and oxygen (O 2 ) that fluctuates independently of water table level. However, because this model (PFLOTRAN TAI ) is based on decomposition rates and organic matter carbon to nitrogen ratios from terrestrial systems, its performance in TAI systems is currently limited by the lack of empirical data to properly parameterize variables. In addition, while PFLOTRAN TAI can simulate movement of O 2 into sediments, it is not currently capable of tracking the movement of CH 4 gas through plant tissues due to both current model structure and lack of available data. Our overall objective of the project was the increase our mechanistic understanding of CH 4 dynamics in response to environmental change, such that we can improve the representation of these dynamics in PFLOTRAN. We installed automated flux chambers in a new field-scale active soil warming experiment and set up marsh organs (mesocosms) to test effects of warming, flooding, and salinity. This resulted in a new dataset consisting of chamber-level CH 4 flux measurements across multiple sites, ecological conditions, and timeframes, as well as corresponding measurements on porewater chemistry, soil carbon quality, plant biomass, and redox reaction rates. Using these data, we improved PFLOTRAN TAI to more accurately model CH 4 dynamics and successfully tested our hypotheses. This grant contributed to the professional development of 3 postdocs, 4 undergraduate interns, 4 teacher interns, 17 technicians, and 15 participatory scientists. The automated chamber technology designed for this grant has also been shared with multiple new projects.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A Short Survey of Current Reactive Burn Model Capabilities as of Mid-Year 2025

The authors were tasked with writing a concise memo (now a short report) that adequately addresses the current ``state-of-the-art'' in the field of reactive flow modeling and burn models for mid-year 2025, along with identifying some of the modeling gaps. It is assumed that the reader has experience with burn models and running hydrocode simulations, and thus every effort is made to ensure brevity, clarity, and utility, so that this document will serve as a quick and useful reference. Note that “mid-year” was used in the title, as it is possible that advancements in artificial intelligence (AI) may significantly impact the field before the year is over.

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

Closure models for the feedback of energetic particles on plasma turbulence

Energetic particles interact with the plasma surrounding them, resonating with certain types of plasma waves to stabilize them while destabilizing others, and changing the character of the background turbulence in ways that have not been fully quantified or understood. Interaction with the turbulent background plasma is key to the acceleration of many types of energetic particles including high-energy cosmic rays, solar energetic particles, and pick-up ions. The acceleration of particles is a process that would ideally be described by a kinetic model, a type of model that follows a probability distribution function (PDF) for all particles in 7-dimensional (x, y, z, v x , v y , v z , t) space. Because of the high dimensionality of a kinetic model, simulations that solve kinetic equations use the largest computational resources currently available, and are yet unable to simulate a realistic number of particles, reach the large scales necessary for astrophysical problems, and use high-precision numerical methods. Two available alternatives to kinetic plasma models have been explored for this problem, with limited success. One is a multi-fluid model produced by a cumulant discarding closure, which evolves coupled equations for the velocity, magnetic field, and internal energy for both the background plasma and the fluid of energetic particles. However, simulations that solve multi-fluid magnetohydrodynamic (MHD) equations are able to include the interaction with energetic particles only in crude ways, typically as an add-on pressure term. The second alternative is to use a hybrid method to couple a fluid description of the background plasma to a kinetic model or a Fokker–Planck model for the energetic particles. These methods are hampered by the physical modeling of the coupling. In this work, we develop a new model, which follows the PDF for all particles; this can be viewed as a step toward physical realism above a multi-fluid MHD model, while also being more computationally efficient than a kinetic model. The equations we develop model both the background plasma and the energetic particles self-consistently. Over the last decade, similar PDF methods have been developed to a high level of sophistication to model reactive flows and turbulent combustion for engineering applications. For treatment of the feedback of the energetic particles on a background plasma, a PDF closure approach should evaluate the mean characteristics, including the density, with better statistical quality than will particle-sampling procedures.

79 ASTRONOMY AND ASTROPHYSICS↗