Home
Research
Molecular Simulations
Fluid Inclusions
Models
Members
Publications
Programs
Gallery
Links

       
  Thermodynamic models

   
         
   
Files
Descriptions
Ref
EOS_CO2_H2O.exe
d_CO2_H2O_NaCl.exe

The program (EOS_CO2_H2O.exe) is used to calculate: (1)-compressibility(Z), (2)-molar volume(V) and (3)-pressure(P) of the CO2-H2O system up to 647 K and 100 MPa. The program (d_CO2_H2O_NaCl.exe) is used to calculate: (1)-desnity(d), (2)-molar volume(V), (3)-apparent molar volume(Vfai) and (4)-apparent molar volume at infinite dilution(Vfai) of the CO2+H2O and CO2+H2O+NaCl systems up to 647 K and 100 MPa.
Zhenhao Duan, Jiawen Hu, Dedong Li and Shide Mao (2008) Density of the CO2-H2O and CO2-H2O-NaCl systems up to 647 K and 100 MPa. Energy & Fuels
viscosity.exe

An accurate viscosity model over a wide range of temperature, pressure and concentration is developed for the LiCl-H2O(273-573 K, 1-400 bar and 0-16.7 m), NaCl-H2O(273-573 K, 1-500 bar and 0-6.0 m) and KCl-H2O(273-543 K, 1-500 bar and 0-4.5 m) systems with the experimental accuracy.
Shide Mao & Zhenhao Duan (2008), The viscosity of aqueous alkali-chloride solutions up to high temperature, pressure and ionic strength, Submitted to Journal of Chemical & Engineering Data.
density_model_on_net.exe
isochore_NaClH2O.exe

An accurate and simple density model over a wide range of temperature, pressure and concentration is developed for binary LiCl-H2O, NaCl-H2O, KCl-H2O, MgCl2-H2O, CaCl2-H2O, SrCl2-H2O and BaCl2-H2O systems within or about the experimental uncertainties. The infinite dilute apparent molar volume, density and isochores of the fluid inclusions can be calculated from this density model.
Shide Mao & Zhenhao Duan (2008) The PVT properties of aqueous chloride fluids up to high temperatures and pressures. Submitted to The Journal of Chemical Thermodynamics
EOS_C_H.exe

This program is used to calculate volumetric properties and activity of the H2O-CH4 673K to 2573K up to 10 Gpa.
Chi Zhang, Zhenhao Duan* and Zhigang Zhang (2006) Molecular simulation of the CH4 and CH4-H2O systems with ab initio potential surface up to 10GPa and 2573K
h2o_cacl2.exe
h2o_mgcl2.exe

These two programs are used to calculate volumetric properties and liquid-vapor phase euiqlibria of the H2O-CaCl2 and H2O-MgCl2 systems from 523K to 973K and from water satuartion pressure up to 1500bar.
Zhenhao Duan, John Weare and Nancy Moller (2006) A high temperature equation of state for the H2O-CaCl2 and H2O-MgCl2 systems. Geochimica Cosmochimica Acta, 70 (15): 3765-3777.[Pdf]
A thermodynamic model is presented to calculate N2 solubility in pure water (273-590 K and 1-600 bar) and aqueous NaCl solutions (273-400 K, 1-600 bar and 0-6 mol kg-1), gas phase composition of the H2O-N2-NaCl system (273-523 K and 1-300 bar), the N2 solubility in seawater, and liquid phase density of H2O-N2-NaCl system at phase equilibrium.
Shide Mao & Zhenhao Duan (2006) A thermodynamic model for calculating nitrogen solubility, gas phase composition and density of the H2O-N2-NaCl system. Fluid Phase Equilibria., 248 (2): 103-114.[Pdf]
EOS for the system of H2O, CO2 and CO2-H2O, which was proposed on the basis of MD simulation results (click here to download) and experimental data. This EOS is expected to show improved accuracy for the PVTx properties of CO2-H2O from 673.15 K up to 2573.15 K and 10.0 GPa. The prediction power of this EOS can be found here.
Zhenhao Duan & Zhigang Zhang (2006) Equation of state of the H2O, CO2, and H2O-CO2 systems up to 10 GPa and 2573.15 K: Molecular dynamics simulations with ab initio potential surface. Geochim. Cosmochim. Acta. , 70 (9): 2311-2324.[Pdf]
This model can calculate water content of gas phase for CH4-H2O-NaCl system; CH4 solubility in pure water and aqueous NaCl solutions ; CH4 solution heat(deltaH) and partial volume(VlCH4) in pure water and aqueous NaCl solutions;The liquid density for CH4-H2O-NaCl system at equilibria and CH4 solubility in seawater and brines. The T-P-m range is: 273-523 K, 1-2000 bar and 0-6 m.
Zhenhao Duan & Shide Mao (2006) A thermodynamic model for calculating methane solubility, density and gas phase composition of methane-bearing aqueous fluids from 273 to 523 K and from 1 to 2000 bar. Geochim. Cosmochim. Acta., 70 (13): 3369-3386.[Pdf]
This program is to calculate ethane solubilities in water and aqueous NaCl solutions. T-P range for ethane-water system is: 273-444 K; 0-1000 bar. T-P-m range for ethane-water-NaCl system is: 273-348 K; 0-16 bar; 0-6.3 m.
Shide Mao, et al. (2005) An accurate model for calculating C2H6 solubility in water and aqueous NaCl solutions. Fluid Phase Equilibria, 238: 77-86 [Pdf]
Programs for the prediction of CH4 and CO2 hydrate phase equilibrium and cage occupancy. All the relevant files (codes and binaries) are included in the .zip compressed file. Separate files can also be downloaded one by one. Instructions can be found in Readme.txt
Rui Sun and Zhenhao Duan (2005) Prediction of CH4 and CO2 hydrate phase equilibrium and cage occupancy from Ab Initio intermolecular potential. Geochim. Cosmochim. Acta, 69: 4411-4424. [Pdf]
Programs for the prediction of CH4 and CO2 hydrate phase equilibrium in aqueous NaCl solution, in seawater, and in brines containing Na+, K+, Ca2+, Mg2+, Cl-, and SO42-.
Zhenhao Duan and Rui Sun, A Model to Predict Phase Equilibrium of CH4 and CO2 Clathrate Hydrate in Aqueous Electrolyte Solutions. American Mineralogist., 91 (8-9): 1346-1354.[Pdf]
Source code in PASCAL and descriptions can be found here. [Duan_Hu_Eos.rar]
Duan Zhenhao, Hu Jiawen, (2004) A new cubic equation of state and its applications to the modeling of vapor-liquid equilibria and volumetric properties of natural fluids. Geochim. Cosmochim. Acta, 68 (14): 2997-3009 [Pdf]
This program calculates the solubility of CO2 in pure water and in aqueous NaCl solution. T-P-X range of this model: 273-533 K, 0-2000 bar, 0-4.5 m NaCl.

