| 查看: 1384 | 回复: 2 | ||
| 【悬赏金币】回答本帖问题，作者shangsiyi2将赠送您 50 个金币 | ||
[求助]
DEFINE_MASS_TRANSFER udf编译求助
|
||
|
各位专家，我在做一个气液相间传质的模拟，自己写了个udf，但是编译不过去，求各位帮忙看看。 代码如下： #include "udf.h" DEFINE_MASS_TRANSFER(o_lbe_transfer,cell,thread,from_index,from_species_index,to_index,to_species_index) { Thread *o, *lbe; o = THREAD_SUB_THREAD(thread, from_index); lbe = THREAD_SUB_THREAD(thread, to_index); real k_pq;/*mass transfer coefficient*/ real sh_q = 20.;/*sherwood number*/ real diff_q;/*diffusivity of o2 in lbe*/ real diam = C_PHASE_DIAMETER(cell, o);/*bubble diameter*/ real vof_o = C_VOF(cell, o); real vof_lbe = C_VOF(cell, lbe); real area_intf;/*interfacial area*/ real m_gl; real YI_SAT;/*soluability,wt%*/ real YI_E;/*mass fraction of o2 in lbe corresponding to partial pressure*/ real R_I; real x_o;/*molar fraction of o2 in ar*/ real ppress;/*partial pressure*/ real R; int i = 0;/*species index*/ /*caculating mass transfer coefficient*/ diff_q = 0.0239*exp(-43073. / 8.314 / C_T(cell, lbe))*0.0001; k_pq = sh_q*diff_q / diam; /*udf for specifying interfacial area*/ area_intf = 6.*vof_o*vof_lbe / diam; /*mass transfer rate*/ R = C_YI(cell, lbe, i)*C_R(cell, lbe);/* mass concentration */ R_I = 0.0; m_gl = 0.0; YI_SAT = POW(10, 2.52 - 4803 / C_T(cell, lbe)); x_o = 5.*C_YI(cell, o, i) / (4. + C_YI(cell, o, i));/*convert mass fraction to molar fraction*/ ppress = x_o*(C_P(cell, o)+101325);/*total pressure*molar fraction*/ YI_E = sqrt(ppress * 100000. / exp(13.558 - 32005 / C_T(cell, lbe))); if (YI_E > YI_SAT) { R_I = YI_SAT*C_R(cell, lbe); } else { R_I = YI_E*C_R(cell, lbe); } m_gl = k_pq*area_intf*(R_I - rho); return (m_gl); } 模型如下：  |
回复此楼» 猜你喜欢
颜宁:基础研究应顶天立地,立地即从需求端作为问题的出发点,立项背景就是面向真实需求
已经有3人回复
-
发现督导在听课，需要跟TA招呼示意吗？有同学看手机、课件有疏漏算大问题吗？
已经有4人回复
-
催稿后秒拒
已经有6人回复
-
前同事拿我以前未中的国自然标书申请今年的国自然，现在基金委发来调查函
已经有10人回复
-
想问一下有没有需要文章但是没时间带学生的研究生导师。
已经有17人回复
-
298求调剂
已经有9人回复
-
申请省自然科学基金，研究区能否是省外区域
已经有13人回复
-
国内科研环境差导致国内科研水平差的根源
已经有3人回复
-
CSC在职申请求问
已经有7人回复
-
家乡二本高校/沿海传统私企，如何抉择？
已经有15人回复
gyctju
金虫 (正式写手)
- 应助: 314 (大学生)
- 金币: 1131.7
- 散金: 227
- 红花: 48
- 沙发: 2
- 帖子: 987
- 在线: 426.1小时
- 虫号: 1569491
- 注册: 2012-01-08
- 专业: 分离过程
2楼2020-04-02 11:29:44
3楼2022-01-03 00:37:55
30