液滴蒸发的简单模型

官官方工学老师·18 页·深入（追求细节与边界）·0 次浏览·2 天前
✦ 液滴蒸发的简单模型工学液滴蒸发的简单模型基本假设
液滴蒸发d²定律传热传质燃烧基础

液滴蒸发的简单模型

看清 d² 定律背后藏着哪些前提——传热、传质与准稳态假设的适用边界

按 空格/→ 演示下一步

1 / 18 页

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

液滴蒸发d²定律传热传质燃烧基础

液滴蒸发的简单模型

看清 d² 定律背后藏着哪些前提——传热、传质与准稳态假设的适用边界

1第 1 页 · 液滴蒸发的简单模型

液滴蒸发无处不在

蜡烛在烧、咖啡在凉、墨水在干、发动机在轰鸣——这四件看似不相关的事，背后都在做同一件事：把液体悄悄变成气体。

蜡烛火焰
熔化的蜡液滴被火焰加热蒸发，气态石蜡与氧气反应放热
咖啡变凉
液面水分子逃逸成蒸汽，带走热量使杯中液体降温
喷墨打印
压电喷出微米级墨滴，溶剂蒸发后染料固定在纸上
航空发动机
燃油雾化成百微米液滴，在燃烧室毫秒级蒸发后点火
同一首歌的不同演奏对应 →液滴蒸发的普遍规律

钢琴、提琴、吉他各奏同一旋律——蜡烛、咖啡、发动机讲同一个蒸发故事

D2(t)=D02−K⋅tD^2(t) = D_0^2 - K \cdot t
2第 2 页 · 液滴蒸发无处不在

为什么需要假设

上一页我们看到液滴蒸发无处不在——但真实场景里，温度、流场、蒸汽浓度、表面张力梯度……全部纠缠在一起。要把这些全写进方程，没人能解得动。我们必须先做一件事。

现实是多场耦合
温度场、速度场、浓度场、表面变形同时演化，方程一旦耦合就几乎无法解析
假设 = 选择性忽略
保留主导物理、丢掉次要效应，把复杂现实压成可写方程的形式
简化后才能求解
方程简化后才有可能给出解析解或低维数值解
假设有适用边界
模型只在假设成立的范围内有效；越界，预测就偏离现实
城市地图对应 →物理简化模型

地图不画每棵树每条裂缝，但保留了街道走向——模型不写所有细节，但保留主导物理过程

3第 3 页 · 为什么需要假设

四大基本假设

建蒸发方程太难，先把真实场景拆成四块可解的拼图。

1
球形
液滴始终保持完美球形，靠表面张力抵抗形变
2
准稳态
温度和蒸汽浓度场瞬时跟上界面退缩速度
3
一维
传热传质只沿径向发生，角向没有差异
4
静止环境
忽略对流与外流场，只有分子扩散搬运质量
4第 4 页 · 四大基本假设

假设的物理图景

由内向外看：从液滴内部到无限远处的空间分层结构，每一层对应一个物理主导过程。

图解渲染中…
a1：球内温度均匀，无内部流动a2：蒸发发生在界面，法向通量为Ja3：蒸气靠扩散向外输运a4：浓度c∞、温度T∞保持不变
5第 5 页 · 假设的物理图景

蒸汽扩散的本质

上页假设蒸汽在静止空气中靠扩散离开液滴——那「扩散」到底是什么？打开香水瓶，几十秒味道飘到几米外，没有风扇，靠的纯粹是分子自己乱跑。

分子随机热运动
每个分子时刻做无规则热运动，朝各方向走的概率均等
浓度梯度产生净通量
单个分子去哪都行，但高浓度侧单位体积分子更多，穿过截面的净分子数从高到低
没有「推力」
扩散不是被什么力推出来的，是高浓度侧分子数更多这一统计后果
菲克第一定律
通量 J = -D∇c，扩散系数 D 反映分子随机运动的快慢
拥挤走廊里的人流对应 →高浓度到低浓度的分子扩散

没人被规定向外走，但人多一侧穿过任何截面的人数天然更多

J=−D∇cJ = -D \nabla c
6第 6 页 · 蒸汽扩散的本质

菲克定律与浓度梯度

上一页我们看到水汽从液滴表面往外散开。'散开有多快'不能靠直觉——菲克定律把通量与浓度梯度用一行等式绑在一起。

