It can be shown that the SPH method is unstable in tension. The instability is shown to
result from an effective stress with a negative modulus (imaginary sound speed)
being produced by the interaction between the constitutive relation and the kernel
function, and is not caused by the numerical time integration algorithm. 1 D.S Balsara 2 states, use special filtering of
velocities (so called conservative smoothing, because momentum
quantities are not modified):図 1.
1Swegle J.W., Hicks D.L., and Attaway S.W., Smoothed Particle
Hydrodynamics: Stability Analysis, Journal of Computational Physics,
Vol. 116, pp. 123-134, 1995.
2Balsara D.S., Von Neumann Stability Analysis of Smoothed Particle
Hydrodynamics Suggestions for Optimal Algorithms, Journal of
Computational Physics, Vol. 121, pp. 357-372, 1995.