关于VOF和熔化/凝固模型的问题

liyanzz

关于VOF和熔化/凝固模型的问题

我做的课题是焊接熔池流动模拟，大概框架为：建立一个长100mm，宽50mm，高4mm的容器，容器壁的材料为钢，使用VOF模型和熔化/凝固模型，定义空气为第一相，钢为第二相，将容器底部高为3mm的区域定义为钢，然后追踪钢和空气的交界面，即找钢的体积分数在0~1之间的网格，然后把高斯热源以能量源相的形式加到追踪到的界面上。我编的源相程序如下：
#include"udf.h"
#define PI 3.1415926
DEFINE_SOURCE(cell_heat_source,cell,mixture_thread,dS,eqn)
{
real x[ND_ND];
real m,n,f,source;
Thread *pri_th,*sec_th;
pri_th = THREAD_SUB_THREAD(mixture_thread, 0);
sec_th = THREAD_SUB_THREAD(mixture_thread, 1);
if(C_VOF(cell,sec_th)<1&&C_VOF(cell,sec_th)>0)
{
m=x[0];
n=x[1];
  f=x[2];
    source=0.65*16*110*6*sqrt(3)/(PI*sqrt(PI)*0.003*0.003*0.004)*exp(-(3*m*m)/(0.003*0.003)-(3*n*n)/(0.003*0.003)-(3*f*f)/(0.004*0.004));
    dS[eqn]=0;
}
else
    source=dS[eqn]=0.;
return source;
}
不知错误到底在哪里，结果根本就不对，很是着急，希望大家能帮帮忙？非常非常感谢！