50 UDF中阿伦尼乌斯公式

Fluent中化学反应速率不是通过源项定义的，对于体积反应速率和多相非均相反应速率可以通过以下定义：

• 体积反应速率

#include "udf.h"
#define K1 2.0e-2
#define K2 5.
DEFINE_VR_RATE(vol_reac_rate,c,t,r,mole_weight,species_mf,rate,rr_t)
{
    real s1 = species_mf[0];
    real mw1 = mole_weight[0];
    if (FLUID_THREAD_P(t) && THREAD_VAR(t).fluid.porous)
    *rate = K1*s1/pow((1.+K2*s1),2.0)/mw1;//可定义指前因子和活化能
    else
    *rate = 0.;
    *rr_t = *rate;
}

• 多相非均相反应

DEFINE_HET_RXN_RATE (name,c,t,r,mw,yi,rr,rr_t)

DEFINE_HET_RXN_RATE(user_evap_condens_react, c, t, hr, mw, yi, rr, rr_t)
{
    Thread **pt = THREAD_SUB_THREADS(t);
    Thread *tp = pt[0];
    Thread *ts = pt[1];
    int i;
    *rr = 0.0;
    .......
    *rr = area_density * flux_evap ;//此处定义相间反应速率
}