Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geochemical Modelling”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Unraveling Hydrogen Induced Geochemical Reaction Mechanisms through Coupled Geochemical Modeling and Machine Learning

Underground hydrogen storage (UHS) provides a promising large-scale, long-term energy storage solution. A reasonable recovery of stored hydrogen is critical for a successful storage scheme. However, in subsurface reservoirs hydrogen is subject to active geochemical reactions that might result in hydrogen loss. In this study, we implemented a geochemical modeling approach coupled with an unsupervised machine learning technique called non-negative matrix factorization (NMF) to unravel the complex brine-rock-H 2 geochemical processes responsible for hydrogen losses, with particular focus on sulfate reduction reactions. NMF is applied to modeled mineral evolution and fluid component profiles to retrieve profiles that can be interpreted to more easily assess competing processes. NMF decouples simulated competing equilibrium reactions. This facilitates separation of overlapping reaction profiles from redox processes, dissolution fronts, and secondary precipitation while considering the effects of simulation parameters such as salinity, temperature, and total H 2 pressure. NMF successfully discriminates these competing effects in nonlinear ways, allowing robust interpretation. In addition, NMF reveals subtle coupled mineral associations and reaction fronts that are invisible to conventional model analysis. This integrated approach strengthens the conceptual understanding of complex nonlinear hydrogen-brine-rock interactions and advances geochemical research on UHS systems to resolve complexities in modeled geochemical systems without the need for direct experiments or prior knowledge. Furthermore, this study highlights the efficacy of combining geochemical modeling with machine learning techniques to enhance the interpretability of the intricate geochemical simulation output through deciphering the overlapping reaction path that cannot be achieved only using conventional analysis of geochemical models alone.

08 HYDROGEN↗
DOE OSTI2023

Geochemical Modeling of the Confining Unit at a CO 2 Storage Site in the Williston Basin, North Dakota

Conference presentation for American Association of Petroleum Geologists (AAPG) Carbon Capture, Utilization, and Storage (CCUS) Conference, Houston, TX, April 25–27, 2023. This geochemical modeling study investigates the evolution of geochemical reactions among CO 2 , formation water, and rocks of the two confining units above and below the CO 2 storage unit at one of the permitted CO 2 storage sites in the Williston Basin. The model is built on PHREEQC-coded chemical equilibrium thermodynamics.

20 FOSSIL-FUELED POWER PLANTS↗
NASA NTRS1997

Hydrothermal hydration of Martian crust: illustration via geochemical model calculations

If hydrothermal Systems existed on Mars, hydration of crustal rocks may have had the potential to affect the water budget of the planet. We have conducted geochemical model calculations to investigate the relative roles of host rock composition, temperature, water-to-rock ratio, and initial fluid oxygen fugacity on the mineralogy of hydrothermal alteration assemblages, as well as the effectiveness of alteration to store water in the crust as hydrous minerals. In order to place calculations for Mars in perspective, models of hydrothermal alteration of three genetically related Icelandic volcanics (a basalt, andesite, and rhyolite) are presented, together with results for compositions based on SNC meteorite samples (Shergotty and Chassigny). Temperatures from 150 degrees C to 250 degrees C, water-to-rock ratios from 0.1 to 1000, and two initial fluid oxygen fugacities are considered in the models. Model results for water-to-rock ratios less than 10 are emphasized because they are likely to be more applicable to Mars. In accord with studies of low-grade alteration of terrestrial rocks, we find that the major controls on hydrous mineral production are host rock composition and temperature. Over the range of conditions considered, the alteration of Shergotty shows the greatest potential for storing water as hydrous minerals, and the alteration of Icelandic rhyolite has the lowest potential.

NASA Discipline Exobiology↗
NASA NTRS1995

A geochemical model for the formation of hydrothermal carbonates on Mars

