[論文レビュー] Cut Finite Element Methods for Partial Differential Equations on Embedded Manifolds of Arbitrary Codimensions

Erik Burman, Peter Hansbo+2 more|arXiv (Cornell University)|Oct 5, 2016

Numerical methods in engineering被引用数 4

ひとこと要約

本稿は、R^k に埋め込まれた任意の余次元の多様体上でラプラス＝ベルトラミ方程式を解くために、安定化されたカット有限要素法の統一的理論的枠組みを提示する。3つの新規安定化形式を導入し、特に表面法線勾配ペナルティ項を新たに提案することで、メッシュと幾何の整合性に依存しない最適収束速度と堅牢な条件数の上限を達成した。数値的にR^3における曲線および曲面に対して検証された。

ABSTRACT

We develop a theoretical framework for the analysis of stabilized cut finite element methods for the Laplace-Beltrami operator on a manifold embedded in $\\mathbb{R}^d$ of arbitrary codimension. The method is based on using continuous piecewise polynomials on a background mesh in the embedding space for approximation together with a stabilizing form that ensures that the resulting problem is stable. The discrete manifold is represented using a triangulation which does not match the background mesh and does not need to be shape-regular, which includes level set descriptions of codimension one manifolds and the non-matching embedding of independently triangulated manifolds as special cases. We identify abstract key assumptions on the stabilizing form which allow us to prove a bound on the condition number of the stiffness matrix and optimal order a priori estimates. The key assumptions are verified for three different realizations of the stabilizing form including a novel stabilization approach based on penalizing the surface normal gradient on the background mesh. Finally, we present numerical results illustrating our results for a curve and a surface embedded in $\\mathbb{R}^3$.

研究の動機と目的

• R^k に埋め込まれた任意の余次元の多様体上で安定化されたカット有限要素法を分析する一般理論的枠組みを構築すること。
• 任意のメッシュ-幾何の交差を許容しても、離散問題における最適な事前誤差推定と条件数の上限を確立すること。
• 新規の表面法線勾配安定化を含む、3つの異なる安定化形式への適用性を検証すること。
• 標準的な最近接点写像が利用できない高余次元における幾何的近似誤差推定と補間誤差の上限を拡張すること。
• 曲線および曲面におけるR^3内の数値実験を通じて、メッシュの配置に依存しない堅牢性と収束性を示すこと。

提案手法

• R^k における背景メッシュ上に定義された連続的な区分多項式有限要素を用い、それが任意にメッシュを貫通する離散的多様体に制限される。
• メッシュと要素の形状が一致する必要がなく、形状正則な三角形分割を要しない安定化形式を導入し、安定性と条件数の制御を保証する。
• 3つの安定化機構を採用する：(1) アクティブメッシュ面をはさんで法線勾配のジャンプをペナルティ化、(2) アクティブ要素上の全勾配をペナルティ化、(3) 背景メッシュ上での表面法線勾配ペナルティを新規に導入。
• 条件数の有界性と最適収束速度を保証する安定化形式に関する抽象的条件を導出する。
• 最近接点写像を明示的に使用しないで、近似と正則性の仮定のみに依存して離散的多様体の幾何的近似誤差推定を確立する。
• 標準的な余次元1を超える設定において、補間誤差推定を一般化し、先行研究を余次元1を超えて拡張する。

実験結果

リサーチクエスチョン

• RQ1 任意の余次元の多様体上でのカット有限要素法に対して、統一的理論的枠組みを構築できるか？
• RQ2 条件数の上限と最適誤差推定を達成するために必要な安定化形式に関する主要仮定が、複数の安定化戦略においても成立するか？
• RQ3 背景メッシュ上での表面法線勾配に基づく新規安定化は、全勾配安定化よりも優れた一貫性と安定性を達成できるか？
• RQ4 最近接点写像が明示的に利用できない高余次元において、幾何的近似誤差はどのように振る舞うか？
• RQ5 条件数はメッシュ-幾何の整合性にどのように依存するか？安定化により、すべての配置において堅牢性を確保できるか？

主な発見

• 提案された枠組みは、新規の表面法線勾配安定化を含む3つの安定化形式すべてに対して、最適なオーダーの事前誤差推定と有界な条件数を保証する。
• 剛性行列の条件数はメッシュサイズ h に対して O(h^{-2}) のスケーリングに従い、数値結果により、さまざまなメッシュ配置でもこの上限が確認された。
• 新規の表面法線勾配安定化は、全勾配安定化よりも著しく小さい一貫性誤差を達成し、より柔軟なパrameter選択と高次近似を可能にする。
• R^3における曲面および曲線に対する数値実験では、安定化が適用された場合、特にペナルティパrameter τが十分に大きい場合には、条件数が堅牢に保たれた。
• 安定化なしの手法では、メッシュ-幾何の整合性の関数として、強い不規則なピークが条件数に現れたが、安定化されたバージョンはすべての配置で良好に条件付けられていた。
• 平均スケーリング条件数 h²κ(𝒜) は、すべてのテストケースおよび安定化形式において有界であり、h→0 の極限で定数に収束することが確認され、理論的予測が裏付けられた。