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.

Diffusion of Gases in Amorphous Polymers: The Monte Carlo Void Method Based on Molecular Dynamics; Application to He Diffusion in Polyethylene

A Monte Carlo Study Of The Glass Transition Of Polymethylene

Diffusivity of Liquid Mixture Confined in Silicate Pore by Molecular Dynamic Simulation

MONTE CARLO SIMULATION OF THE DYNAMICS AND CONFIGURATIONAL DEPENDENT PROPERTIES OF POLYMETHYLENE CHAINS

Measurement of infinite dilution diffusion coefficients of small molecule solvents in nascent polyethylene particles by inverse gas chromatography

A Molecular Dynamics Model Of The Amorphous Regions Of Polyethylene Oxide

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

Molecular Dynamics Simulation Study of Solid Vibration Permeation in Microporous Amorphous Silica Network Voids

Solvent Diffusion In Amorphous Glassy Polymers

Diffusion of small particles in a solid polymeric medium

Application of Molecular Dynamics Simulations in Molecular Property Prediction II: Diffusion Coefficient.

Estimating position-dependent and anisotropic diffusivity tensors from molecular dynamics trajectories: Existing methods and future outlook

Self-diffusion coefficient and viscosity of methane and carbon dioxide via molecular dynamics simulations based on new ab initio-derived force fields

Multiscale Simulations of Graphite-Capped Polyethylene Melts: Brownian Dynamics/Kinetic Monte Carlo Compared to Atomistic Calculations and Experiment

Studying polymer diffusiophoresis with Non-Equilibrium Molecular Dynamics

Gas permeation through a polymer network

Diffusion and Gas Flow Dynamics in Partially Saturated Smectites

All-atom molecular dynamics simulation of solvent diffusion in unentangled polystyrene film

Gas Adsorption and Diffusion Behaviors in Interfacial Systems Composed of a Polymer of Intrinsic Microporosity and Amorphous Silica: A Molecular Simulation Study

Influence of molecular shape on self-diffusion under severe confinement: A molecular dynamics study

Mathematical Models for Estimating Diffusion Coefficients in Concentrated Polymer Solutions from Experimental Data