24小时热门版块排行榜

>论坛更新日志 (3413)
>文献求助 (273)
>虫友互识 (273)
>导师招生 (179)
>硕博家园 (108)
>招聘信息布告栏 (91)
>休闲灌水 (78)
>考博 (76)
>论文道贺祈福 (64)
>博后之家 (51)
>考研 (46)
>教师之家 (45)
>公派出国 (43)
>论文投稿 (43)
>基金申请 (40)
>找工作 (28)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

pointh

金虫 (小有名气)

应助: 7 (幼儿园)
金币: 838.6
红花: 7
帖子: 57
在线: 58.1小时
虫号: 1640109
注册: 2012-02-24
专业: 多相流热物理学

[求助] 如何获取汽液界面的单个网格面积？已有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

回复此楼

1楼2012-05-22 20:29:52

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

就是不买日货

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 18.1
帖子: 10
在线: 3.7小时
虫号: 1505812
注册: 2011-11-23
专业: 机械摩擦学与表面技术

★
臭水沟: 金币+1, 谢谢交流~~ 2012-05-26 11:08:27

您好，我觉得只捕捉C_VOF在0，1之间的cell 得到的分界面不全面吧？因为如果这样的话，对于那种分界面但是正好占满了cell，也就是说volume fraction原本就是1的那部分分界面就不包括在这个循环内了...很纠结...

我的关于分界面疑问的帖子在这：
http://muchong.com/bbs/viewthread.php?tid=4547262