Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Flow chemistry”

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

Initial steady-state core simulation capability or thermal and pool-type molten salt reactors, coupling reactor physics, thermal-hydraulics, and evolving chemistry

This report presents the development and validation of an initial steady-state multiphysics capability for molten salt reactors (MSRs) under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in Fiscal Year 2025. The framework integrates neutronics, thermal-hydraulics, species transport, and thermochemistry to capture the coupled dynamics of liquid-fueled systems. Implementation and testing were performed on two representative designs: the Molten Salt Reactor Experiment (MSRE), a thermal-spectrum, channeled-flow reactor, and the Lotus Molten Salt Reactor (L-MSR), a fast-spectrum, pool-type reactor. The modeling suite employs Griffin for reactor physics and depletion, Pronghorn and SAM for thermal-hydraulics, Thermochimica for chemistry, and Saline for thermophysical properties, with benchmarking and validation carried out against historical MSRE data, experimental flow-loop measurements, and reference depletion calculations from Monte Carlo codes. The framework demonstrated the ability to reproduce key reactor behaviors including temperature feedback, reactivity losses, delayed neutron precursor transport, xenon poisoning, and redox potential evolution. The results confirm the feasibility and accuracy of the coupled models in predicting steady-state and selected transient MSR behaviors. This latter ones are used in this report as a proxy indicating that the steady-state models from which the transient starts are accurate. For MSRE, validation showed good agreement with pump start-up and natural circulation tests, while for the L-MSR, benchmarking confirmed hydraulic calibration and consistency of neutronics–thermal coupling. The tools also provided new insights into species transport, noble metal deposition, and salt solidification dynamics. On the Xenon transport front, the code is validated against the steady state Xenon poisoining measurement and showed good agreement with the experimental value. Identified areas for future work include advanced void transport modeling, three-dimensional simulations, improved alloy corrosion models, and tighter integration with high-fidelity Monte Carlo codes. These developments provide a foundation for high-fidelity MSR simulations that can support reactor design optimization, safety assessments, and long-term operational strategies.

42 - ENGINEERING↗
NASA NTRS2005

Axisymmetric Numerical Modeling of Pulse Detonation Rocket Engines

Pulse detonation rocket engines (PDREs) have generated research interest in recent years as a chemical propulsion system potentially offering improved performance and reduced complexity compared to conventional rocket engines. The detonative mode of combustion employed by these devices offers a thermodynamic advantage over the constant-pressure deflagrative combustion mode used in conventional rocket engines and gas turbines. However, while this theoretical advantage has spurred considerable interest in building PDRE devices, the unsteady blowdown process intrinsic to the PDRE has made realistic estimates of the actual propulsive performance problematic. The recent review article by Kailasanath highlights some of the progress that has been made in comparing the available experimental measurements with analytical and numerical models. In recent work by the author, a quasi-one-dimensional, finite rate chemistry CFD model was utilized to study the gasdynamics and performance characteristics of PDREs over a range of blowdown pressure ratios from 1-1000. Models of this type are computationally inexpensive, and enable first-order parametric studies of the effect of several nozzle and extension geometries on PDRE performance over a wide range of conditions. However, the quasi-one-dimensional approach is limited in that it cannot properly capture the multidimensional blast wave and flow expansion downstream of the PDRE, nor can it resolve nozzle flow separation if present. Moreover, the previous work was limited to single-pulse calculations. In this paper, an axisymmetric finite rate chemistry model is described and utilized to study these issues in greater detail. Example Mach number contour plots showing the multidimensional blast wave and nozzle exhaust plume are shown. The performance results are compared with the quasi-one-dimensional results from the previous paper. Both Euler and Navier-Stokes solutions are calculated in order to determine the effect of viscous effects in the nozzle flowfield. Additionally, comparisons of the model results to performance data from CalTech, as well as experimental flowfield measurements from Stanford University, are also reported.

Morris, Christopher I.↗
DOE OSTI2025

Assessment of flamelet/progress variable methods for supersonic combustion

