| 查看: 6437 | 回复: 16 | |||
[求助]
如何获取汽液界面的单个网格面积？ 已有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 |
回复此楼» 收录本帖的淘帖专辑推荐
 模拟相关 |
» 本帖已获得的红花（最新10朵）
pointh» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Gambit画网格交界面的问题
已经有5人回复
- 汽液界面剪应力的影响
已经有3人回复
- 随着网格数量的增加，VOF模拟的两相流中的液体流动时间变得越来越长？
已经有12人回复
- 关于ANSYS中将带网格的面旋转成带网格的体
已经有3人回复
- 求助如何计算有机实验的产率，新手急需求助，救人一命胜造七级浮屠
已经有5人回复
- 叶面积测量，如何取样？
已经有15人回复
- fluent VOF模型后处理
已经有12人回复
- 各路大神，救人一命，胜造七级浮屠....
已经有8人回复
- 低Re数k-ε模型近壁面区网格加密问题
已经有4人回复
- 【求助】做内标用的正庚烷需要什么纯度，量取时用什么针筒
已经有4人回复
6楼2012-05-26 16:17:29
2楼2012-05-23 22:21:27
3楼2012-05-24 16:24:49
4楼2012-05-25 09:52:45
★
臭水沟: 金币+1, 谢谢交流~~ 2012-05-26 11:08:27
臭水沟: 金币+1, 谢谢交流~~ 2012-05-26 11:08:27
|
您好，我觉得只捕捉C_VOF在0，1之间的cell 得到的分界面不全面吧？因为如果这样的话，对于那种分界面但是正好占满了cell，也就是说volume fraction原本就是1的那部分分界面就不包括在这个循环内了...很纠结... 我的关于分界面疑问的帖子在这： http://muchong.com/bbs/viewthread.php?tid=4547262 |
5楼2012-05-26 10:49:31
★
xiegangmai: 金币-1, 专业版交流贴，请勿灌水！ 2012-05-26 23:33:11
xiegangmai: 金币-1, 专业版交流贴，请勿灌水！ 2012-05-26 23:33:11
 |
7楼2012-05-26 18:26:05
★ ★
xiegangmai: 金币-2, 专业版交流贴，请勿灌水！ 2012-05-26 23:33:26
xiegangmai: 金币-2, 专业版交流贴，请勿灌水！ 2012-05-26 23:33:26
|
|
8楼2012-05-26 18:26:28
9楼2012-05-29 14:49:54
【答案】应助回帖
xiegangmai: 非应助请勿选择“应助回帖” 2012-12-13 23:33:29
| 请教楼主，你的volume和s_g是用什么宏定义的？ 谢谢~~~ 我的邮箱wxw0315@hotmail.com |
10楼2012-12-13 11:30:38
10