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

Simulated Hydrological Dynamics and Coupled Iron Redox Cycling Impact Methane Production in an Arctic Soil

The fate of organic carbon (C) in permafrost soils is important to the climate system due to the large global stocks of permafrost C. Thawing permafrost can be subject to dynamic hydrology, making redox processes an important factor controlling soil organic matter (SOM) decomposition rates and greenhouse gas production. In iron (Fe)-rich permafrost soils, Fe(III) can serve as a terminal electron acceptor, promoting anaerobic respiration of SOM and increasing pH. Current large-scale models of Arctic C cycling do not include Fe cycling or pH interactions. Here, a geochemical reaction model was developed by coupling Fe redox reactions and C cycling to simulate SOM decomposition, Fe(III) reduction, pH dynamics, and greenhouse gas production in permafrost soils subject to dynamic hydrology. In this study we parameterized the model using measured CO 2 and CH 4 fluxes as well as changes in pH, Fe(II), and dissolved organic C concentrations from oxic and anoxic incubations of permafrost soils from polygonal permafrost sites in northern Alaska, United States. In simulations of repeated oxic-anoxic cycles, Fe(III) reduction during anoxic periods enhanced CO 2 production, while the net effect of Fe(III) reduction on cumulative CH 4 fluxes depended on substrate C availability. With lower substrate availability, Fe(III) reduction decreased total CH 4 production by further limiting available substrate. With higher substrate availability, Fe(III) reduction enhanced CH 4 production by increasing pH. Our results suggest that interactions among Fe-redox reactions, pH and methanogenesis are important factors in predicting CH 4 and CO 2 production as well as SOM decomposition rates in Fe-rich, frequently waterlogged Arctic soils.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1985

The relation between the age of the subconducting slab and the recycling of sediments into the mantle

The recycling of sediments into the mantle has become an important issue because recent papers have suggested that the geochemical inverse models of the evolution of radiogenic isotope abundances over the history of the Earth have nonunique solutions. Both the recycling of continent-derived sediments into the mantle and mixing in the mantle could produce similar geochemical effects in the mean isotopic ratios of new igneous material emplaced in continents. Recent models of Archaean heat flow and of plate tectonics during early Earth history have demonstrated that higher internal heat production of the early Earth was mainly dissipated through a higher creation rate of oceanic lithosphere. If the seafloor creation rate was higher on the early Earth, then the residence time of any one piece of oceanic lithosphere on the surface would have been shorter. It is possible that a higher rate of recycling of oceanic lithosphere into the mantle could have resulted in some transport of sediment into the mantle.

Abbott, D.↗
DOE OSTI2021

Drying model of a high salt content cementitious waste form: Effect of capillary forces and salt solution

Highlights: • Drying model for a high salt content cementitious waste form is developed. • Water vapor diffusion and capillary liquid flow are distinguished. • Capillary and salt solution effects are considered in vapor-liquid equilibrium. A water transport model coupling capillary liquid flow with vapor diffusion is developed to describe the drying process for a cementitious waste form with high salinity porewater. Vapor-liquid equilibrium is formulated as the driving force for vapor diffusion and the model accounts for pore capillary and high salinity effects on water thermodynamic activity. Pore filling and porewater surface tension as a function of pore size distribution and water saturation have been quantified for the material. Geochemical speciation modeling is used to simulate porewater activity as a function of composition over the range of saturation. The theoretical relationship between relative humidity and water saturation generally agrees with experimental measurement, and the developed model is capable of predicting drying rates under various external relative humidity conditions. The model was developed to be incorporated into reactive transport models considering the effects of drying such as salt redistribution and efflorescence.

36 MATERIALS SCIENCE↗
NASA NTRS1996

Hydrothermal systems as environments for the emergence of life

Analysis of the chemical disequilibrium provided by the mixing of hydrothermal fluids and seawater in present-day systems indicates that organic synthesis from CO2 or carbonic acid is thermodynamically favoured in the conditions in which hyperthermophilic microorganisms are known to live. These organisms lower the Gibbs free energy of the chemical mixture by synthesizing many of the components of their cells. Primary productivity is enormous in hydrothermal systems because it depends only on catalysis of thermodynamically favourable, exergonic reactions. It follows that hydrothermal systems may be the most favourable environments for life on Earth. This fact makes hydrothermal systems logical candidates for the location of the emergence of life, a speculation that is supported by genetic evidence that modern hyperthermophilic organisms are closer to a common ancestor than any other forms of life. The presence of hydrothermal systems on the early Earth would correspond to the presence of liquid water. Evidence that hydrothermal systems existed early in the history of Mars raises the possibility that life may have emerged on Mars as well. Redox reactions between water and rock establish the potential for organic synthesis in and around hydrothermal systems. Therefore, the single most important parameter for modelling the geochemical emergence of life on the early Earth or Mars is the composition of the rock which hosts the hydrothermal system.