Tabulated chemistry models, including the flamelet/progress variable approach, have been successfully used for a variety of turbulent flame simulations. The progress variable describes the progress of reactions in a system and parameterizes a lookup table of thermochemical variables. This approach reduces the cost of simulations, transporting only one scalar (progress variable) instead of the many species mass fractions required for detailed chemistry. Originally developed for low Mach number flame simulations, recent works have focused on extensions of this approach to compressible flames, supersonic combustion, and detonations, with applications such as scramjet combustors and rotating detonation engines. Unlike low Mach simulations, compressible flow simulations require solving the energy transport equation, which is coupled to the equation of state. This leads to additional modeling challenges regarding the thermodynamics and its impact on the chemistry. The validity of modeling assumptions, for example the relationship between energy and temperature, also varies with the combustion regime. The present work provides a detailed assessment of the existing strategies for chemistry tabulation for compressible/supersonic combustion, including detonations. A priori analysis indicates that approximations which are reasonable for weakly compressible flames may break down for shock-induced combustion. Furthermore, the analysis identifies specific assumptions and approximations that do not hold for detonations, emphasizing that care must be taken when applying tabulated chemistry models outside their intended combustion regimes.

Detonations↗
DOE OSTI2022

Mesoflow: An Open-Source Reacting Flow Solver for Catalysis at Mesoscale

We present the capabilities and software performance metrics of our open-source continuum solver for catalysis, Mesoflow, developed specifically for modeling transport and chemistry at the mesoscale. Our solver utilizes Cartesian block-structured adaptive mesh refinement to resolve complex catalyst surface morphologies directly obtained from X-ray tomography data. An immersed boundary based formulation enables rapid representation of complex geometries prevalent in most mesoporous catalyst interfaces. The solver is developed on top of open-source performance portable library, AMReX, providing parallel execution capabilities on current and upcoming high-performance-computing (HPC) architectures. Our flexible software framework enables integration of complex chemical mechanisms at heterogenous interfaces and time-split algorithms for circumventing highly disparate reaction and flow time-scales. Our current studies indicate a ten-fold performance gain by using graphics-processing-units (GPUs) compared to a single processor for representative problem sizes (2 million cell mesh). We will also present a brief introduction on how to build and use this software for application problems pertaining to catalytic upgrading and gas transport within porous catalyst particles.

adaptive meshing↗
NASA NTRS1993

Towards a generalized computational fluid dynamics technique for all Mach numbers

Currently there exists no single unified approach for efficiently and accurately solving computational fluid dynamics (CFD) problems across the Mach number regime, from truly low speed incompressible flows to hypersonic speeds. There are several CFD codes that have evolved into sophisticated prediction tools with a wide variety of features including multiblock capabilities, generalized chemistry and thermodynamics models among other features. However, as these codes evolve, the demand placed on the end user also increases simply because of the myriad of features that are incorporated into these codes. In order for a user to be able to solve a wide range of problems, several codes may be needed requiring the user to be familiar with the intricacies of each code and their rather complicated input files. Moreover, the cost of training users and maintaining several codes becomes prohibitive. The objective of the current work is to extend the compressible, characteristic-based, thermochemical nonequilibrium Navier-Stokes code GASP to very low speed flows and simultaneously improve convergence at all speeds. Before this work began, the practical speed range of GASP was Mach numbers on the order of 0.1 and higher. In addition, a number of new techniques have been developed for more accurate physical and numerical modeling. The primary focus has been on the development of optimal preconditioning techniques for the Euler and the Navier-Stokes equations with general finite-rate chemistry models and both equilibrium and nonequilibrium thermodynamics models. We began with the work of Van Leer, Lee, and Roe for inviscid, one-dimensional perfect gases and extended their approach to include three-dimensional reacting flows. The basic steps required to accomplish this task were a transformation to stream-aligned coordinates, the formulation of the preconditioning matrix, incorporation into both explicit and implicit temporal integration schemes, and modification of the numerical flux formulae. In addition, we improved the convergence rate of the implicit time integration schemes in GASP through the use of inner iteration strategies and the use of the GMRES (General Minimized Resisual) which belongs to the class of algorithms referred to as Krylov subspace iteration. Finally, we significantly improved the practical utility of GASP through the addition of mesh sequencing, a technique in which computations begin on a coarse grid and get interpolated onto successively finer grids. The fluid dynamic problems of interest to the propulsion community involve complex flow physics spanning different velocity regimes and possibly involving chemical reactions. This class of problems results in widely disparate time scales causing numerical stiffness. Even in the absence of chemical reactions, eigenvalue stiffness manifests itself at transonic and very low speed flows which can be quantified by the large condition number of the system and evidenced by slow convergence rates. This results in the need for thorough numerical analysis and subsequent implementation of sophisticated numerical techniques for these difficult yet practical problems. As a result of this work, we have been able to extend the range of applicability of compressible codes to very low speed inviscid flows (M = .001) and reacting flows.

