24小时热门版块排行榜    

查看: 373  |  回复: 1

sghudng

铜虫 (小有名气)

[交流] 【求助】麻烦哪位大侠帮我看看这个UDF小程序啊~~ 已有1人参与

我这个是体积分数的,需要调用一个面相邻网格上的体积分数,然后求平均赋值给入口边界条件里的第二相体积分数,程序如下:能导进去计算,可是结果好像是没有加进去……

#include "udf.h"
DEFINE_PROFILE(vof,t,i)
{
        Domain *domain=Get_Domain(1);
    Thread *tc_s,*mix_thread,*tc;
        real vof,v,v_tot,vof_ave;
        cell_t c;
    face_t f;
        v_tot=0.0;
        vof_ave=0.0;


        mix_thread=Lookup_Thread(domain,5);

        begin_f_loop(f,mix_thread)
    {
        c=F_C0(f,mix_thread);
        tc=THREAD_T0(mix_thread);
    tc_s=THREAD_SUB_THREAD(tc,1);
        vof=C_VOF(c,tc_s);
        v=C_VOLUME(c,tc);
    v_tot+=v;
    vof_ave+=vof*v;
    }
    end_f_loop(f,mix_thread)
        vof_ave/=v_tot;

    begin_f_loop(f,t)
    {
    F_PROFILE(f,t,i)=vof_ave;
    }
    end_f_loop(f,t)
}
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sghudng

铜虫 (小有名气)

zzuwangshilei:希望解释的更详细点,多谢解释 2010-07-02 22:02:11
算了之后第二相跑了一大半……
2楼2010-07-02 12:21:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 sghudng 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见