（菜鸟求助）关于质量源项问题 请问可以用C_UDSI（c,t,i）和C_UDSI_G（c,t,i）吗

dadadada

先在这里谢谢了
请问这里面的i是什么意思 可以取1、2、3......吗
这两个宏在什么情况下可以用啊
我最近编辑一个关于质量源项的UDF  不知道可不可以这么用
我加了四个scalar

#include "udf.h"


DEFINE_SOURCE(g_source, c, t, dS, eqn)
{
Domain*subdomain_l=Get_Domain(2);
Domain*subdomain_g=Get_Domain(3);
Thread*tl=Lookup_Thread(subdomain_l,2);
Thread*tg=Lookup_Thread(subdomain_g,2);


real mu_t_water;
real density_water=998.2,density_gas=1.225;
real t_prt=0.67;
real x[ND_ND];
real y[ND_ND];
real a,mu_t_gas,A,B,source;

mu_t_water=C_MU_T(c,tl);
a=C_VOF(c,tg);
mu_t_gas=mu_t_water*density_gas/density_water;
C_UDSI(c,t,2)=a;
NV_V(x,=,C_UDSI_G(c,t,2));
A=ND_SUM(x[0],x[1],x[2]);
C_UDSI(c,t,3)=A;
NV_V(y,=,C_UDSI_G(c,t,3));
B=ND_SUM(y[0],y[1],y[2]);
source=a*mu_t_gas*B/t_prt;
  return source;
}



DEFINE_SOURCE(l_source, c, t, dS, eqn)
{
Domain*subdomain_l=Get_Domain(2);
Thread*tl=Lookup_Thread(subdomain_l,2);

real mu_t_water;
real t_prt=1;
real m[ND_ND];
real n[ND_ND];
real a,A,B,source;

mu_t_water=C_MU_T(c,tl);
a=C_VOF(c,tl);
C_UDSI(c,t,0)=a;
NV_V(m,=,C_UDSI_G(c,t,0));
A=ND_SUM(m[0],m[1],m[2]);
C_UDSI(c,t,1)=A;
NV_V(n,=,C_UDSI_G(c,t,1));
B=ND_SUM(n[0],n[1],n[2]);
source=a*mu_t_water*B/t_prt;

  return source;
}

[ 本帖最后由 dadadada 于 2009-4-24 15:37 编辑 ]

dadadada

高手们 急啊 周围也没人可问 比如气含率的梯度该怎么求啊