| 查看: 6401 | 回复: 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文章，我:8.O.55.1.O54,科目全,可伽急
已经有8人回复
-
面上没中，邀请各位路过的虫友分析一下分数
已经有12人回复
-
两块石头
已经有8人回复
-
初秋的晨风
已经有6人回复
-
学科评审组评审是指会评吗？
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有14人回复
-
广西大学-广州大学招聘博士后 欢迎广大优秀人才！！！
已经有7人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有17人回复
-
Ei源刊怎么投
已经有4人回复
-
哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Gambit画网格交界面的问题
已经有5人回复
- 汽液界面剪应力的影响
已经有3人回复
- 随着网格数量的增加，VOF模拟的两相流中的液体流动时间变得越来越长？
已经有12人回复
- 关于ANSYS中将带网格的面旋转成带网格的体
已经有3人回复
- 求助如何计算有机实验的产率，新手急需求助，救人一命胜造七级浮屠
已经有5人回复
- 叶面积测量，如何取样？
已经有15人回复
- fluent VOF模型后处理
已经有12人回复
- 各路大神，救人一命，胜造七级浮屠....
已经有8人回复
- 低Re数k-ε模型近壁面区网格加密问题
已经有4人回复
- 【求助】做内标用的正庚烷需要什么纯度，量取时用什么针筒
已经有4人回复
【答案】应助回帖
xiegangmai: 非应助请勿选择“应助回帖” 2012-12-13 23:33:29
| 请教楼主，你的volume和s_g是用什么宏定义的？ 谢谢~~~ 我的邮箱wxw0315@hotmail.com |
10楼2012-12-13 11:30:38
2楼2012-05-23 22:21:27
3楼2012-05-24 16:24:49
4楼2012-05-25 09:52:45