菲克第一定律
1855 年提出：稳态下通量正比于浓度梯度
质量通量 J
单位时间穿过单位面积的质量 kg/(m²·s)
扩散系数 D
蒸汽的本征属性，约 2.5×10⁻⁵ m²/s
浓度梯度 ∇c
浓度随距离的变化率，负号表方向:高→低
斜坡上的水流对应 →浓度梯度驱动扩散

坡度对应梯度，流速对应通量；都从'高'流向'低'

J=−DdcdxJ = -D \dfrac{dc}{dx}
7第 7 页 · 菲克定律与浓度梯度

浓度分布的建立

从初始阶跃到稳态分布，浓度场经历几步演化才稳定下来。

1
初始阶跃
t=0 时液面为饱和蒸汽浓度 Csat，远场为零，分布是阶跃突变
2
蒸发启动
液面持续向环境释放蒸汽分子，表面浓度得以维持
3
扩散输运
浓度梯度驱动蒸汽由高到低定向迁移，菲克定律发挥作用
4
浓度场建立
随时间推进，表面到远场形成连续递减的浓度分布曲线
5
准稳态达成
分布形状不再变化，但通量仍在持续，日常分析多针对此阶段
8第 8 页 · 浓度分布的建立

d²定律：蒸发速率规律

浓度场稳定了，蒸汽还在不停从液滴表面向外扩散。那液滴自己会怎么缩小？匀速、越来越快、还是越来越慢？菲克定律给出了一个近乎优雅的答案。

d² 线性递减
直径的平方随时间呈直线下降，是蒸发进入准稳态的直接标志
斜率 K 的构成
K = 8D(Cs−C∞)/ρl，由扩散系数、浓度差、液体密度共同决定
球面几何的对冲
dm/dt ∝ d 撞上 m ∝ d³，求导后 d²/dt 自然归一为常数
适用与失效边界
仅在扩散主导、球对称、准稳态下成立；自然对流、Stefan 流、高蒸发通量都会使其偏离
圆柱水杯匀速排水对应 →液滴 d² 线性递减

水杯截面不变→水位匀速下降；液滴球面扩散让 d² 也匀速下降，殊途同归

d2(t)=d02−Kt,K=8D(Cs−C∞)ρld^2(t) = d_0^2 - Kt, \quad K = \dfrac{8D(C_s - C_\infty)}{\rho_l}
9第 9 页 · d²定律：蒸发速率规律

d²定律的物理图像

从左到右读图，看液滴在五个阶段中直径平方逐步递减直至消失。

图解渲染中…
a1：初始时刻，液滴直径为d₀，周围空气几乎无蒸汽a3：扩散层稳定后，进入d²定律适用的准稳态a4：核心阶段：d²=d₀²-Kt，直径平方随时间线性下降a5：d→0时蒸发速率趋于零，液滴完全消失
10第 10 页 · d²定律的物理图像

界面的能量博弈

前面讲的扩散解决了「蒸汽怎么走」，但有个根本问题没解决：液滴变蒸汽需要吸收大量潜热，这股能量从哪来？答案是热传导——气相中的热量必须流到界面来。

热通量抵达界面
气相中靠热传导把热量送到液滴表面
潜热消耗
液体分子挣脱液面变成蒸汽，必须吸收汽化潜热 L
收支平衡
稳态下：到达界面的热通量 = 蒸发带走的潜热速率
温度梯度被锁定
这个平衡条件反过来决定了界面附近的温度梯度大小
月光族的工资对应 →界面的能量收支

工资（热通量）每月花光（潜热），存不下余额——稳态要求收支刚好相抵

kdTdx∣R=L⋅m˙k\left.\dfrac{dT}{dx}\right|_{R} = L \cdot \dot{m}
11第 11 页 · 界面的能量博弈

热流路径分析

热量从环境传到液滴表面要经过几步，每一步都可能成为蒸发速率的瓶颈。

1
环境对流
远处热空气通过对流把热量带到液滴边界层外缘
2
边界层导热
热量穿过气相温度边界层传导到液滴表面
3
界面潜热
表面吸收热量作为汽化潜热，驱动液-气相变
4
内部补充
液滴内部从高温区向表面导热（简单模型中 Td=Tw 而略去）
12第 12 页 · 热流路径分析

蒸发潜热的角色

上一页的热流最终汇向界面，但界面不是终点——它把能量吞了：每一个分子从液态变气态，都要带走一份「潜热」。今天展开这一步：相变消耗的能量到底有多厉害。