Walters, R. W.↗
NASA NTRS2000

Kinetics of the Cl(2)P(J)) + CH4 Reaction: Effects of Secondary Chemistry Below 300 K

Absolute rate data for the Cl(2)P(J) + CH4 yields HCl + CH3 reaction have been obtained from 218 to 298 K by using the discharge-flow resonance fluorescence technique at I Torr total pressure. The result at 298 K is (10.1 +/- 0.6) x 10(exp -14) cu cm/molecule/s. The temperature dependence in Arrhenius form is (6.5 +/- 0.9 ) x 10(exp -12) exp[(-1235 +/- 34 )/T]. The errors given are one standard deviation; overall experimental error is estimated at +/- 15%. Because of the relatively large disagreement among earlier measurements at low temperatures, the results were examined for possible effects of non-Boltzmann spin distribution and vibrational excitation of CH4, secondary chemistry of CH3 radicals, and impurities in the CH4 source. There was no significant change in the observed rate constant when an efficient spin quencher, CF4, was added and estimates indicate that vibrational partitioning in CH4 should be at the ambient reactor temperature before the start of the reaction. The results were also independent of the source of Cl atoms (microwave discharge or thermal decomposition of Cl2) and whether CH4 was purified in-situ. However, the observed rate constant did depend on initial Cl atom concentrations and to a lesser extent on CH4 concentrations. Numerical simulations were used to assess the importance of secondary chemistry over a range of reactant concentrations

Wang, J. J.↗
NASA NTRS2011

Combustion Model of Supersonic Rocket Exhausts in an Entrained Flow Enclosure

This paper describes the Computation Fluid Dynamics (CFD) model developed to simulate the supersonic rocket exhaust in an entrained flow cylinder. The model can be used to study the plume-induced environment due to static firing test of the Taurus II launch vehicle. The finite rate chemistry is used to model the combustion process involving rocket propellant (RP 1) and liquid oxidizer (LOX). A similar chemical reacting model is also used to simulate the mixing of rocket plume and ambient air. The model provides detailed information on the gas concentration and other flow parameters within the enclosed region thus allowing different operating scenarios to be examined in an efficient manner. It is shown that the real gas influence is significant and yields better agreement with the theory.

Vu, Bruce↗
NASA NTRS1977

Shock-tube chemistry. I - The laminar-to-turbulent boundary layer transition

A model is proposed for calculating laminar-turbulent transition in the boundary layer of flows in chemical shock tubes. It is shown that the time to turbulent transition may be treated as the time required by the shocked gas to move over one shock-tube diameter. The model is shown to provide a means of predicting the turbulent transition time in any shock tube experiment, provided the condition of limiting flow is met. Knowing the transition time, it becomes possible to obtain the density and temperature profiles and, hence, to determine precisely the effect of shock tube boundary layers on the measurement of chemical rate constants.

Bander, J. A.↗
NASA NTRS1993

Upwind MacCormack Euler solver with non-equilibrium chemistry

