| 查看: 1579 | 回复: 4 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
求助FLUENT VOF模拟气液分界问题 已有1人参与
|
||
| 在一个封闭旋转容器中，利用VOF模型，材料为水和空气，初始条件patch为气体的体积分数为0.5 ，进行稳态计算。按照定性分析，因为离心力的作用，当容器旋转时应该出现液体在边缘，气体在中间的分布。但是实际的仿真结果是整个容积中体积分数还是在0.4~0.6左右，并没有出现明显的气液分界面。不知道这是什么原因造成的呢？ |
回复此楼» 猜你喜欢
职称评审没过，求安慰
已经有56人回复
-
最近几年招的学生写论文不引自己组发的文章
已经有5人回复
-
26申博自荐
已经有3人回复
-
A期刊撤稿
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助fluent里面的动压计算~两种计算方式结果不一样
已经有1人回复
- Fluent 读入case文件出现receive fatal signal(access-violation)问题 菜鸟求助哈！
已经有1人回复
- fluent udf三维力源项问题求助
已经有14人回复
- 求助大神FLUENT中UDF关于宏函数及返回值的问题
已经有3人回复
- 求助徐子圭博士的fluent udf培训视频！！！！
已经有8人回复
- 求助关于FLUENT压力系数&moment coefficient&adapt的问题
已经有0人回复
- FLUENT 两相流udf求助
已经有0人回复
- 求助fluent动量源项添加及UDFs代码生成 DEFINE_SOURCE
已经有2人回复
- ANSYS FLUENT 14 UDF 编译错误 求助
已经有15人回复
- 【求助】FLUENT UDF 气液两相流 汽蚀 汽化潜热(问题详见内容)
已经有21人回复
- 求助：fluent的沸腾UDF是否可以应用的三维中
已经有0人回复
- 液膜表面蒸发fluent求助！
已经有2人回复
- 【求助】fluent不能启动，重装后还是老样子
已经有5人回复
- 【求助】fluent模拟气固流化床采用欧拉模型并行计算出现问题
已经有12人回复
- 【求助】linux版fluent使用udf并行问题请教
已经有4人回复
- 【求助】fluent计算错误 。求高手指导哪里出错。
已经有8人回复
4楼2015-12-21 23:28:28
卑良2楼2015-12-21 12:16:42
3楼2015-12-21 14:40:33
|
vof格式中，稳态计算只能使用隐格式，而geo-recon是显格式，所以不能选 Fluent中有Degassing边界条件，可以设置某相流体单独流出，但只在eularian两相气液流模型里可以直接使用 对VOF两相流，要使用UDF，在help里的例子贴在下面： 2.3.24.4. Example 2 - Degassing Boundary Condition The following UDFs are used to define the bottom surface as a standard velocity inlet for the gas (primary) phase. The inlet VOF of the droplet phase is 0 and a negative source term for secondary phase mass conservation is set for the layer of cells next to the outlet. The source removes all secondary phase mass in the cell during one time step. The recoil force due to the mass source is also calculated. /*******************************************************************/ This UDF is an implementation of the degassing boundary condition *******************************************************************/ #include "udf.h" #include "sg.h" #include "sg_mphase.h" #include "flow.h" #include "mem.h" #include "metric.h" DEFINE_SOURCE(degassing_source, cell, thread, dS, eqn) { real source; Thread *tm = THREAD_SUPER_THREAD(thread); source = -C_R(cell,thread)*C_VOF(cell,thread)/CURRENT_TIMESTEP ; C_UDMI(cell,tm,0) = source; dS[eqn] = -C_R(cell,thread)/CURRENT_TIMESTEP; return source; } DEFINE_SOURCE(x_prim_recoil, cell, tp, dS, eqn) { real source; Thread *tm = THREAD_SUPER_THREAD(tp); Thread *ts; ts = THREAD_SUB_THREAD(tm,1); source = -C_R(cell,ts)*C_VOF(cell,ts)/CURRENT_TIMESTEP*C_U(cell,tp); dS[eqn] =-C_R(cell,ts)*C_VOF(cell,ts)/CURRENT_TIMESTEP; return source; } DEFINE_SOURCE(x_sec_recoil, cell, ts, dS, eqn) { real source; Thread *tm = THREAD_SUPER_THREAD(ts); source = -C_R(cell,ts)*C_VOF(cell,ts)/CURRENT_TIMESTEP*C_U(cell,ts); dS[eqn] = -C_R(cell,ts)*C_VOF(cell,ts)/CURRENT_TIMESTEP; return source; } DEFINE_SOURCE(y_prim_recoil, cell, tp, dS, eqn) { real source; Thread *tm = THREAD_SUPER_THREAD(tp); Thread *ts; ts = THREAD_SUB_THREAD(tm,1); source = -C_R(cell,ts)*C_VOF(cell,ts)/CURRENT_TIMESTEP*C_V(cell,tp); dS[eqn] =-C_R(cell,ts)*C_VOF(cell,ts)/CURRENT_TIMESTEP; return source; } DEFINE_SOURCE(y_sec_recoil, cell, ts, dS, eqn) { real source; Thread *tm = THREAD_SUPER_THREAD(ts); source = -C_R(cell,ts)*C_VOF(cell,ts)/CURRENT_TIMESTEP*C_V(cell,ts); dS[eqn] = -C_R(cell,ts)*C_VOF(cell,ts)/CURRENT_TIMESTEP; return source; } DEFINE_SOURCE(z_prim_recoil, cell, tp, dS, eqn) { real source; Thread *tm = THREAD_SUPER_THREAD(tp); Thread *ts; ts = THREAD_SUB_THREAD(tm,1); source = -C_R(cell,ts)*C_VOF(cell,ts)/CURRENT_TIMESTEP*C_W(cell,tp); dS[eqn] =-C_R(cell,ts)*C_VOF(cell,ts)/CURRENT_TIMESTEP; return source; } DEFINE_SOURCE(z_sec_recoil, cell, ts, dS, eqn) { real source; Thread *tm = THREAD_SUPER_THREAD(ts); source = -C_R(cell,ts)*C_VOF(cell,ts)/CURRENT_TIMESTEP*C_W(cell,ts); dS[eqn] = -C_R(cell,ts)*C_VOF(cell,ts)/CURRENT_TIMESTEP; return source; } |
5楼2015-12-22 09:40:00