找回密码
 注册
查看: 1213|回复: 0

大家帮忙看看UDF,编译出错阿。

[复制链接]
发表于 2008-3-18 19:54:43 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?注册

x
#include "udf.h"
#define D_h 9.23E-5
#define H 2E-6
real Re(cell_t c, Thread *t)
{
return C_R(c,t)*D_h*C_W(c,t)/C_MU_L(c,t);
}
real Re_h(cell_t c, Thread *t)
{ real c=0.1306*pow(D_h/(2*H),0.3693)*exp(Re(c,t)*(6E-5*D_h/(2*H)-0.0029));
  real u_h=H*C_W_G(c,t);
  return C_R(c,t)*u_h*H/C_MU_L(c,t);
}

DEFINE_ADJUST(turb_adjust, domain)
{
Thread *t;
cell_t c;
        real x[ND_ND];
        C_CENTROID(x, c, t);
        real mu = C_MU_L(c,t);
    real mu_r = c*mu*x[1]*Re_h(c,t)*SQR(1-exp(-Re_h(c,t)*x[1]/(Re(c,t)*H));
thread_loop_c (t, domain)
if (FLUID_THREAD_P(t))
{
begin_c_loop(c,t)
{
    C_MU_T(c,t) = mu+mu_r;
}
end_c_loop(c,t)
}
}
您需要登录后才可以回帖 登录 | 注册

本版积分规则

快速回复 返回顶部 返回列表