Non-NASA Center↗
DOE OSTI2025

Development of a kinetic-thermodynamic model for lime-stabilization of Na-bentonite

This study presents the first kinetic model to predict the solid and pore solution composition of Na-bentonite clay reacting with slaked lime over a period of 720 days. The model successfully accounts for most experimental data using a single kinetic rate constant. The following sequence of reactions was predicted by the model: initial rapid dissolution of portlandite within the first 7 days, leading to a decrease in pH and dissolved calcium, and concurrent formation of calcium silicate hydrates (C-S-H: jennite), calcium aluminate hydrate (C-A-H: C₄AH₁₃), calcium aluminosilicate hydrates (stratlingite) and hydrotalcite. After 7 days, jennite and stratlingite are predicted to transform into tobermorite-II, contributing to strength development up to 28 days. From 28 to 90 days, continued montmorillonite dissolution is predicted, along with minor formation of ettringite, partial tobermorite-II dissolution, and precipitation of secondary phases such as albite and talc. Experimentally, portlandite dissolution was confirmed by TGA and XRD and found to be complete within 7 days, in agreement with model predictions. However, other predicted solid-phase transformations (e.g., tobermorite-II formation and dissolution, ettringite, albite, and talc formation) could not be conclusively verified through experimental techniques. Aqueous phase measurements confirmed that the pH and Ca trends in solution, and that equilibrium was reached by 90 days.

Chemical kinetics↗
DOE OSTI2024

Effect of Anoxic Iron Corrosion on WIPP Brine Geochemistry FY23 Final Report (U)

A 280-day study was completed to evaluate the effect of zero-valent iron (Fe 0 ) on the Waste Isolation Pilot Plant (WIPP) brine geochemistry under anticipated reducing conditions. Hydrogen (H 2 ) gas is expected to be present in the repository after closure due to the anoxic corrosion of a vast quantity of iron contained in the waste forms disposed at WIPP; therefore, a background argon atmosphere containing H 2 was chosen for this study. WIPP groundwater brine pH and E h will impact the mobility and fate of plutonium within the repository. Modeling and laboratory results for Castile WIPP brine indicate that equilibrium fa values relative to the standard hydrogen electrode (SHE) are 40 mV more reducing (i.e., more negative) than those for Salado WIPP brine (-480 mV vs. -440 mV, respectively) because of the higher pH of the Castile brine (pH 9 .3 for Castile vs. pH 8.8 for Salado). The E h and pH data were corrected for the effects of high ionic strength. The experimental results for both brines are consistent with thermodynamic predictions using OLI Systems' Mixed Solvent Electrolyte chemical equilibrium model. The measured and corrected pH and E h data from this study are provided in Table ES-I and Table ES-2, respectively. The experimental study, with four test conditions in triplicate, was performed in a dual glovebox with a nominally 3 vol.% H 2 in argon atmosphere (target H 2 range: 3 ± I vol.%). Simulants containing MgO only ( experimental control) and MgO+Fe 0 (WIPP base case) were prepared for both the Salado and Castile brines. MgO was included in all simulants to account for the use of bulk magnesium oxide in the WIPP repository. Fe 0 was included in some simulants to incorporate the effects of the anoxic corrosion of iron and in-situ hydrogen generation in the study. The brine compositions were developed by Sandia National Laboratory (SNL; Xiong, 2008) and have been used in previous WIPP evaluations. The test method (agitation, etc.) is partially based on ASTM D3987-12. Twelve rounds of periodic measurements of pH and E h were performed over the course of the study. Chemical analysis results for liquids and solids (ICP-MS, ICP-ES, IC Anion, TIC, SEM-EDX) are consistent with the pH, E h , and thermodynamic modeling results. This study included the following conditions that deviate from anticipated post-closure conditions following brine intrusion, but were selected to facilitate bench-scale testing to validate modeling of pH and E h for the post-closure WIP P repository: an anoxic glove box atmosphere containing ≤ 4 vol. % H 2 vs. substantially higher H 2 gas concentrations assumed in the WIPP Performance Assessment (PA); a significantly higher liquid-to-solid test ratio compared to the much lower phase ratio anticipated in the WIP P repository; agitation of the simulant bottles to maximize mass transfer; and finally the use of Fe 0 reagents having a much greater surface area than expected in the WIP P repository. Non-representative conditions were chosen for various reasons such as: to provide bounding conservative results, to provide a margin of safety for testing, or to facilitate simulant sub-sampling and analysis. In a parallel effort, aqueous electrolyte thermodynamic models were developed for the synthetic Salado and Castile brines to inform the experimental design, facilitate laboratory data interpretation, and allow extension of evaluations beyond the parameters tested. Thermodynamic modeling simulations including the MgO and Fe 0 additives that are directly relevant to the experimental measurements (e.g., pH calibration curve, ORP corrections) are included in this report. The measured fa of the simulants was close to the OLI model predictions for both brines and was largely controlled by the background H 2 partial pressure in the vapor phase as well as H 2 generated in situ in the aqueous phase by the Fe 0 corrosion. The H 2 gas-phase concentration tested and thermodynamically evaluated was much lower than is assumed in the WIPP PA; however, H 2 (g) concentrations significantly below this level are still predicted to result in very reducing conditions. In conclusion: • The experimental results are consistent with thermodynamic model predictions for fa, pH, and the effects of high ionic strength. • Evidence to date suggests that the H2 concentration in the glovebox atmosphere ultimately determined the final E h values of the simulants and resulted in highly reducing conditions. As a result, little difference was observed between the control simulants containing only MgO and the WIPP base-case simulants that contained MgO and Fe 0 . • This test methodology is recommended for future studies evaluating WIPP repository conditions. The methodology includes: (1) background H 2 in argon with agitation ( or could alternatively include in-situ-generated H 2 in sealed bottles); (2) carefully measured and corrected ORP data ( with much effort focused on allowing the probes to fully stabilize); and (3) ionic-strength-corrected pH data. Other best practices, such as simulant sparging/handling, ORP probe replacement, etc., should also be considered. • The coupling of experimental studies and thermodynamic modeling is also highly recommended because these methods inform and direct one another leading to greater confidence in and understanding of the results.

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

