E h p 1 (u 2 v2 )
其中 h为混合气体单位质量焓。各组分的质量反应生成率由基元反应的动力学模型给出:
式中, Mi 为第 i 组分的摩尔质量, Rfj 、 Rbf 分别代表摩尔浓度表示的第 j 个基元反应的正反 应和逆反应的速率。
2.3 化学反应模型
本文反应混合物为燃料 H2 和氧化物 Air(O2+3.76N2),化学反应模型为有限速率模型,反 应速率常数通过 Arrhenius 公式计算,考虑碰撞效率和逆反应,采用 7 组分 8 步基元反应,见 表 2.1。采用时间算子分裂算法处理化学反应刚性问题。来`自^751论*文-网www.751com.cn
表 2.1 H2 与 Air(O2+3.76N2)7 组分 8 步基元反应
Number Reaction A n E
1 H2+M=H+H+M 5.5E18 -1.0 432219.918
2 O2+M=O+O+M 7.2E18 -1.0 493352.76
3 H2O+M=H+OH+M 5.2E21 -1.5 493585.552
4 OH+M=O+H+M 8.5E18 -1.0 422600.62
5 H2O+O=OH+OH 5.8E13 0.0 75316.526
6 H2O+H=OH+H2 8.4E13 0.0 84104.424
7 H+O2=OH+O 2.2E14 0.0 70294.87
8 H2+O=OH+H 7.5E13 0.0 46442.004
注:组分:H , O , OH ,H2 , O2 , H2O ,N2
所有基元反应均为可逆反应 单位:Joules,mol,cm3,
其中 M 为第三碰撞体,起催化作用。A 为指前因子,n 为基元反应级数,E 为反应活化能。
2.4 边界条件
该计算域的外圆为预混燃料入口边界,边界上的每个网格由网格单元处的压力 p 决定:
(1)当 p p0
时,流动处于阻塞状态,喷注速度为
Fluent圆盘结构下旋转爆震波传播特性研究(3):http://www.751com.cn/jixie/lunwen_79344.html