[论文解读] A Second-Order Energy Stable Backward Differentiation Formula Method for the Epitaxial Thin Film Equation with Slope Selection

Wenqiang Feng, Cheng Wang+2 more|arXiv (Cornell University)|Jun 6, 2017

Fluid Dynamics and Thin Films参考文献 17被引用 4

一句话总结

该论文提出了一种用于具有斜率选择的外延薄膜方程的二阶能量稳定后向差分公式（BDF）有限差分格式，通过引入稳定化的常系数项以确保无条件能量稳定性。该方法实现了最优收敛性，并采用高效的预条件求解器处理非线性四拉普拉斯系统，在数值模拟中表现出稳健的长时间稳定性与准确性。

ABSTRACT

In this paper, we study a novel second-order energy stable Backward Differentiation Formula (BDF) finite difference scheme for the epitaxial thin film equation with slope selection (SS). One major challenge for the higher oder in time temporal discretization is how to ensure an unconditional energy stability and an efficient numerical implementation. We propose a general framework for designing the higher order in time numerical scheme with unconditional energy stability by using the BDF method with constant coefficient stabilized terms. Based on the unconditional energy stability property, we derive an $L^\infty_h (0,T; H_{h}^2)$ stability for the numerical solution and provide an optimal the convergence analysis. To deal with the 4-Laplacian solver in an $L^{2}$ gradient flow at each time step, we apply an efficient preconditioned steepest descent algorithm and preconditioned nonlinear conjugate gradient algorithm to solve the corresponding nonlinear system. Various numerical simulations are present to demonstrate the stability and efficiency of the proposed schemes and slovers.

研究动机与目标

• 为具有斜率选择的外延薄膜方程开发一种二阶时间离散化方法，以确保无条件能量稳定性。
• 通过在BDF框架中引入稳定化的常系数项，克服在高阶时间格式中保持能量稳定性的挑战。
• 提供严格的收敛性分析，并建立数值解在 $L^{ty}_{h}(0,T;H_{h}^{2})$ 范数下的稳定性。
• 设计高效的迭代求解器，用于求解每时间步产生的非线性四拉普拉斯系统，采用预条件最速下降法与非线性共轭梯度法。
• 通过使用谱方法与快速傅里叶变换的大量数值模拟，验证该格式在长时间尺度下的稳定性与准确性。

提出的方法

• 构建一种带有常系数稳定项的二阶BDF格式，以确保具有斜率选择的外延薄膜方程的无条件能量稳定性。
• 应用一种新颖的能量稳定性框架，利用BDF方法结合稳定化手段处理非线性项，无需显式时间分裂或外推。
• 采用有限差分空间离散化方法，结合周期性边界条件，并使用傅里叶谱方法以实现空间上的高阶精度。
• 在每时间步中，实现预条件最速下降算法与预条件非线性共轭梯度法，以求解由四拉普拉斯项产生的非线性系统。
• 利用快速傅里叶变换（FFT）加速迭代求解器中的线性系统求解，提升计算效率。
• 基于格式的无条件能量稳定性和一致性，推导出在 $L^{ty}_{h}(0,T;H_{h}^{2})$ 范数下的最优收敛速率。

实验结果

研究问题

• RQ1 能否为具有斜率选择的外延薄膜方程构造一种二阶BDF格式，以保证无条件能量稳定性？
• RQ2 如何将稳定化的常系数项融入BDF框架中，以在不依赖显式时间分裂或大稳定化参数的情况下维持能量稳定性？
• RQ3 在能量稳定性条件下，所提格式在 $L^{ty}_{h}(0,T;H_{h}^{2})$ 范数下的最优收敛速率是多少？
• RQ4 预条件迭代求解器在每时间步处理非线性四拉普拉斯系统时的效率如何？
• RQ5 所提格式是否能维持长时间稳定性与准确性，以满足薄膜外延中慢动力学模拟的需求？

主要发现

• 所提出的二阶BDF格式通过新颖的稳定化框架与常系数项，实现了无条件能量稳定性。
• 在 $L^{ty}_{h}(0,T;H_{h}^{2})$ 范数下证明了最优收敛速率，证实了该格式的准确性和鲁棒性。
• 该格式无条件稳定，并在长时间模拟中保持能量衰减，与物理上预期的 $t^{-1/3}$ 能量衰减行为一致。
• 预条件最速下降法与非线性共轭梯度求解器显著提升了每时间步求解四拉普拉斯系统时的效率。
• 数值模拟结果证实了该格式的稳定性、准确性和高效性，即使在长时间模拟下也表现出极小的性能退化。
• 与现有二阶格式相比，该方法在能量稳定性和长时间行为方面表现更优，尤其避免了改进型Crank-Nicolson方法中观察到的性能劣化问题。