Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reacting 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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Feasibility of using reactive silicate particles with temperature-responsive coatings to enhance the security of geologic carbon storage

The large-scale deployment of geologic carbon storage will require minimizing the risk of CO 2 leakage. Here we report on a novel strategy to seal leakage pathways by injecting functionalized reactive calcium silicate (CaSiO 3 ) particles. The particles consist of a mineral silicate core within a temperature-sensitive Poly(N,N-dimethylaminoethyl methacrylate) coating that limits contact between the core and the surrounding CO 2 (aq) above a switching temperature making them relatively unreactive deep in the subsurface. If the CO 2 leaks to shallower depths and cooler temperatures, the coatings would extend and enable the silicate core to react with CO 2(aq) to form solid carbonate precipitates that reduce the permeability of the fracture or flow path. A combined laboratory- and simulation-based investigation was conducted to explore key early-stage research questions related to the feasibility of this leakage management strategy. Sand columns were injected with either coated or uncoated CaSiO 3 particles and the permeability was measured before and after exposure to CO 2 . Experiments were conducted at 35 °C and 70 °C, a range that spans the switching temperature for the coatings. Results show that at the cooler temperature the coating effectively enabled reaction and permeability reduction. At higher temperatures, the coated particles limited mineral carbonation and the permeability decrease was minimal compared to columns with uncoated particles. Here, the morphology and location of the precipitates provides insight about permeability evolution. Model simulations of 2D radial particle transport from an injection point provided insight into the effects of particle size, coating properties, and concentrations of particles in several representative rock formations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Microwave-Assisted Production of Polycarbonate Diols using Carbon Dioxide

Carbon dioxide (CO2) is a cheap and readily available resource that can be converted to value-added chemicals including fuels, polymers, and other products. One lucrative option is polycarbonate diols, which are the precursor to polyurethanes that find wide applications in automotive and aerospace industries. The traditional process for the synthesis of polycarbonate diols involves using hazardous phosgene and alcohols in a strongly basic medium that generates significant amounts of salt. Alternately, CO2 can be reacted with alkanediols over CeO2 catalysts to produce polycarbonate diols. For the current work, a microwave-assisted atmospheric flow system was investigated for converting CO2 into polycarbonate diols and was compared against a thermal system. Various parameters such as type and ratio of solvents, and the addition of dehydrating agents for the microwave system was tested. The different CeO2 catalysts used were characterized using Brunauer–Emmett–Teller (BET) surface area analysis, Raman spectroscopy, X-ray diffraction, and temperature programmed desorption (NH3-TPD and CO2-TPD).

CO2 utilization↗
DOE OSTI2024

CH 3 Radical Generation in Microplasmas for Up-Conversion of Methane

The conversion of methane, CH 4 , into higher value chemicals using low temperature plasmas is challenged by both improving efficiency and selectivity. One path towards selectivity is capturing plasma produced methyl radicals, CH 3 , in a solvent for aqueous processing. Due to the rapid reactions of methyl radicals in the gas phase, the transport distance from production of the CH 3 to its solvation should be short, which then motivates the use of microplasmas. The generation of CH 3 in Ar/CH 4 /H 2 O plasmas produced in nanosecond pulsed dielectric barrier discharge microplasmas is discussed using results from a computational investigation. The microplasma is sustained in the channel of a microfluidic chip in which the solvent flows along one wall or in droplets. CH 3 is primarily produced by electron-impact of and dissociative excitation transfer to CH 4 , as well as CH 2 reacting with CH 4 . CH 3 is rapidly consumed to form C 2 H 6 which, in spite of being subject to these same dissociative processes, accumulates over time, as do other stable products including C 3 H 8 and CH 3 OH. The gas mixture and electrical properties were varied to assess their effects on CH 3 production. CH 3 production is largest with 5% CH 4 in the Ar/CH 4 /H 2 O mixture due to an optimal balance of electron-impact dissociation, which increases with CH 4 percentage, and dissociative excitation transfer and CH 2 reacting with CH 4 , which decrease with CH 4 percentage. Design parameters of the microchannels were also investigated. Increasing the permittivity of the dielectrics in contact with the plasma increased the ionization wave intensity which increased CH 3 production. Increased energy deposition per pulse generally increased CH 3 production as does lengthening pulse length up to a certain point. The arrangement of the solvent flow in the microchannel can also affect the CH 3 density and fluence to the solvent. The fluence of CH 3 to the liquid solvent is increased if the liquid is immersed in the plasma as a droplet or is a layer on the wall where the ionization wave terminates. The solvation dynamics of CH 3 with varying numbers of droplets was also examined. Here, the maximum density of solvated methyl radicals CH 3aq occurs with a large number of droplets in the plasma. However, the solvated CH 3aq density can rapidly decrease due to desolvation, emphasizing the need to quickly react the solvated species in the solvent.

