Deflation-based preconditioning for immersed finite element methods and immersogeometric analysis

この論文「Deflation-based preconditioning for immersed finite element methods and immersogeometric analysis（浸り込み有限要素法および浸り込み幾何解析のためのデフレーションベースの事前条件付け）」は、CAD 幾何形状のトリミング（切り取り）によって生じる浸り込み有限要素法（Immersed FEM）および浸り込み幾何解析（Immersogeometric Analysis, IGA）における極端な数値的不安定性（条件数の悪化）の問題に焦点を当て、既存の事前条件付け手法の限界を指摘し、新しいデフレーション（Deflation）に基づくロバストな事前条件付け手法を提案するものです。

以下に、問題定義、手法、主要な貢献、結果、および意義について詳細な技術的サマリーを記述します。

---

1. 問題定義

浸り込み有限要素法では、複雑な CAD 形状を単純な背景メッシュに埋め込むことでメッシュ生成の負担を軽減しますが、物理領域の境界が背景メッシュの要素を切断（トリミング）することで、**極端に小さな要素（cut elements）**が発生します。これにより、以下の問題が生じます。

• 極端な条件数の悪化: 小さな要素にのみサポートされる基底関数が存在すると、システム行列（剛性行列や質量行列）の最小固有値が極めて小さくなり、条件数が爆発的に増加します。
• 既存手法の限界:
  • 対角スケーリング（Jacobi 事前条件付け）: 基底関数のスケーリング不均衡は解消できますが、近似的な線形依存性（near linear dependencies）を除去できず、条件数が改善されない場合があります。
  • SIPIC（Symmetric Incomplete Permuted Inverse Cholesky）: 対角スケーリングと局所的な直交化を組み合わせますが、特定の幾何配置（カウンター例）において近似的な線形依存性を検出できず、失敗することがあります。
  • Schwarz 事前条件付け: ブロック選択戦略に依存しており、最適なブロック分割が確立されておらず、特定のトリミング構成ではロバスト性が保証されません。
  • 安定化手法（$\alpha$-stabilization など）: 変分形式自体を変更するため、元の解の精度を損なう可能性があります。

2. 手法と分析

著者らは、まず既存の手法の限界を理論的および数値的に分析し、その上で新しいアプローチを提案しました。

2.1 既存手法の限界に関する分析

• 対角スケーリングの解析: 1 次元および 2 次元の様々なトリミング構成（スリバーカット、コーナーカット、ミドルカット）に対して、条件数のスケーリング関係を厳密に導出しました。
  • ラグランジュ基底: 対角スケーリング後も条件数が体積分数 $\eta$ に対して $\eta^{-2(p-1)}$ のように悪化します。
  • B-スプライン基底: 1 次元や特定の 2 次元構成では対角スケーリングのみで条件数が $\eta$ に依存しなくなる（ロバストになる）ことが示されましたが、「ミドルカット」のような特定の 2 次元構成では、基底関数の連続性やノット位置に依存して条件数が再び悪化することが判明しました。
• カウンター例の提示: SIPIC および Schwarz 事前条件付けが失敗する具体的な幾何形状（回転した正方形、複数のリッジを持つ構造など）を構築し、これらの手法が「近似的な線形依存性」を完全に検出・除去できないことを実証しました。

2.2 提案手法：デフレーションベースの事前条件付け

線形代数におけるデフレーションの概念を浸り込み FEM に適用し、**「小さな要素にのみサポートされる基底関数」**が引き起こす小さな固有値を直接除去する手法を提案しました。

• 基本原理:
  • 問題となる基底関数（トリミング領域にのみサポートされるもの）の集合を特定し、それらを張る部分空間（デフレーション部分空間）を定義します。
  • 射影行列 $P = I - AZ(Z^T AZ)^{-1}Z^T$ を用いて、システム行列から小さな固有値を「デフレーション（除去）」します。ここで $Z$ は、問題となる基底関数に対応する単位ベクトルから構成されます。
  • これにより、反復解法（共役勾配法）の収束を阻害する小さな固有値をゼロにシフトさせ、有効条件数を大幅に改善します。
