液滴燃烧的简化模型

官方工学老师·19 页·深入(追求细节与边界)·0 次浏览·2 天前
燃烧学d² 定律扩散火焰理论推导

液滴燃烧的简化模型

从准稳态假设到平方律推导,搞懂直径为何随时间线性消失

按 空格/→ 演示下一步

1 / 19 页

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

燃烧学d² 定律扩散火焰理论推导

液滴燃烧的简化模型

从准稳态假设到平方律推导,搞懂直径为何随时间线性消失

1第 1 页 · 液滴燃烧的简化模型

液滴燃烧的物理场景

液滴在高温环境中蒸发与燃烧的现象

液滴燃烧的物理场景
液滴在高温环境中蒸发与燃烧的现象
2第 2 页 · 液滴燃烧的物理场景

基本假设与简化条件

把真实燃烧过程压缩成可解方程组的四个关键简化前提。

1
准稳态
忽略液滴半径随时间的缓慢变化,每个瞬间按定常场求解
2
球对称
忽略重力和环境流动不对称,问题降为一维径向问题
3
Fick扩散
蒸汽输运仅由浓度梯度驱动,通量形式因此大幅简化
4
界面平衡
液面温度取湿球温度,蒸气分压等于该温度下的饱和压
3第 3 页 · 基本假设与简化条件

气相质量守恒方程

燃料蒸气从液滴表面往外扩散,氧气从外面往里扩散,二者碰面才燃烧。稳态下,每秒从表面蒸发的燃料量,必须等于火焰消耗的量——多一滴少一滴都不行。

球壳控制体
在半径 r 处取一薄层球壳做考察对象,做质量收支
稳态条件
进出球壳的质量流率相等,壳内存量不随时间变
球对称简化
只关心径向方向,其他方向自动消失
守恒结论
通过任意半径球面的燃料质量流率 ṁ 为常数
高速公路稳定车流对应 →球壳间的稳态扩散

没有匝道进出的路段,任何断面车流量都相等;球壳同理,任意 r 处燃料流率都相同

m˙=4πr2ρDdYFdr=const\dot{m} = 4\pi r^2 \rho D \frac{dY_F}{dr} = \text{const}
4第 4 页 · 气相质量守恒方程

质量守恒方程的积分

沿箭头方向,从微分方程逐步积分,最终得到火焰面处的质量通量关系。

图解渲染中…
A1∇·(ρv)=0 球坐标系下展开A3角度项与时间项为零,只剩径向A5ṁ=4πr²ρv_r 穿过任一球面为常数A7把r替换为火焰面位置r_f
5第 5 页 · 质量守恒方程的积分

组分守恒与火焰面条件

从气相组分方程出发,推导火焰面上燃料与氧化剂必须满足的界面条件。

1
组分守恒方程
对燃料与氧化剂各写一个独立的输运方程
2
火焰面假设
把火焰压成无厚度界面,反应只在面上发生
3
代入界面条件
火焰面两侧的通量差等于反应的消耗速率
4
求通量之比
推出燃料通量与氧化剂通量存在固定比值
5
确定火焰位置
火焰面就在燃料与氧化剂恰好全部反应处
6第 6 页 · 组分守恒与火焰面条件

气相能量方程

质量守恒告诉我们ṁ在稳态下是常数,组分守恒把火焰面钉在了某个 r_f。现在切到温度场——气相里这些热量究竟从哪儿来、怎么走、又会从哪儿跑掉?

守恒框架
对流输入 = 导热输出 + 辐射净损失
对流项
ṁ cp dT/dr,由前页ṁ常数结果直接得到
球坐标导热项
(1/r²) d/dr(r²λ dT/dr),Fourier 定律
辐射损失项
q̇_rad,高温或大液滴下不可忽略
边界条件
r→r_s 时 T=T_s,r→∞ 时 T=T_∞
暖气片给房间供热对应 →气相能量方程三项

热水沿管流入=对流项,金属壁向导热=导热项,墙面辐射散失=辐射项

ρucpdTdr=1r2ddr ⁣(r2λdTdr)q˙rad\rho u c_p \dfrac{dT}{dr}=\dfrac{1}{r^2}\dfrac{d}{dr}\!\left(r^2\lambda\dfrac{dT}{dr}\right)-\dot{q}_{rad}
7第 7 页 · 气相能量方程

温度分布的边界条件