03 NATURAL GAS↗
DOE OSTI2025

Dynamic response of a freely rotating butterfly valve in the advanced test reactor − dynamic coefficients modeling

Here, in evaluating the water hammer issue pertaining to the primary-coolant-regulating butterfly valve in the Advanced Test Reactor, the dynamic fluid body interaction (DFBI) approach was implemented in the analysis covered in Part I. Although DFBI modeling accurately and simultaneously solved the dynamic motion of the valve’s disk along with the flow field of the surrounding fluid, it shed little light on the reason behind such motion. For Part II, the reacting torque of the fluid on the disk was decomposed into representations of the dynamic coefficients in terms of stiffness, damping, and added mass. These were evaluated via simulations with steady-state static (stiffness), constant angular speed (damping), and variable angular speed (added mass) disks. Substituting the dynamic coefficients into Newton’s second law enabled the response trajectories to be obtained. Stable (by average) and unstable equilibrium positions and thrust tendencies of the valve were determined based on the stiffness coefficient (or static torque), the response amplitude was dampened or enlarged by the damping coefficient (minorly affected by added mass), and the response frequency was altered by the damping and added mass coefficients. Although the dynamic coefficient approach renders slightly different trajectories, due to the averaging effect of the torque in comparison to the DFBI method, the overall trend of the response aligns with the DFBI simulation, thus confirming the conclusion in Part I that a fix to the current butterfly valve is necessary.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Foot-scale Evaluation of CO2-Responsive Polymer and CO2 Binary Fluid as an Alternative Fracturing Fluid for Enhanced Geothermal Systems

A foot-scale high-temperature true-triaxial fracturing apparatus was used to conduct hydraulic fracture tests on rock samples. The fracturing tests used different fluids including water, CO2, CO2 with water, and CO2 with a 1% polyallylamine (PAA)/water solution. The tests were conducted with both constant pressure injection and constant flow rate injection modes. CO2-based fracturing fluids were found to produce higher breakdown pressures, high transient flow rates, and produce higher-conductivity fractures as compared to water-based fracturing fluids. Additionally, faster pressurization rates with CO2-based fracturing fluids (obtained when fracturing at constant flow rate mode) are found to be associated with higher fracture conductivities. When fracturing with CO2-based fluids in the presence of a PAA aqueous solution, for example if the rock is saturated with 1wt% PAA in water, the volume expansion caused by CO2-induced crosslinking of PAA leads to a faster pressure increase due to the associated volume expansion and increase in viscosity. It was also found that CO2 as a fracturing fluid injected in hot dry rock (HDR) attain the highest fracture conductivity only when injected at very high flow rates, followed very closely by the CO2/PAA fracturing fluid system that generates fractures with, on average, similarly high conductivity values though independently of injection flow rate and using 1/6 of the mass of CO2 as compared to CO2 in HDR. Breakdown pressures were also similar for CO2 stimulation in HDR and CO2/PAA fluid system under identical injection flow rates. Finally, the well-known low viscosity of CO2 phase prevents the efficient transport of proppants while the reacted CO2/PAA fluid system is known to form a high viscosity binary fracturing fluid system with the potential to carry proppant in addition of being a fine fracturing fluid.

Enhanced Geothermal Systems, polyallylamine, Co2-r↗
DOE OSTI2023

Solid Phase Processing at the Extreme

