[論文レビュー] Error analysis of linearized semi-implicit Galerkin finite element methods for nonlinear parabolic equations

Buyang Li, Weiwei Sun|arXiv (Cornell University)|Aug 23, 2012

Advanced Numerical Methods in Computational Mathematics参考文献 19被引用数 108

ひとこと要約

本稿では、非線形放物型方程式、特に3次元ジュール加熱系に適用される線形化半陰的ガレルキン有限要素法に向けた、新しい誤差解析フレームワークを提示する。誤差を時間的および空間的成分に分割し、時間離散系を用いることで、時間ステップ制限なしに、無条件に最適な$L^2$および$H^1$誤差推定を確立した。これは、従来の研究で強いノルムにおける有界性を保証するために制限的な時間ステップ条件を必要としていたという、主な限界を克服したものである。

ABSTRACT

This paper is concerned with the time-step condition of commonly-used linearized semi-implicit schemes for nonlinear parabolic PDEs with Galerkin finite element approximations. In particular, we study the time-dependent nonlinear Joule heating equations. We present optimal error estimates of the semi-implicit Euler scheme in both the $L^2$ norm and the $H^1$ norm without any time-step restriction. Theoretical analysis is based on a new splitting of the error and precise analysis of a corresponding time-discrete system. The method used in this paper can be applied to more general nonlinear parabolic systems and many other linearized (semi)-implicit time discretizations for which previous works often require certain restriction on the time-step size $τ$.

研究の動機と目的

• 非線形放物型PDEの線形化半陰的スキームの誤差解析において一般的に課される時間ステップ制限を解消すること。
• 時間ステップ制限を必要とする逆不等式と$L^∞$-ノルムの有界性に基づく帰納法に依存するのを回避すること。
• 時間ステップサイズ$\tau$に何ら制限を設けずに、$L^2$および$H^1$ノルムにおける最適誤差推定を確立すること。
• ジュール加熱モデルを越えて、より広範な非線形放物型系に適用可能な一般化されたフレームワークを開発すること。

提案手法

• 時間離散放物系を用いて、数値誤差を時間的および空間的成分に分割する新しい手法を導入する。
• 時間離散系に対するガレルキン有限要素誤差を、時間ステップ$\tau$に依存しない空間メッシュサイズ$h$にのみ依存する境界で解析する。
• 帰納法と逆不等式を用いて、事前有界性を仮定せずに$L^\infty$-ノルムの誤差を制御し、時間ステップ制限を回避する。
• グロウォールの不等式を適用して、$L^2$および$H^1$ノルムにおける最適誤差推定を導出する。
• 非線形項$\sigma(u)$を制御するため、$W^{1,p}$推定と補間誤差境界を用いる。
• 帰納法とメッシュ依存の逆不等式を用いて、離散解の強いノルムにおける有界性を確立し、無条件的な誤差制御を可能にする。

実験結果

リサーチクエスチョン

• RQ1 非線形放物型方程式の線形化半陰的ガレルキン有限要素スキームに対して、時間ステップ制限なしに最適誤差推定を導出できるか？
• RQ2 新しい誤差分割戦略を用いることで、誤差解析における数値解の$L^\infty$-ノルムの有界性の必要性を排除できるか？
• RQ3 3次元非線形ジュール加熱系に対して、$H^1$-ノルム誤差推定を無条件に得ることは可能か？
• RQ4 提案されたフレームワークを他の非線形放物型系および高次時間積分スキームに拡張できるか？
• RQ5 新しい解析フレームワーク下で、誤差が空間メッシュサイズ$h$および時間ステップ$\tau$にどのように依存するか？

主な発見

• 時間ステップサイズ$\tau$に何ら制限を設けずに、$O(h^{3/2})$の最適$L^2$誤差推定が得られた。
• 時間ステップ制限なしに、$O(h)$の最適$H^1$誤差推定が無条件に確立された。
• $L^2$誤差境界$\|e_h^n\|_{L^2} \leq Ch^{3/2}$は、帰納法とグロウォールの不等式を用いて証明され、$h$依存の境界にのみ依存する。
• 時間離散系と逆不等式を用いることで、誤差の$L^\infty$-ノルム制御の必要性を回避し、時間ステップ条件の不要性を実現した。
• 解析フレームワークは一般性を持ち、他の非線形放物型系や高次時間積分スキームへの応用が可能である。
• 主な革新点は、誤差の分割と時間離散系の導入による時間的・空間的誤差寄与の分離であり、これにより無条件の安定性と最適性が実現された。