//---------------------------------------------------------------------------
//  Project:         FlowLive
//  Module:          thermo
/*!\file             acentric_factor.h
	\brief            acentric factors
	\version          0.0.2.1
	\date             Created: 2007-09-21
	\date             Last modified: 2007-09-21
	\author           <a href="mailto:xtobias@gmx.net">Tobias Severin, xtobias@gmx.net</a>
*/
//---------------------------------------------------------------------------
#ifndef ACENTRIC_FACTOR_H
#define ACENTRIC_FACTOR_H
//---------------------------------------------------------------------------
#include <thermo/thermo.h>

#include <string>

//---------------------------------------------------------------------------
THERMO_BEGIN

// data from:
// http://de.wikipedia.org/wiki/Azentrischer_Faktor

// Helium-3 (He3)
static const double kAcentricFactor_He3 = -0.473;

// Helium (He)
static const double kAcentricFactor_He = -0.365;

// Argon (Ar)
static const double kAcentricFactor_Ar = 0.0001;

// Xenon (Xe)
static const double kAcentricFactor_Xe = 0.008;

// Wasserstoff (H2)
static const double kAcentricFactor_H2 = -0.216;

// Stickstoff (N2)
static const double kAcentricFactor_N2 = 0.039;

// Sauerstoff (O2)
static const double kAcentricFactor_O2 = 0.025;

// Fluor (F2)
static const double kAcentricFactor_F2 = 0.054;

// Chlor (Cl2)
static const double kAcentricFactor_Cl2 = 0.090;

// Brom (Br2)
static const double kAcentricFactor_Br2 = 0.108;

// Fluorwasserstoff (HF)
static const double kAcentricFactor_HF = 0.329;

// Wasser (H2O)
static const double kAcentricFactor_H2O = 0.344;

// Schweres Wasser (D2O)
static const double kAcentricFactor_D2O = 0.351;

// Ammoniak (NH3)
static const double kAcentricFactor_NH3 = 0.250;

// Methan (CH4)
static const double kAcentricFactor_CH4 = 0.011;

// Ethylen (C2H4)
static const double kAcentricFactor_C2H4 = 0.089;

// Propan (C3H8)
static const double kAcentricFactor_C3H8 = 0.153;

// n-Butan (C4H10)
static const double kAcentricFactor_C4H10 = 0.199;

// Isobutan (C4H10)
static const double kAcentricFactor_i_C4H10 = 0.183;

// Kohlenmonoxid (CO)
static const double kAcentricFactor_CO = 0.066;

// Kohlendioxid (CO2)
static const double kAcentricFactor_CO2 = 0.239;

// Tetrafluormethan (CF4)
static const double kAcentricFactor_CF4 = 0.177;

// Tetrachlorkohlenstoff (CCl4)
static const double kAcentricFactor_CCl4 = 0.193;

// Benzol (C6H6)
//static const double kAcentricFactor_C6H6 = 0.212;

// Toluol (C7H8)
static const double kAcentricFactor_C7H8 = 0.263;

// Methanol (CH4O)
static const double kAcentricFactor_CH3OH = 0.556;

// Ethanol (C2H6O)
static const double kAcentricFactor_C2H5OH = 0.644;

// Aceton (C3H6O)
static const double kAcentricFactor_CH3COCH3 = 0.304;

// Essigs鋟re (C2H4O2 / CH3COOH)
static const double kAcentricFactor_CH3COOH = 0.447;

// Schwefeldioxid (SO2)
static const double kAcentricFactor_SO2 = 0.256;

// Schwefeltrioxid (SO3)
static const double kAcentricFactor_SO3 = 0.481;

// Quecksilber (Hg)
static const double kAcentricFactor_Hg = -0.167;



/*
Data from:
A GENERALIZED VAPOR PRESSURE
EQUATION FOR HEAVY
HYDROCARBONS
Chorng H. Twu*, John E. Coon, and John R.
Cunningham
Simulation Sciences, Inc.
601 South Valencia Avenue
Brea, CA 92621
*/


