Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “HCO”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Isotopically labelled inorganic carbon delivered to algalcultures via bubbler bottle

This protocol describes a method for delivering labelled inorganic carbon as 13 CO 2 to algal cultures by bubbling air through a solution of H 13 CO 3 - , and then into the culture. We developed this method to deliver label to cultures grown under continuous bubbling with air, without the use of 13 CO 2 labelled gas and the necessary equipment to mix labelled gas with air at near-atmospheric levels. Bubbling precludes the more common approach of adding H 13 CO 3 - label directly to the media, because dissolved HCO 3 - is in equilibrium with atmospheric CO 2 . Thus, excess HCO 3 - added to a solution will leave the solution as CO 2 gas as it equilibrates. Bubbling rapidly accelerates this equilibration which is typically diffusion limited. The method described here takes advantage of this aspect of carbonate chemistry, and uses a solution of H 13 CO 3 - —which is less expensive and more convenient than 13 CO 2 gas—to generate a flux of 13 CO 2 that can be bubbled into a culture.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

A bicarbonate-rich liquid condensed phase in non-saturated solutions in the absence of divalent cations

Bicarbonate (HCO 3 - ) and sodium (Na + )-containing solutions contain droplets of a separate, bicarbonate-rich liquid condensed phase (LCP) that have higher concentrations of HCO 3 - relative to the bulk solution in which they reside. The existence and composition of the LCP droplets has been investigated by nanoparticle tracking analysis, nuclear magnetic resonance spectroscopy, refractive index measurements and X-ray pair distribution function analysis. The bicarbonate-rich LCP species is a previously unaccounted-for, ionic phenomenon which occurs even in solutions with solely monovalent cations. Its existence requires re-evaluation of models used to describe and model aqueous solution physicochemistry, especially those used to describe and model carbonate mineral formation.

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

LoTSS Jellyfish Galaxies. IV. Enhanced Star Formation on the Leading Half of Cluster Galaxies and Gas Compression in IC3949

With Mapping Nearby Galaxies at APO integral field spectroscopy, we present a resolved analysis of star formation for 29 jellyfish galaxies in nearby clusters, identified from radio continuum imaging taken by the Low Frequency Array. Simulations predict enhanced star formation on the "leading half" (LH) of galaxies undergoing ram pressure stripping, and in this work we report observational evidence for this elevated star formation. The dividing line (through the galaxy center) that maximizes this star formation enhancement is systematically tied to the observed direction of the ram-pressure-stripped tail, suggesting a physical connection between ram pressure and this star formation enhancement. We also present a case study on the distribution of molecular gas in one jellyfish galaxy from our sample, IC3949, using Atacama Large Millimeter/submillimeter Array CO J = 1 – 0, HCN J = 1 – 0, and HCO + J = 1 – 0 observations from the ALMA MaNGA Quenching and Star Formation Survey. The H2 depletion time (as traced by CO) in IC3949 ranges from ~1 Gyr in the outskirts of the molecular gas disk to ~11 Gyr near the galaxy center. IC3949 shows a clear region of enhanced star formation on the LH of the galaxy where the average depletion time is ~2.7 Gyr, in line with the median value for the galaxy on the whole. Dense gas tracers, HCN and HCO + , are only detected at the galaxy center and on the LH of IC3949. Our results favor a scenario in which ram pressure compresses the interstellar medium, promoting the formation of molecular gas that in turn fuels a localized increase of star formation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Superhydrophobic Phosphonium Modified Robust 3D Covalent Organic Framework for Preferential Trapping of Charge Dispersed Oxoanionic Pollutants

