[論文レビュー] Discrete Hessian complexes in three dimensions

Long Chen, Xuehai Huang|arXiv (Cornell University)|Dec 20, 2020

Black Holes and Theoretical Physics参考文献 37被引用数 5

ひとこと要約

本稿は、多項式テンソル空間の分解を用いて、3次元の単体メッシュ上に適合型仮想要素ヘッシアン複体を構築し、線形化されたアインシュタイン＝ビアンキ系の最適順序離散化を可能にする。この手法は、H²-適合仮想要素とH(div)-適合有限要素を、 Koszul型の分解により組み合わせ、適切に選ばれた自由度と射影に基づく作用素を通じて安定性と計算可能性を確保し、エネルギーノルムにおいて最適収束率を達成する。

ABSTRACT

A family of conforming virtual element Hessian complexes on tetrahedral meshes are constructed based on decompositions of polynomial tensor spaces. They are applied to discretize the linearized time-independent Einstein-Bianchi system with optimal order convergence.

研究の動機と目的

• 3次元の単体メッシュ上に、最小限の自由度で安定かつ適合型の離散ヘッシアン複体を構築すること。
• 仮想要素法を用いて、線形化された時間に依存しないアインシュタイン＝ビアンキ系の最適順序有限要素離散化を可能にすること。
• 従来の有限要素ヘッシアン複体が直面する高次多項式の複雑さを、仮想要素技術を活用することで克服すること。
• 多項式テンソル空間の新規分解を通じて、離散複体の適切な定義と計算可能性を保証すること。
• H(curl)-適合性やL²射影の計算可能性といった重要な構造を保持しつつ、最適収束を達成するフレームワークを確立すること。

提案手法

• 単体上での多項式テンソル空間に対して、Koszul型の分解 Pk(K; S) = ∇²Pk+2(K) ⊕sym(Pk−1(K; T) × x) を用いて離散ヘッシアン複体を構築する。
• 仮想要素空間 Σ(K) = ∇²W(K) ⊕sym(V(K) × x) を定義する。ここで、W(K) はH²-適合仮想要素空間、V(K) = Pk−1(K; T) はH(div)-適合有限要素空間である。
• 射影に基づくアプローチを実装する：QKkσ は内部モーメントから計算可能であり、eQKk+2σ は ΠKσ から再構成されることで、L²射影の計算可能性が保証される。
• 離散双一次形式の定義に、bK と cK の双一次形式を用い、離散インフラス条件を通じて安定性を確保する。
• 変分形式 (91)–(92) に基づく適合型混合仮想要素法を適用し、局所形式のグローバルアセンブリにより解を計算する。
• 一貫性および安定性のある混合有限要素法の定式化を採用し、離散インフラス条件と射影に基づく一貫性推定に依存する。

実験結果

リサーチクエスチョン

• RQ1 仮想要素を用いて、既存の有限要素複体よりも少ない自由度で、3次元の単体メッシュ上に適合型離散ヘッシアン複体を構築できるか？
• RQ2 完全なトレースの特徴付けが得られない多面体領域において、非同次境界条件を伴う局所的H²問題の適切な定義をどのように保証できるか？
• RQ3 どの多項式テンソル空間の分解が、ヘッシアン複体のための安定的かつ適合型仮想要素空間の構築を可能にするか？
• RQ4 この仮想要素フレームワークを用いて、線形化されたアインシュタイン＝ビアンキ系に対して最適収束率を達成できるか？
• RQ5 形状関数の明示的知識がなくても、仮想要素空間のL²射影を多項式空間上に計算可能にするにはどうすればよいか？

主な発見

• 正則性の仮定の下で、提案された仮想要素ヘッシアン複体は、ヘッシアンのH(curl)-ノルムおよびポテンシャルのH²-ノルムにおいて、最適収束次数 hk を達成する。
• 離散インフラス条件は、新規な安定性の議論により確立され、混合定式化の適切な定義が保証される。
• 一貫性推定式 (96)–(98) により、射影に基づく作用素が、それぞれの関連ノルムにおいて hk および hk+1 のオーダーの誤差を導入することが示された。
• L²射影が、内部モーメントと再構成作用素のみを用いて計算可能であることが保証される。
• 収束速度 ∥E − Eh∥H(curl) + ∥σ − σh∥2 ≲ hk(|E|k+1 + ∥σ∥k+2 + |f|k) が証明され、滑らかな解に対して最適順序が確認された。
• 従来の有限要素複体が要請する高次多項式基底関数（例：P9）の必要性を回避し、自由度を顕著に削減した。