Cage occupancy of methane hydrates from Gibbs ensemble Monte Carlo simulations

Paul E. Brumby, Daisuke Yuhara, David T. Wu, Amadeu K. Sum, Kenji Yasuoka

Research output: Contribution to journalArticle

17 Citations (Scopus)

Abstract

Isobaric-isothermal Gibbs ensemble Monte Carlo simulations of sI methane hydrates in equilibrium with bulk methane are performed to calculate large and small cage occupancies. The OPLS united-atom Lennard-Jones potential, and a variation, are used to model methane, while the TIP4P/Ice model represents water. This model system produces systematically higher total cage occupancies than those of simulations by Henley and Lucia [1] using TIP4P-Ew water and TraPPE-UA methane. We also see higher total occupancies than those observed in prior Raman spectrometry studies by Uchida et al. [2], but find good agreement for large cage occupancies. The simulations provide comprehensive predictions of large and small cage occupancies for temperatures and pressures within the ranges 270-290 K and 20-400 bar, respectively. These predictions demonstrate the advantages of the Gibbs ensemble for the simulation of clathrate hydrates; the ability to control the temperature and pressure, while not constraining the chemical composition of the hydrate phase is the key advantage of this ensemble.

Original languageEnglish
Pages (from-to)242-248
Number of pages7
JournalFluid Phase Equilibria
Volume413
DOIs
Publication statusPublished - 2016 Apr 15

Fingerprint

Methane
Hydrates
hydrates
methane
simulation
Lennard-Jones potential
Water
clathrates
Ice
predictions
Spectrometry
water
chemical composition
ice
Atoms
Temperature
temperature
Monte Carlo simulation
Chemical analysis
spectroscopy

Keywords

  • Clathrate hydrate
  • Gibbs ensemble
  • Methane
  • Monte Carlo simulation
  • Occupancy

ASJC Scopus subject areas

  • Chemical Engineering(all)
  • Physical and Theoretical Chemistry
  • Physics and Astronomy(all)

Cite this

Cage occupancy of methane hydrates from Gibbs ensemble Monte Carlo simulations. / Brumby, Paul E.; Yuhara, Daisuke; Wu, David T.; Sum, Amadeu K.; Yasuoka, Kenji.

In: Fluid Phase Equilibria, Vol. 413, 15.04.2016, p. 242-248.

Research output: Contribution to journalArticle

Brumby, Paul E. ; Yuhara, Daisuke ; Wu, David T. ; Sum, Amadeu K. ; Yasuoka, Kenji. / Cage occupancy of methane hydrates from Gibbs ensemble Monte Carlo simulations. In: Fluid Phase Equilibria. 2016 ; Vol. 413. pp. 242-248.
@article{d8dca293e35b4cddabd0d89eeed05e0f,
title = "Cage occupancy of methane hydrates from Gibbs ensemble Monte Carlo simulations",
abstract = "Isobaric-isothermal Gibbs ensemble Monte Carlo simulations of sI methane hydrates in equilibrium with bulk methane are performed to calculate large and small cage occupancies. The OPLS united-atom Lennard-Jones potential, and a variation, are used to model methane, while the TIP4P/Ice model represents water. This model system produces systematically higher total cage occupancies than those of simulations by Henley and Lucia [1] using TIP4P-Ew water and TraPPE-UA methane. We also see higher total occupancies than those observed in prior Raman spectrometry studies by Uchida et al. [2], but find good agreement for large cage occupancies. The simulations provide comprehensive predictions of large and small cage occupancies for temperatures and pressures within the ranges 270-290 K and 20-400 bar, respectively. These predictions demonstrate the advantages of the Gibbs ensemble for the simulation of clathrate hydrates; the ability to control the temperature and pressure, while not constraining the chemical composition of the hydrate phase is the key advantage of this ensemble.",
keywords = "Clathrate hydrate, Gibbs ensemble, Methane, Monte Carlo simulation, Occupancy",
author = "Brumby, {Paul E.} and Daisuke Yuhara and Wu, {David T.} and Sum, {Amadeu K.} and Kenji Yasuoka",
year = "2016",
month = "4",
day = "15",
doi = "10.1016/j.fluid.2015.10.005",
language = "English",
volume = "413",
pages = "242--248",
journal = "Fluid Phase Equilibria",
issn = "0378-3812",
publisher = "Elsevier",

}

TY - JOUR

T1 - Cage occupancy of methane hydrates from Gibbs ensemble Monte Carlo simulations

AU - Brumby, Paul E.

AU - Yuhara, Daisuke

AU - Wu, David T.

AU - Sum, Amadeu K.

AU - Yasuoka, Kenji

PY - 2016/4/15

Y1 - 2016/4/15

N2 - Isobaric-isothermal Gibbs ensemble Monte Carlo simulations of sI methane hydrates in equilibrium with bulk methane are performed to calculate large and small cage occupancies. The OPLS united-atom Lennard-Jones potential, and a variation, are used to model methane, while the TIP4P/Ice model represents water. This model system produces systematically higher total cage occupancies than those of simulations by Henley and Lucia [1] using TIP4P-Ew water and TraPPE-UA methane. We also see higher total occupancies than those observed in prior Raman spectrometry studies by Uchida et al. [2], but find good agreement for large cage occupancies. The simulations provide comprehensive predictions of large and small cage occupancies for temperatures and pressures within the ranges 270-290 K and 20-400 bar, respectively. These predictions demonstrate the advantages of the Gibbs ensemble for the simulation of clathrate hydrates; the ability to control the temperature and pressure, while not constraining the chemical composition of the hydrate phase is the key advantage of this ensemble.

AB - Isobaric-isothermal Gibbs ensemble Monte Carlo simulations of sI methane hydrates in equilibrium with bulk methane are performed to calculate large and small cage occupancies. The OPLS united-atom Lennard-Jones potential, and a variation, are used to model methane, while the TIP4P/Ice model represents water. This model system produces systematically higher total cage occupancies than those of simulations by Henley and Lucia [1] using TIP4P-Ew water and TraPPE-UA methane. We also see higher total occupancies than those observed in prior Raman spectrometry studies by Uchida et al. [2], but find good agreement for large cage occupancies. The simulations provide comprehensive predictions of large and small cage occupancies for temperatures and pressures within the ranges 270-290 K and 20-400 bar, respectively. These predictions demonstrate the advantages of the Gibbs ensemble for the simulation of clathrate hydrates; the ability to control the temperature and pressure, while not constraining the chemical composition of the hydrate phase is the key advantage of this ensemble.

KW - Clathrate hydrate

KW - Gibbs ensemble

KW - Methane

KW - Monte Carlo simulation

KW - Occupancy

UR - http://www.scopus.com/inward/record.url?scp=84957838382&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84957838382&partnerID=8YFLogxK

U2 - 10.1016/j.fluid.2015.10.005

DO - 10.1016/j.fluid.2015.10.005

M3 - Article

AN - SCOPUS:84957838382

VL - 413

SP - 242

EP - 248

JO - Fluid Phase Equilibria

JF - Fluid Phase Equilibria

SN - 0378-3812

ER -