| 查看: 4455 | 回复: 9 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
fluent多孔介质孔隙率是位置的函数问题求助 已有3人参与
|
|||
|
想要给一块板一段流体能通过，一段流体不能通过，我想通过定义孔隙率为0和1来实现 定义了一个流体域，一个多孔域  对多孔域使用了UDF #include "udf.h" DEFINE_PROFILE(Bed_porosity,t,i) { real xc[ND_ND]; cell_t c; begin_c_loop(c,t) { if (xc[0]>-606 && xc[0]<-600,xc[1]>320 && xc[1]<=470,xc[2]>0 &&xc[2]<100) { F_PROFILE(c,t,i)=1; } else if (xc[0]>-606 && xc[0]<-600,xc[1]>470 && xc[1]<620,xc[2]>0 &&xc[2]<100) { F_PROFILE(c,t,i)=0; } } end_c_loop(c,t) } 在fluent里编译也通过了，但是fluent算完后发现多孔域流体全都能通过，并不是一段通过，一段不通过。 求大神指导 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
-
什么时候开奖？
已经有10人回复
-
跳槽后在研项目怎么办？
已经有9人回复
-
2026国自然函评费到账
已经有16人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
让我中一个面上吧！
已经有14人回复
-
今天放榜吗？
已经有16人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有10人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 欧拉模型+多孔介质 求助
已经有4人回复
- 求助关于固定床多孔介质孔隙率的问题
已经有7人回复
- fluent中UDF编译的问题！
已经有11人回复
- 多孔介质孔隙率udf问题
已经有7人回复
- fluent中多空介质的计算
已经有19人回复
- fluent多孔介质设置问题！！求助
已经有10人回复
- 【求助】请问在fluent中怎么设置多孔介质的内热源啊
已经有13人回复
- 多孔介质孔隙率的声明
已经有4人回复
- 急求各位大侠啊~fluent里多孔介质模型参数设置的问题~
已经有5人回复
- 多孔介质在gambit中建模后边界条件设置问题
已经有15人回复
- 急求！！请问fluent里如何利用多孔介质模型模拟两相流吸附过滤？
已经有8人回复
- Fluent多孔介质模拟问题
已经有15人回复
- 【求助】求助啦~双相介质和孔隙介质是一回事吗？
已经有3人回复
- 【求助】大侠支招，另急需FLUENT多孔介质的算例
已经有11人回复
|
谢谢您的解答。 我模拟的过程中发现多孔域内部和外部的速度一样，虽然我确实发现了透过多孔介质有压力降。本来我定义的阻力系数大的地方速度应该小，阻力系数小的得放速度应该大。 我在文献中发现 由于数值模拟中没有表示多孔介质区域的实际几何模型存在，因此 FLUENT 中认为在多孔介质区域内的速度与区域外气流的速度相同。这样无法 模拟多孔介质区域对气流速度增加的影响。————《防风栅气动性能的风洞试验及数值模拟研究》_李兆杨 于是 在标准动量方程中增加动量源项，我们可以通过在 FLUENT 中采用 UDF 编译的方法直接计算多孔介质区域的动量方程源项。由于疏透结构自身对气流 的体积阻塞作用，会对气流产生局部加速作用，如前所述，多孔介质模型未能 考虑这种局部加速作用，我们可以通过考虑单元平均体积疏透率，通过 UDF 编 译程序来求解真实速度或物理速度————《防风栅气动性能的风洞试验及数值模拟研究》_李兆杨 附上他的UDF #include "udf.h" #define C2 20112500e+001 /*porous loss coefficients*/ #define G 8.200000e-001 /*porous stands for the value per unit*/ DEFINE_ADJUST(velocity_adjust,domain) { int ID=16; /* zone ID for porous zone */ Thread *thread=Lookup_Thread[domain,16]; cell_t,c; real x[ND_ND];/*this will hold the position vector*/ begin_c_loop(c,thread) { C_CENTROID(x,c,thread); C_U(c,thread)=G*C_U(c,thread); C_V(c,thread)=G*C_V(c,thread); C_W(c,thread)=G*C_W(c,thread); } end_c_loop(c,thread) } /*source define for x-momentum*/ DEFINE_SOURCE(xmom_source,c,t,dS,eqn) { real x[ND_ND];/*this will hold the position vector*/ real con,source; C_CENTROID(x,c,t); con=c2*0.5*C_R(c,t);/*C_R(c,t) terms density of air*/ source=0.001-con*fabs(C_U(c,t))*C_U(c,t); dS[eqn]=-2*con*fabs(C_U(c,t)); return source |
9楼2014-11-25 10:50:04
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
本帖内容被屏蔽 |
2楼2014-09-28 12:51:21
3楼2014-09-28 16:04:03
4楼2014-10-09 09:18:14
