Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geochemical Modeling”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Modeling Efforts to Gain Insight into Historical Leak Events from the Single-Shell Tank A-105 in the Hanford 241-A Tank Farm - 20113

The 241-A Tank Farm is a single-shell tank (SST) farm constructed to store process waste from Hanford nuclear operations. Millions of gallons of nuclear waste were stored in the 241-A Tank Farm SSTs and some of the SSTs leaked in the past. In addition, spills and pipeline leaks during transfers and storage and intentional discharges to cribs and trenches resulted in releasing waste to the ground. Tank liner leaks are referred to as 'leaks' and all other discharges to the soil are referred to as 'releases.' Liquid waste that could be removed by pumping has been removed from all of the SSTs to reduce the potential for future leaks. The most significant historical leak event in 241-A Tank Farm occurred in 1965 when tank A-105 experienced a pressurized steam event. This event resulted in damage to the inner steel liner flooring of A-105, which was separated from the sidewalls over part of its circumference and buckled up to 2.5 m, vertically. During this event up to 7,570 L of contaminated liquid may have leaked to the adjacent soil. This leak occurred under extreme temperature and pressure conditions. Modeling efforts to better understand the nature and extent of historical leaks from subsurface SSTs that were used to store highly radioactive, self-boiling liquid wastes at 241-A Tank Farm have recently been undertaken. These efforts have compiled the available historical information to support a preliminary non-isothermal, multiphase flow and chemical transport modeling effort to re-create the conditions under which the leaks occurred and to formulate a conceptual model as to the extent and distribution of leaked radioactive contaminants in the adjacent soil material. The focus of this work was on developing an understanding of the key features, processes and bounding conditions related to tank A-105 leak events that occurred in the 1960's. The activities include estimation of leak composition, non-isothermal multiphase flow and transport modeling, and geochemical modeling. The STOMP{sup C} Water-Air-Energy modeling code was utilized to implement the three-dimensional representation of the subsurface tank and the surrounding flow field. Once constructed, the model was used to evaluate the conceptual understanding of leaks originating from different parts of the tank under the elevated temperature boundary conditions that were imposed by surrounding tanks at the time of the historical leaks. The chemical evolution of the liquid waste was also evaluated as it leaked from the tank under transient pressure and temperature gradients. The results of the preliminary evaluation suggested the presence of a heat-pipe effect beneath tank A-105 in which water vapor at an elevated temperature is driven away from the base of the tank to a position where water vapor cools, condenses and is then drawn back toward the tank by the strong capillary attraction of the dry soil. A non-sorbing contaminant (Tc-99) was introduced into this flow field to better understand the potential distribution patterns of leaked contaminants that may have occurred during the historical leak events. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
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↗
DOE OSTI2022

ATS (Advanced Terrestrial Simulator) integrated hydrology and reactive transport model output in Copper Creek, Colorado

This dataset is generated using the ATS (Advanced Terrestrial Simulator) model at Copper Creek, Colorado, the largest catchment in the East River watershed. ATS is an integrated hydrology and reactive transport model to simulate the Concentration-Discharge (C-Q) relationship, and is used to quantify the geochemical export and understand how watersheds respond to climate disturbances. Key hydrogeochemical variables, including discharge, precipitation, evapotranspiration, and geochemical concentrations, were simulated and output at a daily frequency at the Copper Creek outlet. A group of hydrology and geochemical model variables over the spatial domain are output at a bi-monthly frequency. The model used 1-km resolution daily frequency Daymet meteorological forcing, and mineral composition from previous studies. The simulations are performed on NERSC. The datasets include the input files for running the ATS (Advanced Terrestrial Simulator) simulation, including the mesh file, input script and geochemical model files. The output files from the ATS (Advanced Terrestrial Simulator) are the daily time-series discharge and geochemical concentration simulated at the outlet of Copper Creek, and the bi-monthly spatial output over the simulation domain.

54 ENVIRONMENTAL 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↗
NASA NTRS2013

Does Aspartic Acid Racemization Constrain the Depth Limit of the Subsurface Biosphere?

