DEFINE_SPECIFIC_HEAT (cp, T, Tref, h, yi)
{
real cp;
if (T >=540 && T<650)
cp=4.12126E-5*pow(T,4)-0.0947*pow(T,3)+81.65593*pow(T,2)-31302.17859*T+4.5041E6;
else if (T >=650 && T<688)
cp=-0.02116*pow(T,4)+56.4407*pow(T,3)-56449.66828*pow(T,2)+2.50871E7*T-4.18E9;
else if (T >=688 && T<750)
cp=-0.00139*pow(T,4)+3.98394*pow(T,3)-4266.71572*pow(T,2)+2.0283E6*T-3.61059E8;
else
cp=1.97284E-6*pow(T,4)-0.0073*pow(T,3)+10.13277*pow(T,2)-6260.7141*T+1.45669E6;
return cp;
}
|