[論文レビュー] 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}$ ノルムを採用し、高次一致性解析と逆不等式を組み合わせることで、時間方向に最適な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 方程式に適用された1次安定化線形半陽的スキームの厳密な収束推定を確立すること。
• 非局所 Cahn–Hilliard 方程式における非線形項の複雑さに起因する課題を、誤差解析に離散 $H^{-1}$ ノルムを用いることで克服すること。
• 高次一致性解析と逆不等式を組み合わせることで、数値解の一様な $\ell^\infty$ バインドを導出すること。
• 安定化パラメータ $A$ に関する精密な条件の下で、スキームのエネルギー安定性を正当化すること。
• 摂動論法により数値解の $\ell^\infty$ バインドを回復させることで、収束解析の有効性を保証すること。

提案手法

• スキームは、スペクトル精度を達成するため、フーリエ擬スペクトル法を用いた空間離散化と、安定化線形半陽的時間離散化を組み合わせる。
• 非線形項の複雑さを直接取り扱うのを避けるために、誤差関数を離散 $H^{-1}$ ノルムで解析する。
• 時間的および空間的切り捨て誤差を捉えるために、高次一致性展開を用い、逆不等式の適用を可能にする。
• 逆不等式を $H^{-1}$ 誤差推定に適用し、誤差関数の $\ell^\infty$ バインドを導出する。
• 数値解を摂動問題の正確な解とみなすことにより、時間ステップ制約が緩い条件下でも一様な $\ell^\infty$ バインドを導出できる。
• 数値解の $\ell^\infty$ バインドに基づいて、安定化パラメータ $A$ の新たな下界を導出し、エネルギー安定性条件を精緻化する。

実験結果

リサーチクエスチョン

• RQ1 非局所 Cahn–Hilliard 方程式に対する1次安定化線形半陽的スキームは、緩い時間ステップ制約のもとで最適収束を達成できるか？
• RQ2 離散 $H^{-1}$ ノルムを効果的に用いることで、強い正則性仮定を必要とせずに非線形項を誤差解析で扱えるか？
• RQ3 高次一致性解析が、逆不等式を用いて $\ell^\infty$ ノルムにおける誤差をバインドする際に果たす役割は何か？
• RQ4 安定化パラメータ $A$ に対してどのような条件下で、スキームのエネルギー安定性が厳密に正当化されるか？
• RQ5 摂動論法と一致性解析を用いて、数値解の $\ell^\infty$ バインドを確立できるか？

主な発見

• スキームは時間方向に1次収束、空間方向にスペクトル精度を達成し、誤差推定式 $\|\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) を回復することで、収束解析が完了し、帰納法の正当性が保証される。
• 本手法は低次の空間スキームに対してもロバストであり、空間切り捨て誤差を捉えるために高次一致性解析が不可欠である。