24小时热门版块排行榜

返回列表

一只小火车

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1
帖子: 6
在线: 5.6小时
虫号: 5549266
注册: 2017-02-07
性别: GG
专业: 多相流热物理学

[求助] 气液两相流非均相反应的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应该是摩尔分数，这里还是质量分数没有改过来，依然没有结果。

回复此楼