各位老師，我的模型是一段長方形的管子，中間是多孔介質，里面沒有設置reaction，模擬化學反應，反應方程為：4nh3+6NO=5N2+6h20
進口設置：NO：0.0989
                  NH3：0.0989
                  h2o(l):0
出口設置：為設置，全為0
開了能量方程和組分運輸模型，組分運輸模型里面沒有勾選Volumetric，udf編寫如下，反應用源相來表示，請問這段程序哪里有問題呢，有氮氣生成，就是沒有水，還有反應過程不太正常，瞬間就反應完了。謝謝各位老師！
#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 NO   0
#define NH3  1
#define H20  2
#define N2   3
#define NUM_SPECS 4 
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,0)/30;
        ds[eqn]= delt_h*PRE_EXP*exp(-ACTIVE/(UNIVERSAL_GAS_CONSTANT*temp)) *ACTIVE*(1/temp) *(1/temp) *(1/UNIVERSAL_GAS_CONSTANT)*C_R(c,t)*C_YI(c,t,0)/30;
        return source;
}

DEFINE_SOURCE(no_source,c,t,dS,eqn)
{
        real source; 
        source = 6*(-arrhenius_rate(C_T(c,t)))*C_R(c,t)*C_YI(c,t,0)/30;
        ds[eqn]= -6*PRE_EXP*exp(-ACTIVE/(UNIVERSAL_GAS_CONSTANT*temp)) *ACTIVE*(1/temp) *(1/temp) *(1/UNIVERSAL_GAS_CONSTANT)*C_R(c,t)*C_YI(c,t,0)/30;
        return source;
}

DEFINE_SOURCE(nh3_source,c,t,dS,eqn)
{
        real source;
        source = 4*(-arrhenius_rate(C_T(c,t)))*C_R(c,t)*C_YI(c,t,1)/17;
        ds[eqn]= -4*PRE_EXP*exp(-ACTIVE/(UNIVERSAL_GAS_CONSTANT*temp)) *ACTIVE*(1/temp) *(1/temp) *(1/UNIVERSAL_GAS_CONSTANT)*C_R(c,t)*C_YI(c,t,0)/17;
        return source;
}

DEFINE_SOURCE(h20_source,c,t,dS,eqn)
{
        real source;
        source = 6*(arrhenius_rate(C_T(c,t)))*C_R(c,t)*C_YI(c,t,2)/18;
        ds[eqn]= 6*PRE_EXP*exp(-ACTIVE/(UNIVERSAL_GAS_CONSTANT*temp)) *ACTIVE*(1/temp) *(1/temp) *(1/UNIVERSAL_GAS_CONSTANT)*C_R(c,t)*C_YI(c,t,0)/18;
        return source;
}
DEFINE_SOURCE(n2_source,c,t,dS,eqn)
{
        real source;
        source = 5*(arrhenius_rate(C_T(c,t)))*C_R(c,t)*C_YI(c,t,3)/28;
        ds[eqn]= 5*PRE_EXP*exp(-ACTIVE/(UNIVERSAL_GAS_CONSTANT*temp)) *ACTIVE*(1/temp) *(1/temp) *(1/UNIVERSAL_GAS_CONSTANT)*C_R(c,t)*C_YI(c,t,0)/28; 
        return source;
}