Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reaction rates”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Synthesis of wafer-scale SWCNT forests with remarkably invariant structural properties in a bulk-diffusion-controlled kinetic regime

Robust synthesis of vertically aligned carbon nanotubes (VACNT) at large scale is required to accelerate deployment of numerous cutting-edge devices to emerging commercial applications. To address this need, this work demonstrates that the structural characteristics of single-walled CNTs (SWCNTs) produced in a growth regime dominated by bulk diffusion of the gaseous carbon precursor are remarkably invariant over a broad range of process conditions (precursor concentration increased up to 30-fold, catalyst substrate area from 1 cm 2 to 180 cm 2 , growth pressure from 20 to 790 mbar, and gas flowrates up to 8-fold). A simple growth kinetics model is derived that accounts for both reactant bulk diffusion and competing byproduct formation in the presence of excess hydrogen, which quantitatively reproduces the experimental data in the entire isothermal parameter space. The model enables critical predictions for process scale-up optimization including a potential 6-fold increase in CNT production rate, a ~90% carbon conversion efficiency in select conditions, and elimination of reaction rate decay observed at high pressures by using a hydrogen-free growth environment. This model-guided opportunity to substantially improve synthesis throughput and efficiency, in conjunction with the remarkably invariant CNT properties, is appealing for reliable VACNT device fabrication at both small and industrial scales.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Towards ternary binders involving limestone additions — A review

The review summarises literature to examine the transition from portland limestone cement system to composite ternary binder systems involving limestone. Interest in limestone addition as an ideal component in multicomponent binder systems has surged as evident from the large volume of literature published in the recent past. A ternary blended system, with co-substitution of limestone, has the potential to complement the reaction of the supplementary cementitious materials (SCMs). The direct addition of limestone powder helps to attain higher substitution levels of portland cement clinker, improve early age properties, and supplement SCM's reactivity. However, the dilution of hydrates could hamper the long-term benefits. In this review, the interaction of fine limestone is classified and elaborated under two broad umbrellas: physical and chemical interactions. The physical interactions can manifest in three ways, namely, filler action, shearing action and improved packing, which alters reaction rate and extent at early ages. The chemical interactions also modify the reaction kinetics and phase assemblage due to nucleation of C-S-H on calcite surfaces, preservation of the ettringite phase and formation of carboaluminates. Two different forms of carboaluminate hydrates — hemicarboaluminate and mono-carboaluminate can be present in the hydrate matrix depending on the balance between carbonate ions and aluminates in the pore solution. Several factors such as replacement level, particle size, choice of SCM, its reactivity and reactive aluminates content, sulphate levels, curing temperature, and duration of curing can control the carboaluminate formation, reaction degree of SCMs and chemical interaction from limestone additions. A combination of physical and chemical effects makes fine limestone a potential material for co-substitution with aluminosilicate based SCMs, mainly fly ash, slag, and calcined clay. In this review, the factors affecting limestone-SCM composites are summarised based on a detailed literature survey. The effects of SCM-limestone cement composites on hydration kinetics, reaction chemistry, the reactivity of SCMs, the stability of hydrated phases, and contribution to the physical structure development and macroscopic properties by evaluating hydration and mechanical properties are discussed. The importance of AFm (Al{sub 2}O{sub 3}–Fe{sub 2}O{sub 3}-mono) phases in various deterioration mechanisms in concrete and their influence on performance characteristics in different exposure environment is critically appraised.

36 MATERIALS SCIENCE↗
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

Entropic Control of HD Exchange Rates over Dilute Pd-in-Au Alloy Nanoparticle Catalysts

