层流火焰分析

官官方工学老师·26 页·深入（追求细节与边界）·0 次浏览·2 天前
✦ 层流火焰分析工学层流火焰分析简化分析
层流预混火焰传播燃烧分析物理建模

层流火焰分析

从火焰结构、传播速度到失稳机理，掌握完整分析体系

按 空格/→ 演示下一步

1 / 26 页

全部页面点击任意一页，跳回舞台从这页播放

层流预混火焰传播燃烧分析物理建模

层流火焰分析

从火焰结构、传播速度到失稳机理，掌握完整分析体系

1第 1 页 · 层流火焰分析

什么是层流预混火焰

上一节我们谈了层流火焰分析在燃烧学里为何是基础。这一页聚焦其中最经典的一类——层流预混火焰：燃料与氧化剂在燃烧前充分混合，并以层流方式传播。

预混
燃料与氧化剂在进入火焰区前已完成分子级均匀混合
层流
雷诺数低，流线平行有序，无宏观湍流脉动
火焰锋面
极薄的反应区，把已燃区与未燃区截然分隔开
层流燃烧速度 SL
锋面相对未燃气的法向传播速度，是混合物的本征参数
水面推进的薄波对应 →层流预混火焰

未燃混合气=平静水面，火焰锋面=薄波，薄波前进速度=SL

ρuSL=ρbub\rho_u S_L = \rho_b u_b
2第 2 页 · 什么是层流预混火焰

层流 vs 湍流火焰

层流和湍流是燃烧的两大模式，结构差异决定速度与建模方法完全不同，工程场景选错模型代价巨大。

层流火焰
  • 结构：光滑、平整、无限薄的反应层
  • 速度：唯一确定的层流燃烧速度 SL
  • 流动：Re 低，流线有序平行
  • 分析：一维预混火焰模型可直接求解
湍流火焰
  • 结构：皱折、增厚，局部可能熄火
  • 速度：有效湍流燃烧速度 ST 远大于 SL
  • 流动：Re 高，流场强烈随机脉动
  • 分析：LES、PDF、火焰面/小火焰模型
火焰面结构与基础速率选层流分析；工程尺度的燃烧与混合必须用湍流模型。
3第 3 页 · 层流 vs 湍流火焰

为什么要研究层流火焰

上一页我们看到湍流火焰比层流复杂得多，那为什么科研不直接啃硬骨头？因为层流不是湍流的反面，而是它的'基本语法'。

理论基石
湍流燃烧可视为层流小火焰的叠加，层流是可解析求解的简化极限
实验标定
用来验证化学反应机理和输运模型，是燃烧理论的'试金石'
核心参数
层流火焰速度、点火能量、熄灭极限等工程设计的关键数据
工程应用
内燃机、燃气轮机、火焰安全评估都以层流特性为起点
学外语先学字母对应 →研究燃烧先学层流

湍流是'复杂句子'，必须先掌握层流这套'语法'才能拆解

4第 4 页 · 为什么要研究层流火焰

简化分析的基本假设

四条假设，把层流火焰方程从'解不开'砍到'能求解'。

1
稳态
冻结时间——火焰结构不随时间变化
2
准一维
压扁空间——只留传播方向，其它方向均匀
3
忽略辐射
辐射换热不计入——只看传导与对流
4
快速化学
反应瞬时完成——燃烧区近似为无限薄面
5第 5 页 · 简化分析的基本假设

简化分析的控制方程

质量守恒、动量守恒（简化）、能量守恒、组分守恒方程的建立

简化分析的控制方程
质量守恒、动量守恒（简化）、能量守恒、组分守恒方程的建立
6第 6 页 · 简化分析的控制方程

层流火焰速度的定义

前面我们从物理图像出发推到了简化控制方程。本页要把方程里反复出现、又最容易被混淆的那个量单独拎出来讲清楚——层流火焰速度 S_L。

物理含义
未燃混合气进入火焰锋面的法向分速度；方向严格垂直于火焰面
符号与坐标系
记为 S_L；定义在随火焰一起运动的坐标系里，未燃气以 S_L 流入
物性参数
由混合气组分、温度、压力唯一决定，与火焰形状、流场无关
边界提醒
对平面绝热火焰≈实验室系推进速度；对弯曲或拉伸火焰，二者因拉伸率而偏离
雪铲推雪前进对应 →层流火焰速度 S_L

雪铲走过雪面，未压实的新雪以恒速从前方灌入铲刀；铲刀=火焰锋面

ρuSL=ρbSb\rho_u S_L = \rho_b S_b
7第 7 页 · 层流火焰速度的定义

