Reactive transport modeling involving precipitation and dissolution in porous and fractured media
Explore the source record for details and available documents.
Yoon, Hongkyu↗
SEARCH · Search NASA
Results for “precipitation / dissolution models”
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.
Pore‐Scale Modeling of Reactive Transport with Coupled Mineral Dissolution and Precipitation
Abstract We present a new pore‐scale model for multicomponent advective‐diffusive transport with coupled mineral dissolution and precipitation. Both dissolution and precipitation are captured simultaneously by introducing a phase transformation vector field representing the direction and magnitude of the overall phase change. An effective viscosity model is adopted in simulating fluid flow during mineral dissolution‐precipitation that can accurately capture the velocity field without introducing any empirical parameters. The proposed approach is validated against analytical solutions and interface tracking simulations in simplified structures. After validation, the proposed approach is employed in modeling realistic rocks where mineral dissolution and precipitation are dominant at different locations. We have identified three regimes for mineral dissolution‐precipitation coupling: (a) compact dissolution‐precipitation where dissolution is dominant near the inlet and precipitation is dominant near the outlet, (b) wormhole dissolution with clustered precipitation where dissolution generates wormholes in the main flow paths and precipitation clogs the secondary flow paths, and (c) dissolution dominant where all solid grains are gradually dissolved. In the three regimes, the proposed approach provides reliable porosity‐permeability relationships that cannot be described well by traditional macroscale models. We find that the permeability can increase while the overall porosity decreases when the main flow paths are expanded by dissolution and adjacent pore spaces are clogged by precipitation.
58 GEOSCIENCES↗
An Optimization-Based Law of Mass Action Precipitation/Dissolution Model
Rare earth elements (REE) and many other critical minerals are necessary for the manufacturing of modern everyday technologies, including microchips, batteries and electric motors. Recovery of these materials typically involves aqueous systems which can be modeled as chemical equilibrium problems. One common method for solving these problems is the law of mass action approach (LMA), where a system of non-linear equations involving the equilibrium constants is solved. However, despite being theoretically simple, these problems are in practice very difficult to solve. Currently, the use of iterative heuristics based on saturation indices to decide on which species and reactions to include in the calculations is the state of the art to arrive at a solution. Here, we present an optimization-based alternative to solve chemical equilibria problems involving precipitation/dissolution reactions without the need for such heuristics. Our approach is first validated against the LMA software MINTEQ and PHREEQC for a number of case studies, and then applied to a novel REE recovery process reported in the literature. Overall, our approach was found to have close agreement with MINTEQ and PHREEQC, and we were able to successfully replicate the reported yield and purity for the published REE process.
42 ENGINEERING↗
Real-Time Atomic-Scale Structural Analysis Resolves the Amorphous to Crystalline CaCO 3 Mechanism Controversy
Amorphous calcium carbonate (ACC) occurs as a precursor to geological and biogenic calcium carbonate (CaCO 3 ), yet its transformation pathways and reaction mechanisms remain inconsistent and controversial. In this study, we investigated the transformation of ACC to calcite under both solution and dry conditions, in the presence and absence of impurity ions, utilizing operando time-resolved synchrotron X-ray diffraction (TRXRD) and reactive transport modeling. Results demonstrate that TRXRD techniques allow us to differentiate dissolution-reprecipitation versus solid-state transformation mechanisms for amorphous to crystalline phase transitions. Specifically, we observe that in environments with abundant water, ACC transforms to calcite through a dissolution-reprecipitation mechanism. This features an activation energy of 63 ± 2 kJ/mol and unit cell volume contraction during calcite crystal growth. Conversely, under water-limited conditions, ACC to calcite transformation proceeds via a solid-state transformation mechanism, with an activation energy of 210 ± 2 kJ/mol, three times greater than the dissolution-reprecipitation route, and a unit cell expansion during crystalline calcite growth. Further, to illustrate the magnitude of these effects, the rates of calcite growth were similar during dissolution-reprecipitation at 3 °C [0.00207(35) s –1 ] and solid-state transformation at 280 °C [0.00134(11) s –1 ]. Moreover, the incorporation of an impurity, strontium, significantly retards the rate of calcite growth while expanding its unit cell but whose incorporation is history dependent. Reactive transport modeling of the dissolution–precipitation kinetics suggests that ACC must be dissolving as compact aggregates. These various transformation mechanisms drive diverse geological and biological carbonate formations, impacting their use as paleoenvironmental markers and functional materials synthesis.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Reactive-transport modeling of fracture flow to quantify the changes in flow pathways from matrix thermal contraction and mineral precipitation and dissolution; influence of grid resolution
A series of reactive-transport models of Enhanced Geothermal Systems (EGS) were constructed using the reactive transport code PFLOTRAN to examine the effect of matrix thermal contraction and mineral dissolution/precipitation on fracture flow in the context of grid cell size and model complexity. It was found that for thermal drawdown at production well, the impact of fracture zone grid cell size is negligible.
Gatz-Miller, Hannah Sarah↗
Dissolution of Nb(C,N) During Post-weld Heat-Treatment of Electric Resistance Welded X70 Line Pipe
The effect of peak temperature for post-weld heat treatment (PWHT) on the dissolution of Nb-rich carbonitrides, Nb(C,N), in the weldline of an electric resistance welded (ERW) X70 grade pipeline steel was investigated. Selected-area electron diffraction patterns and corresponding electron energy loss spectroscopy mappings revealed both Nb-rich carbonitrides and Ti-rich nitrides, (Ti,Nb)N, in the sample heat-treated to 1080 °C peak temperature. However, only Ti-rich nitrides were revealed in the sample heat-treated to 1220 °C peak temperature. During continuous heating to higher than the Ac3 temperature, i.e., near 1100 °C temperature, the dissolution of Nb(C,N) was observed in dilation and derivative of dilation. Notably, this experimental result agreed with the predictions by a DICTRA model. Furthermore, a phenomenological model was developed to explain the volume shrinkage during the dissolution of Nb(C,N) precipitates.
36 MATERIALS SCIENCE↗
MX precipitate behavior in an irradiated advanced Fe-9Cr steel: Helium effects on phase stability
As part of an ongoing series aimed at optimizing Fe-9Cr reduced activation ferritic/martensitic (RAFM) alloys for fusion energy systems, this study explores MX precipitate behavior under dual-ion irradiations, specifically examining correlations between helium transmutation and irradiation-induced damage. Utilizing single and dual-beam ion irradiation, the research explores the combined effects of helium (10–25 appm He/dpa), temperature (400–600 °C), and damage levels (15–100 dpa) on the microstructural evolution of CNA9 steel, a variant of Castable Nanostructured Alloys (CNAs). The study demonstrates that helium co-implantation hinders radiation-enhanced coarsening of MX-TiC precipitates at 500 and 600 °C, maintaining MX-TiC precipitate stability at moderate damage levels (15 dpa) but failing to prevent complete precipitate dissolution at higher damage levels (≥50 dpa) when irradiated at 500 °C. Here, a generalized precipitate stability model suggests that helium-induced suppression of diffusion alters the balance between recoil resolution and back diffusion for MX-TiC precipitates, enhancing the current understanding of precipitate behavior under damage and transmutation simulated dual-ion irradiation conditions.
Characterization↗
A solvent-targeted recovery and precipitation scheme for the recycling of up to ten polymers from post-industrial mixed plastic waste
Solvent-Targeted Recovery and Precipitation (STRAP TM ) separates polymers within a plastic waste stream by selective dissolution. In this work, the STRAP framework, which combines computational modeling and experiments, was applied to develop a series of steps to separate up to 10 polymers from post-industrial mixed plastic waste (MPW) and the main components recovered were LDPE, HDPE, and PET. The STRAP steps were initially demonstrated with a physical polymer mixture containing LDPE, HDPE, PS, PVC, EVOH, PET, PP, PA6, PA66, and PA66/6, in which recoveries of 89% or higher were achieved for each polymer. Furthermore, this paper demonstrates a solvent selection approach that can be applied to separate unknown plastics materials into purer components.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗
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↗
In-situ ion irradiation of fission products in a spent UO 2 fuel
This study investigates the behavior of fission gas bubbles and five metal precipitates (5MPs) (Mo, Ru, Rh, Tc, Pd) in spent uranium dioxide (UO 2 ) fuel under various ion irradiation doses, temperatures, and flux conditions. Utilizing in-situ ion irradiation and advanced transmission electron microscopy, we analyzed the evolution of fission gas bubbles and 5MPs in UO 2 samples from the Belgium Reactor 3 (BR-3). Our findings reveal significant shrinkage of fission gas bubbles and 5MPs with increasing irradiation dose, accompanied by a decrease in pair density. We demonstrate that ion irradiation induces a homogeneous re-solution process where individual atoms are ejected from bubbles and precipitates, leading to their dissolution and subsequent re-precipitation in the matrix. In conclusion, this study provides critical insights into the dynamic behavior of fission products under irradiation, facilitating the development of predictive models and contributing to the optimization of nuclear fuel performance and safety.
Fission products↗
In-situ ion irradiation of a spent UO2 fuel: evolution of fission products and nanograins with radiation dose
This study investigates the behavior of fission gas bubbles, five metal precipitates (5MPs) (Mo, Ru, Rh, Tc, Pd) and nanograins in spent uranium dioxide (UO2) fuel under various ion irradiation doses, temperatures, and flux conditions. Utilizing in-situ ion irradiation and advanced transmission electron microscopy, we analyzed the evolution of fission gas bubbles, 5MPs and nanograins in UO2 samples from the Belgium Reactor 3 (BR-3). Our findings reveal significant shrinkage of fission gas bubbles and 5MPs with increasing irradiation dose, accompanied by a decrease in pair density. The study also reveals nanograin coarsening under irradiation. We demonstrate that ion irradiation induces a homogeneous re-solution process where individual atoms are ejected from bubbles and precipitates, leading to their dissolution and subsequent re-precipitation in the matrix. This study provides critical insights into the dynamic behavior of fission products and nanograins under irradiation, facilitating the development of predictive models and contributing to the optimization of nuclear fuel performance and safety.
11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Computer simulation of reaction and transport of core-scale serpentinization of Fe-bearing olivine and evolution of geological hydrogen
We present reactive flow and transport simulations of core-scale serpentinization of forsteritic, Fe-bearing olivine under hydrothermal conditions. The model captures fluid flow, evolving porosity–permeability, and redox-controlled H 2 production in a cylindrical core over 55 days at 245 °C and 37 bar. Key processes include olivine dissolution, precipitation of serpentine and magnetite, and oxidation of Fe2+ to Fe 3+ . Here, results show strong coupling among flow velocity, alteration front propagation, and spatially heterogeneous H 2 generation. Elevated H2 concentrations align with inlet-localized magnetite precipitation, consistent with experiments. Secondary mineral formation reduces porosity and permeability, altering transport pathways. Reactive flow at the fluid–solid interface governs H 2 generation rates and distribution. Dissolved SiO2 promotes serpentine and talc formation by suppressing brucite, while bicarbonate extends reaction duration and moderates surface complexation effects.
Geologic hydrogen↗
Geochemical Impact of Acid Spearhead and Slickwater Stimulation on Wolfcamp Shale from the Hydraulic Fracturing Test Site
The Hydraulic Fracturing Test Site 1 (HFTS-1) was a field study performed in the Wolfcamp Formation in the West Texas Permian (Midland) Basin, USA, with a focus on improving the efficiency of hydraulic fracturing. Investigating site-specific rock-fluid geochemical interactions during hydraulic fracturing is an important step to understanding the impact on formation shale porosity, permeability, and long-term shale gas production. During field operations in this region, hydraulic fracturing fluid (HFF) injection usually starts with a concentrated acid spearhead for rapid rock dissolution, followed by the injection of near-neutral pH slickwater containing chemicals and proppants. A multistep sequential injection approach was used to investigate different stages of rock-fluid interactions. The carbonate content in the host rock is important when acid spearhead is considered, as carbonate mineral dissolution is rapid and can result in porosity and permeability changes in the shale matrix. Here, in this study, we designed flow-through experiments using fractured carbonate-rich and clay-rich Wolfcamp shale cores with (1) a short-time acid soaking step and (2) a long-term slickwater flow-through step to simulate the injection method used at HFTS-1. The fluid chemistry was analyzed. A thorough mineralogical progression [e.g., Calcium (Ca) dissolution and iron (Fe) redox progression] in the cores during HFF injection was also characterized and imaged by synchrotron microprobe. Reactive transport modeling was performed based on the experimental setup. The results showed that the acid spearhead is a crucial step in creating a reaction front by mineral dissolution, especially in carbonate-rich shales. A slight layer of ferrihydrite precipitated during the slickwater flow-through period. This study provides insights into potential geochemical impact due to hydraulic fracturing operations in the Permian Basin.
clastic rock↗
Geochemical Assessment for Carbon Sequestration in the Conasauga Group, Northwest Georgia, USA
Sedimentary geological formations are known to be great candidates for geological carbon sequestration. Published studies suggest the southeast of the United States contains many formations suitable for carbon storage. The Cassville 1 Stratigraphic Borehole well could act as a potential carbon reservoir for nearby energy resource facilities in Georgia, United States. Although studies have shown that porous formations are adequate for geological carbon sequestration, it is important to understand possible geochemical reactions between CO 2 and the targeted geological formation before injecting any fluids. In this study, a sandstone sample from the Cassville 1 well is being considered for geological carbon sequestration in the Conasauga Group in Northwest Georgia. Here, the collected sandstone sample, consisting of quartz, K-feldspar, micas, kaolinite, and carbonate minerals such as calcite and dolomite, has a 6% porosity. Leveraging the formation composition and porosity, a one-dimensional continuum reactive transport model was built using CrunchFlow to assess possible geochemical reactions between injected CO 2 and the geological formation. Simulation results show that the carbonate minerals, calcite and dolomite, dissolve during the injection period of 10,000 days, increasing formation porosity from 6% to as much as 30%. The rate and extent of carbonate mineral dissolution and resulting porosity increase are highly sensitive to mineral reactive surface area values. No evidence of mineral precipitation was observed, suggesting that dissolution reactions will control porosity evolution during the CO 2 injection period.
42 ENGINEERING↗
MX precipitate behavior in an irradiated advanced Fe-9Cr steel: Self-ion irradiation effects on phase stability
In an effort to optimize Fe-9Cr reduced activation ferritic/martensitic (RAFM) steels and to inform the design and operation of fusion reactors, this work represents the first in a series of cohesive studies dedicated to the evolution of MX-TiC precipitates under accelerated single and dual ion irradiations. This study investigates CNA9, a simplified Fe-9Cr RAFM steel featuring initial MX-TiC precipitate densities of (2.3±0.3)×10²¹ m⁻³. This material was subjected to single self-ion irradiation at damage levels ranging from 1 to 100 displacements per atom (dpa) over a temperature range of 300 to 600°C, with a nominal dose rate of 7×10⁻⁴ dpa/s. Irradiation-induced coarsening was observed, as evidenced by statistically significant increases in mean diameter sizes, at 15 dpa at both 500°C and 600°C, whereas no coarsening was noted at 300°C or 400°C. Further, complete dissolution of precipitates occurred at damage levels of 50 and 100 dpa across the two temperatures tested (300°C and 500°C) while no significant changes were observed at any doses below 15 dpa at 500°C. Experimentally parameterized recoil resolution modeling suggests that the observed radiation stability of MX-TiC precipitates is intricately linked to diffusional changes of solutes resulting from the co-evolution of microstructural features within the experiments. The findings align with current theoretical perspectives on radiation-induced precipitate stability in complex alloys.
36 MATERIALS SCIENCE↗
Reactive Transport and Mass Balance Modeling of the Stimson Sedimentary Formation and Altered Fracture Zones Constrain Diagenetic Conditions at Gale Crater, Mars
On a planet as cold and dry as present-day Mars, evidence of multiple aqueous episodes offers an intriguing view into very different past environments. Fluvial, lacustrine, and eolian depositional environments are being investigated by the Mars Science Laboratory Curiosityin Gale crater, Mars. Geochemical and mineralogical observations of these sedimentary rocks suggest diagenetic processes affected the sediments. Here, we analyze diagenesis of the Stimson formation eolian parent material, which caused loss of olivine and formation of magnetite. Additional, later alteration in fracture zones resulted in preferential dissolution of pyroxene and precipitation of secondary amorphous silica and Ca sulfate. The ability to compare the unaltered parent material with the reacted material allows constraints to be placed on the characteristics of the altering solutions. In this work we use a combination of a mass balance approach calculating the fraction of a mobile element lost or gained, τ, with fundamental geochemical kinetics and thermodynamics in the reactive transport code CrunchFlow to examine the characteristics of multiple stages of aqueous alteration at Gale crater, Mars. Our model results indicate that early diagenesis of the Stimson sedimentary formation is consistent with leaching of an eolian deposit by a near-neutral solution, and that formation of the altered fracture zones is consistent with a very acidic, high sulfate solution containing Ca, P and Si. These results indicate a range of past aqueous conditions occurring at Gale crater, Mars, with important implications for past Martian climate and environments.
Gale crater↗
Reactive Transport and Mass Balance Modeling of the Stimson Sedimentary Formation and Altered Fracture Zones Constrain Diagenetic Conditions at Gale Crater, Mars
On a planet as cold and dry as present-day Mars, evidence of multiple aqueous episodes offers an intriguing view into very different past environments. Fluvial, lacustrine, and eolian depositional environments are being investigated by the Mars Science Laboratory Curiosity in Gale crater, Mars. Geochemical and mineralogical observations of these sedimentary rocks suggest diagenetic processes affected the sediments. Here, we analyze diagenesis of the Stimson formation eolian parent material, which caused loss of olivine and formation of magnetite. Additional, later alteration in fracture zones resulted in preferential dissolution of pyroxene and precipitation of secondary amorphous silica and Ca sulfate. The ability to compare the unaltered parent material with the reacted material allows constraints to be placed on the characteristics of the altering solutions. In this work we use a combination of a mass balance approach calculating the fraction of a mobile element lost or gained, tau, with fundamental geochemical kinetics and thermodynamics in the reactive transport code CrunchFlow to examine the characteristics of multiple stages of aqueous alteration at Gale crater, Mars. Our model results indicate that early diagenesis of the Stimson sedimentary formation is consistent with leaching of an eolian deposit by a near-neutral solution, and that formation of the altered fracture zones is consistent with a very acidic, high sulfate solution containing Ca, P and Si. These results indicate a range of past aqueous conditions occurring at Gale crater, Mars, with important implications for past martian climate and environments.
Gale↗