[論文レビュー] Adaptive finite element methods for partial differential equations

Rolf Rannacher|ArXiv.org|May 1, 2003

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

ひとこと要約

本稿では、偏微分方程式の有限要素解における後験的誤差推定のためのデュアル・ウェイトド・リジッド（DWR）法を導入し、目的関数に向けた適応的メッシュ・リファインメントを可能にする。残差にデュアル問題の解を用いて重みづけすることで、特定の出力量における誤差を最小化する経済的なメッシュを構築する。流体力学的応用において、抗力計算、抗力低減、安定性解析の各分野で高い効率性と正確性を示した。

ABSTRACT

The numerical simulation of complex physical processes requires the use of economical discrete models. This lecture presents a general paradigm of deriving a posteriori error estimates for the Galerkin finite element approximation of nonlinear problems. Employing duality techniques as used in optimal control theory the error in the target quantities is estimated in terms of weighted `primal' and `dual' residuals. On the basis of the resulting local error indicators economical meshes can be constructed which are tailored to the particular goal of the computation. The performance of this {\it Dual Weighted Residual Method} is illustrated for a model situation in computational fluid mechanics: the computation of the drag of a body in a viscous flow, the drag minimization by boundary control and the investigation of the optimal solution's stability.

研究の動機と目的

• 特定の出力関数への適合したガレルキン有限要素法における後験的誤差推定の一般枠組みの構築を目的とする。
• 目的関数量の所定の精度を満たすために計算コストを最小化するように、効率的かつ目的関数に向けたメッシュ適応を可能にする。
• 抗力計算、抗力低減のための境界制御、最適流れの安定性解析を含む、複雑な流体力学的問題へのこの手法の適用を目的とする。
• 適応的リファインメントを通じて最小限の自由度で高い精度を達成するDWR法の有効性を示す。

提案手法

• 最適制御理論からの双対性技術を用いて、後験的誤差表現を導出し、出力関数の誤差を、プライマルおよびデュアルの残差の重み付き和として表現する。
• 線形化されたデュアル問題の解を用いてデュアル・ウェイトド・リジッド（DWR）推定子を導入し、デュアル解が局所的誤差指標の重みを提供する。
• プライマル残差（離散解からのもの）とデュアル重み（デュアル解からのもの）の積として局所的誤差指標を構築し、局所的メッシュリファインメントを可能にする。
• デュアル解の誤差を近似するために、パッチ単位の高次補間を適用し、要素間での誤差寄与をバランスさせる。
• 誤差バランスを用いて、局所的誤差寄与が最大となる要素をリファインし、メッシュ全体にわたる均一な誤差分布を確保する。
• 必要に応じて、誤差推定子の安定性と収束性を保証するため、出力関数 Functional の正則化を実装する。

実験結果

リサーチクエスチョン

• RQ1 非線形ガレルキン有限要素近似に対して、全解の誤差ではなく、特定の出力関数に適合した後験的誤差推定をどのように導出できるか。
• RQ2 デュアル・ウェイトド・リジッドに基づく適応的メッシュリファインメントは、均一またはヒューリスティックなリファインメント戦略と比較して、より少ない自由度でより高い精度を達成できるか。
• RQ3 DWR法は、境界制御下での抗力計算や抗力低減という実用的流体力学的問題において、どのように性能を発揮するか。
• RQ4 境界制御によって得られた最適化された流れ解の安定性はいかほどか。また、DWR法はこのような解の解析を支援できるか。

主な発見

• 抗力計算において、DWR法は2,244の自由度で抗力係数の誤差1%を達成し、均一リファインメントを上回る効率性を示した。
• Re = 40における抗力低減において、適応的リファインメントは4,264の自由度で抗力係数3.16723を達成したが、均一リファインメントでは164,928の自由度を要して3.11800にまで到達した。
• 抗力低減のDWR適応メッシュは、均一リファインメントと比較して自由度を75%削減したが、同等の精度を維持した。
• 抗力低減の最適解は不安定性の端に位置しており、臨界固有値の実部が0に近づく傾向にあり、ホップ分岐の可能性を示唆した。
• 流れ線の可視化と固有値解析により、最適化された流れにおける不安定性の発生が確認され、臨界固有関数は分岐点近くで振動的挙動を示した。
• DWR法は、正確かつ効率的にデュアル解を計算し、信頼性の高い誤差指標を提供し、すべてのテストケースにおいて効果的なメッシュ適応を可能にした。