[pavankonchada]

I am trying to write udf for variable dynamic viscosity for PCM as shown below
DEFINE_PROPERTY(PCM_Dynamicviscosity,cell,thread)
{
real tempe2,mu,bt;
tempe2 = C_T(cell,thread);
bt=(6722811.795)-(67319.37*tempe2)+(224.7*tempe2*tempe2)-(0.25*tempe2*tempe2*tempe2);
printf("bt=%f \n",bt);
if (tempe2 < 298.6)
mu=180000;
if ((tempe2 >= 298.6) && (tempe2 <= 300.6))
mu=0.001798(1+((100000(1-bt)*(1-bt))/((bt*bt*bt)+0.001)));
else
{
mu=0.001798;
}
return mu;
}
the value of bt is 0.011 and hence mu value is 175632.39
but the value I am getting for bt when initializing is 1.5 and dynamic viscosity is coming wrong
Problem is
0.25*tempe2*tempe2*tempe2=6655940.5 (fluent value)
0.25*tempe2*tempe2*tempe2=6655941.9189(Actual value)