[論文レビュー] Error estimates in balanced norms of finite element methods for higher order reaction-diffusion problems

Sebastian Franz, Hans-G. Roos|arXiv (Cornell University)|Feb 26, 2019

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

ひとこと要約

本稿は、シシュキンメッシュ上での層を多く含む解に特化した新しい非標準的射影作用素を導入することにより、特異摂動反応拡散問題に高次有限要素法を適用した場合のバランスドノルムにおける最適誤差推定を確立する。著者らは、バランスドノルムにおいて $O((N^{-1} \max|\psi'|)^{p+1-m})$ の収束率を達成し、既存の結果を改善することで対数因子を排除し、2次元における高次問題への解析を拡張した。

ABSTRACT

Error estimates of finite element methods for reaction-diffusion problems are often realised in the related energy norm. In the singularly perturbed case, however, this norm is not adequate. A different scaling of the $H^m$ seminorm for $2m$-th order problems leads to a balanced norm which reflects the layer behaviour correctly. We prove error estimates in such balanced norms and improve thereby existing estimates known in literature.

研究の動機と目的

• 特異摂動高次反応拡散問題におけるエネルギーノルムの層挙動の捉えにくさを是正すること。
• 境界層およびコーナー層の漸近的挙動を正しく反映するバランスドノルムにおける最適誤差推定を導出すること。
• 層適応メッシュ上での安定性と精度を維持する上で、標準的射影（例：$L^2$-射影）の限界を克服すること。
• 2次元における2次元問題から高次（$2m$次）問題へのバランスドノルム誤差推定理論の拡張すること。
• 対数因子を含まない収束率の確立。既存の結果に $\mathcal{O}(\ln N)^{1/2}$ 項が含まれていたのを改善すること。

提案手法

• 有限要素空間 $V^N$ に、層領域におけるヘルミート補間と粗い領域におけるリッツ射影を組み合わせた新しい射影作用素 $P$ を導入する。
• 射影を領域ごとに定義：層領域 $\Omega_i^*$ では $Pw_i = Iw_i$、外部では $Pw_i = 0$、外領域では $Pv = Iv$、粗い領域 $\Omega_c$ では $Pv = \pi v$。
• メッシュ生成関数 $\varphi$ と遷移パラメータ $\lambda = \sigma \varepsilon \ln N$ を用い、$\sigma \geq 2m$ として最適解像度を達成するシシュキン型メッシュを定義する。
• 層領域およびコーナー領域における補間誤差 $\|Iv - Pv\|_{W^{1,\infty}}$ の境界を確立し、$\lesssim N^{-(\sigma-1)}$ および $\lesssim (1+hN)N^{-3}$ を示す。
• 双線形形式 $\tilde{a}(\cdot,\cdot)$ に対する安定性および一致性の推定を導出し、$\varepsilon^{1/2}(1+hN)(h + N^{-1}\max|\psi'|)^3$ を含む境界を得る。
• すべての推定を組み合わせ、バランスドノルムにおける最終的な誤差境界を導出：$\sigma \geq 4$ の場合 $\|u - u^N\|_b \lesssim (1+hN)(h + N^{-1}\max|\psi'|)^2$、一般には $\sigma \geq p+1$ の場合 $\lesssim (N^{-1}\max|\psi'|)^{p+1-m}$。

実験結果

リサーチクエスチョン

• RQ1 高次（$2m$次）特異摂動反応拡散問題に対して、バランスドノルムにおける最適誤差推定を導出可能か？
• RQ2 シシュキンメッシュ上での高次FEMにおいて、誤差推定から対数因子を排除することは可能か？
• RQ3 層構造を尊重し、バランスドノルム誤差解析を可能にする安定的かつ高精度な射影作用素を構築する方法は何か？
• RQ4 層適応メッシュ上での $C^{m-1}$-ヘルミート有限要素に対して、バランスドノルムにおける最適収束率は何か？
• RQ5 1次元から2次元への解析拡張は可能か？ その際、最適収束率を維持できるか？

主な発見

• 本稿は、$2m$次反応拡散問題におけるバランスドノルム $\|v\|_b = \varepsilon^{k-1/2}|v|_{H^m} + \|v\|_{H^{m-k}}$ における最適誤差推定を確立した。
• 特に $m=2$, $k=1$, および $\sigma \geq 4$ の場合、バランスドノルムにおける誤差は $\|u - u^N\|_b \lesssim (1+hN)(h + N^{-1}\max|\psi'|)^2$ を満たす。
• $p=3$（双3次ヘルミート要素）の場合、収束次数は $p+1-m = 2$ であり、理論的境界と一致する。
• 従来の推定に含まれていた $\mathcal{O}((\ln N)^{1/2})$ 要因を排除し、純粋な多項式収束を達成した。
• 一般に $m \geq k \geq 1$ および多項式次数 $p \geq 2m$ の場合、$\sigma \geq p+1$ のとき誤差境界は $\|u - u^N\|_b \lesssim (N^{-1}\max|\psi'|)^{p+1-m}$ となる。
• 提案された射影作用素 $P$ がバランスドノルム枠組み内での最適安定性および近似性質を保証することを解析が確認した。