Search NASA⌕ Search

SEARCH · Search NASA

Results for “reaction front”

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

Relating the depth of the water table to the depth of weathering

Weathering of bedrock creates and occludes permeability, affecting subsurface water flow. Often, weathering intensifies above the water table. On the other hand, weathering can also commence below the water table. To explore relationships between weathering and the water table, a simplified weathering model for an eroding hillslope was formulated that takes into account both saturated and unsaturated subsurface water flow (but does not fully account for changes in dissolved gas chemistry). The phreatic line was calculated using solutions to mathematical treatments for both zones. In the model, the infiltration rate at the hill surface sets both the original and the eventual steady-state position of the water table with respect to the weathering reaction front. Depending upon parameters, the weathering front can locate either above or below the water table at steady state. Erosion also affects the water table position by changing porosity and permeability even when other hydrological conditions (e.g. hydraulic conductivity of parent material, infiltration rate at the surface) do not change. The total porosity in a hill (water storage capacity) was found to increase with infiltration rate (all else held constant). Furthermore, this effect was diminished by increasing the erosion rate. We also show examples of how the infiltration rate affects the position of the water table and how infiltration rate affects weathering advance.

58 GEOSCIENCES↗
DOE OSTI2026

Reactive Transport Modeling with Physics-Informed Machine Learning for Critical Minerals Applications

This study presents a physics-informed neural network (PINN) framework for reactive transport modeling for simulating fast bimolecular reactions in porous media. Accurate characterization of cAhemical interactions and product formation in surface and subsurface environments is essential for advancing critical mineral extraction and related geoscience applications. The proposed methodology sequentially addresses the flow and diffusion–reaction subproblems. The flow field is computed using a mixed formulation, while the diffusion–reaction system is modeled via two uncoupled tensorial diffusion equations reformulated in terms of chemical invariants. PINNs are employed to solve the governing equations, enabling data-efficient, mesh-free prediction of chemical concentration fields. The framework is validated through a series of benchmark problems involving flow in heterogeneous porous media. Initial verification is conducted using patch tests for the flow field, followed by validation of the transport problem with emphasis on preserving non-negativity of concentrations. The complete fast bimolecular reaction scenario is then solved, yielding spatial distributions of reactants and product species. Results demonstrate that the PINNs-based approach effectively captures sharp, mixing-limited reaction fronts and dispersive mixing behavior, offering reliable predictions of reactive plume evolution. These capabilities are crucial for evaluating long-term subsurface behavior in applications such as fluid storage, energy extraction, and efficient extraction of critical minerals.

42 ENGINEERING↗
DOE OSTI2026

A Generalized Grain-Scale Model for the Non-Plasma and Plasma-Assisted Hydrogen Direct Reduction of Iron Ore

Direct Reduction of Iron ore using hydrogen (H-DRI) is a promising pathway towards efficient steelmaking and accurate predictive models are a necessity for scale-up and optimization of this technology. However, accurate models of this process remain limited because existing models oversimplify grain-scale phenomena, such as nonlinearity inside grain, self-sufficient porosity, surface reactions, and the role of plasma species. These phenomena are important for flash steelmaking and plasma-assisted H-DRI processes. To address this need, we present a phenomenological model for simulating H-DRI at the scale of a single micron-sized grain of the iron ore. We call this the Transient Reactive Grain Model (TRGM). TRGM incorporates key physical process: gas species transport, a chemical kinetics of material conversion, nanopore structural evolution and, adsorption-desorption surface kinetics at the reactive nanopore surface. The important contribution of this work is that the model provides a dependence on different reductant species, specifically hydrogen atoms versus molecules, so that role of hydrogen plasma reduction can be clarified compared to the use of pure hydrogen gas reduction. TRGM predictions agree well with experimental data for both molecular H2 reduction of Fe2O3 and plasma hydrogen reduction of Fe3O4. Results reveal species concentration gradients with a diffuse reaction zone, and enhanced hydrogen diffusion at the grain outer surface due to evolving porosity. These findings challenge common assumptions in existing models, including sharp reaction fronts, quasi-steady diffusion and kinetics, and the neglect of surface chemistry. As a generalized grain-scale model for H-DRI processes, TRGM has practical applications in flash steelmaking and in-flight reduction using both molecular and plasma hydrogen.