Abstract Charge dispersed and less hydrophilic anionic pollutants are often difficult to be preferentially captured by common cationic framework materials. A gamma radiation modification approach is employed to anchor a superhydrophobic phosphonium‐containing coating on a robust crystalline 3D covalent organic framework (COF). By regulating the surface hydrophobicity, the prepared 3DCOF‐ g ‐VBPPh 3 Cl is endowed with a strong affinity for charge dispersed and less hydrophilic oxoanions, such as MnO 4 − , TcO 4 − and ReO 4 − , to surmount the Hofmeister bias, which much favors remediation of oxoanionic pollutants in complex aqueous systems. Batch and column experiments with 3DCOF‐g‐VBPPh 3 Cl in both tap water and simulated groundwater are performed, and rapid sequestration is achieved with removal efficacy up to 99.995%, record‐high distribution coefficient of 1.0 × 10 8 mL g −1 , and desirable priority over competing anions such as Cl − , SO 4 2− , HCO 3 − , and NO 3 − , confirming promise for remediation of charge dispersed anionic pollutants. Theoretical calculations reveal a mechanism of preferential capture based on electrostatic and dispersion interaction between charge dispersed anions and 3DCOF‐g‐VBPPh 3 Cl that effectively competes with solvation of the aqueous anions.

Wang, Yue↗
DOE OSTI2025

Role of Histidine‐Containing Peptoids in Accelerating the Kinetics of Calcite Growth

Carbonate mineralization, the conversion of CO 2 into stable, thermodynamically favorable carbonate minerals, offers a promising strategy for permanent and environmentally friendly carbon storage, with minimal risk of long-term leakage and minimal monitoring requirements. Drawing inspiration from carbonic anhydrase (CA), a family of zinc-containing metalloenzymes that catalyze the hydration of CO 2 to bicarbonate and promote carbonate precipitation, a class of histidine-containing peptoids was designed that is capable of coordinating with Zn 2+ ions to act as CA mimetics for accelerating calcite step growth. In situ atomic force microscopy (AFM) measurements reveal that these peptoids significantly enhance step advancement, with a more pronounced effect observed when combined with Zn 2+ ions and under higher calcium-to-carbonate activity ratios, indicating that peptoids facilitate the incorporation of CO 3 2− ions at step edges. Solution NMR and 3D atomic force microscopy (3D AFM) analyses show that the coordination of peptoids with Zn 2+ promotes both the deprotonation of HCO 3 − to CO 3 2− and restructures the interfacial hydration layers of calcite, collectively lowering the activation barrier for step growth. These findings establish a design framework for sequence-defined polymers to regulate carbonate mineralization, offering promising applications in CO 2 capture and long-term storage.

CO2 mineralization↗
DOE OSTI2022

Carbon Dioxide Capture at Nucleophilic Hydroxide Sites in Oxidation‐Resistant Cyclodextrin‐Based Metal–Organic Frameworks**

Abstract Carbon capture and sequestration (CCS) from industrial point sources and direct air capture are necessary to combat global climate change. A particular challenge faced by amine‐based sorbents—the current leading technology—is poor stability towards O 2 . Here, we demonstrate that CO 2 chemisorption in γ‐cylodextrin‐based metal–organic frameworks (CD‐MOFs) occurs via HCO 3 − formation at nucleophilic OH − sites within the framework pores, rather than via previously proposed pathways. The new framework KHCO 3 CD‐MOF possesses rapid and high‐capacity CO 2 uptake, good thermal, oxidative, and cycling stabilities, and selective CO 2 capture under mixed gas conditions. Because of its low cost and performance under realistic conditions, KHCO 3 CD‐MOF is a promising new platform for CCS. More broadly, our work demonstrates that the encapsulation of reactive OH − sites within a porous framework represents a potentially general strategy for the design of oxidation‐resistant adsorbents for CO 2 capture.

Zick, Mary E.↗
DOE OSTI2022

Carbon Dioxide Capture at Nucleophilic Hydroxide Sites in Oxidation–Resistant Cyclodextrin–Based Metal–Organic Frameworks

