版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

981308758

新虫 (小有名气)

应助: 3 (幼儿园)
金币: 326.1
散金: 42
红花: 2
帖子: 62
在线: 41.7小时
虫号: 1690698
注册: 2012-03-14
性别: MM
专业: 绿色有机化学

[求助] 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)

回复此楼

» 猜你喜欢

做一只敬业的虫子

1楼 2014-04-23 16:03:05

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

yingwuqiqi

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 880.9
帖子: 2
在线: 1.5小时
虫号: 2066461
注册: 2012-10-16
专业: 流体力学

【答案】应助回帖

★ ★ ★
981308758: 金币+3, ★有帮助, 这个还是知道的依然谢谢你问题已经解决 2014-05-07 15:33:17

编译成功还要在fluent中设置的，在cell condition/fluid/source term中设置，要不然fluent怎么知道你需要使用UDF？

赞一下