FLUENT厂房屋盖数值模拟的设置问题想问问大家
我在做厂房屋盖风压模拟的时候遇到了一些设置问题，大家帮忙一下：
X正向为流体前进方向，建筑物逆时针旋转作为风向角
SST K-W湍流模型，入口流速和湍动能（日本规范建议）采用UDF编写
松弛因子全部默认
离散格式如图：
1，试算的时候出现turbulence viscous rate超过限制的问题，有文献认为是网格质量，Y+值和湍流量的问题，但是前两个问题在SST K-W湍流模型里应该很容易就能达到要求，湍流量是使用UDF编写，文后我附上UDF文档，大家请帮忙检查一下，我实在没办法就在下面这个图片里在Maximum Turb.Viscosity Ratio后面添加5个0，试算完检查Y+值是符合要求的，那么我这样是否可行？
2，我对该数值模拟需要得到的数据是屋盖的风压值（升力）Z轴为竖直向上，以及屋盖弯矩（右手法则，拇指指向和Y轴平行）。那么我下图升力系数的Force Vector和弯矩监视器Moment Center与Moment Axis该怎么设置？
最后，附上UDF文件
#include "udf.h"
DEFINE_PROFILE(velocity_inlet, thread, position) 定义速度入口Vinlet
{
real x[ND_ND];
real z;
face_t f;
begin_f_loop(f, thread)
{
F_CENTROID(x,f,thread);
z=x[2]/10;
F_PROFILE(f, thread, position) =26.83*pow(z,0.15);
end_f_loop(f, thread)
}
}
DEFINE_PROFILE(k_inlet, thread, position) 根据日本规范建议定义湍动能K
{
real x[ND_ND];
real z;
real i;
face_t f;
begin_f_loop(f, thread)
{
F_CENTROID(x,f,thread);
z=x[2];
if(z<5)
{i=0.23;}
if(z>=5,z<=350)
{i=0.1*pow(z/350,-0.2);}
if(z>350)
{i=0.1;}
F_PROFILE(f, thread, position) =1.5*(26.83*pow(z/10,0.15)*i)*(26.83*pow(z/10,0.15)*i);
end_f_loop(f, thread)
}
}
DEFINE_PROFILE(w_inlet, thread, position) 定义比耗散率W
{
real x[ND_ND];
real z;
real i;
real l;
real k;
real e;
face_t f;
begin_f_loop(f, thread)
{
F_CENTROID(x,f,thread);
z=x[2];
if(z<5)
{i=0.23;}
if(z>=5,z<=350)
{i=0.1*pow(z/350,-0.20);}
if(z>350)
{i=0.1;}
l=100*pow(z/30,0.5);
k=1.5*(26.83*pow(z/10,0.15)*i)*(26.83*pow(z/10,0.15)*i);
e=pow(0.09,0.75)*pow(k,1.5)/l;
F_PROFILE(f, thread, position) =e/(0.09*k);
end_f_loop(f, thread)
}
}
