Cohomology of Finite Element Stokes Complexes on Alfeld Splits

技術要約：Alfeld分割における有限要素ストークス複体のコホモロジー

問題設定
本論文は、$\mathbb{R}^3$ 上の多面体領域における混合境界条件を持つストークス複体に対する、適合的な有限要素離散化の構築と解析を扱う。ストークス複体は、関数空間と微分演算子の連鎖として以下のように定義される：
$$0 \to H^2_{\Gamma_0}(\Omega) \xrightarrow{\text{grad}} H^1_{\Gamma_0}(\Omega; \text{curl}) \xrightarrow{\text{curl}} H^1_{\Gamma_0}(\Omega)^3 \xrightarrow{\text{div}} L^2(\Omega) \to 0$$
ここで、$\Gamma_0$ と $\Gamma_1$ は境界 $\partial\Omega$ を分割する。著者らは、Alfeld分割（頂点を重心に接続して細分された四面体）上のメッシュを用いた、多項式空間の離散化に焦点を当てている。

有限要素外微分解析（FEEC）における重要な課題は、離散複体が連続複体のトポロジカルな性質、特にそのコホモロジーを保持することを保証することである。コホモロジー群は、先行する演算子の核（カーネル）を像（レンジ）で割ったものを特徴付ける。離散コホモロジーが連続コホモロジーと同型であることは、頑健な数値手法（例：非圧縮流のための圧力ロバストなスキームや、パラメータ依存の問題に対する前処理）にとって極めて重要である。de Rham 複体や特定の自明な領域上のストークス複体についてはコホモロジーの結果が存在するが、混合境界条件を持つ非自明な領域上の適合的なストークス部分複体のコホモロジーは、本研究以前は文献において完全には扱われていなかった。

手法
著者らは、混合境界条件とストークス複体の特定の正則性要件を扱うために、[26] で開発された手法を拡張し、トレース構造、バブル空間、および骨格複体（skeletal complexes）に基づく枠組みを採用している。

1. 完全多項式空間： 著者らはまず、Alfeld分割上の「完全な」有限要素ストークス複体（空間は、特定の連続性要件（スカラー空間については $C^2$、ベクトル空間については $C^1$）を持つ高次多項式 $p \ge 3$ で構成される）を解析する。彼らは、部分単体（頂点、辺、面）に対するトレース演算子を定義し、大域的な有限要素空間を分解するためのバブル空間を構築する。
2. 境界条件の取り扱い： 混合境界条件を考慮するために、トレース構造を修正することが、本手法における重要な貢献である。著者らは、標準的なトレース条件では、境界 $\Gamma_0$ の「平坦な」部分における微分の消失を自動的に強制できないことを特定した。彼らは、修正された頂点バブル空間（$B_{0,nn}$ および $B_{1,nn}$）を導入し、これによって平坦な境界面に対する特定の法線微分制約を強制し、離散空間が要求される境界条件（$H^k_{\Gamma_0}$）を満たすようにしている。
3. 最小複体： 著者らは、局所的な複体の完全性と大域的なコホモロジーの同型性を維持しつつ、多項式の次数を可能な限り低減できる「最小」の部分複体を追求している。彼らは、Guzmán-Neilan 要素（最低次の $H^1$ 適合速度空間）と、区分的定数関数を利用する。また、最低次の完全多項式を含みつつ、完全な連鎖を形成するために必要な制約を満たす、新しい低次元の $H^2$ 適合および $H^1(\text{curl})$ 適合空間を構築している。
4. コホモロジー解析： 離散空間を骨格成分（頂点ベース）とバブル成分の直和分解として確立することにより、著者らは、離散複体のコホモロジーが相対単体コーチェイン複体のコホモロジーと同型であることを示している。これにより、単体コホモロジーと領域の相対ベッチ数を関連付ける既知の結果を活用することが可能となる。
5. 射影： 著者らは、これらの最小複体に対して、Scott-Zhang 型の有界な可換コーチェイン射影を構築し、これらが境界条件を保持し、微分演算子と可換であることを保証している。

主な貢献と結果

1. コホモロジーの同型性（完全複体）： 著者らは、Alfeld 分割上の完全多項式空間（[19] で定義）によって形成される離散ストークス複体が、任意の有界強リプシッツ領域および混合境界条件に対して、連続ストークス複体と同型のコホモロジーを持つことを証明している。具体的には、$\dim H^{k,h}_{\Gamma_0} = \dim H^k_{\Gamma_0} = b_k(\Omega, \Gamma_0)$ である。これは、[19] の結果を非自明な領域および混合境界条件へと拡張するものである。

2. 最小複体の構築： 本論文は、新しい「最小」適合有限要素複体を構築している。

   • $L^2$ 適合空間は、区分的定数関数（$DG_0$）である。
   • $H^1$ 適合空間は、Guzmán-Neilan 要素である。
   • 完全な複体を完成させるために、新しい低次元の $H^2$ 適合および $H^1(\text{curl})$ 適合空間が構築されている。
   • これらの最小空間は、完全な局所複体を形成し、連続複体と同型の全域コホモロジーを持つことが示されている。

3. 有界な可換射影： 著者らは、これらの最小複体に対する局所的で有界な可換コーチェイン射影を構築している。これらの射影は、境界条件を保持し、勾配（gradient）、回転（curl）、および発散（divergence）の演算子と可換であり、数値スキームの安定性や前処理において極めて重要な特性である。

4. 自由度： 論文では、新しい最小空間のユニソルベントな自由度の集合を明示的に特徴付けている。これらは、Whitney 形式に類似しているが、ストークス複体のより高い正則性に適応させたものであり、頂点の値、一定次数の微分、および辺/面の積分を含む。

意義
本論文は以下の領域において重要性を主張している：

• 理論的完全性： 非自明な領域および混合境界条件における適合的な有限要素ストークス複体のコホモロジーを厳密に確立することで、文献の空白を埋めている。
• 頑健性： コホモロジーの同型性を証明することにより、これらの複体に基づく数値スキーム（非圧縮ナビエ・ストークス方程式や結合応力理論など）が、連続問題のトポロジカルな性質を継承することを保証している。これは、質量保存、圧力ロバスト性、および一様な安定性につながる。
• 効率性： 「最小」複体の構築は、構造的完全性（コホモロジー）を犠牲にすることなく、より低次で計算コストの低い有限要素法への道筋を提供する。
• 前処理： コホモロジーの特性化と有界な可換射影の存在は、ストークス複体から生じるパラメータ依存問題のための、頑健な前処理器を構築するために不可欠である。

本研究は [26] の枠組みに基づき、それをストークス複体および Alfeld 分割の特定の文脈へと拡張しており、流体力学や弾性論においてこれらの特定のトポロジカルな制約が不可欠となる将来的な応用への基礎を提供している。