[B.Hamada]

Hi all.
I work on the phase change, so to modeling this situation, I have to write an UDF for variable Cp.
My UDF :

PHP Code:

#include "udf.h" 
DEFINE_SPECIFIC_HEAT(my_user_cp, T, Tref, h, yi) 
{ 
    real cp;
        real delta=6;
    real Tmelt=299.5;
    real cps=1413;
    real cpl=2705;
    real lf=232000;

    if (T < Tmelt)
    {
        cp = cps;
    }
    else if ((T>Tmelt) && (T<Tmelt+delta))
    {
        cp = ((cps+cpl)/2)+(lf/delta);
    }
    else (T > Tmelt+delta) 
    {
        cp = cpl;
    }
    *h = cp*(T-Tref);
    return cp; 
} 


But when i interpret my UDF, I found errors on line 20 and 23
Please help me to solve the problem.

[upeksa]

#include "udf.h"
DEFINE_SPECIFIC_HEAT(my_user_cp, T, Tref, h, yi)
{
real cp;
real delta=6;
real Tmelt=299.5;
real cps=1413;
real cpl=2705;
real lf=232000;

if (T < Tmelt)
{
cp = cps;
}
else if ((T>=Tmelt) && (T<=(Tmelt+delta)))
{
cp = ((cps+cpl)/2)+(lf/delta);
}
else
{
cp = cpl;
}
*h = cp*(T-Tref);
return cp;
}

//It works now

[B.Hamada]

Originally Posted by upeksa

#include "udf.h"
DEFINE_SPECIFIC_HEAT(my_user_cp, T, Tref, h, yi)
{
real cp;
real delta=6;
real Tmelt=299.5;
real cps=1413;
real cpl=2705;
real lf=232000;

if (T < Tmelt)
{
cp = cps;
}
else if ((T>=Tmelt) && (T<=(Tmelt+delta)))
{
cp = ((cps+cpl)/2)+(lf/delta);
}
else
{
cp = cpl;
}
*h = cp*(T-Tref);
return cp;
}

//It works now

Thank you it works normally