气相能量方程给出一族可能的温度分布,但解不唯一——像只给斜率不给起点,画不出具体曲线。需要三个「锚点」把 T(r) 钉死。

液滴表面 r=r_s
T = T_s,由液-汽相平衡(Clausius-Clapeyron)决定
火焰面 r=r_f
T = T_f,达到绝热火焰温度,燃料在此耗尽
远场 r→∞
T → T_∞,回归未受扰动的环境
三根钉子拉直一根绳子对应 →三个温度边界条件

两端和中间高点固定后,曲线形状唯一

r=rs: T=Ts, r=rf: T=Tf, r: TTr=r_s:\ T=T_s,\ r=r_f:\ T=T_f,\ r\to\infty:\ T\to T_\infty
8第 8 页 · 温度分布的边界条件

温度分布曲线

沿径向 r 看温度:从液滴表面出发,经蒸气区到火焰面达峰,再到远处回落到环境温度。

图解渲染中…
a1T_d:液滴温度,约等于湿球温度或沸点a3T_f:火焰温度,化学反应放热致峰值a5T_∞:环境温度,远离液滴的基准
9第 9 页 · 温度分布曲线

液滴表面能量守恒

前面页画出的温度曲线在 r = r_s 处有个跳变:内部 T_d、表面 T_s、外侧迅速升至 T∞。但 T_s 究竟落到哪一格还没被显式规定——这一步要用表面能量守恒把它锁住。

表面是个能量汇
液体蒸发要吸走潜热 L_vṁ,这部分能量必须从外面传进来,否则表面温度会变
传入项:导热热流
气体侧在 r_s 处的导热通量 k_g(dT/dr)|_{r_s},经球面 4πr_s² 送进液滴
表面能量平衡
导热通量×球面积 = 单位时间蒸发质量×潜热,等式两侧量纲都是 W
气相方程的边界条件
和无穷远处 T∞ 配套,配上气相能量方程可锁定 T_s 与 ṁ;单一方程定不出两者
出汗时皮肤降温对应 →液滴表面蒸发吸热

汗蒸发的潜热来自身体→皮肤的导热,否则皮肤温度会变,液滴完全同构

4πrs2kgdTdrrs=m˙Lv4\pi r_s^2 k_g \left.\frac{dT}{dr}\right|_{r_s}=\dot{m}\,L_v
10第 10 页 · 液滴表面能量守恒

火焰面能量守恒

上一页我们看了液滴表面的能量守恒——蒸发吸热与导热输入的平衡。现在把视线推到火焰面,这里有另一种守恒:反应放出的热量必须等于向两侧导走的热量之和。

反应放热项
燃料蒸气与氧气在 r_f 处反应,每消耗单位质量燃料释放 ΔH_c 热量,构成火焰面的热源强度
内侧导热通量
λ_g(dT/dr)|_{r_f^-} > 0,方向指向液滴,为蒸发和预热燃料蒸气提供能量
外侧导热通量
λ_g(dT/dr)|_{r_f^+} < 0,方向指向远场来流,使氧化剂被预热至 T_f
能量跳跃条件
内、外侧导热之差 = 反应放热,温度梯度在火焰面发生有限跳跃
蜡烛火焰的悬浮位置对应 →火焰面能量守恒

烛焰放热 = 向下熔化蜡 + 向上加热空气;火焰位置由两路热流的平衡决定

m˙YF,fΔHc=λg[(dTdr)rf(dTdr)rf+]\dot{m}'' Y_{F,f}\Delta H_c = \lambda_g\left[\left(\frac{dT}{dr}\right)_{r_f^-} - \left(\frac{dT}{dr}\right)_{r_f^+}\right]
11第 11 页 · 火焰面能量守恒

液气平衡方程联立

把前面所有守恒方程像链条一样串起来,求出液滴的燃烧速率。

1
取温度场解 T(r)
从气相能量方程得到径向温度分布,作为后续热流通量的输入
2
求表面热流密度
对 T(r) 在 r=rs 处求导,得 −λ dT/dr|ᵣₛ
3
代入表面能量守恒
液滴侧能量平衡:导热通量 = ṁ·L,得到蒸发速率 ṁ
4
代入火焰面能量守恒
火焰两侧热流差 = ṁ·Y_{O,∞}·Q,给出第二个约束
5
联立质量守恒求 k
蒸发 ṁ = 扩散 ṁ,与火焰面方程配对解出 k,得 d²/dt(d²)=−k
12第 12 页 · 液气平衡方程联立