08 HYDROGEN↗
DOE OSTI2024

Direct numerical simulation of a high-pressure hydrogen micromix combustor: Flame structure and stabilisation mechanism

A high-pressure hydrogen micromix combustor has been investigated using direct numerical simulation with detailed chemistry to examine the flame structure and stabilisation mechanism. The configuration of the combustor was based on the design by Schefer et al., using numerical periodicity to mimic a large square array. A precursor simulation of an opposed jet-in-crossflow was first conducted to generate appropriate partially-premixed inflow boundary conditions for the subsequent reacting simulation. The resulting flame can be described as an predominantly-lean inhomogeneously-premixed lifted jet flame. Five main zones were identified: a jet mixing region, a core flame, a peripheral flame, a recirculation zone, and combustion products. The core flame, situated over the jet mixing region, was found to burn as a thin reaction front, responsible for over 85% of the total fuel consumption. The peripheral flame shrouded the core flame, had low mean flow with high turbulence, and burned at very lean conditions (in the distributed burning regime). It was shown that turbulent premixed flame propagation was an order-of-magnitude too slow to stabilise the flame at these conditions. Stabilisation was identified to be due to ignition events resulting from turbulent mixing of fuel from the jet into mean recirculation of very lean hot products. Ignition events were found to correlate with shear-driven Kelvin-Helmholtz vortices, and increased in likelihood with streamwise distance. At the flame base, isolated events were observed, which developed into rapidly burning flame kernels that were blown downstream. Further downstream, near-simultaneous spatially-distributed ignition events were observed, which appeared more like ignition sheets. The paper concludes with a broader discussion that considers generalising from the conditions considered here.

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

Chemical timescale effects on detonation convergence

