|
楼主 |
发表于 2005-10-9 17:18:12
|
显示全部楼层
UDF高手请指教
这是源程序
改变湍流模型的粘性的
#include "udf.h"
#define alpha_0 0.025
#define Re_k 6.
#define C_0 13.
real C_s(cell_t c, Thread *t)
{return fabs(C_U(c,t)*(C_U(c,t)*C_DVDX(c,t)+C_V(c,t)*C_DVDY(c,t))-C_V(c,t)*(C_U(c,t)*C_DUDX(c,t)+C_V(c,t)*C_DUDY(c,t)))
/pow((C_U(c,t)*C_U(c,t)+C_V(c,t)*C_V(c,t)),1.5);}
real R_l(cell_t c, Thread *t)
{return C_s(c,t)*sqrt(C_MU_L(c,t)/C_R(c,t)/C_D(c,t));}
real alpha_1(cell_t c, Thread *t)
{return C_0*R_l(c,t)/(0.1+R_l(c,t));}
real Re_t(cell_t c, Thread *t)
{return C_R(c,t)*C_K(c,t)/C_MU_L(c,t)/C_D(c,t);}
real alpha(cell_t c, Thread *t)
{return (alpha_0+alpha_1(c,t)+Re_t(c,t)/Re_k)/(1+Re_t(c,t)/Re_k);}
DEFINE_ADJUST(turb_adjust, domain)
{
Thread *t;
cell_t c;
/* Set the turbulent viscosity */
thread_loop_c (t, domain)
{
if (FLUID_THREAD_P(t))
{
begin_c_loop(c,t)
{
C_MU_T(c,t)=alpha(c,t)*C_R(c,t)*C_K(c,t)/C_D(c,t);
}
end_c_loop(c,t)
}
}
} |
|