Microscopic insights into poly- and mono-crystalline methane hydrate dissociation in Na-montmorillonite pores at static and dynamic fluid conditions

Knowledge on the kinetics of gas hydrate dissociation in clay pores at static and dynamic fluid conditions is a fundamental scientific issue for improving gas production efficiency from hydrate deposits using thermal stimulation and depressurization respectively. Here, molecular dynamics simulations were used to investigate poly- and mono...

Transport properties of mixtures of acid gases with aqueous monoethanolamine solutions: A molecular dynamics study

We investigated the effect of temperature and monoethanolamine (MEA) concentration on the self-diffusivity of acid gases, CO₂, and H₂S in aqueous MEA solutions. For this purpose, we computed densities of pure MEA and 30 wt% MEA/water solutions while scaling the LJ energy (ϵ) parameter and point charges of MEA. Results...

Kinetics of zeolite-catalyzed heptane hydroisomerization and hydrocracking with CBMC-modeled adsorption terms: Zeolite Beta as a large pore base case

A reactor model that deconvolutes thermodynamics of adsorption of hydrocarbon in the pores of zeolite Beta, obtained by Configurational-bias Monte Carlo simulations, from intrinsic, intraporous kinetics of hydroisomerization and hydrocracking reactions, provides a good quantitative description of all significant reactions in the kinetic...

Gibbs Ensemble Monte Carlo Simulation of Fluids in Confinement: Relation between the Differential and Integral Pressures

The accurate description of the behavior of fluids in nanoporous materials is of great importance for numerous industrial applications. Recently, a new approach was reported to calculate the pressure of nanoconfined fluids. In this approach, two different pressures are defined to take into account the smallness of the system: the so-called...

Multiple linear regression and thermodynamic fluctuations are equivalent for computing thermodynamic derivatives from molecular simulation

Partial molar properties are of fundamental importance for understanding properties of non-ideal mixtures. Josephson and co-workers (Mol. Phys. 2019, 117, 3589–3602) used least squares multiple linear regression to obtain partial molar properties in open constant-pressure ensembles. Assuming composition-independent partial molar properties...

Improving the accuracy of computing chemical potentials in CFCMC simulations

The CFCMC simulation methodology considers an expanded ensemble to solve the problem of low insertion/deletion acceptance probabilities in open ensembles. It allows for a direct calculation of the chemical potential by binning of the coupling parameter λ and using the probabilities p(λ = 0) and p(λ = 1), which require extrapolation. Here, we...

The dynamic behavior of gas hydrate dissociation by heating in tight sandy reservoirs: A molecular dynamics simulation study

Knowledge on the kinetics of gas hydrate dissociation in microporous sediments is very important for developing safe and efficient approaches to gas recovery from natural gas hydrate (NGH) deposits. Herein, molecular dynamics (MD) simulations are used to study the dissociation kinetics in microporous sediments. The hydrate phase occupies a...

Computation of partial molar properties using continuous fractional component Monte Carlo

An alternative method for calculating partial molar excess enthalpies and partial molar volumes of components in Monte Carlo (MC) simulations is developed. This method combines the original idea of Frenkel, Ciccotti, and co-workers with the recent continuous fractional component Monte Carlo (CFCMC) technique. The method is tested for a system...

Solubility of sulfur compounds in commercial physical solvents and an ionic liquid from Monte Carlo simulations

Natural gas, synthesis gas, and flue gas typically contain a large number of impurities (e.g., acidic gases), which should be removed to avoid environmental and technological problems, and to meet customer specifications. One approach is to use physical solvents to remove the acidic gases. If no experimental data are available, the solubility...

Product shape selectivity of MFI-type, MEL-type, and BEA-type zeolites in the catalytic hydroconversion of heptane

The influence of product shape selectivity on the bifunctional conversion of n-C₇ by zeolite catalysts is investigated. Three different zeolite catalysts with different pore sizes (MFI-type, MEL-type, and BEA-type zeolites) have been investigated experimentally. For all three catalysts, n-C₇ is isomerized to monobranched...

Computing equation of state parameters of gases from Monte Carlo simulations

Monte Carlo (MC) simulations in ensembles with a fixed chemical potential or fugacity, for example the grand-canonical or the osmotic ensemble, are often used to compute phase equilibria. Chemical potentials can be computed either with an equation of state (EoS) or from molecular simulations. The accuracy of the computed chemical potentials...

Solubilities of CO₂, CH₄, C₂H₆, and SO₂ in ionic liquids and Selexol from Monte Carlo simulations

Monte Carlo simulations are used to calculate the solubility of natural gas components in ionic liquids (ILs) and Selexol, which is a mixture of poly(ethylene glycol) dimethyl ethers. The solubility of the pure gases carbon dioxide (CO₂), methane (CH₄), ethane (C₂H₆), and sulfur dioxide (SO<sub>2<...

Computing bubble-points of CO₂/CH₄ gas mixtures in ionic liquids from Monte Carlo simulations

Computing bubble-points of multicomponent mixtures using Monte Carlo simulations is a non-trivial task. A new method is used to compute gas compositions from a known temperature, bubble-point pressure, and liquid composition. Monte Carlo simulations are used to calculate the bubble-points of carbon dioxide (CO₂) and methane (CH<sub...