Dilute Pd-in-Au alloy catalysts are promising materials for selective hydrogenation catalysis. Extensive surface science studies have contributed mechanistic insight on the energetic aspect of hydrogen dissociation, migration and recombination on dilute alloy systems. Yet, translating these fundamental concepts to the kinetics and free energy of hydrogen dissociation on nanoparticle catalysts operating at ambient pressures and temperatures remains challenging. Here, the effect of the Pd concentration and Pd ensemble size on the catalytic activity, apparent activation energy and rate limiting process is addressed by combining experiment and theory. Experiments in a flow reactor show that a compositional change from 4 to 8 atm% Pd of the Pd-in-Au alloy catalyst leads to strong increase in activity, even exceeding the activity per Pd atom of monometallic Pd under the same conditions, albeit with an increase in apparent activation energy. First-principles calculations show that the rate and apparent activation enthalpy for HD exchange increase when increasing the Pd ensemble size from single Pd atoms to Pd trimers in a Au surface, suggesting that the ensemble size distribution shifts from mainly single Pd atoms on the 4 atm% Pd alloy to larger Pd ensembles of at least three atoms for the 8 atm% Pd/Au catalyst. The DFT studies also indicated that the rate-controlling process is different: H 2 (D 2 ) dissociation determines the rate for single atoms whereas recombination of adsorbed H and D determines the rate on Pd trimers. Both experiment and theory suggest that the increased reaction rate with increasing Pd content and ensemble size stems from an entropic driving force. Lastly, our results support hydrogen migration between Pd sites via Au and indicate that the dilute alloy design prevents the formation of subsurface hydrogen, which is crucial in achieving high selectivity in hydrogenation catalysis.

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

First principles reaction discovery: from the Schrodinger equation to experimental prediction for methane pyrolysis

Our recent success in exploiting graphical processing units (GPUs) to accelerate quantum chemistry computations led to the development of the ab initio nanoreactor, a computational framework for automatic reaction discovery and kinetic model construction. In this work, we apply the ab initio nanoreactor to methane pyrolysis, from automatic reaction discovery to path refinement and kinetic modeling. Elementary reactions occurring during methane pyrolysis are revealed using GPU-accelerated ab initio molecular dynamics simulations. Subsequently, these reaction paths are refined at a higher level of theory with optimized reactant, product, and transition state geometries. Reaction rate coefficients are calculated by transition state theory based on the optimized reaction paths. The discovered reactions lead to a kinetic model with 53 species and 134 reactions, which is validated against experimental data and simulations using literature kinetic models. We highlight the advantage of leveraging local brute force and Monte Carlo sensitivity analysis approaches for efficient identification of important reactions. Both sensitivity approaches can further improve the accuracy of the methane pyrolysis kinetic model. The results in this work demonstrate the power of the ab initio nanoreactor framework for computationally affordable systematic reaction discovery and accurate kinetic modeling.

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

Steam-created grain boundaries for methane C–H activation in palladium catalysts

Defects may display high reactivity because the specific arrangement of atoms differs from crystalline surfaces. We demonstrate that high-temperature steam pretreatment of palladium catalysts provides a 12-fold increase in the mass-specific reaction rate for carbon-hydrogen (C–H) activation in methane oxidation compared with conventional pretreatments. Through a combination of experimental and theoretical methods, we demonstrate that an increase in the grain boundary density through crystal twinning is achieved during the steam pretreatment and oxidation and is responsible for the increased reactivity. The grain boundaries are highly stable during reaction and show specific rates at least two orders of magnitude higher than other sites on the palladium on alumina (Pd/Al 2 O 3 ) catalysts. Theoretical calculations show that strain introduced by the defective structure can enhance C–H bond activation. Introduction of grain boundaries through laser ablation led to further rate increases.

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

Proton branching ratios of 23 Mg levels

