| 查看: 1359 | 回复: 2 | ||
[求助]
气液两相流非均相反应的UDF
|
|
参考之前的一篇帖子，自己写的udf，计算结果始终是没有发生反应。 反应式是A(l)=B(l)+C(g),VOF模型，A和B是混合相（第一相，液相），B第二项气相。希望高手指点。 附程序： #include "udf.h" #include "mem.h" #include "sg_mphase.h" #define FLUID_ZONE_ID 7 /* zone ID of reaction */ #define R 8.31434 /* gas constant J/K-mol */ #define Pre 3.97e+12 /* pre-exponential factors units 1/s */ #define E1 228000 /* Activation energies (J/mol) */ DEFINE_HET_RXN_RATE(my_rate, c, t, r, mw, yi, rr, rr_t) { int zone_ID = 0; Thread **pt = THREAD_SUB_THREADS(t); Thread *prim_t = pt[0]; /* thread for primary phase */ Thread *sec_t = pt[1]; /* thread for secondary phase */ real T_prim = C_T(c,prim_t); /* primary phase temperature (K) */ real C2=C_YI(c, prim_t, 0); /* mass fraction of P in primary phase*/ zone_ID = THREAD_ID(t); if (zone_ID == FLUID_ZONE_ID) *rr = C2*Pre*exp(-E1/(R*T_prim)); else *rr=0.; } PS：C2应该是摩尔分数，这里还是质量分数没有改过来，依然没有结果。 |
回复此楼» 猜你喜欢
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有9人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有9人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有9人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有8人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有10人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有12人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有11人回复
2楼2018-05-30 16:43:32
dzsz99185
木虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1641.1
- 散金: 10
- 帖子: 74
- 在线: 91.1小时
- 虫号: 1847697
- 注册: 2012-06-05
- 性别: GG
- 专业: 有色金属冶金
3楼2018-06-13 08:55:50
