2. 晶界热力学与迁移驱动力
晶界迁移是再结晶、晶粒长大与相变组织演化的核心动力学过程。要把它写进相场模型，第一步不是直接写演化方程，而是先回答一个更基础的问题：晶界为什么会动、朝哪个方向动、动得多快。这取决于晶界的热力学状态与作用在其上的驱动力。本节从晶界能出发，逐层建立曲率驱动、储存能驱动、化学自由能差三类驱动力，并说明温度场如何通过迁移率耦合进速度表达式。
一、晶界能：迁移的“价格标签”
晶界是相邻晶粒之间的二维缺陷面，原子排列失配使其具有超额自由能，即晶界能 γ（单位 J/m²）。γ 的大小取决于两侧晶粒的取向差、晶界面法向与晶界类型：小角度晶界由位错阵列构成，γ 随取向差增大而升高；大角度晶界趋于饱和，典型金属约 0.3～1.0 J/m²；共格孪晶界能量极低，可低至普通大角度晶界的十分之一。晶界能不仅决定界面形状，还通过界面刚度（γ 加二阶导数项）影响曲率驱动的强弱。
在相场模型中，晶界能通常以梯度能系数 κ 与势垒高度 w 的形式进入自由能泛函：γ ∝ √(κ·w)。这意味着调 γ 实际上是同时调 κ 与 w，只改一个会连带改变界面宽度与迁移率，这是建模时最容易踩的坑。
二、四类驱动力与迁移通量
晶界迁移的广义驱动力是界面两侧自由能密度差，可写成压力量纲（Pa）。常见来源有四类：
- 曲率驱动：弯曲界面存在 Laplace 压力差 ΔP = γ·κ_c（κ_c 为曲率，凸向一侧为正）。它使晶界向曲率中心移动，是晶粒长大与晶界网络粗化的主因。
- 储存能驱动：冷变形后位错密度 ρ 带来储存能 E_s ≈ α·μ·b²·ρ（α≈0.5，μ 剪切模量，b 柏氏矢量）。再结晶时，变形基体与新生无畸变晶粒之间的储存能差可达 10⁶～10⁷ J/m³，是再结晶的主导驱动力。
- 化学自由能差：两相或两成分间摩尔自由能差 ΔG_m 提供驱动力 P_chem = ΔG_m / V_m（V_m 摩尔体积），主导相变与析出伴随的界面迁移。
- 温度场耦合：温度不直接产生力，但通过迁移率 M(T) 与扩散系数强烈调制速度，并可通过热梯度引入额外界面力。
把上述驱动力统一记为总驱动力 P = P_curv + P_store + P_chem + P_ext，则晶界迁移速度与通量可写成：
v = M(T) · P
J = -M(T) · ∇(δF/δφ) （相场通量形式）
M(T) = M0 · exp(-Q_m / (R·T)) （Arrhenius 迁移率）
P_curv = γ · κ_c
P_store = α · μ · b² · (ρ_def - ρ_rec)
P_chem = ΔG_m / V_m其中 M0 为指前因子，Q_m 为迁移激活能，R 为气体常数。速度 v 与驱动力 P 在低驱动力下近似线性（迁移率控制区）；当驱动力极高时，界面可能进入扩散控制或溶质拖曳区，线性关系失效，此时需引入拖曳因子或溶质拖曳模型。
三、驱动力与温度场的耦合关系
温度对迁移的影响是双重的：一方面 M(T) 随温度指数上升，升温显著加快晶界迁移；另一方面，储存能与化学自由能差本身也随温度变化（如 ΔG_m 随过冷度增大而增大）。因此实际速度是“驱动力随温度变化”与“迁移率随温度变化”的乘积。在焊接热影响区、增材制造熔池等存在陡峭温度梯度的场景中，同一晶界在不同位置经历不同 T，迁移速度差异可达数个数量级，这正是热影响区晶粒异常长大的根源。
下面这张图把四类驱动力、迁移率与速度表达式的耦合关系集中呈现，便于对照建模时各参数的输入位置。
四、建模时的参数取值与量级参考
把表达式落到相场代码里，需要给出可用的量级。以下为一组教学示例参数（纯铜类金属，仅用于说明量级，非实测值）：
| 参数 | 符号 | 示例取值 | 说明 |
|---|---|---|---|
| 晶界能 | γ | 0.5 J/m² | 大角度晶界典型值 |
| 储存能差 | ΔE_s | 5×10⁶ J/m³ | 中等变形量 |
| 化学自由能差 | ΔG_m | 10³ J/mol | 小过冷度 |
| 摩尔体积 | V_m | 7×10⁻⁶ m³/mol | 铜 |
| 迁移激活能 | Q_m | 1.2×10⁵ J/mol | 晶界扩散控制 |
| 指前因子 | M0 | 10⁻² m⁴/(J·s) | 需标定 |
以 T=700 K 估算：M(T)=10⁻²·exp(-1.2×10⁵/(8.314×700))≈10⁻²·exp(-20.6)≈1.1×10⁻¹¹ m⁴/(J·s)。若仅由储存能驱动，P≈5×10⁶ Pa，则 v≈5.5×10⁻⁵ m/s，即约 55 μm/s，与再结晶初期晶界迁移的合理量级相符。若同时叠加曲率驱动（κ_c≈10⁶ m⁻¹，P_curv≈5×10⁵ Pa），总驱动力提高约 10%，速度相应增加。这说明再结晶早期储存能主导，后期曲率驱动逐渐接管。
五、从表达式到可执行判断
实际建模时可按以下顺序检查：先确认 γ、ΔE_s、ΔG_m 的量级是否与合金体系匹配；再确认 M(T) 的激活能与指前因子是否有文献或标定依据；然后检查温度场是否随时间更新，避免用恒温近似处理焊接或增材过程；最后用单界面算例验证 v 是否落在预期区间。若速度明显偏大，优先排查 M0 是否取值过高或界面宽度过窄；若速度几乎不随温度变化，则检查 M(T) 是否真的接入了温度场，而不是被固定为常数。
掌握这四类驱动力与迁移率的关系后，下一层工作才是把它们嵌入相场演化方程，让多晶组织在储存能与曲率共同作用下自然粗化与再结晶。