[論文レビュー] Maximum bound principles for a class of semilinear parabolic equations and exponential time differencing schemes

Qiang Du, Lili Ju+2 more|arXiv (Cornell University)|May 23, 2020

Advanced Mathematical Modeling in Engineering参考文献 93被引用数 12

ひとこと要約

本稿は、半線形放物型方程式における最大境界原理（MBP）の理論的枠組みを確立し、初期条件および境界条件がそのような境界を満たす場合に、解が絶対値で一様に有界に保たれることを証明する。時間に依存しないMBPを保証する指数型時間差分（ETD）スキームを構築し、時間ステップサイズに依存しない形でMBPを保証する。安定性と高次精度を両立させ、スカラー、ベクトル、行列値系において、厳密な誤差推定とエネルギー安定性を証明する。

ABSTRACT

The ubiquity of semilinear parabolic equations has been illustrated in their numerous applications ranging from physics, biology, to materials and social sciences. In this paper, we consider a practically desirable property for a class of semilinear parabolic equations of the abstract form $u_t=\mathcal{L}u+f[u]$ with $\mathcal{L}$ being a linear dissipative operator and $f$ being a nonlinear operator in space, namely a time-invariant maximum bound principle, in the sense that the time-dependent solution $u$ preserves for all time a uniform pointwise bound in absolute value imposed by its initial and boundary conditions. We first study an analytical framework for some sufficient conditions on $\mathcal{L}$ and $f$ that lead to such a maximum bound principle for the time-continuous dynamic system of infinite or finite dimensions. Then, we utilize a suitable exponential time differencing approach with a properly chosen generator of contraction semigroup to develop first- and second-order accurate temporal discretization schemes, that satisfy the maximum bound principle unconditionally in the time-discrete setting. Error estimates of the proposed schemes are derived along with their energy stability. Extensions to vector- and matrix-valued systems are also discussed. We demonstrate that the abstract framework and analysis techniques developed here offer an effective and unified approach to study the maximum bound principle of the abstract evolution equation that cover a wide variety of well-known models and their numerical discretization schemes. Some numerical experiments are also carried out to verify the theoretical results.

研究の動機と目的

• 抽象的半線形放物型方程式に対して、時間に依存しない最大境界原理（MBP）を保証するための、線形散乱作用素 𝒮 および非線形作用素 f に関する十分な解析的条件を確立すること。
• 時間ステップサイズに依存しない形でMBPを保証する、第一および第二の指数型時間差分（ETD）スキームを構築すること。
• 提案されたETDスキームのスカラーおよびベクトル／行列値設定における、厳密な誤差推定およびエネルギー安定性を導出・証明すること。
• ディリクレ境界条件、ノイマン境界条件、周期的境界条件を含むさまざまな境界条件に適応可能な抽象的枠組みを拡張し、異なる行列ノルム（2ノルムおよびフロベニウスノルム）に関するMBPを分析すること。
• 積分因子（IF）スキームにおいて、時間に依存しないMBPを保証する可能性を検討し、ETD法と比較すること。

提案手法

• 収縮半群理論に基づき、線形作用素 𝒮 および非線形作用素 f に関する解析的条件を導出し、解がすべての時間において絶対値で有界に保たれることを保証する。
• 収縮半群の生成子を用いてETDスキームを構築し、時間離散化においてMBPを時間に依存せずに保証する。
• 指数積分法を適用して、時間ステップサイズに制限なく第一および第二の時間精度を達成しつつ、MBPを維持する。
• エネルギー安定性解析および誤差推定を用いて、さまざまな境界条件下でのETDスキームの収束性および頑健性を検証する。
• ベクトルおよび行列値系に枠組みを拡張し、ℓ²、ℓ∞、行列2ノルムおよびフロベニウスノルムに関するMBPを分析する。
• ETD法と積分因子（IF）スキームを比較し、IFスキームは時間ステップ制限のもとでMBPを条件付きで保証できるが、ETDスキームは時間に依存しない形でMBPを保証できることを示す。

実験結果

リサーチクエスチョン

• RQ1 線形作用素 𝒮 および非線形作用素 f がどのような条件下に置かれれば、半線形放物型方程式が時間連続設定において最大境界原理（MBP）を満たすか？
• RQ2 時間ステップサイズに依存しない形でMBPを保証する指数型時間差分（ETD）スキームを構築できるか？
• RQ3 提案されたETDスキームのスカラー、ベクトル、行列値系における誤差推定およびエネルギー安定性の性質は何か？
• RQ4 行列値系において、2ノルムおよびフロベニウスノルムに関してMBPはどのように振る舞うか？
• RQ5 積分因子（IF）スキームは時間に依存しない形でMBPを保証できるか？もしそうでない場合、どのような条件下で条件付きで保証できるか？

主な発見

• 本稿は、初期および境界データが境界を満たす限り、抽象的半線形放物型方程式の解がすべての時間において絶対値で有界に保たれることを保証するための十分条件を確立した。
• 第一および第二のETDスキームが構築され、時間ステップサイズに依存せずMBPを保証する（安定性制限なし）。
• 理論的誤差推定が導出され、解およびデータの適切な正則性仮定のもとで、ETDスキームの最適収束率が示された。
• ETDスキームのエネルギー安定性が証明され、離散エネルギーが時間とともに増大しないことが保証された。これは、長期的な数値的信頼性を支持する。
• 行列値系において、2ノルムに関するMBPはETDスキームによって時間に依存せず保証されるが、フロベニウスノルムに関しては非線形項のバインディングに課題があるため、未解決の問題のままである。
• IFスキームは時間ステップ制限のもとでMBPを条件付きで保証できるが、適切な安定化子を用いることでIFスキームでも時間に依存しないMBP保証が達成可能である可能性が示唆され、今後の研究の方向性として提示された。