Uniform inf–sup stability of quartic and quintic Scott–Vogelius elements on Freudenthal meshes: a protected raw edge-star lifting

技術要約：Freudenthalメッシュ上の4次および5次Scott–Vogelius要素の一様inf–sup安定性

問題設定
本論文は、一様な3次元Freudenthal三角形分割におけるストークス方程式に対するScott–Vogelius有限要素の、一様なinf–sup安定性を扱っている。具体的には、多項式次数 $k=4$ および $k=5$ を対象としている。Zhangは以前、$k \ge 6$ におけるメッシュ一様安定性を確立したが、その構成は、エッジ・トレースを構築した後に要素ごとの発散平均を修復するために、6次のバブルを挿入することに依存していた。特定された決定的な欠落は、エッジ・トレース関数自体がすでに次数 $k \ge 4$ を持っているにもかかわらず、既存の構成では、高次で使用される逐次的な「借用」ロジックに頼ることなく、これらのトレースを扱うための真に局所的かつ次数に依存しない持ち上げ（lifting）メカニズムが欠けていたことである。不足している要素は、平均修復のための独立した戦略と組み合わせることで、$k=4, 5$ においても安定性を証明できる、局所的なエッジ段階の持ち上げである。

手法
著者は、厳密な有理数演算と整数行列の恒等式に基づいたコンピュータ支援証明戦略を用いて、「保護された生の（protected raw）エッジ・スター持ち上げ」を構築している。手法は以下の段階を経て進行する：

1. 幾何学的列挙と分類：
   著者は、Freudenthalメッシュにおけるすべての可能なエッジ構成を分類している。接する四面体と立方体を分析することにより、正確に7つの幾何学的クラスと37の方向/境界構成を特定した。この列挙は、1辺あたり $N \ge 3$ 個の立方体を持つメッシュにおけるすべてのエッジタイプ（内部、境界、面対角線、体対角線）を網羅している。

2. 厳密な局所的持ち上げの構築：
   各37の構成および次数 $k=4, 5$ に対して、著者は局所的な線形写像を構築する。

• ソース空間： 連続性の制約（特異エッジおよび境界の適合関係）に従う、対象エッジ上の発散トレースの空間。
• ターゲット空間： エッジに接する立方体の和集合（「スター」）上の速度係数の空間。
• 保護（Protection）： この写像は、ターゲットのエッジ・トレースを再現しつつ、すべての非ターゲットのエッジ・トレースおよび頂点値がゼロであることを保証するように設計されている。
• 検証： これらの写像は、非公式な「自由度の借用」によって導かれたものではなく、特定の代数的恒等式（例：$HB=0, DB=s_0 I$）を満たす整数行列 $B, D, C$ を用いて、厳密な整数行列によって証明されている。ソース・トレース空間は、整数行列 $H$ の核として特徴付けられる。

3. ノルム推定：
   厳密な有理数Bernstein質量および剛性計算を用いて、著者は持ち上げ写像の $H^1$ ノルムを計算している。彼らは、すべての構成に対して一様な参照パッチ・バウンド $C_{ref} < 385$ を確立している。このバウンドは、アフィン・スケーリング特性により、メッシュサイズ $h$ に依存しない。

4. 組み立てとオーバーラップ制御：
   グローバルな解を組み立てるために、著者はエッジの方向と端点の座標を3を法として用いた彩色スキームを定義している。これにより、同じ色のパッチが互いに素となる189色彩色が得られる。最大サポート・オーバーラップは19に制限される。組み立ての依存深度はゼロであり、これはすべてのエッジ持ち上げが逐次的依存関係なしに同時に適用可能であることを意味する。

5. 平均修復：
   生の（raw）エッジ持ち上げは、必ずしも個々の四面体における発散のゼロ平均を保持しない。著者は、保護された生の（protected raw）エッジ・スター・レマと、4次（$P_4$）マクロパッチ演算子を組み合わせた手法を導入している。この演算子は、2つの立方体パッチを用いて要素ごとの発散平均を修復し、エッジ・トレースを変更することなく、最終的な場が各四面体上でゼロ平均であることを保証する。この4次修復は、Zhangの元の構成で使用されていた6次のバブルを、平均補正段階において置き換えるものである。一方で、新しいエッジ段階の持ち上げは、$k=4, 5$ におけるトレースの一致を局所的に処理する。

主要な貢献

• $k=4, 5$ のギャップの解決： 本論文は、Freudenthalメッシュ上のScott–Vogelius要素について、多項式次数4および5に対する一様なinf–sup安定性の最初の証明を提供し、すべての $k \ge 4$ に対する結果を完結させた。
• 保護されたエッジ・スター持ち上げ： 局所的で次数に依存せず、厳密な整数行列によって証明された「保護された」持ち上げメカニズムの導入。これにより、Zhangのエッジ段階の次数要求に関する観察を、逐次的な「1つの自由度を借りる」ロジックに頼ることなく、低次に対する厳密な証明へと転換した。
• 完全な幾何学的センサス： Freudenthalパターンに対して網羅的であることを検証された、37のエッジ構成と117の境界装飾付きソース・エンベロープの厳密な列挙。
• コンピュータ支援検証： 構成は、厳密な有理数証明を含む再現性パッケージに依存している。著者は、代数的恒等式、ノルム・バウンド、および組み立て特性を検証するために、独立したPython検証器およびLean形式化（Mathlibを使用）を利用している。Leanコードには、コアとなる代数的チェックにおいて sorry や未証明の公理は含まれていない。

結果

• 定理 1.1： 固定された整数 $k \ge 4$ に対して、メッシュサイズ $h$ に依存しない定数 $\beta_k > 0$ が存在し、Freudenthalメッシュ上のScott–Vogeliusペア $(V^k_h, Q^k_h)$ に対してinf–sup条件が成立する。
• 一様バウンド： 構築された右逆写像 $R_h: Q^k_h \to V^k_h$ は、$|R_h q_h|_1 \le C_k \|q_h\|_0$ を満たす。ここで $C_k$ は $h$ に依存しない。
• 具体的な定数： エッジ持ち上げの参照パッチ・バウンドは $C_{ref} < 385$ である。グローバルな組み立ては、最大オーバーラップ19を伴う189彩色を用いる。

意義と主張
本論文は、「真に局所的で次数に依存しない持ち上げ」を提供し、エッジ段階の次数要求に関する観察を、低次に対する厳密な証明へと変容させることを主張している。エッジ段階における「1つの自由度を借りる」という逐次的なロジックを、エッジ段階のための有限な厳密な行列に置き換え、また平均修復のために（Zhangの仕事で使用された）6次のバブルの代わりに4次マクロパッチ演算子を使用することで、著者は$k=4, 5$におけるエッジ段階でのアドホックな次数上昇の必要性を排除している。本研究は、$k=4, 5$ の安定性が、前述の $k \ge 6$ の解析で必要とされた次数6のバブル補正に頼ることなく、それぞれ4次および5次の要素のみを用いて達成可能であることを示している。厳密な有理数証明とLeanによる形式化の使用は、浮動小数点誤差なしに、安定性定数と構造的特性が数学的に検証されていることを保証している。この結果は、Freudenthalメッシュ上のすべての $k \ge 4$ に対する安定性理論を統一するものである。