[論文レビュー] Fully Adaptive Newton-Galerkin Time Stepping Methods for Singularly Perturbed Parabolic Evolution Equations

Mario Amrein, Thomas P. Wihler|arXiv (Cornell University)|Oct 2, 2015

Differential Equations and Numerical Methods参考文献 13被引用数 3

ひとこと要約

本稿では、境界層やスパイクを示す特異摂動半線形放物型方程式に対して、$\varepsilon$-ロバストな誤差制御を実現する完全に適応的なニュートン・ガレルキン時間刻み法を提案する。ニュートン線形化、適応的後退オイラー時間刻み、空間的有限要素の細分化を組み合わせることで、境界層やスパイクが存在する中でも$\varepsilon$-ロバストな収束性を達成し、変動する摂動パラメータを用いた数値実験により検証された。

ABSTRACT

In this paper we develop an adaptive procedure for the numerical solution of semilinear parabolic problems, with possible singular perturbations. Our approach combines a linearization technique using Newton's method with an adaptive discretization-which is based on a spatial finite element method and the backward Euler time stepping scheme-of the resulting sequence of linear problems. Upon deriving a robust a posteriori error analysis, we design a fully adaptive Newton-Galerkin time stepping algorithm. Numerical experiments underline the robustness and reliability of the proposed approach for various examples.

研究の動機と目的

• 微小な拡散係数（$\varepsilon \ll 1$）を有する半線形放物型問題を効率的かつ信頼性高く解くための数値法を開発すること。このような問題は境界層、スパイク、またはショックを示すことがある。
• 標準的手法が安定性または精度の喪失により失敗する特異摂動領域におけるロバストな数値近似の課題に取り組むこと。
• 線形化、時間刻み、空間離散化誤差のバランスを最適化するための適応的アルゴリズムを設計すること。
• すべての摂動領域で信頼性のある誤差監視が可能となるよう、$\varepsilon \to 0$ においても誤差推定器がロバストであることを保証すること。
• 既知の特異挙動を示すベンチマーク問題（非線形源項や爆発解を含む）に対して、手法の性能を検証すること。

提案手法

• 半線形放物型問題に対してニュートン法を適用し、線形化することで、線形時間発展問題の系列に変換する。
• 線形化された問題は、後退オイラースキームによる時間離散化と$\mathbb{P}_1$-有限要素による空間離散化を用いて離散化する。
• 合計残差を線形化誤差、時間離散化誤差、空間離散化誤差の3つの計算可能な成分に分解する、後験的誤差解析を導出する。
• 特異摂動楕円型問題および線形放物型理論の技術を活用することで、誤差推定器を$\varepsilon$-ロバストに構築する。
• 各反復で支配的となる誤差成分に基づき、時間ステップまたは空間メッシュを細分化する完全に適応的なアルゴリズムを設計する。
• 残差指標に基づき、ニュートンステップサイズ、時間ステップ、空間メッシュを動的に調整することで、効率性とロバスト性を確保する。

実験結果

リサーチクエスチョン

• RQ1 半線形放物型問題における特異摂動パラメータ$\varepsilon$に対してロバスト性を維持する完全に適応的なニュートン・ガレルキン法を構築可能か？
• RQ2 効果的な適応的制御のため、合計誤差を線形化誤差、時間離散化誤差、空間離散化誤差に明確に分解できるか？
• RQ3 提案された適応的アルゴリズムは、$\varepsilon \to 0$ の特異摂動領域において境界層や内部スパイクを信頼性高く解像できるか？
• RQ4 誤差推定器は、$\varepsilon$ の異なる値においても、有界かつ安定した効率インデックスを提供できるか？
• RQ5 爆発や急激な動的挙動を示す問題において、推定誤差の時間的挙動はどのように変化するか？

主な発見

• 提案された適応的アルゴリズムは$\varepsilon$-ロバストな収束性を達成しており、$p = 1,\dots,5$ に対して$\varepsilon = 10^{-p}$ のすべてのテスト値で効率インデックスが有界に保たれている。
• 推定誤差は時間とともに$\mathcal{O}(t_n^{1/2})$ の割合で増大し、誤差推定器の理論的予想と整合的である。
• $\varepsilon = 10^{-5}$ の場合に鋭い境界層、$\varepsilon = 10^{-3}$ の場合に内部スパイクを効果的に解像している。
• 爆発問題（$\beta = 4$）において、時間の経過に伴い$x=2$付近の発展するスパイクを適応的手法が正確に捉えている。
• 推定誤差のlog/logプロットでは勾配が$1/2$となり、理論的誤差増大率を確認し、誤差推定器の信頼性を裏付けた。
• 適応戦略は線形化、時間、空間の細分化を効果的にバランスさせ、最適な計算複雑性とロバストな性能を達成している。