[論文レビュー] Two-level nonlinear Schwarz methods - a parallel implementation with application to nonlinear elasticity and incompressible flow problems

Kyrill Ho, Axel Klawonn+1 more|arXiv (Cornell University)|Mar 25, 2026

Advanced Numerical Methods in Computational Mathematics被引用数 0

ひとこと要約

論文は、GDSW型粗空間（TrilinosのFROSch）を用いた初の並列二レベル非線形Schwarz実装を提示し、非線形Navier–StokesおよびNeo-Hookeanねじれ弾性問題に対して強いスケーラビリティとロバスト性を示し、多くの設定で標準のNewton-Krylov-Schwarz法より優れていることを実証する。

ABSTRACT

Nonlinear Schwarz methods are a type of nonlinear domain decomposition method used as an alternative to Newton's method for solving discretized nonlinear partial differential equations. In this article, the first parallel implementation of a two-level nonlinear Schwarz method leveraging the GDSW-type coarse spaces from the Fast and Robust Overlapping Schwarz (FROSch) framework in Trilinos is presented. This framework supports both additive and hybrid two-level nonlinear Schwarz methods and makes use of modifications to the coarse spaces constructed by FROSch to further enhance the robustness and convergence speed of the methods. Efficiency and excellent parallel performance of the software framework are demonstrated by applying it to two challenging nonlinear problems: the two-dimensional lid-driven cavity problem at high Reynolds numbers, and a Neo-Hookean beam deformation problem. The results show that two-level nonlinear Schwarz methods scale exceptionally well up to 9\,000 subdomains and are more robust than standard Newton-Krylov-Schwarz solvers for the considered Navier-Stokes problems with high Reynolds numbers or, respectively, for the nonlinear elasticity problems and large deformations. The new parallel implementation provides a foundation for future research in scalable nonlinear domain decomposition methods and demonstrates the practical viability of nonlinear Schwarz techniques for large-scale simulations.

研究の動機と目的

• 離散化された非線形偏微分方程式に対するニュートン法の代替として、スケーラブルな非線形領域分解ソルバーを動機づけ開発する。
• ロバストな粗空間（GDSW/RGDSW）を用いて一段階の非線形Schwarz法を二段階スキームへ拡張し、収束性とスケーラビリティを改善する。
• FROSch（Trilinos）内に並列二段階非線形Schwarzフレームワークを実装し、FEDDLibへの統合を通じて性能評価を行う。
• 高レイノルズ数のリード付きキャビティ流れとNeo-Hookeanビーム変形といった挑戦的な問題でのロバスト性と効率を示す。
• 粗空間設計とそれが非線形収束性とスケーラビリティに与える影響に関する洞察を提供する。

提案手法

• ASPIN/RASPEN一段階非線形Schwarzフレームワークとその接線計算を説明する。
• 局所非線形補正と粗非線形補正から構成される二段階非線形Schwarzのバリアント（加法的・ハイブリッド）を導入する。
• 非線形ヌル空間様性を満たすようなGDSW型粗空間（RGDSWおよびMsFEM系）を構築・適用し、スケーラビリティを向上させる。
• モジュラーなクラス構造と幽霊層を介した並列組み立てに焦点を置き、FROSchフレームワーク（Trilinos）内にソルバーを実装する。
• 内側Newton反復で局所非線形問題を解き、外側Newton反復の接線を組み立てる。
• Navier–Stokesおよび弾性問題で速度と圧力成分を結合するモノリシックな粗空間を使用し、Newton法で粗非線形問題を解く。

実験結果

リサーチクエスチョン

• RQ1 二段階非線形Schwarz法は、弾性と非圈縮性流れの非線形PDEで一段階法より収束性とロバスト性を改善するか。
• RQ2 粗空間の選択（RGDSW対GDSW対MsFEM）は、並列設定でのスケーラビリティ、収束性、解出時間にどのような影響を与えるか。
• RQ3 二段階非線形Schwarzの最大サブドメイン数千規模までの並列スケーラビリティはどの程度達成可能で、Newton-Krylov-Schwarzと比較してどうか。
• RQ4 非線形の二段階アプローチは、標準ソルバーより高いレイノルズ数と大きな変形をよりロバストに扱えるか。

主な発見

• 二段階非線形Schwarz法は最大9,000サブドメインまで非常に良くスケールする。
• 非線形Schwarzの実装は、線形化系を線形Schwarz前処理によって解く古典的なNewtonアプローチより一般に優れている。
• 二段階法は多くのケースで弱いスケーラビリティが強化され、解出時間が速く、レイノルズ数の上昇や大変形に対してもロバスト性が高い。
• RGDSW/MsFEM/GDSW粗空間は、非線形設定での収束を改善する堅牢な粗補正を提供する。
• 複数の非線形一段階Schwarz法と比較して、二段階アプローチは対象としたNavier–Stokesおよび弾性問題で性能が向上する。