It is often argued that substantially more carbon dioxide and water were degassed from the martian interior than can be found at present in the atmosphere, polar caps and regolith. Calculations have shown that atmospheric escape cannot account for all of the missing volatiles. Suggestions that carbon dioxide is stored as marine or lacustrine deposits, are challenged by Earth-based and spacecraft remote-sensing data. Moreover, recent modelling of the martian atmosphere suggests that rainfall or open bodies of water are in any case unlikely to have persisted for extended periods of time. Hydrothermal carbonates therefore provide a possible solution to this dilemma. Using an accessible terrestrial system (Iceland) as a guide to the underlying processes, and a host rock composition inferred from the least-altered martian meteorite, we present a geochemical model for the formation of carbonates in possible martian hydrothermal systems. Our results suggest that an extensive reservoir of carbonate minerals--equivalent to an atmospheric pressure of carbon dioxide of at least one bar--could have been sequestered beneath the surface by widespread hydrothermal activity in the martian past.

NASA Discipline Exobiology↗
NASA NTRS1980

Equilibrium geochemical modeling of a seasonal thermal energy storage aquifer field test

A geochemical mathematical modeling study designed to investigate the well plugging problems encountered at the Auburn University experimental field tests is summarized. The results, primarily of qualitative interest, include: (1) loss of injectivity was probably due to a combination of native particulate plugging and clay swelling and dispersion; (2) fluid-fluid incompatibilities, hydrothermal reactions, and oxidation reactions were of insignificant magnitude or too slow to have contributed markedly to the plugging; and (3) the potential for and contributions from temperature-induced dissolved gas solubility reductions, capillary boundary layer viscosity increases, and microstructural deformation cannot be deconvolved from the available data.

Stottlemyre, J. S.↗
DOE OSTI2023

Geochemical Modeling Using OLI Systems Mixed Solvent Electrolyte Model in Support of the Waste Isolation Pilot Plant

Aqueous electrolyte thermodynamic models were developed for synthetic Salado (GWB) and Castile (ERDA-6) brines to inform the experimental design and to facilitate laboratory data interpretation in geochemical studies supporting disposal in the Waste Isolation Pilot Plant (WIPP). The thermodynamic models, which include the relevant vapor-liquid-solid equilibria for the oxidation-reduction conditions of interest, enable computation of system parameters such as activity-based pH, E h [vs. standard hydrogen electrode (SHE)], ionic strength, total dissolved solids, and solution density; the chemical speciation of all phases; saturation tendencies for thermodynamically feasible solids; and changes in brine chemistry because of the addition of H 2 (g), argon, MgO, Fe(0), and other waste form components. The impacts of oxidation-reduction potential (ORP) and pH on aqueous speciation and dominant solid phases were also assessed with the aid of real-solution Pourbaix diagrams for key elements of interest. In addition to the OLI modeling results, limited measurement data for pH and ORP standards are provided, analyzed, and discussed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1979

A global geochemical model for the evolution of the mantle

It is proposed that the upper mantle transition region, 220 to 670 km, is composed of eclogite which has been derived from primitive mantle by about 20 percent partial melting and that this is the source and sink of oceanic crust. The remainder of the upper mantle is garnet peridotite which is the source of continental basalts and hotspot magmas. This region is enriched in incompatible elements by hydrous and CO2 rich metasomatic fluids which have depleted the underlying layers in the L.I.L. elements and L.R.E.E. The volatiles make this a low-velocity, high attenuation, low viscosity region. The eclogite layer is internally heated and its controls the convection pattern in the upper mantle. Plate tectonics is intermittent. The continental thermal anomaly at a depth of 150-220 km triggers kimberlite and carbonatite activity, alkali and flood basalt volcanism, vertical tectonics and continental breakup. Hot spots remain active after the continents leave and build the oceanic islands. Mantle plumes rise from a depth of about 220 km. Midocean ridge basalts rise from the depleted layer below this depth. Material from this layer can also be displaced upwards by subducted oceanic lithosphere to form back-arc basins.

Anderson, D. L.↗
NASA NTRS2007

The Aqueous Alteration of CR Chondrites: Experiments and Geochemical Modeling

CR carbonaceous chondrites are of major interest since they contain some of the most primitive organic matter known. However, aqueous alteration has more or less overprinted their original features in a way that needs to be assessed. This study was initiated by comparing the mineralogy and modal abundances of the most altered CR1 chondrite, GRO 95577, to a less altered CR2. Calculated element distributions imply that GRO 95577 may result from aqueous alteration of Renazzo by an isochemical process on their parent asteroid, whose mineralogical composition was estimated ( Unaltered CR shown included table).

Perronnet, M.↗
DOE OSTI2022