层流火焰速度的推导

从简化假设出发，一步步推到层流火焰速度的解析式。

1
建立简化模型
一维稳态、绝热、Lewis=1、单步不可逆反应
2
写出控制方程
连续、能量、组分三个守恒方程
3
解耦合并
能量与组分通过放热参数耦合为一个方程
4
求解特征值
绝热火焰温度确定，火焰速度即为特征值
5
得到解析式
S_L 与热扩散率、反应速率的方根成正比
8第 8 页 · 层流火焰速度的推导

简化分析的温度分布

沿火焰传播方向看温度变化：从未燃气体 T₀，经预热区渐升，到反应区陡升至绝热火焰温度 T_ad。

图解渲染中…
a1：T₀ 为未燃气体的初始温度，≈环境温度a2：预热区宽度较大，仅靠热传导升温a3：反应区极薄，化学反应集中释热a4：T_ad 是热力学决定的火焰最高温度
9第 9 页 · 简化分析的温度分布

火焰厚度的定义与估算

上页看到火焰区温度从常温陡升到绝热火焰温度——这个'陡坡'有多陡？反应区到底有多宽？我们需要给它一个量化的长度尺度。

最大温度梯度定义
用温升除以反应区内温度梯度的最大值，得到等效厚度
0.1~1 毫米量级
甲烷空气约 0.5 mm；氢气更快，可降至 0.3 mm；稀释弱火焰可逼近 1 mm
与 S_L 的标度律
δ_T ~ α/S_L，热扩散越快、火焰速度越小，反应区越厚
反应与热扩散交叠
化学放热和热传导共同支撑的薄层，是把火焰作薄面处理的判据
头发丝对应 →火焰反应区厚度

头发直径约 0.07 mm，火焰厚度约 0.1~1 mm——刚好是肉眼刚能分辨、却又必须凑近才看清的尺度

δT  =  Tb−Tumax⁡∣dT/dx∣  ∼  αSL\delta_T \;=\; \frac{T_b - T_u}{\max |dT/dx|} \;\sim\; \frac{\alpha}{S_L}
10第 10 页 · 火焰厚度的定义与估算

简化与详细分析的差异

两种思路的前提假设根本不同，适用场景也截然分开。

简化分析
  • 假设反应无限快
  • 火焰面无限薄
  • 只解流动和能量方程
  • 温度由热力学唯一确定
详细分析
  • 有限化学反应速率
  • 火焰有有限厚度
  • 流动能量组分方程联立
  • 温度与组分耦合求解
工程估算用简化，机理研究和污染物预测用详细。
11第 11 页 · 简化与详细分析的差异

化学动力学的作用

上页看到，温度在很薄的一层内陡升，这正说明燃烧不是瞬间完成。反应速率有限，使放热和组分转化逐步发生，火焰便显出预热区与反应区。

有限反应速率
燃料消耗和放热在有限空间内完成，因此火焰具有可分辨的厚度。
预热区
反应区向未燃气体传热，主要发生升温并伴随有限的低温反应。
核心反应区
燃料与氧化剂快速消耗，温度和大部分放热集中发生。
连续过渡
两区之间没有物理界面，温度、组分和速率均连续变化。
时间尺度竞争
化学时间与流动、输运时间相当时形成分区；化学极快时才接近瞬时反应面。
闸机限流的长队对应 →有限速率的火焰结构

闸机处理能力对应反应速率；等待区像预热区，集中通过区像反应区。越慢队伍越长，火焰厚度还受输运影响。

Da≡τflow/τchem\mathrm{Da}\equiv\tau_{\mathrm{flow}}/\tau_{\mathrm{chem}}
12第 12 页 · 化学动力学的作用

预热区的温度变化

预热区里，未燃气体怎么被一步步加热到接近绝热火焰温度？

1
未燃气体逼近火焰
在火焰坐标系下，冷气流以 S_L 速度稳定流向反应区
2
热传导逆流加热
高温反应区向上游冷气体释放热通量，方向与流动相反
3
温度指数式上升
能量守恒导出 T 沿 x 指数增长，特征长度约 α/S_L
4
临近 T_ad 过渡
温度逼近绝热火焰温度，化学反应开始主导
13第 13 页 · 预热区的温度变化

反应区的温度变化

反应区是一层极薄的化学工厂，未燃气在此快速跨越到燃烧产物，温度陡升至平衡值。

