|
定义一个y+,带入到udf子程序的公式中进行计算。
假如是圆管紊流,y+定义为,
Ypuls=y*u/mu,
其中,y为到壁面的距离,宏定义为C_WALL_DIST(c,t);
u为摩擦速度,定义u=power((Tw/den),0.5), 摩擦速度好像也没有宏定义
Tw为壁面切应力=0.25D*(压力梯度),D为管直径,压力梯度有宏定义为C_P_G(c,t)
den为密度,有宏定义C_R(c,t);
mu为运动粘度,有宏定义。
udf程序如下:
/* Yplus definitions */
real Y_plus(cell_t c, Thread *t)
{
return C_R(c,t)*pow(0.25*D0*C_P_G(c,t)/C_R(c,t),0.5)*C_WALL_DIST(c,t)/C_MU_L(c,t);
}
|
|