Numerical simulations of detonation-containing flows have emerged as crucial tools for designing next- generation power and propulsion devices. As these tools mature, it is important for the combustion community to properly understand and isolate grid resolution effects when simulating detonations. To this end, the objective of this work is to provide a comprehensive analysis of the numerical convergence of unsteady detonation simulations, with focus on isolating the impacts of chemical timescale modifications on convergence characteristics in the context of operator splitting. With the aid of an AMReX-based adaptive mesh refinement flow solver-which enables resolutions up to ($\mathcal{O}$ (1000) cells-per-induction length-the convergence analysis is conducted using two kinetics configurations: (1) the simplified three-step Arrhenius-based model mechanism of Short and Quirk (1997), where chemical timescales in the detonation are modified by adjusting activation energies in the initiation and branching reactions, and (2) a detailed hydrogen- air mechanism, where the chemical timescales are adjusted by varying the ambient pressure. The convergence of unsteady self-sustained detonations in one-dimensional channels is then analyzed with reference to steady-state theoretical baseline solutions using these mechanisms. The goal of the analysis is to provide a detailed comparison of the effects of grid resolution on both macroscopic (peak pressures and wave speeds) and microscopic (wave structure) quantities of interest, drawing connections between the deviations from steady-state baselines and minimum chemical timescales. In particular, chemical timescale reductions were found to have minimal impact on the convergence of macroscopic properties. Furthermore, analyses of microscopic convergence trends, particularly in the reaction front location, revealed a key insight: maintaining the induction time while eliminating prohibitive chemical timescales through mechanism simplifications and combustion modeling can significantly enhance detonation convergence properties. Ultimately, this work uncovers resolution-dependent unsteady detonation convergence regimes and highlights the important role played by not only the chemical timescales, but also the ratio between the chemical timescale and induction time on the numerical convergence of the detonation wave structure.

Adaptive mesh refinement↗
DOE OSTI2020

Real Time 3D Observations of Portland Cement Carbonation at CO2 Storage Conditions

Depleted oil reservoirs are considered a viable solution to the global challenge of CO 2 storage. A key concern is whether the wells can be suitably sealed with cement to hinder the escape of CO 2 . Under reservoir conditions, CO 2 is in its supercritical state, and the high pressures and temperatures involved make real-time microscopic observations of cement degradation experimentally challenging. Here, we present an in situ 3D dynamic X-ray micro computed tomography (μ-CT) study of well cement carbonation at realistic reservoir stress, pore-pressure, and temperature conditions. The high-resolution time-lapse 3D images allow monitoring the progress of reaction fronts in Portland cement, including density changes, sample deformation, and mineral precipitation and dissolution. By switching between flow and nonflow conditions of CO 2 -saturated water through cement, we were able to delineate regimes dominated by calcium carbonate precipitation and dissolution. For the first time, we demonstrate experimentally the impact of the flow history on CO 2 leakage risk for cement plugging. In-situ μ-CT experiments combined with geochemical modeling provide unique insight into the interactions between CO 2 and cement, potentially helping in assessing the risks of CO 2 storage in geological reservoirs.

58 GEOSCIENCES↗
DOE OSTI2019

Moisture Corrosion of LiH: A Kinetic Investigation by DRIFT Spectroscopy

Lithium hydride (LiH) is a unique, ionic compound with applications in a variety of industries. Unfortunately, LiH is very reactive toward H 2 O even at ppm levels, forming oxide (Li 2 O) and hydroxide (LiOH) corrosion layers while outgassing H 2 . An effective means to eliminate unwanted outgassing is vacuum-heating to convert LiOH into Li 2 O, although subsequent re-exposure to moisture during transport/handling reconverts some Li 2 O back to LiOH. A corrosion growth model for previously vacuum-baked LiH is necessary for long-term prediction of the hydrolysis of LiH. In this work, a para-linear hydroxide corrosion growth model is proposed for the reaction of previously vacuum-baked LiH samples with moisture. This model, composed of two competing diffusion reaction fronts at the LiOH/Li 2 O and Li 2 O/LiH interfaces, is validated experimentally by subjecting a previously vacuum-baked polycrystalline LiH sample to 35 ppm of H 2 O at room temperature while monitoring the corrosion growth as a function of time with diffuse-reflectance infrared Fourier transform (DRIFT) spectroscopy. The para-linear growth model for the hydrolysis of previously vacuum-baked LiH proposed in this report can also serve as a template for the hydrolysis of other hygroscopic oxides grown on metal or metal hydride substrates.

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

Revealing High-Temperature Reduction Dynamics of High-Entropy Alloy Nanoparticles via In Situ Transmission Electron Microscopy

Understanding the behavior of high-entropy alloy (HEA) materials under hydrogen (H 2 ) environment is of utmost importance for their promising applications in structural materials, catalysis, and energy-related reactions. Herein, the reduction behavior of oxidized FeCoNiCuPt HEA nanoparticles (NPs) in atmospheric pressure H-2 environment was investigated by in situ gas-cell transmission electron microscopy (TEM). The reduction reaction front was maintained at the external surface of the oxide. During reduction, the oxide layer expanded and transformed into porous structures where oxidized Cu was fully reduced to Cu NPs while Fe, Co, and Ni remained in the oxidized form. In situ chemical analysis showed that the expansion of the oxide layer resulted from the outward diffusion flux of all transition metals (Fe, Co, Ni, Cu). Revealing the H-2 reduction behavior of HEA NPs facilitates the development of advanced multicomponent alloys for applications targeting H 2 formation and storage, catalytic hydrogenation, and corrosion removal.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Replacement of Calcium Carbonate Polymorphs by Cerussite

Calcium carbonate (CaCO 3 ) polymorphs, calcite, aragonite, and vaterite, serve as a major sink to retain various metal ions in natural and engineered systems. Here, we visualize the systematic trends in reactivities of calcite, vaterite, and aragonite to Pb 2+ dissolved in acidic aqueous solutions using in situ optical microscopy combined with ex situ scanning electron and transmission X-ray microscopies. All three polymorphs undergo pseudomorphic replacement by cerussite (PbCO 3 ) but with distinct differences in the evolution of their morphologies. The replacement of calcite and aragonite occurs through the formation of a pseudomorphic cerussite shell (typically 5–10 μm thick) followed by a slower inward propagation of reaction fronts through a thin solution gap (~0.1 μm wide) between the shell and the CaCO 3 substrate. The replacement of vaterite is characterized by the formation of a thinner cerussite shell (≤1 μm thick) and a larger cavity between the shell and the host mineral. These systematic differences in cerussite morphology for different CaCO 3 polymorphs are explained by the relative dissolution and precipitation rates of the reactant and product minerals, coupled with the role of ion transport through the cerussite shells. We also find that the replacement of calcite by cerussite is the slowest when all three polymorphs coexisted. Furthermore, our results provide mechanistic insights into the growth mode of cerussite on dissolving calcium carbonate and demonstrate these CaCO 3 polymorphs as promising substrate materials for removal and recycling of Pb from acidic polluted water and industrial effluents.

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

Operando Neutron Imaging of Reaction Extent and Particle Swelling Informs Limiting Factors for Salt Hydrate Thermochemical Energy Storage

Salt hydrates are a promising thermochemical energy storage medium that stores heat through the reversible uptake (hydration) and release (dehydration) of water vapor. Our study deploys operando neutron imaging to investigate salt hydrate performance with high spatial resolution (42 μm pixels). For flow over a packed bed with diffusion-driven transport, measurements reveal the formation of a solid diffusion layer due to particle swelling for the pure SrBr2 salt. In contrast, the SrBr2–vermiculite composite exhibits significantly less swelling and more than a 2-fold increase in the apparent water vapor diffusivity. For axial flow through a packed bed, neutron imaging confirms theoretically predicted transitions from a moving reaction front to a homogeneous profile with an increase in humid air flow rate. Our study establishes neutron imaging as a powerful technique to advance fundamental understanding of thermochemical systems and help guide composite material design.

Kinzer, Bryan [ORNL] (ORCID:0000000337804910)↗
DOE OSTI2023

Asynchronous domain dynamics and equilibration in layered oxide battery cathode

To improve lithium-ion battery technology, it is essential to probe and comprehend the microscopic dynamic processes that occur in a real-world composite electrode under operating conditions. The primary and secondary particles are the structural building blocks of battery cathode electrodes. Their dynamic inconsistency has profound but not well-understood impacts. In this research, we combine operando coherent multi-crystal diffraction and optical microscopy to examine the chemical dynamics in local domains of layered oxide cathode. Our results not only pinpoint the asynchronicity of the lithium (de)intercalation at the sub-particle level, but also reveal sophisticated diffusion kinetics and reaction patterns, involving various localized processes, e.g., chemical onset, reaction front propagation, domains equilibration, particle deformation and motion. These observations shed new lights onto the activation and degradation mechanisms of state-of-the-art battery cathode materials.

25 ENERGY STORAGE↗
DOE OSTI2023

Grain boundary widening controls siderite (FeCO3) replacement of limestone (CaCO3)

Abstract The microstructure of minerals and rocks can significantly alter reaction rates. This study focuses on identifying transport paths in low porosity rocks based on the hypothesis that grain boundary widening accelerates reactions in which one mineral is replaced by another (replacement reaction). We conducted a time series of replacement experiments of three limestones (CaCO 3 ) of different microstructures and solid impurity contents using FeCl 2 . Reacted solids were analyzed using chemical imaging, small angle X-ray and neutron scattering and Raman spectroscopy. In high porosity limestones replacement is reaction controlled and complete replacement was observed within 2 days. In low porosity limestones that contain 1–2% dolomite impurities and are dominated by grain boundaries, a reaction rim was observed whose width did not change with reaction time. Siderite (FeCO 3 ) nucleation was observed in all parts of the rock cores indicating the percolation of the solution throughout the complete core. Dolomite impurities were identified to act as nucleation sites leading to growth of crystals that exert force on the CaCO 3 grains. Widening of grain boundaries beyond what is expected based on dissolution and thermal grain expansion was observed in the low porosity marble containing dolomite impurities. This leads to a self-perpetuating cycle of grain boundary widening and reaction acceleration instead of reaction front propagation.

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

Operando neutron radiography validates a parameter-free transport–kinetics model for thick solid-state battery cathodes

Tortuosity-weighted interfacial flux for lithium (TWIF-Li) predicts through-thickness Li gradients in thick composite all-solid-state cathodes without fitted parameters. Image-derived microstructures, GITT-derived concentration-dependent solid diffusion, and tortuosity-weighted interfacial kinetics reproduce operando neutron radiography across practical rates, delivering transferable design rules to suppress transport-limited reaction fronts.

Adam, Andre [ORNL] (ORCID:0000000245023033)↗
DOE OSTI2025

Two-dimensional mapping of absolute OH densities in an atmospheric pressure plasma effluent via planar laser-induced fluorescence: effects of He/H 2 O and He/O 2 mixtures in N 2 and air, with and without solid targets

Planar laser-induced fluorescence (LIF) was employed to measure the absolute density of hydroxyl radicals (OH) in the effluent of the COST Reference Microplasma Jet for two feed gas mixtures: He/H 2 O and He/O 2 . Experiments were conducted with the effluent propagating into air and N 2 environments. For the He/H 2 O case, measurements were also performed with the effluent impinging on a solid target at varying distances from the jet nozzle. Calibration of the OH-LIF signal from the COST-Jet was achieved by comparing it to a reference signal generated by the photofragmentation of H 2 O 2 . Results demonstrated that OH densities were sustained longer when the effluent propagates in a nitrogen environment compared to air, particularly with water added to the feed gas. The broader OH distribution in N 2 suggests slower consumption due to the absence of oxygen, which accelerates OH depletion in air via reactions involving O 2 and HO 2 . Even when water was not added to the feed, as in the He/O 2 case, appreciable OH densities were observed, due to gas impurities and reactive species interactions with atmospheric humidity, forming reaction fronts that delineate the gas flow. Two-dimensional fluid dynamics simulations elucidated the influence of atmospheric gas entrainment and solid targets on the OH distribution. Experimental trends were further compared with a zero-dimensional chemistry model to explore OH production and consumption mechanisms in air and nitrogen environments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Impact of aqueous solution pH on network structure of corrosion–induced surface layers of boroaluminosilicate glass

Understanding of the long-term stability of vitrified nuclear waste glass is critical to its safe storage, but knowledge gaps remain regarding the corrosion mechanism, especially with respect to surface alteration layers whose presence changes the observed corrosion rates. Here, we report the aqueous corrosion of a model boroaluminosilicate glass at ground water relevant pHs as investigated by vibrational spectroscopy to understand the network structure of the alteration layer that we hypothesize controls the corrosion rates. Both specular reflection infrared spectroscopy and infrared variable angle spectroscopic ellipsometry were used to elucidate details of the structure of the alteration layers formed at pH 7 and pH 9. After 2312 days, the thickness of the alteration layer is found to be ~2000 nm for corrosion at pH 7 and ~300 nm for corrosion at pH 9. Here, the interpretation of the vibrational spectra suggests that alteration layers formed at pH 7 and pH 9 both exhibit a porous network with "silica-like" structures, but with differences in the bond lengths and bond angles in the glass network. Moreover, the relative abundance of the hydrous species within the porous network appears to be dependent on pH. Additionally, the spectroscopic data suggest that the B and Ca trapped in the alteration layer during the secondary corrosion of the pH 9 pre-corroded sample in pH 3 solution are hydrated within the nano-porous layer and not bound to the silicate network. However, the precise chemical form remains to be determined. Based on these findings, we propose that the combination of the network structure of the alteration layer and the relative abundance of hydrous species within the alteration layer act to control the transport of ionic species dissolved from the reaction front to the aqueous solution outside the glass.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Geochemical Impact of Acid Spearhead and Slickwater Stimulation on Wolfcamp Shale from the Hydraulic Fracturing Test Site

The Hydraulic Fracturing Test Site 1 (HFTS-1) was a field study performed in the Wolfcamp Formation in the West Texas Permian (Midland) Basin, USA, with a focus on improving the efficiency of hydraulic fracturing. Investigating site-specific rock-fluid geochemical interactions during hydraulic fracturing is an important step to understanding the impact on formation shale porosity, permeability, and long-term shale gas production. During field operations in this region, hydraulic fracturing fluid (HFF) injection usually starts with a concentrated acid spearhead for rapid rock dissolution, followed by the injection of near-neutral pH slickwater containing chemicals and proppants. A multistep sequential injection approach was used to investigate different stages of rock-fluid interactions. The carbonate content in the host rock is important when acid spearhead is considered, as carbonate mineral dissolution is rapid and can result in porosity and permeability changes in the shale matrix. Here, in this study, we designed flow-through experiments using fractured carbonate-rich and clay-rich Wolfcamp shale cores with (1) a short-time acid soaking step and (2) a long-term slickwater flow-through step to simulate the injection method used at HFTS-1. The fluid chemistry was analyzed. A thorough mineralogical progression [e.g., Calcium (Ca) dissolution and iron (Fe) redox progression] in the cores during HFF injection was also characterized and imaged by synchrotron microprobe. Reactive transport modeling was performed based on the experimental setup. The results showed that the acid spearhead is a crucial step in creating a reaction front by mineral dissolution, especially in carbonate-rich shales. A slight layer of ferrihydrite precipitated during the slickwater flow-through period. This study provides insights into potential geochemical impact due to hydraulic fracturing operations in the Permian Basin.

clastic rock↗
DOE OSTI2020

Exploring an ‘ideal hill': how lithology and transport mechanisms affect the possibility of a steady state during weathering and erosion

We introduce a model of chemical reaction within hills to explore how evolving porosity (and by inference, permeability) affects flow pathways and weathering. The model consists of hydrologic and reactive-transport equations that describe alteration of ferrous minerals and feldspar. These reactions were chosen because previous work emphasized that oxygen- and acid-driven weathering affects porosity differently in felsic and mafic rocks. A parameter controlling the order of the fronts is presented. In the absence of erosion, the two reaction fronts move at different velocities and the relative depths depend on geochemical conditions and starting composition. In turn, the fronts and associated changes in porosity drastically affect both the vertical and lateral velocities of water flow. For these instances, estimates of weathering advance rates based on simple models that posit unidirectional constant-velocity advection do not apply. In the model hills, weathering advance rates diminish with time as the Darcy velocities decrease with depth. The system can thus attain a dynamical steady state at any erosion rate where the regolith thickness is constant in time and velocities of both fronts become equal to one another and to the erosion rate. The slower the advection velocities in a system, the faster it attains a steady state. For example, a massive rock with relatively fast-dissolving minerals such as diabase -- where solute transport across the reaction front mainly occurs by diffusion -- can reach a steady state more quickly than granitoid rocks in which advection contributes to solute transport. The attainment of a steady state is controlled by coupling between weathering and hydrologic processes that force water to flow horizontally above reaction fronts where permeability changes rapidly with depth and acts as a partial barrier to fluid flow.

54 ENVIRONMENTAL SCIENCES↗