| 查看: 705 | 回复: 1 | ||
[求助]
udf指点 已有1人参与
|
|
我这是从网上看的一段udf的代码 为啥我把它弄进fluent 中 成功编译了但是 在算的时候为啥都没有反应，生成物的含量都为0，反应物没有变化都 囊中羞涩求大神指点 #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)); } #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; /*下面是甲烷生成焓 单位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) |
回复此楼» 猜你喜欢
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
物理学I论文润色/翻译怎么收费?
已经有181人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有0人回复
2楼2014-05-06 22:00:08