// Argon
//static const double kAcentricFactor_Ar = 0.0012; // method 1
//static const double kAcentricFactor_Ar = 0.0015; // method 2

// Methane
//static const double kAcentricFactor_CH4 = 0.0133;
// 0.0136

// Ethane
static const double kAcentricFactor_C2H6 = 0.0992;
// 0.0994

// Propane
//static const double kAcentricFactor_C3H8 = 0.1534;
// 0.1536

// n-Butane
//static const double kAcentricFactor_C4H10 = 0.1993;
// 0.1994

// n-Pentane
static const double kAcentricFactor_C5H12 = 0.2509;
// 0.2509

// n-Hexane
static const double kAcentricFactor_C6H14 = 0.2980;
// 0.2980

// n-Heptane
static const double kAcentricFactor_C7H16 = 0.3487;
// 0.3487

// n-Octane
static const double kAcentricFactor_C8H18 = 0.3977;
// 0.3977

// n-Nonane
static const double kAcentricFactor_C9H20 = 0.4396;
// 0.4396

// n-Decane
static const double kAcentricFactor_C10H22 = 0.4872;
// 0.4872

// n-Undecane
static const double kAcentricFactor_C11H24 = 0.5376;
// 0.5377

// n-Dodecane
static const double kAcentricFactor_C12H26 = 0.5765;
// 0.5766

// n-Tridecane
static const double kAcentricFactor_C13H28 = 0.6250;
// 0.6251

// n-Tetradecan
static const double kAcentricFactor_C14H30 = 0.6659;
// 0.6661

// n-Pentadecan
static const double kAcentricFactor_C15H32 = 0.7126;
// 0.7128

// n-Hexadecane
static const double kAcentricFactor_C16H34 = 0.7480;
// 0.7483

// n-Heptadecane
static const double kAcentricFactor_C17H36 = 0.7721;
// 0.7483

// n-Octadecane
static const double kAcentricFactor_C18H38 = 0.7961;
// 0.7964

// n-Nonadecane
static const double kAcentricFactor_C19H40 = 0.8075;
// 0.8079


// extrapolated values
static const double kAcentricFactor_C20H42 = 0.8367;
static const double kAcentricFactor_C21H44 = 0.8614;
static const double kAcentricFactor_C22H46 = 0.8843;
static const double kAcentricFactor_C23H48 = 0.9054;
static const double kAcentricFactor_C24H50 = 0.9247;
static const double kAcentricFactor_C25H52 = 0.9422;
static const double kAcentricFactor_C26H54 = 0.9579;



// iso-Butane
//static const double kAcentricFactor_i_C4H10 = 0.1832;
// 0.1833

// iso-Pentane
static const double kAcentricFactor_i_C5H12 = 0.2266;
// 0.2266

// 2-m-Pentane
static const double kAcentricFactor_2m_C5H12 = 0.2776;
// 0.2776

// Ethylene
//static const double kAcentricFactor_C2H4 = 0.0897;
// 0.0899

// Propylene
static const double kAcentricFactor_C3H6 = 0.1419;
// 0.1420

// 1-Butene
static const double kAcentricFactor_C4H8 = 0.1900;
// 0.1901

// Cyclopentane
static const double kAcentricFactor_cy_C5H10 = 0.1945;
// 0.1946

// m-Cyclopentane
//static const double kAcentricFactor_??? 0.2304;
// 0.2305

// Cyclohexane
static const double kAcentricFactor_cy_C6H12 = 0.2108;
// 0.2109

// Benzene
static const double kAcentricFactor_C6H6 = 0.2111;
// 0.2112

// Toluene
//static const double kAcentricFactor_??? 0.2624;
// 0.2625

//---------------------------------------------------------------------------
//! Returns the acentric factor of a given specie
THERMO_EXPORT double AcentricFactor(const std::string &aSpecieName);


//---------------------------------------------------------------------------
THERMO_END

#endif