Integrated Experimental and Modeling Study of Geochemical Reactions of Simple Fracturing Fluids with Caney Shale

Interactions between rock minerals and hydraulic fluids directly impact the geochemical and geomechanical properties of shale formations. However, the mechanisms of geochemical reactions in shale unconventional reservoirs remain poorly understood. Here, to investigate the geochemical reactions between shale and hydraulic fracturing fluids, a series of batch reactor experiments were undertaken. Three rock samples with different mineralogical compositions and three fluid samples of different compositions (deionized water, deionized water + 2% potassium chloride (KCl), and deionized water + 0.5% choline chloride (C 5 H 14 ClNO) were used. Experiments were undertaken at reservoir temperature and atmospheric pressure. Elemental compositions of effluents after 1, 3, 7, 14, 28 days were analyzed using Inductively Couple Plasma Mass Spectrometry (ICP-MS). Medical Computed Tomography (CT) scan and X-ray Fluorescence (XRF) were conducted on the entire core run to help upscale results obtained from rock-fluid interaction experiments. Geochemical modeling using a reactive simulator, TOUGHREACT, was undertaken to corroborate experimental results. Results show that lower pH triggered high dissolution rates in the rock samples, especially the carbonate components. As pH increased, the rate of dissolution declined significantly, though for most cases dissolution still continued. Observed dissolved silica concentrations were much higher than the quartz solubility, suggesting that much of the silica originates from more soluble silica polymorphs and possibly desorption from clay mineral exchange sites. Concentration of most elemental species in solution increased but aluminium and magnesium concentrations declined rapidly following initial entry into solution. Geochemical modeling corroborated the conclusions regarding mineral dissolution and precipitation observed from experiments, notably; the dissolution of calcite and oxidation of pyrite in reacted shale samples, the likely presence of silica polymorphs such as opal, chalcedony or amorphous silica in these samples, and the reduction of Al and Mg concentrations in solution by precipitation of secondary aluminosilicate phases. The de-flocculation of clay minerals during reaction implies fines migration after hydraulic fracturing. This is detrimental to reservoir productivity as clay fines are displaces and lodged within the micro and nano-fractures created during fracturing. The immediate consumption of aluminium and magnesium also has implications on blockage of hydrocarbon pathways due to precipitation of new minerals in these locations.

04 OIL SHALES AND TAR SANDS↗
NASA NTRS1992

A carbonate-silicate aqueous geochemical cycle model for Mars

A model for the carbonate-silicate geochemical cycle of an early, wet Mars is under development. The results of this study will be used to constrain models of the geochemical history of Mars and the likely mineralogy of its present surface. Although Mars today is a cold, dry planet, it may once have been much warmer and wetter. Values of total outgassed CO2 from several to about 10 bars are consistent with present knowledge (Pollack et al. 1987), and this amount of CO2 implies an amount of water outgassed at least equal to an equivalent depth of 500-1000 meters (Carr 1986). Pollack et al. (1987), in addition, estimate that a thick CO2 atmosphere may have existed for an extended period of time, perhaps as long as a billion years. The greenhouse effect of such an atmosphere would permit the presence of liquid water on the surface, most likely in the form of a shallow sea in the lowest regions of the planet, such as the northern plains (Schaefer 1990). The treatment of geochemical cycles as complex kinetic chemical reactions has been undertaken for terrestrial systems in recent years with much success (Lasaga 1980, 1981; Berner et al. 1983; Lasaga et al. 1985). Although the Martian system is vastly less well understood, and hence less well-constrained, it is also a much simpler system, due to the lack of biogenic reactions that make the terrestrial system so complex. It should be possible, therefore, to use the same techniques to model the Martian system as have been used for terrestrial systems, and to produce useful results. A diagram of the carbonate-silicate cycle for Mars (simplified from the terrestrial system) is given.

Schaefer, M. W.↗
NASA NTRS1976

Chronology of the early lunar crust

The Rb-Sr systematics of lunar 'anorthositic' rocks was studied in an effort to understand the evolution of the lunar crust. The contamination by KREEP of the anorthositic crust is corrected using chemical mixing models; cumulate geochemical models are used to establish the chronology of the early lunar crust; a method is presented to derive an estimate of the Rb/Sr ratio of the moon. A mathematical model for the Rb-Sr evolution of the lunar crust that includes mare basalt genesis is presented.

Schonfeld, E.↗
NASA NTRS2024

Revision 1: Characterizing Basalt-Atmosphere Interactions on Venus: A Review of Thermodynamic and Experimental Results

The surface of Venus is in contact with a hot (~470° C), high pressure (92 bars), and caustic (CO 2 with S, but little H 2 O) atmosphere, which should cause progressive alteration of the crust in the form of sulfate and iron-oxide coatings; however, the exact rate of alteration and mineral species are not well constrained. Different experimental approaches, each with its own limitations, are currently being used to constrain mineralogy and alteration rate. One note is that no experimental approach has been able to fully replicate the necessary conditions and sustain them for a significant length of time. Further, geochemical modeling studies can also constrain surface alteration mineralogy, again with different assumptions and limitations. Here we review recent geochemical modeling and experimental studies to constrain the state of the art for alteration mineralogy, rate of alteration, open questions about the surface mineralogy of Venus, and what can be constrained before the fleet of missions arrive later this decade. Combining the new results confirm that basalt on the surface of Venus should react quickly and form coatings of sulfates and iron-oxides; however, the mineralogy and rate of alteration is dependent on physical properties of the protolith (including bulk composition, mineralogy, and crystallinity), as well as atmospheric composition, and surface temperature. Importantly, the geochemically modelling results show that the mineralogy is largely controlled by atmospheric oxygen fugacity, which is not well constrained for the near surface environment on Venus. Therefore, alteration experiments run over a range of oxygen and sulfur fugacities are needed across a wide range of Venus analog materials with varying mineralogy and crystallinity.

Justin Filiberto↗
DOE OSTI2022

Stable carbon isotope values of syndepositional carbonate spherules and micrite record spatial and temporal changes in photosynthesis intensity

Marine and lacustrine carbonate minerals preserve carbon cycle information, and their stable carbon isotope values (δ 13 C) are frequently used to infer and reconstruct paleoenvironmental changes. However, multiple processes can influence the δ 13 C values of bulk carbonates, confounding the interpretation of these values in terms of conditions at the time of mineral precipitation. Co-existing carbonate forms may represent different environmental conditions, yet few studies have analyzed δ 13 C values of syndepositional carbonate grains of varying morphologies to investigate their origins. Here, we combine stable isotope analyses, metagenomics, and geochemical modeling to interpret δ 13 C values of syndepositional carbonate spherules (>500 μm) and fine-grained micrite (<63 μm) from a ~1600-year-long sediment record of a hypersaline lake located on the coral atoll of Kiritimati, Republic of Kiribati (1.9°N, 157.4°W). Petrographic, mineralogic, and stable isotope results suggest that both carbonate fractions precipitate in situ with minor diagenetic alterations. The δ 13 C values of spherules are high compared to the syndepositional micrite and cannot be explained by mineral differences or external perturbations, suggesting a role for local biological processes. We use geochemical modeling to test the hypothesis that the spherules form in the surface microbial mat during peak diurnal photosynthesis when the δ 13 C value of dissolved inorganic carbon is elevated. In contrast, we hypothesize that the micrite may precipitate more continuously in the water as well as in sub-surface, heterotrophic layers of the microbial mat. Both metagenome and geochemical model results support a critical role for photosynthesis in influencing carbonate δ 13 C values. The down-core spherule–micrite offset in δ 13 C values also aligns with total organic carbon values, suggesting that the difference in the δ13C values of spherules and micrite may be a more robust, inorganic indicator of variability in productivity and local biological processes through time than the δ 13 C values of individual carbonate forms.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2002

Emergence of RNA and Early Life; Experimental Modeling and Geochemical Constraints

Support from this grant continues to fund our research on two related problems. One involves attempts to model the abiotic formation of simple source compounds for functional biomolecules, their concentration from dilute state in the hydrosphere and in several cases surface induced reactions to form precursor monomers for bioactive end products. The second area is a search for the earliest traces of life on Earth and the early environments on Earth and Mars. For this purpose we have analyzed the isotopic compositions of carbon and nitrogen in graphitic residues that have been thought to represent the remains of microorganisms. The results of these studies have been presented in a number of publication, listed as references.

Arrhenius, Gustaf↗