| 查看: 5640 | 回复: 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» 猜你喜欢
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？
已经有6人回复
-
带资进组求博导收留
已经有10人回复
-
最近几年招的学生写论文不引自己组发的文章
已经有11人回复
-
需要合成515-64-0，50g，能接单的留言
已经有3人回复
-
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）
已经有3人回复
-
临港实验室与上科大联培博士招生1名
已经有8人回复
-
想换工作。大多数高校都是 评职称时 认可5年内在原单位取得的成果吗？
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Gambit画网格交界面的问题
已经有5人回复
- 汽液界面剪应力的影响
已经有3人回复
- 随着网格数量的增加，VOF模拟的两相流中的液体流动时间变得越来越长？
已经有12人回复
- 关于ANSYS中将带网格的面旋转成带网格的体
已经有3人回复
- 求助如何计算有机实验的产率，新手急需求助，救人一命胜造七级浮屠
已经有5人回复
- 叶面积测量，如何取样？
已经有15人回复
- fluent VOF模型后处理
已经有12人回复
- 各路大神，救人一命，胜造七级浮屠....
已经有8人回复
- 低Re数k-ε模型近壁面区网格加密问题
已经有4人回复
- 【求助】做内标用的正庚烷需要什么纯度，量取时用什么针筒
已经有4人回复
4楼2012-05-25 09:52:45
2楼2012-05-23 22:21:27
3楼2012-05-24 16:24:49
★
臭水沟: 金币+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