[論文レビュー] Scalable smoothing strategies for a geometric multigrid method for the immersed boundary equations

Amneet Pal Singh Bhalla, Matthew G. Knepley+3 more|arXiv (Cornell University)|Dec 7, 2016

Lattice Boltzmann Simulation Studies参考文献 34被引用数 5

ひとこと要約

本稿では、ストークスおよびナビエ＝ストークス流れの条件下で浸漬境界（IB）方程式に適用される幾何的多重グリッド・ソルバーに、スケーラブルで並列化可能な平滑化戦略を開発する。加法的ヴァンカ平滑化法とスシュール補間に基づく平滑化法を導入し、低〜中程度のレイノルズ数および高材料剛性において有効な性能を示し、複雑な流体構造挙動問題においてもスケーラビリティを実現する。

ABSTRACT

The immersed boundary (IB) method is a widely used approach to simulating fluid-structure interaction (FSI). Although explicit versions of the IB method can suffer from severe time step size restrictions, these methods remain popular because of their simplicity and generality. In prior work (Guy et al., Adv Comput Math, 2015), some of us developed a geometric multigrid preconditioner for a stable semi-implicit IB method under Stokes flow conditions; however, this solver methodology used a Vanka-type smoother that presented limited opportunities for parallelization. This work extends this Stokes-IB solver methodology by developing smoothing techniques that are suitable for parallel implementation. Specifically, we demonstrate that an additive version of the Vanka smoother can yield an effective multigrid preconditioner for the Stokes-IB equations, and we introduce an efficient Schur complement-based smoother that is also shown to be effective for the Stokes-IB equations. We investigate the performance of these solvers for a broad range of material stiffnesses, both for Stokes flows and flows at nonzero Reynolds numbers, and for thick and thin structural models. We show here that linear solver performance degrades with increasing Reynolds number and material stiffness, especially for thin interface cases. Nonetheless, the proposed approaches promise to yield effective solution algorithms, especially at lower Reynolds numbers and at modest-to-high elastic stiffnesses.

研究の動機と目的

• 浸漬境界（IB）法における幾何的多重グリッド・ソルバーにおいて、既存のヴァンカ型平滑化法の並列スケーラビリティに限界がある問題に対処する。
• ストークスおよびナビエ＝ストークス条件の下で、半陽的IB定式化に適合する並列高性能コンピューティングに適した新たな平滑化戦略を開発する。
• 変動する材料剛性およびレイノルズ数を有する流体構造挙動問題におけるソルバー効率とスケーラビリティを向上させる。
• 半陽的時間離散化によるIB方程式から生じる線形系のロバストでスケーラブルな解法を可能にする。
• 薄いおよび厚い構造的モデルを含む多様な物理的状態において、新規平滑化法の性能を調査する。

提案手法

• ストークス-IB系における収束のロバスト性を維持しつつ、並列性を向上させるために、ヴァンカ平滑化法の加法的バージョンを採用する。
• 線形系のブロック構造を活用して収束速度を向上させる、スシュール補間に基づく平滑化法を導入する。
• 局所的に細分化されたカルテジアングリッド上で、標準的なガラーキンのコarseningおよび prolongation 演算子を用いた幾何的多重グリッドを実装する。
• EulerianおよびLagrangian方程式の空間離散化に2次精度の有限差分法を用いる。
• 正確で安定したLagrangian-Eulerian結合を実現するため、ペスキンの4点正則化デルタ関数を適用する。
• Krylov部分空間ソルバー用の幾何的多重グリッド予約子に平滑化法を統合し、鞍部問題のスケーラブルな解法を可能にする。

実験結果

リサーチクエスチョン

• RQ1 加法的ヴァンカ平滑化法は、浸漬境界方程式の幾何的多重グリッドフレームワークにおいて、効果的に並列化を可能にしながらも、有効に使用可能か？
• RQ2 スシュール補間に基づく平滑化法は、従来の平滑化法と比較して、ストークス-IB系において優れた収束特性を示すか？
• RQ3 浸漬境界定式化において、レイノルズ数および材料剛性の増加に伴い、ソルバー性能はどのように低下するか？
• RQ4 異なる構造的厚さおよび物理的パラメータにおいて、提案された平滑化法のスケーラビリティおよびロバスト性はいかがなものか？
• RQ5 新規平滑化法は、流体構造挙動シミュレーションにおいて、低レイノルズ数から中〜高レイノルズ数の両方の流れにおいても、効率を維持できるか？

主な発見

• 加法的ヴァンカ平滑化法は、従来の点別ヴァンカ平滑化法の代替として効果的かつ並列化可能であり、分散メモリアーキテクチャにおけるスケーラビリティを著しく向上させる。
• スシュール補間に基づく平滑化法は、特に高材料剛性および複雑な界面幾何形状を有する問題において、強力な収束性能を示す。
• 線形ソルバーの性能はレイノルズ数および材料剛性の増加に伴い低下するが、特に薄い界面ケースでは顕著である。一方、新規平滑化法は低〜中程度のレイノルズ数においては依然として有効である。
• 提案されたソルバーは、ストークスおよび非ストークス流れを含む広範な物理的パラメータ範囲でスケーラブルな性能を達成する。
• 平滑化法は、厚いおよび薄い構造的モデルの両方においてロバストで効率的であり、変動する弾性剛性下でも収束率が安定を保つ。
• パフォーマンスプロファイルの結果、RASおよびRMS平滑化法は適切に最適化されており、ハイパフォーマンスコンputングプラットフォームにおける大規模シミュレーションに適している。