Prediction of Adsorption and Diffusion of Shale Gas in Composite Pores Consisting of Kaolinite and Kerogen using Molecular Simulation.

Natural gas production from shale formations is one of the most recent and fast growing developments in the oil and gas industry. The accurate prediction of the adsorption and transport of shale gas is essential for estimating shale gas production capacity and improving existing extractions. To realistically represent heterogeneous shale formations, a composite pore model was built from a kaolinite slit mesopore hosting a kerogen matrix. Moreover, empty slabs (2, 3, and 4 nm) were added between the kerogen matrix and siloxane surface of kaolinite. Using Grand-Canonical Monte Carlo (GCMC) and molecular dynamics (MD) simulations, the adsorption and diffusion of pure methane, pure ethane, and a shale gas mixture were computed at various high pressures (100, 150, and 250 atm) and temperature of 298.15 K. The addition of an inner slit pore was found to significantly increase the excess adsorption of methane, as a pure component and in the shale gas mixture. The saturation of the composite pore with methane was observed to be at a higher pressure compared to ethane. The excess adsorption of carbon dioxide was not largely affected by pressure, and the local number density profile showed its strong affinity to kerogen micropores and the hydroxylated gibbsite surface under all conditions and pore widths. Lateral diffusion coefficients were found to increase with increasing the width of the empty slab inside the composite pore. Statistical errors of diffusion coefficients were found to be large for the case of shale gas components present at low composition. A larger composite pore configuration was created to investigate the diffusion of methane in different regions of the composite pore. The calculated diffusion coefficients and mean residence times were found to be indicative of the different adsorption mechanisms occurring inside the pore.

Solubilities and Transport Properties of CO2, Oxalic Acid, and Formic Acid in Mixed Solvents Composed of Deep Eutectic Solvents, Methanol, and Propylene Carbonate.

Recently, deep eutectic solvents (DES) have been considered as possible electrolytes for the electrochemical reduction of CO2 to value-added products such as formic and oxalic acids. The applicability of pure DES as electrolytes is hindered by high viscosities. Mixtures of DES with organic solvents can be a promising way of designing superior electrolytes by exploiting the advantages of each solvent type. In this study, densities, viscosities, diffusivities, and ionic conductivities of mixed solvents comprising DES (i.e., reline and ethaline), methanol, and propylene carbonate were computed using molecular simulations. To provide a quantitative assessment of the affinity and mass transport of CO2 and oxalic and formic acids in the mixed solvents, the solubilities and self-diffusivities of these solutes were also computed. Our results show that the addition of DES to the organic solvents enhances the solubilities of oxalic and formic acids, while the solubility of CO2 in the ethaline-containing mixtures are in the same order of magnitude with the respective pure organic components. A monotonic increase in the densities and viscosities of the mixed solvents is observed as the mole fraction of DES in the mixture increases, with the exception of the density of ethaline-propylene carbonate which shows the opposite behavior due to the high viscosity of the pure organic component. The self-diffusivities of all species in the mixtures significantly decrease as the mole fraction of DES approaches unity. Similarly, the self-diffusivities of the dissolved CO2 and the oxalic and formic acids also decrease by at least 1 order of magnitude as the composition of the mixture shifts from the pure organic component to pure DES. The computed ionic conductivities of all mixed solvents show a maximum value for mole fractions of DES in the range from 0.2 to 0.6 and decrease as more DES is added to the mixtures. Since for most mixtures studied here no prior experimental measurements exist, our findings can serve as a first data set based on which further investigation of DES-containing electrolyte solutions can be performed for the electrochemical reduction of CO2 to useful chemicals.

How sensitive are physical properties of choline chloride-urea mixtures to composition changes: Molecular dynamics simulations and Kirkwood-Buff theory.