Carbon capture and sequestration (CCS) from industrial point sources and direct air capture are necessary to combat global climate change. A particular challenge faced by amine-based sorbents—the current leading technology—is poor stability towards O 2 . Here, we demonstrate that CO 2 chemisorption in γ-cylodextrin-based metal-organic frameworks (CD-MOFs) occurs via HCO 3 – formation at nucleophilic OH – sites within the framework pores, rather than via previously proposed pathways. Here, the new framework KHCO 3 CD-MOF possesses rapid and high-capacity CO 2 uptake, good thermal, oxidative, and cycling stabilities, and selective CO 2 capture under mixed gas conditions. Because of its low cost and performance under realistic conditions, KHCO 3 CD-MOF is a promising new platform for CCS. More broadly, our work demonstrates that the encapsulation of reactive OH – sites within a porous framework represents a potentially general strategy for the design of oxidation-resistant adsorbents for CO 2 capture.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Chelator‐Assisted Precipitation‐Based Separation of the Rare Earth Elements Neodymium and Dysprosium from Aqueous Solutions

The rare earth elements (REEs) are critical resources for many clean energy technologies, but are difficult to obtain in their elementally pure forms because of their nearly identical chemical properties. Here, an analogue of macropa, G-macropa, was synthesized and employed for an aqueous precipitation-based separation of Nd 3+ and Dy 3+ . G-macropa maintains the same thermodynamic preference for the large REEs as macropa, but shows smaller thermodynamic stability constants. Molecular dynamics studies demonstrate that the binding affinity differences of these chelators for Nd 3+ and Dy 3+ is a consequence of the presence or absence of an inner-sphere water molecule, which alters the donor strength of the macrocyclic ethers. Leveraging the small REE affinity of G-macropa, we demonstrate that within aqueous solutions of Nd 3+ , Dy 3+ , and G-macropa, the addition of HCO 3 − selectively precipitates Dy 2 (CO 3 ) 3 , leaving the Nd 3+ −G-macropa complex in solution. With this method, remarkably high separation factors of 841 and 741 are achieved for 50 : 50 and 75 : 25 mixtures. Further studies involving Nd 3+ :Dy 3+ ratios of 95 : 5 in authentic magnet waste also afford an efficient separation as well. Finally, G-macropa is recovered via crystallization with HCl and used for subsequent extractions, demonstrating its good recyclability.

Gao, Yangyang↗
DOE OSTI2024

Structure Activity Relationships for Second‐Coordination Sphere Functional Group Dependent CO 2 Reduction by Manganese Bipyridyl Electrocatalysts

Abstract A series of twelve second coordination sphere (SCS) functionalized manganese tricarbonyl bipyridyl complexes are investigated for their electrocatalytic CO 2 reduction properties in acetonitrile. A qualitative and quantitative assessment of the SCS functional groups is discussed with respect to the catalysts’ thermodynamic and kinetic efficiencies, and their product selectivities. In probing a broad scope of functional groups, it is clear that only the aprotic ortho‐arylester SCS is capable of promoting the highly desired low‐overpotential proton‐transfer electron‐transfer (PT‐ET) pathway for selective CO production. The ortho‐phenolic analogues cause an increase in overpotential with a product selectivity favoring H 2 evolution, consistent with a high‐overpotential pathway via the anionic [Mn−H] − intermediate. Alternative aprotic Lewis base functional groups such as trifluoromethyl, morpholine and acetamide are shown to also be capable of intermediate manganese hydride generation. The tertiary amine substituent, 2‐morpholinophenyl, exhibits a desirable product distribution characteristic of syn‐gas (CO : H 2 =30 : 48) with an impressive turnover frequency, while the secondary amine group, 2‐acetamidophenyl, induces a notable shift in selectivity with a faradaic yield of 55 % for the formate (HCO 2 − ) product. In addition to their catalytic properties, cyclic voltammetry and infrared spectroelectrochemistry (IR‐SEC) studies are presented to probe pre‐catalyst electronic properties and the two‐electron reduction activation pathway.

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

Resolving the Chemical Formula of Nesquehonite via NMR Crystallography, DFT Computation, and Complementary Neutron Diffraction

