Thermal Properties of Supercritical Carbon Dioxide by Monte Carlo Simulations

We present simulation results for the volume expansivity, isothermal compressibility, isobaric heat capacity, Joule-Thomson coefficient and speed of sound for carbon dioxide (CO 2 ) in the supercritical region, using the fluctuation method based on Monte Carlo simulations in the isothermal-isobaric ensemble. We model CO 2 as a quadrupolar two-center Lennard-Jones fluid with potential parameters reported in the literature, derived from vapor-liquid equilibria (VLE) of CO 2 . We compare simulation results with an equation of state (EOS) for the two-center Lennard-Jones plus point quadrupole (2CLJQ) fluid and with a multiparametric EOS adjusted to represent CO 2 experimental data. It is concluded that the VLE-based parameters used to model CO 2 as a quadrupolar two-center Lennard-Jones fluid (both simulations and EOS) can be used with confidence for the prediction of thermodynamic properties, including those of industrial interest such as the speed of sound or Joule-Thomson coefficient, for CO 2 in the supercritical region, except in the extended critical region.

[1]  J. Ilja Siepmann,et al.  Vapor–liquid equilibria of mixtures containing alkanes, carbon dioxide, and nitrogen , 2001 .

[2]  Kwong H. Yung,et al.  Carbon Dioxide's Liquid-Vapor Coexistence Curve And Critical Properties as Predicted by a Simple Molecular Model , 1995 .

[3]  Jie Yang,et al.  COMPASS Force Field for 14 Inorganic Molecules, He, Ne, Ar, Kr, Xe, H2, O2, N2, NO, CO, CO2, NO2, CS2, and SO2, in Liquid Phases , 2000 .

[4]  W. Wagner,et al.  A New Equation of State for Carbon Dioxide Covering the Fluid Region from the Triple‐Point Temperature to 1100 K at Pressures up to 800 MPa , 1996 .

[5]  J. Fischer,et al.  Determination of an effective intermolecular potential for carbon dioxide using vapour-liquid phase equilibria from NpT + test particle simulations , 1994 .

[6]  T. Kristóf,et al.  Heat capacities of two-centre Lennard–Jones fluids from NpT ensemble Monte Carlo simulations , 1999 .

[7]  J. Ilja Siepmann,et al.  Direct Gibbs Ensemble Monte Carlo Simulations for Solid−Vapor Phase Equilibria: Applications to Lennard−Jonesium and Carbon Dioxide , 2001 .

[8]  J. Fischer,et al.  Vapor-liquid equilibria of binary mixtures containing methane, ethane, and carbon dioxide from molecular simulation , 1996 .

[9]  Hans Hasse,et al.  A Set of Molecular Models for Symmetric Quadrupolar Fluids , 2001 .

[10]  A. D. Santis,et al.  The potential model dependence of the neutron radial and partial distribution functions for liquid CO2 , 1987 .

[11]  I. Szalai,et al.  Monte Carlo simulation of some thermophysical properties of two-centre Lennard-Jones fluids along the vapour-liquid equilibrium curve , 1998 .

[12]  K. Shukla,et al.  Effect of quadrupole moment on excess properties of binary 2CLJQ mixtures from molecular dynamics simulation: model and real mixtures , 1997 .

[13]  Antoinette J. Taylor,et al.  Thermodynamic and structural properties of liquids modelled by ‘2-Lennard-Jones centres’ pair potentials , 1977 .

[14]  H. Hasse,et al.  Comprehensive study of the vapour–liquid equilibria of the pure two-centre Lennard–Jones plus pointquadrupole fluid , 2001, 0904.3413.

[15]  K. Gubbins,et al.  Accurate CO2 Joule-Thomson inversion curve by molecular simulations , 2002 .

[16]  Johann Fischer,et al.  An equation of state for two-center Lennard-Jones fluids , 1997 .

[17]  I. R. Mcdonald,et al.  Electrostatic interactions in molecular crystals , 1983 .

[18]  Jeffrey J. Potoff,et al.  Molecular simulation of phase equilibria for mixtures of polar and non-polar components , 1999 .

[19]  J. Fischer,et al.  Construction and application of physically based equations of state: Part II. The dipolar and quadrupolar contributions to the Helmholtz energy , 1992 .

[20]  Philippe Ungerer,et al.  Prediction of thermodynamic derivative properties of fluids by Monte Carlo simulation , 2001 .