研发埠

标题: 化学udf，没有发生反应 [打印本页]

作者: orgring 时间: 2018-3-19 10:05
标题: 化学udf，没有发生反应
各位老师，我的模型是一段长方形的管子，中间是多孔介质，里面没有设置reaction，模拟化学反应，反应方程为：4nh3+6NO=5N2+6h20
进口设置：NO：0.0989
               NH3：0.0989
               h2o(l):0
出口设置：为设置，全为0
开了能量方程和组分运输模型，组分运输模型里面没有勾选Volumetric，udf编写如下，反应用源相来表示，请问这段程序哪里有问题呢，有氮气生成，就是没有水，还有反应过程不太正常，瞬间就反应完了。谢谢各位老师！
#include "udf.h"

#define PRE_EXP 1.8e+08
#define ACTIVE  1.3e+08
#define BETA 0.0
real arrhenius_rate(real temp)
{
return PRE_EXP*pow(temp,BETA)*exp(-ACTIVE/(UNIVERSAL_GAS_CONSTANT*temp));
}
#define NO 0
#define NH3  1
#define H20  2
#define N2 3
#define NUM_SPECS 4
DEFINE_SOURCE(energy_source,c,t,dS,eqn)
{
      real delt_h, source;
      delt_h = -7.489518e+07;
      source = -delt_h*(-arrhenius_rate(C_T(c,t)))*C_R(c,t)*C_YI(c,t,0)/30;
      ds[eqn]= delt_h*PRE_EXP*exp(-ACTIVE/(UNIVERSAL_GAS_CONSTANT*temp)) *ACTIVE*(1/temp) *(1/temp) *(1/UNIVERSAL_GAS_CONSTANT)*C_R(c,t)*C_YI(c,t,0)/30;
      return source;
}

DEFINE_SOURCE(no_source,c,t,dS,eqn)
{
      real source;
      source = 6*(-arrhenius_rate(C_T(c,t)))*C_R(c,t)*C_YI(c,t,0)/30;
      ds[eqn]= -6*PRE_EXP*exp(-ACTIVE/(UNIVERSAL_GAS_CONSTANT*temp)) *ACTIVE*(1/temp) *(1/temp) *(1/UNIVERSAL_GAS_CONSTANT)*C_R(c,t)*C_YI(c,t,0)/30;
      return source;
}

DEFINE_SOURCE(nh3_source,c,t,dS,eqn)
{
      real source;
      source = 4*(-arrhenius_rate(C_T(c,t)))*C_R(c,t)*C_YI(c,t,1)/17;
      ds[eqn]= -4*PRE_EXP*exp(-ACTIVE/(UNIVERSAL_GAS_CONSTANT*temp)) *ACTIVE*(1/temp) *(1/temp) *(1/UNIVERSAL_GAS_CONSTANT)*C_R(c,t)*C_YI(c,t,0)/17;
      return source;
}

DEFINE_SOURCE(h20_source,c,t,dS,eqn)
{
      real source;
      source = 6*(arrhenius_rate(C_T(c,t)))*C_R(c,t)*C_YI(c,t,2)/18;
      ds[eqn]= 6*PRE_EXP*exp(-ACTIVE/(UNIVERSAL_GAS_CONSTANT*temp)) *ACTIVE*(1/temp) *(1/temp) *(1/UNIVERSAL_GAS_CONSTANT)*C_R(c,t)*C_YI(c,t,0)/18;
      return source;
}
DEFINE_SOURCE(n2_source,c,t,dS,eqn)
{
      real source;
      source = 5*(arrhenius_rate(C_T(c,t)))*C_R(c,t)*C_YI(c,t,3)/28;
      ds[eqn]= 5*PRE_EXP*exp(-ACTIVE/(UNIVERSAL_GAS_CONSTANT*temp)) *ACTIVE*(1/temp) *(1/temp) *(1/UNIVERSAL_GAS_CONSTANT)*C_R(c,t)*C_YI(c,t,0)/28;
      return source;
}

: 登录/注册后可看大图

: 登录/注册后可看大图

: 登录/注册后可看大图

: 登录/注册后可看大图

: 登录/注册后可看大图

: 登录/注册后可看大图

欢迎光临研发埠 (http://bbs.yanfabu.com/)