| 查看: 1501 | 回复: 3 | ||
| 【悬赏金币】回答本帖问题，作者鹿1213将赠送您 5 个金币 | ||
[求助]
Fluent模拟CH4和O2的化学反应出现temperature limit，缩小松弛因子也不行，求大佬指导
|
||
|
用的EDM和编译UDF源项，恳请前辈们看一下会不会是udf的问题？ 我的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) /*定义化学反应速率函数*/ { return PRE_EXP*pow(temp, BETA)*exp(-ACTIVE / (UNIVERSAL_GAS_CONSTANT*temp)); /*pow(a,b)函数是求a的b次方*/ } //定义化学反应物成分以及组分数量 #define CH4 0 #define O2 1 #define CO2 2 #define H2O 3 #define N2 4 #define NUM_SPECS 5 //定义能量方程源项 DEFINE_SOURCE(energy_source, c, t, dS, eqn) { real delt_h, source; /*下面是CH4生成焓 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; /*计算CH4的放热速度 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); /*计算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_source, 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_source, 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; }  捕获.JPG |
回复此楼» 猜你喜欢
法国博士招聘
已经有0人回复
-
法国博士招聘
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有149人回复
-
法国博士招聘
已经有0人回复
-
法国博士招聘
已经有0人回复
-
超导三元凸包图
已经有0人回复
-
最近遇到这样一个问题
已经有10人回复
-
论文接收，散金币了！
已经有175人回复
-
vasp算居里温度
已经有0人回复
-
郑州大学光电器件方向招收一名博士
已经有6人回复
2楼2019-06-26 11:30:44
|
本帖内容被屏蔽 |
3楼2019-11-12 08:25:18
4楼2021-04-07 12:44:46
10