[求助]高手们指点一下吧！

草原鹿

用udf程序实现入口处的速度，湍动能。为什么后处理时发现入口湍流强度很大（200％～300％）呢？
下面是我的udf程序，高手帮我看看：
#include "udf.h"
#define UFREE 15
#define A 0.175
#define B 0.16
DEFINE_PROFILE(inlet_velocity,thread,zys)
{real z,x[3];
face_t f;
begin_f_loop(f,thread);
{F_CENTROID(x,f,thread);
z=x[2];
F_PROFILE(f,thread,zys)=1.0*UFREE*pow(z/10,B);
}
end_f_loop(f,thread)
}
DEFINE_PROFILE(inlet_K_profile,thread,zys)
{real z,x[3];
face_t f;
begin_f_loop(f,thread);
{F_CENTROID(x,f,thread);
z=x[2];
F_PROFILE(f,thread,zys)=1.5*UFREE*UFREE*A*A*pow(z/10,-0.1);
}
end_f_loop(f,thread)
}
上面的湍动能由k=(3/2)*(U*I)^2计算的，U=UFREE*(z/10)^B,I=A*(z/10)^(-B-0.05),z为大气边界层的高度！B为B类场地的幂指数！
谢谢！