表面张力添加
做的是焊接模拟，只建了一个方块模型，想在上表面添加表面张力，我是在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);
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); 我也是从论坛里下的，可以看看
页:
[1]