[論文レビュー] A modified SPH approach for fluids with large density differences

F. Ott, Erik Schnetter|arXiv (Cornell University)|Mar 27, 2003

Fluid Dynamics Simulations and Interactions参考文献 9被引用数 21

ひとこと要約

本稿では、質量密度の代わりに粒子密度を離散化することで、密度差が大きな多相流のシミュレーションを安定化する修正されたスムージング粒子法（SPH）を提案する。粒子密度を用いた密度再構築に再定義することで、不自然な圧力振動が解消され、衝撃管、音波伝播、および移流問題において精度が向上する。特に、粒子の質量や空間的解像度が著しく異なる場合に顕著である。

ABSTRACT

We introduce a modified SPH approach that is based on discretising the particle density instead of the mass density. This approach makes it possible to use SPH particles with very different masses to simulate multi-phase flows with large differences in mass density between the phases. We test our formulation with a simple advection problem, with sound waves encountering a density discontinuity, and with shock tubes containing a contact discontinuity between air and Diesel oil. For all examined problems where particles have different masses, the new formulation yields better results than standard SPH. This is also the case for problems in which different spatial resolutions are needed while the mass density does not change.

研究の動機と目的

• 密度差が著しい多相流をシミュレートする際の標準SPHの不安定性および不正確さを是正すること。特に、流体界面における問題を対象とする。
• 標準SPHにおける根本的問題、すなわち、不連続部近傍で粒子間の質量差が大きくなると密度および圧力計算が不正確になるという点を克服すること。
• 粒子の質量や空間的解像度が異なる場合でも、質量密度が各相で一定のままでも、安定したシミュレーションを可能にすること。
• 圧縮性流体、特に液体-気体または液体-液体界面を含む、剛性の高い状態方程式を伴う問題において、安定性と精度を維持する定式化を開発すること。
• 接触不連続、衝撃波、音波を伴う高密度対比問題に対して、標準SPHの代替として数値的に安定した代替手法を提供すること。

提案手法

• 質量密度に基づく標準SPHの密度再構築を、粒子体積の逆数として定義される粒子密度に基づくものに置き換える。
• SPHカーネルの和を再解釈し、質量密度の代わりに粒子密度を計算する式を用いる。式は $ \tilde{\rho}_{\text{part}}(\mathbf{x}_i) = \sum_j V_j W_{ij} $ であり、$ V_j $ は粒子体積を表す。
• 粒子質量 $ m_i $ と粒子密度 $ \rho_{\text{part},i} $ を用いて、実際の流体の質量密度 $ \rho_i = m_i / V_i $ を計算し、粒子運動と整合性を保つ。
• 運動方程式（運動量およびエネルギーの時間発展）に標準SPH式を適用するが、ラグランジュ的性質を維持する。
• 衝撃波の安定化のため、$ \alpha = 0.5 $ の人工粘性を用い、約8個のスムージング長さにわたって分散させる。初期不連続性の増幅は生じない。
• 初期データにのみ平滑化を適用する必要がある場合に限り適用するが、不連続な初期条件が新しい定式化を不安定化させないことを示す。

実験結果

リサーチクエスチョン

• RQ1 質量密度の代わりに粒子密度に基づく修正SPH定式化が、密度対比が著しい多相流のシミュレーションにおいて安定性と精度を向上させることができるか？
• RQ2 標準SPHと比較して、密度不連続の移流をシミュレートする際の新定式化の性能は？
• RQ3 修正SPH法は、空気とディーゼル油（密度比約1:80）を含む高密度対比衝撃管において、衝撃波および膨張波を正確に捉えることができるか？
• RQ4 粒子の質量や空間的解像度が異なる場合でも、不自然な圧力振動を引き起こさずに処理できるか？
• RQ5 初期データに不連続性や高周波数モードが含まれる場合でも、新定式化は頑健か？

主な発見

• 修正SPH定式化は、密度不連続を伴う1次元移流問題を安定的にシミュレートでき、不自然な圧力振動が発生せず、界面の維持も正確に再現された。
• 音波テストでは、新手法は正しく透過・反射係数を捉えており、標準SPHは誤った波動挙動と圧力スパイクを生じさせた。
• 空気とディーゼル油（密度比約1:80）を含む衝撃管において、接触不連続は安定に保たれ、時間経過に伴い圧力振動が増幅されることもなかった。
• 衝撃波速度（$ v_s = 1056 \pm 20 \, \text{m/s} $）、衝撃後速度（$ v_D = 5.94 \pm 0.02 \, \text{m/s} $）、圧力上昇（$ \Delta p_D = 4.94 \pm 0.02 \, \text{MPa} $）のシミュレーション結果は、解析解とよく一致した。
• 新手法は解析的衝撃関係と非常に良好に一致し、すべての主要な衝撃パラメータについて相対誤差が2%未満であった。これは、非常に高い精度を示している。
• 標準SPHは、剛性の高い流体では物理的に意味のない圧力誤差が生じ、実際の圧力に対して数個のオーダーも超える場合があり、同じ条件下では安定なシミュレーションを生成できなかった。