Nesquehonite is a magnesium carbonate mineral relevant to carbon sequestration envisioned for carbon capture and storage of CO 2 . Its chemical formula remains controversial today, assigned as either a hydrated magnesium carbonate [MgCO 3 · 3H 2 O], or a hydroxy- hydrated- magnesium bicarbonate [Mg(HCO 3 )OH · 2H 2 O]. Here, the resolution of this controversy is central to understanding this material‘s thermodynamic, phase, and chemical behavior. In an NMR crystallography study, using rotational-echo double-resonance 13 C{ 1 H} (REDOR), 13 C- 1 H distances are determined with precision, and the combination of 13 C static NMR lineshapes and density functional theory (DFT) calculations are used to model different H atomic coordinates. [MgCO 3 · 3H 2 O] is found to be accurate, and evidence from neutron powder diffraction bolsters these assignments. Refined H positions can help understand how H-bonding stabilizes this structure against dehydration to MgCO 3 . More broadly, these results illustrate the power of NMR crystallography as a technique for resolving questions where X-ray diffraction is inconclusive.

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

Influence of the double bond position in combustion chemistry of methyl butene isomers: A shock tube and laser absorption study

Shock tube experiments have been carried out on 2–methyl–1–butene (2M1B), 2–methyl–2–butene (2M2B), and 3–methyl–1–butene (3M1B)—the three isomers of methyl butene compound. Here, carbon monoxide (CO) time–histories and ignition delay times are obtained behind reflected shockwaves over the temperature range of 1350–1630 K and pressures of 8.3–10.5 atm with stoichiometric mixtures of 0.075% fuel in O 2 /Ar. Comparative ignition study reveals that 3M1B ignites significantly faster than the other two isomers, while 2M1B dissociates earlier but ignites later than 2M2B. Possible mechanisms for this behavior are discussed with ignition delay time sensitivity and reaction path analysis. In addition, time–resolved CO measurements are compared with three different reaction mechanisms from the literature. Sensitivity analyses have been carried out to identify important reactions that need attention to accurately predict the chemistry of these isomers. Further investigation into the rates of unimolecular fuel decomposition reactions and C 3 H 3 + O 2 = CH 2 CO + HCO reaction are suggested based on the current investigation.

09 BIOMASS FUELS↗
DOE OSTI2021

Extended thermodynamic model for high salinity produced waters

To support the development of desalination processes for high salinity produced waters, we developed a comprehensive thermodynamic model that reliably predicts the thermodynamic properties and phase behavior of produced water. Here this study extends a previously developed thermodynamic model for the aqueous hexary oceanic salt system by incorporating the Sr 2+ and Ba 2+ ions. Based on the electrolyte nonrandom two-liquid theory, the model requires two adjustable binary interaction parameters for each water-electrolyte and electrolyte-electrolyte pair that has a common ion. The binary interaction parameters for the electrolyte-electrolyte pairs involving Sr 2+ and Ba 2+ ions were identified using thermodynamic and salt solubility data. The model was validated for temperatures from 273.15 K to 473.15 K and electrolyte concentrations from infinite dilution to salt saturation. Extension of the model to include HCO 3 − , CO 3 2− , and CO 2 is currently in progress.

42 ENGINEERING↗
DOE OSTI2023

Formation of hydrocarbons favored by high pressure at subduction zone conditions