Previous studies of the subsurface biosphere have deduced average cellular doubling times of hundreds to thousands of years based upon geochemical models. We have directly constrained the in situ average cellular protein turnover or doubling times for metabolically active micro-organisms based on cellular amino acid abundances, D/L values of cellular aspartic acid, and the in vivo aspartic acid racemization rate. Application of this method to planktonic microbial communities collected from deep fractures in South Africa yielded maximum cellular amino acid turnover times of approximately 89 years for 1 km depth and 27 C and 1-2 years for 3 km depth and 54 C. The latter turnover times are much shorter than previously estimated cellular turnover times based upon geochemical arguments. The aspartic acid racemization rate at higher temperatures yields cellular protein doubling times that are consistent with the survival times of hyperthermophilic strains and predicts that at temperatures of 85 C, cells must replace proteins every couple of days to maintain enzymatic activity. Such a high maintenance requirement may be the principal limit on the abundance of living micro-organisms in the deep, hot subsurface biosphere, as well as a potential limit on their activity. The measurement of the D/L of aspartic acid in biological samples is a potentially powerful tool for deep, fractured continental and oceanic crustal settings where geochemical models of carbon turnover times are poorly constrained. Experimental observations on the racemization rates of aspartic acid in living thermophiles and hyperthermophiles could test this hypothesis. The development of corrections for cell wall peptides and spores will be required, however, to improve the accuracy of these estimates for environmental samples.

cellular↗
DOE OSTI2024

UV-Vis spectrophotometric determination of rare earth elements (REE) speciation at near-neutral to alkaline pH. Part I: m-cresol purple properties from 25-75 °C and Er hydrolysis

The speciation and mobility of rare earth elements (REE) strongly depends on pH which controls the formation of charged aqueous hydroxyl species. The latter potentially play an important role in controlling heavy REE adsorption on clay minerals in near-neutral to alkaline waters such as in regolith-hosted REE mineral deposits. However, accurate REE hydrolysis constants are needed for developing geochemical models that can predict the role of these charged species in natural systems. Here, we develop a robust experimental UV-Vis spectrophotometric method using m-cresol purple to determine in situ pH from 25 to 75 °C. This method is used to derive the average ligand number and hydrolysis constants of erbium (Er) at 25 °C in aqueous solutions with low ionic strength (≤ 0.001 mol/L) at pH from ~7 to 9.5 and in the presence of Er concentrations from 0 to 0.057 mM. The average ligand number ranges between 1 and 3 indicating that Er(OH) 2+ , Er(OH) 2 + and Er(OH) 3 0 control speciation in the experiments. The logarithm of the Er hydrolysis constants (log*β n ° , n= 1 to 3) derived at infinite dilution for the reaction Er 3+ + nH 2 O = Er(OH) n 3-n + nH + are: *β 1 ° = –7.22 ± 0.10, *β 2 ° = –14.52 ± 0.08, *β 3 ° = –23.24 ± 0.04. Implementation of these experimental data into a geochemical model indicates that the Er(OH) 2 + and Er(OH) 3 0 species are both stable in a much wider pH range than previously predicted. Consequently, the positively charged REE hydroxyl complexes can potentially control the fractionation of light vs. heavy REE via adsorption as observed in the formation of certain regolith-hosted REE deposits.

58 GEOSCIENCES↗
DOE OSTI2024

Hydrothermal solubility of Dy hydroxide as a function of pH and stability of Dy hydroxyl aqueous complexes from 25 to 250 °C

The rare earth elements (REE) have important applications in green energy technologies. The formation of mineral deposits in geologic systems commonly involves hydrothermal fluids which can mobilize the REE. However, the REE speciation is not well known as a function of pH. The thermodynamic properties of REE hydroxyl complexes used in geochemical models are based on the Helgeson-Kirkham-Flowers (HKF) equation of state parameters which were derived by extrapolation of low temperature experimental and estimated data. In this study, Dy hydroxide solubility experiments are combined with available literature data to improve these models from 25 to 250 °C and optimize the thermodynamic properties of Dy 3+ and Dy hydroxyl complexes using GEMSFITS. Batch-type solubility experiments were conducted from 150 to 250 °C and at saturated water vapor pressure in perchloric acid solutions with initial pH values of 2 to 5 in 0.5 pH unit increments. The measured solubility of Dy hydroxide is retrograde with temperature and decreases with pH. The logarithm of total dissolved Dy molality ranges from –2.3 to –5.3 at 150 °C (pH 4.7–5.5), from –2.4 to –5.6 at 200 °C (pH 3.9–5.1), and from –3.7 to –6.9 at 250 °C (pH of 3.4 and 5.0). The optimized standard partial molal Gibbs energies of formation (Δ f G° T ) derived for Dy 3+ and DyOH 2+ display a close to linear relationship with temperature, fitting with previous optimizations based on DyPO 4 solubility data in the literature. A comparison of the optimized ΔfG°T values for aqueous Dy species with predictions from available HKF parameters indicates significant differences ranging from +11 to –26 kJ/mol between 25 and 250 °C. The experimental fits are used to derive the Dy hydroxide solubility products (K s0 ) and formation constants for the hydrolysis of Dy (β n with n = 1 to 3; Dy 3+ + nOH – = DyOH n 3-n ) as a function of temperature. The optimization method presented yields accurate thermodynamic properties for the Dy 3+ aqua ions and the DyOH 2+ species at the acidic to mildly acidic pH studied whereas more experimental work is needed at near-neutral and alkaline conditions to better constrain the other hydroxyl complexes. Furthermore, the optimized thermodynamic data have a significant impact on geochemical modeling of the mobility and solubility of REE minerals in acidic hydrothermal fluids.

