稳态火焰小尺度模型
稳态火焰小尺度的理论基础
概述
老师，火焰小尺度模型是什么？
火焰小尺度模型是Peters在1984年提出的非预混合湍流燃烧模型。将湍流扩散火焰视为"众多薄层流火焰（火焰小尺度）的集合体"，将火焰的内部结构归结为混合分数 $Z$ 的一维问题。
可以将3D湍流燃烧简化为1D吗？
是的。火焰小尺度的假设是"火焰厚度小于湍流的最小尺度（Kolmogorov尺度）"。在这种情况下，火焰的内部结构局部上等同于1D层流对向流火焰。
火焰小尺度方程
请告诉我火焰小尺度方程。
化学种 $Y_i$ 的火焰小尺度方程可写为：
其中 $\chi$ 是标量耗散率，表示 $Z$ 空间中的扩散。
稳态火焰小尺度中的$\partial/\partial t = 0$。
是的。在稳态火焰小尺度中，时间导数项为零，耗散率 $\chi_{st}$（在化学计量面处的值）作为参数的常微分方程。当 $\chi_{st}$ 增大时，反应跟不上，最终导致消焰（quenching dissipation rate $\chi_q$）。
S型曲线（S-curve）
什么是S型曲线？
稳态火焰小尺度解的最高温度对 $\chi_{st}$ 作图会形成S形曲线。上支是燃烧状态，下支是未燃状态，中间支是不稳定解。
$\chi_q$ 大约是多少？
在甲烷/空气情况下 $\chi_q \approx 20-50$ s$^{-1}$，氢气/空气时 $\chi_q \approx 1000$ s$^{-1}$ 非常大。这意味着氢气火焰不容易消焰。
火焰小尺度库
将稳态火焰小尺度解以 $\chi_{st}$ 为参数计算多个，并存储为 $(Z, \chi_{st})$ 的二维表。如果提前应用乱流的 $\beta$-PDF平均，就成为 $(\widetilde{Z}, \widetilde{Z''^2}, \widetilde{\chi_{st}})$ 的三维查询表。
火焰小尺度模型的核心是"1D层流火焰的表 + 乱流PDF"的组合吧。
完全正确。包含详细化学反应的火焰结构是离线计算的，3D CFD中只需查表，所以计算成本非常低。
把火焰变成"数据库"的人——Norbert Peters的火焰小尺度思想
火焰小尺度模型的发起人Norbert Peters在1984年提出的想法简单说就是"乱流中的火焰最后是小层流火焰的集合"。那么就可以提前计算层流火焰并存储在表中，本计算中只需参考一个混合分数即可。这种"提前计算→查表"的思想，戏剧性地降低了燃烧CFD的计算成本。Peters在亚琛工业大学发展了这个概念，后来衍生出FGM（Flamelet-Generated Manifolds）和FPV（Flamelet Progress Variable）模型。如今涡轮机设计离不开他的遗产。
稳态火焰小尺度的数值计算方法
数值方法详细说明
火焰小尺度库怎样制作？
库的构建有两个步骤。(1) 对向流扩散火焰的计算，(2) PDF积分进行表的生成。
对向流扩散火焰的计算
具体怎样计算？
用Cantera、FlameMaster、OPPDIF（CHEMKIN-PRO）等求解1D对向流扩散火焰。逐步增加应变率（应变速率）$a$，追踪到消焰点。$a$ 与 $\chi_{st}$ 的关系为 $\chi_{st} \approx a \cdot f(Z_{st})$。
典型的库构建参数：
| 参数 | 推荐值 | 备注 |
|---|---|---|
| $Z$ 方向格点 | 128-256 | 在 $Z_{st}$ 附近集中 |
| $\chi_{st}$ 的分割 | 30-50点 | 对数等距 |
| 应变速率范围 | 1 - $a_q$ s$^{-1}$ | 到消焰 |
| 反应机理 | GRI-Mech 3.0等 | 采用详细机理 |
Fluent中的实现
Fluent里怎样设置？
在Fluent的Non-Premixed Combustion模型内进行以下设置。
1. 选择Flamelet Model（与Equilibrium Chemistry切换）
2. 以CHEMKIN格式导入反应机理
3. 设置Number of Flamelet数量（默认20，推荐30-50）
4. 设置PDF表的分辨率
5. 计算开始后，$\widetilde{Z}$、$\widetilde{Z''^2}$ 的传输方程被求解，通过查表确定温度·化学种
FGM（Flamelet Generated Manifold）的关系
FGM与火焰小尺度模型有什么区别？
FGM是van Oijen等人（2000年）提出的方法，从火焰小尺度解构造低维多样体（manifold）。在稳态火焰小尺度基础上引入进度变量 $C$（Progress Variable），能表达着火·消焰的过渡过程。
稳态火焰小尺度与FGM的比较：
| 特性 | 稳态火焰小尺度 | FGM |
|---|---|---|
| 表维度 | 2-3D ($Z$, $Z''$, $\chi$) | 3-4D ($Z$, $Z''$, $C$, $C''$) |
| 消焰再现 | 通过S曲线可能 | 通过Progress Variable自然再现 |
| 自着火 | 困难 | 可对应 |
| 部分预混合 | 困难 | 可对应 |
| 计算成本 | 非常低 | 低（查表） |
FGM的通用性更好啊。
是的。在STAR-CCM+和OpenFOAM中FGM已成为主流。Fluent在R2以后也强化了FGM选项。不过稳态火焰小尺度作为最简单、最稳定的方法，仍被广泛应用。
火焰小尺度模型的关键是"制作优质的表"吧。
完全正确。表的分辨率和反应机理的合理性直接决定模型精度。
PDF的"形状"选择决定结果——β-PDF与Clipped Gaussian的比较
在火焰小尺度/PDF模型的数值实现中避免不了"混合分数PDF形状的选择"问题。最常见的是β-PDF（贝塔分布），ANSYS Fluent默认采用。但β-PDF虽然数学上易处理，却只能表示单峰形态，在完全燃料·完全氧化剂边界处行为过于陡峭。Clipped Gaussian（截断高斯分布）在许多情况下与实验数据拟合较好，但表生成耗时。使用哪一种会导致最高温度相差50～100K，这是实务技能差异的所在。
稳态火焰小尺度的实务应用
实践指南
请教我火焰小尺度模型的实务使用方法。
在燃气轮机燃烧室、锅炉的设计中，火焰小尺度/FGM模型是最有效的选择。
适用范围的判断
请教我火焰小尺度模型能应用的条件。
根据燃烧图（Borghi/Peters diagram）上的位置进行判断。
| 机制 | $Da$ | $Ka$ | 火焰小尺度适用 |
|---|---|---|---|
| 层流火焰 | -- | $Ka < 1$ | 可应用（可能不必要） |
| Flamelet制度 | $Da > 1$ | $1 < Ka < 100$ | 最优 |
| Thin reaction zone | $Da > 1$ | $Ka > 100$ | 有限 |
| Well-stirred reactor | $Da < 1$ | -- | 不适（EDC更合适） |
用Damkohler数和Karlovitz数进行判断啊。
是的。大多数燃气轮机燃烧室和锅炉都属于Flamelet制度，所以火焰小尺度模型可以应用。不过在稀释燃烧条件下Da接近1时，火焰小尺度假设就变得可疑了。
网格要求
火焰小尺度模型采用查表，所以化学反应的网格要求较松，但混合场的分辨率很重要。
| 要求 | RANS | LES |
|---|---|---|
| 总单元数（燃气轮机燃烧室） | 200万-500万 | 2000万-5000万 |
| 喷射孔附近 | 喷射孔直径/10 | 喷射孔直径/20 |
| 火焰区域 | 1-2 mm | 0.5-1 mm |
| y+ | 30-100 | < 1 (壁面分辨) |
常见失败与对策
| 现象 | 原因 | 对策 |
|---|---|---|
| 温度与Equilibrium相同 | $\chi_{st}$的效果弱 | 确认$\chi$的传输方程，调整$C_\chi$常数 |
| 消焰未被再现 | 稳态火焰小尺度的极限 | 改用FGM + Progress Variable |
| 表外的$Z$值 | 燃料入口Z设置错误 | 确认边界条件Z值（0和1） |
| CO排放过高 | 表分辨率不足 | 增加$Z_{st}$附近表分割 |
火焰小尺度模型的优点是计算轻快和详细化学反应兼顾吧。
完全正确。即使采用GRI-Mech 3.0的53化学种，3D CFD的运行时只是查表，成本不到EDC的1/100，可获得详细的化学种分布。
火焰小尺度表的"大小"问题——吃内存太多的失败事例
火焰小尺度模型实践中初学者容易遇坑的是查询表的大小设置。混合分数Z、标量耗散率χ、进度度变量C三维表中，各轴各设100点就变成100×100×100=100万条目。每张表几百MB。在大规模并行计算中，各计算节点都会复制表，在单节点2GB RAM上限下，仅表就会撑爆内存——真实发生过这样的悲剧。现场规则是"表最初从50点开始，确认网格收敛后→最终80～100点"。由于表粗糙引起的数值误差常被忽视，要留意。
稳态火焰小尺度的软件比较
商业工具比较
请告诉我火焰小尺度/FGM模型的工具支持情况。
| 工具 | 稳态火焰小尺度 | FGM | 表生成 | LES支持 |
|---|---|---|---|---|
| Ansys Fluent | 标准配置 | R2以后强化 | 自动（CHEMKIN） | 有 |
| STAR-CCM+ | 配置 | 标准推荐 | DARS联动 | 有 |
| OpenFOAM | 社区版 | flameletFoam等 | 手动（Cantera） | 有 |
| CONVERGE | 配置 | 配置 | SAGE联动 | 有 |
| FlameMaster | 表生成专用 | -- | 1D火焰小尺度计算 | -- |
Ansys Fluent
Fluent火焰小尺度模型的特点？
Fluent在Non-Premixed Combustion模型内，可在Equilibrium / Steady Flamelet / Unsteady Flamelet间切换。表生成完全自动化是最大优点。不过相比STAR-CCM+，FGM的实现稍显后发。
STAR-CCM+
STAR-CCM+ FGM的特点？
STAR-CCM+的Flamelets/FGM模型可灵活定制Progress Variable的定义。通过DARS库联动，从1D火焰小尺度计算到表生成无缝衔接。向Part Premixed模型的拓展也容易，在燃气轮机燃烧室LES中有多项实绩。
FlameMaster
什么是FlameMaster？
这是Peters课题组（RWTH亚琛）开发的1D火焰小尺度计算专用代码。许多研究者把它与OpenFOAM组合使用。对向流扩散火焰计算精度高，S曲线追踪稳定。可免费利用。
表生成工作流
OpenFOAM + FlameMaster/Cantera 的典型工作流：
1. 用Cantera/FlameMaster计算1D对向流扩散火焰（应变率扫描）
2. 将火焰小尺度解转换到Z空间
3. 通过$\beta$-PDF积分生成$(\widetilde{Z}, \widetilde{Z''^2})$表
4. 转换为OpenFOAM格式的查询表
5. 用OpenFOAM的flameletFoam等读取进行计算
火焰小尺度模型工具选择多样，能根据用途灵活组合啊。
是的。商业工具的自动表生成很方便，但OpenFOAM+Cantera的组合能完全控制表的构建过程，有优势。
火焰小尺度模型工具——Ansys Fluent Flamelet vs StarCCM+
火焰小尺度模型是扩散燃烧代表性CFD方法，Ansys FluentFlamelet Generated Manifold（FGM）与Siemens StarCCM+的Composed FGM是商业工具主流。FGM通过1D火焰计算（Cantera等）生成查询表，导入3D CFD，模拟详细化学反应，同时高速化计算。StarCCM+用户自定义反应模型定制性强，特殊燃料（氨·氢·SAF）燃烧模型适应快。氢混焼喷枪设计中，2023年后两工具的GRI-Mech 3.0支持和氢火焰lift解析成为竞争焦点。
稳态火焰小尺度的前沿研究
前沿话题和研究动向
请教我火焰小尺度模型的最新动向。
火焰小尺度模型是成熟技术，但其应用范围拓展和精度提高的研究在继续。
非稳态火焰小尺度模型
什么是非稳态火焰小尺度？
稳态火焰小尺度令$\partial/\partial t = 0$，非稳态火焰小尺度模型保留这一项，对时间变化的$\chi_{st}(t)$进行响应追踪。能再现消焰后的再着火以及自着火过程。Fluent中可作为"Unsteady Flamelet"选择。
多维FGM
FGM还在进一步发展吗？
最近FGM研究方向是增加表的维度以提高精度。
| 表变量 | 维数 | 能再现的现象 |
|---|---|---|
| $Z$, $C$ | 2D | 基本消焰·着火 |
| $Z$, $C$, $h$ | 3D | 热损失效果 |
| $Z$, $C$, $h$, $Y_{NO}$ | 4D | NOx预报改善 |
| $Z$, $Z''$, $C$, $C''$ | 4D | 乱流完全PDF效果 |
4D表内存没问题吗？
各维各100点时4D达到 $10^8$ 条目。压缩表和稀疏网格的实用化在推进。用机器学习代替表也有研究。
LES-Flamelet
LES与火焰小尺度的组合？
LES-Flamelet是乱流燃烧LES中计算效率最高的手法。用$\beta$-PDF模型化亚网格尺度的$Z$方差。GE航空和赛峰等燃气轮机厂商正式导入设计LES-Flamelet/FGM。
氢气和氨气燃烧的扩展
新燃料的应用如何？
氢火焰因为Lewis数效应需要修正火焰小尺度假设（差扩散火焰小尺度）。氨气燃烧的化学时间尺度长，Da值小，火焰小尺度假设难以成立。NH3燃烧用EDC或PaSR更合适的情况多。
火焰小尺度模型并非万能，需根据燃烧制度灵活选择啊。
完全正确。根据Peters图判断是否属于Flamelet制度是使用火焰小尺度模型前必须的基本检查。
FGM与FPV——火焰小尺度拓展的两大流派目标
稳态火焰小尺度（SFM）的先端形式FGM（Flamelet-Generated Manifolds）与FPV（Flamelet Progress Variable）都用"混合分数+进度度变量"的2D表表达非平衡燃烧，但表的制作方法不同。FGM由亚琛工大van Oijen在2000年提出，能向预混合火焰拓展，SFM难以处理。FPV由斯坦福Pierce和Moins在2004年提出，在再着火和消焰模型化上强。两者都基于"火焰用2个变量记述"假设，在高压·超稀薄条件下该假设破裂，这就是至今仍有研究继续的理由。
稳态火焰小尺度的故障排除
故障排除
请教我火焰小尺度模型特有的故障。
1. 火焰小尺度表生成失败
症状：1D火焰小尺度计算在到达消焰点前发散。
对策：
- 用0D计算（着火延迟、层流燃烧速度）确认反应机理的妥当性
- 细化应变速率步长（特别是消焰附近）
- 用Cantera/FlameMaster单独生成表，导入CFD
2. 平衡与差异不出现
稳态火焰小尺度却与Equilibrium结果相同…
原因：乱流场的$\chi_{st}$仅到达S曲线上支（不发生消焰条件）。此时稳态火焰小尺度实际上与Equilibrium等价。
确认方法：后处理对比$\widetilde{\chi_{st}}$最大值与$\chi_q$。如果$\widetilde{\chi_{st}} << \chi_q$则消焰效果可忽视。要看消焰效果的话，在应变强的条件（高速喷流等）下验证。
3. 进度变量（FGM）的问题
| 症状 | 原因 | 对策 |
|---|---|---|
| 进度变量$C$不收敛 | $C$的定义不当 | 改用$C = Y_{CO_2} + Y_{H_2O}$等稳定定义 |
| 表查询出现clipping | $C$超出表范围 | 对$C$设min/max限制 |
| 着火延迟与实验不符 | 表的$C$方向分辨率不足 | 增加表点数 |
4. 部分预混合条件
不完全是非预混合的条件（部分预混合）怎样？
如有试点火焰附带的主预混合气等部分预混合条件，纯火焰小尺度模型就不够。要改用Fluent的Partially Premixed Combustion模型（火焰小尺度+进度变量）或STAR-CCM+的FGM模型。
一般的调试步骤
1. 独立确认1D火焰小尺度计算（S曲线、消焰点）
2. 对Z的函数将表温度·化学种分布可视化
3. 确认CFD的$\widetilde{Z}$场与$\widetilde{Z''^2}$场（边界条件整合性）
4. 确认$\widetilde{\chi_{st}}$场（是否到达消焰条件）
火焰小尺度模型的调试是表和CFD的切分要点啊。
完全正确。表品质差了再好的CFD也改善不了结果。先在1D世界确认物理，再进3D才是正道。
火焰小尺度表的"洞"——补间误差引发的谜之发散
火焰小尺度模型故障排除中难的是"表补间误差"。混合分数Z和标量耗散率χ网格粗的区域，CFD计算点落在表网格点间隙，线性补间精度急剧恶化。特别是消焰附近（χ大的区域）的表有极端变化，网格疏则补间误差导致温度局部急变，引发发散。经验上"看表生成日志，重新细化χ方向网格（对数轴等距）就解决了"的例子多。发散原因归罪于求解器前，要先怀疑表。这是燃烧诊断的铁则。
相关话题
价值
更详细
错误