The effective application of high-flow-stress materials in the realm of Solid Phase Processing frequently necessitates operating at elevated temperatures and substantially high force. These rigorous processing conditions give rise to intricate predicaments concerning the durability of the tools and dies involved, as well as the efficient utilization of process energy. This project is centered on the investigation of the fundamental scientific aspects inherent to these challenges. We aim to delve into the details of how materials react when subjected to intense shear forces at die interfaces, determine the metallurgical interactions occurring at the interfaces with the dies, and investigate the consequences of extreme temperatures on material flow, microstructural evolution, and resultant properties.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Emergent Behavior at the Calcite–Water Interface during Reactive Transport in a Simple Microfluidic Channel

Geochemical reactive transport processes in natural mineral-fluid systems may produce a wide array of emergent phenomena that are difficult to predict from basic principles and to reproduce in model systems. In this study, we present experimental results obtained from a simple microfluidic system with which we explored the consequences of reacting the calcite (104) cleavage surface with an acidic Pb-bearing solution (pH = 3.5, [Pb] total = 5 mM) as a function of flow rate. This system is relevant to passive remediation systems for Pb-rich acid mine drainage. We observed periodic banding in the amounts of Pb sorption at flow velocities ≥ 926 μm s -1 , where the band spacing was spatially correlated with the amount of calcite dissolution and the development of micropyramidal topography on the calcite (104) surface. The equivalent coverage of Pb deposited in these Pb-rich bands was at least several monolayers per unit cell, yet there was no evidence for precipitation of any secondary Pb phase implying incorporation of Pb within the near-surface calcite lattice. We also observed spatial variations in nucleation and growth of euhedral secondary Pb-carbonate minerals hydrocerusite and cerussite at flow rates ≤ 278 μm s -1 ,. These findings demonstrate potential for exploiting the rich phenomenology afforded by the interplay among transport phenomena and chemical kinetics in experimental systems designed to yield deeper insights into geochemical self-organization.

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

Ammonia Permeation Through Inorganic Membranes

Moisture probes are used in the tritium facility to measure the moisture content of the gas stream. When tritium reacts with nitrogen to form ammonia, the ammonia gas damages the sensitive moisture probe. The gas stream tested during this experiment for moisture probe protection is composed of nitrogen and argon gas in conjunction with ammonia and water vapor. The gas stream is filtered through an inorganic membrane due to the high reactivity of tritium and a residual gas analyzer (RGA) is used to measure the ammonia and water vapor in the gas stream, as well as the permeability of nitrogen and argon through the inorganic membranes. The membranes tested were boron nitride and silver. The permeability was compared for each membrane at different pressures and flow rates. Fourier-transformed infrared spectroscopy (FT-IR) was used to compare the membranes before and after ammonia permeation. Ammonia gas is a byproduct formed from tritium gas reacting with nitrogen in the gas stream, similar to the Haber Process which occurs when hydrogen reacts with nitrogen [2]. The corrosive nature of ammonia damages metals and sensitive instruments such as moisture probes. Moisture probes are essential in various gas and liquid processes for measuring the moisture content of a gas. Membrane research has found that graphene oxide membranes successfully uptake ammonia from the gas stream at a humidity of approximately 30% at which point the pores pin closed, preventing gas permeation. Further research is being conducted to find an inorganic membrane that will prevent ammonia permeation thus circumventing damage to the sensitive moisture probe. The permeation measurements for argon and nitrogen gas and water and ammonia vapor were measured at various pressures of 800, 1000, 1200, 1500, 1800, and 2000 torr. The permeability average (PA) was calculated for each filter using a Barrer calculation. The average permeations for each filter were then compared. The gas composition was measured using the RGA Nitrogen was permeated through the membranes after ammonia permeation and a FT-IR spectrum was obtained. Boron Nitride: As the nitrogen pressure increased, the permeability of the membrane decreased from 800-1200 torr and remained relatively unchanged from 1200-2000 torr. The membrane exhibited a slight decrease in argon permeability from 800-1200 torr and a slight increase from 1500-2000 torr. The boron nitride membrane when permeated with ammonia and nitrogen prevented less ammonia permeation compared to the blank filter. The boron nitride filter demonstrated a decrease in permeability after ammonia permeation. Silver Membrane: When permeated with nitrogen alone, the membrane permeability decreased from 800 torr to 1200 torr. The permeability of the membrane showed little difference as pressure increased from 1200 torr. When permeated with argon, the same phenomenon occurred but on a much smaller scale. When compared to the blank and the boron nitride membranes the silver membrane exhibited significantly lower permeability. This may be in part due to the coating on the silver membrane and differences in pore size of the membranes. During permeation with ammonium hydroxide vapor and nitrogen, the partial pressure of the ammonia increased as compared to the blank membrane. The permeability before and after permeation with ammonia and nitrogen were compared. Permeability values suggested that after being permeated with the ammonia that the permeability for lower pressures increased slightly. There was no significant change in the FT-IR spectrums before and after permeation. Future work: Continue testing various inorganic membranes for ammonia permeation. Further research chemical separation of ammonia from the gas stream, which may prove difficult due to the high reactivity of tritium.

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

