24小时热门版块排行榜    

查看: 440  |  回复: 1

ereklin

金虫 (职业作家)

[求助] 固固反应模拟,关于转化率怎么定义,怎么用UDF写反应源项

如图所示,加热固体散料图,固固之间反应,温度从外向内一层层传,温度达到一定值就反应吸收反应热,怎么定义转化率;同时怎么用UDF写这个反应源项
固固反应模拟,关于转化率怎么定义,怎么用UDF写反应源项
1.jpg
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ereklin

金虫 (职业作家)

刚接触UDF,不太懂。这是我的UDF程序,计算时中间温度高于了周围(源项有问题),哪儿出错了!
#include "udf.h"
#include "math.h"
#define rq 8314
#define k0  0.1018419
#define ml  236
#define qf  -9.53*pow(10,6)
DEFINE_SOURCE(energy_source, c, t, dS, eqn)        
{
real source;
real Et;
real temp;
real tb;
real zhzhi;
temp=C_T(c,t);
tb=CURRENT_TIME;
if(tb<=900)
        {Et=66691.31;
         }
else
     {  if(tb>12600)
             {Et=79065.39;
             }
         else
             {Et=64605.507+2.478*tb-1.844*pow(10,-4)*pow(tb,2)+6.251*pow(10,-9)*pow(tb,3);
              }
}
if(temp<=1473&&temp>=1273)
  {zhzhi=k0*exp(-Et/(rq*temp));
     source=ml*qf*(-zhzhi*pow(1-(zhzhi*tb)/3,2));
    dS[eqn] = 0;
    }
else
  {source=dS[eqn]=0;
   }
return source;
}
2楼2013-06-07 18:04:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ereklin 的主题更新
信息提示
请填处理意见