[mabdolahim]

Hi everyone

I am supposed to model a material that its phase changes with temperature; its thermal conductivity behaves as a function of temperature for three temperature ranges. I write a UDF but it is not working as expected. Can anyone help me with this problem? the functions are attached and there is my code:


#include "udf.h"

#define A 0.23152914
#define B 2.4105618
#define C 0.0014546952
#define D 0.21946513
#define E 0.0000168736358
#define F 0.005086847

DEFINE_PROPERTY(thermal_npcm,c,t)
{
real k;
real temp = C_T(c,t);
if(temp<302.)
k = A;
else if (temp>309.)
k = A+E*pow(temp,1.5)-F*pow(temp,0.5);
else
k = A+(B*pow(temp,2.5))-C*pow(temp,1.5)+D*pow(temp,0.5);

return k;
}

[Cooper24]

What is the error you are receiving??

[mabdolahim]

Originally Posted by Cooper24

What is the error you are receiving??

Dear Cooper24,
The solving proceeds without any error but the problem is the outputs that are so wrong. By constant thermal conductivity, the material melts very well, while by UDF thermal conductivity the material state remains unchanged.