菲克第一定律：稳态扩散的数学描述
各位，咱们今天聊点实在的。扩散这个东西，说白了就是原子在材料里"乱跑"。但你别小看这个"乱跑"，它决定了你做的合金能不能均匀、半导体掺杂能不能到位、甚至电池用久了会不会衰减。
我个人习惯，讲扩散一定先从菲克第一定律入手。为啥？因为它最简单、最直观，而且很多工程问题都能用它来估算。你想想看，如果连稳态扩散都搞不明白，那后面那些非稳态、反应扩散就更别谈了。
浓度梯度与扩散通量
先问大家一个问题：为什么气体会从高压区往低压区流动？为什么墨水会在水里散开？
答案很简单——因为有"势"的差异。在扩散里，这个"势"就是浓度梯度。
菲克第一定律的数学形式，我直接写出来：
J = -D · (dC/dx)这里：
- J —— 扩散通量，单位是 atoms/(m²·s) 或 kg/(m²·s)。说白了就是每秒钟有多少物质穿过单位面积。
- dC/dx —— 浓度梯度，单位是 atoms/(m⁴) 或 kg/(m⁴)。它描述的是浓度在空间上的变化快慢。
- D —— 扩散系数，单位是 m²/s。这个后面重点讲。
- 负号 —— 表示扩散方向是从高浓度往低浓度走。这是热力学第二定律的体现，别问我为什么，大自然就是这么安排的。
我在项目中遇到过一件事：有次做铜铝复合板的界面分析，发现界面附近成分分布不对称。一开始以为是工艺问题，后来一算，原来是忽略了浓度梯度在界面处的突变。嗯，这里要注意，浓度梯度不是连续的，尤其是在相界面处。
核心要点：菲克第一定律只适用于稳态扩散。什么叫稳态？就是浓度分布不随时间变化，即 ∂C/∂t = 0。说白了，你测出来的浓度曲线今天和明天一模一样。
扩散系数的物理意义
扩散系数 D，这个参数太重要了。我经常跟年轻工程师说：你要是能把这个 D 搞明白，扩散这块你就通了八成。
D 的物理意义是什么？它衡量的是原子在基体中"跳来跳去"的难易程度。D 越大，原子跑得越快；D 越小，原子就像被粘住了一样。
从微观机制来看，D 可以写成：
D = D₀ · exp(-Q/RT)其中：
- D₀ —— 频率因子，跟原子振动频率和跳跃距离有关
- Q —— 扩散激活能，单位 J/mol。这是原子要"翻越"的能量壁垒
- R —— 气体常数，8.314 J/(mol·K)
- T —— 绝对温度，K
我曾经犯过一个低级错误：做高温合金的扩散实验时，忘了考虑温度对 D 的指数影响。结果算出来的扩散深度差了三个数量级。后来我学乖了，每次算扩散之前，先确认温度范围，再查 D₀ 和 Q 的文献值。
工程小技巧：实际工作中，D 值很难直接测量。我建议的做法是：先做一组扩散退火实验，测出浓度分布曲线，然后反推 D。这叫"扩散偶法"，很实用。
典型扩散系数数值参考
给大家一个直观感受。不同体系中 D 的典型值差别很大：
| 扩散体系 | 温度 (℃) | 扩散系数 D (m²/s) | 备注 |
|---|---|---|---|
| Cu 在 Al 中 | 500 | ~10⁻¹⁴ | 常见合金体系 |
| C 在 γ-Fe 中 | 1000 | ~10⁻¹¹ | 渗碳工艺参考 |
| Ni 在 Ni 中（自扩散） | 1000 | ~10⁻¹⁶ | 自扩散很慢 |
| H 在 Pd 中 | 25 | ~10⁻¹⁰ | 氢扩散极快 |
你看，同样是扩散，氢在钯中的扩散系数比镍自扩散快了 10⁶ 倍。这就是为什么氢脆问题那么棘手——氢原子跑得太快了，你根本拦不住。
知识体系结构图
下面这张图，是我自己总结的菲克第一定律知识框架。你把它存下来，以后遇到扩散问题，先对着这张图捋一遍思路：
避坑指南：我曾经在计算扩散通量时，直接把浓度差除以距离当作梯度。这是错的！梯度是微分，不是差分。只有当浓度分布是线性时，才能用 ΔC/Δx 近似。如果曲线很陡，误差会非常大。
好了，菲克第一定律就讲到这里。记住三件事：通量正比于梯度、负号代表方向、D 是温度的函数。下次你看到扩散数据，先问自己三个问题——是稳态吗？梯度算对了吗？D 值查准了吗？
公众号:蓝海资料掘金营，微信deep3321