Development of a defensible source-term model (STM), usual ly 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+—Nd3+—NO3−(—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(NO3)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(NO3)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)-NO3 association schema in the OWC model suffers from two shortcomings. Firstly, its second stepwise stability constant for Nd(NO3)2+ (log K2) is much higher than the first stepwise stability constant for NdNO32+ (log K1), thus violating the general rule of (log K2–log K1) < 0, or (Formula presented.). Secondly, the OWC model predicts abnormally high activity coefficients for Nd(NO3)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(NO3)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.
Oakes et al. (2023) published a review article in this journal. In that paper, Oakes et al. (2023) developed thermodynamic models to describe electrolyte solutions for HClO4–NaClO4–H2O and HBr–NaBr–H2O systems, based on literature data. In their paper, previously published work from researchers in the field was criticized; some of it is ours. Here, in this brief Comment, we first comment on their models, and then we briefly provide a technical response to that criticism.
Understanding pure H2 and H2/CH4 adsorption and diffusion in earth materials is one vital step toward a successful and safe H2 storage in depleted gas reservoirs. Despite recent research efforts such understanding is far from complete. In this work we first use Nuclear Magnetic Resonance (NMR) experiments to study the NMR response of injected H2 into Duvernay shale and Berea sandstone samples, representing materials in confining and storage zones. Then we use molecular simulations to investigate H2/CH4 competitive adsorption and diffusion in kerogen, a common component of shale. Our results indicate that in shale there are two H2 populations, i.e., free H2 and adsorbed H2, that yield very distinct NMR responses. However, only free gas presents in sandstone that yields a H2 NMR response similar to that of bulk H2. About 10 % of injected H2 can be lost due to adsorption/desorption hysteresis in shale, and no H2 loss (no hysteresis) is observed in sandstone. Our molecular simulation results support our NMR results that there are two H2 populations in nanoporous materials (kerogen). The simulation results also indicate that CH4 outcompetes H2 in adsorption onto kerogen, due to stronger CH4-kerogen interactions than H2-kerogen interactions. Nevertheless, in a depleted gas reservoir with low CH4 gas pressure, about ∼30 % of residual CH4 can be desorbed upon H2 injection. The simulation results also predict that H2 diffusion in porous kerogen is about one order of magnitude higher than that of CH4 and CO2. This work provides an understanding of H2/CH4 behaviors in deleted gas reservoirs upon H2 injection and predictions of H2 loss and CH4 desorption in H2 storage.
The objective of the crystalline disposal work packages is to advance our understanding of long-term disposal of used fuel in crystalline rocks and to develop necessary experimental and computational capabilities to evaluate various disposal concepts in such media.
Machine learning methodologies can provide insight into Brønsted-Guggenheim-Scatchard specific ion interaction theory (SIT) parameter values where experimental data availability may be limited. This study develops and executes machine learning frameworks to model the SIT interaction coefficient, ε. Key findings include successful estimations of ε via artificial neural networks using clustering and value prediction approaches. Applicability to other chemical parameters is also assessed briefly. Models developed here provide support for a use-case of machine learning in geologic nuclear waste disposal research applications, namely in predictions of chemical behaviors of high ionic strength solutions (i.e., subsurface brines).
Gas intercalation into clay interlayers may result in hydrogen loss in the geological storage of hydrogen; a phenomenon that has not been fully understood and quantified. Here we use metadynamics molecular simulations to calculate the free energy landscape of H2 intercalation into montmorillonite interlayers and the H2 solubility in the confined water; in comparison with results obtained for CO2. The results indicate that H2 intercalation into hydrated interlayers is thermodynamically unfavorable while CO2 intercalation can be favorable. H2 solubility in hydrated clay interlayers is in the same order of magnitude as that in bulk water and therefore no over-solubility effect due to nanoconfinement is observed - in striking contrast with CO2. These results indicate that H2 loss and leakage through hydrated interlayers due to intercalation in a subsurface storage system, if any, is limited.
Neodymium (Nd), a rare earth element (REE), is critical to numerous industries. Neodymium can be extracted from ore concentrates, waste materials, or recycled materials such as recycled Nd-Fe-B permanent magnets. In a standard process, concentrated sulfuric acid (H2SO4) is used as an extraction/leaching agent. Therefore, knowledge of Nd(III)–sulfate interaction at high ionic strengths is important for optimization of the extraction process. In addition, sulfate is also a major species in natural surface waters and present in nuclear waste streams. Nd(III) has been used a chemical analog to trivalent actinides in nuclear waste research and development. Consequently, knowledge of Nd(III)-sulfate interactions is also impactful to the field of nuclear waste management. In this study, we have developed a thermodynamic model that can describe the interaction of Nd(III) with sulfate to ionic strengths up to ~ 16.5 mol·kg–1 and to temperatures up to 100 °C. The model adopts the Pitzer formulation to describe activity coefficients of aqueous species. This model can be used to design and optimize a chemical process for REE recovery from ore concentrates, recycled materials, and acid mine drainage (AMD) and to understand the mobility of REEs and actinides in the environment.
Numerous experimental investigations indicated that expansive clays such as montmorillonite can intercalate CO2 preferentially into their interlayers and therefore potentially act as a material for CO2 separation, capture, and storage. However, an understanding of the energy-structure relationship during the intercalation of CO2 into clay interlayers remains elusive. Here, we use metadynamics molecular dynamics simulations to elucidate the energy landscape associated with CO2 intercalation. Our free energy calculations indicate that CO2 favorably partitions into nanoconfined water in clay interlayers from a gas phase, leading to an increase in the CO2/H2O ratio in clay interlayers as compared to that in bulk water. CO2 molecules prefer to be located at the centers of charge-neutral hydrophobic siloxane rings, whereas interlayer spaces close to structural charges tend to avoid CO2 intercalation. The structural charge distribution significantly affects the amount of CO2 intercalated in the interlayers. These results provide a mechanistic understanding of CO2 intercalation in clays for CO2 separation, capture, and storage.
This research presents a simple method to additively manufacture Cone 5 porcelain clay ceramics by using the direct ink-write (DIW) printing technique. DIW has allowed the application of extruding highly viscous ceramic materials with relatively high-quality and good mechanical properties, which additionally allows a freedom of design and the capability of manufacturing complex geometrical shapes. Clay particles were mixed with deionized (DI) water at different ratios, where the most suitable composition for 3D printing was observed at a 1:5 w/c ratio (16.2 wt.%. of DI water). Differential geometrical designs were printed to demonstrate the printing capabilities of the paste. In addition, a clay structure was fabricated with an embedded wireless temperature and relative humidity (RH) sensor during the 3D printing process. The embedded sensor read up to 65% RH and temperatures of up to 85 °F from a maximum distance of 141.7 m. The structural integrity of the selected 3D printed geometries was confirmed through the compressive strength of fired and non-fired clay samples, with strengths of 70 MPa and 90 MPa, respectively. This research demonstrates the feasibility of using the DIW printing of porcelain clay with embedded sensors, with fully functional temperature- and humidity-sensing capabilities.
Understanding the formation of H2CO3 in water from CO2 is important in environmental and industrial processes. Although numerous investigations have studied this reaction, the conversion of CO2 to H2CO3 in nanopores, and how it differs from that in bulk water, has not been understood. We use ReaxFF metadynamics molecular simulations to demonstrate striking differences in the free energy of CO2 conversion to H2CO3 in bulk and nanoconfined aqueous environments. We find that nanoconfinement not only reduces the energy barrier but also reverses the reaction from endothermic in bulk water to exothermic in nanoconfined water. Also, charged intermediates are observed more often under nanoconfinement than in bulk water. Stronger solvation and more favorable proton transfer with increasing nanoconfinement enhance the thermodynamics and kinetics of the reaction. Here our results provide a detailed mechanistic understanding of an important step in the carbonation process, which depends intricately on confinement, surface chemistry, and CO2 concentration.
Clays are known for their small particle sizes and complex layer stacking. We show here that the limited dimension of clay particles arises from the lack of long-range order in low-dimensional systems. Because of its weak interlayer interaction, a clay mineral can be treated as two separate low-dimensional systems: a 2D system for individual phyllosilicate layers and a quasi-1D system for layer stacking. The layer stacking or ordering in an interstratified clay can be described by a 1D Ising model while the limited extension of individual phyllosilicate layers can be related to a 2D Berezinskii–Kosterlitz–Thouless transition. This treatment allows for a systematic prediction of clay particle size distributions and layer stacking as controlled by the physical and chemical conditions for mineral growth and transformation. Clay minerals provide a useful model system for studying a transition from a 1D to 3D system in crystal growth and for a nanoscale structural manipulation of a general type of layered materials.
Heavy metals released from kerogen to produced water during oil/gas extraction have caused major enviromental concerns. To curtail water usage and production in an operation and to use the same process for carbon sequestration, supercritical CO2 (scCO2) has been suggested as a fracking fluid or an oil/gas recovery agent. It has been shown previously that injection of scCO2 into a reservoir may cause several chemical and physical changes to the reservoir properties including pore surface wettability, gas sorption capacity, and transport properties. Using molecular dynamics simulations, we here demonstrate that injection of scCO2 might lead to desorption of physically adsorbed metals from kerogen structures. This process on one hand may impact the quality of produced water. On the other hand, it may enhance metal recovery if this process is used for in-situ extraction of critical metals from shale or other organic carbon-rich formations such as coal.
Structural health monitoring of an engineered component in a harsh environment is critical for multiple DOE missions including nuclear fuel cycle, subsurface energy production/storage, and energy conversion. Supported by a seeding Laboratory Directed Research & Development (LDRD) project, we have explored a new concept for structural health monitoring by introducing a self-sensing capability into structural components. The concept is based on two recent technological advances: metamaterials and additive manufacturing. A self-sensing capability can be engineered by embedding a metastructure, for example, a sheet of electromagnetic resonators, either metallic or dielectric, into a material component. This embedment can now be realized using 3-D printing. The precise geometry of the embedded metastructure determines how the material interacts with an incident electromagnetic wave. Any change in the structure of the material (e.g., straining, degradation, etc.) would inevitably perturbate the embedded metastructures or metasurface array and therefore alter the electromagnetic response of the material, thus resulting in a frequency shift of a reflection spectrum that can be detected passively and remotely. This new sensing approach eliminates complicated environmental shielding, in-situ power supply, and wire routing that are generally required by the existing active-circuit-based sensors. The work documented in this report has preliminarily demonstrated the feasibility of the proposed concept. The work has established the needed simulation tools and experimental capabilities for future studies.
To mitigate adverse effects from molten corium following a reactor pressure vessel failure (RPVF), some new reactor designs employ a core catcher and a sacrificial material (SM), such as ceramic or concrete, to stabilize the molten corium and avoid containment breach. Existing reactors cannot easily be modified to include these SMs but could be modified to allow injectable cooling materials. Current reactor designs are limited to using water to stabilize the corium, but this can create other issues such as reaction of water with the concrete forming hydrogen gas. The novel SM proposed here is a granular carbonate mineral that can be used in existing light water reactor plants. The granular carbonate will decompose when exposed to heat, inducing an endothermic reaction to quickly solidify the corium in place and producing a mineral oxide and carbon dioxide. Corium spreading is a complex process strongly influenced by coupled chemical reactions, including decay heat from the corium, phase change, and reactions between the concrete containment and available water. A recently completed Sandia National Laboratories laboratory directed research and development (LDRD) project focused on two research areas: experiments to demonstrate the feasibility of the novel SM concept, and modeling activities to determine the potential applications of the concept to actual nuclear plants. Small-scale experiments using lead oxide (PbO) as a surrogate for molten corium demonstrate that the reaction of the SM with molten PbO results in a fast solidification of the melt due to the endothermic carbonate decomposition reaction and the formation of open pore structures in the solidified PbO from CO2 released during the decomposition. A simplified carbonate decomposition model was developed to predict thermal decomposition of carbonate mineral in contact with corium. This model was incorporated into MELCOR, a severe accident nuclear reactor code. A full-plant MELCOR simulation suggests that by the introduction of SM to the reactor cavity prior to RPVF ex-vessel accident progression, e.g., core-concrete interaction and core spreading on the containment floor, could be delayed by at least 15 h; this may be enough for additional accident management to be implemented to alleviate the situation.
Swelling clay hydration/dehydration is important to many environmental and industrial processes. Experimental studies usually probe equilibrium hydration states in an averaged manner and thus cannot capture the fast water transport and structural change in interlayers during hydration/dehydration. Using molecular simulations and thermogravimetric analyses, we observe a two-stage dehydration process. The first stage is controlled by evaporation at the edges: water molecules near hydrophobic sites and the first few water molecules of the hydration shell of cations move fast to particle edges for evaporation. The second stage is controlled by slow desorption of the last 1-2 water molecules from the cations and slow transport through the interlayers. The two-stage dehydration is strongly coupled with interlayer collapse and the coordination number changes of cations, all of which depend on layer charge distribution. This mechanistic interpretation of clay dehydration can be key to the coupled chemomechanical behavior in natural/engineered barriers.
Upon extraction/injection of a large quantity of gas from/into a subsurface system in shale gas production or carbon sequestration, the gas pressure varies remarkably, which may significantly change the wettability of porous media involved. Mechanistic understanding of such changes is critical for designing and optimizing a related subsurface engineering process. Using molecular dynamics simulations, we have calculated the contact angle of a water droplet on various solid surfaces (kerogen, pyrophyllite, calcite, gibbsite, and montmorillonite) as a function of CO2 or CH4 gas pressure up to 200 atm at a temperature of 300 K. The calculation reveals a complex behavior of surface wettability alteration by gas pressure variation depending on surface chemistry and structure, and molecular interactions of fluid molecules with surfaces. As the CO2 gas pressure increases, a partially hydrophilic kerogen surface becomes highly hydrophobic, while a calcite surface becomes more hydrophilic. Considering kerogen and calcite being the major components of a shale formation, we postulate that the wettability alteration of a solid surface induced by a gas pressure change may play an important role in fluid flows in shale gas production and geological carbon sequestration.
Actinide oxalates are chemical compounds important to nuclear industry, ranging from actinide separation in waste reprocessing, to production of specialty actinides, and to disposal of high level nuclear waste (HLW) and spent nuclear fuel (SNF). In this study, the solubility constants for Pr2(C2O4)3·10H2O and Nd2(C2O4)3·10H2O by performing solubility experiments in HNO3 and mixtures of HNO3 and H2C2O4 at 23.0 ± 0.2 °C have been determined. The targeted starting materials, Pr2(C2O4)3·10H2O and Nd2(C2O4)3·10H2O, were successfully synthesized at room temperature using PrCl3, NdCl3 and oxalic acid as the source metrials. Then, we utilized the targeted solubility-controlling phases to conduct solubility measurements. There was no phase change over the entire periods of experiments, demonstrating that Pr2(C2O4)3·10H2O and Nd2(C2O4)3·10H2O were the solubility-controlling phases in our respective experiments. Based on our experimental data, we have developed a thermodynamic model for Pr2(C2O4)3·10H2O and Nd2(C2O4)3·10H2O in the mixtures of HNO3 and H2C2O4 to high ionic strengths. The model for Pr2(C2O4)3·10H2O reproduces well the reported experimental data for Pu2(C2O4)3·10H2O, which are not utilized for the model development, demonstrating that Pr(III) is an excellent analog for Pu(III). Similarly, the model for Nd2(C2O4)3·10H2O reproduces the solubility of Am2(C2O4)3·10H2O and Cm2(C2O4)3·10H2O. The Pitzer model was used for the calculation of activity coefficients. Based on the published, well established model for dissociation constants for oxalic acid and stability constants for actinide-oxalate complexes [i.e., AmC2O4+, and Am(C2O4)2−] to high ionic strengths, we have obtained the solubility constants (log10K0) for the following reactions at 25 °C,Pr2(C2O4)3·10H2O ⇌ 2Pr3+ + 3C2O42− + 10H2O(l)Nd2(C2O4)3·10H2O ⇌ 2Nd3+ + 3C2O42− + 10H2O(l) to be −30.82 ± 0.30 (2σ), and −31.14 ± 0.35 (2σ), respectively. These values can be directly applied to Pu2(C2O4)3·10H2O, Am2(C2O4)3·10H2O and Cm2(C2O4)3·10H2O. The model established for actinide oxalates by this study provides the needed knowledge with regard to solubilities of actinide/REE oxalates at various ionic strengths, and is expected to find applications in many fields, including the geological disposal of nuclear waste and the mobility of REE under the surface conditions, as Pr2(C2O4)3·10H2O and Nd2(C2O4)3·10H2O can be regarded as the pure Pr and Nd end-members of deveroite, a recently discovered natural REE oxalate with the following stoichiometry, (Ce1.01Nd0.33La0.32Pr0.11Y0.11Sm0.01Pb0.04U0.03Th0.01Ca0.04)2.01(C2O4)2.99·9.99H2O. Regarding its importance in the geological disposal of nuclear waste, Am2(C2O4)3·10H2O/Pu2(C2O4)3·10H2O/Cm2(C2O4)3·10H2O can be the source-term phase for actinides, as demonstrated by the instance in the disposal in clay/shale formations. This is exemplified by the stability of Am2(C2O4)3·10H2O in comparison with Am(OH)3(am), Am(OH)3(s) and AmCO3(OH)(s) under the relevant geological repository conditions.