请教DEFINE_UDS_FLUX

momotto

我想第一个和第二个uds的对流通量为：密度×速度，而其他uds的对流通量为0。
下面这个（帮助文件里的）的udf哪里出错了？怎么用不了？
知道的指点我一下，谢谢了。
DEFINE_UDS_FLUX(UDS_flux,f,t,i)
{
Thread *t0,*t1=NULL;
cell_t c0,c1=-1;
real NV_VEC(psi_vec),NV_VEC(A);
t0=F_C0_THREAD(f,t);
c0=F_C0(f,t);
if(NULL!=F_C1_THREAD(f,t))
{
t1=F_C1_THREAD(f,t);
c1=F_C1(f,t);
}
else
{
t1=NULL;
c1=-1;
}
if(NULL==t1)
{
NV_D(psi_vec, =,F_U(f,t),F_V(f,t),F_W(f,t));
NV_S(psi_vec, *=,F_R(f,t));
}
else
{
NV_D(psi_vec, =,C_U(c0,t0),C_V(c0,t0),C_W(c0,t0));
NV_D(psi_vec,+=,C_U(c1,t1),C_V(c1,t1),C_W(c1,t1));
NV_S(psi_vec,/=,2.);
NV_S(psi_vec,*=,(C_R(c0,t0)+C_R(c1,t1))/2.);
}
F_AREA(A,f,t);
if(i==0)
    return NV_DOT(psi_vec,A);
else if(i==1)
return NV_DOT(psi_vec,A);
else
return 0;
}

gallongallon

没有什么错误，使用有什么问题？

weiyj83

我也是遇到同样的问题，期待解答
/******************************************************************/
DEFINE_UDS_FLUX(sed_convective, f, t, i)
{
  Thread *t0, *t1 = NULL;
  cell_t  c0,  c1 = -1;
  real NV_VEC(psi_vec), NV_VEC(A);
  real ws;
  
  ws = 10 * nu / s * (sqrt(1 + 0.01 * D_non * D_non * D_non) - 1);
  t0 = THREAD_T0(f,t);
  c0 = F_C0(f,t);
  if (NULL != THREAD_T1(f,t))
    {
      t1 = THREAD_T1(f,t);
      c1 = F_C1(f,t);
    }
  else
    {
      t1 = NULL;
      c1 = -1;
    }
  if (NULL == t1)
    {
      NV_D(psi_vec,  =, F_U(f,t), F_V(f,t), (F_W(f,t) - ws));
      NV_S(psi_vec, *=, F_R(f,t));
    }
  else
    {
      NV_D(psi_vec,  =, C_U(c0,t0), C_V(c0,t0), (C_W(c0,t0) - ws));
      NV_D(psi_vec, +=, C_U(c1,t1), C_V(c1,t1), (C_W(c1,t1) - ws));
      NV_S(psi_vec, /=, 2.);     /* averaging. */
      NV_S(psi_vec, *=, (((C_R(c0,t0) + C_R(c1,t1)) / 2.)));
    }
  F_AREA(A, f, t);
  return NV_DOT(psi_vec, A);
}
/******************************************************************/

cpp -ID:\Fluent.Inc\fluent6.1.22/src -ID:\Fluent.Inc\fluent6.1.22/cortex/src -ID:\Fluent.Inc\fluent6.1.22/client/src -ID:\Fluent.Inc\fluent6.1.22/multiport/src -I. -DUDFCONFIG_H="<udfconfig.h>" D:\Fluent.Inc
tbin
tx86\pile01\pile.c
D:\Fluent.Inc\ntbin\ntx86\pile01\pile.c:137: macro `THREAD_T0'; used with too many (2) args
`THREAD_T0'; used with too many是怎么回事啊[br][br][以下内容由 weiyj83 在 2007年07月25日 07:33pm 时添加] [br]
我是在帮助里面的算例基础上修改一下的。
期待解答，谢谢
qq：35813407

tchllc

THREAD_T0(f,t);改写成THREAD_TO(t)