Publications

The US Department of Energy Office of Nuclear Energy is conducting a brine availability heater test to characterize the thermal, mechanical, hydrological and chemical response of salt at elevated temperatures. In the heater test, brines will be collected and analyzed for chemical compositions. In order to support the geochemical modeling of chemical evolutions of the brines during the heater test, we are recalibrating and validating the solubility models for the mineral constituents in salt formations up to 100°C, based on the solubility data in multiple component systems as well as simple systems from literature. In this work, we systematically compare the model-predicted values based on the various solubility models related to the constituents of salt formations, with the experimental data. As halite is the dominant constituent in salt formations, we first test the halite solubility model in the Na-Mg-Cl dominated brines. We find the existing halite solubility model systematically over-predict the solubility of halite. We recalibrate the halite model, which can reproduce halite solubilities in Na-Mg-Cl dominated brines well. As gypsum/anhydrite in salt formations controls the sulfate concentrations in associated brines, we test the gypsum solubility model in NaCl solutions up to 5.87 mol•kg-1 from 25°C to 50°C. The testing shows that the current gypsum solubility model reproduces the experimental data well when NaCl concentrations are less than 1 mol•kg-1. However, at NaCl concentrations higher than 1, the model systematically overpredicts the solubility of gypsum. In the Na - Cl - SO4 - CO3 system, the validation tests up to 100°C demonstrate that the model excellently reproduces the experimental data for the solution compositions equilibrated with one single phase such as halite (NaCl) or thenardite (Na2SO4), with deviations equal to, or less than, 1.5 %. The model is much less ideal in reproducing the compositions in equilibrium with the assemblages of halite + thenardite, and of halite + thermonatrite (Na2CO3•H2O), with deviations up to 31 %. The high deviations from the experimental data for the multiple assemblages in this system at elevated temperatures may be attributed to the facts that the database has the Pitzer interaction parameters for Cl - CO3 and SO4 - CO3 only at 25°C. In the Na - Ca - SO4 - HCO3 system, the validation tests also demonstrate that the model reproduces the equilibrium compositions for one single phase such as gypsum better than the assemblages of more than one phase.

Abstract not provided.

Abstract not provided.

Aqueous dissolution of silicate glasses and minerals plays a critical role in global biogeochemical cycles and climate evolution. The reactivity of these materials is also important to numerous engineering applications including nuclear waste disposal. The dissolution process has long been considered to be controlled by a leached surface layer in which cations in the silicate framework are gradually leached out and replaced by protons from the solution. This view has recently been challenged by observations of extremely sharp corrosion fronts and oscillatory zonings in altered rims of the materials, suggesting that corrosion of these materials may proceed directly through congruent dissolution followed by secondary mineral precipitation. Here we show that complex silicate material dissolution behaviors can emerge from a simple positive feedback between dissolution-induced cation release and cation-enhanced dissolution kinetics. This self-accelerating mechanism enables a systematic prediction of the occurrence of sharp dissolution fronts (vs. leached surface layers), oscillatory dissolution behaviors and multiple stages of glass dissolution (in particular the alteration resumption at a late stage of a corrosion process). Our work provides a new perspective for predicting long-term silicate weathering rates in actual geochemical systems and developing durable silicate materials for various engineering applications.