简化理论描述

官方工学老师·19 页·深入(追求细节与边界)·0 次浏览·2 天前
层流燃烧非预混火焰简化建模边界分析

层流非预混火焰的简化理论

能判断简化模型的适用边界,看清每一步近似背后的代价

按 空格/→ 演示下一步

1 / 19 页

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

层流燃烧非预混火焰简化建模边界分析

层流非预混火焰的简化理论

能判断简化模型的适用边界,看清每一步近似背后的代价

1第 1 页 · 层流非预混火焰的简化理论

层流非预混火焰是什么

上一节我们把火焰当成一张薄薄的反应面。但真实燃烧里,燃料和氧气不是一开始就混好的——这叫非预混火焰。

不预先混合
燃料与氧气各走各的输运通道,到相遇时才发生反应
扩散混合区
极薄的反应面,由燃料与氧化剂扩散相遇形成
层流
流场有序平滑,无湍流脉动,可写出解析解
典型场景
蜡烛、油灯、未预混的射流扩散火焰
两人从两侧走近碰面对应 →层流非预混火焰

燃料从一侧走来,氧气从另一侧走来,相遇瞬间才燃烧

2第 2 页 · 层流非预混火焰是什么

为什么要用简化理论

上一页我们定义了层流非预混火焰——燃料与氧化剂分别输运、在界面混合燃烧的薄层。听起来结构清晰,但真要在计算机里把这个薄层「算」出来,会撞上一连串硬骨头。

多场强耦合
流体力学、传热传质、化学反应同时求解,每个子问题已是棘手的非线性方程
化学维度灾难
甲烷机理已含50+组分、300+基元反应,大分子燃料可达上千组分,刚性极强
多尺度跨度
从分子碰撞纳秒级到火焰结构毫米级,时间与空间跨度近9个量级
抓主要矛盾
简化不是粗糙,而是剥离次要细节、保留燃烧主导机制的策略性取舍
工程落地桥梁
以可接受精度换取计算可行性,让工业燃烧器设计与优化真正可行
地图测绘对应 →简化燃烧理论

真实地理事无巨细但无法使用,地图只保留路网与关键地标;简化理论只保留主导物理,却能算出火焰整体行为

3第 3 页 · 为什么要用简化理论

层流火焰面假设

从左到右看:燃料与氧化剂分别从两侧向中心扩散,在火焰面汇合反应。

图解渲染中…
S无限薄的反应界面,燃料与氧化剂在此相遇并瞬时反应完毕Z假设火焰面厚度为零,把三维反应区简化成一张几何面P反应生成的产物(如CO₂、H₂O)从火焰面两侧向外扩散Q反应释放的热量通过传导传给未燃气体
4第 4 页 · 层流火焰面假设

无限快速反应假设

第二个基本假设及其意义

无限快速反应假设
第二个基本假设及其意义
5第 5 页 · 无限快速反应假设

混合分数的定义

上一页说无限快反应下,化学来不及思考,瞬间把混合物烧光。那「烧掉的」到底是什么比例?我们需要一把尺子,描述每个微元里燃料和氧化剂各占多少——这就是混合分数。

混合分数 ξ
无量纲标量场,描述当地微元来自燃料流的质量比例
取值范围 [0, 1]
ξ=0 处全为氧化剂,ξ=1 处全为燃料,中间是两者掺混
守恒标量
ξ 在非反应输运中只被对流与扩散搬运,无源项
基于元素定义
用某元素(如 C)在燃料流与氧化剂流的质量分数差归一化
解耦关键
无限快反应下,温度与组分都退化为 ξ 的单值函数
调鸡尾酒兑比例对应 →混合分数 ξ 跟踪燃料/氧化剂比例

ξ=1 是纯燃料瓶,ξ=0 是纯氧化剂瓶,ξ=0.3 就是三成燃料兑七成氧化剂

ξ=ZiZi,oxZi,fuelZi,ox\xi = \frac{Z_i - Z_{i,\mathrm{ox}}}{Z_{i,\mathrm{fuel}} - Z_{i,\mathrm{ox}}}
6第 6 页 · 混合分数的定义

混合分数的物理意义

两股燃料流和氧化剂流在炉内相遇后,每个流体微元都带着“来自哪一股”的质量来源信息。Z 的物理意义,就是用 0 到 1 记录这份来源比例;它衡量混合偏向,不衡量反应完成度。

来源份额
采用燃料流为 1 的常用约定,Z=0 对应纯氧化剂流,Z=1 对应纯燃料流;中间值表示两边来源的质量份额。
守恒的混合标记
在等扩散、Lewis 数为 1 的近似下,Z 只由对流和扩散输运,化学反应不改变其值。
局部偏向
Z 越接近 1,局部组成越偏燃料流;越接近 0,越偏氧化剂流。中间值只表示两者共存,不等于完全均匀。
有限信息
Z 不是燃料质量分数本身;给定流股组成后可由它换算燃料质量分数,但不能单凭 Z 给出温度或反应进度。
两种原液配饮料对应 →混合分数