Integrated Disposal Facility FY2011 Glass Testing Summary Report [Erratum]

Pacific Northwest National Laboratory was contracted by Washington River Protection Solutions, LLC to provide the technical basis for estimating radionuclide release from the engineered portion of the disposal facility (e.g., source term). Vitrifying the low-activity waste at Hanford is expected to generate over 1.6 x 10 5 m 3 of glass (Certa and Wells 2010). The volume of immobilized low-activity waste (ILAW) at Hanford is the largest in the DOE complex and is one of the largest inventories (approximately 8.9 x 10 14 Bq total activity) of long-lived radionuclides, principally 99 Tc (t 1/2 = 2.1 x 10 5 ), planned for disposal in a low-level waste (LLW) facility. Before the ILAW can be disposed, DOE must conduct a performance assessment (PA) for the Integrated Disposal Facility (IDF) that describes the long-term impacts of the disposal facility on public health and environmental resources. As part of the ILAW glass testing program PNNL is implementing a strategy, consisting of experimentation and modeling, in order to provide the technical basis for estimating radionuclide release from the glass waste form in support of future IDF PAs. The purpose of this report is to summarize the progress made in fiscal year (FY) 2011 toward implementing the strategy with the goal of developing an understanding of the long-term corrosion behavior of low-activity waste glasses.

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

Solubility controls on plutonium and americium release in subsurface environments exposed to acidic processing wastes

To identify the role of waste composition and sediment interactions in controlling Pu and Am mobility in contaminated sediments at Hanford, a legacy nuclear site, Pu and Am concentrations in solutions equilibrated with contaminated sediments from beneath the 216-Z-9 (Z-9) Trench were compared to the solubilities of PuO 2 materials, synthesized by methods representative of the disposed wastes, in the absence of sediments. Furthermore, this work shows that the solubilities of PuO 2 materials synthesized by different methods, and with varying particle sizes, agree with PuO 2 (am,hyd), although dissolution kinetics differed between materials. According to saturation index (SI) calculations, PuO 2 (am,hyd) is likely also controlling Pu release from sediments under conditions where phosphate concentrations are low. However, both Pu-phosphate and Am-phosphate phases, identified in SI calculations and by high resolution transmission electron microscopy, play roles in controlling release in low pH, high phosphate, shallow sediments just below the Z-9 Trench. The elevated phosphate is likely due to decomposition of tributyl phosphate from waste solutions over time. Sediments from deeper in the subsurface beneath the Z-9 Trench are less acidic and contain less phosphate, with Pu solubility likely controlled by PuO 2 (am,hyd) that precipitated following neutralization of the acidic waste stream. Controls on Am concentrations in deeper sediments are more complex and potentially involve sediment adsorption and/or release from Pu (1-x) Am x O 2 following Am in-growth. The concentrations of both Pu and Am were elevated in the colloidal fraction associated with shallow sediments, but not in PuO 2 experiments, suggesting the presence of Pu/Am pseudocolloids (e.g., Pu/Am associated with mineral colloids). However, Pu and Am association with the colloidal size fraction was not observed in deeper sediments, suggesting transport of Pu and Am to these depths beneath the Z-9 Trench was not due to colloidal transport.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

