颗粒尺寸设置的UDF问题

tchllc

我现在用欧拉多相流计算气-液两相流,考虑在入口边界的气泡大小为5MM,在其他流体部分为2MM,下面是我的UDF程序,计算显示入口边界的气泡大小也2MM,请帮忙看什么地方有问题.
DEFINE_PROPERTY(bubble_diameter,f,t)
{
     real d;
    begin_f_loop(f,t)
{
if(THREAD_ID(t)==4)/*inlet ID=4*/
     d=C_PHASE_DIAMETER(f,t)=0.005;
else
         d=C_PHASE_DIAMETER(f,t)=0.002;
}
end_f_loop(f,t)
return d;  
}