[論文レビュー] A multi-resolution SPH method for fluid-structure interactions

Chi Zhang, Massoud Rezavand+1 more|arXiv (Cornell University)|Nov 21, 2019

Fluid Dynamics Simulations and Interactions参考文献 38被引用数 5

ひとこと要約

本稿では、流体と構造物の相互作用（FSI）のためのマルチスケール平滑化粒子法（SPH）を提示する。流体領域と固体領域で空間的・時間的解像度を異なるものにすることで、計算効率を向上させる。位置ベースのVerletスキームと、固体粒子の時間平均速度・加速度を用いることで、運動量保存と正確な力の一致を確保し、生物学的力学系を含む複雑なFSI問題において高い精度と頑健性を示す。

ABSTRACT

In this paper, we present a multi-resolution smoothed particle hydrodynamics (SPH) method for modeling fluid-structure interaction (FSI) problems. By introducing different smoothing lengths and time steps, the spatio-temporal discretization is applied with different resolutions for fluid and structure. To ensure momentum conservation at the fluid-structure coupling, a position-based Verlet time integration scheme is introduced. Furthermore, the time-averaged velocity and acceleration of solid particles are introduced to enhance force matching in the fluid and solid equations. A set of numerical examples including several bio-mechanical problems are considered to demonstrate the efficiency, accuracy and robustness of the present method. A open-source code for all the examples is also provided.

研究の動機と目的

• 大規模な構造物変形を伴うFSI問題における単一解像度SPHシミュレーションの高い計算コストを低減すること。
• 構造物の動態が流体の動態よりも細かい時間解像度を必要とする場合に、均一な時間ステップでは非効率になる問題を克服すること。
• 複数の時間ステップと解像度を用いる際の、流体-構造界面における運動量保存と正確な力の伝達を保証すること。
• フラッピングする魚のような体や静脈弁といった複雑な生物学的力学的FSI系の効率的かつ正確なシミュレーションを可能にすること。
• 再現可能性とFSI研究分野への広範な応用を目的として、堅牢でオープンソースの実装を提供すること。

提案手法

• 流体領域と固体領域で異なるスムージング長と粒子間隔を適用することで、空間的解像度の適応を可能にする。
• 複数の時間ステップ（固体では短く、流体では長く）を用いることで、安定性を保ちつつ計算効率を向上させる。
• 流体-構造界面における運動量保存を強制するために、位置ベースのVerlet時間積分スキームを実装する。
• 各流体時間ステップ内で固体粒子の速度および加速度を時間平均化することで、流体および固体の式における力の一致を向上させる。
• 統一されたSPHフレームワーク内において、流体動力学には更新ラグランジュSPH、固体動力学には全ラグランジュSPHを適用する。
• 弾性構造物にはネオ＝フーケの材料モデルを適用し、ベンチマークケースではスリップ禁止および流入/流出境界条件を強制する。

実験結果

リサーチクエスチョン

• RQ1 可変時間ステップを用いたマルチスケールSPH法は、流体-構造相互作用シミュレーションにおいて計算コストを顕著に低減できるか？
• RQ2 流体と固体領域で異なる時間ステップと空間解像度を用いる場合、運動量保存はどのように維持されるか？
• RQ3 固体粒子の速度および加速度を時間平均化することで、力の一致とシミュレーションの正確性はどの程度向上するか？
• RQ4 本手法は、柔軟体の受動的フラッピングといった複雑な大変形生物学的現象を正確に捉えることができるか？
• RQ5 破砕ダムの崩壊と弾性ゲート、および流れ誘発振動といったベンチマークFSI問題において、本手法はどの程度の性能を示すか？

主な発見

• 本手法は、流体の解像度を粗くし、時間ステップを長くすることで、顕著な計算効率の向上を達成しながらも、高い精度を維持している。
• 位置ベースのVerletスキームにより、異なる時間ステップを用いても、流体-構造界面における運動量保存が成功裏に保持された。
• 固体粒子の時間平均速度および加速度により、より安定的かつ正確な力の伝達が実現され、振動の低減と結合の忠実度向上が達成された。
• 本手法は、交互に発生する渦の対と、周波数比が約2.25:1のリサージュ曲線に類似した尾の軌道を含む、柔軟な魚のような体の受動的フラッピングを正確に再現した。
• 尾の垂直振幅は体幅の72.5％に達し、生物学的力学系に特徴的な大変形を扱える能力を示した。
• 弾性ゲート付きの破砕ダムおよび流れ誘発振動のベンチマーク結果は、本手法の頑健性、正確性、および複雑なFSI応用への適性を確認した。