[论文解读] Extending phase field models of grain boundaries to three dimensions

Ryosuke Kobayashi, James A. Warren|arXiv (Cornell University)|Feb 16, 2005

Solidification and crystal growth phenomena被引用 4

一句话总结

本文通过用 3D 旋转矩阵 P ∈ SO(3) 替代标量取向场 θ，将二维晶界相场模型扩展至三维，实现了多晶微结构演化的逼真模拟。作者采用测地线或基于矩阵范数的取向梯度度量方法构建自由能泛函，推导出一致的演化方程，并在具有完整 3D 取向动力学的薄膜状系统中展示了三维晶粒生长、碰撞与粗化过程。

ABSTRACT

In this letter we describe a method of extending an existing phase field model of polycrystalline solidification from two to three dimensions (3D).

研究动机与目标

• 将现有的二维晶界相场模型扩展至三维，实现具有完整取向动力学的三维多晶微结构模拟。
• 解决在三维晶体学取向表示中的数学挑战，将二维角度 θ 替换为 3×3 旋转矩阵 P ∈ SO(3)。
• 在三维中定义物理上一致的取向梯度度量方法，采用矩阵距离或 SO(3) 上的测地线长度。
• 推导出相场 φ 和取向场 P 的演化方程，确保时间演化过程中保持 SO(3) 约束。
• 展示该模型在模拟薄膜状系统中三维晶粒生长、粗化和碰撞方面的能力。

提出的方法

• 使用 SO(3) 的元素表示三维晶体取向，即行列式为 1 的 3×3 正交矩阵，嵌入于 R^9 中。
• 采用两种候选度量方法定义取向梯度：I 型（R^9 中的 ||P - Q||）和 II 型（通过 PQ⁻¹ 的迹计算的测地线距离）。
• 通过在二维模型中以 |∇P| 替代 |∇θ|，构建自由能泛函，保留相变和晶界能相关项。
• 通过对自由能泛函关于 φ 和 P 求变分导数，推导出两组场的演化方程。
• 采用 SO(3) 坐标上的约束动力学（如 Rodrigues 向量）或 R^9 中的投影方法实现演化方程，以保持 SO(3) 成员身份。
• 在时间积分过程中使用数值投影方法强制执行 SO(3) 约束，确保取向演化具有物理一致性。

实验结果

研究问题

• RQ1 如何在保持晶界能和取向差效应物理真实性的同时，将二维晶界相场模型推广至三维？
• RQ2 二维中 |∇θ| 的适当三维类比是什么？哪种几何度量最能准确捕捉三维中的取向差？
• RQ3 如何在时间积分过程中数值演化取向场 P ∈ SO(3) 的动力学，同时在整个时间演化过程中保持 SO(3) 约束？
• RQ4 三维取向动力学对多晶微结构演化（包括晶粒生长和粗化）有何影响？
• RQ5 不同取向梯度度量方法（I 型与 II 型）在三维模拟结果中产生何种影响？

主要发现

• 该模型通过用 SO(3) 旋转矩阵替代标量取向，成功将二维晶界相场模型推广至三维，实现了完整的三维微结构模拟。
• II 型测地线度量方法在三维中对取向差的表示更具物理准确性，尤其在不连续区域表现更优。
• 变分公式推导出的 φ 和 P 演化方程具有一致性，通过投影或参数化方法正确施加了约束。
• 模拟结果展示了薄膜系统中逼真的三维晶粒生长、碰撞与粗化过程，彩色编码的欧拉角可视化了三维取向演化。
• 该模型在数值上稳定，能够模拟复杂的微结构演化，未来工作可进一步引入各向异性和晶体对称性。