1 前言
2 建模与网格
3 求解设置
#include "udf.h"double x[ND_ND];double z;double uH=4.3;double h0=0.16;DEFINE_PROFILE(velocity_normal, thread, position){face_t f;begin_f_loop(f,thread){F_CENTROID(x,f,thread);z=x[2];F_PROFILE(f,thread,position)=pow(z/h0,0.25)*uH;}end_f_loop(f,thread)}DEFINE_PROFILE(kinetic, thread, position){face_t f;begin_f_loop(f,thread){F_CENTROID(x,f,thread);z=x[2];F_PROFILE(f,thread,position)=pow(uH,2)*0.033*exp(-0.32*z/h0);}end_f_loop(f,thread)}
#include "udf.h"
double x[ND_ND];
double z;
double uH=4.3;
double h0=0.16;
DEFINE_PROFILE(velocity_normal, thread, position)
{
face_t f;
begin_f_loop(f,thread)
F_CENTROID(x,f,thread);
z=x[2];
F_PROFILE(f,thread,position)=pow(z/h0,0.25)*uH;
}
end_f_loop(f,thread)
DEFINE_PROFILE(kinetic, thread, position)
F_PROFILE(f,thread,position)=pow(uH,2)*0.033*exp(-0.32*z/h0);
4 计算结果
参考文献