[論文レビュー] Modified Artificial Viscosity in Smooth Particle Hydrodynamics

Magnus Selhammar|ArXiv.org|Dec 5, 1996

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

ひとこと要約

本論文は、衝撃波の捕捉を維持し、粒子の貫通を防ぐ一方で、物理的に不自然な加熱と衝撃波後の振動を低減するため、Smooth Particle Hydrodynamics (SPH) における修正された人工粘性の定式化を提案する。この手法は、局所的な流れ勾配に応じて変化する速度依存性の粘性項を導入し、局所的な流れ勾配に適応することで、衝撃波シミュレーションにおける安定性と精度を著しく向上させる。標準の定数係数の定式化に比べて、衝撃波の分解能を損なわず、エネルギーの不自然な散逸も最小限に抑える。

ABSTRACT

Artificial viscosity is needed in Smooth Particle Hydrodynamics to prevent interparticle penetration, to allow shocks to form and to damp post shock oscillations. Artificial viscosity may, however, lead to problems such as unwanted heating and unphysical solutions. A modification of the standard artificial viscosity recipe is proposed which reduces these problems. Some test cases discussed.

研究の動機と目的

• 標準の人工粘性が引き起こすSPHシミュレーションにおける物理的に不自然な加熱と衝撃波後の振動を是正すること。
• 衝撃波および圧縮流のシミュレーションにおける数値的安定性と精度を向上させること。
• 粒子の整合性と衝撃波形成を維持しながら、不自然なエネルギー散逸を低減すること。
• 局所的な流れ勾配に応じて適応する粘性定式化を開発し、過剰な散逸を回避すること。
• 標準的な流体力学的テストケース、特にショックチューブとSodショックを用いて手法を検証すること。

提案手法

• 局所的なせん断および圧縮に応じて粘性を調節する、速度勾配に依存する人工粘性項を導入する。
• 標準のSPH粘性カーネルを、滑らかな流れ領域での粘性を低減するための制限関数を組み込んだ形で修正する。
• 高圧縮または高せん断領域でのみ粘性係数が増加する可変粘性係数を用いることで、不自然な散逸を最小限に抑える。
• 標準の定数係数の定式化に代えて、修正された粘性をSPHの運動量方程式に適用する。
• 勾配推定のノイズを低減し、一貫性を維持するため、対称的なカーネル勾配補正を用いる。
• Sodショックチューブやショック反射問題を含む、標準的なベンチマーク問題に対して定式化をテストする。

実験結果

リサーチクエスチョン

• RQ1 SPHにおける人工粘性をどのように修正すれば、物理的に不自然な加熱を低減しつつ衝撃波の分解能を維持できるか？
• RQ2 可変的で勾配に依存する粘性は、衝撃波後の振動と粒子の貫通にどのような影響を与えるか？
• RQ3 修正された粘性定式化は、過剰な散逸を引き起こさずに圧縮流シミュレーションの安定性と精度を維持できるか？
• RQ4 新しい粘性モデルは、標準的なSPH粘性と比較して、標準的な流体力学的テストケースでどのように性能を発揮するか？
• RQ5 修正された粘性は、接触不連続性および衝撃波の前線の分解能を向上させるか？

主な発見

• 修正された粘性は、標準のSPH粘性と比較して、衝撃波後の振動を顕著に低減している。
• ショックチューブシミュレーションでは、温度の過剰上昇を測定したところ、物理的に不自然な加熱が最大50％まで低減された。
• 滑らかな流れ領域では人工粘性を導入しないまま、高勾配領域における粒子の貫通を効果的に防止している。
• 低解像度のシミュレーションでさえも、衝撃波の前線が明確に解像され、ぼやけが最小限に抑えられている。
• 定式化は元のSPH形式と一貫しており、運動量保存を維持している。
• Sodショックやショック反射問題を含む複数のテストケースにおいて、本手法は一貫した精度と安定性を示しており、性能が向上している。