Reactive Transport Modeling of Hydrogen Production from Serpentinization of Olivine

Hydrogen production from serpentinization of ultramafic rocks represents a promising natural pathway for generating carbon-free energy, yet its kinetics and controlling factors remain incompletely understood. A key challenge in advancing serpentinization research lies in the heterogeneity of porosity and permeability in rocks, which leads to nonuniform fluid velocity fields, as well as uncertainties in estimating reactive surface area and identifying appropriate mineral reaction equilibria. Additional complexities arise from the role of dissolved SiO 2 , Fe 2+ /Fe 3+ partitioning, and the limited effect of pH variations within the strongly alkaline regime on hydrogen yields. These challenges hinder straightforward extrapolation from laboratory tests to practical applications of hydrogen production from natural rocks. Here, in this work, we address these questions using a simulation-based reactive transport modeling framework calibrated against controlled laboratory experiments reported elsewhere. The model couples geochemical kinetics, multiphase flow, and mineralogical feedbacks, enabling systematic evaluation of how surface area, dissolved silica concentration, Fe redox state, temperature, and pressure govern serpentinization and H2 generation. We find that surface area exerts the strongest control on reaction rates and hydrogen yields, while Fe 2+ /Fe 3+ ratios act as secondary modulators. Elevated dissolved silica concentrations suppress hydrogen production but accelerate serpentine precipitation, whereas increasing pH beyond 12 within the strongly alkaline regime produces only marginal gains. Finally, we demonstrate that integrating targeted experiments with calibrated simulations offers a powerful and efficient approach for predicting hydrogen yields and assessing parameter trade-offs in industrial-scale applications. This integration can substantially reduce the experimental burden while improving predictive capability, thereby enhancing both the mechanistic understanding and the practical feasibility of hydrogen production from serpentinization.

08 HYDROGEN↗
DOE OSTI2021

Thin Water Films Enable Low-Temperature Magnesite Growth Under Conditions Relevant to Geologic Carbon Sequestration

Injecting supercritical CO2 (scCO2) into basalt formations for long-term storage is a promising strategy for mitigating CO2 emissions. Mineral carbonation can result in permanent entrapment of CO2; however, carbonation kinetics in thin H2O films in humidified scCO2 is not well understood. We investigated forsterite (Mg2SiO4) carbonation to magnesite (MgCO3) via amorphous magnesium carbonate (AMC; MgCO3·xH2O, 0.5 < x < 1), with the goal to establish the fundamental controls on magnesite growth rates at low H2O activity and temperature. Experiments were conducted at 25, 40, and 50 °C in 90 bar CO2 with a H2O film thickness on forsterite that averaged 1.78 ± 0.05 monolayers. In situ infrared spectroscopy was used to monitor forsterite dissolution and the growth of AMC, magnesite, and amorphous SiO2 as a function of time. Geochemical kinetic modeling showed that magnesite was supersaturated by two to three orders of magnitude and grew according to a zero-order rate law. The results indicate that the main drivers for magnesite growth are sustained high supersaturation coupled with low H2O activity, a combination of thermodynamic conditions not attainable in bulk aqueous solution. This improved understanding of reaction kinetics can inform subsurface reactive transport models for better predictions of CO2 fate and transport.

Kerisit, Sebastien N.↗
DOE OSTI2025

Geo-Responsive Chemomechanics in Aluminum Oxyhydroxide via Alkali-Driven Dehydroxylation for Supercritical Geothermal Systems

Widespread use of enhanced geothermal systems can revolutionize global renewable electrical power access, yet its advancement is hindered by the inherent instability of Portland cement-based chemistries for geothermal well construction under high-temperature corrosive conditions. Here, in this work, we demonstrate the tunable mechanical performance of aluminum oxyhydroxides as cementitious materials through an alkali-controlled dehydroxylation reaction pathway for long-term applications under supercritical geothermal environments. Notably, the synthesized aluminum oxyhydroxides demonstrate remarkable stability, maintaining superior mechanical performance under supercritical conditions for over 30 days. Synchrotron X-ray diffraction, spectroscopy measurements and geochemical thermodynamic modeling uncover that the gibbsite dehydroxylation pathway functions as a key dial for tuning the chemomechanics, rendering the aluminum oxyhydroxide a strong cementitious material. By uncovering the mechanistic role of alkali-driven dehydroxylation, this work proposes a cementitious chemistry distinct from conventional Portland cement and geopolymer-dominated alkali-activated systems, laying the groundwork for developing next-generation cementitious materials for supercritical geothermal energy exploitation.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

