Prediction of the phase equilibria of methane hydrates using the direct phase coexistence methodology
- Environmental Research Laboratory, National Center for Scientific Research NCSR “Demokritos,” Aghia Paraskevi, Attiki GR-15310 (Greece)
The direct phase coexistence method is used for the determination of the three-phase coexistence line of sI methane hydrates. Molecular dynamics (MD) simulations are carried out in the isothermal–isobaric ensemble in order to determine the coexistence temperature (T{sub 3}) at four different pressures, namely, 40, 100, 400, and 600 bar. Methane bubble formation that results in supersaturation of water with methane is generally avoided. The observed stochasticity of the hydrate growth and dissociation processes, which can be misleading in the determination of T{sub 3}, is treated with long simulations in the range of 1000–4000 ns and a relatively large number of independent runs. Statistical averaging of 25 runs per pressure results in T{sub 3} predictions that are found to deviate systematically by approximately 3.5 K from the experimental values. This is in good agreement with the deviation of 3.15 K between the prediction of TIP4P/Ice water force field used and the experimental melting temperature of ice Ih. The current results offer the most consistent and accurate predictions from MD simulation for the determination of T{sub 3} of methane hydrates. Methane solubility values are also calculated at the predicted equilibrium conditions and are found in good agreement with continuum-scale models.
- OSTI ID:
- 22416051
- Journal Information:
- Journal of Chemical Physics, Vol. 142, Issue 4; Other Information: (c) 2015 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA); ISSN 0021-9606
- Country of Publication:
- United States
- Language:
- English
Similar Records
Phase equilibrium of liquid water and hexagonal ice from enhanced sampling molecular dynamics simulations
Computer simulation of alkali metal cation-montmorillonite hydrates