相变是能量黑洞
每蒸发 1 kg 液体，界面必须吸走 L 焦耳热量，相当于这公斤物质欠下的能量债
潜热数值惊人
水的 L≈2.4 MJ/kg，相当于把 1 kg 0℃ 水一路烧到 100℃，再烧到近 700℃
蒸发＝自冷却
能量被不断抽走，界面温度降到湿球温度 T_s<T_∞，温差反过来又驱动向内传热
潜热耦合两类方程
没有这部分能量支出，扩散方程与导热方程各管各的、互不相干，模型不再自洽
炉子上的烧水壶对应 →界面能量收支

火给壶的热功率决定水烧多快；环境给界面的热通量决定蒸发多快。火撤了水壶就停，潜热没补上蒸发就停

m˙ L  =  4πR2 k (T∞−Ts)\dot{m}\,L \;=\; 4\pi R^{2}\,k\,(T_\infty - T_s)
13第 13 页 · 蒸发潜热的角色

液滴寿命公式推导

从质量守恒出发，五步推导出液滴寿命 τ = r₀²/K。

1
质量守恒
dm/dt = −4πr²·J，质量守恒的微分形式
2
用半径表示
m = 4πr³ρ_l/3，把对 m 的方程换成对 r 的
3
联立化简
两式消去 4πr²，得到 ρ_l·dr/dt = −J
4
代入扩散通量
把扩散模型的 J 代入，得到 d(r²)/dt = −K
5
积分得寿命
r² 从 r₀² 线性衰减到 0，τ = r₀²/K
14第 14 页 · 液滴寿命公式推导

寿命与初始直径的关系

从 d² 定律出发，追随初始直径流向寿命公式，看直径变化如何非线性放大寿命。

图解渲染中…
a2：直径平方随时间线性衰减a3：令 d=0 得总寿命，与 d₀² 成正比
15第 15 页 · 寿命与初始直径的关系

模型边界与局限性

四条假设搭建的理想模型，遇到真实燃料液滴几乎条条偏离——看清边界在哪。

理想模型预测
  • d²定律严格成立，蒸发速率恒定
  • 液滴内部温度均匀，无内部流动
  • 球对称扩散，无变形与Marangoni效应
  • 单组分蒸发，无相间平衡偏移
实际液滴偏差
  • 多组分蒸发，轻组分先挥发，d²不再线性
  • 内部环流导致温度不均，加速蒸发
  • 变形、破裂、Marangoni对流破坏球对称
  • 相变温度与活度修正，蒸汽压非平衡
理想模型是定性分析利器；工程量化须补多组分、Marangoni、非平衡相变三块拼图。
16第 16 页 · 模型边界与局限性

液滴蒸发模型全貌

  • ✓四假设勾勒理想图景，d²定律由此自然涌现
  • ✓质量扩散与热量传递在界面耦合，缺一不可
  • ✓寿命正比于初始直径平方，是模型标志性预言
  • ✓简化既赋予模型力量，也划定其适用范围
延伸主题：考虑马兰戈尼对流的蒸发多组分液滴的蒸发竞争环境湿度对寿命的影响
17第 17 页 · 液滴蒸发模型全貌

课后思考

先合上笔记，对每道题想两分钟，再看参考答案——重点不是对错，而是你思考时卡在了哪里。

1d² 定律为什么是线性的？回忆推导过程，什么物理量决定了这条直线的斜率？

参考答案根源在扩散方程的解：蒸汽通量与剩余直径成正比，所以 d² 随时间线性变化。斜率 K 由扩散系数 D、饱和浓度差和液体密度共同决定——K 越大，蒸发越快。

2如果把环境从氮气换成氦气（氦气中蒸汽扩散系数更大），液滴寿命会怎样变化？哪些量也跟着变？

参考答案氦气中 D 更大，d² 定律里的斜率 K 跟着变大，液滴寿命缩短。但前提是准稳态、均匀温度等假设仍成立；若变化太剧烈，这些假设会率先被打破。

3对于直径只有几十纳米的液滴，这个模型还成立吗？四条基本假设中哪一条最先失效？

参考答案纳米尺度下连续介质假设失效——分子数太少，浓度概念失去意义，克努森效应显著，蒸发通量不再遵循菲克定律。需要分子动力学或过渡区理论来修正。

18第 18 页 · 课后思考