燃料流对应其中一瓶原液;Z 记录它在这份饮料中的质量份额,不记录饮料是否已放热或反应。

7第 7 页 · 混合分数的物理意义

守恒标量方程的推导

从组分方程出发,加权求和并做扩散近似,推导出混合分数的守恒输运方程。

1
组分方程出发
列出每个物种的标准质量守恒方程:含对流、扩散、化学源项
2
加权求和
按燃料元素原子质量分数组合各物种方程,得到新标量 Z
3
源项相消
化学反应不增减元素,所有化学源项对加权和自动为零
4
等扩散近似
假设所有物种质量扩散系数相等(Le=1),让扩散项并成一项
5
得到守恒方程
得到与反应无关的标准标量输运方程,只取决于混合过程
8第 8 页 · 守恒标量方程的推导

混合分数的概率密度函数

守恒标量方程告诉你某瞬时某位置上 ξ 是几。但真实湍流里 ξ 一直在抖、每个瞬间都不一样——它变成了随机变量。本节引入 PDF,用统计语言刻画这种脉动。

湍流中 ξ 是随机变量
瞬时场每个空间点的 ξ 都不同,必须用统计方法描述整体行为
PDF 的定义 P(ξ)
曲线下面积恒为 1;ξ 落在某区间的概率=该段曲线下的面积
统计矩来自 P(ξ)
均值 ~ξ、方差等都由 P(ξ) 积分得到,刻画整体脉动水平
Favre 平均写法
密度加权 ~ξ=ρξ/ρ̄,使守恒方程避开密度脉动项,形式更简洁
燃烧中的应用
反应速率、标量场均对 P(ξ) 加权积分,免去逐点求解详细机理
全班同学的身高分布对应 →混合分数的 PDF

身高是随机变量,分布曲线告诉你 170cm 附近人多、两端人少;P(ξ) 同样告诉你哪个 ξ 最常出现

P(ξ)0,01P(ξ)dξ=1ξ~=01ξP(ξ)dξP(\xi) \geq 0,\quad \int_0^1 P(\xi)\, d\xi = 1 \\ \widetilde{\xi} = \int_0^1 \xi\, P(\xi)\, d\xi
9第 9 页 · 混合分数的概率密度函数

状态关系式的概念

上一节我们拿到了 Z 的概率密度函数。但 PDF 只是统计描述,真正预测温度场和浓度场,还得把 Z「翻译」成温度和组分——这层翻译就是状态关系式。

状态量只是 Z 的函数
温度 T 与各组分 Y_k 都可写成 T(Z)、Y_k(Z),空间坐标与时间消失
三个假设撑起唯一性
无限快速反应 + 绝热壁面 + 等扩散率(Le=1),三者缺一不可
反应动力学被绕过
关系式中不出现反应速率 ω,化学细节完全抹除
火焰塌缩为 Z_st 面
Z>Z_st 富燃、Z<Z_st 贫燃,火焰位于 Z=Z_st 等值面
配方比例与瞬时烹饪对应 →状态关系式

食材比例(Z)固定、烹调无限快→成品味道(T,Yi)只由比例决定,不看过程(反应动力学)

T=T(Z),Yk=Yk(Z)T = T(Z), \quad Y_k = Y_k(Z)
10第 10 页 · 状态关系式的概念

完全燃烧的温度与组分

从左上读起:先由 Z 算出组分与当量比 φ;再按燃烧区域三分支;最终汇成 Tad-Z 曲线,峰值落在化学计量点。

图解渲染中…
In输入:燃料与氧化剂种类 + 当地混合分数 ZStoichZ = Z_st 时 φ = 1,燃料与氧化剂恰好同时耗尽PeakTad 在 Z_st 处取最大,两侧因多余物质稀释而下降Burn完全燃烧假设:贫燃侧氧化剂有剩,富燃侧燃料有剩
11第 11 页 · 完全燃烧的温度与组分

有限反应速率的影响

左路:理想假设推出的状态关系。右路:现实里有限速率引发的四类偏离,汇成同一结论。

图解渲染中…
D2H、O、OH等自由基在反应区之外仍有非零浓度D3辐射散热或不完全燃烧使局部温度低于绝热值D4富燃侧氧化剂耗尽,部分燃料来不及燃尽就离开反应区
12第 12 页 · 有限反应速率的影响

Zeldovich变换

方程虽然写出来了,但仍是三维 PDE——输运和化学反应纠缠在一起。Zeldovich 的聪明之处:既然守恒标量本身就是 ξ,能不能把它直接当成新的自变量?

