求大神帮忙修改适合SST K-W湍流模型入口边界的UDF
本帖最后由 zf392222982 于 2013-6-28 15:07 编辑这是我正在使用的适用于K-E湍流模型的UDF，哪位大神能帮我改写成适合SST K-W湍流模型的UDF？
以下的数据含义如下：
26.83——根据《建筑荷载规范》取50年一遇的基本风压换算成的10米高处的风速；
350——《建筑荷载规范》规定的梯度风高度；
0.16——《建筑荷载规范》规定的B累场地地面粗糙度
#include "udf.h"
DEFINE_PROFILE(velocity_inlet, thread, position)
{
real x;
real z;
face_t f;
begin_f_loop(f, thread)
{
F_CENTROID(x,f,thread);
z=x/10;
F_PROFILE(f, thread, position) =26.83*pow(z,0.16);
end_f_loop(f, thread)
}
}
DEFINE_PROFILE(k_inlet, thread, position)
{
real x;
real z;
real i;
face_t f;
begin_f_loop(f, thread)
{
F_CENTROID(x,f,thread);
z=x;
if(z<5)
{i=0.23;}
if(z>=5,z<=350)
{i=0.1*pow(z/350,-0.21);}
if(z>350)
{i=0.1;}
F_PROFILE(f, thread, position) =1.5*(26.83*pow(z/10,0.16)*i)*(26.83*pow(z/10,0.16)*i);
end_f_loop(f, thread)
}
}
DEFINE_PROFILE(e_inlet, thread, position)
{
real x;
real z;
real i;
real l;
real k;
face_t f;
begin_f_loop(f, thread)
{
F_CENTROID(x,f,thread);
z=x;
if(z<5)
{i=0.23;}
if(z>=5,z<=350)
{i=0.1*pow(z/350,-0.21);}
if(z>350)
{i=0.1;}
l=100*pow(z/30,0.5);
k=1.5*(26.83*pow(z/10,0.16)*i)*(26.83*pow(z/10,0.16)*i);
F_PROFILE(f, thread, position) =pow(0.09,0.75)*pow(k,1.5)/l;
end_f_loop(f, thread)
}
}
页:
[1]