Background: The anomalous 22 Ne abundance measured in certain presolar graphite grains is thought to arise from the decay of 22 Na that was synthesized at high temperatures in core-collapse supernovae. To better interpret this abundance anomaly, the primary destruction mechanism of 22 Na, the 22 Na(p,γ) 23 Mg reaction, must be better understood. Purpose: Determine proton branching ratios of several 23 Mg excited states that play a role in the high-temperature 22 Na(p,γ) 23 Mg reaction rate. Methods: Particle decays of 23 Mg excited states populated with the previously reported 24 Mg(p,d) 23 Mg transfer reaction measurement [Kwag et al., Eur. Phys. J. A 56, 108 (2020)] were analyzed to extract proton branching ratios. The reaction was studied using a 31-MeV proton beam from the Holifield Radioactive Ion Beam Facility of Oak Ridge National Laboratory and 24 Mg solid targets. Results: Proton branching ratios of several 23 Mg excited states in the energy range of Ex = 8.044 - 9.642 MeV were experimentally determined for the first time for the p0 and p1'(p1+p2+p3) decay channels. Conclusions: These new branching ratios for 23 Mg levels can provide an experimental foundation for an improved high-temperature rate of the 22 Na(p,γ) 23 Mg reaction needed to understand production of anomalously high 22 Ne abundance in core-collapse supernovae.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Multiple redox mechanisms for water-gas shift reaction on Fe 3 O 4 (1 1 1) surface: A density functional theory and mean-field microkinetic modeling study

In this work, we investigated redox mechanisms for the high temperature (HT) water–gas shift reaction (WGSR) on the Fe 3 O 4 (1 1 1) surface using density functional theory (DFT) calculations and mean-field microkinetic modeling (MKM). The redox pathways branched into three Langmuir-Hinshelwood processes (LH1, LH2 and LH3) and one Mars-van-Krevelen (MVK) process (in the presence of oxygen vacancy) from multidentate binding of CO and CO 2 over four available reactive sites (Fe oct2 , Fe tet1 , Fe bridge, and O1). We found that the LH1 and LH2 processes have CO adsorption at a single iron site (Fe oct2 or Fe tet1 ), while the LH3 and MVK processes have stronger chemisorption of CO or CO 2 by both the Fe oct2 site and the O 1 site. From the mechanistic study of these reaction path, we recognized that availability of O 1 sites was key to proceed to either CO oxidation by the single site (LH1, LH2) or by the dual site (LH3, MVK). We observed that the reaction energetics were significantly different in the CO oxidation steps, where the single or dual site results in markedly different apparent activation energy (73–281 kJ/mol) and reaction rates (10 -3 –10 -10 mol∙m -2 ∙s -1 ) among the four reaction mechanisms. The utilization of mean field MKM with DFT reaction energetics helps explain the experimental debates for the catalytic reaction details as well as providing a possible direction to engineer the catalyst with higher activity.

42 ENGINEERING↗
DOE OSTI2022

Operability thresholds for thermally damaged EBW detonators

Operability thresholds that differentiate between functional RP-87 exploding bridge wire (EBW) detonators and nonfunctional RP-87 EBW detonators (duds) were determined by measuring the time delay between initiation and early wall movement (function time). The detonators were inserted into an externally heated hollow cylinder of aluminum and fired with current flow from a charged capacitor using an exploding bridge wire (EBW initiated). Functioning detonators responded like unheated pristine detonators when the function time was 4 μs or less. The operability thresholds of the detonators were characterized with a simple decomposition cookoff model calibrated using a modified version of the Sandia Instrumented Thermal Ignition (SITI) experiment. These thresholds are based on the calculated state of the PETN when the detonators fire. The operability threshold is proportional to the positive temperature difference (ΔT) between the maximum temperature within the PETN and the onset of decomposition (~406 K). Furthermore, the temperature difference alone was not sufficient to define the operability threshold. The operability threshold was also proportional to the time that the PETN had been at elevated temperatures. That is, failure was proportional to both temperature and reaction rate. The reacted gas fraction is used in the current work for the reaction correlation. Melting of PETN also had a significant effect on the operability threshold. Detonator failure occurred when the maximum temperature exceeded the nominal melting point of PETN (414 K) for 45±5 s or more.

42 ENGINEERING↗
DOE OSTI2025

Constraining (n,γ) cross sections via surrogate measurements with hyperion