Zhenhao Duan and Rui Sun (2003) An improved model calculating CO2 solubility in pure water and aqueous NaCl solutions from 273 to 533 K and from 0 to 2000bar. Chem. Geol., 193: 253-271 [Pdf]

Duan ZH, Sun R, Zhu C, et al. (2006) An improved model for the calculation of CO2 solubility in aqueous solutions containing Na+, K+, Ca2+, Mg2+, Cl-, and SO42-. Marine Chem., 98 (2-4): 131-139. [Pdf]

This program calculates molar volume and the composition of CH4-H2O system at phase equilibrium. T-P range of this model: 273-383 K, 0 -1000 bar.
Rui Sun and Zhenhao Duan (2003) A new equation of state and Fortran 77 program to calculate vapor-liquid phase equilibria of CH4-H2O system at low temperatures. Computers & Geosciences 29: 1291-1299 [Pdf]
EOS for pure water, which was proposed on the basis of MD simulation results and experimental data. This EOS is expected to show improved accuracy for the volumetric properties of liquid water from 100.0 MPa and 273.15 K up to substantial ionization limit (about 20 GPa and 2000 K)
Zhigang Zhang and Zhenhao Duan. (2005) Isothermal-isobaric molecular dynamics simulations of the PVT properties of water over wide range of temperatures and pressures. Phys. Earth Planet Interiors, 149: 335-354 [Pdf]
 
         
  Molecular Simulation Packages


   
   
Files
Descriptions
Ref
Isothermal-Isobaric molecular dynamics simulation program for versatile purposes. (Readme, Examples)
Zhigang Zhang and Zhenhao Duan. (2005) Isothermal-isobaric molecular dynamics simulations of the PVT properties of water over wide range of temperatures and pressures. Phys. Earth Planet Interiors, 149: 335-354 [Pdf]
Gibbs Ensemble Monte Carlo simulation for spherical models (CH4-C2H6) (Examples)
Zhigang Zhang and Zhenhao Duan. (2002) Phase equilibria of the system methane-ethane from temperature scaling Gibbs ensemble Monte Carlo simulation. Geochim. Cosmochim. Acta. 66 (19): 3431-3439 [Pdf]
 
         
   
Caution: the programs listed above are provided "as is" - their correctness and performance for any particular application cannot be guaranteed.
 
         

 

WWW www.geochem-model.org

Zhenhao Duan Research Group
State Key Laboratory of Lithospheric Evolution
Institute of Geology and Geophysics
Chinese Academy of Sciences
Beijing, China, 100029
Tel: 86-10-82998377
Fax: 86-10-82024251
E-mail: duanzhenhao@yahoo.com