单颗粒碳的燃烧-单膜模型

官方工学老师·26 页·深入(追求细节与边界)·0 次浏览·2 天前
燃烧机理单膜模型传质扩散颗粒燃烧

单颗粒碳的燃烧:单膜模型

用一张图看懂氧扩散、反应界面与颗粒缩小的燃烧链条

按 空格/→ 演示下一步

1 / 26 页

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

燃烧机理单膜模型传质扩散颗粒燃烧

单颗粒碳的燃烧:单膜模型

用一张图看懂氧扩散、反应界面与颗粒缩小的燃烧链条

1第 1 页 · 单颗粒碳的燃烧:单膜模型

什么是单膜模型

想象一颗碳粒静静悬浮在高温空气中,氧气从四面八方赶来。但颗粒附近空气几乎不流动,氧气只能靠分子扩散穿过一层薄气膜到达表面——这就是单膜模型的核心:三维燃烧被压成「径向穿过薄膜」。

球对称假设
颗粒视为球体,所有物理量只随径向 r 变化,θ 和 φ 方向均匀
表面反应
反应只发生在颗粒外表面(异相反应),内部无燃烧
单层气膜
颗粒外存在一层极薄气膜,所有传质阻力都集中于此
准稳态
燃烧过程中气膜内的浓度场不随时间变化,仅为 r 的函数
一维扩散
膜内氧气只沿径向 r 方向扩散,三维问题简化为一维
饺子在沸水里煮熟对应 →单膜模型

热量先穿过饺子表面被冷却的水膜才能传进内部;氧气也先穿过碳粒外被消耗的气膜才能反应

NO2=DO2rp(CO2,CO2,s)N_{O_2}=\frac{D_{O_2}}{r_p}(C_{O_2,\infty}-C_{O_2,s})
2第 2 页 · 什么是单膜模型

单膜模型的物理图像

前页说到单膜模型假设反应极快、速率由扩散控制。这一页,我们把视线放到碳球本身——O2 究竟是怎么穿过周围气体、最终抵达球面的?

碳球外表面
碳粒是实心球,燃烧反应只发生在外表面一层
浓度边界层
球面外侧 O2 浓度从远处 C∞ 逐渐降到表面 0
等效「膜」
把渐变的浓度梯度压成厚度 δ 的虚拟薄壳,便于计算
膜外对流、膜内扩散
膜外由气流补充 O2,膜内仅靠分子扩散穿过 δ
围圈递饭的人墙对应 →碳球外的扩散膜

最内圈碗空相当于 O2≈0,最外圈饭足相当于 C∞,饭穿过人墙的速度即扩散速率

NO2=kcC,kc=D/δN_{O_2} = k_c \cdot C_\infty,\quad k_c = D/\delta
3第 3 页 · 单膜模型的物理图像

单膜模型的几何与变量

沿径向从中心往外看,分五层嵌套;箭头即 r 增大的方向。

图解渲染中…
C反应界面:O₂ 在此被碳消耗,浓度降为 0D气膜=边界层,O₂ 在内只靠分子扩散输送
4第 4 页 · 单膜模型的几何与变量

质量守恒的基本框架

上一页画出了单膜模型的几何——一颗碳球、一层贴附的气膜、外围远场气流。但'几何'只能告诉我们'长什么样',要回答'烧多快'还需要方程。质量守恒就是这个方程的起点。

控制体 CV
在颗粒外围画一个虚拟边界,包住颗粒和整层气膜,把这块空间当作'账本'来记账
累积项
控制体内总质量随时间的变化率 dmcv/dt,非稳态时不为零
流动项
O2、N2 等随气流从远场流入,CO2、CO 等产物从表面离开
源汇项
表面反应消耗碳为'汇',同时生成 CO2/CO 等气体为'源'
稳态假设
通常假设颗粒缓慢燃烧,dm/dt≈0,方程简化为流入+生成=流出+消耗
家里的水桶对应 →单颗粒燃烧的控制体

水龙头进水=空气流入;桶底出水=产物气流出;舀水喝掉=碳消耗;反应产气=额外冒泡

