Thermodynamics and kinetics of methane hydrate formation in nanoporous media : theory and molecular simulation

Methane hydrate is a non-stoichiometric crystal in which water molecules form hydrogen-bonded cages that entrap methane molecules. Abundant methane hydrate resources can be found on Earth, especially trapped in mineral porous rocks (e.g., clay, permafrost, seafloor, etc.). For this reason, understan...

Full description

Bibliographic Details
Main Author: Jin, Dongliang
Other Authors: Laboratoire Interdisciplinaire de Physique Saint Martin d’Hères (LIPhy), Centre National de la Recherche Scientifique (CNRS)-Université Grenoble Alpes 2016-2019 (UGA 2016-2019 ), Université Grenoble Alpes, Benoît Coasne
Format: Doctoral or Postdoctoral Thesis
Language:English
Published: HAL CCSD 2018
Subjects:
Online Access:https://tel.archives-ouvertes.fr/tel-02181782
https://tel.archives-ouvertes.fr/tel-02181782/document
https://tel.archives-ouvertes.fr/tel-02181782/file/JIN_2018_diffusion.pdf
Description
Summary:Methane hydrate is a non-stoichiometric crystal in which water molecules form hydrogen-bonded cages that entrap methane molecules. Abundant methane hydrate resources can be found on Earth, especially trapped in mineral porous rocks (e.g., clay, permafrost, seafloor, etc.). For this reason, understanding the thermodynamics and formation kinetics of methane hydrate confined in porous media is receiving a great deal of attention. In this thesis, we combine computer modeling and theoretical approaches to determine the thermodynamics and formation kinetics of methane hydrate confined in porous media. First, the state-of-the-art on the thermodynamics and formation kinetics of methane hydrate is presented. Second, different molecular simulation strategies, including free energy calculations using the Einstein molecule approach, the direct coexistence method, and the hyperparallel tempering technique, are used to assess the phase stability of bulk methane hydrate at various temperatures and pressures. Third, among these strategies, the direct coexistence method is chosen to determine the shift in melting point upon confinement in pores, ∆Tm=Tm^{pore}-Tm^{bulk} where Tm^{pore} and Tm^{bulk} are the melting temperatures of bulk and confined methane hydrate. We found that confinement decreases the melting temperature, Tm^{pore}≺Tm^{bulk}. The shift in melting temperature using the direct coexistence method is consistent with the Gibbs-Thompson equation which predicts that the shift in melting temperature linearly depends on the reciprocal of pore width, i.e., ∆Tm/Tm^{bulk}∼k{GB}/Dp. The quantitative validity of this classical thermodynamic equation to describe such confinement and surface effects is also addressed. The surface tensions of methane hydrate-substrate and liquid water-substrate interfaces are determined using molecular dynamics to quantitatively validate the Gibbs-Thompson equation. Molecular dynamics simulations are also performed to determine important thermodynamic properties of bulk and confined methane ...