Abstract:This paper explores the three-phase equilibrium of structure-I methane hydrate by utilizing the extracted Mie potential parameters for non-bonded molecular interactions in the LAMMPS molecular dynamics simulator. An in-house code was developed to numerically solve the SAFT-VR Mie equation of state (EOS). To reduce numerical errors, unlike common practices, we incorporated all partial derivatives analytically in the code. The microscopic simulation parameters for methane and water molecules were selected so that saturation pressure and liquid density data from EOS were fitted on experimental data. Using this approach, the three-phase ( V−Lw−H ) equilibrium temperatures for pressure values of 30, 70, 100, and 300 bar were calculated. The solubility of methane in the aqueous phase at three-phase equilibrium was also determined. The results were validated against available laboratory data. To evaluate the performance of Mie potential parameters, the MD simulations were repeated using three well-known force fields for methane molecules, namely, Optimized Potentials for Liquid Simulations-United Atom (OPLS-UA), Universal force field (UFF), and General Amber force field (GAFF). The TIP4P/Ice force field was used for water molecules in all cases. The comparison with experimental data reveals that while the accuracy of the UFF force field is the highest, our simulation approach outperforms OPLS-UA and GAFF force fields in predicting methane solubility in the aqueous phase. Besides molecular dynamics simulation, the van der Waals-Platteeuw model was combined with the SAFT VR Mie EOS to account for the hydrate phase. In our case study, the SAFT-VR Mie + vdW-P model proved to be more accurate in predicting the three-phase equilibrium data compared to the MD approach.

Simulation of tetrahydrofuran-water hydrate formation with a combined cellular automata method and finite difference

Simulation of the THF hydrate-water interfacial free energy from computer simulation

Experimental study on the formation and agglomeration of tetrahydrofuran hydrate under flowing condition

Prediction of the univariant two-phase coexistence line of the tetrahydrofuran hydrate from computer simulation

Effect of the Temperature and Tetrahydrofuran (THF) Concentration on THF Hydrate Formation in Aqueous Solution

Microscopic morphology evolution of the crystal structure of tetrahydrofuran hydrate under flowing condition

Study of THF Hydrate Crystallization Based on In Situ Observation with Atomic Force Microscopy

Coupled thermal-hydrodynamic-mechanical–chemical numerical simulation for gas production from hydrate-bearing sediments based on hybrid finite volume and finite element method

Molecular dynamics simulation of methane hydrates: Prediction of the phase equilibria using extracted microscopic parameters from SAFT-VR Mie EOS

Kinetic studies of the secondary hydrate formation in porous media based on experiments in a cubic hydrate simulator and a new kinetic model

Experiment design as a tool for hydrate phase equilibria prediction of CO 2 + CH 4 gas mixture in the presence of Tera-Butylammonium Bromide in L W –H–V systems

Effects of Additive Mixture (THF/SDS) on the Thermodynamic and Kinetic Properties of CO2/H2 Hydrate in Porous Media

Three-phase equilibria of hydrates from computer simulation. I. Finite-size effects in the methane hydrate

Differential Scanning Calorimetry Studies of Clathrate Hydrate Formation

Hydrate Formation Conditions Of Methane Plus Ethylene Plus Tetrahydrofuran Plus Water Systems

Experimental measurement and thermodynamic modeling of hydrate stability conditions in the methane + cyclopentane + tetra-n-butyl ammonium chloride (TBAC) aqueous solution system

Study of Vapor–hydrate Two-Phase Equilibria

A Simulator Based on Coupling of Reaction Transport Model and Multiphase Hydrate Simulator and Its Application to Studies of Methane Transportation in Marine Sediments