fluent关于化学反应的udf编写
UDF的代码如下#include "udf.h"
/*动力学参数常量*/
#define PRE_EXP 1.8e+08 /*指前因子 单位1/s*/
#define ACTIVE1.3e+08 /*活化能 J/kmol*/
#define BETA 0.0 /*温度指数*/
real arrhenius_rate(real temp)
{
return PRE_EXP*pow(temp,BETA)*exp(-ACTIVE/(UNIVERSAL_GAS_CONSTANT*temp));
}
#define CH40
#define O2 1
#define CO22
#define H2O3
#define N2 4
#define NUM_SPECS 5
/*定义能量方程源项*/
DEFINE_SOURCE(energy_source,c,t,dS,eqn)
{
real delt_h, source;
/*下面是甲烷生成焓 单位J/kmol*/
delt_h = -7.489518e+07;
source = delt_h*(-arrhenius_rate(C_T(c,t)))*C_R(c,t)*C_YI(c,t,CH4)/16;/*计算甲烷的放热速度 单位j/(m3 s)*/
return source;
}
/*定义CH4组分方程源项*/
DEFINE_SOURCE(ch4_source,c,t,dS,eqn)
{
real source;
source = (-arrhenius_rate(C_T(c,t)))*C_R(c,t)*C_YI(c,t,CH4);/*计算甲烷的消耗速度 单位kg/(m3 s)*/
return source;
}
/*定义O2组分方程源项*/
DEFINE_SOURCE(o2_source,c,t,dS,eqn)
{
real source;
source = 2*(-arrhenius_rate(C_T(c,t)))*C_R(c,t)*C_YI(c,t,CH4);/*计算O2的消耗速度 单位kg/(m3 s)*/
return source;
}
/*定义CO2组分方程源项*/
DEFINE_SOURCE(co2,c,t,dS,eqn)
{
real source;
source = (arrhenius_rate(C_T(c,t)))*C_R(c,t)*C_YI(c,t,CH4);/*计算CO2的生成速度 单位kg/(m3 s)*/
return source;
}
/*定义H2O组分方程源项*/
DEFINE_SOURCE(h2o,c,t,dS,eqn)
{
real source;
source = 2*(arrhenius_rate(C_T(c,t)))*C_R(c,t)*C_YI(c,t,CH4);/*计算H2O的生成速度 单位kg/(m3 s)*/
return source;
}
页:
[1]