| 查看: 4101 | 回复: 7 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
FLUENT UDF化学反应 已有1人参与
|
|||
|
请各位大神指导一下，初学UDF，我想通过UDF编写化学反应速率以及各方程源项，来实现CH4和O2的燃烧反应，编写的代码如下： #include "udf.h" #include "math.h" #define PRE_EXP 2.119e+011 #define ACTIVE 2.027e+08 #define BETA 0.0 real arrhenius_rate(real temp) { return PRE_EXP*pow(temp,BETA)*exp(-ACTIVE/(UNIVERSAL_GAS_CONSTANT*temp)); } 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)/16; return source; } DEFINE_SOURCE(ch4_source,c,t,dS,eqn) { real source; real ss; source=-arrhenius_rate(C_T(c,t))*C_R(c,t)*C_YI(c,t,0); return source; } DEFINE_SOURCE(o2_source,c,t,dS,eqn) { real source; source=-arrhenius_rate(C_T(c,t))*C_R(c,t)*C_YI(c,t,1); return source; } 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,2); return source; } DEFINE_SOURCE(h2o_source,c,t,dS,eqn) { real source; source=arrhenius_rate(C_T(c,t))*C_R(c,t)*C_YI(c,t,3); return source; } 程序中的指前因子，活化能等参数都是fluent中的材料库中数据。选用的是无反应输运模型，编译UDF后，分别添加各方程源相。但是Fluent运行后并没有发生反应，查看发现arrhenius_rate=0，但是这些数据全部都是来自fluent?那么Fluent又是如何计算arrhenius_rate的呢？后来我认为将arrhenius_rate=1000.0，但是还是没有反应？不知道是怎么回事？请大神指点一下啊？模型说明在附件里面。 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : tut13.pdf
2015-11-18 16:27:09, 1.08 M
» 收录本帖的淘帖专辑推荐
 Fluent技巧 |
» 猜你喜欢
284求调剂
已经有10人回复
-
一志愿山东大学药学学硕求调剂
已经有4人回复
-
07化学280分求调剂
已经有4人回复
-
298-一志愿中国农业大学-求调剂
已经有12人回复
-
求材料，环境专业调剂
已经有3人回复
-
335求调剂
已经有5人回复
-
求调剂
已经有7人回复
-
一志愿吉大化学322求调剂
已经有4人回复
-
环境学硕288求调剂
已经有8人回复
-
341求调剂(一志愿湖南大学070300)
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- fluent udf
已经有4人回复
- fluent udf
已经有0人回复
- fluent udf
已经有2人回复
- FLUENT udf
已经有2人回复
- FLUENT udf
已经有0人回复
- fluent UDF
已经有2人回复
5楼2015-11-24 13:29:22
nwj0502
至尊木虫 (知名作家)
- 应助: 5 (幼儿园)
- 金币: 11135.3
- 散金: 6315
- 红花: 38
- 帖子: 9032
- 在线: 1320.5小时
- 虫号: 771979
- 注册: 2009-05-16
- 性别: GG
- 专业: 物理化学
2楼2015-11-23 16:25:24
3楼2015-11-24 10:26:43
4楼2015-11-24 10:34:53
