fluent关于化学反应的udf编写
合集下载
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
放热速度 单位 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);/*计算甲烷的消耗速度
udf的代码如下includeudfh动力学参数常量definepreexp18e08指前因子单位1sdefineactive13e08活化能jkmoldefinebeta00温度指数realarrheniusraterealtempreturnpreexppowtempbetaexpactiveuniversalgasconstanttemp
单位 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;
}
/*定义 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);/*计算 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;/*计算甲烷的
位 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; }
{
return PRE_EXP*pow(temp,BETA)*exp(-ACTIVE/(UNIVERSAL_GAS_CONSTANT*temp));
}
#define CH4 0
#define O2 1
#define CO2 2
#define H2O 3
#define N2 4
#define NUM_SPECS 5
UDF 的代码如下
#include "udf.h"
/*动力学参数常量*/
#define PRE_EXP 1.8e+08 /*指前因子 单位 1/s*/
#define ACTIVE 1.3e+08 /*活化能 J/kmol*/
#define BETA 0.0
/*温度指数*/
real arrhenius_rate(real temp)