[論文レビュー] Revisiting finite element accuracy by a mixed functional-probabilistic approach.

Joël Chaskalovic, Franck Assous|arXiv (Cornell University)|Mar 26, 2018

Probabilistic and Robust Engineering Design参考文献 12被引用数 3

ひとこと要約

本稿では、$P_k$ と $P_m$ 要素（$k < m$）間の相対的精度を確率変数としてモデル化することにより、有限要素法の精度を混合的機能的・確率的アプローチで評価する手法を提案する。ブラームブ・ヒルベルト補題の幾何的解釈を用いて、あるメッシュサイズ $h$ において一方の要素が他方よりもより正確である確率がどのようになるかを示す2つの確率分布を導出する。$m-k \to \infty$ の漸近的挙動も分析している。主な貢献は、$h$ に応じた要素の優位性を予測する確率的フレームワークを提供し、有限要素選択に関する新たな知見をもたらす点である。

ABSTRACT

The aim of this paper is to provide a new perspective on finite element accuracy. Starting from a geometrical reading of the Bramble-Hilbert lemma, we derive two probability distributions that estimates the relative accuracy, considered as a random variable, between two finite elements $P_k$ and $P_m$, ($k < m$). We establish mathematical properties of these laws, particularly studying their asymptotic relation when the difference $m-k$ goes to infinity. Then, we get new insights which, among others, show that $P_k$ or $P_m$ is more likely accurate than the other, depending on a the value of the mesh size $h$.

研究の動機と目的

• 決定的誤差境界を越えて、確率的視点から有限要素法の精度を再評価すること。
• 2つの有限要素 $P_k$ と $P_m$（$k < m$）間の相対的精度を確率変数としてモデル化すること。
• 一方の要素が他方よりも正確である確率を推定する確率分布を導出すること。
• 多項式次数の差 $m-k$ が増加する際のこれらの分布の漸近的挙動を分析すること。
• $P_k$ や $P_m$ がより正確である可能性が高いメッシュサイズ $h$ の条件を特定すること。

提案手法

• 有限要素近似誤差を関数空間内でモデル化するため、ブラームブ・ヒルベルト補題の幾何的解釈を導出する。
• $P_k$ と $P_m$ 要素間の相対的精度を推定する2つの確率分布を構築する。
• 関数解析を用いて、近似ノルムに基づく相対的精度を表す確率変数を定義する。
• 分布の漸近的性質を $m-k \to \infty$ の下で分析し、収束の傾向を明らかにする。
• モーメントや尾部挙動を含む導出法則の数学的性質を確立し、確率的推論を支援する。
• 得られた確率的モデルをメッシュサイズ $h$ と関連づけ、一方の要素が他方を上回る条件を特定する。

実験結果

リサーチクエスチョン

• RQ1 有限要素 $P_k$ と $P_m$ 間の相対的精度を確率変数としてどのようにモデル化できるか？
• RQ2 $P_k$ が $P_m$ よりも正確である確率を記述する確率分布は何か？
• RQ3 メッシュサイズ $h$ は、精度の観点から一方の要素が他方を上回る確率にどのように影響するか？
• RQ4 多項式次数の差 $m-k \to \infty$ の下で、精度の確率分布の漸近的挙動はいかなるものか？
• RQ5 どのような $h$ の条件下で $P_k$ が $P_m$ よりも正確である可能性が高くなり、逆に $P_m$ が優位になるか？

主な発見

• 本稿では、$P_k$ と $P_m$ 有限要素間の相対的精度を確率変数としてモデル化する2つの確率分布を導出している。
• これらの分布は、$P_k$ や $P_m$ がより正確である確率が、メッシュサイズ $h$ に強く依存することを示している。
• 多項式次数の差 $m-k$ が大きくなるにつれて、導出された確率法則は特定の漸近的収束パターンを示す。
• 小さな $h$ の場合、$P_m$ がより正確である可能性が高くなるが、大きな $h$ の場合、問題の正則性に応じて $P_k$ が優位になる可能性がある。
• 確率的フレームワークにより、従来の誤差推定を超えて、$h$ に基づく有限要素選択の新たな基準が提供される。