| 查看: 6455 | 回复: 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» 猜你喜欢
师弟论文见刊大半年才想起来申请专利，还能抢救一下吗？
已经有3人回复
-
上海工程技术大学 激光智能制造课题组 2027级博士研究生招生
已经有5人回复
-
课题组招2027级博士 上海工程技术大学 激光智能制造方向
已经有5人回复
-
有机合成以后会不会被AI改变？做科研的虫友怎么看
已经有8人回复
-
求合成方法
已经有7人回复
-
现代”学阀”该如何界定
已经有12人回复
-
我的奶奶
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Gambit画网格交界面的问题
已经有5人回复
- 汽液界面剪应力的影响
已经有3人回复
- 随着网格数量的增加，VOF模拟的两相流中的液体流动时间变得越来越长？
已经有12人回复
- 关于ANSYS中将带网格的面旋转成带网格的体
已经有3人回复
- 求助如何计算有机实验的产率，新手急需求助，救人一命胜造七级浮屠
已经有5人回复
- 叶面积测量，如何取样？
已经有15人回复
- fluent VOF模型后处理
已经有12人回复
- 各路大神，救人一命，胜造七级浮屠....
已经有8人回复
- 低Re数k-ε模型近壁面区网格加密问题
已经有4人回复
- 【求助】做内标用的正庚烷需要什么纯度，量取时用什么针筒
已经有4人回复
3楼2012-05-24 16:24:49
2楼2012-05-23 22:21:27
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