A computer code, designated UMPIRE, is currently under development to solve the Euler equations in two dimensions with non-equilibrium chemistry. UMPIRE employs an explicit MacCormack algorithm with dissipation introduced via Roe's flux-difference split upwind method. The code also has the capability to employ a point-implicit methodology for flows where stiffness is introduced through the chemical source term. A technique consisting of diagonal sweeps across the computational domain from each corner is presented, which is used to reduce storage and execution requirements. Results depicting one dimensional shock tube flow for both calorically perfect gas and thermally perfect, dissociating nitrogen are presented to verify current capabilities of the program. Also, computational results from a chemical reactor vessel with no fluid dynamic effects are presented to check the chemistry capability and to verify the point implicit strategy.

Sherer, Scott E.↗
DOE OSTI2023

Paragenetic controls on CO 2 -fluid-rock interaction and weakening in a macroporous-dominated sandstone

The injection and storage of anthropogenic CO 2 in the subsurface is being deployed as a climate change mitigation tool; however, diagenetic-paragenetic heterogeneity in sandstone reservoirs often contributes to interval specific chemomechanical changes that affect injection and can increase leakage risk.In this paper we address reservoir heterogeneities’ impact on chemomechanical changes in a macroporous-dominated lithofacies of Morrow B sandstone, a formation containing several diagenetically-distinct hydraulic facies while undergoing enhanced oil recovery (EOR) and carbon dioxide (CO 2 ) sequestration. We performed three flow-through experiments using a CO 2 -charged or uncharged formation water combined with four indirect tensile strength tests per post-test sample. We then used the microstructure and paragenetic sequence to understand chemomechanical weakening with key observations as follows: dissolution of carbonates and feldspars changed porosity; increased permeability led to reclassifying each sample in a different hydraulic flow unit; decreased ultrasonic velocity; and did not lead to a loss of tensile strength. Tensile strength maintenance occurred due to the low abundance and minor dissolution of siderite, the stability of quartz, and the relative position of diagenetic ankerite within feldspar. This macroporous-dominated lithofacies is the primary reservoir for the Morrow B Sandstone, and is analogous to other porous sandstone reservoirs. It represents an end-member of a chemomechanically low-risk siliceous CO 2 sequestration and CO 2 -EOR reservoir.

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

Incompressible spatially-developing free-shear flows

The first goal of this research is to understand the dynamics of three dimensional vortical structures, especially the interaction of the ribs and two dimensional rollers. Because of vortex stretching effects, this requires accurate and well resolved simulations. Another goal is to distinguish temporally- and spatially-growing shear layers in terms of the dynamics and observed vortical structures. This is important since the former is much easier to compute, but the error incurred in using temporal results for certain applications is unknown. The third goal is to investigate the effects of asymmetry on passive scalar mixing and fast-chemistry product formation. In this report a brief summary of a new numerical method is presented, followed by results for mixing layers and wake flows.

Buell, J. C.↗
DOE OSTI2025

Bayesian calibration of bubble size dynamics applied to CO 2 gas fermenters

To accelerate the scale-up of gaseous CO 2 fermentation reactors, computational models need to predict gas-to-liquid mass transfer which requires capturing the bubble size dynamics, i.e. bubble breakup and coalescence. However, the applicability of existing models beyond air–water mixtures remains to be established. Here, an inverse modeling approach, accelerated with a neural network surrogate, calibrates the breakup and coalescence closure models, that are used in class methods for population balance modeling (PBM). The calibration is performed based on experimental results obtained in a CO 2 -air–water-coflowing bubble column reactor. Bayesian inference is used to account for noise in the experimental dataset and bias in the simulation results. To accurately capture gas holdup and interphase mass transfer, the results show that the breakage rate needs to be increased by one order of magnitude. In conclusion, the inferred model parameters are then used on a separate configuration and shown to also improve bubble size distribution predictions.

09 BIOMASS FUELS↗
DOE OSTI2021

Towards chemical equilibrium in thermochemical water splitting. Part 1: Thermal reduction

The efficiency of many processes strongly depends on their thermodynamic reversibility, i.e., proximity to equilibrium throughout the process. In thermochemical cycles for water and/or carbon dioxide splitting, thermochemical air separation, and thermochemical energy storage, operating near equilibrium means that the oxygen chemical potential of the solid and gas phases must not differ significantly. Herein we show that approaching this ideal is possible in thermal reduction only if the reaction step occurs at a specific, reaction coordinate- and material-dependent temperature. The resulting thermal reduction temperature profile also depends on the ratio of gas and solid flows.