Indirect measurement techniques, such as the surrogate method, are necessary to constrain cross sections and reaction rates that are inaccessible for direct measurement. The surrogate method uses an alternate reaction channel to populate a nucleus of interest, and combines resulting experimental data with theory to constrain the desired (n,γ) cross section. Individual γ-decay probabilities are experimentally extracted using γ-ray spectroscopy. This project utilizes surrogate analysis of 64,70 Zn(p,d) and 64,70 Zn(p,t) reactions measured with the Hyperion array at Texas A&M University's Cyclotron Facility to probe the relevant neutron capture reactions 61,62,67,68 Zn(n,γ). As an example, preliminary γ-decay probabilities are presented for 69 Zn.

Neutron capture↗
DOE OSTI2020

Nuclear spin conservation enables state-to-state control of ultracold molecular reactions

Quantum-state control of reactive systems has enabled microscopic probes of underlying interaction potentials and the alteration of reaction rates using quantum statistics. However, extending such control to the quantum states of reaction outcomes remains challenging. Here in this work, we realize this goal by utilizing the conservation of nuclear spins throughout the reaction. Using resonance-enhanced multiphoton ionization spectroscopy to investigate the products formed in bimolecular reactions between ultracold KRb molecules we find that the system retains a near-perfect memory of the reactants’ nuclear spins, manifested as a strong parity preference for the rotational states of the products. We leverage this effect to alter the occupation of these product states by changing the coherent superposition of initial nuclear spin states with an external magnetic field. In this way, we are able to control both the inputs and outputs of a reaction with quantum-state resolution. The techniques demonstrated here open up the possibilities to study quantum entanglement between reaction products, and ultracold reaction dynamics at the state-to-state level.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Lasers for the observation of multiple order nuclear reactions

Nuclear reaction rates become nonlinear with respect to flux (cm −2 s −1 ) in extreme environments such as those found during stellar nucleosynthesis and terrestrial nuclear detonations. To observe these effects directly in the laboratory, extremely high particle fluences (cm −2 ) are necessary but not sufficient. Reactor-based neutron sources, such as the Institut Laue-Langevin’s high-flux neutron reactor, were previously the closest to meeting this challenge, albeit over ∼hour time scales. In ultra-high flux environments, where multiple reactions occur on picosecond time scales, nuclei are unable to return to their ground states between reactions; consequently, reactions take place on excited nuclei. To accurately model high-flux environments, data on the cross-sections of excited nuclear states are required, which differ significantly from those of ground states due to spin/parity effects. In order to replicate these effects in the laboratory, short high-fluence pulses on the order of the lifetime of a typical nuclear excited state (generally ≲1 ns) are required. Particle beams generated by high-intensity lasers are uniquely positioned to meet this need with the potential to produce fluences of 10 17 protons/cm 2 and 10 22 neutrons/cm 2 over a few pico-seconds or less. In addition to providing a quantitative analysis of the rates of multiple rapid reactions in general, the present work examines a number of laser-based experiments that could be conducted in the near future to observe multiple rapid reactions for laboratory-based astrophysics and the measurement of exotic cross-sections.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Effects of oligomerization and decomposition on the nanoparticle growth: a model study

Abstract. The rate at which freshly formed secondary aerosol particles grow is an important factor in determining their climate impacts. The growth rate of atmospheric nanoparticles may be affected by particle-phase oligomerization and decomposition of condensing organic molecules. We used the Model for Oligomerization and Decomposition in Nanoparticle Growth (MODNAG) to investigate the potential atmospheric significance of these effects. This was done by conducting multiple simulations with varying reaction-related parameters (volatilities of the involved compounds and reaction rates) using both artificial and ambient measured gas-phase concentrations of organic vapors to define the condensing vapors. While our study does not aim at providing information on any specific reaction, our results indicate that particle-phase reactions have significant potential to affect the nanoparticle growth. In simulations in which one-third of a volatility basis set bin was allowed to go through particle-phase reactions, the maximum increase in growth rates was 71 % and the decrease 26 % compared to the base case in which no particle-phase reactions were assumed to take place. These results highlight the importance of investigating and increasing our understanding of particle-phase reactions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Understanding the transport and break up of reactive ejecta

