4. 相场序参量设计

相场模拟晶界迁移与再结晶，本质上是在一个连续场框架内同时描述“晶粒长什么样”“哪里已经再结晶”“驱动力从哪来”“热激活有多强”四类信息。这四类信息分别由多晶取向场、再结晶分数场、位错密度场和温度场承载。序参量设计一旦不合理，后续演化方程再精细也难以得到物理可信的结果。本节围绕这四个场的选取原则、物理约束和相互耦合关系展开，并给出可操作的参数化思路。

一、多晶取向场：用连续变量区分晶粒

多晶取向场用于区分不同晶粒及其取向差。常见做法有两种：一是为每个晶粒分配一个局部序参量 η_i（i=1…N），晶粒内部 η_i≈1、其余≈0，晶界处多个 η 同时取中间值；二是采用连续取向角 θ(r) 或四元数场，直接表示局部晶格取向。前者适合晶粒数较少、拓扑变化剧烈的再结晶场景，后者适合描述取向梯度与织构演化。

选取原则：序参量必须满足对称性与归一化约束。若采用多序参量方案，应保证 Ση_i=1 或 Ση_i²=1，避免出现“无主晶粒”的真空区。晶界能 γ 通常写成 η 梯度项与势垒项的组合，例如 γ ∝ ∫[κ|∇η|² + f(η)]dx，其中势垒高度决定界面宽度。界面宽度 l 与梯度系数 κ、势垒高度 w 满足 l ∝ √(κ/w)，因此 κ 和 w 不能独立随意取值，必须与网格尺寸 Δx 匹配：一般要求 l ≥ 4Δx，否则界面无法被网格分辨，演化会出现虚假钉扎。

取向差角 Δθ 与晶界能、迁移率的关系需要显式引入。低角度晶界（Δθ < 15°）能量和迁移率都显著低于高角度晶界，可用 Read-Shockley 型关系 γ(Δθ)=γ_m(Δθ/θ_m)(1−ln(Δθ/θ_m)) 描述。若取向场只记录“晶粒编号”而不记录取向差，就无法区分小角与大角晶界，再结晶形核与长大行为会失真。

二、再结晶分数场：区分变形基体与新晶粒

再结晶分数场 X(r,t) 用于标记材料是否已完成再结晶，通常取 0 表示变形基体、1 表示再结晶晶粒。它不是一个独立演化的序参量，而是由取向场和位错密度场共同决定的“状态标签”。工程上常写成 X = 1 − exp(−k tⁿ) 的 JMAK 形式，但在相场中更推荐用局部判据：当某区域位错密度低于临界值 ρ_c 且取向场已切换到新晶粒时，X 由 0 向 1 平滑过渡。

物理约束：X 必须单调、有界，且与取向场保持几何一致——不能出现“X=1 但取向仍为变形基体”的矛盾。为避免数值振荡，X 的演化宜采用弛豫形式 ∂X/∂t = −L_X(∂f/∂X)，其中 f 是双阱势，保证 X 在 0 和 1 附近稳定。再结晶形核位置通常由位错密度场的高值区、晶界三叉节点或第二相粒子界面提供，因此 X 的初始分布应结合这些位置设定，而不是均匀随机撒点。

三、位错密度场：提供再结晶驱动力

位错密度场 ρ(r,t) 是再结晶的“燃料”。储存能 E_s ≈ αGb²ρ（α≈0.5，G 为剪切模量，b 为伯氏矢量），它构成再结晶的驱动力 P_D = E_s。相场中通常把 ρ 作为独立标量场演化，或通过储存能密度间接耦合。

选取原则：ρ 的初始值应反映变形工艺。例如冷轧 50% 的铜，位错密度约 10¹⁵–10¹⁶ m⁻²；退火再结晶时 ρ 在再结晶区迅速降至 10¹² m⁻² 量级。演化方程可写为 ∂ρ/∂t = −k_rρ + k_g|∇v|，其中第一项为回复、第二项为变形增殖。关键约束是：ρ 不能为负，且再结晶区 ρ 的下降速率应与 X 的上升速率匹配，否则会出现“再结晶已完成但储存能仍很高”的非物理状态。

四、温度场：控制热激活与迁移率

温度场 T(r,t) 通过 Arrhenius 关系控制晶界迁移率 M = M₀exp(−Q/RT) 和形核率。选取原则：若模拟等温退火，T 可取常数；若模拟焊接、激光退火等存在温度梯度的过程，T 必须作为空间场参与演化，并满足热传导方程 ρ_mc_p∂T/∂t = ∇·(k∇T) + q。此时 T 与再结晶场双向耦合：再结晶释放潜热，温度又反过来加速迁移。

物理约束：T 必须高于绝对零度，且 M 的指数项不能因 T 过小导致数值下溢。工程上常把 Q/RT 控制在 10–40 范围内，对应 M 变化 4–17 个数量级，既保证物理真实性又避免数值刚性。

五、四场耦合关系与参数一致性检查

四个场不是并列关系，而是有明确的驱动链：温度场决定迁移率和形核率；位错密度场提供储存能驱动力；取向场决定晶界位置与取向差；再结晶分数场标记状态并反馈到位错密度演化。设计时必须做一致性检查：界面宽度与网格匹配、驱动力与储存能量级匹配、迁移率与温度匹配、再结晶分数与位错密度下降同步。

六、参数选取的实操建议

以纯铜等温退火为例，给出可参考的示例参数（属于合理假设，非实测值）：网格 Δx = 0.1 μm，界面宽度 l = 0.5 μm，梯度系数 κ 与势垒 w 按 l ∝ √(κ/w) 反推；初始位错密度 ρ₀ = 5×10¹⁵ m⁻²，对应储存能约 1.5×10⁶ J/m³；退火温度 T = 673 K，取 Q = 120 kJ/mol，则 M 约为 10⁻¹⁴ m⁴/(J·s) 量级；再结晶临界位错密度 ρ_c = 10¹⁴ m⁻²。时间步长 Δt 需满足 Δt ≤ Δx²/(4Mκ) 的稳定性条件，通常取 0.01–0.1 s。

综上，序参量设计的核心不是“选几个场”，而是让每个场承担明确的物理角色，并在界面宽度、驱动力、迁移率和状态切换四个层面保持量级一致。只有四场耦合自洽，后续的晶界迁移与再结晶演化才具有可解释性和可复现性。