Factors governing the behaviour of aqueous methane in narrow pores

Molecular Density Fluctuations Control Solubility and Diffusion for Confined Aqueous Hydrogen

Hydrogen's contribution to a sustainable energy transformation requires intermittent storage technologies, e.g., underground hydrogen storage (UHS). Toward designing UHS sites, atomistic molecular dynamics (MD) simulations are used here to quantify thermodynamic and transport properties for confined aqueous H2. Slit-shaped pores of width 10 and 20 Å are carved out of kaolinite. Within these pores, water yields pronounced hydration layers. Molecular H2 distributes along these hydration layers, yielding solubilities up to ∼25 times those in the bulk. Hydrogen accumulates near the siloxane surface, where water density fluctuates significantly. On the contrary, a dense hydration layer forms on the gibbsite surface, which is, for the most part, depleted of H2. Although confinement reduces water mobility, the diffusion of aqueous H2 increases as the kaolinite pore width decreases, also a consequence of water density fluctuations. These results relate to H2 permeability in underground hydrogen storage sites.

Breakthrough pressure of oil displacement by water through the ultra-narrow kerogen pore throat from the Young-Laplace equation and molecular dynamic simulations

As one common unconventional reservoir, shale plays a pivotal role to compensate the depletion of conventional oil resources. There are numerous nanoscale pores and ultra-narrow pore throats (sub 2 nm) in shale media. To displace oil through ultra-narrow pore throats by water, one needs to overcome excessively-high capillary pressure. Understanding the water-oil two-phase displacement process through pore throats is critical to numerical simulation on tight/shale oil exploitation and ultimate oil recovery estimation. In this work, we use molecular dynamics simulations to investigate oil (represented by n-octane) displacement by water through a ~2 nm kerogen (represented by Type II-C kerogen) pore throat. Besides, the applicability of the Young-Laplace equation to the ultra-narrow kerogen pore throat has been assessed. We find that although the Type II-C kerogen is generally oil-wet, water has an excellent displacement efficiency without the oil film on the substrate, thanks to the hydrogen bonding formed between water and heteroatoms (such as O, N, and S) on the kerogen surface. Unlike previous studies, the capillary pressure obtained from the widely used Young-Laplace equation shows a good agreement with the breakthrough pressure obtained from MD simulations for the ∼2 nm kerogen pore throat. Our work indicates that explicitly considering intermolecular interactions as well as atomistic and molecular level characteristics is imperative to study the two-phase displacement process through ultra-narrow pore throats.

Homogenization of Continuum-Scale Transport Properties from Molecular Dynamics Simulations: An Application to Aqueous-Phase Methane Diffusion in Silicate Channels

Silica-based materials including zeolites are commonly used for wide-ranging applications including separations and catalysis. Substrate transport rates in these materials can significantly influence the efficiency of such applications. Two factors that contribute to transport rates include (1) the porosity of the silicate matrix and (2) nonbonding interactions between the diffusing species and the silicate surface. These contributions generally emerge from disparate length scales, namely, "microscopic" (roughly nanometer-scale) and "macroscopic" (roughly micron-scale), respectively. Here, we develop a simulation framework to estimate the simultaneous impact of these factors on methane mass transport in silicate channels. Specifically, we develop a model of methane transport using homogenization theory to obtain transport parameters valid at length scales of hundreds to thousands of nanometers. These parameters implicitly reflect interactions taking place at fractions of a nanometer. The inputs to the homogenization analysis are data from extensive molecular dynamics simulations that incorporate atomistic-scale interactions, processed to yield local diffusion coefficients and mean force potentials. With this model, we demonstrate how nuances in silicate hydration and silica/methane interactions impact methane diffusion rates in silicate materials, including the effects of silicate surface chemistry such as the presence of silanol groups. The molecular dynamics simulations indicate that methane diffusivity at the silica surface is lower than the bulk-like rates observed at the center of channels of sufficient width. However, potentials of mean force generally evidence attractive methane/silica interactions that enhance diffusion overall. By simultaneously accounting for both of these effects, we show that the effective diffusion coefficient for the nanoporous silicate can be approximately double the value of estimates assuming fully bulk-like behavior in the channel. This study therefore demonstrates the importance of determining diffusion coefficients and potentials of mean force at an atomistic level when estimating transport properties in bulk materials. Importantly, we provide a simple homogenization framework to incorporate these molecular-scale attributes into bulk material transport estimates. This hybrid homogenization/molecular dynamics approach will be of general use for describing small-molecule transport in materials with detailed molecular interactions.

