[論文レビュー] A note on robust preconditioners for monolithic fluid-structure interaction systems of finite element equations

Ulrich Langer, Huidong Yang|arXiv (Cornell University)|Dec 21, 2014

Advanced Numerical Methods in Computational Mathematics参考文献 55被引用数 12

ひとこと要約

本稿では、結合有限要素行列の近似 $LDU$ 因数分解を用いて、モノリシック流体-構造相互作用（FSI）システムのための頑健な前処理行列を提案する。$L$、$D$、$U$ を近似因子で構築し、スチェル補完を近似するための特化した代数多重グリッド（AMG）サイクルを用いることで、メッシュ、流体密度、時間刻み幅の変動に対しても頑健性を達成。従来の前処理行列と比較して、GMRES反復回数を約10倍削減する。

ABSTRACT

In this note, we consider preconditioned Krylov subspace methods for discrete fluid-structure interaction problems with a nonlinear hyperelastic material model and covering a large range of flows, e.g, water, blood, and air with highly varying density. Based on the complete $LDU$ factorization of the coupled system matrix, the preconditioner is constructed in form of $\hat{L}\hat{D}\hat{U}$, where $\hat{L}$, $\hat{D}$ and $\hat{U}$ are proper approximations to $L$, $D$ and $U$, respectively. The inverse of the corresponding Schur complement is approximated by applying one cycle of a special class of algebraic multigrid methods to the perturbed fluid sub-problem, that is obtained by modifying corresponding entries in the original fluid matrix with an explicitly constructed approximation of the exact perturbation coming from the sparse matrix-matrix multiplications.

研究の動機と目的

• 既存のモノリシック FSI システムにおける前処理付き Krylov 法の流体密度の変動に対する頑健性の欠如を是正すること。
• ブロック Gauss-Seidel や ILU タイプの前処理行列が示すメッシュ依存性および時間刻み幅に敏感な特性を克服すること。
• 水や空気を含む広範な流体特性の範囲で、効率性と収束の頑健性を維持する前処理行列を開発すること。
• 完全に結合されたモノリシック多重グリッドソルバに依存せずに、実装および計算のオーバーヘッドを低減しながら頑健性を達成すること。

提案手法

• 前処理行列は、結合 FSI システム行列の近似 $LDU$ 因数分解により構築され、$L$、$D$、$U$ の正確な因子の近似として $\t{L}$、$\t{D}$、$\t{U}$ を用いる。
• スチェル補完の逆行列は、摂動を加えた流体部分行列に特化した代数多重グリッド（AMG）法を1サイクル適用することで近似する。
• 摂動は、スパース行列同士の積を用いて明示的に導出され、構造が流体系に与える結合効果を捉える。
• 完全な $LDU$ 因数分解を避けるために、低ランク近似と構造化されたスムージング手法を用いることで、計算効率を維持する。
• 前処理行列は GMRES およびフレキシブルな GMRES ソルバに統合され、流体密度およびメッシュの細分化の変動に応じた性能が評価された。

実験結果

リサーチクエスチョン

• RQ1 近似 $LDU$ 因数分解に基づく前処理行列は、流体密度の変動に対してもモノリシック FSI システムで頑健な収束を達成できるか？
• RQ2 提案された前処理行列は、メッシュの細分化および時間刻み幅の変動に対しても頑健性を維持するか？
• RQ3 反復回数とスケーラビリティの観点から、ブロック Gauss-Seidel、SSOR、および ILU(0) タイプと比較して、新しい前処理行列の性能はどの程度か？
• RQ4 摂動を加えた流体行列に1回の AMG サイクルを適用することで、スチェル補完の逆行列を効率的に近似できるか？

主な発見

• 新しい前処理行列は、$\t{P}_{ILU}$ や $\t{P}_{SSOR}$ のような古い前処理行列と比較して、GMRESの反復回数を約10倍削減する。特に、空気（$\rho_f = 0.0011$ g/cm³）のような低密度流体において顕著な改善が得られる。
• $\rho_f = 1.1$ g/cm³（水に近い）の場合、追加質量効果という FSI ソルバで知られる課題に対しても、反復回数が低く安定していることから、本手法の頑健性が確認された。
• 時間刻み幅 $\Delta t = 0.00125$ ms の場合、3つのメッシュレベルにわたり反復回数がほとんど変化せず、メッシュ細分化に対する頑健性が示された。
• 流体密度 $\rho_f \in \{1.1, 0.11, 0.011, 0.0011\}$ g/cm³ の広い範囲で安定した性能を維持し、密度が低下するに従い反復回数が安定または減少する傾向を示した。
• フレキシブルな GMRES は標準 GMRES よりわずかに性能が劣るが、両者とも従来の手法と比較して著しく効率的であった。