预混合火焰模型
预混合火焰的理论基础
概要
老师，预混合火焰与扩散火焰有什么区别？
预混合火焰是燃料和氧化剂在燃烧前充分混合，火焰传播的形态。对应于汽油发动机、燃气轮机的贫燃燃烧器、家用燃气灶。火焰面具有明确的边界，分隔未燃混合气和已燃气体。
火焰传播速度是重要参数呢。
是的。层流燃烧速度 $S_L$ 是预混合火焰的基本参数。甲烷/空气（当量比1.0、常温常压）中 $S_L \approx 0.36$ m/s，氢气/空气中 $S_L \approx 2.1$ m/s。
进行变量 $c$ 和支配方程
预混合火焰的CFD中使用什么变量？
预混合火焰使用进行变量 $c$ 追踪火焰面。$c=0$ 为未燃混合气，$c=1$ 为已燃气体。
其中 $\dot{\omega}_c$ 是反应源项，仅在火焰面附近非零。
源项 $\dot{\omega}_c$ 如何建模？
这是预混合火焰建模的核心。主要有3种方法。
主要的预混合燃烧模型
| 模型 | 原理 | 优点 | 缺点 |
|---|---|---|---|
| G-equation (Level Set) | 用 $G=0$ 等值面追踪火焰面 | 几何学直观 | 火焰面内部结构缺失 |
| TFC (Turbulent Flame Closure) | Zimont 模型。$S_T = A(u'/S_L)^n S_L$ | 实现容易 | $S_T$ 依赖经验公式 |
| FSD (Flame Surface Density) | 火焰面积密度 $\Sigma$ 的输送方程 | 基于物理 | $\Sigma$ 方程的模型常数 |
| c-equation + 反应速度 | $\dot{\omega}_c = \rho_u S_L | \nabla c |$ | 直接 | 火焰厚度的数值问题 |
乱流燃烧速度
乱流中的火焰速度如何变化？
乱流燃烧速度 $S_T$ 随乱流强度 $u'$ 增加。Zimont 关联式被广泛使用。
其中 $\alpha$ 是温度扩散系数，$l_t$ 是乱流的积分尺度，$A$ 是模型常数（$A \approx 0.52$）。
乱流越强，火焰面就变得越皱，见的燃烧速度就增大，是这个意思吗？
正确。乱流使火焰面产生皱褶（wrinkle），增加面积，从而增大单位截面的燃烧速度。这是Damkohler（1940）的经典描述，现代CFD模型都基于这个思想。
预混合火焰与扩散火焰的核心区别就是"火焰面追踪"呢。
正是这样。用进行变量代替混合分率来描述火焰，这是根本的不同。
"火焰厚度"小于1mm——测量层流燃烧速度的难度
预混合火焰理论中的基础量"层流燃烧速度SL"实际上很难测量。有多种实验方法（悬挂化焰法、对向流法、球形传播法等），但同一燃气用不同方法测量时，往往有10～20%的偏差。原因是"如何修正拉伸效应"这一问题上各研究者的看法不同。再者，火焰厚度仅有0.1～1mm，插入温度计测量时计测仪器本身会干扰火焰。这种"测量本身破坏对象"的问题与量子力学的不确定性原理性质不同，但从测量难度的角度看是相似的结构。GRI-Mech 3.0的反应常数验证中使用了这些谨慎得到的测量数据。
预混合火焰的数值计算方法
数值方法的详细
用CFD解预混合火焰的主要数值问题是什么？
最大的问题是火焰厚度的分析。层流预混合火焰的厚度 $\delta_L \approx \alpha/S_L$，甲烷/空气约0.5 mm，氢气/空气约0.2 mm。用RANS的数mm级别网格直接分析是不可能的。
加厚火焰模型（TFM）
怎么解决这个问题？
LES中广泛使用的是加厚火焰模型（Thickened Flame Model，Colin et al., 2000）。把火焰人工加厚，使其能被网格分析。
增加扩散系数 $F$ 倍，减少反应速度 $1/F$ 倍。
这样火焰厚度增加到 $F\delta_L$，但 $S_L$ 保持不变。通常 $F = 5-20$。
但是，加厚火焰会改变与乱流的相互作用吧？
问得好。加厚的火焰无法分析小尺度乱流皱褶。因此引入效率函数 $E$ 进行修正。
Charlette效率函数是代表性模型，形式为 $E = E(\Delta/\delta_L, u'/S_L)$。
Fluent中的实现
在Fluent中怎么设置预混合火焰？
Fluent提供以下模型。
1. 预混合燃烧 (Zimont TFC 模型)：基于c方程。用于RANS
2. 部分预混合燃烧：预混合+非预混合混合模型
3. FGM（Flamelet Generated Manifold）：进行变量 + 混合分率
Zimont TFC模型的设置：
- Models > Species > Premixed Combustion
- Turbulent Flame Speed model：Zimont
- Laminar Flame Speed：输入值或计算值（依赖当量比）
- Flame Stretch Factor：默认0.26
OpenFOAM中的实现
OpenFOAM中呢？
XiFoam 是预混合燃烧求解器。求解火焰皱褶系数 $\Xi$（= $S_T/S_L$）的输送方程。
| 求解器 | 对象 | 模型 |
|---|---|---|
| XiFoam | 预混合压缩性 | $\Xi$-equation |
| reactingFoam + PaSR | 预混合/部分预混合 | Species Transport |
| fireFoam | 火灾 | EDM/扩散火焰 |
加厚火焰模型在OpenFOAM中是标准配置吗？
标准版本中没有，但社区版本（如TFM4OpenFOAM）可用。在燃气轮机LES研究中被广泛应用。
预混合火焰的数值解法核心就是"加厚火焰"这个大胆的思想呢。
正是。TFM在物理上是精妙的技巧，已成为LES预混合燃烧的实际标准。
"进行度变量c"的本质——预混合火焰模型中最困扰的变量
预混合火焰模型的数值实现中，许多人会在反应进行度变量c的定义上出现问题。c表示0（未燃）到1（已燃）的变量，但"用哪种化学物质的质量分率定义c"因模型和研究者而异。有用CO2定义的学派，有用规格化温度的学派，有用多个成分线性组合定义的学派——各种定义会产生微妙的差异。Fluent的默认设置使用生成物的组合，但根据燃料和当量比，使用其他定义有时会获得更高精度。"改变c的定义后突然吻合了"的经验谈在CAE学会中经常被提及。
预混合火焰的实务应用
实践指南
请教预混合火焰分析的实际操作流程。
燃气轮机贫燃燃烧器或SI发动机的预混合火焰分析流程。
分析流程
1. 层流燃烧速度 $S_L$ 的验证 -- 用Cantera计算1D层流预混合火焰，确认当量比、压力、温度依存性
2. 燃烧制度的确认 -- 用Peters图判定条件是否在皱褶火焰还是薄反应区制度
3. 模型选择 -- RANS：Zimont TFC或FGM，LES：TFM或FGM
4. 网格设计 -- 火焰带有充分分辨率（使用TFM时，$F\delta_L$ 至少5个点）
5. 点火方法 -- 补丁点火（设定$c=1$ 区域）或火花点火模型
6. 后处理 -- 评估火焰面位置、$S_T$ 值、回火/吹熄风险
$S_L$ 的设置
$S_L$ 的设置很关键吧？
是的。$S_L$ 偏差10%会导致 $S_T$ 偏差超过10%。准确设置 $S_L$ 作为温度、压力、当量比三变量的函数至关重要。
| 燃料 | $\phi=1.0$, 1atm, 300K 时 $S_L$ [cm/s] | 出处 |
|---|---|---|
| CH4/air | 36 | GRI-Mech 3.0 |
| C3H8/air | 39 | USC Mech II |
| H2/air | 210 | Li et al. |
| NH3/air | 7 | Otomo et al. |
氨气是7 cm/s吗。是甲烷的1/5呢。
因此单独NH3很难维持稳定火焰，主流策略是混合10-30%的H2来改善$S_L$。
常见失败和对策
| 症状 | 原因 | 对策 |
|---|---|---|
| 火焰回火 | $S_T$ 过大 | 降低Zimont常数$A$，细化网格 |
| 火焰吹熄 | $S_T$ 过小或混合不良 | 确认$S_L$设置，检查入口乱流强度 |
| 火焰面数值过厚 | 网格分辨率不足 | 降低TFM的$F$ 值（需结合网格细化） |
| 压力振动（热声） | 火焰-声学耦合 | 设定非反射边界条件，推荐LES |
热声不稳定性也能用CFD预测吗？
用LES可以直接评估压力变动与发热变动的相关（Rayleigh判据）。燃气轮机燃烧振动问题是LES预混合燃烧的重要应用领域。Fluent的压缩性求解器或OpenFOAM的XiFoam可以处理。
预混合火焰分析的出发点就是正确设置$S_L$呢。
不过言。"掌握$S_L$者掌握预混合燃烧CFD"。
"不着火"问题的80%是初始条件——预混合火焰实践中的常见陷阱
预混合火焰CFD实践中最常见的咨询是"计算收敛但一直不着火"问题。Fluent、STAR-CCM+等都会以"无反应流动"状态收敛，用户往往察觉不到。原因大多是两种之一：初始化（Patch）中没有在燃烧区设置高温气体区域，或点火单元的温度没有超过活性化能閾值。现场的定式是"在燃烧区贴2000K以上的小补丁→运行数百步确认发热→本计算"的两阶段点火程序。掌握这一点就能大大缩短初期调试时间。
预混合火焰的软件比较
商用工具比较
请告诉我预混合火焰模型的工具支持情况。
| 工具 | TFC/Zimont | TFM | FGM | G-equation | LES支持 |
|---|---|---|---|---|---|
| Ansys Fluent | 标准配置 | UDF实现 | R2及以后 | 无 | 有 |
| STAR-CCM+ | 配置 | 配置 | 标准推荐 | 无 | 有 |
| OpenFOAM | XiFoam($\Xi$-eq) | 社区版 | 社区版 | 社区版 | 有 |
| CONVERGE | 配置 | 配置 | 配置 | G-equation配置 | 有 |
G-equation只有CONVERGE支持吗？
CONVERGE专为SI发动机配置了G-equation（Level Set法）作为标准配置。用G=0等值面追踪火焰，用于爆燃预测和循环变动评估。其他工具的G-equation是研究级实现。
燃气轮机LES
燃气轮机的LES预混合燃烧用哪种合适？
SI发动机
SI发动机的火焰传播分析呢？
许可证和成本
| 工具 | 预混合模型追加成本 | LES计算规模 |
|---|---|---|
| Fluent | 基础包含 | 500万-2000万单元 |
| STAR-CCM+ | 反应流Token | 1000万-5000万单元 |
| CONVERGE | 基础含 | AMR实效500万-2000万单元 |
| OpenFOAM | 免费 | 网格生成手动 |
预混合燃烧因应用而最优工具差别很大呢。
正确。SI发动机选CONVERGE，燃气轮机选STAR-CCM+或OpenFOAM，这是现状的分工。
ANSYS vs STAR-CCM+ 预混合火焰比较——"哪个是对的"这个问题的难度
商用工具的预混合火焰模型对比中经常出现"同一燃烧器Fluent和STAR-CCM+的结果相差10%"的问题。原因大多是模型的默认参数（乱流Schmidt数、火焰拉伸模型系数等）在各厂商间有微妙差异。"哪个正确？"这个问题实际上没有意义，正确的问法是"哪个工具的验证数据集与我们的条件接近？"。现场的行家一致的说法是"先对同一实验数据检查哪个工具契合，然后再信任默认值"。工具的优劣不是关键，而是针对企业设计条件的验证积累才是要点。
预混合火焰的前沿研究
前沿课题和研究动向
请介绍预混合火焰的最前沿研究。
重点关注3个方向。(1) 氢气预混合火焰的LES，(2) 热声不稳定性的预测，(3) 基于DNS数据的模型改良。
氢气预混合火焰
氢气的预混合火焰有什么特殊问题？
氢气的Lewis数 $Le \approx 0.3$ 极小，因此产生强的差扩散效应，在火焰面局部产生温度升高（热扩散不稳定性）。火焰面形成细胞结构，通常的预混合燃烧模型假设破裂。
H2预混合火焰的LES需要采取以下措施。
- TFM效率函数中包含Lewis数效应
- FGM表格用$Le \neq 1$重新计算
- 可能的话用DNS进行验证
热声不稳定性
燃烧振动的CFD预测能做到什么程度？
用LES获取压力和发热变动的时序，计算Rayleigh指数。
$RI > 0$ 的区域驱动不稳定性。这种方法已能定量预测燃气轮机燃烧振动的频率和幅度。CERFACS（法）和TU Munich的研究组走在前沿。
DNS-informed 模型
DNS的成果被用于模型改良吗？
| DNS知见 | 对模型的影响 |
|---|---|
| 高Ka数火焰面断裂 | FSD方程消滞项修正 |
| 差扩散导致局部$S_L$变动 | TFM效率函数修正 |
| 火焰-乱流逆相关 | $\Xi$-equation 生成项修正 |
DNS被当作"数值实验"来改良模型呢。
氢气燃烧速度340m/s让发动机设计者烦恼的原因
汽油的层流燃烧速度约0.4m/s，而氢气约3m/s，在希薄条件下也能燃烧。某些研究还报告接近爆邻域时达到340m/s。这样异常高的速度导致氢气在预混合发动机中易出现"早期着火（预点火）"和"逆火（回火）"。火焰逆流入吸气口的现象。目前的氢发动机开发中为了防止这些，主要采用直喷+高EGR策略，但氢燃烧速度过快的处理仍然是尖端研究课题。
预混合火焰故障排除
故障排除
请教预混合火焰计算的故障。
1. 火焰向上游传播（回火）
症状：火焰从燃烧室向混合区逆流。
原因和对策：
- 乱流燃烧速度 $S_T$ 评估过大 → 降低Zimont常数$A$（0.52 → 试试0.4）
- 网格太粗，数值扩散导致火焰前线向上游扩散 → 细化包含火焰面的区域
- 边界层内低速区域火焰传播 → 检查壁面附近网格分辨率
2. 火焰不稳定（吹熄）
症状：火焰流出燃烧室并熄灭。
对策：
- 检查 $S_L$ 设置值是否过小（特别是高压条件下 $S_L$ 减小）
- 试点火焰或保炎器的再循环区是否充分分析
- 乱流模型是否正确预测再循环区（k-$\varepsilon$ 有时过度扩散，推荐SST）
3. 进行变量 $c$ 的数值问题
| 症状 | 原因 | 对策 |
|---|---|---|
| $c$ 超出 [0,1] 范围 | 数值方案的溅出 | 用有界方案（避免QUICK） |
| $c$ 前线过厚 | 网格分辨率不足 | 审视TFM的$F$值，细化网格 |
| $c$ 噪声大 | 压力-密度耦合不良 | 密度比补正，Under-Relaxation调整 |
4. 压力振动（热声）
计算中出现非物理的压力振动怎么办？
5. 当量比变动的情况
燃料浓度不均匀的情况怎么处理？
这不是完全预混合，而是部分预混合（Partially Premixed）条件。用Fluent的部分预混合燃烧（火焰片元+进行变量）或STAR-CCM+的FGM处理。纯粹预混合模型适用会导致局部过浓、希薄区域出现非物理结果。
预混合火焰的麻烦根源就是"$S_L$的准确性"和"火焰面的数值分辨率"呢。
正确。先用Cantera验证$S_L$，再检查网格分辨率。跳过这两步调试就会陷入泥沼，务必执行。
预混合火焰"振动"——热声不稳定的可怕现实
预混合火焰故障排除中容易被忽视的是"热声不稳定"。燃烧室内热释放率振动，与压力波共鸣会导致振动增幅。低NOx燃烧器采用希薄预混合时，稳定裕度降低，这种振动易发生。实际上1990年代燃气轮机新型燃烧器因热声振动导致燃烧室壁面损伤的事故有多起。CFD预测需要非定常LES，定常RANS看不到这个问题。"定常计算无问题"→"实机因燃烧器振动部件破损"这样的陷阱特别需要在设计评审时警戒。
相关课题
价值
详细
错误