dmcvdt=m˙inm˙out+m˙genm˙cons\frac{dm_{cv}}{dt}=\dot m_{in}-\dot m_{out}+\dot m_{gen}-\dot m_{cons}
5第 5 页 · 质量守恒的基本框架

一维稳态扩散方程

从稳态质量守恒出发,按流程一步步推出扩散-反应微分方程。

图解渲染中…
A稳态:浓度场不随时间变化,流入+生成=流出+消耗B在膜内任意位置取厚度dx的薄控制体FTaylor展开:把x+dx处的通量在x处展成级数Hd/dx(D dC/dx) = -k C,即一维稳态扩散-反应ODE
6第 6 页 · 一维稳态扩散方程

氧气与CO₂的双组分分析

上一页我们写下了单个组分的一维稳态扩散方程。但碳燃烧涉及两种气体——O₂ 进去、CO₂ 出来——它们在膜内同时反向运动,互相牵制。这就是双组分耦合问题。

O₂ 向内扩散
氧气从气相主体穿过气膜,输送到颗粒表面供反应消耗
CO₂ 向外扩散
反应产物CO₂从颗粒表面穿过气膜,排入气相主体
等摩尔反向耦合
C+O₂→CO₂,每耗1 mol O₂生成1 mol CO₂,两股摩尔流大小相等方向相反
无总体流动
正反两向摩尔流相互抵消,膜内无净的 Stefan 流动,纯靠分子扩散传质
双向自动扶梯对应 →O₂与CO₂的反向扩散

上楼的人=O₂进燃烧面,下楼的人=CO₂离开,单位时间人数相等,净人流为零

$N_{O_2} = -\,N_{CO_2}$
7第 7 页 · 氧气与CO₂的双组分分析

组分守恒的边界条件

从远场到表面,看O₂与CO₂两类边界条件如何衔接。

图解渲染中…
a1无穷远处未受扰动的主体气相a4颗粒-气体界面,反应在此发生a5化学计量1:1:消耗1 O₂生成1 CO₂a6O₂进入通量等于CO₂离开通量
8第 8 页 · 组分守恒的边界条件

燃烧速率的表达式

扩散方程描述氧气怎么穿过气膜到表面,边界条件告诉我们表面反应多快——现在要把这两端扣在一起:扩散送来的氧必须刚好够反应消耗,稳态才成立。

耦合条件
稳态下,扩散到表面的氧通量等于表面反应消耗的氧通量
双驱动力
浓度差 (C∞-Cs) 驱动扩散,Cs 又被表面动力学反过来钉住
两种极限
k_s·δ/D ≫ 1 时反应控制;≪ 1 时扩散控制
速率表达式
总速率对两种阻力取调和平均,正比于外部浓度 C∞
工厂供料与产线节拍对应 →扩散-反应耦合

原料运输速度决定车间最大产能;二者必须匹配才能稳态运行

$\dot{r}_{O_2} = \dfrac{D/\delta \cdot k_s}{D/\delta + k_s}\, C_\infty = \dfrac{C_\infty}{\delta/D + 1/k_s}$
9第 9 页 · 燃烧速率的表达式

碳-氧表面反应机理

C + O₂ → CO₂ 的表面反应步骤

碳-氧表面反应机理
C + O₂ → CO₂ 的表面反应步骤
10第 10 页 · 碳-氧表面反应机理

反应速率的阿伦尼乌斯表达

前页我们说到 C + O₂ → CO₂ 这条机理路线——但反应到底有多快?这就要靠表面反应速率常数 k_s 来刻画。阿伦尼乌斯公式把温度 T 与 k_s 用一条指数关系牢牢锁住。

阿伦尼乌斯公式
k_s = A·exp(−E/RT),温度与反应速率的桥梁
指前因子 A
反映分子有效碰撞频率与取向概率,碳燃烧典型 10~10⁴ m/s
活化能 E
反应必须越过的能量门槛,碳氧反应约 100~180 kJ/mol
指数敏感性
E 越大、T 越高,k_s 对 T 变化越敏感,决定反应处于化学还是扩散控制区
登山翻越垭口对应 →阿伦尼乌斯反应速率

