[論文レビュー] Virtual Element Method for Second Order Elliptic Eigenvalue Problems

Francesca Gardini, Giuseppe Vacca|arXiv (Cornell University)|Oct 12, 2016

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

ひとこと要約

本稿では、多角形／多面体メッシュ上での2階楕円型固有値問題（ラプラス固有値問題を含む）を解くための仮想単位要素法（VEM）を提示する。本手法は、射影に基づく安定化を用いて安定かつ一貫性のある離散双線形形式を構築することで、固有値および固有関数の両方において最適収束率を達成する。数値実験により、さまざまなメッシュタイプおよび特異解を含む理論的予測が確認された。

ABSTRACT

We introduce the Virtual Element Method (VEM) for elliptic eigenvalue problems. The main result of the paper states that VEM provides an optimal order approximation of the eigenmodes. A wide set of numerical tests confirm the theoretical analysis.

研究の動機と目的

• 2階楕円型固有値問題（ラプラス固有値問題を含む）のための仮想単位要素法（VEM）を開発すること。
• VEMフレームワークにおいて、固有値および固有関数の両方の最適収束率を確立すること。
• 固有値問題における安定化および非安定化の両定式化を分析すること。
• 多角形メッシュ（非凸領域および不連続係数を含む）を用いた広範な数値的検証を通じて理論的結果を検証すること。
• VEMの適用範囲を正則性が低い解および特異解を有する固有値問題へ拡張すること。

提案手法

• 本手法は、一般の多角形／多面体要素上で定義された仮想単位要素空間を採用し、非多項式関数の複雑な積分を避けるために任意の多角形メッシュを許容する。
• エネルギー（勾配－勾配）形式とL²内積の両方の離散双線形形式を、射影に基づく近似を用いて構築する。
• 離散エネルギー形式は、多項式空間への勾配の射影を用い、一貫性および安定性を保証する。
• エネルギー形式に安定化項を導入し、スカラーおよび非スカラーの安定化オプションを備えることで、最適収束を保証する。
• 弱ガレルキン形式として定式化され、固有値問題は離散空間上での一般化固有値問題として解かれる。
• 局所的射影および多項式の一貫性を保つ度数自由度の使用に依存した実装であり、非多項式基底関数の明示的計算を回避する。

実験結果

リサーチクエスチョン

• RQ1 仮想単位要素法は、2階楕円型固有値問題において、固有値および固有関数の両方で最適収束率を達成できるか？
• RQ2 安定化パラメータの選択が、VEM近似のロバスト性および精度にどのように影響するか？
• RQ3 非凸領域における特異固有関数を有する問題に対しても、VEMは最適収束を維持するか？
• RQ4 PDEにおける不連続な拡散係数に対して、この手法はどのように動作するか？
• RQ5 VEMは、固有値問題における混合型定式化および適応的グリッド・リファインメントに拡張可能か？

主な発見

• VEMは、局所的仮想単位要素空間の多項式次数 $k$ と一致する収束次数を有する、固有値および固有関数の両方で最適収束率を達成する。
• L字型領域における最初の固有値では、正則性の低下により収束率が最適でない（$r=2/3 - \varepsilon$）が、解析的であるとされる3番目の固有関数は、完全な最適収束率を達成する。
• 数値的テストにより、安定化パラメータ $\tau_P$ に対して本手法はロバストであることが示され、$\tau_P \sim h^{-2}$ の場合にのみ失敗が観察された。
• 不連続な拡散係数（$\delta = 0.50, 0.10, 0.01, 10^{-8}$）を有する問題に対しても、最適収束が維持され、係数のジャンプに対してもロバストであることが確認された。
• テスト7.3および7.4における $k=1,2$ の収束プロットは、予測された収束率を確認し、理論的誤差推定の妥当性を検証した。
• 本手法はボロノイメッシュおよび一般の多角形メッシュでも良好に動作し、複雑な幾何形状における柔軟性およびロバスト性を示した。