A Geo‐Structurally Based Correction Factor for Apparent Dissolution Rates in Fractured Media

Abstract Field measurements of apparent geochemical weathering reaction rates in subsurface fractured porous media are known to deviate from laboratory measurements by multiple orders of magnitude. To date, there is no geologically based explanation for this discrepancy that can be used to predict reaction rates in field systems. Proposed correction factors are typically based on ad hoc characterizations related to geochemical kinetic models. Through a series of high‐fidelity reactive transport simulations of mineral dissolution within explicit 3D discrete fracture networks, we are able to link the geo‐structural attributes with reactive transport observations. We develop a correction factor to linear transition state theory for the prediction of the apparent dissolution rate based on measurable geological properties. The modified rate law shows excellent agreement with numerical simulations, indicating that geological structure could be a primary reason for the discrepancy between laboratory and field observations of apparent dissolution rates in fractured media.

58 GEOSCIENCES↗
DOE OSTI2024

Weathering Incongruence in Mountainous Mediterranean Climates Recorded by Stream Lithium Isotope Ratios

Lithium isotope ratios (δ 7 Li) of rivers are increasingly serving as a diagnostic of the balance between chemical and physical weathering contributions to overall landscape denudation rates. Here, we show that intermediate weathering intensities and highly enriched stream δ 7 Li values typically associated with lowland floodplains can also describe small upland watersheds subject to cool, wet climates. This behavior is revealed by stream δ 7 Li between +22.4 and +23.5‰ within a Critical Zone observatory located in the Cévennes region of southern France, where dilute stream solute concentrations and significant atmospheric deposition otherwise mask evidence of incongruence. The water-rock reaction pathways underlying this behavior are quantified through a multicomponent, isotope-enabled reactive transport model. Using geochemical characterization of soil profiles, bedrock, and long-term stream samples as constraints, we evolve the simulation from an initially unweathered granite to a steady state weathering profile which reflects the balance between (a) fluid infiltration and drainage and (b) bedrock uplift and soil erosion. Enriched stream δ 7 Li occurs because Li is strongly incorporated into actively precipitating secondary clay phases beyond what prior laboratory experiments have suggested. Chemical weathering incongruence is maintained despite relatively slow reaction rates and moderate clay accumulation due to a combination of two factors. First, reactive primary mineral phases persist across the weathering profile and effectively “shield” the secondary clays from resolubilization due to their greater solubility. Second, the clays accumulating in the near-surface profile are relatively mature weathering byproducts. These factors promote characteristically low total dissolved solute export from the catchment despite significant input of exogenous dust.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1984

An integrated dynamical and geochemical approach to lunar origin modelling

The three major categories of models of lunar origin which explain the Moon's properties are complete and a more general scenario is presented. The model presented is as the Earth grew by planetesimal bombardment, a circumterrestrial cloud of particles was created from a combination of impact ejected mantle material and planetismals captured directly into orbit around the Earth. The compositional properties are explained two ways: (1) a few big late planetismals of diverse composition are captured in orbit and/or hit the Earth; and (2) the circumterrestrial swarm acts as a filter, preferentially capturing small weak silicate bodies, while passing large iron planetismals.

Greenberg, R.↗
NASA NTRS1985

Seafloor hydrothermal activity and spreading rates - The Eocene carbon dioxide greenhouse revisited

A suggestion has been made that enhanced rates of hydrothermal activity during the Eocene could have caused a global warming by adding calcium to the ocean and pumping CO2 into the atmosphere (Owen and Rea, 1984). This phenomenon was purported to be consistent with the predictions of the CO2 geochemical cycle model of Berner, Lasaga and Garrels (1983) (henceforth BLAG). In fact, however, the BLAG model predicts only a weak connection between hydrothermal activity and atmospheric CO2 levels. By contrast, it predicts a strong correlation between seafloor spreading rates and pCO2, since the release rate of CO2 from carbonate metamorphism is assumed to be proportional to the mean spreading rate. The Eocene warming can be conveniently explained if the BLAG model is extended by assuming that the rate of carbonate metamorphism is also proportional to the total length of the midocean ridges from which the spreading originates.

Kasting, J. F.↗