尝试人数对应 A,垭口海拔对应 E,登山者平均体力对应 T,翻过去的人次就是 k_s

ks=Aexp(E/RT)k_s = A \cdot \exp(-E/RT)
11第 11 页 · 反应速率的阿伦尼乌斯表达

动力学控制 vs 扩散控制

看图方法:从判断框架出发,左右对比两种极限机制与中间过渡区。

图解渲染中…
Jks随温度指数变化,是判断核心KC2表观活化能≈真实活化能DC2表观活化能≈0,温度不再敏感
12第 12 页 · 动力学控制 vs 扩散控制

为什么要用电路比拟

上一页已经用'阻力'来描述扩散和反应的快慢对比。把'阻力'再推一步:它和电路里的电阻在数学上完全是一回事。

形式同构
扩散通量公式与欧姆定律数学结构完全一致
三要素映射
通量↔电流,浓度差↔电压,阻力↔电阻
串联叠加
多层扩散阻力像电阻一样相加
直观优势
复杂传质可视化,便于识别限速环节
适用范围
仅稳态一维线性扩散成立;非线性需修正
串联电阻电路对应 →多层扩散传质

电压降累加对应浓度降累加;总电阻决定总电流对应总阻力决定总通量

NA=cbcsRdI=ΔVRN_A=\dfrac{c_b-c_s}{R_d}\quad\longleftrightarrow\quad I=\dfrac{\Delta V}{R}
13第 13 页 · 为什么要用电路比拟

扩散阻力 ≡ 欧姆电阻

左右两条平行的因果链:上为菲克定律,下为欧姆定律。三组节点一一对应——这就是「电路比拟」的本体论根基。

图解渲染中…
f1膜外 C_g 与表面 C_s 之差,扩散的驱动力f21/(k_c·A),传质阻力的量化形式o1电路两端电压差,电流的驱动力o2R_elec,材料的本征导电阻力
14第 14 页 · 扩散阻力 ≡ 欧姆电阻

串联电阻模型

O₂ 依次穿过两层阻力,对应两个串联的电阻。

图解渲染中…
a2外扩散阻力 R_ext = 1/k_ext,气膜越厚越大a4表面反应阻力 R_surf = 1/k_s,温度越低越大a3表面浓度 C_s 由两段阻力共同决定
15第 15 页 · 串联电阻模型

总燃烧速率的电路解

上一页我们把扩散阻力和反应阻力画成了两个串联的电阻。现在问题来了:怎么像解电路那样,一口气求出总燃烧速率?

串联电阻求和
总阻力等于扩散阻力加反应阻力:R总 = Rd + Rr
驱动浓度差
无穷远处到表面的氧气浓度差,就是电路里的「电压」
欧姆定律形式
摩尔通量 N = 浓度差 / 总阻力
统一速率表达式
N = C_O2,∞ / (1/kd + 1/kr),kd、kr分别是传质与反应速率常数
极限行为一致
哪个阻力大,总速率就由哪一边决定——与上一页极限分析吻合
串联电阻电路对应 →燃烧速率求解

电压↔浓度差,电流↔摩尔通量,电阻↔传质+反应阻力;总电流=总电压/总电阻

NO2=CO2,1/kd+1/kr=ΔCRd+RrN_{O_2}=\dfrac{C_{O_2,\infty}}{1/k_d+1/k_r}=\dfrac{\Delta C}{R_d+R_r}
16第 16 页 · 总燃烧速率的电路解

不同控制机制的物理图像

两种控制机制都限制了燃烧速率,但瓶颈位置不同——电路里要看哪个电阻在主导总电流。

动力学控制
  • 慢在表面反应这一步
  • R_r 占总电阻绝大部分
  • 速率由 k_s(T) 主导
  • 升温明显加快
扩散控制
  • 慢在边界层气膜扩散
  • R_d 占总电阻绝大部分
  • 速率与 k_s 无关
  • 升温几乎不变
低温小颗粒或慢反应→动力学;高温大颗粒或快反应→扩散。判据:直接比较 R_r 与 R_d 谁更大。
17第 17 页 · 不同控制机制的物理图像

能量守恒方程的形式

