FLUENT中的多相流问题

ljzzhet

最近在研究风作用下的雪漂移运动，看文献中说采用多相流模型进行求解，其中空气相采用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 编辑 ]

lijiajun_Euler

你好，我想请问一下，这几句语句是什么作用：
if(F_C1_THREAD(f,t)!=NULL)
     {
      t1=F_C1_THREAD(f,t);
      c1=F_C1(f,t);
     }
   else
     {
      t1=NULL;
      c1=-1;
     }