• ランク削減技術（Rank-reduction）:
  • 対角スケーリングだけで既に改善される基底関数を除外し、本当に問題となる（近似的な線形依存性を生む）基底関数のみを $Z$ に含めることで、デフレーション行列のランクを削減し、計算コストを低減するヒューリスティックなアルゴリズムを提案しました。
• アルゴリズム: 提案された**デフレーション事前条件付き共役勾配法（DPCG）**を実装し、停止基準として残差と有効条件数を監視する方式を採用しました。

3. 主要な貢献

1. 対角スケーリングの限界の明確化: 既存の文献で「対角スケーリングは滑らかなスプライン基底では十分である」という直観的議論に対し、基底の選択、カット構成、ノット位置が条件数に決定的な影響を与えることを理論的に証明し、2 次元における複雑な依存関係を解明しました。
2. 既存手法のロバスト性欠如の証明: SIPIC や Schwarz 事前条件付けが、一見単純な幾何形状であっても失敗する具体的なカウンター例を提示し、これらの手法が「完全にロバスト」ではないことを示しました。
3. 新しいロバストな事前条件付け手法の提案: 浸り込み FEM の特性（小さな要素に起因する少数の小さな固有値）に特化したデフレーション手法を初めて適用し、対角スケーリングと組み合わせることで、トリミング構成に依存しないロバストな事前条件付けを実現しました。
4. ランク削減アルゴリズムの開発: 不要な基底関数を除外する効率的な手法を開発し、デフレーション行列のサイズを縮小しつつ、条件数の改善効果を維持することを示しました。

4. 数値結果

多様なテストケース（$L^2$ 射影、ポアソン問題、波動方程式の陰的時間積分など）において、提案手法を既存手法（Jacobi, SIPIC, Schwarz）と比較しました。

• 条件数の改善: 提案されたデフレーション手法（ランク削減あり・なし）は、トリミング要素のサイズ（$\delta$）が小さくなっても、条件数が一定に保たれることを示しました。一方、Jacobi や SIPIC は $\delta$ の減少に伴い条件数が急増しました。
• 反復回数と収束性:
  • 提案手法は、トリミング構成や多項式次数に依存せず、安定して少ない反復回数で収束しました。
  • 特に、SIPIC や Schwarz が失敗するケース（例：複数のリッジを持つ構造や、特定の回転角度を持つ格子構造）において、提案手法のみが成功して解を求めました。
  • ランク削減を適用しても、条件数の改善効果はほぼ維持され、計算効率の向上が期待できました。
• 時間依存問題への適用: 波動方程式の陰的時間積分シミュレーションにおいても、時間ステップごとの線形方程式求解において、提案手法が安定した収束を示し、他の手法では 1000 反復以上で発散または停滞するケースでも成功しました。

5. 意義と結論

• 実装の非侵襲性: この手法は、変分形式や近似空間そのものを変更せず、単に反復解法に追加される「事前条件付け」として機能します。したがって、既存の FEM ソフトウェアへの統合が容易で、安定化手法が引き起こす精度低下のリスクがありません。
• 理論的・実用的なブレイクスルー: 浸り込み FEM/IGA における「トリミングによる悪条件化」に対する、理論的に裏付けられ、数値的にロバストな解決策を提供しました。
• 今後の展望: 本手法は対称正定行列に焦点を当てていますが、非対称や不定行列（混合要素法など）への拡張や、$h$-$p$ ロバストなマルチグリッド法との組み合わせが今後の課題として挙げられています。

総じて、この論文は、浸り込み有限要素法の実用化における最大の障壁の一つである「悪条件化」に対し、既存の複雑な安定化手法に代わる、シンプルかつ極めて強力なデフレーションベースの解決策を提示した画期的な研究と言えます。