A remedy for numerical oscillations in weakly compressible smoothed particle hydrodynamics

Weakly Compressible Smoothed Particle Hydrodynamics (WCSPH) can lead to non-physical oscillations in the pressure and density fields when simulating incompressible flow problems. This in turn may result in tensile instability and sometimes divergence. In this paper, it is shown that this difficulty originates from the specific form of spatial discretization used for the pressure term when solving the mass conservation equation. After describing the pressure–velocity decoupling problem associated with the so-called colocated grid methods, a modified approach is presented that overcomes this problem using a different discretization scheme for the second derivative of pressure. The modified scheme is employed for solving a number of benchmark problems including both single-phase and two-phase test cases. Copyright © 2010 John Wiley & Sons, Ltd.

[1]  Nikolaus A. Adams,et al.  A multi-phase SPH method for macroscopic and mesoscopic flows , 2006, J. Comput. Phys..

[2]  Holger Wendland,et al.  Piecewise polynomial, positive definite and compactly supported radial functions of minimal degree , 1995, Adv. Comput. Math..

[3]  C. Rhie,et al.  Numerical Study of the Turbulent Flow Past an Airfoil with Trailing Edge Separation , 1983 .

[4]  Rui Xu,et al.  Comparisons of weakly compressible and truly incompressible algorithms for the SPH mesh free particle method , 2008, J. Comput. Phys..

[5]  J. Morris,et al.  Modeling Low Reynolds Number Incompressible Flows Using SPH , 1997 .

[6]  J. Monaghan,et al.  SPH elastic dynamics , 2001 .

[7]  Paul W. Cleary,et al.  Modelling confined multi-material heat and mass flows using SPH , 1998 .

[8]  J. Monaghan.  SPH without a Tensile Instability , 2000 .

[9]  S. Shao,et al.  INCOMPRESSIBLE SPH METHOD FOR SIMULATING NEWTONIAN AND NON-NEWTONIAN FLOWS WITH A FREE SURFACE , 2003 .

[10]  J. Monaghan.  On the problem of penetration in particle methods , 1989 .

[11]  J. Monaghan.  Simulating Free Surface Flows with SPH , 1994 .

[12]  Siamak Kazemzadeh Hannani,et al.  A fully explicit three‐step SPH algorithm for simulation of non‐Newtonian fluid flow , 2007 .

[13]  S. Cummins,et al.  An SPH Projection Method , 1999 .

[14]  J. Monaghan.  Smoothed particle hydrodynamics , 2005 .

[15]  S. Miyama,et al.  Numerical Simulation of Viscous Flow by Smoothed Particle Hydrodynamics , 1994 .

[16]  F. Harlow,et al.  Numerical Calculation of Time‐Dependent Viscous Incompressible Flow of Fluid with Free Surface , 1965 .

[17]  Pep Español,et al.  Incompressible smoothed particle hydrodynamics , 2007, J. Comput. Phys..

[18]  A. Colagrossi,et al.  Numerical simulation of interfacial flows by smoothed particle hydrodynamics , 2003 .

[19]  Joe J. Monaghan,et al.  SPH simulations of swimming linked bodies , 2008, J. Comput. Phys..

[20]  Nikolaus A. Adams,et al.  An incompressible multi-phase SPH method , 2007, J. Comput. Phys..