Subduction zones enable carbon transport to the deep Earth. Inorganic and organic carbon in fluids released from the subducted slab are important in water-rock interactions. However, the potential role of reduced carbon species in fluids during subduction is still unclear. In particular, the pressure, temperature, and oxidation state conducive to the formation of fluid hydrocarbons are poorly understood. Here we focused on the role of pressure influencing the formation of hydrocarbons using diamond anvil cell experiments starting with 0.5 M Ca-acetate solutions over a broad range of pressures, from 1.6 to 4.6 GPa at temperatures of 300 and 350°C. Droplets of immiscible hydrocarbon fluid coexisting with aqueous carbonate or bicarbonate and crystals of calcite or aragonite were formed at peak pressure and temperature conditions in the experiments. Furthermore, greater quantities of hydrocarbons formed at higher pressure as revealed by measuring the in situ volume percent of the droplets. Measurement of the in situ Raman spectra of the co-existing aqueous solutions and analysis of the peak area ratios of acetate to HCO 3 – showed faster decrease at higher pressure, which indicated that the reactions forming hydrocarbons were accelerated by pressure. Compared to previous Na-acetate experiments at 300°C and 3.0 GPa, this study found hydrocarbons that were more methane-rich and had more propane and less isobutane. In experiments at 350°C, isobutane transformed to one or more aromatic hydrocarbons in an immiscible fluid. The air-dried aromatic hydrocarbons were measured using a UV laser, and show peaks of bitumen. Overall, our study supports the hypothesis that high pressures and high temperatures facilitate the possible occurrence of immiscible hydrocarbon fluids in the deep carbon cycle.

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

Automated identification and calculation of prompt effects in kinetic mechanisms using statistical models

The kinetics of prompt dissociation involves rovibrationally excited species (generally formed by exothermic reactions) which may dissociate or isomerize prior to thermalization via collisions with the bath gas. Treating such rovibrationally excited species (so-called "hot" species) with standard kinetic phenomenology may result in incorrect macroscopic representation of their reactivity. Here this work presents the first fully automated methodology for the calculation of prompt effects of a chosen species in a kinetic mechanism, including (i) reaction selection; (ii) theoretical calculation of rate constants and prompt branching fractions; and (iii) final rate constant fitting. The energy partition between hot fragments is estimated using a variety of statistical models, including a new physically sound microcanonical statistical model based on the rovibrational density of states of the fragments. The methodology is validated against literature data for the prompt dissociations of HCO and C 3 H 7 radicals. The microcanonical statistical model is in better agreement with trajectory simulations for larger species and is thus applicable for practical systems that typically involve large molecules, for which direct dynamics calculations are impractical. The automated workflow is applied to the evaluation of the effects of prompt dissociation for two isomeric radicals C 4 H 7 1-3 (1-methylallyl) and C 4 H 7 1-4 (3-buten-1-yl). Twelve H-atom abstraction reactions are selected and the corresponding rate constants are computed with first principles theory. The microcanonical statistical model predicts that prompt dissociations of C 4 H 7 1-3 and C 4 H 7 1-4 are already significant at 1000K, resulting in differences of up to an order of magnitude at 2000K with respect to the phenomenological thermal rate constants. To illustrate the effects of prompt dissociation on simulations of experimental data, the calculated prompt rate constants are implemented in both CRECK and C3MechV3.3 kinetic mechanisms. Simulations of experimental flame data illustrate the noticeable impact of prompt dissociation kinetics on the high-temperature combustion reactivity of C 4 H 8 -1 and C 4 H 8 -2.

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

Reduced chemical kinetic model for CH 4 -air non-premixed flames including excited and charged species

Electric fields can impact small laminar flames by changing their shape and overall behavior by acting on charged species produced in combustion. However, no reduced chemical kinetic model has been developed considering both major species and minor species related to flame characterization and flame behavior in the presence of an electric field. Here, this study presents a reduced chemical kinetic model for methane-air combustion which includes minor excited species (CH*and OH*) and charged species (H 3 O + , HCO + , C 2 H 3 O + , CH 5 O + , O$_{2}^{–}$, OH – , e – , CO$_{3}^{–}$, CHO$_{2}^{–}$, O – , CHO$_{3}^{–}$). The results employing the reduced chemistry model have been validated for a two-dimensional flame geometry by comparison with (i) detailed chemistry simulation results for species location and peak values, and (ii) experimental CH* chemiluminescence location, considering the self-repulsion of charges yet without externally applied electric field to the flame. This reduced chemical kinetic model, with 45 species and 216 reactions, shows a computational demand one-third that of employing its equivalent detailed chemistry (83 species and 394 reactions). The reduction is modest but significant considering that high fidelity is needed to capture the behavior of the chemi-ion and chemiluminescent species. Future works will involve the use of this model for simulations predicting flame behavior with applied electric field (i.e., field strength ≠0 kV/cm).

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

