[論文レビュー] Guaranteed error control bounds for the stabilised space-time IgA approximations to parabolic problems

Ulrich Langer, Svetlana Matculevich+1 more|arXiv (Cornell University)|Dec 16, 2017

Advanced Numerical Methods in Computational Mathematics参考文献 9被引用数 6

ひとこと要約

本稿では、ストレクチャライズド変換を用いた空間時間等値的有限要素法（IgA）近似における放物型初期境界値問題の、保証付きで完全に計算可能な後験的誤差境界を提示する。関数解析と積分恒等式に基づき、メッシュに依存しない誤差上限（majorants）が得られ、エネルギー空間内の任意の近似に対して有効である。これにより、局所誤差指標を用いた信頼性高く効率的な適応的メッシュ・エンリッチメントが可能となり、数値結果により、複雑で特異的かつ非自明な幾何形状においても高い効率性と信頼性が確認された。

ABSTRACT

The paper is concerned with space-time IgA approximations of parabolic initial-boundary value problems. We deduce guaranteed and fully computable error bounds adapted to special features of IgA approximations and investigate their applicability. The derivation method is based on the analysis of respective integral identities and purely functional arguments. Therefore, the estimates do not contain mesh-dependent constants and are valid for any approximation from the admissible (energy) class. In particular, they provide computable error bounds for norms associated with stabilised space-time IgA approximations as well as imply efficient error indicators enhancing the performance of fully adaptive solvers. The last section of the paper contains a series of numerical examples where approximate solutions are recovered by IgA techniques. The mesh refinement algorithm is governed by a local error indicator generated by the error majorant. Numerical results discussed in the last section illustrate both reliability, as well as the quantitative efficiency of the error estimates presented.

研究の動機と目的

• 安定化された空間時間IgA近似に対する信頼性があり、完全に計算可能な後験的誤差推定を確立すること。
• 誤差境界がメッシュサイズに依存しない、許容可能なエネルギー空間内の任意の近似に対して有効であることを保証すること。
• 誤差上限（majorants）から導出される局所誤差指標を用いて、効率的な適応的メッシュ・エンリッチメントを可能にすること。
• 複雑で非自明な幾何形状および特異的解を有する状況において、提案された誤差推定の効率性と信頼性を示すこと。

提案手法

• 積分恒等式と純粋な関数的議論を用いた、関数型の後験的誤差推定の導出。メッシュ依存定数を回避する。
• ストリームライン拡散および時間上向き（time-upwind）のテスト関数に基づく、安定化された空間時間変分形式の適用。
• エネルギー誤差ノルムの保証付き上界（majorants）および下界（minorants）の構築。
• IgAにおける適応的メッシュ・エンリッチメントのための誤差指標として、上限の局所的寄与を用いる。
• 非構造的で局所的なメッシュ・エンリッチメントを実現するため、THBスプラインを用いた実装。
• 上限の最小化により、非定数の効率インデックスが1に近い非常に正確な上界を得る。

実験結果

リサーチクエスチョン

• RQ1 安定化された空間時間IgA近似に対する、メッシュ依存定数を含まない保証付きで完全に計算可能な誤差境界を導出できるか？
• RQ2 既存の後験的境界と比較して、提案された関数型誤差推定は信頼性および効率性においてどのように性能を発揮するか？
• RQ3 上限の局所的寄与は、IgAにおける適応的メッシュ・エンリッチメントのための効果的かつ信頼性のある誤差指標として機能できるか？
• RQ4 空間時間IgAにおける解の構築と比較して、誤差境界の評価に要する計算コストはどの程度か？
• RQ5 誤差上限に基づく適応的メッシュ・エンリッチメント戦略は、特異的解および複雑な幾何形状を有する問題に対してどの程度効果的か？

主な発見

• 提案された関数的後験的誤差推定は、保証付きで完全に計算可能であり、メッシュ依存定数なしにエネルギー空間内の任意の近似に対して有効である。
• 誤差上限（majorant）の効率インデックスは常に1に近く、高品質な上界であることを示している。
• 局所的上限寄与に基づく適応的メッシュ・エンリッチメントは、最適収束率を達成し、均一メッシュ・エンリッチメントと比較して収束順序が向上している。
• テスト例において、誤差境界の計算コストはIgA解の計算コストのわずか一部にとどまった。
• 空間時間THBスプラインに基づく適応的スキームは、特異性および複雑な幾何形状を有する問題を効率的に処理でき、さまざまなテストケースにおいて頑健な性能を示した。