上一页用电阻网络求出了总燃烧速率,但还差一块:颗粒温度 $T_s$ 到底是多少?要回答这个问题,必须给颗粒再立一条方程——能量守恒,把气体导热与反应放热合成一笔账,与温度变化平账。

三方能量收支
气体加热 + 反应放热 = 颗粒内能积累,是骨架方程
加热项
$hA_p(T_\infty-T_s)$:膜内导热给颗粒的热流,温差驱动
反应热项
$(-\Delta H)\dot{r}_s A_p$:表面异相反应放热,$\dot{r}_s$ 由前几节定
准稳态化简
$dT_s/dt=0$:积累项消失,$T_s$ 由代数方程唯一确定
蓄水池对应 →颗粒内能

进水(加热+反应热)=水位(温度)变化;稳态时水位锁定,自动平衡

ρpcp,pVpdTsdt=hAp(TTs)+(ΔH)r˙sAp\rho_p c_{p,p} V_p \frac{dT_s}{dt} = hA_p(T_\infty-T_s) + (-\Delta H)\dot{r}_s A_p
18第 18 页 · 能量守恒方程的形式

表面反应热与热释放

上页能量守恒里那个 Q 还裹着面纱——它就是反应焓。这一页拆开:C + O₂ → CO₂ 释放多少热、在哪释放、强度几何。

反应焓 ΔH
标准反应 C + O₂ → CO₂ 释放 ΔH ≈ -393.5 kJ/mol,是典型强放热反应
表面释放
反应只在颗粒外表面进行,热也在表面放出,气相主体无内热源
能量方程源项
单位面积放热率 q = Rc·(-ΔH),是颗粒升温的能源
数量级震撼
-393.5 kJ/mol 极大;Rc = 10⁻⁶ mol/(cm²·s) 即可释放 ~0.4 W/cm²
瞬时释放假设
化学时间 << 扩散时间,热视为在表面瞬时放出,无须额外动力学建模
瀑布水位落差对应 →反应焓 ΔH

水从高势能落到低势能,差额变热;原子从高能键重排到低能键,差额也变热

C+O2CO2,ΔH=393.5  kJ/mol\text{C} + \text{O}_2 \to \text{CO}_2, \quad \Delta H^\circ = -393.5\;\text{kJ/mol}
19第 19 页 · 表面反应热与热释放

颗粒温度分布特征

上一页讲到反应放热集中在颗粒外表面——这股热量除了传给气体,还会往哪走?答案是颗粒内部。表面是热源,同时向气体和核心两个方向传热。

表面是热源
反应仅在颗粒外表面进行,热量同时向气体与颗粒内部两个方向传导
表面比核心热
稳态下 T_s > T_c,温差由内部导热热阻与外部对流热阻共同决定
Biot 数刻画温差
Bi = h·r_p/k_s,Bi 越小颗粒越接近等温,Bi→0 时 T_s≈T_c
碳颗粒近似等温
碳的 k_s 很高,Bi<<1,ΔT 极小,单膜模型可用 T_p 一个温度描述
燃烧的火柴对应 →燃烧的碳颗粒

火柴也是表面先烧,火焰在外、热量往里传,未燃部分(核心)始终比表面凉

Bi=hrpks\text{Bi} = \dfrac{h \, r_p}{k_s}
20第 20 页 · 颗粒温度分布特征

温度与燃烧速率的耦合

温度与燃烧速率互为因果,形成自催化正反馈环。

1
颗粒温度升高
受外部供热或前期反应放热,颗粒表面温度T_s开始上升
2
速率指数跃升
k∝exp(-E/RT),T_s略升即让表面反应速率急剧放大
3
表面氧更贫
动力学加快使碳表面更'饿',O₂摩尔分数Y_O₂,s显著降低
4
扩散梯度加大
Y_O₂,∞与Y_O₂,s差值拉大,扩散驱动力随之增强
5
总速率加快
电路解q_total=ΔY/(R_diff+R_kin)显著增大
6
更多热量释放
ΔH_rxn·q_total更大,反哺颗粒使T_s继续上升,回到st1
21第 21 页 · 温度与燃烧速率的耦合

