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。采用时间算子分裂算法处理化学反应刚性问题。来`自^优尔论*文-网www.youerw.com

表 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

时，流动处于阻塞状态，喷注速度为