1
进入反应区
未燃气体已被预热到接近自燃温度，仅差临门一脚
2
链式反应爆发
自由基浓度指数累积，反应速率爆炸式增长
3
集中释放热量
单位体积释热率远高于预热区，构成火焰最热层
4
跃升至绝热火焰温度
温度从近着火温度陡升至T_ad，达到热力学平衡
14第 14 页 · 反应区的温度变化

详细火焰结构图解

从左到右看层流火焰的空间结构：未燃混合气经预热升温、剧烈反应，最终变为已燃产物。

图解渲染中…
a1：常温T₀，燃料浓度最高，无化学反应发生a2：温度由T₀升至点火温度T_ign，组分几乎不变a3：火焰前锋所在，燃料被快速消耗，温度接近T_ada4：反应终止，温度为绝热火焰温度T_ad，全为产物
15第 15 页 · 详细火焰结构图解

温度分布：简化 vs 详细

简化分析给"阶跃式"温升，详细分析是"渐变式"——同一物理问题，两条曲线的差异背后是根本不同的假设。

简化分析温升
  • 曲线形状：阶跃突变，温度从 T_u 直跳到 T_ad
  • 预热区：无厚度，瞬间完成热传导
  • 反应区：无穷大速率，火焰厚度趋于零
  • 终温：精确等于绝热火焰温度
详细分析温升
  • 曲线形状：S 型渐变，光滑上升趋近 T_ad
  • 预热区：导热先于反应，温度逐步攀升
  • 反应区：有限速率放热，温升出现明显曲率
  • 终温：渐近趋近 T_ad，存在余温尾巴
简化模型抓"终点"便于估算，详细模型显"过程"才能解释火焰厚度与传播——机理研究必须用详细模型。
16第 16 页 · 温度分布：简化 vs 详细

速度分布与火焰传播

从左到右是气流穿过火焰锋面的空间推进；下方连续性方程把速度与温度分布耦合在一起。

图解渲染中…
a1：未燃气入口，速度即层流火焰速度 u0，温度 T0a3：反应区，集中放热，速度梯度 du/dx 最大a4：已燃气出口，热膨胀使 uf/u0 = Tf/T0e1：一维稳态连续性方程，是 v–T 耦合之根
17第 17 页 · 速度分布与火焰传播

层流火焰传播机制

前面把火焰拆成了几个区、给出了 S_L 的公式。问题来了：这一薄层凭什么能自己前进？靠什么机制把'火'传下去？

传播本质：自维持行波
在火焰坐标系下火焰前沿是稳态行波，无外力推动也能持续向前推进
热量反向传导至预热区
高温反应区向尚未燃烧的上游冷气导热，温差本身就是传播的引擎
预热触发化学加速
预热气体温度越过临界后，阿伦尼乌斯反应速率呈指数级上升，开始剧烈反应
传热=放热：稳态平衡
预热区吸热速率等于反应区放热速率，这个平衡把传播速度 S_L 锁定下来
燃着的导火索对应 →层流火焰前锋

燃烧段把热量传给下一段，使其达到着火温度后再燃烧——火头一节接一节往下传

ρuSL cpΔT  ∼  λ ΔT/δL⇒SL∼α/δL\rho_u S_L \, c_p \Delta T \;\sim\; \lambda \,\Delta T / \delta_L \quad\Rightarrow\quad S_L \sim \alpha / \delta_L
18第 18 页 · 层流火焰传播机制

详细分析的控制方程

四分支同构展开，每个守恒方程都是『瞬态+对流=扩散+源』的统一形式。

图解渲染中…
B：∂ρ/∂t+∇(ρv)=0，无反应源C：ρDv/Dt=-∇p+∇τ，含粘性D：ρcp DT/Dt=∇(k∇T)+ω̇E：ρDY_i/Dt=∇(ρD∇Y_i)+ω̇_i
19第 19 页 · 详细分析的控制方程

反应速率的严格处理

上一页我们看到简化分析把反应压成一张无限薄曲面，详细分析却画出一个有厚度的反应区。这中间的差距，就藏在 Arrhenius 形式的温度指数敏感性里。

Arrhenius 速率形式
ω̇ ~ A·exp(-E/RT)，对温度是指数级敏感
有限速率 → 有限厚度
化学时间 τ_chem 有限，反应区尺度 ~ α/S_L
Damköhler 数与极限
Da = τ_flow/τ_chem；Da→∞ 时趋近薄火焰极限
反应区厚度的量级
δ_r ~ α/S_L，化学时间与扩散时间的平衡
热水泡茶对应 →Arrhenius 速率与反应区厚度

