ÎÒÓõļ×ÍéµÄȼÉÕUDF
ζȱ仯·¶Î§Ì«Ð¡ÁË26.5µ½26.8¶È£¬´íÎóÔÚÄÜÁ¿·½³ÌµÄ±àд£¿£¿
#include "udf.h"
#define PRE_EXP 1.8e+08
#define ACTIVE 1.3e+08
#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;
delt_h = -7.489518e+07;
source = delt_h*(-arrhenius_rate(C_T(c,t)))*C_R(c,t)*C_YI(c,t,CH4)/16;
return source;
}
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);
return source;
}
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);
return source;
}
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);
return source;
}
DEFINE_SOURCE(h2o,c,t,dS,eqn)
{
real source;
source = 2*(arrhenius_rate(C_T(c,t)))*C_R(c,t)*C_YI(c,t,CH4);
return source;
}
![UDFζÈÔÆÍ¼±ê³ß²»ºÏÀíÈçºÎ´¦Àí]()
UDF.png
![UDFζÈÔÆÍ¼±ê³ß²»ºÏÀíÈçºÎ´¦Àí-1]()
mm.png |