Implementation of a High-Fidelity Interface Resolving Method in Nek5000

The development and utilization of computational fluid dynamics (CFD) models for large, high-temperature electric melters in the Waste Treatment and Immobilization Plant (WTP) in eastern Washington State have proven to be valuable for various purposes. These models allow for a better understanding of the physio-chemical processes occurring within the melter vessels and can contribute to improving operational efficiency, throughput, and addressing operational issues related to vitrification. The CFD models employed for these melter vessels incorporate multiphase fluid flow and heat transfer simulations in different regions, including the plenum, cold cap, and molten glass regions. As the tank waste and glass formers are introduced into the melter, a reacting batch layer known as the cold cap forms on top of the molten glass. To enhance the melt rate, forced convection bubblers located at the bottom of the melters generate convection currents that help homogenize the molten glass and provide heat to the cold cap. As the bubbles rise through the highly viscous glass, they adopt a spherical-cap shape [1]. Meanwhile, the conversion of the batch to glass generates significant amounts of gases (such as water vapor, carbon dioxide, sulfur dioxide, and NOx) due to thermal decomposition [2]. These gases become trapped between the cold cap and molten glass, forming a foam layer [3]. For modeling multiphase flow in CFD and heat transfer simulations of waste glass melters at different scales, efforts are underway to augment the capabilities of the Nek5000 [4] and NekRS [5] open-source codes [6]. Nek5000/NekRS is a scalable and efficient spectral element code that has been successfully applied to a wide range of fluid dynamics problems. By leveraging the Nek5000/NekRS software, it becomes possible to model the melter systems more affordably and with lower computational requirements compared to currently utilized commercial CFD software. The specific objectives of this ongoing effort include: 1. Implementation of a level set method in Nek5000/NekRS: The level set method is a numerical technique commonly employed in CFD simulations to track and represent the interface between different phases or materials accurately. By incorporating this method into Nek5000/NekRS, the ability to simulate multiphase flows in waste glass melters at a high level of fidelity can be achieved. 2. Demonstration of capability for air bubbling through molten glass: As part of the development process, a specific case of air bubbling through molten glass will be simulated using the augmented Nek5000/NekRS code. This demonstration aims to showcase the ability of the software to accurately capture and analyze the complex phenomena involved in the multiphase flow within waste glass melters. By achieving these objectives, the improved Nek5000/NekRS code will offer a powerful computational tool for simulating and analyzing waste glass melter systems, enabling better understanding, optimization, and troubleshooting of these vitrification processes. The ability to accurately model and simulate multiphase flows has broad relevance across many industries and scientific domains, and the improved functionality can contribute to advancements in various fields beyond waste glass melter simulations.

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

Limitations of cetane number to predict transient combustion phenomena in high-pressure fuel sprays