水温越高溶出越快（指数敏感），但总要泡一段有限时间——不会瞬间无味

ω=Aexp⁡ ⁣(−ERT)[F]a[O]b,Da=τflowτchem\omega = A\exp\!\left(-\frac{E}{RT}\right)[F]^a[O]^b, \quad Da = \frac{\tau_{flow}}{\tau_{chem}}
20第 20 页 · 反应速率的严格处理

质量分数分布

燃料、氧化剂、产物浓度在火焰锋面内的变化规律

质量分数分布
燃料、氧化剂、产物浓度在火焰锋面内的变化规律
21第 21 页 · 质量分数分布

层流火焰速度的影响因素

逐一改变条件，看火焰速度如何响应——这是理解影响因素的标准方法。

1
建立基准
1 atm、298 K、φ=1的甲烷-空气，SL约0.4 m/s，作为对比零点
2
改变温度
SL ∝ T^n，n≈2，每升100 K近翻倍，Arrhenius指数放大效应
3
改变压力
SL ∝ p^m，m为负（甲烷约−0.5），高压反而减慢
4
更换燃料
H2可达3 m/s，甲烷仅0.4 m/s，分子结构决定基础速度量级
5
调当量比
φ≈1.05–1.1附近出现峰值，过贫过富都下降，呈非单调曲线
6
综合关联
实验拟合SL ∝ T^α p^β f(φ)，多参数耦合，无统一解析式
22第 22 页 · 层流火焰速度的影响因素

Markstein长度定义

前面我们看到火焰有'厚度'——预热区让气体升温、反应区集中放热。但真实火焰前缘往往不是平面：它会弯曲、会受应变。这种拉伸的'影响能传到多深'，要用Markstein长度来表征。

定义
火焰前缘对曲率和应变响应的特征长度
物理意义
拉伸扰动在火焰结构内'传播'的距离
符号与稳定性
正→抑制扰动；负→放大扰动（稀氢火焰）
微观根源
反映热扩散与质量扩散（Lewis数）的相对强弱
涟漪扩散距离对应 →Markstein长度

石子入水，涟漪传多远；Markstein长度就是拉伸信号在火焰中的'扩散半径'

Sb=SL−LM⋅κS_b = S_L - L_M \cdot \kappa
23第 23 页 · Markstein长度定义

火焰拉伸与稳定性

从拉伸出现到熄火发生，火焰经历五步临界演化，每一步都在逼近稳定性的边界。

1
拉伸的来源
火焰面受曲率和流体应变率共同拉伸，单位时间内面积变化
2
热质传递变化
拉伸改变未燃气体进入反应区的速率，破坏局部热平衡
3
火焰响应判断
用Markstein长度判别火焰速度随拉伸是增大还是减小
4
临界拉伸率
存在临界拉伸率，超过后火焰传播速度无法补偿热损失
5
熄火发生
反应区释热小于散热，火焰失去自维持能力，熄灭
24第 24 页 · 火焰拉伸与稳定性

层流火焰分析要点

  • ✓简化假设一步无限快反应；详细保留完整有限速率化学
  • ✓火焰速度由释热与扩散竞争决定，是多参数的函数
  • ✓预热区与反应区尺度共同决定火焰厚度量级
  • ✓Markstein 长度刻画拉伸响应，决定火焰稳定性边界
延伸主题：湍流火焰模型与燃烧模化点火、熄火与可燃极限火焰-声学/流场耦合
25第 25 页 · 层流火焰分析要点

课后思考

先自己想，再看参考答案。三问分别落在核心机制、工程意义和模型边界三个层次。

1为什么反应越快、热扩散越慢，层流火焰速度反而越大？背后是什么竞争机制？

参考答案火焰传播是「反应放热→热扩散预热未燃气体→未燃气体被点燃」的接力。反应快则放热快，热扩散慢则热量集中在薄层内预热效率高，两者都加速这场接力。

2既然发动机里几乎都是湍流燃烧，为什么我们仍要花大量精力研究层流火焰？

参考答案层流火焰是湍流燃烧模型的基元参数——火焰面模型、输运方程方法等都需要层流火焰速度和厚度作为输入。没有可靠的层流数据，湍流模型无从校准。

3简化分析忽略了辐射热损失和详细输运效应，加入这两项后温度分布会如何偏离？

参考答案辐射从高温区额外带走热量，压低峰值温度并使分布变平缓；多组分扩散与热扩散的差异让各组分扩散速率不同，反应区位置和火焰速度都会偏离简化预测。

26第 26 页 · 课后思考
层流火焰分析 · 知识图解