A Stabilized Incompressible SPH Method by Relaxing the Density Invariance Condition

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.

[1]  J. Monaghan,et al.  Smoothed particle hydrodynamics: Theory and application to non-spherical stars , 1977 .

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

[3]  Jacek Pozorski,et al.  SPH computation of incompressible viscous flows , 2002 .

[4]  S. Koshizuka,et al.  Moving-Particle Semi-Implicit Method for Fragmentation of Incompressible Fluid , 1996 .

[5]  C. K. Thornhill,et al.  Part IV. An experimental study of the collapse of liquid columns on a rigid horizontal plane , 1952, Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences.

[6]  Min-Cheol Ryu,et al.  Numerical simulation of impact loads using a particle method , 2010 .

[7]  L. Lucy A numerical approach to the testing of the fission hypothesis. , 1977 .

[8]  Hitoshi Gotoh,et al.  Enhancement of stability and accuracy of the moving particle semi-implicit method , 2011, J. Comput. Phys..

[9]  Hitoshi Gotoh,et al.  ENHANCED PREDICTIONS OF WAVE IMPACT PRESSURE BY IMPROVED INCOMPRESSIBLE SPH METHODS , 2009 .

[10]  Timothy Nigel Phillips,et al.  A consistent reflected image particle approach to the treatment of boundary conditions in smoothed particle hydrodynamics , 2009 .

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

[12]  S. Shao,et al.  Corrected Incompressible SPH method for accurate water-surface tracking in breaking waves , 2008 .

[13]  Arthur Veldman,et al.  A Volume-of-Fluid based simulation method for wave impact problems , 2005 .

[14]  Ronald Fedkiw,et al.  Two-Way Coupled SPH and Particle Level Set Fluid Simulation , 2008, IEEE Transactions on Visualization and Computer Graphics.

[15]  Masayuki Tanaka,et al.  Stabilization and smoothing of pressure in MPS method by Quasi-Compressibility , 2010, J. Comput. Phys..

[16]  J. Meijerink,et al.  An iterative solution method for linear systems of which the coefficient matrix is a symmetric -matrix , 1977 .

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

[18]  Guirong Liu,et al.  Smoothed Particle Hydrodynamics: A Meshfree Particle Method , 2003 .

[19]  I. J. Schoenberg Contributions to the problem of approximation of equidistant data by analytic functions. Part A. On the problem of smoothing or graduation. A first class of analytic approximation formulae , 1946 .

[20]  Hitoshi Gotoh,et al.  SPH-LES Model for Numerical Investigation of Wave Interaction with Partially Immersed Breakwater , 2004 .

[21]  HighWire Press Philosophical Transactions of the Royal Society of London , 1781, The London Medical Journal.

[22]  Damien Violeau,et al.  Numerical modelling of complex turbulent free‐surface flows with the SPH method: an overview , 2007 .

[23]  Nikolaus A. Adams,et al.  A constant-density approach for incompressible multi-phase SPH , 2009, J. Comput. Phys..

[24]  I. J. Schoenberg Contributions to the Problem of Approximation of Equidistant Data by Analytic Functions , 1988 .

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

[26]  Nikolaus A. Adams,et al.  Angular-momentum conservative smoothed particle dynamics for incompressible viscous flows , 2006 .

[27]  Behzad Ataie-Ashtiani,et al.  Numerical simulation of landslide impulsive waves by incompressible smoothed particle hydrodynamics , 2008 .

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

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

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

[31]  Rui Xu,et al.  Accuracy and stability in incompressible SPH (ISPH) based on the projection method and a new approach , 2009, J. Comput. Phys..

[32]  J. C. Martin An experimental study of the collapse of liquid column on a rigid horizontal plane , 1952 .

[33]  Hitoshi Gotoh,et al.  Turbulence particle models for tracking free surfaces , 2005 .

[34]  A. Chorin Numerical solution of the Navier-Stokes equations , 1968 .

[35]  Stéphane Ploix,et al.  Application of weakly compressible and truly incompressible SPH to 3-D water collapse in waterworks , 2010 .

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

[37]  Seiichi Koshizuka,et al.  Improvement of stability in moving particle semi‐implicit method , 2011 .

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