[論文レビュー] Asymptotically compatible energy of variable-step fractional BDF2 formula for time-fractional Cahn-Hilliard model

Hong-lin Liao, Nan Liu+1 more|arXiv (Cornell University)|Oct 22, 2022

Fractional Differential Equations Solutions被引用数 6

ひとこと要約

本稿は、Caputo微分の次数 $\alpha \in (0,1)$ を持つ時間分数 Cahn-Hilliard モデルに対して、変化する時間ステップを用いた分数階 BDF2 スキームを提案する。分数階 BDF2 公式の局所的・非局所的分割を用いて、新たな離散勾配構造を導入し、弱いステップ比制約の下で離散エネルギー散逸則を確立する。$\alpha \to 1^{-}$ の極限において、古典的 Cahn-Hilliard エネルギー則と漸近的に整合性を保ち、適応的時間刻み法を用いることで長時間安定性と精度を確保する。

ABSTRACT

A new discrete energy dissipation law of the variable-step fractional BDF2 (second-order backward differentiation formula) scheme is established for time-fractional Cahn-Hilliard model with the Caputo's fractional derivative of order $α\in(0,1)$, under a weak step-ratio constraint $0.4753\le τ_k/τ_{k-1}

研究の動機と目的

• 時間分数 Cahn-Hilliard モデルに対して、一般の非一様時間メッシュを用いた安定的で2次精度の時間離散化を構築すること。
• 変化する時間ステップ下でも物理的エネルギー減少性を保つ離散エネルギー散逸則を確立すること。
• $\alpha \to 1^{-}$ の極限において、古典的 Cahn-Hilliard エネルギー則と漸近的に整合性を保つこと。
• 初期に急激な変化が生じる領域と、後に粗大化する段階を的確に捉える適応的時間刻み法を用いて、長時間シミュレーションを効率的に行うこと。

提案手法

• 分数階 BDF2 公式を、古典的2段階後退差分に類似する局所的成分と、Caputo微分のL1公式に類似する非局所的成分に分離する局所的・非局所的分割を導入する。
• この分割に基づき、一貫性があり安定なエネルギー散逸則を導出する新たな離散勾配構造を構築する。
• 安定性およびエネルギー散逸のための弱いステップ比制約 $0.4753 \leq \tau_k / \tau_{k-1} < r^*(\alpha)$ を導出する。ここで $r^*(\alpha) \geq 4.660$ であり、$\alpha \in (0,1)$ に対して成り立つ。
• 初期特異性と動的変化を解像するために、$\tau_{\text{ada}} = \max\{\tau_{\min}, \tau_{\max}/\Pi_\phi\}$ で与えられる適応的時間刻み法を実装する。ここで $\Pi_\phi = \sqrt{1 + \eta \|\partial_\tau \phi^n\|^2}$ である。
• 空間離散化にフーリエ擬スペクトル法（$128 \times 128$ グリッド）を用い、$\varepsilon = 10^{-12}$ を用いて分数階 BDF2 公式の高速評価に SOE 技法を適用する。

実験結果

リサーチクエスチョン

• RQ1 時間分数 Cahn-Hilliard モデルに対して、離散エネルギー散逸則を保つ2次精度の変化する時間ステップの分数階 BDF2 スキームを構築可能か？
• RQ2 離散エネルギー則は、$\alpha \to 1^{-}$ の極限において、古典的 Cahn-Hilliard エネルギー散逸則と漸近的に整合性を保つか？
• RQ3 重み付きメッシュおよび適応的メッシュ上でも、スキームが最適な時間精度 $O(\tau^{\min\{\gamma, 3 - \alpha\}})$ を維持できるか？
• RQ4 初期特異性および長時間にわたる粗大化ダイナミクスを解像するための適応的時間刻み法は、どの程度有効か？

主な発見

• 提案された変化する時間ステップの分数階 BDF2 スキームは、重み付きメッシュおよび適応的メッシュ上でも2次精度 $O(\tau^{\min\{\gamma, 3 - \alpha\}})$ を達成しており、$\alpha = 0.4$ および $\alpha = 0.7$ に対して数値的に確認された。
• 弱いステップ比制約 $0.4753 \leq \tau_k / \tau_{k-1} < r^*(\alpha)$ の下で、新たな離散エネルギー散逸則が確立された。ここで $r^*(\alpha) \geq 4.660$ であり、$\alpha \in (0,1)$ に対して成り立つ。
• 離散エネルギーおよびエネルギー散逸則は、$\alpha \to 1^{-}$ の極限において、古典的 Cahn-Hilliard モデルと漸近的に整合性を保ち、整数階極限での一貫性を保証する。
• $\eta = 10^3$ を用いた適応的時間刻み法により、初期段階の急激なエネルギー散逸と長時間にわたる粗大化ダイナミクスが的確に捉えられ、時間ステップが解の勾配に応じて動的に調整された。
• 数値結果から、粗大化速度および相分離の進化は $\alpha$ に強く依存しており、$\alpha \to 1^{-}$ の極限に近づくにつれて、古典的 $O(t^{-1/3})$ の速度に近づくことが示された。