The swelling of clay at high temperature and pressure is important for applications including nuclear waste storage but is not well understood. A molecular dynamics study of the swelling of Na montmorillonite in water at several temperatures (T = 298, 400, and 500 K) and water environment pressures (Pe = 5 and 100 MPa) is reported here. Adopting a rarely used setup that enables swelling pressure to be resolved with an accuracy of ~1 MPa, the swelling pressure was computed at basal spacings 1.6–2.6 nm (or 2–5 water layers between neighboring clay sheets), which has not been widely studied before. At T = 298 K and Pe = 5 MPa, swelling pressure Ps oscillates at d-spacing d smaller than 2.2 nm and decays monotonically as d increases. Increasing T to 500 K but keeping Pe at 5 MPa, Ps remains oscillatory at small d, but its repulsive peak at d = 2.2 nm shifts to ~2.0 nm and Ps at different d-spacings can grow more attractive or repulsive. At d > 2.0 nm, Ps is weakened greatly. Keeping T at 500 K and increasing Pe to 100 MPa, Ps recovers toward that at T = 298 K and Pe = 5 MPa, however, the repulsive peak at d = 2.0 nm remains the same. The opposite effects of increasing temperature and pressure on the density and dielectric screening of water, which control ion correlations and thus double layer repulsion, are essential for understanding the observed swelling pressure at elevated temperatures and its response to environment pressures.
Water flow in nanometer or sub-nanometer hydrophilic channels bears special importance in diverse fields of science and engineering. However, the nature of such water flow remains elusive. Here, we report our molecular-modeling results on water flow in a sub-nanometer clay interlayer between two montmorillonite layers. We show that a fast advective flow can be induced by evaporation at one end of the interlayer channel, that is, a large suction pressure created by evaporation (∼818 MPa) is able to drive the fast water flow through the channel (∼0.88 m/s for a 46 Å-long channel). Scaled up for the pressure gradient to a 2 μm particle, the velocity of water is estimated to be about 95 μm/s, indicating that water can quickly flow through a μm-sized clay particle within seconds. The prediction seems to be confirmed by our thermogravimetric analysis of bentonite hydration and dehydration processes, which indicates that water transport at the early stage of the dehydration is a fast advective process, followed by a slow diffusion process. The possible occurrence of a fast advective water flow in clay interlayers prompts us to reassess water transport in a broad set of natural and engineered systems such as clay swelling/shrinking, moisture transport in soils, water uptake by plants, water imbibition/release in unconventional hydrocarbon reservoirs, and cap rock integrity of supercritical CO2 storage.
Due to a stress redistribution after the excavation of an underground tunnel for radioactive waste disposal, an Ed/DZ (excavation disturbed/damaged zone) will be generated in the near field of the opening, resulting in significant changes in the hydraulic and mechanical properties of the rock mass in the zone. Initially more or less randomly distributed hydrocarbons at grain boundaries in rock salt, which sometimes can only be observed with ultraviolet light, can then be mobilised and migrate at a potentially significant rate towards the excavation. Within the international cooperative project DECOVALEX 2019, the migration mechanism of such fluid inclusions in rock salt is being studied intensively. A multi-scale modelling strategy has been developed. A macroscale coupled hydro-mechanical modelling of an underground excavation was performed to determine hydraulic and time-dependent deviatoric stress conditions, by taking into account the rock salt creep behaviour. Under the obtained macro-scale constraints, micro-scale modelling of a pathway dilation along halite grain boundaries was performed using different model strategies: a) coupled hydromechanical modelling with a consideration of hydraulic pressure-induced dilatant deformation, b) nonlinear dynamic model taking account of fluid migration, stress-dependent grain boundary wetting and shear-induced dilatancy of salt, and c) phase-field modelling of flow pathway propagation. The permeability increase resulting from the pathway dilation is estimated to be as high as two orders of magnitude. Based on the permeability determined, a series of pressure build-ups measured from a borehole with a high hydrocarbon release rate, a total of 430 build-ups within a monitoring time of 938 days, can be simulated with a macro-scale compressible flow model accounting for different zones around the opening.
Disposal of large, heat-generating waste packages containing the equivalent of 21 pressurized water reactor (PWR) assemblies or more is among the disposal concepts under investigation for a future repository for spent nuclear fuel (SNF) in the United States. Without a long (>200 years) surface storage period, disposal of 21-PWR or larger waste packages (especially if they contain high-burnup fuel) would result in in-drift and near-field temperatures considerably higher than considered in previous generic reference cases that assume either 4-PWR or 12-PWR waste packages (Jové Colón et al. 2014; Mariner et al. 2015; 2017). Sevougian et al. (2019c) identified high-temperature process understanding as a key research and development (R&D) area for the Spent Fuel and Waste Science and Technology (SFWST) Campaign. A two-day workshop in February 2020 brought together campaign scientists with expertise in geology, geochemistry, geomechanics, engineered barriers, waste forms, and corrosion processes to begin integrated development of a high-temperature reference case for disposal of SNF in a mined repository in a shale host rock. Building on the progress made in the workshop, the study team further explored the concepts and processes needed to form the basis for a high-temperature shale repository reference case. The results are described in this report and summarized..
In shale gas production, gas composition may vary over time. To understand this phenomenon, we use molecular dynamics simulations to study the permeation of CH4, C2H6 and their mixture from a source container through a pyrophyllite nanopore driven by a pressure gradient. For a pure gas, the flow rate of CH4 is always higher than that of C2H6, regardless of pore size. For a 1:1 C2H6: CH4 mixture, however, C2H6:CH4 flow rate ratio is higher than the compositional ratio in the container (i.e., 1:1) when the pore size is smaller than ~1.8 nm. The selective transport is caused by the competitive adsorption of C2H6 over CH4 in the nanopore. The selectivity is also determined by the interplay between the surface diffusion of the adsorbed molecules and the viscous flow in the center of the pore, and it diminishes as the viscous flow becomes to dominate the surface diffusion when the pore size becomes larger than 1.8 nm. Our work shows that compositional differentiation of shale gas in production is a consequence of nanopore confinement and therefore a key characteristic of an unconventional reservoir. The related compositional information can potentially be used for monitoring the status of a production well such as its recovery rate.
In high-level radioactive waste disposal, a heat-generating waste canister is generally encased with a layer of bentonite-based buffer material acting as an engineered barrier to limit water percolation and radionuclide release. The low thermal conductivity of bentonite (~0.5 W/m∙K) combined with a high thermal loading waste package may result in a high surface temperature on the package that can potentially impact the structural integrity of the package itself as well as the surrounding buffer material. We show here that the thermal conductivity of bentonite can be effectively enhanced by embedding copper wires/meshes across the buffer layer to form fully connected high heat conduction pathways. A simple calculation based on Rayleigh’s model shows that a required thermal conductivity of 5 W/m∙K for effective heat dissipation can be achieved simply by adding ~1 vol % of copper wires/meshes into bentonite. As a result, the peak surface temperature on a large waste package such as a dual-purpose canister can be reduced by up to 300°C, leading to a significant reduction in the surface storage time for waste cooling and therefore the overall cost for direct disposal of such waste packages. Because of the ensured full thermal percolation across the buffer layer, copper wires/meshes turn out to be much more effective than any other materials currently suggested (such as graphene or graphite) in enhancing the thermal conductivity of buffer material. Furthermore, the embedded copper wires/meshes can help reinforce the mechanical strength of the buffer material, thus preventing the material from a potential erosion by a possible intrusion of dilute groundwater.
The U.S. Department of Energy Office of Spent Fuel Waste Disposition (SFWD) established in fiscal year 2010 (FY10) the Spent Fuel Waste Science & Technology (SFWST) Program (formerly the Used Fuel Disposition Campaign - UFDC) program to conduct the research and development (R&D) activities related to storage, transportation and disposal of used nuclear fuel and high level nuclear waste. The Mission of the SFWST is: To identify alternatives and conduct scientific research and technology development to enable storage, transportation and disposal of used nuclear fuel and wastes generated by existing and future nuclear fuel cycles. Significant progress has been made in FY20 in both experimental and modeling arenas in evaluation of used fuel disposal in crystalline rocks, especially in model demonstration using field data. The work covers a wide range of research topics identified in the R&D plan.
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.