[論文レビュー] Numerical study of goal-oriented error control for stabilized finite element methods

Marius Paul Bruchhäuser, Kristina Schwegler+1 more|arXiv (Cornell University)|Mar 27, 2018

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

ひとこと要約

本稿では、安定化有限要素法と組み合わせることで、対流優勢問題における目的関数指向の適応型有限要素法を提案する。Dual Weighted Residual (DWR) 法を用い、目的関数のためのロバストで計算可能な誤差制御を達成するため、双対問題の高次離散化を採用している。2次元および3次元のベンチマークにおいて、不自然な振動が著しく減少し、効果性指数が1に近い結果が得られた。

ABSTRACT

The efficient and reliable approximation of convection-dominated problems continues to remain a challenging task. To overcome the difficulties associated with the discretization of convection-dominated equations, stabilization techniques and a posteriori error control mechanisms with mesh adaptivity were developed and studied in the past. Nevertheless, the derivation of robust a posteriori error estimates for standard quantities and in computable norms is still an unresolved problem and under investigation. Here we combine the Dual Weighted Residual (DWR) method for goal-oriented error control with stabilized finite element methods. By a duality argument an error representation is derived on that an adaptive strategy is built. The key ingredient of this work is the application of a higher order discretization of the dual problem in order to make a robust error control for user-chosen quantities of interest feasible. By numerical experiments in 2D and 3D we illustrate that this interpretation of the DWR methodology is capable to resolve layers and sharp fronts with high accuracy and to further reduce spurious oscillations.

研究の動機と目的

• 標準の後験的誤差推定法がしばしば失敗する、対流優勢問題における信頼性と効率性に優れた誤差制御の課題に取り組む。
• 安定化有限要素法における既存の後験的誤差推定法の限界である、ロバストでない、または計算不能な誤差ノルムの問題を克服する。
• グローバルノルムではなく、ユーザー定義の目的関数に焦点を当てた目的関数指向の適応戦略を開発する。
• 双対問題の近似を改善することで、層や鋭いフロントを的確に捉えるロバスト性と正確性を確保する。
• 複雑な速度場を有する2次元および3次元の対流優勢移流問題において、本手法の有効性を示す。

提案手法

• 安定化有限要素解の目的関数指向の後験的誤差表現を導出するために、Dual Weighted Residual (DWR) 法を適用する。
• 誤差推定子の精度を向上させるために、双対問題に高次有限要素空間を用いることで、低次双対近似の欠陥を回避する。
• プライマル問題に残差に基づく安定化（特にSUPG）を実装し、対流優勢領域における不自然な振動を抑制する。
• DWR誤差表現から、局所的残差を双対解の値で重み付けした計算可能な誤差指標を構築する。
• 誤差指標に基づいてメッシュの適応的な細分化および粗化を実行し、自由度を目的関数に最も影響を与える領域に集中させる。
• プライマル問題よりも高次要素を用いて双対問題を解くことで、特に層が豊富な領域でも信頼性とロバスト性に優れた誤差推定を確保する。

実験結果

リサーチクエスチョン

• RQ1 高次双対離散化は、対流優勢問題における安定化有限要素法の目的関数指向誤差推定のロバスト性と信頼性を向上させるか？
• RQ2 DWRに基づく適応戦略は、鋭い層やフロントを的確に解像しつつ、不自然な振動を最小限に抑えるか？
• RQ3 提案手法は、2次元および3次元の設定において、目的関数の効果性指数を1に近づけるか？
• RQ4 SUPG安定化と目的関数指向の適応戦略の組み合わせは、グローバル誤差制御に比べてより効率的なメッシュを生成するか？
• RQ5 変動する対流場および複雑な幾何形状を有する3次元問題において、本手法はどのように性能を発揮するか？

主な発見

• テストされたすべての目的関数量において、効果性指数が1に非常に近い値を達成しており、高精度な誤差制御が実現している。
• 特に境界層および内部層において、数値解の不自然な振動が顕著に減少しており、これは双対問題の近似精度向上に起因する。
• 2次元および3次元のテストケースにおいて、適応的メッシュ細分化が、流れ方向に沿った層に特に関連する領域に要素を集中させることに成功した。
• 点値および平均値の目的関数に対して、粗い初期メッシュでも、ターゲット領域における誤差が最小限に抑えられた高精度な近似が得られた。
• 変動する対流場を有する3次元問題においても、本手法はロバスト性を維持しており、複雑で現実的な輸送問題への応用可能性を示した。
• 高次双対離散化を用いることで、低次双対近似に比べて誤差推定がより正確になり、誤差表現におけるソース誤差が低減された。