|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?注册
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)
}
} |
|