|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
关于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;
}
不知错误到底在哪里,结果根本就不对,很是着急,希望大家能帮帮忙?非常非常感谢! |
|