3. 扩散的宏观定律

前两节我们从原子跳跃的微观图像出发，建立了扩散的物理直觉：原子在浓度梯度、化学势梯度或应力场驱动下发生定向迁移。但工程中真正要回答的问题往往是宏观的——渗碳层在给定温度和时间下能长到多深？两块异种金属焊接后界面附近的成分如何随时间演化？这些问题的答案不能靠逐个原子追踪，而必须借助连续介质框架下的宏观扩散方程。本节围绕Fick第一定律、Fick第二定律以及扩散方程在典型边界条件下的解展开，重点说明每个方程的物理含义、适用条件以及工程中如何用它们做定量估算。

一、Fick第一定律：通量与浓度梯度的线性响应

Fick第一定律描述的是稳态扩散条件下，扩散通量与浓度梯度之间的定量关系。其数学形式为：

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

其中 J 为扩散通量，单位是 mol/(m²·s) 或 kg/(m²·s)，表示单位时间通过单位面积的物质量；D 为扩散系数，单位 m²/s；∂C/∂x 为沿 x 方向的浓度梯度；负号表示扩散方向与浓度梯度方向相反，即物质从高浓度区向低浓度区净迁移。

这里需要特别强调参照系的约定：本文及后续图示中，x 轴正方向统一取为浓度降低的方向，因此 ∂C/∂x 为负值，负号使 J 为正，表示通量沿 x 正方向。若读者采用相反的坐标约定，负号的位置会变化，但物理结论不变。工程计算中务必先明确坐标方向，再代入符号。

Fick第一定律的适用前提是稳态扩散，即浓度场不随时间变化（∂C/∂t = 0）。严格满足稳态的扩散过程在工程中并不多见，但许多场景可以近似为稳态：例如薄板两侧维持恒定浓度差、扩散偶在长时间保温后中间区域浓度分布趋于线性、气体通过薄壁的渗透等。在这些情况下，第一定律可以直接用于估算通量和扩散速率。

一个典型的教学示例：纯铁在 900°C 渗碳，表面碳浓度维持在 1.0 wt%，心部碳浓度为 0.1 wt%，渗碳层厚度约 1 mm。若取该温度下碳在奥氏体中的扩散系数 D ≈ 1.5×10⁻¹¹ m²/s，并假设浓度沿渗碳层近似线性分布，则浓度梯度约为 (0.1−1.0)/0.001 = −900 wt%/m，通量 J ≈ 1.5×10⁻¹¹ × 900 ≈ 1.35×10⁻⁸ wt%·m/s。这个数值虽然不能直接给出渗层深度，但能帮助判断碳原子向内部的供给速率量级。需要说明的是，该算例中的 D 和浓度分布均为合理假设，用于说明计算方法，不代表实测结果。

二、Fick第二定律：浓度场随时间的演化

实际扩散过程中，浓度场几乎总是随时间变化的。要描述这种非稳态扩散，需要在第一定律的基础上引入质量守恒。考虑一维微元体，进入的通量减去流出的通量等于微元体内物质的积累，结合第一定律消去通量，得到Fick第二定律：

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

当扩散系数 D 不随浓度和位置变化时，方程简化为：

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

Fick第二定律的物理含义非常明确：某一点浓度随时间的变化率，正比于该点浓度分布曲线的曲率。曲率越大，浓度变化越快；曲率为零（浓度线性分布）时，浓度不随时间变化，这正是稳态扩散的特征。因此，第一定律可以看作第二定律在 ∂C/∂t = 0 时的特例。

第二定律是一个偏微分方程，求解需要初始条件和边界条件。初始条件描述扩散开始前浓度随位置的分布，边界条件描述扩散过程中试样表面或端部的浓度约束。不同的初始—边界条件组合对应不同的工程场景，解的形式也各不相同。

三、典型边界条件与扩散方程的解

工程中最常见的边界条件有三类。第一类是恒定表面浓度边界条件：试样表面浓度在扩散过程中保持恒定，例如渗碳时炉内碳势稳定、表面碳浓度维持在某一固定值。第二类是恒定通量边界条件：表面通量保持恒定，例如某些表面渗氮或离子注入过程。第三类是有限源或瞬时源条件：表面沉积一层有限量的扩散物质，随后向内部扩散，例如薄膜扩散偶或表面涂覆后的退火。

对于恒定表面浓度、半无限长试样的经典情形，初始条件为 C(x,0) = C₀（x > 0），边界条件为 C(0,t) = C_s（t > 0），C(∞,t) = C₀。其解为：

(C(x,t) − C₀) / (C_s − C₀) = erfc( x / (2√(Dt)) )

其中 erfc 为余误差函数。这个解是渗碳、渗氮等表面强化工艺计算渗层深度的核心工具。工程上常把达到某一特征浓度（如渗碳层通常取 0.4 wt% 或某一规定值）的位置定义为渗层深度 x_d，则 x_d 与 √(Dt) 成正比。这意味着扩散时间增加四倍，渗层深度才增加一倍，这是扩散工艺中“时间换深度”效率递减的根本原因。

下面这张图把Fick第一定律、第二定律以及恒定表面浓度解之间的关系和典型边界条件组织在一起，帮助读者建立从通量到浓度场演化的整体认识。

四、工程应用中的关键判断

在实际使用这些方程时，有几个容易出错的点需要特别注意。第一，扩散系数 D 对温度极其敏感，通常服从 Arrhenius 关系 D = D₀·exp(−Q/RT)。温度升高几十度，D 可能翻倍甚至更多。因此渗碳工艺中温度的控制精度直接决定渗层深度的重复性。第二，Fick第二定律的解析解大多建立在 D 为常数的假设上。当浓度变化范围很大时（如渗碳从 1.0 wt% 到 0.1 wt%），D 随浓度变化不可忽略，此时需要用数值方法求解变系数扩散方程。第三，半无限长假设要求试样厚度远大于扩散特征长度 √(Dt)。对于薄件或长时间扩散，必须改用有限尺寸的解，否则会高估渗层深度。

判断一个扩散问题该用哪个方程、哪类边界条件，可以遵循以下思路：先看浓度场是否随时间变化，若基本不变则用第一定律估算通量；若明显随时间演化，则用第二定律。再看表面条件：表面浓度是否恒定、是否有外部供给、是否只有有限量的扩散物质。最后看试样尺寸与 √(Dt) 的关系，决定采用半无限解还是有限尺寸解。这套判断逻辑在渗碳、渗氮、扩散焊、涂层退火等工艺中反复出现，是宏观扩散计算的基本功。

本节建立的宏观扩散方程体系，是后续讨论相变动力学、组织演化以及扩散型相变的基础工具。无论是析出相的长大、珠光体的层片间距，还是焊接热影响区的成分偏析，其定量描述最终都要回到 Fick 定律及其在各种边界条件下的解。掌握这些方程的物理含义、适用条件和典型解的形式，是进行材料扩散与相变工程计算的第一步。