A Novel Modeling Approach to Stochastically Evaluate the Impact of Pore Network Geometry, Chemistry and Topology on Fluid Transport

Fine-grained sandstones, siltstones, and shales have become increasingly important to satisfy the ever-growing global energy demands. Of particular current interest are shale rocks, which are mudstones made up of organic and inorganic constituents of varying pore sizes. These materials exhibit high heterogeneity, low porosity, varying chemical composition and low pore connectivity. Due to the complexity and the importance of such materials, many experimental, theoretical and computational efforts have attempted to quantify the impact of rock features on fluids diffusivity and ultimately on permeability. In this study, we introduce a stochastic kinetic Monte Carlo approach developed to simulate fluid transport. The features of this approach allow us to discuss the applicability of 2D vs 3D models for the calculation of transport properties. It is found that a successful model should consider realistic 3D pore networks consisting of pore bodies that communicate via pore throats, which however requires a prohibitive amount of computational resources. To overcome current limitations, we present a rigorous protocol to stochastically generate synthetic 3D pore networks in which pore features can be isolated and varied systematically and individually. These synthetic networks do not correspond to real sample scenarios but are crucial to achieve a systematic evaluation of the pore features on the transport properties. Using this protocol, we quantify the contribution of the pore network’s connectivity, porosity, mineralogy, and pore throat width distribution on the diffusivity of supercritical methane. A sensitivity analysis is conducted to rank the significance of the various network features on methane diffusivity. Connectivity is found to be the most important descriptor, followed by pore throat width distribution and porosity. Based on such insights, recommendations are provided on possible technological approaches to enhance fluid transport through shale rocks and equally complex pore networks. The purpose of this work is to identify the significance of various pore network characteristics using a stochastic KMC algorithm to simulate the transport of fluids. Our findings could be relevant for applications that make use of porous media, ranging from catalysis to radioactive waste management, and from environmental remediation to shale gas production.

Structure and dynamics of ethane confined in silica nanopores in the presence of CO2

Fundamental understanding of the subcritical/supercritical behavior of key hydrocarbon species inside nano-porous matrices at elevated pressure and temperature is less developed compared to bulk fluids, but this knowledge is of great importance for chemical and energy engineering industries. This study explores in detail the structure and dynamics of ethane (C2H6) fluid confined in silica nanopores, with a focus on the effects of pressure and different ratios of C2H6 and CO2 at non-ambient temperature. Quasi-elastic neutron scattering (QENS) experiments were carried out for the pure C2H6, C2H6:CO2 = 3:1, and 1:3 mixed fluids confined in 4-nm cylindrical silica pores at three different pressures (30 bars, 65 bars, and 100 bars) at 323 K. Two Lorentzian functions were required to fit the spectra, corresponding to fast and slow translational motions. No localized motions (rotations and vibrations) were detected. Higher pressures resulted in hindrances of the diffusivity of C2H6 molecules in all systems investigated. Pore size was found to be an important factor, i.e., the dynamics of confined C2H6 is more restricted in smaller pores compared to the larger pores used in previous studies. Molecular dynamics simulations were performed to complement the QENS experiment at 65 bars, providing supportive structure information and comparable dynamic information. The simulations indicate that CO2 molecules are more strongly attracted to the pore surface compared to C2H6. The C2H6 molecules interacting with or near the pore surface form a dense first layer (L1) close to the pore surface and a second less dense layer (L2) extending into the pore center. Both the experiments and simulations revealed the role that CO2 molecules play in enhancing C2H6 diffusion ("molecular lubrication") at high CO2:C2H6 ratios. The energy scales of the two dynamic components, fast and slow, quantified by both techniques, are in very good agreement. Herein, the simulations identified the fast component as the main contributor to the dynamics. Molecule motions in the L2 region are mostly responsible for the dynamics (fast and slow) that can be detected by the instrument.

Understanding methane/carbon dioxide partitioning in clay nano- and meso-pores with constant reservoir composition molecular dynamics modeling

CRC-MD simulations show that nanopores in shales bounded by clay minerals have a strong preference for CO₂ relative to CH₄.

Reduction in the Repulsive Forces between Two Charged Surfaces in Aqueous Solutions Containing Salts by a Liquid Flow

