[Cooper24]

Hi. I have written a UDF for my source terms in energy equation. The part of the UDF inside DEFINE_SOURCE is:


if (r2 < R2 && myid == C_PART(c,sec_th)) // line 1
{
l = round((x[0] - sp[0] + R)/dz);
m = round((x[1] - sp[1] + R)/dz);
source = 1.0;
dS[eqn] = 0.0;
C_UDMI(c,t,7) = 1.0;
C_UDMI(c,t,8) = source;
}
else
{
source = 0.0;
C_UDMI(c,t,8) = source;
C_UDMI(c,t,7) = 0.0;
dS[eqn] = 0.0;
}
return source;


When I am using the if condition of line 1 and using C_UDMI(c,t,7) and C_UDMI(c,t,8) values in CFD-Post for plotting, it works fine. The values of these two UDMIs go from 0 to 1.


But when I am using the following UDF, the maximum values of UDMIs 7 and 8 is reduced by half in CFD-Post. They go from 0 to 0.5. I am unable to understand what is the issue.



if (r2 < R2 && myid == C_PART(c,sec_th)) // line 1
{
l = round((x[0] - sp[0] + R)/dz);
m = round((x[1] - sp[1] + R)/dz);

if( x[2] - (sp[2] - Hlaser[l][m]) <= 0.05*dz && x[2] - (sp[2] - Hlaser[l][m]) >= -0.1*dz) // line 2
{
source = 1.0;
dS[eqn] = 0.0;
C_UDMI(c,t,7) = 1.0;
C_UDMI(c,t,8) = source;
}
else
{
source = 0.0;
C_UDMI(c,t,8) = source;
C_UDMI(c,t,7) = 0.0;
dS[eqn] = 0.0;
}
}
else
{
source = 0.0;
C_UDMI(c,t,8) = source;
C_UDMI(c,t,7) = 0.0;
dS[eqn] = 0.0;
}
return source;


Can someone please help!! Thanks in advance.