58 GEOSCIENCES↗
NASA NTRS1987

Mutual interactions of redox couples via electron exchange in silicate melts - Models for geochemical melt systems

The mutual interactions via electron exchange of redox couples in glass-forming melts were investigated both theoretically and experimentally. A thermodynamic approach for considering the mutual interactions leads to conclusion that the degree of mutual interaction in the melt should be proportional in part to the difference in relative reduction potentials of the interacting redox couples. Experimental studies verify this conclusion for numerous redox couples in several composition/temperature/oxygen fugacity regimes. Geochemical systems simultaneously possess many potentially multivalent elements; the stabilized redox states in the resulting magmas can be explained in part by mutual interactions and by redox buffering through the central Fe(III)- Fe(II) couples in the melts. The significance of these results for basaltic magmas of the earth, moon, and meteorites is addressed.

Schreiber, Henry D.↗
DOE OSTI2021

Solubility Product of Vanadinite Pb 5 (VO 4 ) 3 Cl at 25 °C—A Comprehensive Approach to Incongruent Dissolution Modeling

Although vanadinite (Pb 5 (VO 4 ) 3 Cl occurs in abundance in various terrestrial geochemical systems of natural and anthropogenic origin and is seriously considered as a potential nuclear waste sequestering agent, its actual application is severely limited by a lack of understanding of its basic thermodynamic parameters. In this regard, the greatest challenge is posed by its incongruent dissolution, which is a pivotal hurdle for effective geochemical modeling. Our paper presents an universal approach for geochemical computing of systems undergoing incongruent dissolution which, along with unique, long-term experiments on vanadinites’ stability, allowed us to determine the mineral solubility constant. The dissolution experiments were carried out at pH = 3.5 for 12 years. Vanadinite has dissolved incongruently, continuously re-precipitating into chervetite (Pb 2 V 2 O 7 ) with the two minerals remaining in mutual equilibrium until termination of the experiments. The empirically derived solubility constant K sp,V,298 = 10 –91.89 ± 0.05 of vanadinite was determined for the first time. The proposed modeling method is versatile and can be adopted to other mineral systems undergoing incongruent dissolution.

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

Coupled Lattice Boltzmann Modeling Framework for Pore-Scale Fluid Flow and Reactive Transport

In this paper, we propose a modeling framework for pore-scale fluid flow and reactive transport based on a coupled lattice Boltzmann model (LBM). We develop a modeling interface to integrate the LBM modeling code parallel lattice Boltzmann solver and the PHREEQC reaction solver using multiple flow and reaction cell mapping schemes. The major advantage of the proposed workflow is the high modeling flexibility obtained by coupling the geochemical model with the LBM fluid flow model. Consequently, the model is capable of executing one or more complex reactions within desired cells while preserving the high data communication efficiency between the two codes. Meanwhile, the developed mapping mechanism enables the flow, diffusion, and reactions in complex pore-scale geometries. We validate the coupled code in a series of benchmark numerical experiments, including 2D single-phase Poiseuille flow and diffusion, 2D reactive transport with calcite dissolution, as well as surface complexation reactions. The simulation results show good agreement with analytical solutions, experimental data, and multiple other simulation codes. In addition, we design an AI-based optimization workflow and implement it on the surface complexation model to enable increased capacity of the coupled modeling framework. Compared to the manual tuning results proposed in the literature, our workflow demonstrates fast and reliable model optimization results without incorporating pre-existing domain knowledge.

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

Studies on the Encapsulation of Brine Generated from a Process Using Selective Electrodialysis Reversal

This topical report summarizes the laboratory testing, geochemical modeling, and mineralogical testing conducted to examine the efficacy of wastewater encapsulation for the brine generated from the Saltworks Flex electrodialysis reversal technology. Within the laboratory testing, the short- and long-term physical properties of seven separate mix designs were examined along with the leaching of each mix design. Additionally, geochemical thermodynamic modeling was used along with mineralogical analyses to help gain a fundamental understanding of the chemistry involved in the encapsulation process. Encapsulation is a process that results in the immobilization of constituents of interest by creating an engineered material that contains the waste within its matrix, both physically and chemically, in such a manner that the receiving environment is protected from adverse effects. Within the context of this study, encapsulation mix designs involved mixing the Saltworks brine with fly ash and various additives, principally lime (calcium oxide) and cement, to form a paste that eventually sets and hardens into a monolith upon placement in a landfill. Moisture conditioning fly ash with wastewater sources for on-site landfill disposal is a process that is in use at some sites and has been presented as a viable option to dispose of liquid wastes. However, because chemical reactions that can lead to sequestration of constituents of interest are not optimized, leaching can occur. Therefore, in this study, mixes were created to mimic moisture conditioning, as well as engineered pastes of varying mix designs. The study found that engineered encapsulation mixes provide significant resistance to leaching compared to moisture conditioning alone. The best performing encapsulation mix used lime as the additive, which proved superior to the use of cement. After 90 days of curing, this encapsulated mix was 5.5 times stronger than the moisture conditioned mix. Additionally, this mix exhibited an approximate 2.5 order of magnitude drop in hydraulic conductivity compared to the moisture conditioned mix. When comparing leaching, nearly all the chloride was found to evolve from the moisture conditioned mix, whereas only 6.7% of the chlorine contained in the best performing encapsulated mix leached out. Mineralogical analysis via X-ray diffraction (XRD) showed a significant concentration of the chloride sequestering mineral, Friedel’s salt, in the best performing encapsulation mix, whereas none was observed in the moisture conditioned mix, pointing to it being the primary chemical mechanism for chloride retention, along with the physical impermeability of the material. With respect to the modeling study, the model was successful in predicting the presence of most mineral phases observed in the XRD analysis.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2025

AquaMEND: Reconciling multiple impacts of salinization on soil carbon biogeochemistry

Soil salinization, exacerbated by climate change, poses a global threat to coastal ecosystem function and soil quality. Salinity influences carbon cycling through direct effects on microbial activity and indirect alterations to soil physicochemical properties including cation exchange, pH, and soil organic carbon availability. Current models inadequately represent these complexities, relying on linear reduction functions that overlook specific physicochemical changes induced by salinity. To address this gap, we propose an integrated model framework, AquaMEND, that combines microbial-explicit carbon decomposition and geochemical models. This model allows cation exchange and surface complexation processes to capture solute chemistry and nutrient availability in soils upon saltwater intrusion. Using response functions that capture salinity impacts on both salt-sensitive and salt-resistant microbial processes, AquaMEND simulates how the abiotic and biotic mechanisms work individually and collectively to regulate organic and inorganic pools and fluxes. Here, the parallel structure of aqueous and non-aqueous phases, together with microbial functions, result in a versatile model for solving dynamic coupling of organics, minerals and microbes under various environmental settings.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Assessment of chemo-mechanical impacts of CO 2 sequestration on the caprock formation in Farnsworth oil field, Texas

This study evaluates the chemo-mechanical influence of injected CO 2 on the Morrow B sandstone reservoir and the upper Morrow shale caprock utilizing data from the inverted 5-spot pattern centered on Well 13-10A within the Farnsworth unit (FWU). This study also seeks to evaluate the integrity of the caprock and the long-term CO 2 storage capability of the FWU. The inverted 5-spot pattern was extracted from the field-scale model and tuned with the available field observed data before the modeling work. Two coupled numerical simulation models were utilized to continue the study. First, a coupled hydro-geochemical model was constructed to simulate the dissolution and precipitation of formation minerals by modeling three intra-aqueous and six mineral reactions. In addition, a coupled hydro-geomechanical model was constructed and employed to study the effects of stress changes on the caprock’s porosity, permeability, and ground displacement. The Mohr–Coulomb circle and failure envelope were used to determine caprock failure. In this work, the CO 2 -WAG injection is followed by the historical field-observed strategy. During the forecasting period, a Water Alternating Gas (WAG) injection ratio of 1:3 was utilized with a baseline bottom-hole pressure constraint of 5500 psi for 20 years. A post-injection period of 1000 years was simulated to monitor the CO 2 plume and its effects on the CO 2 storage reservoir and caprock integrity. The simulation results indicated that the impacts of the geochemical reactions on the porosity of the caprock were insignificant as it experienced a decrease of about 0.0003% at the end of the 1000-year post-injection monitoring. On the other hand, the maximum stress-induced porosity change was about a 1.4% increase, resulting in about 4% in permeability change. It was estimated that about 3.3% of the sequestered CO 2 in the formation interacted with the caprock. Despite these petrophysical property alterations and CO 2 interactions in the caprock, the caprock still maintained its elastic properties and was determined to be far from its failure.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2025

Guiding Principles for Geochemical/Thermodynamic Model Development and Validation in Nuclear Waste Disposal: A Close Examination of Recent Thermodynamic Models for H + —Nd 3+ —NO 3 - (—Oxalate) Systems

Development of a defensible source-term model (STM), usually a thermodynamical model for radionuclide solubility calculations, is critical to a performance assessment (PA) of a geologic repository for nuclear waste disposal. Such a model is generally subjected to rigorous regulatory scrutiny. In this article, we highlight key guiding principles for STM model development and validation in nuclear waste management. We illustrate these principles by closely examining three recently developed thermodynamic models with the Pitzer formulism for aqueous H + —Nd 3+ —NO 3 - (—oxalate) systems in a reverse alphabetical order of the authors: the XW model developed by Xiong and Wang, the OWC model developed by Oakes et al., and the GLC model developed by Guignot et al., among which the XW model deals with trace activity coefficients for Nd(III), while the OWC and GLC models are for concentrated Nd(NO 3 ) 3 electrolyte solutions. The principles highlighted include the following: (1) Principle 1. Validation against independent experimental data: A model should be validated against experimental data or field observations that have not been used in the original model parameterization. We tested the XW model against multiple independent experimental data sets including electromotive force (EMF), solubility, water vapor, and water activity measurements. The results show that the XW model is accurate and valid for its intended use for predicting trace activity coefficients and therefore Nd solubility in repository environments. (2) Principle 2. Testing for relevant and sensitive variables: Solution pH is such a variable for an STM and easily acquirable. All three models are checked for their ability to predict pH conditions in Nd(NO 3 ) 3 electrolyte solutions. The OWC model fails to provide a reasonable estimate for solution pH conditions, thus casting serious doubt on its validity for a source-term calculation. In contrast, both the XW and GLC models predict close-to-neutral pH values, in agreement with experimental measurements. (3) Principle 3. Honoring physical constraints: Upon close examination, it is found that the Nd(III)-NO 3 association schema in the OWC model suffers from two shortcomings. Firstly, its second stepwise stability constant for Nd(NO 3 ) 2+ (log K 2 ) is much higher than the first stepwise stability constant for NdNO 3 2+ (log K 1 ), thus violating the general rule of (log K 2 –log K 1 ) < 0, or $\frac{K1}{K2}$>1. Secondly, the OWC model predicts abnormally high activity coefficients for Nd(NO 3 ) 2 + (up to ~900) as the concentration increases. (4) Principle 4. Minimizing degrees of freedom for model fitting: The OWC model with nine fitted parameters is compared with the GLC model with five fitted parameters, as both models apply to the concentrated region for Nd(NO 3 ) 3 electrolyte solutions. The latter appears superior to the former because the latter can fit osmotic coefficient data equally well with fewer model parameters. The work presented here thus illustrates the salient points of geochemical model development, selection, and validation in nuclear waste management.

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

Investigation of Coupled Processes in Fractures and the Bordering Matrix via a Micro‐Continuum Reactive Transport Model

Abstract In multi‐mineral fractured rocks, the altered porous layer on the fracture surface resulting from preferential dissolution of the fast‐reacting minerals can have profound impacts on subsequent chemical‐physical alteration of the fractures. This study adopts the micro‐continuum approach to provide further understanding of reactive transport processes in the altered layer (AL), and mass exchanges with the bordering matrix and fracture. The modeling framework couples the Darcy‐Brinkman‐Stokes (DBS) solver in COMSOL Multiphysics and the geochemical modeling capability of CrunchFlow. Three‐dimensional steady state simulations with systematically varied chemical‐physical parameters of the AL were performed to examine the impacts of individual factors and processes. Our simulation results confirm previous observations that dissolution of the fast‐reacting mineral (i.e., calcite) is largely controlled by diffusion across the AL. We also show that dissolution of the slow‐reacting mineral (i.e., dolomite), which controls AL development and fracture enlargement, increases with surface area and has a complex dependence on different local rate‐limiting processes. In particular, advection can result in evident spatial variations in the local dissolution rates of dolomite, although it does not affect the bulk chemistry significantly. The difference in the spatial patterns between simulations with and without advection in the AL is more noticeable in the locations with smaller apertures, with up to 20% difference in local reaction rates. Therefore, it is important to include a full depiction of advection, diffusion, and reactions for accurately capturing local dynamics that control long‐term fracture evolution.

58 GEOSCIENCES↗