[论文解读] An Efficient Finite Difference-based Implicit Solver for Phase-Field Equations with Spatially and Temporally Varying Parameters

Zirui Mao, G. R. Liu+1 more|arXiv (Cornell University)|Jan 22, 2024

Solidification and crystal growth phenomenaMaterials Science被引用 3

一句话总结

本论文提出了一种基于有限差分法的隐式求解器，用于求解具有空间和时间变化参数的相场方程，利用Eyre定理确保无条件稳定性。该方法在非均匀且随时间变化的条件下模拟微观结构演化时表现出高精度和高效率，尽管目前仍需使用均匀网格和特定边界条件。

ABSTRACT

The phase field method is an effective tool for modeling microstructure evolution in materials. Many efficient implicit numerical solvers have been proposed for phase field simulations under uniform and time-invariant model parameters. We use Eyre's theorem to develop an unconditionally stable implicit solver for spatially non-uniform and time-varying model parameters. The accuracy, unconditional stability, and efficiency of the solver is validated against benchmarking examples. In its current form, the solver requires a uniform mesh and may only be applied to problems with periodic, Neumann, or mixed periodic and Neumann boundary conditions.

研究动机与目标

• 开发一种稳定且高效的相场模拟数值求解器，适用于非均匀且随时间变化的模型参数。
• 将现有隐式求解器扩展至相场建模中通常使用的均匀且时间不变参数之外的范围。
• 利用Eyre定理确保在空间和时间参数变化情况下仍具有无条件稳定性。
• 通过与既有的基准问题对比，验证求解器的精度和效率。
• 识别当前方法在网格类型和边界条件要求方面的局限性，以利于实际应用。

提出的方法

• 求解器基于具有空间和时间变化参数的相场方程的有限差分离散化。
• 应用Eyre定理构建凸分解格式，以保证无条件能量稳定性。
• 通过牛顿-拉夫森法迭代求解所得非线性系统，以实现稳健收敛。
• 该方法假设使用均匀计算网格，以维持稳定性和效率。
• 边界条件被限制为周期性、诺伊曼型或周期性-诺伊曼混合型。
• 求解器以半隐式方式实现，以在精度和计算成本之间取得平衡。

实验结果

研究问题

• RQ1 当模型参数在空间和时间上变化时，隐式有限差分求解器是否能保持无条件稳定性？
• RQ2 在非均匀参数条件下，所提出的求解器在精度和效率方面与现有求解器相比如何？
• RQ3 所提出方法在网格类型和边界条件适用性方面存在哪些限制？
• RQ4 在复杂的相场模拟中，使用Eyre定理在多大程度上能确保能量稳定性？
• RQ5 该求解器在具有变化参数的标准基准问题上的表现如何？

主要发现

• 所提出的求解器在具有空间和时间变化参数的相场方程中表现出无条件稳定性，这由Eyre定理保证。
• 在基准模拟中，求解器实现了高精度，与已知参考解非常接近。
• 通过采用半隐式、基于牛顿法的迭代求解策略，保持了计算效率。
• 该方法受限于均匀网格和特定边界条件，限制了其在复杂几何结构中的直接应用。
• 与标准基准的对比验证了求解器在微观结构演化模拟中的可靠性。
• 该方法能够实现非均匀且随时间变化参数条件下的稳定模拟，而这类条件在真实材料加工过程中十分常见。