我这个UDF编写的对不对啊？

zf392222982

湍流模型采用 SST K-omiga
比如我现在以B类地貌，0.45基本风压为基础（参考处流速26.83）
#include "udf.h"               
DEFINE_PROFILE(velocity_inlet, thread, position)               #定义速度入口
{
  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)                     #定义湍动能
{
  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)                 #定义比耗散率
{
  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)
    }
}

以上我这样编写的UDF文件对吗？

小七工作室

好多文献都有相关这个例子  风环境的