传质理论基础

官官方工学老师·17 页·深入（追求细节与边界）·0 次浏览·2 天前
✦ 传质理论基础工学传质理论基础菲克扩散定律
传质扩散Fick定律边界条件

传质理论基础

掌握扩散的物理图景、Fick定律推导与边界条件三层逻辑

按 空格/→ 演示下一步

1 / 17 页

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

传质扩散Fick定律边界条件

传质理论基础

掌握扩散的物理图景、Fick定律推导与边界条件三层逻辑

1第 1 页 · 传质理论基础

扩散是什么

把一滴墨水滴进一杯清水，墨滴会慢慢散开，整杯水最后变成均匀的淡灰色——这个过程没人推、没人搅，它自己就完成了。

微观本质
分子/粒子因热运动发生位置迁移，整体呈现从高浓度区向低浓度区的定向移动
驱动力：浓度梯度
∇c ≠ 0 才有扩散；浓度差越大，单位时间通过的物质量越多
自发过程
不需要外界做功，熵增驱动，沿梯度方向自发进行，最终趋于浓度均匀
与对流的区别
扩散靠分子随机热运动；对流靠外力造成的流体整体宏观流动
墨水滴入清水慢慢散开对应 →扩散

墨水分子从浓度高的墨滴处，自发迁移到浓度低的清水里，最后整杯颜色均匀

J=−D∇cJ = -D \nabla c
2第 2 页 · 扩散是什么

分子的随机热运动

上一页说扩散是「分子从高浓度跑到低浓度」，好像分子自带方向感。其实单个分子只是在乱撞——本页看它怎么乱、为何乱中有规律。

热运动随机性
室温下分子速度每秒几百米，但方向下一秒就变，瞬时路径完全不可预测
布朗运动
悬浮微粒被周围分子从各方向撞击，因瞬时受力不平衡而做无规则抖动
统计规律
单个轨迹混乱，但大量分子的均方位移⟨x²⟩与时间成正比
通向扩散
正是这种随机碰撞的统计结果，产生了宏观上可预测的扩散通量
醉汉随机踉跄对应 →布朗运动

每步方向随机不可预测，但均方位移与步数成正比

⟨x2⟩=2Dt\langle x^2 \rangle = 2Dt
3第 3 页 · 分子的随机热运动

微观运动→宏观浓度

从单个分子的无规则跳跃出发，经统计平均，涌现出宏观可测的浓度与通量关系。

图解渲染中…
a1：单个分子做随机游走，没有方向偏好a3：净通量：单位时间穿过单位面积的分子净数a5：Fick定律形式：J = -D∇c，梯度与通量线性挂钩
4第 4 页 · 微观运动→宏观浓度

菲克第一定律

上页看到无数分子的随机运动汇集成了宏观扩散。但扩散有多快、朝哪个方向走？菲克用一个简洁的关系式回答了这个问题。

稳态条件
浓度分布不随时间变化，每个截面通过的物质量恒定
扩散通量 J
单位时间通过单位截面的物质量，单位 mol/(m²·s)
浓度梯度 dC/dx
浓度沿空间方向的变化率，反映不均匀程度
扩散系数 D
物质的本征扩散能力，单位 m²/s，受温度、黏度等影响
负号的方向意义
扩散方向与梯度方向相反——始终从高浓度流向低浓度
走廊里一侧挤满人对应 →稳态扩散

人越挤走得越快，朝人少方向流动——通量正比于密度差，方向相反

J=−DdCdxJ = -D \dfrac{dC}{dx}
5第 5 页 · 菲克第一定律

浓度梯度的方向

上页我们看到菲克第一定律里有个 ∂C/∂x，这就是浓度梯度。它描述浓度随位置变化的快慢，但方向究竟指向哪里？这页讲清楚。

浓度梯度的定义
浓度随空间位置的变化率 ∂C/∂x，越大表示浓度变化越剧烈
梯度指向增大方向
数学上梯度永远指向函数值增加最快的方向，从低浓度指向高浓度
扩散方向与梯度相反
物质从高浓度向低浓度迁移，所以扩散方向 = −梯度方向
热力学本质
从高到低使化学势降低，系统自发趋于更稳定状态
山坡上滚动的小球对应 →浓度梯度与扩散

