3. 相场模型的热力学基础

相场模型之所以能描述晶粒长大，根本原因在于它把“系统往哪儿演化”这件事交给了一个标量泛函——总自由能。只要写出总自由能关于序参量场的表达式，再让系统沿自由能下降方向演化，晶粒的粗化、界面迁移、曲率驱动等行为就会自然涌现。因此，理解相场模型的第一步，不是写动力学方程，而是把自由能泛函的每一项“是什么、为什么需要、参数如何影响结果”讲清楚。

一、总自由能泛函：把能量写成场的函数

在单相多晶系统中，常用一组连续场变量 \(\eta_1(\mathbf{r}),\eta_2(\mathbf{r}),\dots,\eta_p(\mathbf{r})\) 表示不同取向的晶粒。每个 \(\eta_i\) 在属于第 \(i\) 个晶粒的区域接近 1，在其它区域接近 0，界面处则在 0 与 1 之间连续过渡。系统的总自由能泛函通常写成：

F = ∫ [ f(η_1,…,η_p) + Σ_i (κ_i/2)|∇η_i|² ] dV

其中第一项 \(f\) 是体自由能密度，第二项是梯度能密度，\(\kappa_i\) 为梯度能系数。这个积分把局部能量密度在空间上累加，得到整个系统的总能量。相场演化的驱动力来自 \(F\) 对 \(\eta_i\) 的变分导数，即 \(\delta F/\delta\eta_i\)。因此，泛函中每一项的形式和参数，都会直接改变晶粒形貌、界面宽度和长大速率。

二、体自由能：决定“稳定相”与“界面势垒”

体自由能密度 \(f\) 只与局部序参量有关，不含梯度。对多晶系统，常用双阱势或四阶多项式构造：

f = Σ_i [ a η_i² + b η_i⁴ ] + c Σ_{i≠j} η_i² η_j²

当 \(a<0\)、\(b>0\) 时，单个 \(\eta_i\) 在 0 和 1 附近各有一个极小值，分别对应“不属于该晶粒”和“属于该晶粒”。交叉项 \(c η_i² η_j²\) 则惩罚两个晶粒同时占据同一位置，保证界面处不同 \(\eta_i\) 不会同时接近 1。体自由能决定了界面两侧的平衡值，也决定了界面能中“势垒”部分的大小。若 \(a\)、\(b\) 取值不当，可能出现界面过宽、界面能异常或数值刚性过强等问题。

三、梯度能与界面能：界面不是数学突变，而是有限宽度

梯度能 \(\frac{\kappa}{2}|\nabla\eta_i|^2\) 惩罚序参量的空间变化。若没有梯度能，界面会退化为无限薄的数学面，无法用连续场描述。梯度能系数 \(\kappa\) 越大，界面越宽，界面能也越高。在简单一维双阱情形下，界面能可近似写成：

γ ≈ ∫ [ Δf(η) + (κ/2)(dη/dx)² ] dx

其中 \(\Delta f\) 是体自由能相对平衡态的增量。界面能由“体自由能势垒”和“梯度能”共同贡献，二者比例决定界面宽度：\(\kappa\) 越大或势垒越低，界面越宽。工程上常通过目标界面能 \(\gamma\) 和界面宽度 \(w\) 反推 \(\kappa\) 与势垒高度，而不是随意给参数。例如，若希望 \(w\approx 1\) nm、\(\gamma\approx 0.5\) J/m²，就需要让 \(\kappa\) 与势垒高度匹配到相应量级。

四、多晶各向异性界面能的引入方式

真实晶粒的界面能依赖晶界取向差和界面法向。若所有界面能相同，晶粒长大只会得到近似等轴的统计形貌，无法再现某些晶粒的择优长大、刻面或织构演化。相场中引入各向异性，常见有三种做法：

• 梯度能系数各向异性：令 \(\kappa\) 依赖界面法向 \(\mathbf{n}=\nabla\eta_i/|\nabla\eta_i|\)，即 \(\kappa(\mathbf{n})\)。这样不同法向的界面具有不同梯度能，从而产生各向异性界面能。
• 高阶梯度项：在泛函中加入 \((\nabla^2\eta_i)^2\) 或 \(\nabla^2\eta_i\) 相关项，使界面能依赖曲率和法向，适合描述强各向异性或刻面。
• 取向差依赖的界面能函数：对多晶系统，令界面能显式依赖相邻晶粒的取向差 \(\Delta\theta\)，例如 \(\gamma(\Delta\theta)=\gamma_0[1-\delta\cos(m\Delta\theta)]\)，再通过相场变量耦合到演化方程中。

第一种方式实现简单、数值稳定，适合大多数晶粒长大模拟；第二种能刻画刻面，但需要更高阶离散和更小时间步；第三种物理意义最直接，但需要额外追踪晶粒取向或取向差场。选择哪种方式，取决于你关心的是统计粗化、织构演化，还是界面形貌细节。

五、参数选择与结果判读：从泛函到可观察形貌

把泛函写出来只是第一步，真正影响模拟可信度的是参数标定。一个实用的做法是：先确定目标界面能 \(\gamma\) 和界面宽度 \(w\)，再反推体自由能势垒高度与梯度能系数。以一维平衡界面为例，若取双阱势 \(f = W\eta^2(1-\eta)^2\)，则界面能近似为 \(\gamma \propto \sqrt{\kappa W}\)，界面宽度 \(w \propto \sqrt{\kappa/W}\)。因此，增大 \(\kappa\) 会同时增大 \(\gamma\) 和 \(w\)，而增大 \(W\) 会增大 \(\gamma\)、减小 \(w\)。工程上通常先固定 \(w\) 为 3～5 个网格间距，再调整 \(W\) 使 \(\gamma\) 匹配目标值。

各向异性引入后，结果判读要关注三点：第一，界面是否出现择优方向，例如某些晶界迁移更快；第二，晶粒形貌是否从等轴变为拉长或刻面；第三，织构是否随演化增强。若模拟中所有晶界仍以相同速率迁移，说明各向异性没有真正耦合到演化方程中，需要检查 \(\kappa(\mathbf{n})\) 是否随法向变化，或 \(\gamma(\Delta\theta)\) 是否被正确代入。若界面出现非物理的网格钉扎或锯齿，则可能是各向异性过强、网格分辨率不足或时间步过大。

总结来说，相场模型的热力学基础可以概括为：体自由能决定“稳定相”和势垒，梯度能决定界面宽度和界面能，各向异性则决定界面能的方向与取向差依赖。三者共同构成总自由能泛函，而泛函的变分导数就是晶粒长大的驱动力。把这三项的参数与目标物理量对应起来，才能让模拟结果既有物理意义，又具备工程参考价值。