public static double CalculateWaterSaturationPressure(double temperature)
{
double pSaturation = 1.0;
//if (temperature >= 253.15 && temperature <= 283.15) {
// pSaturation = 8.8365E-10 * Math.Pow((temperature -193.03), 6.2145);
//}
//else if (temperature <= 453.15) {
// pSaturation = Math.Exp(23.197 - 3816.44/(temperature - 46.13));
//}
if (temperature <= 273.33) //temperature belown which steam table can not deal
//double t = temperature - 273.15;
//pSaturation = 100.0 * 6.1121 * Math.Exp(17.502 * t/(240.97 + t));
//To Do: When system is not air-water, need to change!!!
{
pSaturation = ThermalPropCalculator.CalculateWaterSaturationPressureBelowFreezingPoint(temperature);
}
else
{
//formulation coming from Perry's Chemical Engineer's Handbook
//pSaturation = Math.Exp(73.649 - 7258.2/temperature - 7.3037 * Math.Log(temperature)
// + 4.1653e-6 * temperature * temperature);
//pSaturation = ThermalPropCalculator.CalculateVaporPressure(temperature, vapPressureCoeffs);
SteamTable steamTable = SteamTable.GetInstance();
pSaturation = steamTable.GetSaturationPressure(temperature);
}
return(pSaturation);
}