小球从山顶滚向山脚（扩散方向），坡度指向山顶（梯度方向），两者正好相反

J=−D∂C∂xJ = -D \frac{\partial C}{\partial x}
6第 6 页 · 浓度梯度的方向

菲克第二定律推导

从第一定律出发，对一个微元做质量守恒，就能推出浓度随时间变化的方程。

1
取微元
在 x 与 x+dx 之间取薄片，左右两侧通量不同
2
质量守恒
进入通量 − 离开通量 = 浓度时间变化率 × dx
3
泰勒展开
J(x+dx) ≈ J(x) + ∂J/∂x · dx，把进出差量化
4
代入第一定律
把 J = −D ∂C/∂x 代入 ∂J/∂x，得到 −D ∂²C/∂x²
5
得到结果
化简得 ∂C/∂t = D ∂²C/∂x²，即非稳态扩散方程
7第 7 页 · 菲克第二定律推导

边界条件与初始条件

求解扩散方程必需的三类条件

边界条件与初始条件
求解扩散方程必需的三类条件
8第 8 页 · 边界条件与初始条件

随机行走模型

上一页说分子在做无规则热运动，但菲克定律描述的扩散却是有方向的——中间缺一座桥。爱因斯坦的随机行走模型就是这座桥：从单个分子的'瞎走'出发，推出宏观扩散规律。

随机游走假设
分子每撞一次换一次方向，两次碰撞间走一步 λ（平均自由程），方向完全随机
均方位移
走 N 步后位移平均为 0，但位移平方的平均 ⟨x²⟩ = Nλ² 与步数成正比
扩散系数的微观身世
把 N = t/τ 代入，得 ⟨x²⟩ = 2Dt，D = λ²/2τ，从此 D 不再是经验参数
桥接菲克定律
对大量分子做统计平均，随机行走自然导出菲克第一定律
醉汉走夜路对应 →分子随机游走

每一步随机迈出，单步完全无序；但 N 步后位置分布可预测，方差正比于步数

⟨x2⟩=Nλ2=λ2τ t=2Dt\langle x^2 \rangle = N\lambda^2 = \frac{\lambda^2}{\tau}\,t = 2Dt
9第 9 页 · 随机行走模型

扩散系数的统计推导

从分子随机跳跃出发，用统计平均一步步推出扩散系数 D。

1
单步位移
每次跳跃方向随机、步长固定为 l
2
均方位移⟨r²⟩
对大量粒子位置平方取平均，刻画扩散的统计量
3
N步独立求和
⟨r²⟩=N⟨l²⟩，交叉项因方向对称相消
4
用时Nτ导出D
N步耗时 t=Nτ，得⟨r²⟩=6Dt，D=l²/6τ
10第 10 页 · 扩散系数的统计推导

爱因斯坦关系式

前页由随机行走得到三维扩散系数 D=〈r²〉/(6t)。若给同一种分子一个极小的外力，它还会沿受力方向缓慢漂移；爱因斯坦关系把这两个看似不同的运动参数连起来。

核心等式
对同一粒子与扩散坐标，D=μkBT；μ 是单位外力产生的漂移速度。
长时随机行走
三维各向同性随机行走在长时极限有 D=〈r²〉/(6t)，与前页统计结果一致。
球形粒子近似
球形、连续介质且无滑移时，ζ=6πηr，故 D=kBT/(6πηr)。
适用边界
这里要求线性响应、同一扩散坐标且记忆效应可忽略；团簇、界面流体会使比值偏离。
汽车自由滑行与轻推加速对应 →扩散与外力漂移

自由滑行反映热涨落，轻推后的速度响应反映迁移率；同一摩擦环境由 kBT 把二者关联。

D=μkBT=kBTζ,μ=vdF,D=⟨r2⟩6t,ζ=6πηr⇒D=kBT6πηrD=\mu k_{\mathrm B}T=\frac{k_{\mathrm B}T}{\zeta},\quad \mu=\frac{v_d}{F},\quad D=\frac{\langle r^2\rangle}{6t},\quad \zeta=6\pi\eta r\Rightarrow D=\frac{k_{\mathrm B}T}{6\pi\eta r}
11第 11 页 · 爱因斯坦关系式

微元体分析方法

