| 查看: 442 | 回复: 1 | ||
[求助]
固固反应模拟,关于转化率怎么定义,怎么用UDF写反应源项
|
|
刚接触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













回复此楼