| 查看: 3914 | 回复: 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算完后发现多孔域流体全都能通过，并不是一段通过，一段不通过。 求大神指导 |
回复此楼» 猜你喜欢
读博
已经有5人回复
-
博士申请都是内定的吗？
已经有6人回复
-
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有5人回复
-
博士读完未来一定会好吗
已经有29人回复
-
投稿精细化工
已经有4人回复
-
高职单位投计算机相关的北核或SCI四区期刊推荐，求支招！
已经有4人回复
-
导师想让我从独立一作变成了共一第一
已经有9人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
-
小论文投稿
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
欧拉模型+多孔介质 求助
已经有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 (小学生)
- 金币: 6628.4
- 散金: 52
- 红花: 29
- 帖子: 542
- 在线: 283.8小时
- 虫号: 1776059
- 注册: 2012-04-24
- 性别: GG
- 专业: 化学工程及工业化学
- 管辖: 仿真模拟
10楼2015-02-06 15:49:35