In spite of the fact that a flow is often present in the liquid in which charged particles are dispersed, the effect of a flow on the forces controlling the dispersion is not clear. Here, we used a combined atomic force microscope-peristaltic pump system to determine the effect of a flow in aqueous solutions between a negatively charged silica particle and a negatively charged silicon wafer on the forces in the system. The effect of a flow on the forces in water or aqueous solutions of NaCl or MgCl₂·6H₂O was studied for salt concentrations lower than the concentrations needed to invert the charge of the silica and silicon surfaces. This was done to prevent the formation of a reversed flow in the system due to a charge inversion of the silica surface. A flow was seen to decrease the intersurface repulsive forces, if the water contained salt (NaCl or MgCl₂·6H₂O). An increased bulk salt concentration was also seen to decrease the repulsive forces further in the presence of a liquid flow. The surface potentials and effective ionic concentrations of the systems were determined by comparing the experimental curves with the theoretically calculated ones. The surface potentials and effective ionic concentrations were seen to decrease and increase, respectively, as the flow rate and bulk salt concentrations were increased. This change was explained by the shrinking of the diffuse layers by the liquid flow, due to part of the diffuse layer being washed away by the flowing liquid.

Confinement Effects on Carbon Dioxide Methanation: A Novel Mechanism for Abiotic Methane Formation

An important scientific debate focuses on the possibility of abiotic synthesis of hydrocarbons during oceanic crust-seawater interactions. While on-site measurements near hydrothermal vents support this possibility, laboratory studies have provided data that are in some cases contradictory. At conditions relevant for sub-surface environments it has been shown that classic thermodynamics favour the production of CO2 from CH4, while abiotic methane synthesis would require the opposite. However, confinement effects are known to alter reaction equilibria. This report shows that indeed thermodynamic equilibrium can be shifted towards methane production, suggesting that thermal hydrocarbon synthesis near hydrothermal vents and deeper in the magma-hydrothermal system is possible. We report reactive ensemble Monte Carlo simulations for the CO2 methanation reaction. We compare the predicted equilibrium composition in the bulk gaseous phase to that expected in the presence of confinement. In the bulk phase we obtain excellent agreement with classic thermodynamic expectations. When the reactants can exchange between bulk and a confined phase our results show strong dependency of the reaction equilibrium conversions, [Formula: see text], on nanopore size, nanopore chemistry, and nanopore morphology. Some physical conditions that could shift significantly the equilibrium composition of the reactive system with respect to bulk observations are discussed.

Confined Water Determines Transport Properties of Guest Molecules in Narrow Pores

We computed the transport of methane through 1 nm wide slit-shaped pores carved out of solid substrates. Models for silica, magnesium oxide, and alumina were used as solid substrates. The pores were filled with water. The results show that the methane permeability through the hydrated pores is strongly dependent on the solid substrate. Detailed analysis of the simulated systems reveals that local properties of confined water, including its structure, and more importantly, evolution of solvation free energy and hydrogen bond structure are responsible for the pronounced differences observed. The simulations are extended to multicomponent systems representative of natural gas, containing methane, ethane, and H2S. The results show that all pores considered have high affinity for H2S, moderate affinity for methane, and low affinity for ethane. The H2S/methane transport selectivity through the hydrated alumina pore is comparable, or superior, to that reported for existing commercial membranes. A multiscale approach was then implemented to demonstrate that a Smoluchowski one-dimensional model is able to reproduce the molecular-level results for short pores when appropriate values for the local self-diffusion coefficients are used as input parameters. We propose that the model can be extended to predict methane transport through uniform hydrated pores of macroscopic length. When verified by experiments, our simulation results could have important implications in applications such as natural gas sweetening and predictions of methane migration through hydraulically fractured shale formations.

The Carbon-Water Interface: Modeling Challenges and Opportunities for the Water-Energy Nexus

Providing clean water and sufficient affordable energy to all without compromising the environment is a key priority in the scientific community. Many recent studies have focused on carbon-based devices in the hope of addressing this grand challenge, justifying and motivating detailed studies of water in contact with carbonaceous materials. Such studies are becoming increasingly important because of the miniaturization of newly proposed devices, with ubiquitous nanopores, large surface-to-volume ratio, and many, perhaps most of the water molecules in contact with a carbon-based surface. In this brief review, we discuss some recent advances obtained via simulations and experiments in the development of carbon-based materials for applications in water desalination. We suggest possible ways forward, with particular emphasis on the synergistic combination of experiments and simulations, with simulations now sometimes offering sufficient accuracy to provide fundamental insights. We also point the interested reader to recent works that complement our short summary on the state of the art of this important and fascinating field.