[論文レビュー] A Study of the Stability Properties of SPH

Joseph P. Morris|arXiv (Cornell University)|Apr 3, 1995

Fluid Dynamics Simulations and Interactions被引用数 22

ひとこと要約

本稿は、磁気力学（MHD）や非圧縮性流れのように応力が負になる状況において、運動量保存型のスムージング粒子法（SPH）に生じる不安定性を調査する。高次のスプラインカーネル、特にガウス関数に近いものでは、不正な横方向モードの抑制と不安定性の成長率の低減が顕著に見られ、これにより安定性が著しく向上することが示された。一方、圧力差を用いる代替形式では、これらの問題が完全に回避される。

ABSTRACT

When using a formulation of Smooth Particle Hydrodynamics (SPH) which conserves momentum exactly the motion of the particles is observed to be unstable to negative stress. It is also found that under normal circumstances a lattice of SPH particles is potentially unstable to transverse waves. This document is a summary of a detailed report (Morris 1994) investigating the nature of these and other instabilities in depth. Approaches which may be used to eliminate these instabilities are suggested. It is found that the stability properties of SPH in general improve as higher order spline interpolants, approximating a Gaussian, are used as kernels.

研究の動機と目的

• 運動量保存型SPH形式において応力が負になる場合に生じる数値的不安定性の根本的要因を特定すること。
• 特に1次元および格子構造に基づく設定におけるSPHの安定性を分析すること。
• 高次のスプラインカーネル（特にガウス関数に近いもの）が不安定モードの成長率に与える影響を評価すること。
• 負の応力や調整可能な背景圧力が関与する問題においても運動量保存を維持しながら不安定性を回避できる代替形式を提案すること。
• MHDや非圧縮性流れを含む複雑な流体動力学問題に適した安定的かつ高精度なSPHカーネルおよび形式の選定ガイドラインを提供すること。

提案手法

• 対称的カーネル（ガウス関数および立方スプラインを含む）を用いた1次元SPHの線形化分散関係を導出し、波動伝播および安定性を分析する。
• 離散的粒子和をフーリエ級数に変換するためのポisson和公式を用いて、周期的格子上での安定性解析を実施する。
• 標準的な運動量保存型形式（背景圧力に依存）と、圧力差を用いる代替形式（背景圧力に依存しない）の2通りの形式を比較する。
• 数値解析を用いて、長方形およびヘキサゴナル格子における不安定な横方向モードの成長率を計算する。
• カーネルのコンパクトサポートおよびフーリエ変換の減衰特性が不安定性の成長に与える影響を評価し、ガウス関数に近いカーネルでは不正なモード成長が抑制されることを示す。
• 応力テンソルを正の成分（運動量保存型SPHで処理）と残差（差分法により処理）に分割することで、安定性を維持する手法を提案する。

実験結果

リサーチクエスチョン

• RQ1 運動量保存型SPHが負の応力下で不安定化する理由は何か？ その背後にある根本的メカニズムは何か？
• RQ2 特に高次のスプラインカーネル（ガウス関数に近いもの）が、1次元および格子構造に基づく流れにおけるSPHの安定性にどのように影響を与えるか？
• RQ3 運動量方程式における作用・反作用ペアの力の回避を図る代替形式は、負の応力に起因する不安定性を完全に排除できるか？
• RQ4 カーネルのフーリエ変換の減衰率は、SPHにおける不安定な横方向モードの抑制に果たす役割は何か？
• RQ5 粒子格子の幾何構造（長方形対ヘキサゴナル）および次元（1次元、2次元、3次元）は、不安定性の発生および成長にどのように影響を与えるか？

主な発見

• 標準的な運動量保存型SPH形式は、粒子同士のペアリングと凝集が生じ、粒子が互いのポテンシャル井戸に落ち込むことによって負の応力下で不安定化する。
• 高次のスプラインカーネル、特にガウス関数に近いものでは、横方向モードの不安定性が顕著に抑制され、実用的なスムージング長さでは成長率が無視できる程度にまで低下する。
• 圧力差形式（式9）は背景圧力依存性を完全に排除し、不安定性を完全に回避するが、運動量の正確な保存は失われる。
• フーリエ変換の減衰が速いカーネル（ガウス関数や高次のスプラインカーネル）を用いることで安定性が著しく向上し、不正なモード成長が抑制される。
• コンパクトサポートを持つカーネルを用いた場合、長方形格子およびヘキサゴナル格子の両方で横方向モードの不安定性が生じるが、高次のカーネルでは成長率が著しく低下する。
• 複雑な状態方程式や強い衝撃波が関与する場合、応力テンソルを正の成分と残差に分割し、差分法を用いることで安定性を維持しつつ精度を保つ戦略が有効である。