[论文解读] Numerical Energy Dissipation for Time-Fractional Phase-Field Equations

Chaoyu Quan, Tao Tang+1 more|arXiv (Cornell University)|Sep 14, 2020

Solidification and crystal growth phenomena参考文献 31被引用 11

一句话总结

本文利用一阶稳定半隐式 L1 格式，为时间分数阶相场方程建立了离散能量律，证明了在非均匀时间步长下能量耗散与有界性。该分析进一步扩展至 (2−α)-阶 L1-SAV 格式，表明即使在 α 较小时，修正能量仍保持有界且耗散，且经严格数值验证。

ABSTRACT

The numerical integration of phase-field equations is a delicate task which needs to recover at the discrete level intrinsic properties of the solution such as energy dissipation and maximum principle. Although the theory of energy dissipation for classical phase field models is well established, the corresponding theory for time-fractional phase-field models is still incomplete. In this article, we study certain nonlocal-in-time energies using the first-order stabilized semi-implicit L1 scheme. In particular, we will establish a discrete fractional energy law and a discrete weighted energy law. The extension for a $(2-α)$-order L1 scalar auxiliary variable scheme will be investigated. Moreover, we demonstrate that the energy bound is preserved for the L1 schemes with nonuniform time steps. Several numerical experiments are carried to verify our theoretical analysis.

研究动机与目标

• 为时间分数阶相场方程建立离散分数阶与加权能量律。
• 分析一阶稳定半隐式 L1 格式在非均匀时间步长下的能量稳定性。
• 将分析扩展至更高阶的 (2−α)-阶 L1-SAV 格式。
• 验证能量有界性与耗散性在数值上是否得以保持，尤其是当 α 较小时。
• 为适用于时间分数阶偏微分方程的能量稳定数值格式提供一个框架。

提出的方法

• 使用一阶稳定半隐式 L1 格式来近似相场方程中的 Caputo 时间分数阶导数。
• 应用一种特殊的 Cholesky 分解技术，以分析能量稳定性并推导离散能量律。
• 为 L1 格式推导出离散分数阶能量律与离散加权能量律。
• 将分析扩展至 (2−α)-阶 L1-SAV 格式，确保能量有界性与耗散性。
• 采用非均匀时间网格（ graded ）以提高在 t=0 附近的分辨率，尤其适用于 α 较小时。
• 通过 Allen–Cahn、Cahn–Hilliard 与 MBE 方程的数值实验验证理论结果。

实验结果

研究问题

• RQ1 能否为应用于时间分数阶相场方程的一阶 L1 格式建立离散分数阶能量律？
• RQ2 在非均匀时间步长下，L1 格式是否能保持能量有界性与耗散性？
• RQ3 能否使 (2−α)-阶 L1-SAV 格式在时间分数阶相场模型中保持能量稳定性与有界性？
• RQ4 当 α 减小时，L1 格式中的能量耗散与数值分辨率如何变化？
• RQ5 Cholesky 分解技术在多大程度上能够实现时间分数阶问题的能量稳定性分析？

主要发现

• 为一阶 L1 格式建立了离散分数阶能量律，确保在半离散层面实现能量耗散。
• 在非均匀时间步长下，L1 格式保持了能量有界性，即使 α 较小，经数值验证成立。
• 当 α=0.3 时，标准 L1 格式因分辨率不足导致能量耗散被破坏，但能量有界性仍保持。
• (2−α)-阶 L1-SAV 格式确保了能量有界性，并满足离散分数阶能量律。
• 数值实验表明，修正能量在所有测试的 α 值与网格类型下，始终有界于初始能量之上。
• Cholesky 分解技术在时间分数阶相场问题的能量稳定性分析中被证明是有效的。