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 ;//此处定义相间反应速率
}