燃烧速率常数

K从前述方程组解出,看它如何归结到d²定律

图解渲染中…
B1半径方向的温度场T(r)C1半径方向的浓度场Y(r)D1火焰面无限薄反应假设E1单位时间液滴蒸发量
13第 13 页 · 燃烧速率常数

液滴寿命计算

上一页我们得到了燃烧速率常数 K,现在把它代回去积分,就能算出液滴从初始半径一路烧完要多久。一个反直觉的事实:不是直径本身、而是直径的平方随时间线性下降——这就是著名的 d² 定律。

d² 线性递减
d² 随时间以恒定斜率 −K 减小,叫 d² 定律
积分表达式
d²(t) = d₀² − Kt,标准直线方程
燃尽时间
令 d = 0,得 t_b = d₀²/K
平方放大效应
初始直径翻倍,寿命变 4 倍;用半径 r 同样形式,仅 K 差 4 倍系数
蜡烛燃烧对应 →液滴燃尽寿命

粗蜡烛比细蜡烛耐用得多——直径放大一点,寿命平方级延长

d2(t)=d02Kt,tb=d02/Kd^2(t) = d_0^2 - Kt,\quad t_b = d_0^2/K
14第 14 页 · 液滴寿命计算

对流效应的引入

前面所有公式都默认液滴静止、四周空气不流动——这就是纯扩散模型。现实中液滴要么在下落,要么周围热气在上升,对流不可避免。

纯扩散假设的边界
球对称+球壳火焰只在静止空气中成立
自然对流
火焰把空气加热→密度变小上升→冷空气从下方补入
强迫对流
液滴运动或外加风场强制带走产物、送来氧气
对流的双重影响
输运增强→燃烧更快;火焰被拉长倾斜、不再是球壳
无风中点香 vs 有风中点香对应 →纯扩散燃烧 vs 对流燃烧

无风烟笔直火苗稳;有风烟被吹歪、燃烧加快

15第 15 页 · 对流效应的引入

对流对d²定律的修正

沿箭头先看平方律如何从K₀变为Kc,再检查Kc的依赖与适用边界。

图解渲染中…
a5把流场效应折算进等效燃烧速率常数a10C为对流修正因子,不是新的普适常数a11低速、小温差时可用近似,但需评估误差a13准稳态、球形等前提仍须同时满足
16第 16 页 · 对流对d²定律的修正

d²定律应用例题

用一道典型例题走一遍 d² 定律求液滴寿命的全流程。

1
确定工况
燃料种类、初直径 d₀、环境温度与氧浓度
2
求燃烧常数 K
由 Spalding 数 B 与气液物性算出 K
3
写出 d² 定律
核心方程:d²(t) = d₀² − K·t
4
定终止条件
取 d=0 完全燃尽,或按工程要求取临界直径
5
解出寿命
t_b = d₀²/K,代入数值即得结果
17第 17 页 · d²定律应用例题

核心要点回顾

  • 球对称准稳态假设是全部推导的起点
  • 液面与火焰面边界条件耦合定解温度与浓度场
  • d²定律源于蒸发与燃烧的速率匹配
  • 燃烧常数K仅依赖物性,与初始条件无关
  • 对流效应通过输运数B修正纯扩散结果
延伸主题:多组分液滴的内部扩散高压条件下偏离球对称火焰液滴群燃烧的相互作用
18第 18 页 · 核心要点回顾

深入思考

先合上书自己想,每题想三分钟再看参考答案。

1d²定律说液滴寿命与初始直径平方成正比,这一结论的关键假设有哪些?

参考答案依赖准稳态、球对称、无强迫对流、蒸发系数为1等。这些假设忽略了液滴运动与传质阻力,故d²定律只是零级近似。

2实际喷雾中液滴大小不一,d²定律能直接用于估算整团液雾的燃烧时间吗?

参考答案不能直接套用。喷雾粒径分布很宽(Rosin-Rammler 分布),需按尺寸加权积分,总燃烧时间由最小液滴决定。

3若液滴在强剪切流中被拉成椭球,d²定律还成立吗?为什么?

参考答案基本不成立。模型基于球对称假设,椭球液滴蒸发面积与扩散路径都变了,火焰面也不再是球面,需重建立坐标系。

19第 19 页 · 深入思考
液滴燃烧的简化模型 · 知识图解