//cp unit is J/kg.k
private static double CalculateCp(ArrayList materialComponents, double temperature, SubstanceState state)
{
double heatCapacity = 0.0;
double cp = 0;
double molarWeight;
double massFrac;
foreach (MaterialComponent mc in materialComponents)
{
Substance s = mc.Substance;
molarWeight = s.MolarWeight;
ThermalPropsAndCoeffs tpc = s.ThermalPropsAndCoeffs;
massFrac = mc.GetMassFractionValue();
if (state == SubstanceState.Gas)
{
cp = ThermalPropCalculator.CalculateGasHeatCapacity1(temperature, tpc.GasCpCoeffs) / molarWeight;
}
else if (state == SubstanceState.Liquid)
{
cp = ThermalPropCalculator.CalculateLiquidHeatCapacity1(temperature, tpc.LiqCpCoeffs) / molarWeight;
}
heatCapacity += cp * massFrac;
}
return(heatCapacity);
}
public double GetSpecificHeatOfLiquid(double temperature)
{
//this correlation comes from Perry's Chemical Engineers Handbook
//double cp = 2.7637e5 - 2.0901e3 * temperature + 8.125 * temperature * temperature
// -1.4116e-2 * Math.Pow(temperature, 3) + 9.3701e-6 * Math.Pow(temperature, 4);
double cp = ThermalPropCalculator.CalculateLiquidHeatCapacity1(temperature, liqCpCoeffs);
return(cp / molarWeight);
}