A Jacobian-free Newton-Krylov method for cell-centred finite volume solid mechanics

技術要約：セル中心型有限体積固体力学のためのJacobian-free Newton-Krylov法

問題提起
従来の有限体積（FV）法は、主に分離型インプリシット手法（ピカード反復）または完全なエクスプリシット手法に依存している。分離型手法は実装の簡便さとメモリ効率において優れているが、幾何学的非線形性、大きな回転、または複雑な結合を伴う問題において、収束が遅くなる傾向がある。一方、有限要素法のコミュニティで一般的なニュートン・ラプソン型の手法は、二次収束性を提供するが、フル・ヤコビ行列の明示的な組み立て、格納、および反転を必要とする。この要件は計算コストとメモリ消費を増大させる。また、既存の分離型FVフレームワークにブロック結合型ニュートン法を適応させることは、真のヤコビの導出と格納の複雑さから容易ではない。したがって、大規模なコードのリファクタリングを行うことなく、既存の分離型FVフレームワークに統合できる、堅牢かつ効率的なブロック結合型ソルバーが必要とされている。

手法
著者らは、セル中心型有限体積固体力学に特化したJacobian-free Newton-Krylov（JFNK）アプローチを提案している。核となる手法は、有限差分を用いてベクトル（$v$）に対するヤコビ行列（$J$）の作用を近似することで、ヤコビ行列の明示的な形成を回避するものである：
$$Jv \approx \frac{R(u + \epsilon v) - R(u)}{\epsilon}$$
ここで、$R$は残差関数、$u$は解ベクトル、$\epsilon$は摂動スカラーである。この近似は、ニュートン反復内の線形システムを解くために、Generalised Minimal Residual（GMRES）クライロ・部分空間ソルバーと組み合わされる。

効率性と堅牢性を確保するために、本手法は前処理戦略を採用している。正確なヤコビ行列（利用不可能なもの）を使用する代わりに、慣性項と拡散に似た項（分離型ソルバで使用されるRhie-Chow安定化に基づく）から導出されたコンパクト・ステンシル近似ヤコビ行列を利用する。この前処理行列（$P$）は、概念的には従来の分離型アプローチで使用される行列と同様であり、本手法がOpenFOAMのsolids4foamツールボックスのようなフレームワーク内の既存の行列組み立ておよび格納ルーチンを活用することを可能にしている。

実装には、非線形ソルバ（SNES）および線形代数のためにPETScライブラリを使用し、これをsolids4foamFVソルバとインターフェースさせている。主なアルゴリズムの特徴は以下の通りである：

• 前処理器（Preconditioners）: LU直接解法（MUMPS）、不完全LU分解（ILU）、および代数マルチグリッド（HYPRE BoomerAMG）の比較。
• グローバル化（Globalisation）: 解を二次収束領域へと導くためのラインサーチおよび過渡継続（予測子・修正子法）の使用。
• 安定化（Stabilisation）: チェッカーボード現象（ゼロエネルギーモード）を抑制するためにRhie-Chow型の項を残差に追加し、収束速度のために調整されたグローバルスケーリング因子（$\alpha$）を用いる。

主な貢献

1. 有限体積固体力学への初の適用: 本研究は、有限体瘤力学に特化したJFNK法を初めて提案したものであり、流体力学への応用と固体力学の間の溝を埋めるものである。
2. シームレスな統合: コンパクト・ステンシル近似ヤコビ行列を前処理器として使用することで、最小限のコード修正で既存の分離型FVフレームワーク内でニュートン型収束を採用することを可能にした。
3. 包括的なベンチマーク: 次元（2D/3D）、材料挙動（線形弾性、超弾性、弾塑性）、および物理現象（静的、動的、流体構造連成）のバリエーションにわたる一連のケースを用いて、手法を厳密にテストした。

結果
本研究では、計算時間、メモリ使用量、および堅牢性の指標を用いて、JFNK法を分離型アプローチと比較検証した。

• 効率性: JFNK法は、すべての線形および非線形弾性ケースにおいて、分離型アプローチを凌駕した。多くの事例で1〜2桁の高速化が観察され、2D線形弾性Cook's membraneケースでは最大341倍の高速化が記録された。
• メモリ: 自由度が$10^4$未満の問題では、メモリ使用量は分離型手法と同等であった。より大規模な問題では、GMRESソルバおよびマルチグリッド前処理器の格納要件により、JFNKのメモリ使用量は増加し、球状空洞ケースにおいて分離型手法の最大3.6倍に達した。
• 堅牢性と限界:
  • JFNK法は、分離型手法が収束に失敗した超弾性ケース（大きな歪み/回転）において、優れた堅牢性を示した。
  • しかし、JFNK法は弾塑性ケース（塑性を含むCook's membrane）において発散した。著者らは、小歪み弾性をベースとしたコンパクト・ステンシル前処理器が、大きな塑性歪みに対してヤコビ行列の不適切な近似となっているためであると分析している。
  • 前処理器の選択:
    • LU直接解法は、小規模から中規模のサイズおよび2Dケースにおいて最も高速であった。
    • **マルチグリッド（AMG）**は、より大規模な3D問題においてより効果的であり、LUよりも優れたスケーラビリティと低いメモリ増加率を示した。
    • **ILU(k)**は、一般的に最も低速で堅牢性に欠けていた。
  • 安定化: Rhie-Chow安定化因子（$\alpha$）を大きくすると、収束速度と堅牢性は一般に向上したが、漸近的精度の低下という代償を伴った。

意義と主張
本論文は、Jacobian-free Newton-Krylov法が、最小限の変更で導入が可能であることから、既存の分離型フレームワークにとって、特に有限体積固体力学における有望な進歩であることを主張している。本研究は、フル・ヤコビ行列の組み立てに伴う過度なメモリコストを回避しつつ、JFNKが弾性問題において大幅な高速化と堅牢性の向上を実現できることを示している。

著者らは現在の手法の限界についても謙虚であり、現在の前処理戦略では弾塑性に対してはまだ堅牢ではないことを明記している。彼らは、塑性に対する高度な前処理器の開発や、非直交補正の導入が必要な今後の課題として特定している。実装は、コミュニティによる探索と拡張を促進するため、solids4foamツールボックスを通じて公開されている。