This manuscript investigates reactive- versus hydrodynamic-breakup processes of ejecta. For this paper, the reactive metal is cerium (Ce) and the nonreactive metal is tin (Sn), the nonreactive gas is helium (He) and the reactive gas is deuterium (D 2 ) or hydrogen (H 2 ). Experiments were performed in vacuum and the reactive- and nonreactive-gases at various pressures, where we endeavored to match the post-shock gas densities to differentiate between reactive- versus hydrodynamic-breakup processes. Hydrodynamic breakup sensitively links to the Weber number (gas density, liquid fragment diameter, surface tension, and the square of the relative velocity between the fragment and the gas), whereas reactive breakup links to the reactive dynamics which includes two processes. In one case the reactive metal breaks up into smaller fragments as rapidly as the reaction rate, and in the other a crust grows on the liquid fragments as the reactions occur, a diffusion limited process. In the latter case, the particle diameters increase with time as the crust grows. In this process, which is indicated by the data, particles breakup as the CeD 2 loses strength with increasing temperature, leaving an exponentially increasing diameter.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Bridging Scales to Model Reactive Diffusive Transport in Porous Media

Published by ECS. Two novel scale-bridging algorithms to model reaction-diffusion transport in porous media are presented. The algorithms are based on direct numerical simulations and couple the information of a micro-scale model, which accounts for the large field of view provided by micro X-ray computed tomography (X-ray CT), and a nano-scale model, which locally resolves transport in the fine structure extracted from nano X-ray CT. The micro-scale model is discretized in the through-plane direction into a 1D grid, where effective properties and internal boundaries are determined based on the results from the nano-scale model. The validated algorithms are used to examine transport of oxygen in precious group metal-free electrodes considering both zero- and first-order kinetics. Unlike conventional methods, the results show that the effective diffusivity is not a passive property but increases in regions where the reaction-rate coefficient is large. The proposed algorithms account for the multiscale coupling of reaction-diffusion transport and material microstructure, thus improving the predictions compared to conventional methods.

42 ENGINEERING↗
DOE OSTI2020

Improved Performance of Cu(InGa)(SeS) 2 PV Modules Using the Reaction of Metal Precursors. Final Report

This project “Improved Performance of Cu(InGa)(SeS) 2 PV Modules using the Reaction of Metal Precursors” was a partnership led by the Institute of Energy Conversion (IEC) at the University of Delaware with Columbia University and the Molecular Foundry at the Lawrence Berkeley National Laboratory. The aim was to develop pathways to improve Cu(InGa)(SeS) 2 (CIGSS) thin film photovoltaic modules using processes compatible with low manufacturing cost. The CIGSS approach investigated was a two-step process including deposition of metal precursor films following by reaction in hydride gases utilizing IEC’s novel reactor. The process was similar to that under commercial development by the project’s industry partner Stion. When Stion went out of business mid-project the focus changed to a rapid thermal process considered more commercially viable. Approaches to improve the performance of solar cells using the reacted films focused on two material innovations. First, the overall Ga content was increased to increase the operating voltage, which is desirable for scale-up to commercial modules. Second, the processing and performance advantages arising from Ag alloying were investigated. Advanced characterization guided process and material development including control of relative composition gradients. Research on the formation of Cu-Ga-In metal precursors utilized sputtering deposition which is normally used in commercial applications. The work resulted in processes for deposition of precursor stacks with increased relative Ga content and effects of deposition parameters on morphology and phase composition were established. It was shown that the metal precursor films have comparable phase composition and morphology so subsequent reaction follows from the same starting point. The addition of Ag to the metal precursors gave more uniform morphology and improved adhesion of reacted films which enable higher reaction temperature for faster processing. A novel outcome was the discovery of a previously undocumented material phase in sputter-deposited and evaporated Ag-Cu-In-Ga thin films. Hydride gas reaction processes including time-temperature-concentration profiles were developed for different precursor compositions. This enables control of composition profiles to engineer through-film gradients for solar cell optimization with characterization and simulations used to correlate measured film composition profiles to measurements of devices. In particular, the gradient of sulfur at the front of the CIGSS film was found to be critical. The simulations guided process development leading to improved reproducibility of devices improved performance with higher Ga content and higher voltage. With Ag-alloyed precursors, the reaction pathways leading were determined. A significant finding was that Ag-alloying increases the reaction rate to completely convert precursor films to the final chalcopyrite which could enable reduced reaction time to benefit manufacturability. To maintain potential commercial viability, the process under investigation was refocused to a rapid thermal process that could potentially be incorporated into an in-line process for manufacturing. Precursors with different composition were capped with an extra selenium layer and reacted in hydrogen sulfide 5-15 minutes, compared to typically 2 hours in the previous multi-step batch process. Critical RTP parameters were identified to control the reaction. Further optimization would be needed for high efficiency solar cells but pathways to high quality devices with further optimization and improved heating uniformity were developed. The project also developed new optoelectronic characterization approaches with a focus on development and application of spatial- and time-resolved photoluminescence and a custom mapping photoluminescence microscope built. It was shown how critical electronic transport properties strongly depend on the chemical composition of the material and that a wide range of samples show inhomogeneity on a length scale larger than the grains in the films. Additionally, two-photon excitation capability was developed to distinguish bulk vs surface losses. The project advances the state-of-the -art for precursor reaction processes in several ways that could impact manufacturing. This includes validation of approaches to increase voltage and establishment of model-guided control to form optimal composition profiles. The application of process control approaches with knowledge of phase formation and reaction pathways can be critically valuable in designing a large-scale process.