上一页我们用随机行走与统计平均推导出爱因斯坦关系,那是微观视角。现在切换到连续介质视角——把扩散场切成无数微元体,在每个微元上列质量平衡,就能反推出菲克第二定律。

微元体选取
取dx为边长的小立方体,足够小可作连续介质,又含足够多分子
质量平衡通式
累积率 = 流入率 - 流出率 + 反应生成率,守恒律的微分表达
净流入通量
x与x+dx两面通量之差取极限,得到对流通量项∂J/∂x
代入菲克定律
把J=-D∂C/∂x代入通量项,即可从平衡过渡到扩散方程
水管上切一小段对应 →扩散场中取微元体

水流累积=流入-流出,与浓度累积在数学上完全同构

∂C∂t=−∂J∂x+RA\frac{\partial C}{\partial t} = -\frac{\partial J}{\partial x} + R_A
12第 12 页 · 微元体分析方法

组分的流入与流出

考察微元体两侧的通量差异，看清累积到底是怎么产生的。

1
选定微元体
沿 x 方向取厚度 dx、面积 A 的薄片作为控制体
2
左侧流入
左侧通量为 J，单位时间带入 J·A 的物质量
3
右侧流出
右侧通量为 J+∂J/∂x·dx，因梯度不同而与左侧有差
4
差值即累积
流入减流出等于微元体内部浓度随时间的变化率
13第 13 页 · 组分的流入与流出

组分守恒方程的推导

沿箭头看：守恒定律如何拆成对流项与扩散项，再拼成完整对流-扩散方程。

图解渲染中…
A2：单位时间净流入加源项生成量，等于积累量A4：流体宏观运动携带的组分通量，含速度场A5：浓度梯度驱动的扩散通量，由菲克定律给出A7：传质学核心方程，含时间项、对流项、扩散项、源项
14第 14 页 · 组分守恒方程的推导

有化学反应的守恒方程

守恒方程的框架已搭好，但反应会凭空生成或消耗物质——这要靠源项描述。对比看源项引入前后方程发生了什么变化。

无反应守恒方程
  • 物质仅靠扩散与对流进出微元体
  • 方程右端为零（无源无汇）
  • ∂c/∂t + ∇·N_i = 0
  • 适用于纯物理传质过程
有反应守恒方程
  • 物质除流动外，还由反应生成或消耗
  • 方程右端新增反应源项 R_i
  • ∂c/∂t + �·N_i = R_i
  • 适用于反应-扩散耦合过程
工程中绝大多数传质都伴随反应，应使用含源项的形式；无反应方程仅是理想化特例。
15第 15 页 · 有化学反应的守恒方程

核心要点回顾

  • ✓扩散是随机热运动在宏观尺度上的统计表现
  • ✓菲克定律是线性响应近似，浓差大时即失效
  • ✓守恒方程是骨架，源项与边界条件决定具体问题
  • ✓三大支柱层层嵌套，构成分析传质的完整链路
延伸主题：非菲克扩散与异常扩散多组分耦合传质对流传质与边界层
16第 16 页 · 核心要点回顾

课后思考

先独立思考再看参考答案，把卡住的地方记下来——这才是课后学习的真正起点。

1菲克第一定律说扩散通量正比于浓度梯度——但扩散的真正驱动力是什么？浓度梯度还是化学势梯度？

参考答案严格说是化学势梯度。浓度梯度只是它在理想稀溶液中的线性近似。当浓度高或存在强相互作用时，会出现"上坡扩散"——物质从低浓度向高浓度迁移。

2常温下氢气扩散系数约是氧气的 3 倍（6×10⁻⁵ vs 2×10⁻⁵ m²/s）。这个比值背后的物理本质是什么？

参考答案根源在分子量。气体动力学给出 D ∝ 1/√M：氢气 M=2，氧气 M=32，理论比值 √16 = 4，与实验值 3 倍接近。分子越小、热运动越剧烈，扩散越快。

3随机行走导出 ⟨x²⟩ = 2Dt。如果把布朗粒子从一维水箱换到二维或三维管道里，这个公式还成立吗？

参考答案仍成立，但形式变为 ⟨r²⟩ = 2nDt，n 是维数。三个方向各贡献 Dt，总位移随维数线性增加——这就是"为什么三维空间里迷路更容易"。

17第 17 页 · 课后思考
传质理论基础 · 知识图解