Sherwood数的物理意义

前几页的单膜模型默认气体整体静止,氧气靠分子扩散抵达碳表面。但炉膛里有气流吹过颗粒,对流把氧气更快推过来。Sherwood数就是这种加速相对于纯扩散的比值。

Sh定义
Sh = k_c·L/D_AB,对流传质系数×特征长度÷分子扩散系数
物理意义
实际传质通量(含对流)与纯扩散通量的比值
无对流极限
Re→0时Sh→2,是球体在静止流体中纯扩散的精确解
工程关联式
Sh = 2 + 0.6·Re^0.5·Sc^1/3(Ranz-Marshall)
厨房油烟扩散对应 →Sherwood数的物理含义

关油烟机=纯扩散(Sh≈2),开油烟机=对流扩散(Sh>2),加号右侧即加速量

Sh=kcLDAB=J实际J扩散Sh = \dfrac{k_c \cdot L}{D_{AB}} = \dfrac{J_{\text{实际}}}{J_{\text{扩散}}}
22第 22 页 · Sherwood数的物理意义

Damköhler数的控制分区

上一页我们把反应与扩散画成两根串联电阻,总阻力由两者共同决定。但工程师必须追问:瓶颈到底在哪一根?这就需要 Damköhler 数这个无量纲判据。

Da 的定义
反应速率与扩散速率之比,无量纲判据
Da ≪ 1 动力学控制
反应慢、扩散快,瓶颈在表面反应
Da ≫ 1 扩散控制
反应快、扩散慢,表面浓度趋于零
Da ∼ 1 过渡区
两根阻力相当,必须同时求解
对应电路模型
Da = R_diff / R_rxn,比值即哪根主导
漏斗倒水的瓶颈对应 →Damköhler 数判据

漏斗口径 = 反应速率;上方补给 = 扩散速率;Da 大则出口宽补给慢,瓶颈在补给

Da=ksdpDABSh=RdiffRrxnDa = \dfrac{k_s \, d_p}{D_{AB} \, Sh} = \dfrac{R_{diff}}{R_{rxn}}
23第 23 页 · Damköhler数的控制分区

单膜模型核心检验

点击作答

在单膜模型的Damköhler数分析中,当Da远小于1时,燃烧过程的控制机制和颗粒表面氧浓度分别是?

24第 24 页 · 单膜模型核心检验

单膜模型知识地图

  • 串联电阻是统一框架:总阻力 = 扩散阻力 + 反应阻力
  • Da 数分区判据:Da→0 纯动力学,Da→∞ 纯扩散
  • 温度正反馈:升温加速反应同时抬高Da,使控制机制漂移
  • Sherwood 数让结论无量纲化,可推广到不同粒径与流场
  • 模型边界:忽略孔内扩散,不适用于多孔碳(需用缩核模型)
延伸主题:多孔碳的缩核模型颗粒群燃烧的相互作用实际煤粉燃烧的边界层修正
25第 25 页 · 单膜模型知识地图

延伸思考与开放问题

三道题,先自己思考再看参考答案。覆盖核心回顾、参数推演、模型边界。

1为什么高温下燃烧由扩散控制,低温下由反应动力学控制?这个转变的物理本质是什么?

参考答案Da数揭示本质:Da>>1时表面反应极快,氧在抵达前耗尽,速率受扩散限制;Da<<1时表面反应是瓶颈。转变温度由活化能与颗粒尺寸共同决定。

2如果把颗粒直径缩小到原来的1/10,扩散控制与动力学控制区的边界温度会如何变化?

参考答案扩散路径缩短、氧气供应增强,但表面反应速率基本由温度决定。结果是Da增大,动力学控制区间扩大,转变温度向高温方向移动。

3单膜模型假设球对称稳态扩散。若考虑气流冲刷导致的对流,模型还成立吗?需要怎么修正?

参考答案不再完全成立。需用Sh=2+0.6·Re^0.5·Sc^0.33取代Sh=2,把强制对流的增强项显式纳入扩散阻力,Re很大时进入对流控制新极限。

26第 26 页 · 延伸思考与开放问题