[論文レビュー] Goal-oriented a posteriori error control for nonstationary convection-dominated transport problems

Kristina Schwegler, Marius Paul Bruchhäuser+1 more|arXiv (Cornell University)|Jan 25, 2016

Advanced Numerical Methods in Computational Mathematics参考文献 18被引用数 4

ひとこと要約

本稿では、非定常な対流優勢輸送問題に対して、安定化有限要素法と組み合わせたデュアル重み付き残差（DWR）法を用いた、目的関数指向の後験的誤差制御フレームワークを提示する。高次近似を用いた双対問題の近似と「双対化してから安定化する」手法を採用することで、誤差推定のロバスト性が向上し、不自然な振動が顕著に低減され、目的関数の誤差推定の効果的指標は1に非常に近い値を示す。

ABSTRACT

The numerical approximation of convection-dominated problems continues to remain subject of strong interest. Families of stabilization techniques for finite element methods were developed in the past. Adaptive techniques based on a posteriori error estimates offer potential for further improvements. However, there is still a lack in robust a posteriori error estimates in natural norms of the discretizations. Here we combine the dual weighted residual method for goal-oriented error control with stabilized finite element approximations. By a duality argument an error representation is derived on that a space-time adaptive approach is built. It differs from former works on the dual weighted residual method. Numerical experiments illustrate that our schemes are capable to resolve layers and sharp fronts with high accuracy and to further reduce spurious oscillations of approximations.

研究の動機と目的

• 安定化有限要素法を用いた対流優勢問題に適用する際、自然なノルムにおけるロバストな後験的誤差推定器の不足を解消すること。
• 急峻な層や波-front を有する領域での精度を向上させるとともに、不自然な振動を最小限に抑える適応戦略の開発。
• 対流優勢領域において極めて重要な役割を果たす小拡散係数 ε に対してロバストであること。
• 双対解の近似品質を向上させることで、誤差推定の信頼性を強化すること。
• 安定化離散化と組み合わせたDWR法を用いて、時間依存輸送問題に対する効率的で目的関数指向の適応型メッシュ生成を可能にすること。

提案手法

• 空間時間領域における目的関数指向誤差制御のための誤差表現を導出するために、デュアル重み付き残差（DWR）法を適用する。
• 「双対化してから安定化する」戦略を採用：双対問題の弱形式を導出した後、安定化有限要素法（SUPGおよびショックキャプチャリング）を用いて解く。
• 層領域における精度を向上させるために、双対解に対して高次有限要素近似（例：区分的2次または4次多項式）を用いる。
• プライマル問題および双対問題の両方に、ストリームラインアップウィンドPetrov–ガレルキン（SUPG）安定化と追加のショックキャプチャリング（SOLD）を組み合わせ、振動を制御する。
• 目的関数指向誤差推定器に従って、誤差寄与度が高い領域を特定し、空間時間適応型メッシュを生成する。
• 対流優勢問題の物理的挙動を反映するノルムにおける誤差推定を導出し、一貫性およびロバスト性を確保する。

実験結果

リサーチクエスチョン

• RQ1 安定化有限要素法を用いた非定常な対流優勢輸送問題に対して、目的関数指向の後験的誤差推定がロバストに実現可能か。
• RQ2 双対問題の近似方法（例：多項式次数、安定化手法）の選択が、誤差推定器の精度および信頼性にどのように影響するか。
• RQ3 DWR法は、対流優勢問題の数値解における不自然な振動をどの程度低減できるか。
• RQ4 双対解に対して高次近似を用いることで、標準的な補間手法に比べて誤差推定器の効果的指標が向上するか。
• RQ5 提案された適応戦略は、最小限の自由度で急峻なフロントおよび境界層をどの程度正確に解明できるか。

主な発見

• 提案手法は、目的関数の誤差推定において、効果的指標が1に非常に近い値を示しており、誤差推定器の高い信頼性を示している。
• 特に高次双対近似を用いた場合、数値解における不自然な振動が顕著に低減されている。
• 双対問題にショックキャプチャリング安定化を適用しても、明確な改善効果は認められず、SUPGのみで十分であることが示唆された。
• 双対解の高次近似（例：区分的2次または4次多項式）により、低次近似に比べて収束性が向上し、誤差低減効果も顕著に改善された。
• 本手法は、3次元および時間依存設定においても、層および急峻なフロントを高い精度で解明できた。
• DWR法と安定化有限要素法の組み合わせにより、対流優勢問題に対してロバストで目的関数指向の適応型メッシュ生成が可能となり、ナビエ–ストークス流れのような複雑な応用分野への強い応用可能性を示した。