• Volume/Page
  • Keyword
  • DOI
  • Citation
  • Advanced
   
 
 
 

You Tube Flickr Twitter UniPHY Group iResearch App Facebook

J. Chem. Phys. 124, 234710 (2006); doi:10.1063/1.2207138 (12 pages)

Multiscale approach to CO2 hydrate formation in aqueous solution: Phase field theory and molecular dynamics. Nucleation and growth

György Tegze1, Tamás Pusztai1, Gyula Tóth1, László Gránásy1, Atle Svandal2, Trygve Buanes2, Tatyana Kuznetsova2, and Bjørn Kvamme2

1Research Institute for Solid State Physics and Optics, P.O. Box 49, H-1525 Budapest, Hungary
2Institute of Physics and Technology, University of Bergen, Allégaten 55, N-5007 Bergen, Norway

View MapView Map

(Received 14 December 2005; accepted 27 April 2006; published online 20 June 2006)

A phase field theory with model parameters evaluated from atomistic simulations/experiments is applied to predict the nucleation and growth rates of solid CO2 hydrate in aqueous solutions under conditions typical to underwater natural gas hydrate reservoirs. It is shown that under practical conditions a homogeneous nucleation of the hydrate phase can be ruled out. The growth rate of CO2 hydrate dendrites has been determined from phase field simulations as a function of composition while using a physical interface thickness (0.85±0.07 nm) evaluated from molecular dynamics simulations. The growth rate extrapolated to realistic supersaturations is about three orders of magnitude larger than the respective experimental observation. A possible origin of the discrepancy is discussed. It is suggested that a kinetic barrier reflecting the difficulties in building the complex crystal structure is the most probable source of the deviations.

© 2006 American Institute of Physics

Article Outline

  1. INTRODUCTION
  2. MULTISCALE APPROACH TO GAS HYDRATE FORMATION
    1. Phase field theory
      1. Equilibrium interface
      2. Calculation of the nucleation barrier
      3. Phase field simulation of single crystal growth
      4. Numerical implementation
    2. Molecular dynamics
  3. PHYSICAL PROPERTIES
  4. RESULTS AND DISCUSSION
    1. The planar hydrate-solution interface
      1. Molecular dynamics
      2. Phase field theory
    2. Nucleation of CO2 hydrate
    3. Growth of CO2 hydrate dendrites in aqueous solution
  5. SUMMARY

RELATED DATABASES

To view database links for this article, you need to log in.

KEYWORDS and PACS

PACS

  • 64.60.Q-

    Nucleation

  • 61.50.Ah

    Theory of crystal structure, crystal symmetry; calculations and modeling

  • 68.70.+w

    Whiskers and dendrites (growth, structure, and nonelectronic properties)

PUBLICATION DATA

ISSN:

0021-9606 (print)  
1089-7690 (online)

For access to fully linked references, you need to log in.
    Y. C. Shen and D. W. Oxtoby, J. Chem. Phys. 105, 6517 (1996)JCPSA6000105000015006517000001.

    J. W. Cahn and J. E. Hilliard, J. Chem. Phys. 28, 258 (1958)JCPSA6000028000002000258000001.

    J. W. Cahn and J. E. Hilliard, J. Chem. Phys. 31, 688 (1959)JCPSA6000031000003000688000001.

    A. Masuhr, T. A. Waniuk, R. Busch, and W. L. Johnson, Phys. Rev. Lett. 82, 2290 (1999)
    K. L. Ngai, J. H. Magill, and D. J. Plazek, J. Chem. Phys. 112, 1887 (2000)JCPSA6000112000004001887000001;, S. F. Swallen, P. A. Bonvallet, R. J. McMahon, and M. D. Ediger, Phys. Rev. Lett. 90, 015901 (2003).

    M. Conti, Phys. Rev. E 56, 3197 (1997).

    (a) A. Karma and W. J. Rappel, Phys. Rev. E 53, R3017 (1996)
    (b) 57, 4323 (1998)
    (c) A. Karma, Phys. Rev. Lett. 87, 115701 (2001)
    (d) B. Etchebarria, R. Folch, A. Karma, and M. Plapp, Phys. Rev. E 70, 061604 (2004).

    G. J. Martyna, D. J. Tobias, and M. L. Klein, J. Chem. Phys. 101, 4177 (1994)JCPSA6000101000005004177000001.

    W. G. Hoover, Phys. Rev. A 31, 1695 (1985).

    T. Bryk and A. D. J. Haymet, J. Chem. Phys. 117, 10258 (2002)JCPSA6000117000022010258000001.

    D. Turnbull, J. Appl. Phys. 21, 1022 (1950)JAPIAU000021000010001022000001;, J. Chem. Phys. 18, 768 (1950)JCPSA6000018000005000768000002.

    Various theoretical studies indicate that the induction time of nucleation theta is proportional to 1/D, where D is the self-diffusion coefficient in the liquid phase [see, e.g., K. F. Kelton, A. L. Greer, and C. V. Thompson, J. Chem. Phys. 79, 6261 (1983)JCPSA6000079000012006261000001
    henceforth referenced as KGT]. Here, we use the expression by Kashchiev [see Eq. (31) of KGT] to estimate the induction time, while for a lithium-disilicate glass of viscosity of ~4.5×109  Pa  s at 750  K, one obtains an induction time of theta=62.6  s (KGT). In our reference aqueous solution, the induction time after which the steady state nucleation (of extremely low rate) is established is only theta=0.33  ns. This ~11 orders of magnitude difference in the induction time originates mainly from an almost 14 orders of magnitude difference in the self-diffusion coefficient.

    W. W. Mullins and R. F. Sekerka, J. Appl. Phys. 34, 323 (1963)JAPIAU000034000002000323000001;, 35, 444 (1964)JAPIAU000035000002000444000001.

    A. Karma and W.-J. Rappel, Phys. Rev. E 60, 3614 (1998).

    K. Vetsigian and N. Goldenfeld, Phys. Rev. E 68, 060601 (2003).


For access to citing articles, you need to log in.


Figures (12) Tables (3)

Access to article objects (figures, tables, multimedia) requires a subscription; log in to view available files.
(Access to supplementary files, where available, is free for this journal.)

Access to article objects (figures, tables, multimedia) requires a subscription; log in to view available files.
(Access to supplementary files, where available, is free for this journal.)



Close
Google Calendar
ADVERTISEMENT

close