08 HYDROGEN↗
DOE OSTI2023

An effective air–liquid contactor for CO 2 direct air capture using aqueous solvents

Here, the development of a cost-effective, corrosion resistant, high-flux direct air capture (HiDAC) contactor for a solvent-based direct air capture (DAC) process is reported. Literature technoeconomic analyses suggest that the air–liquid contactor can cost over 20 % of the overall DAC plant’s annualized capital costs. To bring down the overall cost of DAC, it is imperative that an effective contactor is developed. A hybrid contactor consisting of a commercial polyvinylchloride structured packing enhanced with stainless-steel 410 random packing has been developed to provide a high surface area for air–liquid contact. The contactor geometry, wettability, corrosion resistance, pressure drop, along with its CO 2 uptake efficiency, CO 2 uptake rates, and extended loading potential using potassium sarcosinate solutions are investigated. Results show that the HiDAC contactor has a relatively high specific surface area (885 m 2 /m 3 ), which allows for CO 2 uptake efficiencies of up to 75 % and capture rates of up to 550 g of CO 2 per day for a 0.3 × 0.25 × 0.3 m contactor. The contactor also exhibits high levels of wettability and corrosion resistance with amino acid-based DAC solvents. A CO 2 uptake model was developed, and modeling results are compared to experimental data to simulate and predict the performance of the contactor in a DAC process using potassium sarcosinate solvent. The system’s overall mass transfer coefficient and theoretical pressure drop were also calculated, and these results were compared to DAC data found in literature. The results presented indicate that the HiDAC contactor is well suited for DAC due to its high specific surface area, resistance to corrosion, and high degree of wettability.

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

Systems Analysis Approach to Polyethylene Terephthalate and Olefin Plastics Supply Chains in the Circular Economy: A Review of Data Sets and Models

The environmental and economic impacts of implementing a circular economy in plastic waste supply chains are not well understood. The proposed systems analysis framework assesses environmental, social, and economic impacts of plastic waste supply chains in a circular economy. The first objective of this article is to identify datasets, models and knowledge gaps associated with waste plastic supply chain processes, mainly in the U.S. Our literature review indicated that the best datasets exist for virgin plastic resin production, mechanical recycling, landfilling, and incineration, with materials recovery facility being intermediate, and with chemical recycling the lowest. The second objective of this article is to develop an illustrative application of the framework by conducting a preliminary systems analysis of PET bottles with closed-loop recycling. Here, the preliminary systems analysis of PET bottles utilized a linear programming optimization method. Our optimization model indicated that both chemical and mechanical recycling processes are needed to achieve a true circular economy of PET bottles with the least greenhouse gas emissions, specifically reductions of 24% when compared with the linear economy. Good quality and standardized life cycle assessment and techno-economic analysis studies are needed to better understand the environmental, economic, and social impacts of advanced sorting and chemical recycling technologies.

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

Laser-spectroscopy testbed for impurity monitoring in liquid metal-cooled fast reactors

A significant challenge in sodium-cooled fast reactors is controlling impurities, in particular oxygen impurities, within the sodium coolant, as they can accelerate corrosion and indicate leaks. Optical methods offer the potential to rapidly detect small concentrations of both gaseous and metal impurities that accelerate corrosion, plug coolant channels, and lead to increased activation of isotopes in the coolant. We present the design and performance of an apparatus designed to enable the application of multiple optical analytical techniques, such as laser-induced breakdown spectroscopy, to detect elemental impurities in the sodium melt with high sensitivity. Here, we experimentally demonstrate the detection of characteristic sodium and oxygen spectral lines in liquid sodium, which sets the stage for the optimization of its analytical sensitivity. A robust sensor of this type integrated with the sodium cooling loop has the potential to significantly improve the safety and operational efficiency of generation IV nuclear reactors.

47 OTHER INSTRUMENTATION↗