| 查看: 4069 | 回复: 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算完后发现多孔域流体全都能通过，并不是一段通过，一段不通过。 求大神指导 |
回复此楼» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有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人回复
★
小木虫: 金币+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
5楼2014-10-09 17:00:20
6楼2014-10-10 14:29:07
送红花一朵 |
您好，我一直调试，看了一些别的例子，似乎要加上动量源项。。下面是我看到的一个二维多孔介质用的udf。。。 #include "udf.h" DEFINE_SOURCE(xmom_source,cell,thread,dS,eqn) { const real c2=100; real x[ND_ND]; 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; } |
7楼2014-11-25 09:26:46
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
本帖内容被屏蔽 |
8楼2014-11-25 10:16:01
|
谢谢您的解答。 我模拟的过程中发现多孔域内部和外部的速度一样，虽然我确实发现了透过多孔介质有压力降。本来我定义的阻力系数大的地方速度应该小，阻力系数小的得放速度应该大。 我在文献中发现 由于数值模拟中没有表示多孔介质区域的实际几何模型存在，因此 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
xingfuww
专家顾问 (正式写手)
-
专家经验: +78 - 仿真EPI: 2
- 应助: 46 (小学生)
- 金币: 6630.4
- 散金: 52
- 红花: 29
- 帖子: 542
- 在线: 283.8小时
- 虫号: 1776059
- 注册: 2012-04-24
- 性别: GG
- 专业: 化学工程及工业化学
- 管辖: 仿真模拟
10楼2015-02-06 15:49:35