Deep eutectic solvents (DESs) have emerged as a cheaper and greener alternative to conventional organic solvents. Choline chloride (ChCl) mixed with urea at a molar ratio of 1:2 is one of the most common DESs for a wide range of applications such as electrochemistry, material science, and biochemistry. In this study, molecular dynamics simulations are performed to study the effect of urea content on the thermodynamic and transport properties of ChCl and urea mixtures. With increased mole fraction of urea, the number of hydrogen bonds (HBs) between cation-anion and ion-urea decreases, while the number of HBs between urea-urea increases. Radial distribution functions (RDFs) for ChCl-urea and ChCl-ChCl pairs shows a significant decrease as the mole fraction of urea increases. Using the computed RDFs, Kirkwood-Buff Integrals (KBIs) are computed. KBIs show that interactions of urea-urea become stronger, while interactions of urea-ChCl and ChCl-ChCl pairs become slightly weaker with increasing mole fraction of urea. All thermodynamic factors are found larger than one, indicating a non-ideal mixture. Our results also show that self- and collective diffusivities increase, while viscosities decrease with increasing urea content. This is mainly due to the weaker interactions between ions and urea, resulting in enhanced mobilities. Ionic conductivities exhibit a non-monotonic behavior. Up to a mole fraction of 0.5, the ionic conductivities increase with increasing urea content and then reach a plateau.

Solubility of Carbon Dioxide, Hydrogen Sulfide, Methane, and Nitrogen in Monoethylene Glycol; Experiments and Molecular Simulation.

Knowledge on the solubility of gases, especially carbon dioxide (CO2), in monoethylene glycol (MEG) is relevant for a number of industrial applications such as separation processes and gas hydrate prevention. In this study, the solubility of CO2 in MEG was measured experimentally at temperatures of 333.15, 353.15, and 373.15 K. Experimental data were used to validate Monte Carlo (MC) simulations. Continuous fractional component MC simulations in the osmotic ensemble were performed to compute the solubility of CO2 in MEG at the same temperatures and at pressures up to 10 bar. MC simulations were also used to study the solubility of methane (CH4), hydrogen sulfide (H2S), and nitrogen (N2) in MEG at 373.15 K. Solubilities from experiments and simulations are in good agreement at low pressures, but deviations were observed at high pressures. Henry coefficients were also computed using MC simulations and compared to experimental values. The order of solubilities of the gases in MEG at 373.15 K was computed as H2S > CO2 > CH4 > N2. Force field modifications may be required to improve the prediction of solubilities of gases in MEG at high pressures and low temperatures.

Kirkwood-Buff Integrals Using Molecular Simulation: Estimation of Surface Effects.

Kirkwood-Buff (KB) integrals provide a connection between microscopic properties and thermodynamic properties of multicomponent fluids. The estimation of KB integrals using molecular simulations of finite systems requires accounting for finite size effects. In the small system method, properties of finite subvolumes with different sizes embedded in a larger volume can be used to extrapolate to macroscopic thermodynamic properties. KB integrals computed from small subvolumes scale with the inverse size of the system. This scaling was used to find KB integrals in the thermodynamic limit. To reduce numerical inaccuracies that arise from this extrapolation, alternative approaches were considered in this work. Three methods for computing KB integrals in the thermodynamic limit from information of radial distribution functions (RDFs) of finite systems were compared. These methods allowed for the computation of surface effects. KB integrals and surface terms in the thermodynamic limit were computed for Lennard-Jones (LJ) and Weeks-Chandler-Andersen (WCA) fluids. It was found that all three methods converge to the same value. The main differentiating factor was the speed of convergence with system size L. The method that required the smallest size was the one which exploited the scaling of the finite volume KB integral multiplied by L. The relationship between KB integrals and surface effects was studied for a range of densities.

Brick-CFCMC: Open Source Software for Monte Carlo Simulations of Phase and Reaction Equilibria Using the Continuous Fractional Component Method.

We present a new molecular simulation code, Brick-CFCMC, for performing Monte Carlo simulations using state-of-the-art simulation techniques. The Continuous Fractional Component (CFC) method is implemented for simulations in the NVT/NPT ensembles, the Gibbs Ensemble, the Grand-Canonical Ensemble, and the Reaction Ensemble. Molecule transfers are facilitated by the use of fractional molecules which significantly improve the efficiency of the simulations. With the CFC method, one can obtain phase equilibria and properties such as chemical potentials and partial molar enthalpies/volumes directly from a single simulation. It is possible to combine trial moves from different ensembles. This enables simulations of phase equilibria in a system where also a chemical reaction takes place. We demonstrate the applicability of our software by investigating the esterification of methanol with acetic acid in a two-phase system.