Fundamental understanding of in-cylinder processes in diesel engines is important to screen emerging biofuels and advanced combustion modes that can reduce greenhouse gas emissions and regulated pollutants including soot. In this study, the role of fuel properties on spray development and combustion is investigated by systematically isolating chemical and thermophysical effects. Three different fuels are considered, two with similar chemical properties and two with similar thermophysical properties with one fuel common to both groups. Experiments are performed in a constant-pressure flow chamber de-signed to provide stable test conditions and facilitate acquisition of at least 150 injections in quick succession for each fuel under reacting conditions at high-pressure, high-temperature ambient conditions using a modified conventional diesel engine injector. Further, high speed optical diagnostics including rainbow schlieren deflectometry, OH* chemiluminescence, and a two-color pyrometry system are employed to simultaneously image the transient spray and reacting jets. Image analysis is performed to determine liquid length, vapor penetration length, timing and location of first stage and main ignition events, lift-off location, total soot mass, and more. Results show that fuels with similar chemical properties or cetane number (CN) exhibit similar delay times for first stage and main ignition events as may be expected, but very different liquid length, first stage and main ignition locations, lift-off length, apparent turbulent flame speed, and soot formation. As such, the ability to characterize candidate biofuels with CN or other parameters derived from simple flame configurations is called into question. In this study, thermo-physical properties controlling the liquid length are identified as the main contributing factor for the observed differences.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Superior performance biodiesel from biomass-derived fusel alcohols and low grade oils: Fatty acid fusel esters (FAFE)

We demonstrated production of a superior performance biodiesel referred to here as fatty acid fusel alcohol esters (FAFE) – by reacting fusel alcohols (isobutanol, 3-methyl-1-butanol, and (S)-(-)-2-methyl-1-butanol) with oil (glyceryl trioleate) using lipase from Aspergillus oryzae. Reaction conditions corresponding to a molar ratio of 5:1 (fusel alcohols to oil), enzyme loading of 2% w/w, reaction temperature of 35 °C, shaking speed of 250 rpm, and reaction time of 24 h achieved >97% conversion to FAFE. Further, FAFE obtained from reacting a fusel alcohol mixture with corn oil were evaluated for use as a fuel for diesel engines. FAFE mixtures showed superior combustion and cold-flow properties, with the derived cetane numbers up to 4.8 points higher, cloud points up to -6 °C lower, and the heat of combustion up to 2.1% higher than the corresponding FAME samples, depending on the fusel mixture used. This represents a significant improvement for all three metrics, which are typically anti-correlated. Finally, FAFE provides a new opportunity for expanded usage of biodiesel by addressing feedstock limitations, fuel performance, and low temperature tolerance.

09 BIOMASS FUELS↗
DOE OSTI2025

Direct numerical simulations for hybrid rocket boundary layers: Performance modeling and scaling

This paper presents a comprehensive performance and scaling analysis of direct numerical simulations for reacting boundary layers, focusing on slab burner configurations. Using a PETSc-based finite volume CFD framework, the study evaluates the scalability and computational cost of flow, chemistry, and radiation evaluations across 2D and 3D simulations. Polymethyl methacrylate (PMMA) is the fuel with pure O 2 as the oxidizer, modeled using a detailed chemical kinetics mechanism with 113 species and 660 reactions. A ray-tracing-based radiation solver, designed for distributed memory applications, is implemented to model radiation heat transfer. Parallel scalability is analyzed for the coupled flow, chemistry, and radiation heat transfer processes. Weak and strong scaling studies are conducted on up to 15,000 computational ranks, revealing robust performance when flow cells exceed 200 per rank. Chemistry evaluations dominate the computational cost in large 3D simulations, accounting for approximately 40% of the total runtime, while flow processes contribute around 35%, and radiation solver contributions remain below 10% due to reduced evaluation frequencies. GPU accelerated chemistry evaluation, implemented with Zero-RK, demonstrates significant promise, achieving up to a 4x speedup for workloads exceeding 30,000 cells per GPU. However, diminishing returns are observed for smaller workloads due to CPU-GPU communication overhead. This study identifies key challenges, including memory bottlenecks and the effects of domain partitioning on flow scalability, while highlighting the potential of GPU-accelerated chemistry to reduce computational costs. In conclusion, these findings provide realizable run configurations for 2D, 3D, and GPU-accelerated cases, offering insights for optimizing reactive flow solvers.

