风工程湍流粘度超限

壬生浪人

我用fluent在做一个标准caarc高层建筑的风荷载模拟下面有我的网格模型采用混合网格，内部为非结构四面体网格，外部为结构网格。

边界层网格用y+为200，在NASA的y+计算器上算得第一层网格为0.008（参考点取建筑顶面182.88m，速度11.7m/s，参考长度为迎风面宽度45.7m）

边界条件采用速度入口，udf格式（用指数风剖面）。
出口为outflow，四周为sym边界，混合网格界面设置interface，其余壁面用不可滑移壁面

问题:计算初始化的时候出现全单元湍流粘度过大

以下为udf内容:

/*profile.c*/
/*UDF for specifying steady-state velocity profile boundary condition*/   
#include"udf.h"   
#include"math.h"
#define zo 182.88
#define uo 11.7
#define zb 5
#define zg 350
#define tur 0.23
#define cmu 0.09

DEFINE_PROFILE(inlet_y_velocity,thread,index)
{
   real x[ND_ND];   
   real z;     
   face_t f;   
   begin_f_loop(f,thread)     
     {
      F_CENTROID(x,f,thread);     
      z=x[2];   
      F_PROFILE(f,thread,index)=uo*pow(z/zo,0.16);
      }
    end_f_loop(f,thread)
}
/*PROFILE FOR KINETIC ENERGY*/
DEFINE_PROFILE(k_profile,thread,index)
{
     real x[ND_ND];
     real z,lr,um;
     face_t f;
     begin_f_loop(f,thread)
    {
       F_CENTROID(x,f,thread);
       z=x[2];
   if(z<zb)
  lr=tur*pow(1.0,1.0);
     else
  lr=0.1*pow(z/zo,-0.20);
     um=uo*pow(z/zo,0.16);
     F_PROFILE(f,thread,index)=1.5*pow(um*lr,2.0);
      }
    end_f_loop(f,thread)
}
/*PROFILE FOR DISSIPATION RATE*/
DEFINE_PROFILE(E_profile,thread,index)
{
  real x[ND_ND];
  real z,lr,um,l,k;
   face_t f;
   begin_f_loop(f,thread)
  {
     F_CENTROID(x,f,thread);
     z=x[2];
  if(z<zb)
  lr=tur*pow(1.0,1.0);
     else
     lr=0.1*pow(z/zo,-0.20);
     um=uo*pow(z/zo,0.16);
     k=1.5*pow(um*lr,2.0);
     l=100*pow(z/30,0.5);
     F_PROFILE(f,thread,index)=pow(cmu,3/4)*pow(k,3/2)/l;
  }
  end_f_loop(f,thread)
}

壬生浪人

顶一下/。。。。。。。。。。。。

壬生浪人

我自己回答一下，y+值设置有误