14 SOLAR ENERGY↗
DOE OSTI2021

Modeling of Thermal Decomposition of TATB-Based Explosive for Safety Analysis

We investigate and model the cook-off behavior of LX-17 to understand the response of explosive systems in abnormal thermal environments. Decomposition has been explored via conventional ODTX (One-Dimensional Time-to-eXplosion), PODTX (ODTX with pressure-measurement), TGA (Thermo-Gravimetric Analysis), and DSC (Differential Scanning Calorimetry) experiments under isothermal and ramped temperature profiles. The data were used to fit reaction rate parameters for proposed schemes in an ALE3D computational model. This model includes chemical reactions, thermo- and hydro-dynamics, and material properties, including thermal expansion, compressibility, and strength. These parameterizations were carried out utilizing a Python evolutionary optimization method on LLNL’s high-performance computing clusters. Additional experiments are being developed to further characterize and monitor decomposition intermediates to improve the model. Once experimentally validated, this model will be scalable to several applications involving LX-17. Furthermore, the optimization methodology developed herein should be applicable to other high explosive materials.

Chemistry - Chemical explosives↗
DOE OSTI2023

Reactive Transport Modeling of Anthropogenic Carbon Mineralization in Stacked Columbia River Basalt Reservoirs

Numerical simulation of CO2 storage in basalts and related reactive lithologies requires modeling complex, coupled hydrologic and chemical processes, including multi-phase flow and transport, partitioning of CO2 into the aqueous phase, and chemical interactions with aqueous fluids and rock minerals. We conducted reactive transport simulations of the Wallula pilot-scale CO2 injection into the flow tops of the Grande Ronde Basalt using our PNNL STOMP-CO2 simulator with the ECKEChem reactive module. Our mineralization simulation of the ~1,000 tons of injected CO2 into the interflow zones was based on the hydrologic transport model we previously developed. For this work, the simulations considered geochemical reactions involving the basalt components, precipitates, formation brine, and injected CO2. In our benchmark case, carbonate minerals precipitated, resulting in ~20% of the CO2 being mineralized in 10 years. Increasing the reaction rate of a single primary mineral phase (clinopyroxene) by an order of magnitude resulted in a carbon mineralization reaction extent of ~90% over the same time interval. Based on these initial sensitivity analysis results, it is clear that a thorough understanding of primary mineral dissolution rates is required for accurately predicting long-term fate and transport of injected CO2 into basalt formations. Our reactive transport numerical simulations will be key components of commercial-scale CO2 storage operation permitting, de-risking, and optimization in mafic and ultramafic reservoirs.

Cao, Ruoshi↗