[論文レビュー] Further investigation of the spurious interface fragmentation in multiphase Smoothed Particle Hydrodynamics

Kamil Szewc, Michał T. Lewandowski|arXiv (Cornell University)|Feb 25, 2016

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

ひとこと要約

この論文は、多相流スムージング粒子法（SPH）における人工的界面補正項の広範な使用を疑問視し、表面張力がない場合に生じるいわゆる「誤った界面破壊」が、物理的に正しいケルビン＝ヘルムホルツ不安定性の現れ方であると主張する。安定性解析とSPH/VOFシミュレーションを用いて、こうした補正項が非物理的解を引き起こす可能性があることを示し、補正パラメータεと平滑化長さhの間の経験的に観察された逆比例関係が、流体力学的安定性理論の直接的結果であることを説明する。

ABSTRACT

This article presents results of further investigation of the problem of spurious interface fragmentation in the multiphase SPH. In order to remove arising instabilities, many authors introduced the artificial interface correction procedure. In the present paper we show that the interface instabilities are physical and the introduction of the interface correction procedure may leads to non-physical solutions. We also explain the puzzling relation between the parameters $\varepsilon$ and $h$. The analysis is performed on the basis of the stability analysis and numerical calculations using SPH and Volume Of Fluid (VOF) approach.

研究の動機と目的

• 表面張力を無視した場合の多相SPHにおける界面不安定性の物理的起源を調査すること。
• 一般的に使用される界面補正項の有効性を疑問視すること。
• 先行研究で観察されたε_min ∝ 1/hという経験的関係を明確にすること。
• サスマンのバブルテストケースにおいて、表面張力が小さくても重要な安定化役を果たしていることを示すこと。
• SPHシミュレーションにおける界面補正手順の使用を慎重に行うべきであると提言すること。

提案手法

• 一様な速度および密度対比を持つ二相流界面に対して線形化されたケルビン＝ヘルムホルツ安定性解析を実施する。
• 界面補正（ε項）をモデル化するための人工的法線応力項を含む分散関係を導出する。
• 先行ベンチマークと整合性を保つために、次元なしの数値単位（R = 1 cm、g = 1 cm/s²、ρ = 1 g/cm³）を用いる。
• 界面補正項を含む・含まないSPHシミュレーションを実施し、VOF法の出力と比較する。
• 数値実験を用いて、最小安定ε（ε_min）が平滑化長さhにどのように依存するかを分析する。
• 安定性条件から理論的にε_min ∝ 1/hの関係を導出し、数値的観察と一致させることを確認する。

実験結果

リサーチクエスチョン

• RQ1 表面張力が無視できる場合にSPHシミュレーションで観察される界面破壊は、数値的誤差であるのか、物理的不安定性であるのか。
• RQ2 人工的界面補正項（ε項）がSPH界面を安定化させる真の物理的役割は何か。
• RQ3 最小安定ε値が、経験的に観察されたように平滑化長さhと逆比例する理由は何か。
• RQ4 サスマンテストケースにおいて、その大きさが小さいにもかかわらず、表面張力が顕著な安定化役を果たしているのか。
• RQ5 観察されたε_min ∝ 1/hの関係を、流体力学的安定性理論から解析的に導出できるか。

主な発見

• 表面張力を無視したSPHシミュレーションで観察される界面不安定性は、物理的に正しいものであり、ケルビン＝ヘルムホルツ不安定性の始まりに対応している。
• 人工的界面補正項（ε項）は物理的欠陥を是正するのではなく、実際の不安定性を隠蔽しており、非物理的解を引き起こす可能性がある。
• 経験的に観察されたε_min ∝ 1/hのスケーリング関係は、ケルビン＝ヘルムホルツ解析から導出された安定性条件によって解析的に説明できる。
• サスマンテストケースにおいては、表面張力の大きさが小さくても、界面の破壊を防ぐために十分な安定化効果を発揮している。
• 界面補正項は主に数値的「美化」の役割を果たしており、物理的安定化要因ではなく、使用には注意が必要である。