[論文レビュー] Toward a 2D SPH Multiphysic Code with Solid-Solid & Fluid Interactions for Industrial Related Problems

David Cébron, Jean‐François Sigrist|arXiv (Cornell University)|Sep 30, 2008

Fluid Dynamics Simulations and Interactions参考文献 15被引用数 2

ひとこと要約

本論文では、多相流および固体相互作用を含む産業スケールのシミュレーションを想定した、拡張されたオープンソース2次元滑らか粒子動力学（SPH）コードSPHYSICS2を提示する。固体同士の接触検出と流体-固体カップリングを統合することで、複雑な衝突および滑り運動の正確なモデル化が可能となり、エネルギー保存則に整合した結果が得られたが、わずかな数値的減衰が生じる場合がある。

ABSTRACT

In the present study, applications of the SPH method to industrial related issues are considered by starting from an existing open source 2D SPH code, namely the SPHYSICS code, which offers an effective ground for numerical developments, which are performed in order to bring an answer to industrial problems, such as simulations of solid/fluid coupling in a free surface flow context. The purpose of the present paper is therefore to expose the numerical developments which yield an enhanced version (referred to as "SPHYSIC2") of the initial code. Firstly, the different features added to obtain the operational code needed for engineering applications are described, and so are the problems raised on this way, offering a kind of review of SPH methods for engineers. Secondly, the validation of the proposed code is partially presented with two well known but difficult test cases, namely the classical "dam break" and "wedge entry" problems. Thirdly, principles of a method to solve solid/solid contacts, frequently present in realistic configuration, are exposed and applied to achieve more complex simulations. Finally, perspectives for new features of the SPHYSIC2 code are exposed and discussed.

研究の動機と目的

• 産業用事前設計応用を想定し、オープンソースのSPHYSICSコードを多相流および多固体相互作用をサポートするように拡張すること。
• 計算コストを低く抑えつつ、効率的に動作する堅牢な固体-固体接触検出アルゴリズムを実装すること。
• 流体-構造および固体-固体相互作用を含むベンチマーク問題に対して、拡張されたコードを検証すること。
• タンク内の水面上で衝突し、滑りながら移動するシリンダーのような、複雑で現実的なシミュレーションの実現可能性を示すこと。
• 将来の産業分野におけるSPHの安全性評価および動的衝突解析への応用基盤を提供すること。

提案手法

• 粒子の近接度に基づく接触検出アルゴリズムをSPHYSICSに統合し、粒子と固体境界との間の単純な距離チェックを用いる。
• 非弾性衝突および滑り挙動を模倣するための法線および接線方向の接触力モデルを実装した。
• 安定性および物理的リアリズムを向上させるために、SPH運動方程式に人工的粘性項および拡散項を追加した。
• 正確な粒子補間と収束性を保証するため、コンパクトサポートおよび正規化を備えたカーネル関数を用いた。
• 流体動力学に弱圧縮性SPH式法と三次スプラインカーネルを適用した。
• 壁粒子法を用いて重力および境界条件をモデル化し、固体境界を表現した。

実験結果

リサーチクエスチョン

• RQ1 SPHYSICSは、産業的意義を持つ固体-固体および流体-固体相互作用を効果的にシミュレートするように拡張可能か？
• RQ2 拡張されたSPHYSICS2モデルは、落下するシリンダーが水に衝突する際の運動およびエネルギー散逸をどの程度正確に再現できるか？
• RQ3 接触アルゴリズムによって導入される数値的減衰の程度は何か？ また、エネルギー保存にどのような影響を及えるか？
• RQ4 自由落下、流体への衝突、境界面での滑りといった複雑な過程を含むシーケンスを、一貫性のある力学的挙動でシミュレート可能か？
• RQ5 造船および機械工学分野における事前設計解析において、SPH手法がどの程度実用的であるか？

主な発見

• SPHYSICS2コードは、水で満たされたタンクにシリンダーが落下する現象を的確にシミュレートでき、衝突、減速、タンク壁への滑りといった主要な力学的挙動を捉えている。
• 自由落下中の速度履歴は理論的予測と良好に一致しており、わずかな数値的減衰が観察されるにとどまる。
• 自由落下中はエネルギー保存が維持されており、衝突時の運動エネルギー損失は主に流体抵抗によるものであり、数値的散逸によるものではない。
• 短時間の接触イベントにおいて、接触検出アルゴリズムが無視できるほどのわずかな数値的減衰しか導入しないため、非弾性衝突のモデル化に適していることが確認された。
• タンク壁に沿ったシリンダーの滑りおよび水面への跳ね返りのシミュレーションは、期待される物理的挙動と定性的・定量的に整合している。
• コードは、複雑なマルチフィジックス産業用シミュレーションの実現可能性に前向きな兆しを示しているが、さらなる検証が継続中である。