Journal of Applied Mathematics
Volume 2012 (2012), Article ID 139583, 24 pages
Research Article

A Stabilized Incompressible SPH Method by Relaxing the Density Invariance Condition

1Department of Civil Engineering, Kyushu University, 744 Motooka, Nishi–ku, Fukuoka 819-0395, Japan
2Faculty of Education and Human Science, Yokohama National University, 79-1 Tokiwadai, Hodogaya-ku, Yokohama 240-8501, Japan

Received 5 January 2012; Accepted 16 March 2012

Academic Editor: Hiroshi Kanayama

Copyright © 2012 Mitsuteru Asai et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


A stabilized Incompressible Smoothed Particle Hydrodynamics (ISPH) is proposed to simulate free surface flow problems. In the ISPH, pressure is evaluated by solving pressure Poisson equation using a semi-implicit algorithm based on the projection method. Even if the pressure is evaluated implicitly, the unrealistic pressure fluctuations cannot be eliminated. In order to overcome this problem, there are several improvements. One is small compressibility approach, and the other is introduction of two kinds of pressure Poisson equation related to velocity divergence-free and density invariance conditions, respectively. In this paper, a stabilized formulation, which was originally proposed in the framework of Moving Particle Semi-implicit (MPS) method, is applied to ISPH in order to relax the density invariance condition. This formulation leads to a new pressure Poisson equation with a relaxation coefficient, which can be estimated by a preanalysis calculation. The efficiency of the proposed formulation is tested by a couple of numerical examples of dam-breaking problem, and its effects are discussed by using several resolution models with different particle initial distances. Also, the effect of eddy viscosity is briefly discussed in this paper.