[論文レビュー] A posteriori error estimators for hierarchical B-spline discretizations

Annalisa Buffa, Eduardo M. Garau|arXiv (Cornell University)|Nov 23, 2016

Advanced Numerical Analysis Techniques参考文献 11被引用数 5

ひとこと要約

本稿では、線形第二順序楕円型問題の階層的Bスプライン離散化に対する信頼性があり関数に基づく事後誤差推定法を提示する。Bスプライン重み付きPoincaré不等式を活用し、基底関数の切断を回避することで、メッシュ構成の制限なしにエネルギー誤差のグローバルな上界を提供する。これにより、多項式次数にかかわらず最適収束率を達成する効率的な適応的グリッド・エンリッチメントが可能となる。

ABSTRACT

In this article we develop function-based a posteriori error estimators for the solution of linear second order elliptic problems considering hierarchical spline spaces for the Galerkin discretization. We prove a global upper bound for the energy error. The theory hinges on some weighted Poincaré type inequalities, where the B-spline basis functions are the weights appearing in the norms. Such inequalities are derived following the lines in [Veeser and Verfürth, 2009], where the case of standard finite elements is considered. Additionally, we present numerical experiments that show the efficiency of the error estimators independently of the degree of the splines used for discretization, together with an adaptive algorithm guided by these local estimators that yields optimal meshes and rates of convergence, exhibiting an excellent performance.

研究の動機と目的

• 階層的Bスプラインガレルキン離散化に対する、メッシュ構成の制限なしに信頼性のある事後誤差推定法を開発すること。
• Bスプライン重み付きPoincaré不等式から導出された関数ベースの推定法を用いて、エネルギー誤差のグローバルな上界を確立すること。
• 基底関数の切断を回避し、凸かつ単純なサポートを維持することで、アイソジオメトリック解析における効率的かつ安定的な適応的グリッド・エンリッチメントを可能にすること。
• 数値実験を通じて、多項式次数に応じて最適収束率と推定法の効率性を示すこと。
• Bスプライン基底関数を重み関数として用いることで、高次スプラインに対するPoincaré型不等式を一般化すること。

提案手法

• Bスプライン基底関数を重み関数として用いる、新規のBスプライン重み付きPoincaré型不等式を導出。重み関数としてのBスプライン基底関数が、ノルムに組み込まれる。
• 切断なしの階層的スプライン基底を用い、数値的安定性と効率性を高めるために、凸かつ連結したサポートを維持する。
• 局所的な残差に基づく誤差指標を個々の基底関数に局所化することで、関数ベースの事後誤差推定法を構築する。
• 導出された重み付きPoincaré不等式を用いてエネルギー誤差を上界で抑え、グローバル誤差推定法の信頼性を証明する。
• 階層的メッシュの細分化において推定法の縮小性を導入し、収束を保証する。
• 局所誤差推定法に従う適応的アルゴリズムを採用。誤差の均等分配を実現し、最適収束率を達成する。

実験結果

リサーチクエスチョン

• RQ1 階層的Bスプライン離散化に対して、階層的メッシュ構成の制限なしに信頼性のある事後誤差推定法を構築できるか？
• RQ2 エネルギー誤差を階層的スプライン空間で上界で抑えるために、どのようにBスプライン重み付きPoincaré不等式を導出できるか？
• RQ3 提案された関数ベース誤差推定法は、多項式次数が変化しても効率性と信頼性を維持できるか？
• RQ4 これらの推定法に従う適応的アルゴリズムは、実際の計算において最適収束率を達成できるか？
• RQ5 推定法の効率性と計算コストの面で、標準的な残差ベースの手法と比較して、性能にどのような差があるか？

主な発見

• 提案された事後誤差推定法は、階層的メッシュ構成に制限がなく、エネルギー誤差のグローバルな上界を提供する。
• 実験的に高い効率性を示し、すべてのテストケースにおいて効率指数が0.8から1.2の間で保たれ、真のエネルギー誤差とよく一致している。
• 推定法に従う適応的グリッド・エンリッチメントは、二乗、立方、四次スプラインにおいても、解の特異性が存在する場合でさえ最適収束率を達成する。
• 例7.5（四分の1円環）では、二乗スプラインを用いた適応的グリッド・エンリッチメントで、均一なグリッド・エンリッチメントに比べて必要な自由度の3％で同等の精度が達成された。
• 基底関数の切断を回避することで、単純で凸なサポートを維持し、非凸または非連結なサポートに伴う計算コストを排除した。
• 数値結果から、高次スプラインを含むすべてのテストされた多項式次数において、推定法が有効かつ効率的であることが確認された。