急急！！关于udf的问题，请大侠帮助。。。

hongfu2233

我做得是风吹建筑模拟，剖面风已有udf控制，地面风速应该为0，但输入flunt观察之后不为0，请哪位大侠看看下面的udf（关于剖面风的）哪里出错了，在下不胜感激！！！

#include "udf.h"

#define a 0.22
#define height 10
#define v_height 5.194
#define Cu 0.09
#define Zg 300
#define Zb 5
#define VKC 0.41
#define scale 1

DEFINE_PROFILE(inlet_v_profile,thread,index)
{
       
        real x[ND_ND];  /*this will hold the position vector*/
        real z;
        face_t f;
               
        begin_f_loop(f,thread)
        {      
                F_CENTROID(x,f,thread);
                z=scale*x[2];/*x=[ND_ND]is one-dimensioned array;x=x[0];y=x[1];z=x[2]*/
                if(z<=Zg)
                  F_PROFILE(f,thread,index)=v_height*pow((z/height),a);
                else
                  F_PROFILE(f,thread,index)=v_height*pow((Zg/height),a);
        }
        end_f_loop(f,thread)
}

2005011lcj

第一次看见有人这么定义的
#define a 0.22
#define height 10
#define v_height 5.194
#define Cu 0.09
#define Zg 300
#define Zb 5
#define VKC 0.41
#define scale 1

liuoseulb

udf没有问题，你说的“地面风速不为0”，在fuent里面显示的是哪里的值？是指地面所在的边界面上，还是地面附近？

hongfu2233

地面上不为0，我取了一个z=0的面，就是地面，用area average-weighted求解了一下地面加权平均风速竟然为4m/s，另外看了一下背面风压线，竟然为正值就是压力为正。请大侠指示。。。

hongfu2233

不对吗？？！！

2005011lcj

对，只是我用不习惯

liuoseulb

你地面边界条件设的是wall吗?

hongfu2233

对 是wall  类型no slip

liuoseulb

那我觉得可能是算法不同造成的，就是节点值和单元值。因为本身z=0的位置归第一层网格管，如果你显示的是单元值，它就不是0，应该是第一层网格范围内都是这个值；如果你用的是节点值，通过插值得出的就应该是0
你可以用display--->contour，去显示z=0的平面风速，然后比较勾上node value和不勾的结果

hongfu2233

display--->contour中没有显示速度这个选项，display--->vector中有显示速度选项但是没有node value选项，大侠，主要是我做的模拟 背风面风压是负值有悖常理，能否指导一二？多谢了。。。

liuoseulb

contour中怎么可能没有velocity选项？？？你仔细看看
迎风面正压，背风面负压，这就是常理

hongfu2233

哦 不好意思，说错了 背风面也是正压，是哪里错了呢？

uesoft

原帖由 2005011lcj 于 2011-9-1 11:04 发表

                               
登录/注册后可看大图

对，只是我用不习惯

估计他是非计算机专业的，能把C和Fluent搞到这么熟练的程度还是值得鼓励的。
不过他最好还是按照编程规范来写代码，便于代码评审、沟通交流。

[ 本帖最后由 uesoft 于 2011-9-2 18:30 编辑 ]