FLUENT中的多相流问题
最近在研究风作用下的雪漂移运动,看文献中说采用多相流模型进行求解,其中空气相采用FLUENT中默认的控制方程即质量守恒、N-S方程和湍流方程进行求解,而雪相控制方程如附件图:我想请问大家,求解雪相的控制方程是否需要将雪相方程采用UDF的方式编入到FLUENT中呢?
我试着编了下雪相的对流项,部分程序如下:
DEFINE_UDS_FLUX(snow_uds_flux,f,t,i)
{
Thread*t0,*t1=NULL;
cell_t c0,c1=-1;
real NV_VEC(snowfraction_vec),NV_VEC(A);
t0=F_C0_THREAD(f,t);
c0=F_C0(f,t);
if(F_C1_THREAD(f,t)!=NULL)
{
t1=F_C1_THREAD(f,t);
c1=F_C1(f,t);
}
else
{
t1=NULL;
c1=-1;
}
if(t1==NULL)
{
NV_D(snowfraction_vec,=,F_U(f,t),F_V(f,t),F_W(f,t));
NV_S(snowfraction_vec,*=,F_R(f,t));
}
}
其中: NV_D(snowfraction_vec,=,F_U(f,t),F_V(f,t),F_W(f,t));
NV_S(snowfraction_vec,*=,F_R(f,t));如何判断所取的速度和密度是空气相还是雪相呢?新手上路,请大家多多指教!!!
[ 本帖最后由 ljzzhet 于 2012-12-10 14:56 编辑 ] 你好,我想请问一下,这几句语句是什么作用:
if(F_C1_THREAD(f,t)!=NULL)
{
t1=F_C1_THREAD(f,t);
c1=F_C1(f,t);
}
else
{
t1=NULL;
c1=-1;
}
页:
[1]