Methyl formate oxidation kinetics up to 100 atm

Methyl formate (MF, CH3OCHO), the simplest ester, is a representative oxygenated fuel with high oxygen content, and low sooting tendency. However, its oxidation behavior under high-pressure and intermediate-temperature conditions remains insufficiently understood, especially where low-temperature peroxy radical chemistry, methanol chemistry, and pressure-dependent reaction pathways play a critical role. In this study, MF oxidation experiments were conducted in the Princeton supercritical-pressure jet-stirred reactor (SP-JSR) at 20 and 100 atm over the temperature range of 400–950 K under both fuel-lean and fuel-rich conditions. Based on the experimental results, an updated HP-Mech was developed by incorporating previous MF sub-mechanisms, expanded low-temperature peroxy pathways, and evaluated pressure-dependent decomposition kinetics. The newly updated HP-Mech shows greatly improved performance in predicting the onset temperature, the key intermediate species fractions, methanol formation, and the progression of MF oxidation across all the experimental conditions. Path flux analysis indicates that MF consumption at the onset stage is dominated by H-abstraction at the methyl site, forming CH2OCHO radicals that lead to the formation and isomerization of O2CH2OCHO, driving low-temperature chain propagation. Moreover, H-abstraction at the formate site forms CH3OCO radicals that preferentially decompose to CH3, initiating the methanol formation pathway linked to CH3O2 and HO2 chemistry. At the same time, HO2 formation is strongly coupled to MF oxidation through multiple MF-derived radical pathways. HCO originates from MF oxidation and acts as a key coupling species linking fuel consumption to HO2 buildup, especially under high-pressure and intermediate-temperature conditions. In addition to this dominant channel, supplementary HO2 formation pathways involving CH3, CH3O, CH2OH, and CH3O2 reacting with O2 further connect methanol chemistry and oxygenated radical chemistry to the HO2 pool, indicating the central role of HO2 in governing MF oxidation. Sensitivity analysis identifies MF with OH/HO2/CH3O2 reactions and the HO2/H2O2/OH sequence as the key factors controlling reactivity in the high-pressure and intermediate-temperature regime. MF directly reacts with OH/HO2/CH3O2 to consume the fuel and produce reactive radicals like CH2OCHO and CH3OCO that undergo subsequent oxidation pathways. Moreover, HO2 recombination suppresses oxidation at lower temperatures, while thermal decomposition of H2O2 accelerates OH production and promotes fuel consumption as temperature increases. The direct formation of active OH from HO2 radicals further completes the mechanism, improving its prediction especially during the oxidation onset stage.

Low-temperature Chemistry↗
DOE OSTI2024

Coupled feldspar dissolution and secondary mineral precipitation in batch systems: 6. Labradorite dissolution, calcite growth, and clay precipitation at 60 °C and pH 8.2–8.4

