表面张力添加

wenqi73

做的是焊接模拟，只建了一个方块模型，想在上表面添加表面张力，我是在fluent中壁面边界条件中的shear stress添加的.读入udf，再打开VOF之后会报错；关闭VOF的话可以读取，但是计算式x动量方向不收敛！，udf如下，请大神求解。  

DEFINE_PROFILE(surface_tension_g,t,i)
{cell_t c;

begin_c_loop(c,t)
     C_UDMI(c,t,0)=C_T_G(c,t)[0];
end_c_loop(c,t)     
return 0;
}

DEFINE_PROFILE(surface_tension,t,i)
{
  cell_t c;

   begin_c_loop(c,t)               
     {
           real temp=C_T(c,t);
       real Kseg=k1*exp(-dH/R0/temp);
       real beta=1+Kseg*a;
       real sur_co,sur_f;
             
       sur_co=-Ar-R0*gamma*log(beta)-Kseg*gamma*dH/beta/temp;
       F_PROFILE(c,t,i)=sur_co*C_UDMI(c,t,0);

lee_corn

我也是从论坛里下的，可以看看