3. 燃烧反应动力学基础

在Fluent中做燃烧仿真，热力学决定“能不能烧、烧到多热”，而动力学决定“烧得多快、在哪儿烧、火焰有多厚”。如果只给热力学平衡而不给反应速率，求解器只能得到“最终会变成什么”，却无法得到“在多长时间、多长空间内变成什么”。本节围绕反应速率、Arrhenius定律、总包反应与详细机理的差别，以及Fluent中反应模型的选择依据展开。

一、反应速率：从基元反应到宏观消耗

对任一基元反应，例如 H + O2 → OH + O，其反应速率遵循质量作用定律：

r = k · [H] · [O2]        （单位：kmol/(m³·s)）

其中 [H]、[O2] 为体积摩尔浓度，k 为反应速率常数。若反应可逆，净速率还要减去逆反应项 k_b·[OH]·[O]。对非基元的总包反应，如 CH4 + 2O2 → CO2 + 2H2O，其速率不能由化学计量式直接写出，必须通过实验拟合得到，且反应级数往往不是整数。

在Fluent中，反应速率最终以源项形式进入组分输运方程：

∂(ρYi)/∂t + ∇·(ρvYi) = -∇·Ji + Ri

其中 Ri 是第 i 组分的净生成/消耗速率，由所选反应模型计算。因此，反应模型的选择，本质上就是选择“如何计算 Ri”。

二、Arrhenius定律：温度如何控制反应快慢

反应速率常数 k 随温度的变化由Arrhenius定律描述：

k = A · T^n · exp(-Ea / (R·T))

式中 A 为指前因子，n 为温度指数（多数反应取0），Ea 为活化能，R 为通用气体常数（8.314 J/(mol·K)）。指数项 exp(-Ea/RT) 是核心：Ea 越大，温度对速率的敏感性越强。例如 Ea 从 100 kJ/mol 提高到 200 kJ/mol，在 1500 K 下速率常数会下降约 6 个数量级。

在Fluent的Arrhenius设置中，需要为每个反应输入 A、n、Ea。常见做法是：先查文献或机理文件（如GRI-Mech、CHEMKIN格式），再核对单位。Fluent默认 A 的单位与反应级数有关，例如二级反应为 m³/(kmol·s)，若直接照搬文献中的 cm³/(mol·s)，必须换算，否则速率会差 10⁶ 倍以上。

三、总包反应与详细反应机理

总包反应（Global Reaction）把复杂燃烧过程压缩为一到几步，例如常用的两步甲烷总包机理：

CH4 + 1.5 O2 → CO + 2 H2O
CO + 0.5 O2 → CO2

它的优点是计算量小、参数少、容易收敛，适合工程上关注温度场、速度场和主要组分分布的场合。缺点是它无法预测CO、NOx、中间自由基的准确生成，也无法反映点火延迟和熄火极限。

详细反应机理（Detailed Mechanism）则包含数十到数千个基元反应，例如GRI-Mech 3.0含53种组分、325个反应，能较好预测点火、火焰传播和污染物。代价是刚性极强、计算成本高，对网格和求解器稳定性要求高。

介于两者之间的是骨架机理（Skeletal）和简化机理（Reduced），通过保留关键路径、剔除次要反应来折中精度与成本。选择时可用“组分数量—反应数量—目标精度”三者权衡。

四、Fluent中反应模型的选择依据

Fluent提供的主要燃烧模型包括：Species Transport（含有限速率、EDC、层流有限速率）、Non-Premixed（PDF）、Partially Premixed、Premixed以及PDF输运等。选择时可按以下逻辑判断：

• 火焰类型：燃料与氧化剂预先混合再进入燃烧区，选预混或部分预混模型；两者分开进入并靠扩散混合燃烧，选Non-Premixed或PDF。
• 精度目标：只关心温度场和主要组分，可用总包反应+有限速率；需要CO、NOx、点火延迟，必须用详细机理+EDC或PDF输运。
• 湍流—化学相互作用：湍流时间尺度与化学时间尺度之比决定用EDC、PDF还是有限速率。Da数大、化学快，PDF更合适；Da数小、化学慢，有限速率更合适。
• 计算资源：详细机理+EDC在三维网格上计算量极大，工程上常先用总包或骨架机理试算，再局部加密或改用简化机理。

五、一个可复现的示例判断

以甲烷/空气非预混射流火焰为例，若目标是观察火焰长度和温度分布，可先采用Non-Premixed + 平衡PDF，快速得到温度场；若发现CO预测偏低，再切换到Species Transport + EDC + GRI-Mech 3.0，并适当加密火焰区网格。判断是否收敛，除了残差，还应监测出口平均温度、CO质量分数和火焰面位置是否稳定。若温度场在100步内波动小于1%，且组分分布不再整体漂移，可认为统计稳定。

总之，反应速率是燃烧仿真的“发动机”，Arrhenius定律是它的“油门曲线”，总包与详细机理是不同精度的“燃料”，而Fluent模型选择则是把这三者匹配到具体火焰类型和工程目标上的决策过程。理解这四者的关系，才能在后续设置中做出有依据的选择。