Here, we conducted experiments on concurrent labradorite dissolution, calcite precipitation, and clay precipitation in batch reactor systems and tracked reaction processes using multiple isotope tracers. Labradorite was chosen for its role as a major and reactive component in basalt; the experiments thus directly impact our understanding of CO 2 storage in basalt aquifers and enhanced rock weathering. We doped initial solutions with 29 Si, 43 Ca, and Ca 13 CO 3 (s). Experiments were conducted at 60 °C and pH ~ 8.3 for up to 840 h, with isotope ratios in the experimental aqueous solutions measured using MC-ICP-MS. Unidirectional rates of labradorite dissolution near equilibrium were approximately two orders of magnitude slower than far-from-equilibrium rates reported in the literature. Calcite growth occurred near equilibrium and the rates were limited by the labradorite dissolution rates. In the steady state phase, the interplay of these three heterogeneous reactions—labradorite dissolution, calcite growth, and clay precipitation—results in a coupled system that approaches a near-equilibrium state. The system does not reach true equilibrium because labradorite continues to dissolve, albeit at a much slower rate near equilibrium. The overall reaction can be approximated as, Na 0.4 Ca 0.6 Al 1.6 Si 2.4 O 8 + 0.6HCO 3 - + 1·.7H 2 O + 0.4H + → 0.4Na + + 0.6CaCO 3(s) + 0.5Al 2 Si 2 O 5 (OH) 4(s) + 0.6Al(OH) 4 - + 1.4SiO 2 o (aq). The experimental results show that using short-term far-from-equilibrium rate constants would lead to an overestimation of feldspar weathering rates at the Earth’s surface (e.g., basalt weathering and enhanced rock weathering) and CO 2 mineralization in basalt aquifers.

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

Machine learning to discover mineral trapping signatures due to CO 2 injection

Mineral trapping is pursued as a geological CO 2 sequestration (GCS) mechanism because it permanently stores CO 2 in solid phases or minerals. However, CO 2 mineral-trapping mechanisms are poorly understood due to (1) lack of sufficient field and laboratory data characterizing these complex processes, and (2) challenges to develop site-specific reactive-transport models coupling fluid flow and geochemical reactions occurring at various temporal (from milliseconds to years) and spatial (from pore (millimeters) to field (kilometers)) scales. Reactive transport with additional complexities such as heterogeneity can make the simulation outputs even more difficult to interpret because of complex nonlinearity and multi-scale interdependencies. Furthermore, the values of model outputs such as concentrations can vary by several orders of magnitude, making it harder to correlate and characterize the impact of the variables via traditional data interpretation techniques such as exploratory data analyses. Recently, machine learning (ML) has shown promise in feature discovery and in highlighting hidden mechanisms that cannot be obtained by existing data-analytics and statistical methods. In this study, we applied an unsupervised ML approach, non-negative matrix factorization with custom -means clustering (NMF) to the data generated by reactive-transport simulations of GCS. The reactive-transport data consisted of 19 attributes, including four physio-chemical variables (pH, porosity, aqueous CO 2 , and sequestered CO 2 ), six chemical species (K + , Na + , HCO, Ca 2+ , Mg 2+ , Fe 2+ ), and four carbonate minerals (calcite, dolomite, siderite, and ankerite), a feldspar mineral (albite), and four clay minerals (illite, clinochlore, kaolinite, and smectite) over a period of 200 years of simulation time. Furthermore, the simulation data used was for Morrow B sandstone at the Farnsworth hydrocarbon unit in Texas. Data are sampled at two locations within the model domain: (1) at the injection well and (2) 200 m west of the injection well. The injection was performed for a period of 10 years. Using NMF, we estimated the temporal interdependencies among the 19 attributes over a span of 200 years. We found that NMF was able to identify four reaction stages and their dominant attributes; these cannot be directly discerned through traditional visualization (e.g., line plots, Pareto analysis, Glyph-based visualization methods) or exploratory data analysis tools of the simulation data. The four stages were: reactions in the injection phase followed by short-, mid-, and long-term reactions. The NMF analysis also revealed that 10 among the 19 attributes are dominant. These dominant attributes for mineral trapping include calcite, dolomite at injection well, siderite at 200 m away from the injection well, clinochlore, kaolinite, Na + , K + , Ca 2+ , Mg 2+ , pH, and aqeuous CO 2 . Finally, at late times (65–200 years), our results showed that calcite plays a major role in mineral trapping with insignificant contribution from siderite, ankerite, and clay minerals. These findings make the proposed unsupervised ML-model attractive for reactive-transport sensing towards real-time GCS monitoring.

54 ENVIRONMENTAL SCIENCES↗