CFD Scalability↗
DOE OSTI2022

Comprehensive uncertainty quantification (UQ) for full engineering models by solving probability density function (PDF) equation

This report details a new method for propagating parameter uncertainty (forward uncertainty quantification) in partial differential equations (PDE) based computational mechanics applications. The method provides full-field quantities of interest by solving for the joint probability density function (PDF) equations which are implied by the PDEs with uncertain parameters. Full-field uncertainty quantification enables the design of complex systems where quantities of interest, such as failure points, are not known apriori. The method, motivated by the well-known probability density function (PDF) propagation method of turbulence modeling, uses an ensemble of solutions to provide the joint PDF of desired quantities at every point in the domain. A small subset of the ensemble is computed exactly, and the remainder of the samples are computed with approximation of the driving (dynamics) term of the PDEs based on those exact solutions. Although the proposed method has commonalities with traditional interpolatory stochastic collocation methods applied directly to quantities of interest, it is distinct and exploits the parameter dependence and smoothness of the dynamics term of the governing PDEs. The efficacy of the method is demonstrated by applying it to two target problems: solid mechanics explicit dynamics with uncertain material model parameters, and reacting hypersonic fluid mechanics with uncertain chemical kinetic rate parameters. A minimally invasive implementation of the method for representative codes SPARC (reacting hypersonics) and NimbleSM (finite- element solid mechanics) and associated software details are described. For solid mechanics demonstration problems the method shows order of magnitudes improvement in accuracy over traditional stochastic collocation. For the reacting hypersonics problem, the method is implemented as a streamline integration and results show very good accuracy for the approximate sample solutions of re-entry flow past the Apollo capsule geometry at Mach 30.

42 ENGINEERING↗
DOE OSTI2021

Membrane-Induced Vanadium Crossover-Blocking Polybenzimidazole Copolymer with Exceptional Proton Selectivity

Broad deployment of redox flow batteries is hindered, in part, due to the lack of highly selective ion-exchange membranes with high proton conductivity. We report a post-functionalization strategy for polybenzimidazole copolymer (PBI-co) membranes by a reacting mixture of concentrated sulfuric and phosphoric acids that led to improved membrane selectivity while retaining high proton conductivity. Fourier-transform infrared spectroscopy confirmed successful sulfonation and protonation of benzimidazole functionalities. PBI-co exhibits a very high room-temperature proton conductivity of 0.138 S cm -1 , as measured by electrochemical impedance spectroscopy. PBI-co also demonstrates electromigration crossover mitigation under applied current densities up to 500 mA cm -2 , and voltage pulses up to 1.4 V, while also inhibiting diffusion-driven vanadium crossover for over 20 days. Single-cell vanadium flow battery testing on post-functionalized PBI-co further confirms state-of-the-art battery performance with an ~99.5% Coulombic efficiency at 100 mA cm -2 and an area-specific resistance of ~20 mO cm 2 lower than Nafion-212 membrane.

25 ENERGY STORAGE↗
DOE OSTI2023

Refinery gas processing method

A process for hydrogen recovery from refinery gas system comprising supplying the refinery gas to an inlet manifold fluidly coupled to a conditioning stage, the conditioning stage comprising a reactor having a reforming catalyst deposited on an ultra-short-channel-length metal substrate; supplying oxidant to the conditioning stage via the inlet manifold; supplying steam from a steam generator to the conditioning stage via the inlet , manifold; reacting the refinery gas in the conditioning stage; and discharging a product through a discharge outlet fluidly coupled to the conditioning stage, the discharge outlet configured to flow the product for use by a downstream reformer. The process allows to either increase the H2 production rate or lower the firing rate while maintaining a constant H2 production rate for the downstream steam reformer, independent of the feed compositional variability of the refinery or still gas.

Vilekar, Saurabh↗
DOE OSTI2022

High-Fidelity Simulations of Supercritical CO2 Based Oxy-Combustion with Non-Ideal Equation of State