变量替换
把自变量从物理坐标 xᵢ 换成混合分数 ξ
变换形式
空间导数拆成 ξ 的梯度乘以沿 ξ 的导数
输运-化学解耦
∇ξ 负责输运,源项只跟 ξ 的局部值相关
PDE 退化为 ODE
多维偏微分方程降为沿 ξ 的一维常微分方程
三维地形里找位置对应 →沿海拔 ξ 这一个轴

x、y、z 还在,但描述标量场只需一个标尺——海拔 ξ

ϕxi=ξxidϕdξ\frac{\partial \phi}{\partial x_i}=\frac{\partial \xi}{\partial x_i}\frac{\mathrm{d}\phi}{\mathrm{d}\xi}
13第 13 页 · Zeldovich变换

变换的具体步骤

用一组数学步骤把 ξ 空间中的非线性燃烧问题,改造成可以分区解析求解的形式。

1
写出原方程
混合分数方程含指数型 Arrhenius 反应源项,难以直接求解
2
识别反应薄层
高 Damköhler 数下反应集中在 ξ=ξ_st 附近极薄区域内
3
施加 Zeldovich 变换
重新定义坐标与温度尺度,让源项在薄层外几乎为零
4
分区求解并匹配
薄层外按无限速率求解,薄层内按激活能渐进展开,再边界衔接
14第 14 页 · 变换的具体步骤

Burke-Schumann解析解

从对流扩散方程到温度组分剖面:左侧是起点,右侧是终点,中间每一步都是一次线性化简化。

图解渲染中…
B3Sh 取常数,扩散系数与速度的耦合被拆开B5用流函数把 PDE 化为 ODE,可严格积分B6Z=Zst 处即火焰面,把流场一分为二B7由状态关系式给出,不依赖流场细节
15第 15 页 · Burke-Schumann解析解

不同火焰结构的解

BS理论在三种几何下的特化。从中央框架分出三条分支,分别对应扩散项的不同形态。

图解渲染中…
B1反向流火焰,仅Z方向有梯度,扩散项消失C1圆柱坐标下,径向热质扩散起作用D1平面坐标下,y方向扩散决定分布
16第 16 页 · 不同火焰结构的解

Shvab-Zeldovich方法

前面我们用Zeldovich变换把反应项藏起来了。但藏起来之后,方程怎么解?Shvab-Zeldovich方法给了一个通用框架:先解守恒标量,再把化学装回去。

耦合函数
把组分与温度写成混合分数的线性组合,所有耦合函数都服从同一个方程
反应项抵消
化学计量比关系让非线性反应源项在耦合函数方程中自动消去
线性化求解
守恒标量方程化为纯对流-扩散方程,与化学反应解耦
双曲函数通解
对向流等典型构型下,Z的解自然呈现双曲函数结构
解多元方程组先消元再回代对应 →S-Z方法先求Z再装回化学

消去非线性反应项,只剩一个简单的Z方程;求出Z后,所有组分都按线性关系算出来

ρDZDt=(ρDZ)\rho \frac{DZ}{Dt} = \nabla \cdot (\rho D \nabla Z)
17第 17 页 · Shvab-Zeldovich方法

无限薄火焰面的近似

处理无限薄火焰面状态关系时,线性与分段线性两种思路的核心分歧在哪?

线性近似
  • T与Y_i对Z呈单一斜率直线
  • 隐含Le=1、无热损失假设
  • 燃料/氧化剂侧共用同一斜率
  • 形式简洁,参数少
分段线性近似
  • Z_st为折点,分两段直线
  • 可容纳Le≠1与热损失
  • 燃料/氧化剂侧斜率各自独立
  • 形式略繁但更贴近实际
理想绝热、Le=1条件下用线性近似更简洁;热损失或Le≠1场合必须升级到分段线性近似。
18第 18 页 · 无限薄火焰面的近似

课后思考

先独立思考,再对照参考答案;不必一次答全,重在培养物理直觉。

1混合分数的输运方程为何不含反应源项?这一性质在数值计算中带来什么好处?

参考答案输运方程无源项意味着 ξ 是纯被动标量;数值上只需解一个标量方程,就能确定流场中任意点的燃料与氧化剂配比。

2射流工况下燃料与氧化剂速度不等,Burke-Schumann 推导中的哪些假设需要修改?

参考答案等压等物性假设仍近似成立;需补充入口混合分数分布,并考虑流线曲率与拉伸对火焰面形状的修正。

3无限快速反应假设会在什么物理场景下严重失真?为什么?

参考答案点火延迟、淬熄、强拉伸等场景;此时有限速率效应不可忽略,火焰面不再是无厚度的几何面。

19第 19 页 · 课后思考