2. 扩散理论基础

扩散是材料中原子、离子或分子在浓度梯度、化学势梯度或外场作用下发生定向迁移的过程。理解扩散，不能只停留在“高浓度向低浓度跑”的直觉层面，而需要建立从通量定义、驱动势、守恒关系到适用边界的完整链条。本节围绕 Fick 定律、扩散系数、化学势驱动和质量守恒方程展开，并说明它们在材料工程中的适用条件。

一、通量与 Fick 第一定律

扩散通量 J 定义为单位时间通过单位面积的物质量，常用单位为 mol·m⁻²·s⁻¹ 或 kg·m⁻²·s⁻¹。Fick 第一定律描述稳态扩散：

J = -D · (∂C/∂x)

其中 D 为扩散系数（m²·s⁻¹），C 为浓度（mol·m⁻³），x 为空间坐标。负号表示扩散方向与浓度梯度方向相反，即物质从高浓度区向低浓度区迁移。这里约定：坐标 x 向右为正，浓度梯度 ∂C/∂x 为正表示浓度沿 +x 方向升高，则通量 J 为负，表示物质沿 −x 方向流动。全文及图示均采用这一约定。

Fick 第一定律适用于稳态或准稳态扩散，即浓度分布不随时间显著变化。在材料工程中，渗碳初期、薄膜扩散偶的短时扩散、以及某些稳态渗透实验，都可以用第一定律近似描述。

二、Fick 第二定律与质量守恒

当浓度随时间变化时，需要引入质量守恒方程。在一维情况下，取微元体，根据“进入通量 − 流出通量 = 内部浓度变化”可得：

∂C/∂t = -∂J/∂x

将 Fick 第一定律代入，若 D 为常数，则得到 Fick 第二定律：

∂C/∂t = D · (∂²C/∂x²)

该方程说明：某点浓度随时间的变化率，正比于该点浓度分布曲线的曲率。曲率越大，浓度变化越快。若 D 随浓度或位置变化，则更一般的形式为：

∂C/∂t = ∂/∂x [ D(C) · ∂C/∂x ]

在三维情况下，质量守恒方程可写为 ∂C/∂t = −∇·J，结合 J = −D∇C 得到 ∂C/∂t = ∇·(D∇C)。这是扩散问题数值求解和解析求解的共同出发点。

三、扩散系数及其物理含义

扩散系数 D 是扩散动力学的核心参数，通常遵循 Arrhenius 关系：

D = D₀ · exp(-Q / (R·T))

其中 D₀ 为频率因子（m²·s⁻¹），Q 为扩散激活能（J·mol⁻¹），R 为气体常数 8.314 J·mol⁻¹·K⁻¹，T 为绝对温度（K）。激活能越高，扩散对温度越敏感。例如，碳在 γ-Fe 中的扩散激活能约为 135～145 kJ·mol⁻¹，而在 α-Fe 中约为 80 kJ·mol⁻¹，因此碳在铁素体中的扩散比在奥氏体中快得多。

扩散系数还受晶体结构、缺陷、合金元素和晶界影响。晶界扩散系数通常比晶内扩散高几个数量级，因此在低温或细晶材料中，晶界扩散可能成为主导机制。实际工程中，若只给出单一 D 值，往往隐含了“以晶内扩散为主”的假设。

四、化学势驱动：从浓度梯度到势梯度

严格来说，扩散的真正驱动力不是浓度梯度，而是化学势梯度。化学势 μ 可写为：

μ = μ₀ + R·T·ln(a)

其中 a 为活度。对于理想溶液，a = C，化学势梯度与浓度梯度一致；对于非理想体系，活度系数偏离 1，扩散方向可能由活度梯度决定，甚至出现“上坡扩散”——原子从低浓度区向高浓度区迁移。例如，在调幅分解中，溶质原子在化学势驱动下发生上坡扩散，形成周期性成分起伏。

更一般的通量表达式为：

J = -M · C · ∇μ

其中 M 为迁移率。当 μ 与 C 的关系已知时，可将其转化为 Fick 形式，得到互扩散系数或本征扩散系数。因此，Fick 定律是化学势驱动在理想或近理想条件下的简化形式。

五、适用条件与工程判断

Fick 定律在材料中的适用需要满足若干条件：

• 连续介质假设：扩散尺度远大于原子间距，浓度可定义为连续函数。
• 局部平衡假设：每个微小区域内部达到热力学平衡，化学势有定义。
• 无对流或对流可忽略：若存在流体流动、应力驱动或电场驱动，需在通量方程中增加相应项。
• 各向同性或已知各向异性：立方晶系可近似各向同性；六方、四方等结构需用扩散张量描述。
• 稀溶液或已知活度系数：高浓度合金中必须考虑活度修正。

在工程计算中，若扩散距离远大于晶粒尺寸、温度较高、时间较长，通常可用 Fick 定律配合有效扩散系数进行估算。若存在晶界快速通道、应力梯度或强化学势耦合，则需采用更复杂的模型。

六、教学示例：渗碳扩散估算

以低碳钢渗碳为例，假设 900 ℃ 下碳在奥氏体中的扩散系数 D ≈ 1.5×10⁻¹¹ m²·s⁻¹，表面碳浓度维持 1.0 wt%，心部初始碳浓度 0.2 wt%。若要求渗层深度（定义为碳浓度达到 0.4 wt% 的位置）约为 0.8 mm，可用误差函数解估算时间：

C(x,t) = C_s - (C_s - C_0)·erf( x / (2√(D·t)) )

代入 C_s = 1.0，C_0 = 0.2，C = 0.4，x = 0.8 mm，可反解出 erf 参数约为 0.75，对应自变量约 0.81，进而得到 t ≈ (x²)/(4·D·0.81²) ≈ 1.6×10⁴ s，约 4.5 小时。该结果属于示例估算，实际工艺还需考虑炉气碳势、表面反应、晶界扩散和温度均匀性。

综上，Fick 定律给出了扩散通量与浓度梯度的定量关系，质量守恒方程将其推广到非稳态；扩散系数决定动力学快慢，化学势揭示真正驱动来源。在材料工程中，只有同时检查连续介质、局部平衡、无外场和各向同性等条件，才能合理使用这些方程，并在必要时引入晶界扩散、活度修正和耦合场模型。