[論文レビュー] Two-side a posteriori error estimates for the DWR method

Bernhard Endtmayer, Ulrich Langer+1 more|arXiv (Cornell University)|Nov 19, 2018

Advanced Numerical Methods in Computational Mathematics参考文献 12被引用数 5

ひとこと要約

本稿では、複数の目的関数 functional を持つ非線形 PDE における双対重み付き残差 (DWR) 法について、両側 a posteriori 評価を提示する。離散化誤差と非線形反復誤差の両方を含む。飽和仮定の下で、1に近い効果的係数を有する効率的誤差推定子を確立し、数値的検証により収束速度の向上と無視可能な剰余項の存在を示した。

ABSTRACT

In this work, we derive two-sided a posteriori error estimates for the dual-weighted residual (DWR) method. We consider both single and multiple goal functionals. Using a saturation assumption, we derive lower bounds yielding the efficiency of the error estimator. These results hold true for both nonlinear partial differential equations and nonlinear functionals of interest. Furthermore, the DWR method employed in this work accounts for balancing the discretization error with the nonlinear iteration error. We also perform careful studies of the remainder term that is usually neglected. Based on these theoretical investigations, several algorithms are designed. Our theoretical findings and algorithmic developments are substantiated with some numerical tests.

研究の動機と目的

• 非線形 PDE における目的関数に特化した適応型有限要素法の信頼性と効率性のある a posteriori 評価を開発すること。
• DWR 法を拡張し、離散化誤差と非線形反復誤差の両方を考慮することで、ニュートン法などの反復ソルバーの適応的停止を可能にすること。
• 標準的な DWR 評価子で無視されがちな剰余項を厳密に解析し、その高次項としての寄与が効率性を損なわないことを保証すること。
• 飽和仮定の下で、一般の非線形 PDE と非線形目的関数 functional に対して、誤差推定子の理論的下界（効率性）を確立すること。
• 離散化誤差と反復誤差をバランスさせる強固な適応アルゴリズムを設計・検証し、最適収束速度を達成すること。

提案手法

• 後退積分を避けるために、階層的拡張とユニオン分割局在化を用いた双対重み付き残差法により、a posteriori 誤差推定子を導出する。
• 非線形反復誤差を誤差推定子の計算可能な成分として組み込む修正された DWR フレームワークを導入し、ニュートン法の適応的停止を可能にする。
• 飽和仮定を用いて誤差推定子の下界を導出し、一般の非線形 PDE と目的関数 functional に対してその効率性を証明する。
• 誤差推定子における剰余項を解析し、それが高次項であることを示し、実用上無視可能であることを示し、標準的定式化における省略が正当化されることを示す。
• 効果的係数を用いて推定子の品質を監視しながら、メッシュの適応的グレーディングと反復誤差の制御を同時に実行する適応アルゴリズムを設計する。
• 誤差推定子の変分的定式化を用い、p-Laplacian モデルに複数の目的関数 functional を適用した数値実験により、手法を検証する。

実験結果

リサーチクエスチョン

• RQ1 非線形 PDE と複数の目的関数 functional を持つ DWR 法について、離散化誤差と反復誤差の両方を含む二重側 a posteriori 評価を導出できるか？
• RQ2 誤差推定子が効率的である、すなわち定数が1に近い下界を有するための条件は何か？
• RQ3 DWR 誤差推定子における剰余項の寄与はどの程度顕著であり、厳密に評価可能か、あるいは無視可能か？
• RQ4 DWR 法を、離散化誤差と非線形反復誤差をバランスさせる形に拡張し、反復ソルバーの適応的停止を可能にすることができるか？
• RQ5 提案された推定子と適応アルゴリズムは、実際の非線形的・非対称的問題においても最適収束速度を達成するのか？

主な発見

• 数値実験において、提案された二重側誤差推定子は効果的係数が1に極めて近い値を達成し、その信頼性と効率性を確認した。
• DWR 推定子における剰余項が高次項であることが示され、標準的定式化における省略が推定子の品質を損なわないことを裏付けた。
• 飽和仮定により下界が導出可能であり、単一および複数の目的関数 functional に対して誤差推定子の効率性が証明された。
• 完全な誤差推定子に基づく適応的グレーディングにより収束速度が向上し、結合目的関数 $J_{\mathfrak{E}}$ に対して観測された収束次数は $\mathcal{O}(\text{DOFs}^{-1})$ であった。
• 反復誤差を推定子に組み込むことで、ニュートン法のアフィン不変停止基準が可能となり、ロバスト性と効率性が向上した。
• 数値結果から、プライマル部のみでは信頼性のある誤差推定が不十分であり、随伴部と反復部を含む完全な推定子が正確な効果的係数を得るために不可欠であることが確認された。