fluent关于化学反应的udf编写

cora

UDF的代码如下

1. #include "udf.h"
   
2. /*动力学参数常量*/
   
3. #define PRE_EXP 1.8e+08 /*指前因子 单位1/s*/
   
4. #define ACTIVE  1.3e+08    /*活化能 J/kmol*/
   
5. #define BETA 0.0         /*温度指数*/
   
6. real arrhenius_rate(real temp)
   
7. {
   
8.         return PRE_EXP*pow(temp,BETA)*exp(-ACTIVE/(UNIVERSAL_GAS_CONSTANT*temp));
   
9. }
   
10. #define CH4  0  
    
11. #define O2   1  
    
12. #define CO2  2  
    
13. #define H2O  3  
    
14. #define N2   4
    
15. #define NUM_SPECS 5  
    
16. /*定义能量方程源项*/
    
17. DEFINE_SOURCE(energy_source,c,t,dS,eqn)
    
18. {
    
19.         real delt_h, source;
    
20.         /*下面是甲烷生成焓 单位J/kmol*/
    
21.         delt_h = -7.489518e+07;
    
22.         source = delt_h*(-arrhenius_rate(C_T(c,t)))*C_R(c,t)*C_YI(c,t,CH4)/16;/*计算甲烷的放热速度 单位j/(m3 s)*/
    
23.         return source;
    
24. }
    
25. /*定义CH4组分方程源项*/
    
26. DEFINE_SOURCE(ch4_source,c,t,dS,eqn)
    
27. {
    
28.         real source;
    
29.         source = (-arrhenius_rate(C_T(c,t)))*C_R(c,t)*C_YI(c,t,CH4);/*计算甲烷的消耗速度 单位kg/(m3 s)*/
    
30.         return source;
    
31. }
    
32. /*定义O2组分方程源项*/
    
33. DEFINE_SOURCE(o2_source,c,t,dS,eqn)
    
34. {
    
35.         real source;
    
36.         source = 2*(-arrhenius_rate(C_T(c,t)))*C_R(c,t)*C_YI(c,t,CH4);/*计算O2的消耗速度 单位kg/(m3 s)*/
    
37.         return source;
    
38. }
    
39. /*定义CO2组分方程源项*/
    
40. DEFINE_SOURCE(co2,c,t,dS,eqn)
    
41. {
    
42.         real source;
    
43.         source = (arrhenius_rate(C_T(c,t)))*C_R(c,t)*C_YI(c,t,CH4);/*计算CO2的生成速度 单位kg/(m3 s)*/
    
44.         return source;
    
45. }
    
46. /*定义H2O组分方程源项*/
    
47. DEFINE_SOURCE(h2o,c,t,dS,eqn)
    
48. {
    
49.         real source;
    
50.         source = 2*(arrhenius_rate(C_T(c,t)))*C_R(c,t)*C_YI(c,t,CH4);/*计算H2O的生成速度 单位kg/(m3 s)*/
    
51.         return source;
    
52. }

复制代码