[论文解读] Convergence analysis for a stabilized linear semi-implicit numerical scheme for the nonlocal Cahn-Hilliard equation

Xiao Li, Zhonghua Qiao+1 more|arXiv (Cornell University)|Feb 13, 2019

Solidification and crystal growth phenomena参考文献 37被引用 8

一句话总结

本文针对非局部 Cahn–Hilliard 方程，采用傅里叶伪谱空间离散化，对一种稳定化的线性半隐式格式进行了收敛性分析。通过采用离散 $H^{-1}$ 范数进行误差估计，并结合高阶一致性分析与逆不等式，作者建立了时间方向最优一阶收敛与空间方向谱精度，同时在对稳定化参数的改进条件下，证明了均匀 $\ell^\infty$ 有界性与能量稳定性。

ABSTRACT

In this paper, we provide a detailed convergence analysis for a first order stabilized linear semi-implicit numerical scheme for the nonlocal Cahn-Hilliard equation, which follows from consistency and stability estimates for the numerical error function. Due to the complicated form of the nonlinear term, we adopt the discrete $H^{-1}$ norm for the error function to establish the convergence result. In addition, the energy stability obtained in [Du et al., J. Comput. Phys., 363:39--54, 2018] requires an assumption on the uniform $\ell^\infty$ bound of the numerical solution and such a bound is figured out in this paper by conducting the higher order consistency analysis. Taking the view that the numerical solution is indeed the exact solution with a perturbation, the error function is $\ell^\infty$ bounded uniformly under a loose constraint of the time step size, which then leads to the uniform maximum-norm bound of the numerical solution.

研究动机与目标

• 为应用于非局部 Cahn–Hilliard 方程的首阶稳定化线性半隐式格式建立严格的收敛性估计。
• 通过在误差分析中使用离散 $H^{-1}$ 范数，克服非局部 Cahn–Hilliard 方程中非线性项带来的挑战。
• 通过结合高阶一致性分析与逆不等式，推导出数值解的均匀 $\ell^\infty$ 有界性。
• 在对稳定化参数 $A$ 的改进条件下，严格证明该格式的能量稳定性。
• 通过扰动论证恢复数值解的 $\ell^\infty$ 有界性，确保收敛性分析的有效性。

提出的方法

• 该格式采用稳定化的线性半隐式时间离散化与傅里叶伪谱空间离散化，以实现谱精度。
• 通过在离散 $H^{-1}$ 范数下分析误差函数，避免了对非线性项复杂性的直接处理。
• 使用高阶一致性展开以捕捉时间与空间的截断误差，从而支持逆不等式的应用。
• 将逆不等式应用于 $H^{-1}$ 误差估计，以推导出误差函数的 $\ell^\infty$ 有界性。
• 将数值解视为一个扰动问题的精确解，从而在较宽松的时间步长约束下推导出均匀 $\ell^\infty$ 有界性。
• 通过基于解的 $\ell^\infty$ 有界性推导稳定化参数 $A$ 的新下界，从而改进能量稳定性条件。

实验结果

研究问题

• RQ1 对于非局部 Cahn–Hilliard 方程的首阶稳定化线性半隐式格式，是否能在放宽时间步长约束的条件下实现最优收敛？
• RQ2 如何有效利用离散 $H^{-1}$ 范数在收敛性分析中处理非线性项，而无需强正则性假设？
• RQ3 高阶一致性分析在使逆不等式可用于 $\ell^\infty$ 范数误差有界性推导中起到何种作用？
• RQ4 在何种条件下，该格式的能量稳定性可被严格证明？
• RQ5 能否通过扰动论证与一致性分析建立数值解的 $\ell^\infty$ 有界性？

主要发现

• 该格式在时间方向实现一阶收敛、空间方向实现谱精度，误差估计为 $\|\hat{e}^{n+1}\|_{-1,N} + \bigl(\gamma_0 \Delta t \sum_{k=1}^{n+1} \|\hat{e}^k\|_2^2\bigr)^{1/2} \leq C^* (\Delta t^3 + h^m)$。
• 在时间步长满足线性缩减路径 $\Delta t \leq Ch$ 的条件下，若 $h \leq 1/(2C_4 C^*)$，则误差函数的 $\ell^\infty$ 有界性为 $\|\hat{e}^{n+1}\|_{\infty} \leq C_4 C^* h$。
• 推导出数值解的均匀 $\ell^\infty$ 有界性，确保在对 $A$ 的改进假设下能量稳定性条件的有效性。
• 能量稳定性条件被改进为 $A \geq \max\bigl\{ \frac{18M_0^4}{\gamma_0}, \frac{3}{2}M_0^2 - \frac{1}{2} \bigr\}$，其中 $M_0 = 1 + \max_{1 \leq k \leq N_k} \|u_N^k\|_{\infty}$。
• 通过在每一步时间推进中恢复先验假设 (68)，完成收敛性分析，确保归纳法成立。
• 该方法对低阶空间格式同样具有鲁棒性，此时高阶一致性分析对于捕捉空间截断误差至关重要。