Enhancement of stability and accuracy of the moving particle semi-implicit method

As a Lagrangian gridless particle method, the MPS (Moving Particle Semi-implicit) method has been proven useful in a wide range of engineering applications. Up to now, most of MPS applications have been limited to problems with compressive stress states. This paper investigates the performance and stability of MPS method in simulation of more general hydrodynamic problems, including those characterized by tensile stresses or by changes in the stress states. It is shown that MPS-based simulations are prone to become destabilized in presence of attractive interparticle forces. This instability appears to be similar to the so-called tensile instability in SPH method. Two new modifications, namely, a modified Poisson Pressure Equation and a corrective matrix for pressure gradient model, are proposed to resolve this shortcoming. These two new modifications together with two previously proposed improvements are shown to stabilize and enhance the performance of MPS method.