24小时热门版块排行榜    

查看: 5780  |  回复: 16
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

pointh

金虫 (小有名气)

[求助] 如何获取汽液界面的单个网格面积? 已有3人参与

Sample Text
       本人现在正在用VOF模型模拟单个气泡的生长,现已知通过汽液界面的总的质量流量(单位:kg/(m2.s)),由于Fluent是基于有限容积法,想转化为单个网格单元的源项(kg/(m3.s)),但在写UDF时,单个网格的面积不知该怎么处理。我现在有两种想法,但好像结果都不对。部分代码如下:
1. begin_c_loop (c,mixture)
          {     if(0<=C_VOF(c,t)&&C_VOF(c,t)<1.0)           //考虑气液界面处网格,有人说参考VOF分相模型,应该选0.5~1之间,我一直没找到,为什么?
                {        c_face_loop(c,t,n)                                //循环网格里的面
                        {        f=C_FACE(c,t,n);         
                                tf=C_FACE_THREAD(c,t,n);        //获取面指针
                                F_AREA(A,f,tf);
                                area=NV_MAG(A);                           //获取网格某个面的面积
                        }
                       mass_gas=G*area*C_VOF(c,t)/volume;    //将总的质量流量转为单个网格的质量流率
               }       
          }
      end_c_loop (c,mixture)
2.  begin_c_loop (c,mixture)
          {           if(0.<=C_VOF(c,t)&&C_VOF(c,t)<1.0)
                {         num=num+1;                               //获取气液界面处的网格个数
                }
          }
      end_c_loop (c,mixture)
      begin_c_loop (c,mixture)
          {        if(0.<=C_VOF(c,t)&&C_VOF(c,sp)<1.0)
                {         area=s_g/num;                           //气泡的表面积除以网格数得每个网格的面积
                        mass_gas=G*area*C_VOF(c,t)/volume;
                }       
           }
      end_c_loop (c,mixture)
        对于第一种方法,F_AREA获取的是网格指向法线外的那个面吗?循环了一下不知有没有关系?(这是因为我只看到了用c_face_loop循环来获取网格内的面指针,不知有没其他方法?)而且气液界面一定是在网格外吗?(对于界面重构是不是可能出现在网格内)
        对于第二种方法,由于需要局部加密,所以网格单元大小不一样,且每个网格内气液界面的面积也不一样,这个我认为问题很大。
       另外,我想在监视窗口显示通过整个气泡壁面的质量流率,但Report Type那一项好像始终没有(单位应该是kg/(m2.s),但那个Mass Flow Rate监测出来的是kg/s)。这个又该怎么解决呢?

气泡模型



公式



网格内的气液界面



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

limeng390

银虫 (初入文坛)

最近看了这样的帖子,有很多启发,我想vof在边界上做一个空化
有容乃大,无欲则刚
16楼2016-05-22 11:22:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 17 个回答

pointh

金虫 (小有名气)

都没有高手看一下吗?
2楼2012-05-23 22:21:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pointh

金虫 (小有名气)

我看了一下,有很多人都对这个问题有疑问,俗话说“救人一命,胜造七级浮屠”,这可能会救很多条人命啊,为什么没人回答呢?
3楼2012-05-24 16:24:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pointh

金虫 (小有名气)

送鲜花一朵
别让帖子沉到地板上啊!这个应该可以为很多做VOF两相流模拟的打开思路啊,怎么都没高手帮忙解答一下啊?辛苦一人,幸福千万家啊!
4楼2012-05-25 09:52:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +5 pnpwoqbg8f 2026-02-16 5/250 2026-02-18 18:25 by 7pdbc5bdo2
[博后之家] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 18:12 by 7pdbc5bdo2
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 18:10 by 7pdbc5bdo2
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 17:57 by 7pdbc5bdo2
[基金申请] 情人节自我反思:在爱情中有过遗憾吗? +5 瞬息宇宙 2026-02-15 6/300 2026-02-18 12:51 by 月下雪林
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 11:09 by lqtl9djx19
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:54 by lqtl9djx19
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:39 by lqtl9djx19
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:53 by lqtl9djx19
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 07:55 by lotyj5cz79
[基金申请] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:40 by lotyj5cz79
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:38 by lotyj5cz79
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:23 by lotyj5cz79
[论文投稿] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 00:40 by tk2gfblvuz
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 00:23 by tk2gfblvuz
[基金申请] 基金正文30页指的是报告正文还是整个申请书 +3 successhe 2026-02-16 4/200 2026-02-17 20:56 by successhe
[基金申请] 今年春晚有几个节目很不错,点赞! +5 瞬息宇宙 2026-02-16 6/300 2026-02-17 12:49 by jymy19840415
[微米和纳米] 球磨粉体时遇到了大的问题,请指教! 10+3 6sbiam 2026-02-12 15/750 2026-02-16 15:03 by tgzxzqj
[基金申请] 过年走亲戚时感受到了所开私家车的鄙视链 +3 瞬息宇宙 2026-02-15 5/250 2026-02-16 14:23 by aspect3000
[硕博家园] 江汉大学解明教授课题组招博士研究生/博士后 +3 cleverlyy 2026-02-12 3/150 2026-02-12 21:02 by qsdf1
信息提示
请填处理意见