在fluent中如何设置指数型速度边界条件

wjzsdu

各个大虾，请教一下，在Fluent中如何把入口速度边界条件设置为“指数型速度边界条件”？

不胜感激！

凌白日

自己写个UDF不可以吗？

vincentsuen

/* 定义大气边界层速度，湍流动能，湍流耗散率（alpha=0.16，10m高度风速40m/s，边界层高度350m） */

#include "udf.h"
#define UMEAN 40
#define Alpha 0.16
#define ZB 10.0
#define ZG 350.0

DEFINE_PROFILE(x_velocity,thread,position)
{
  real z,u,x[ND_ND];
  face_t f;
  begin_f_loop(f,thread)
  {  
    F_CENTROID(x,f,thread);
    z=x[1];
    if(z<=ZB)
      u=UMEAN*pow(ZG/10,0.16)*pow(ZB/ZG,Alpha);
    else if(z>ZB&&z<=ZG)
      u=UMEAN*pow(ZG/10,0.16)*pow(z/ZG,Alpha);
    else
      u=UMEAN*pow(ZG/10,0.16)*pow(ZG/ZG,Alpha);
    F_PROFILE(f,thread,position)=u;
   }
  end_f_loop(f,thread)
}


DEFINE_PROFILE(k_profile,thread,position)
{
  real z,ii,u,x[ND_ND];
  face_t f;
  begin_f_loop(f,thread)
  {  
    F_CENTROID(x,f,thread);
    z=x[1];
    if(z<=ZB)
      {
      u=UMEAN*pow(ZG/10,0.16)*pow(ZB/ZG,Alpha);
      ii=0.1*pow(ZB/ZG,-0.05-Alpha);
      }
    else if(z>ZB&&z<=ZG)
      {
      u=UMEAN*pow(ZG/10,0.16)*pow(z/ZG,Alpha);
      ii=0.1*pow(z/ZG,-0.05-Alpha);
      }
    else
      {
      u=UMEAN*pow(ZG/10,0.16)*pow(ZG/ZG,Alpha);
      ii=0.1*pow(ZG/ZG,-0.05-Alpha);
      }
    F_PROFILE(f,thread,position)=0.3*pow(2.5*ii,2);
   }
  end_f_loop(f,thread)
}


DEFINE_PROFILE(dissip_profile,thread,position)
{
  real z,ii,u,x[ND_ND],lx,kay;
  face_t f;
  begin_f_loop(f,thread)
  {  
    F_CENTROID(x,f,thread);
    z=x[1];
    if(z<=ZB)
      {
       u=UMEAN*pow(ZG/10,0.16)*pow(ZB/ZG,Alpha);
      ii=0.1*pow(ZB/ZG,-0.05-Alpha);
      }
    else if(z>ZB&&z<=ZG)
      {
      u=UMEAN*pow(ZG/10,0.16)*pow(z/ZG,Alpha);
      ii=0.1*pow(z/ZG,-0.05-Alpha);
      }
    else
      {
      u=UMEAN*pow(ZG/10,0.16)*pow(ZG/ZG,Alpha);
      ii=0.1*pow(ZG/ZG,-0.05-Alpha);
      }
    kay=0.3*pow(2.5*ii,2);
    lx=100*pow(z/30,0.5);
    F_PROFILE(f,thread,position)=pow(0.09,0.75)*pow(kay,1.5)/(0.07*lx);
   }
  end_f_loop(f,thread)
}

wjzsdu

非常感谢啊~！！！

lcw

3楼朋友能否画个简单的示意图？这样更清楚一些。谢谢

vincentsuen

原帖由 lcw 于 2011-1-8 10:51 发表

                               
登录/注册后可看大图

3楼朋友能否画个简单的示意图？这样更清楚一些。谢谢

大气边界层指数率入口模型示意

eclipse

对楼主的感激之情犹如滔滔江水连绵不绝。小弟这里调用了！