|
|
刚接触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;
} |
|