24小时热门版块排行榜    

查看: 539  |  回复: 0
【悬赏金币】回答本帖问题,作者1294219827将赠送您 30 个金币

1294219827

新虫 (初入文坛)

[求助] VOF模型,管道流动,我想求截面气相所占面积,下面是UDF

VOF模型,我想求截面气相所占面积,下面是UDF,有错误,不知道哪错了,请各位同仁指正!

DEFINE_ON_DEMAND(VOID)
{
    real air_area=0;                                                //存储气相所占面积
    Thread *tf;
    Domain *d;
    face_t f;
    cell_t c0;                                                            //存储与f相邻的单元
    Thread *tc0 = NULL;                                         //存储与f相邻的单元线指针
    d = Get_Domain(3);                                            //得到气相域指针
    tf= Lookup_Thread(d, 12);                                  //获取截面所在的线指针

begin_f_loop(f, tf)
{
        tc0 = Lookup_Thread(d, 26);
        c0 = F_C0(f, tf);                                                                                                                                                                                    
        if (C_VOF(c0, tc0) <= 1 && C_VOF(c0, tc0) >= 0.5)                                                                                   
        {
                F_AREA(A, f, tf);                                                      //获得面向量              
                air_area+ = NV_MAG(A);                                         //求得面的面积                    
        }
}
  end_f_loop(f, t)
  printf("air_area=%g,avg=%g\n", air_area);
}
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 1294219827 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
信息提示
请填处理意见