The supercritical carbon dioxide (sCO2) Allam-Fetvedt cycle employs high-temperature oxy-combustion using sCO2 as the working fluid, leading to a unique and challenging environment for materials. In this cycle, pressurized gaseous fuel is combusted in the presence of a hot oxidant flow containing a mixture of CO2 and oxygen and a hot CO2 diluent recycle stream at approximately 300 bar under lean combustion conditions. We present high-fidelity reacting simulations of this oxy-combustor environment and compare the results when performing simulations with ideal and non-ideal equations of state. By doing so, we identify if the combustor behavior at the regime of interest is sensitive enough to the equation of state to justify the added computational cost incurred with a non-ideal equation of state.

equation of state↗
DOE OSTI2021

Core-flood Effluent and Shale Surface Chemistries in Predicting Interaction between Shale, Brine, and Reactive Fluid

We report field and laboratory observations to date indicate that the efficiency of hydraulic fracturing, as it relates to hydrocarbon recovery, depends significantly on geochemical alterations to rock surfaces that diminish accessibility by partial or total plugging of the pore and fracture networks. This is caused by mineral scale deposition such as coating of fracture surfaces with precipitates, particle migration, and/or crack closure due to dissolution under stress. In reactive flow-through experiments, mineral reactions in response to acidic fluid injection significantly reduced system porosity and core permeability. The present study focuses on changes to fluid chemistry and shale surfaces (inlet and fracture walls) resulting from shale-fluid interactions and integrating these findings for an improved estimate of transport-related consequences. The reacted shale surfaces were examined by spatially-resolved scanning electron microscopy - energy dispersive spectroscopy (SEM-EDS) analysis. Importantly, inductively coupled plasma - mass spectrometry/optical emission spectroscopy (ICP-MS/OES) was utilized to probe the chemical evolution of the core-flood effluents. The three study cores selected from the Marcellus formation represent different mineralogies and structural features. In flow-through experiments, lab-generated brine and HCl-based fracture fluid (pH=2) were injected sequentially under effective stress (up to 500 psi) at reservoir temperature (80°C). SEM-EDS results confirmed by the ICP concentration trends showed significant Fe hydroxide precipitates in clay- and pyrite-rich outcrop samples due to partial oxidation of Fe-bearing phases in the case of intrusion of low salinity water-based fluids. Porosity reduction in the MSEEL (Marcellus Shale Energy and Environmental Laboratory) carbonate-rich sample is related to compaction of cores under stress due to matrix softening with substantial dissolution, and pore-filling by hydroxides, as well as barite and salts. Despite the same fluid compositions and experimental conditions used for both MSEEL samples, barite precipitation was much more intense in the MSEEL clay-rich sample due to its greater sorption capacity and additional sulfate source as well as fissile nature with multiple lengthwise cracks. ICP tests revealed time-resolved concentration trends in produced brine and reactive fluids that in turn complemented the pre-/post-reaction SEM-EDS observations.

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

Reservoir-scale model of geologic hydrogen production from serpentinization: Cyclic injection in a dual-permeability fracture-matrix system

Geologic hydrogen (GeoH 2 ) from serpentinization is a promising low-carbon resource, but its reservoir-scale behavior remains poorly understood. We develop a dual-permeability reactive transport model for an injector–producer well pair in ultramafic rock that couples multiphase flow, heat transfer, geochemistry, and porosity–permeability evolution. Here, the model is calibrated to olivine flow-through experiments and upscaled to two-year long simulations with continuous injection and cyclic injection with shut-in-to-injection ratios (SIR = 1, 0.5, 0.1). Olivine reacts along high-flux pathways to form lizardite and magnetite, increasing pH; H 2 (aq) and H 2 (g) peak early and then decline as exsolution and advective export outpace local generation. Continuous injection yields the highest cumulative H 2 but the lowest water-use efficiency (1.138 x 10 –6 mol/kgw). Cyclic injection increases this ratio to 1.35 x 10 –6 , 1.377 x 10 –6 , and 1.182 x 10 –6 mol/kgw for SIR = 1, 0.5, and 0.1, respectively; SIR = 0.5 provides a ~20% improvement over continuous injection and the best compromise for pilot design.

08 HYDROGEN↗