No abstract

Understanding the thermal and mechanical properties of CH4 and CO2 hydrates is essential for the replacement of CH4 with CO2 in natural hydrate deposits as well as for CO2 sequestration and storage. In this work, we present isothermal compressibility, isobaric thermal expansion coefficient and specific heat capacity of fully occupied single-crystal sI-CH4 hydrates, CO2 hydrates and hydrates of their mixture using molecular dynamics simulations. Eight rigid/nonpolarisable water interaction models and three CH4 and CO2 interaction potentials were selected to examine the atomic interactions in the sI hydrate structure. The TIP4P/2005 water model combined with the DACNIS united-atom CH4 potential and TraPPE CO2 rigid potential were found to be suitable molecular interaction models. Using these molecular models, the results indicate that both the lattice parameters and the compressibility of the sI hydrates agree with those from experimental measurements. The calculated bulk modulus for any mixture ratio of CH4 and CO2 hydrates varies between 8.5 GPa and 10.4 GPa at 271.15 K between 10 and 100 MPa. The calculated thermal expansion and specific heat capacities of CH4 hydrates are also comparable with experimental values above approximately 260 K. The compressibility and expansion coefficient of guest gas mixture hydrates increase with an increasing ratio of CO2-to-CH4, while the bulk modulus and specific heat capacity exhibit the opposite trend. The presented results for the specific heat capacities of 2220-2699.0 J kg(-1) K(-1) for any mixture ratio of CH4 and CO2 hydrates are the first reported so far. These computational results provide a useful database for practical natural gas recovery from CH4 hydrates in deep oceans where CO2 is considered to replace CH4, as well as for phase equilibrium and mechanical stability of gas hydrate-bearing sediments. The computational schemes also provide an appropriate balance between computational accuracy and cost for predicting mechanical and thermal properties of gas hydrates in the high temperature range (≥260 K), and the schemes may be useful for the study of other complex hydrate systems.

In this paper we apply the general analysis described in our first paper to a binary mixture of cyclohexane and n-hexane. We use the square gradient model for the continuous description of a nonequilibrium surface and obtain numerical profiles of various thermodynamic quantities in various stationary state conditions. Details of the numerical procedure are given and discussed. In the second part of this paper we focus on the verification of local equilibrium of the surface as described with excess densities. We give a definition of the temperature and chemical potentials of the surface and verify that these quantities are independent of the choice of the dividing surface. We verify numerically that the surface in a stationary state of the mixture can be described in terms of Gibbs excess densities, which are found to be in good approximation equal to their equilibrium values at the stationary state temperature and chemical potentials of the surface.

We present a theory that describes the transport properties of the interfacial region with respect to heat and mass transfer. Postulating the local Gibbs relation for a continuous description inside the interfacial region, we derive the description of the Gibbs surface in terms of excess densities and fluxes along the surface. We introduce overall interfacial resistances and conductances as the coefficients in the force-flux relations for the Gibbs surface. We derive relations between the local resistivities for the continuous description inside the interfacial region and the overall resistances of the surface for transport between the two phases for a mixture. It is shown that interfacial resistances depend among other things on the enthalpy profile across the interface. Since this variation is substantial, the coupling between heat and mass flow across the surface is also substantial. In particular, the surface puts up much more resistance to the heat and mass transfer than the homogeneous phases over a distance comparable to the thickness of the surface. This is the case not only for the pure heat conduction and diffusion but also for the cross effects such as thermal diffusion. For the excess fluxes along the surface and the corresponding thermodynamic forces, we derive expressions for excess conductances as integrals over the local conductivities along the surface. We also show that the curvature of the surface affects only the overall resistances for transport across the surface and not the excess conductivities along the surface.