扩散火焰和混合分数
扩散火焰和混合分数的理论基础
概要
老师，扩散火焰是什么样的燃烧形态？
扩散火焰（非预混合火焰）是燃料和氧化剂分别供应，混合与燃烧同时进行的燃烧形态。蜡烛火焰、柴油发动机、燃气轮机燃烧器均属此类。火焰形成于燃料和氧化剂以化学当量比相遇的面（化学计量面）。
预混合火焰在本质上是不同的。
是的。在预混合火焰中，燃烧速度是控制参数，但在扩散火焰中，混合速度成为限速步骤。利用这个特性，用混合分数 $Z$这个变量来描述火焰结构是Burke-Schumann理论的核心。
混合分数
请告诉我混合分数 $Z$ 的定义。
混合分数是一个保守标量，表示"流体单元中含有多少来自燃料流的质量"。广泛使用Bilger定义。
其中$Y_F$是燃料质量分数，$Y_O$是氧化剂质量分数，下标1表示燃料流，下标2表示氧化剂流的值。$s$是质量基础的化学计量比。
$Z=0$是纯氧化剂，$Z=1$是纯燃料，是吧？
正是这样。化学计量混合分数$Z_{st}$定义为$Z_{st} = Y_{O,2}/(s\,Y_{F,1} + Y_{O,2})$，对于甲烷/空气，$Z_{st} \approx 0.055$。火焰位于$Z = Z_{st}$的等值面上。
Burke-Schumann解
Burke-Schumann解是什么？
当化学反应速度无限快（$Da \to \infty$）时，燃料和氧化剂在火焰面上瞬时反应，温度和化学物种只是混合分数的函数。
这就是Burke-Schumann解，温度在$Z_{st}$处达到峰值，呈三角形轮廓。这个解也是小尺度涡模型的出发点。
混合分数的输运方程
混合分数 $Z$ 遵循无源项（不依赖化学反应）的输运方程。
因为没有源项，所以不需要知道化学反应的细节就能求解$Z$的场。
这正是混合分数方法的主要优势。我们将化学反应的复杂性压入到$Z$空间的小尺度涡方程中，CFD只需要输运$Z$及其方差$\widetilde{Z''^2}$。
扩散火焰的理论基于混合分数概念的全部基础。
没错。混合分数可以说是非预混合燃烧CFD中最重要的变量。
蜡烛火焰是"扩散炎"的教科书——Burke-Schumann解的局限性
1928年Burke和Schumann提出的扩散炎分析解基于"无限反应速度"的假设，即燃料和氧化剂的混合瞬间完成。蜡烛火焰在这个模型下的再现效果非常好。然而，当在实际燃气轮机燃烧器中使用同一假设时，NOx生成量可能比实测高10倍以上。虽然混合应该是"限速步骤"，但在高温区域，反应速度也是有限的。这就是为什么现场工程师说"Burke-Schumann解对概念理解最好，但对设计没用"的原因。
扩散火焰和混合分数的数值计算方法
数值方法的详细说明
如何在CFD中实现基于混合分数的燃烧模型？
在非预混合燃烧模型中，求解$Z$和$\widetilde{Z''^2}$（混合分数方差）的输运方程，从查找表中获取化学反应信息。
输运方程
乱流场中Favre平均混合分数$\widetilde{Z}$和方差$\widetilde{Z''^2}$的方程如下。
方差方程最后的项$\widetilde{\chi}$是什么？
标量消散率（scalar dissipation rate）。定义为$\chi = 2D|\nabla Z|^2$，表示混合分数的微观结构消散得有多快。在乱流条件下，建模为$\widetilde{\chi} = C_\chi \frac{\varepsilon}{k}\widetilde{Z''^2}$（$C_\chi \approx 2.0$）。
PDF（概率密度函数）
请解释PDF的作用。
在乱流场中，$Z$在单元内波动，仅用平均值无法正确表达火焰结构。假设$Z$的概率密度函数$P(Z)$，通过积分求平均温度和平均化学物种。
PDF形状普遍采用$\beta$函数分布。由$\widetilde{Z}$和$\widetilde{Z''^2}$这两个参数唯一确定$\beta$分布的形状。
查找表的构建
表是怎样生成的？
预先求解小尺度涡方程或化学平衡，生成以$Z$和$\chi_{st}$（化学计量散逸率）为参数的温度和化学物种查找表。这个表用PDF积分后，在CFD运行时参照。
| 表变量 | 维数 | 用途 |
|---|---|---|
| $\widetilde{Z}$, $\widetilde{Z''^2}$ | 2D | 平衡化学、薄火焰 |
| $\widetilde{Z}$, $\widetilde{Z''^2}$, $\widetilde{\chi_{st}}$ | 3D | 定常小尺度涡 |
| $\widetilde{Z}$, $\widetilde{Z''^2}$, $\widetilde{C}$ | 3D | FGM (加进度变量) |
在Fluent中如何设置？
在Fluent中选择Models > Species > Non-Premixed Combustion，以CHEMKIN格式导入反应机制生成PDF表。表的分辨率（$Z$方向分割数）最低64点，建议128点。
混合分数+PDF表的方法巧妙地在事前求解了化学反应。
是的。3D CFD的运行时不必求解化学反应，所以即使使用详细化学反应，计算成本基本不变，这是最大的优势。
如何求解标量消散率χ——扩散炎模型实现中最有争议的变量
当实现扩散炎的小尺度涡模型时，最令实现者困扰的是标量消散率χ（chi）的计算。χ表示炎面的混合激烈程度（乱流对燃料-氧化剂的拉伸率），χ越高，火焰越容易熄灭。但在乱流CFD中直接计算χ的公式有多种，根据Favre平均和模型常数的选择，结果会改变。特别是当再循环区（喷嘴后流）χ接近0的区域和剪切层中χ高的区域共存时，空间分布的计算精度直接关系到NOx和烟尘的预测。"如何计算χ"已经成为供应商之间的商业秘密。
扩散火焰和混合分数的实务应用
实践指南
请教我扩散火焰CFD分析的实务步骤。
这是燃气轮机燃烧器、锅炉等非预混合燃烧分析的典型流程。
分析流程
1. 定义燃料成分 -- 天然气、煤油、氢气等。多成分燃料的情况下，定义代用混合物
2. 反应机制选择和PDF表生成 -- Fluent的Non-Premixed Combustion模型中导入CHEMKIN格式
3. 网格生成 -- 在混合区（燃料喷射附近）确保充分的分辨率
4. 乱流模型选择 -- k-$\varepsilon$系在旋转流中精度不足。推荐SST k-$\omega$或RSM
5. 边界条件设置 -- 燃料入口($Z=1$)和空气入口($Z=0$)。流量从总当量比反推
6. 辐射模型 -- DO（离散纵标法）或P-1。用WSGM设置吸收系数
网格设计要点
混合分数模型在网格上有什么注意事项？
混合精度决定了模型的成败。特别要细化以下区域。
| 区域 | 推荐单元大小 | 原因 |
|---|---|---|
| 燃料喷嘴出口 | 喷嘴直径的1/10以下 | 喷流初期混合的分辨 |
| 混合层（$0.01 < Z < 0.1$） | 1-2 mm | 火焰带的分辨 |
| 旋流涡的核心 | 2-3 mm | 再循环区的捕捉 |
| 壁面附近 | y+ < 1 (LES) / 30-300 (RANS) | 壁面热传导 |
如果燃料喷嘴是1 mm，就需要0.1 mm网格。
是的。所以燃气轮机燃烧器的严格LES通常有数千万个网格。RANS可以用数百万网格，但涡旋诱导涡破裂（VBB）的预测精度就会受限。
常见故障和对策
| 现象 | 原因 | 对策 |
|---|---|---|
| 火焰温度超过绝热值 | PDF表Z方向分割太粗 | 表分辨率改为128点以上 |
| 排出口存在Z>0的未燃燃料 | 混合不良（网格粗糙） | 细化喷射附近网格 |
| NOx是实验值的2倍以上 | thermal NOx的温度依赖性陡峭 | 先改善温度场精度（辐射模型） |
| 火焰长度不匹配 | 乱流Sc数不合适 | $Sc_t = 0.7$（改变默认的0.85） |
乱流Schmidt数可以改吗？
$Sc_t$是模型常数，在燃烧流中通常0.7左右与实验相符。Fluent的默认值0.85扩散量偏少。自由喷流火焰中，将其改为0.7可以改善火焰长度的一致性。
虽然混合分数模型设置简单，但混合分辨率决定了成败。
这就是"混合控制一切（Mixing controls everything）"这句非预混合燃烧的铁律。
燃气轮机燃烧器的"稀薄燃烧"——从扩散炎的历史转变
20世纪60~80年代，燃气轮机燃烧器主流是强调稳定性的扩散炎设计。虽然在燃烧过程中混合燃料和空气确实稳定，但高温区NOx激增。面对排放规制强化，从1990年代开始，"预混合稀薄燃烧（Lean Premixed, LPM）"成为主流。GE的DLN（Dry Low NOx）喷嘴和三菱重工的DLEC喷嘴就是典型，设计需要预混合小尺度涡模型不可或缺。但稀薄预混合容易吹灭，部分负荷下的稳定性确保仍是困扰工程师的实际课题。
扩散火焰和混合分数的软件对比
商业工具对比
请告诉我扩散火焰模型的工具支持情况。
非预混合燃烧（混合分数方法）已被主要CFD工具几乎全部支持。
| 工具 | 模型名称 | PDF表生成 | 多燃料流支持 | LES支持 |
|---|---|---|---|---|
| ANSYS Fluent | Non-Premixed Combustion | 自动（CHEMKIN输入） | 支持2个燃料流 | 支持 |
| STAR-CCM+ | Flamelet/FGM | 自动（DARS联动） | 支持多流 | 支持 |
| OpenFOAM | flameletFoam等 | 手动（用Cantera等事前生成） | 定制 | 支持 |
| CONVERGE | Flamelet model | SAGE联动 | 支持 | 支持 |
ANSYS Fluent
请讲解Fluent的设置步骤。
1. 选择Species > Non-Premixed Combustion
2. 用CHEMKIN格式定义燃料成分（Rich Flammability Limit一侧）
3. 自动计算PDF表（Number of Grid Points: 推荐128）
4. 选择绝热/非绝热（非绝热情况下，加入Enthalpy offset表）
5. 在边界条件中为各Inlet设置$Z$值（燃料入口: Z=1，空气入口: Z=0）
STAR-CCM+
在STAR-CCM+中怎样操作？
STAR-CCM+将其实现为Flamelets/FGM模型。与DARS库联动自动生成小尺度涡表。能够添加Progress Variable的FGM扩展功能强大，可以重现局部熄火。
OpenFOAM
OpenFOAM是标准配置吗？
标准发行版中有基础的小尺度涡模型，但实务级别广泛使用社区版flameletFoam（TU Darmstadt开发）。PDF表用Cantera事前生成，作为OpenFOAM的lookup table导入。
选择指南
各工具分别适用于哪些场景？
不同用途有最优的工具。
是的。非预混合燃烧模型本身在各工具中以相同的理论基础实现，所以网格生成的易用性、LES支持、喷雾模型联动等周边功能是差异所在。
Fluent的"Non-Premixed Combustion"——扩散炎模型实现成为黑箱的原因
使用ANSYS Fluent进行扩散炎计算时使用的"Non-Premixed Combustion"模型，内部运行PDF（概率密度函数）表和小尺度涡计算。但这个表生成流程在默认设置下，散逸率范围和网格点数在用户意识之外决定。特别是在高压条件（航空发动机、高压锅炉）下，实际燃烧条件可能超出表的范围。STAR-CCM+在同等功能中参数可视化程度更高，但因为参数多，初学者容易困惑。无论哪个，"理解默认值的含义后再使用"是前提。
扩散火焰和混合分数的前沿研究
前沿话题和研究动向
扩散火焰的研究发展到什么阶段了？
扩散火焰的CFD建模是成熟领域，但活跃研究仍在继续。特别是(1) 局部熄火和重新着火、(2) 差压喷雾火焰、(3) 氢气扩散火焰这3个领域受关注。
局部熄火和重新着火
混合分数模型无法处理熄火现象吗？
标准的平衡PDF模型无法重现熄火。当标量消散率超过临界值$\chi_q$时火焰局部消灭的现象需要用非定常小尺度涡模型或FGM（小尺度涡生成流形）加进度变量来处理。
熄火条件可用下列Damkohler数表示。
当$Da < Da_{\text{crit}} \approx 1$时发生熄火。
条件矩闭合法 (CMC)
CMC是什么？
CMC是以混合分数为条件变量，输运化学物种条件平均的方法。由Klimenko和Bilger独立提出。不依赖于PDF表假设，对局部熄火和抬升火焰的预测精度更高。
| 方法 | 局部熄火 | 计算成本 | 实现难度 |
|---|---|---|---|
| 平衡PDF | 不可能 | 低 | 容易 |
| 定常小尺度涡 | 有限 | 低 | 容易 |
| FGM | 可能 | 中等 | 中等 |
| 非定常小尺度涡 | 可能 | 中等 | 中等 |
| CMC | 良好 | 高 | 高 |
| 输运PDF | 最优 | 最高 | 高 |
氢气扩散火焰
氢气的扩散火焰有特殊考虑吗？
氢气的质量扩散系数比其他气体大4~5倍。因此等Le数假设（$Le=1$）不成立，差扩散（differential diffusion）效果明显。普通混合分数模型隐含假设$Le=1$，所以氢气火焰需要修正。
具体如何处理？
面向氢社会，氢气扩散火焰的建模变得越来越重要。
正是这样。氢的Le数效应影响火焰稳定性和NOx排放，所以直接用甲烷模型会很危险。
激光计测击毁CFD——LIF和PIV揭示扩散炎的"意外真实"
从2000年代平面激光诱导荧光（PLIF）和粒子图像测速（PIV）普及以来，扩散炎CFD的验证方式发生了戏剧性变化。以前"基本相符"的甲烷-空气扩散炎的OH自由基分布实测与CFD相比在形状和峰值上都偏差20~30%。根本原因几乎总是"由于火焰伸长导致的熄灭的低估"。这直接导致了非定常小尺度涡和熄火重着火模型研究的加速。实验计测揭示CFD的缺陷，进而催生下一代模型——这个循环推动着燃烧CFD的最前沿。
扩散火焰和混合分数的故障处理
故障排除
请教我混合分数模型常见的故障。
让我整理一下非预混合燃烧模型特有的故障。
1. PDF表生成错误
现象: Fluent/STAR-CCM+的PDF表计算失败。
原因和对策:
- 反应机制有问题: CHEMKIN格式语法错误。用Cantera转为YAML后验证
- 小尺度涡无法收敛: 散逸率过高导致解发散。降低$\chi_{st}$上限
- 热力学数据不一致: NASA7系数格式错误。检查温度范围
2. 非物理温度
温度有时会超过绝热火焰温度…
原因: 当PDF表$Z$方向离散化太粗时，$\beta$ PDF积分产生插值误差。$Z_{st}$附近特别重要。
对策:
- 提高表分辨率（最低64点，推荐128点）
- 在$Z_{st}$附近集中网格点（Fluent自动执行）
- 非绝热模型中，焓缺陷方向也要充分分割数（最低20点）
3. 火焰形状与实验不符
| 现象 | 可能原因 | 对策 |
|---|---|---|
| 火焰过长 | 乱流混合弱 | 降低$Sc_t$(0.7)，细化网格 |
| 火焰过短 | 数值扩散导致混合过度 | 使用二阶精度格式 |
| 火焰不抬升 | 平衡模型无法再现熄火 | 改用FGM或部分预混合模型 |
| 火焰径向扩展过大 | 乱流模型扩散过度 | 从k-$\varepsilon$改为SST或RSM |
4. 多燃料流问题
燃料从多个地方进入怎样处理？
Fluent的Non-Premixed Combustion模型基本上是1个燃料流+1个氧化剂流的2流体系统。如果有导流燃料和主燃料成分不同，需要定义Secondary Stream或改用部分预混合模型。STAR-CCM+的FGM模型原生支持多燃料流。
混合分数模型的故障大多集中在表的质量和混合分辨率。
正是这样。既然化学反应委托给表，表的质量和CFD中$Z$场的精度就是全部。遵守先在0D/1D验证表的行为后再进3D的步骤。
"混合分数超过1.0了！"——扩散炎CFD的标配错误原因与对策
在扩散炎模拟中，初学者常遇到"混合分数Z超出0~1范围"的问题。Z按定义应该在0（纯氧化剂）~1（纯燃料）范围内，但数值计算时，边界条件不一致或梯度陡峭会导致1.01或-0.02等非物理值出现。这样的话，小尺度涡表的参考范围会超出，温度和密度会急剧变化导致计算发散。处理方法除了再确认边界条件外，还要启用标量或Z的截断处理。Fluent有"Bounded scalar"选项，OpenFOAM有`min/max`约束设置。
相关话题
更详细
错误