Abstract:Clays and clay rocks are considered good natural and engineered barriers for deep geological disposal of nuclear waste worldwide. Metal corrosion and organic waste degradation in underground repositories generate significant amounts of gas that should be able to migrate through the multibarrier system to avoid potential pressure buildup, which could be compromising the integrity of the barriers and host rocks. The gas is expected to accumulate in larger pores and eventually form an interconnected network. Under such conditions, the migration of gas molecules takes place both in pore water films and gas-filled macropores. Therefore, mass fluxes depend on the distribution of gas molecules between the water-rich and gas-rich phases and their mobility in both compartments. Classical molecular dynamics (MD) simulations were employed to investigate the mobilities of He, H2, CO2, Ar, and CH4 in a Na-montmorillonite mesopore as a function of the degree of saturation, as well as evaluate the hydrodynamic behavior of the pore fluid in partially saturated clays. The diffusivity of the gas molecules was determined by observing the asymptotic behavior of the mean square displacement in the gas-rich phase and at the gas-water interface. The partition coefficient and Gibbs free energy were analyzed to investigate the transfer of gas molecules between the gas-rich and water-rich phases by observing the molecular trajectories as they cross the vapor-liquid interface. The results revealed that the diffusion coefficient in the gas phase increased with increasing gas-filled pore width and converged asymptotically toward the diffusion coefficient in the bulk state. It could be shown that the diffusion coefficient of gas molecules dissolved in the water films remained constant as long as the interacting water surface was in the bulk-liquid-like phase. This behavior changes in very thin water films. It was observed that the partitioning coefficient of gas molecules at the solid-liquid interface is nearly the same as that in the bulk-liquid-like phase. Partitioning is observed to be strongly dependent on the temperature and gas molecular weights. In the second part of the study, nonequilibrium molecular dynamics (NEMD) simulations were performed to investigate the mobility of gases in pressure-driven decoupled gas-phase dynamics (DGPD) and coupled gas and water phase dynamics (CGWPD) in a partially saturated Na-montmorillonite slit mesopore. The dynamic viscosity of the gas phase was calculated from NEMD simulations and indicated that the viscosity of the gas phase was almost the same in both methods (DGPD and CGWPD). The average slip length for gas molecules at the gas-water interface was also calculated, revealing that the slip-free boundary condition assumed in continuum models is generally invalid for microfluidics and that a slip boundary condition exists at the microscale for specific surface interactions. Finally, a Bosanquet-type equation was developed to predict the diffusion coefficient and dynamic viscosity of gas as a function of the average pore width, gas mean-free path, geometric factor, and thickness of the adsorbed water film.

Temperature dependence study of water dynamics in Fluorohectorite clays using Molecular dynamics simulations

Investigating Temperature Effects on the Shear Behavior of Clays: Molecular Dynamics Simulations

A Molecular Dynamics Study on the Structural and Mechanical Properties of Hydrated Kaolinite System under Tension

Molecular Dynamics Simulation of Diffusion and Electrical Conductivity in Montmorillonite Interlayers

The Affection to Water Clusters on Temperatures and Electric Fields by Molecular Dynamics Simulation

Molecular Simulation Study of Montmorillonite in Contact with Water

Molecular Dynamics Simulations of Carbon Dioxide Intercalation in Hydrated Na-Montmorillonite

Nanoscale Prediction of the Thermal, Mechanical, and Transport Properties of Hydrated Clay on 106- and 1015-Fold Larger Length and Time Scales

Dynamics of Water in LiCl and CaCl2 Aqueous Solutions Confined in Silica Matrices: A Backscattering Neutron Spectroscopy Study

Molecular Dynamics Simulation of Self-Diffusion Coefficient and Its Relation with Temperature Using Simple Lennard-Jones Potential

Diffusion and Gas Flow Dynamics in Partially Saturated Smectites

Molecular dynamics simulation of hydrogen diffusion in water-saturated clay minerals; implications for Underground Hydrogen Storage (UHS)

Elevated temperature effects (T > 100 °C) on the interfacial water and microstructure swelling of Na-montmorillonite

A First-Principles Study of the Effects of Temperature on the Atomic and Electronic Structures of Mg-montmorillonite with H2O Molecule Adsorption

Molecular Dynamics Modeling of Hydraulic Conductivity of Soil Considering Variable Viscosity and Adsorbed Water

Molecular Dynamics Simulation of Illite: from Particle Associations to Hydration Properties

Experimental Study of the Stability and Phase Relations of Clays at High Temperature in a Thermal Gradient

Molecular Simulation of Interlayer Structure and Dynamics Properties in Lithium Montmorillonite Hydrates

Molecular Dynamic Simulation of Interlayer Water Structure in Na-montmorillonite

MOLECULAR DYNAMICS SIMULATIONS FOR HYDRATION SWELLING AND INTERLAYER STRUCTURE OF Na-MONTMORILLONITE