定常框架模型
理论与物理
概述
老师，什么是小火焰模型？
小火焰模型是Peters于1984年提出的非预混湍流燃烧模型。它将湍流扩散火焰视为“众多薄层流火焰（小火焰）的集合体”，并将火焰内部结构归结为混合分数 $Z$ 的一维问题。
也就是说，可以将3D的湍流燃烧简化成1D问题吗？
是的。小火焰的假设是“火焰厚度小于湍流的最小尺度（Kolmogorov尺度）”。此时火焰的内部结构在局部上与一维层流对向流火焰等同。
小火焰方程
请讲解一下小火焰方程。
化学组分 $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）模型等系列谱系。现在的燃气轮机设计离不开他的遗产。
各项的物理意义
- 时间项 $\partial(\rho\phi)/\partial t$：想象一下拧开水龙头的瞬间。最初水流会不稳定地喷溅，过一会儿才会变成稳定的水流，对吧？描述这个“变化过程中”的就是时间项。心脏跳动时血流脉动，发动机阀门每次开闭时流动变化，这些都是非稳态现象。那么稳态分析是什么？就是只看“经过足够时间流动稳定之后”——也就是将此项设为零。计算成本会大幅下降，所以先用稳态求解是CFD的基本策略。
- 对流项 $\nabla \cdot (\rho \mathbf{u} \phi)$：把落叶丢进河里会怎样？会被水流带着运往下游，对吧？这就是“对流”——流体运动搬运物质的效果。暖风的暖气能送到房间的另一头，也是因为空气这个“搬运工”通过对流输送了热量。这里有趣的是——这项包含“速度×速度”，因此是非线性的。也就是说，流速变快这项会急剧增强，变得难以控制。这就是湍流的根本原因。常见的误解：“对流和传导差不多”→ 完全不一样！对流是流动搬运，传导是分子传递。效率有天壤之别。
- 扩散项 $\nabla \cdot (\Gamma \nabla \phi)$：有过在咖啡里倒入牛奶后放置不管的经历吗？即使不搅拌，过一会儿也会自然混合。那就是分子扩散。那么下一个问题——蜂蜜和水，哪个更容易流动？当然是水对吧。因为蜂蜜的粘性（$\mu$）高，不容易流动。粘性越大扩散项越强，流体的运动就变得“粘稠”。雷诺数小的流动（缓慢、粘稠）中扩散占主导。相反，Re数大的流动中对流占压倒性优势，扩散则成为配角。
- 压力项 $-\nabla p$：按下注射器的活塞，液体就会从针头有力地射出，对吧？为什么呢？因为活塞侧是高压，针头侧是低压——这个压力差产生了推动流体的力。大坝放水也是同样的原理。天气图上等压线密集的地方会怎样？没错，会刮强风。“有压力差的地方就会产生流动”——这就是纳维-斯托克斯方程压力项的物理意义。这里的误解点：CFD的“压力”通常指表压而非绝对压力。切换到可压缩分析时结果突然变得奇怪，原因可能就是混淆了绝对压力/表压。
- 源项 $S_\phi$：被加热的空气会上升——为什么呢？因为比周围空气轻（密度低），被浮力推上去了。这个浮力作为源项添加到方程中。此外，燃气灶的火焰产生化学反应热，工厂的电磁泵对金属熔液施加洛伦兹力……这些都是“从外部向流体注入能量或力”的作用，都用源项表示。忘记源项会怎样？自然对流分析中忘记加入浮力，流体就完全不动——冬天在房间里开了暖气，暖空气却不上浮，这种物理上不可能的结果就会出现。
假设条件与适用范围
- 连续介质假设：克努森数 Kn < 0.01（分子平均自由程 ≪ 特征长度）时成立
- 牛顿流体假设：剪切应力与应变速率呈线性关系（非牛顿流体需要粘度模型）
- 不可压缩假设（Ma < 0.3时）：将密度视为常数。马赫数0.3以上需考虑可压缩性效应
- Boussinesq近似（自然对流）：仅在浮力项中考虑密度变化，其他项使用恒定密度
- 不适用的情形：稀薄气体（Kn > 0.1）、超音速/高超音速流动（需要捕捉激波）、自由表面流动（需要VOF/Level Set等）
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 速度 $u$ | m/s | 入口条件中从体积流量换算时，注意截面积单位 |
| 压力 $p$ | Pa | 区分表压和绝对压力。可压缩分析中使用绝对压力 |
| 密度 $\rho$ | kg/m³ | 空气: 约1.225 kg/m³@20°C、水: 约998 kg/m³@20°C |
| 粘性系数 $\mu$ | Pa·s | 注意与运动粘性系数 $\nu = \mu/\rho$ [m²/s] 混淆 |
| 雷诺数 $Re$ | 无量纲 | $Re = \rho u L / \mu$。层流/湍流转换的判定指标 |
| CFL数 | 无量纲 | $CFL = u \Delta t / \Delta x$。直接关系到时间步长的稳定性 |
数值解法与实现
数值方法的细节
小火焰库是怎么制作的？
库的构建有两个步骤。(1) 对向流扩散火焰的计算，(2) 通过PDF积分创建表格。
对向流扩散火焰的计算
具体是怎么计算的呢？
使用Cantera、FlameMaster、OPPDIF（CHEMKIN-PRO）等求解一维对向流扩散火焰。逐步提高应变率（应变速率）$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 Flamelets（默认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，“理所当然地尝试改变PDF形状”的经验会成为实务技能的差距。
迎风格式（Upwind）
一阶迎风：数值扩散大但稳定。二阶迎风：精度提高但有振荡风险。高雷诺数流动中必须使用。
中心差分（Central Differencing）
二阶精度，但Pe数 > 2时会产生数值振荡。适用于低雷诺数的扩散主导流动。
TVD格式（MUSCL、QUICK等）
通过限制器函数抑制数值振荡同时保持高精度。对捕捉激波或陡峭梯度有效。
有限体积法 vs 有限元法
なった
詳しく
報告