[Bamse]

Hi,

Im trying to create a UDF for the specific heat being dependent on the temperature in a solid material. As you can see the we have one equation for T >= 973.15 and another for below that. I have complied this UDF and haven't gotten any error. The simulation I'm making is a transient one. The question I have though is the calculation of "*h = cp_all*(T-Tref)" correct? Should this be an integration or not?
Do I even need to have "*h = cp_all*(T-Tref)" included in the UDF?
Another question that I have is what should I set the Tref to? I didn't see any guidelines on this in the UDF manual.

Thanks,

#include "udf.h"
DEFINE_SPECIFIC_HEAT(my_user_cp, T, Tref, h, yi)
{
real cp_all;
if (T >= 973.15)
{
cp_all = 0.05 * T + 102.72;
*h = cp_all*(T-Tref);
}
else
{	
cp_all = -0.3547 * T + 254.;
*h = cp_all*(T-Tref);
}
return cp_all;
}