An SPH Projection Method

A new formulation is introduced for enforcing incompressibility in Smoothed Particle Hydrodynamics (SPH). The method uses a fractional step with the velocity field integrated forward in time without enforcing incompressibility. The resulting intermediate velocity field is then projected onto a divergence-free space by solving a pressure Poisson equation derived from an approximate pressure projection. Unlike earlier approaches used to simulate incompressible flows with SPH, the pressure is not a thermodynamic variable and the Courant condition is based only on fluid velocities and not on the speed of sound. Although larger time-steps can be used, the solution of the resulting elliptic pressure Poisson equation increases the total work per time-step. Efficiency comparisons show that the projection method has a significant potential to reduce the overall computational expense compared to weakly compressible SPH, particularly as the Reynolds number, Re, is increased. Simulations using this SPH projection technique show good agreement with finite-difference solutions for a vortex spin-down and Rayleigh?Taylor instability. The results, however, indicate that the use of an approximate projection to enforce incompressibility leads to error accumulation in the density field.

[1]  George Keith Batchelor,et al.  An Introduction to Fluid Dynamics. , 1969 .

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

[3]  W. Roger Briley,et al.  Numerical method for predicting three-dimensional steady viscous flow in ducts , 1974 .

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

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

[6]  Robert L. Lee,et al.  The cause and cure (!) of the spurious pressures generated by certain fem solutions of the incompressible Navier‐Stokes equations: Part 2 , 1981 .

[7]  J. Monaghan,et al.  Kernel estimates as a basis for general particle methods in hydrodynamics , 1982 .

[8]  John C. Strikwerda,et al.  Finite Difference Methods for the Stokes and Navier–Stokes Equations , 1984 .

[9]  R. Sani,et al.  On pressure boundary conditions for the incompressible Navier‐Stokes equations , 1987 .

[10]  S. Abdallah.  Numerical solutions for the pressure poisson equation with Neumann boundary conditions using a non-staggered grid, 1 , 1987 .

[11]  Joseph John Monaghan,et al.  The dynamics of interstellar cloud complexes , 1988 .

[12]  P. Colella,et al.  A second-order projection method for the incompressible navier-stokes equations , 1989 .

[13]  H. O. Nordmark,et al.  Rezoning for higher order vortex methods , 1991 .

[14]  Philip M. Gresho,et al.  On the theory of semi‐implicit projection methods for viscous incompressible flow and its implementation via a finite element method that also introduces a nearly consistent mass matrix. Part 1: Theory , 1990 .

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

[16]  J. Bell,et al.  A Second-Order Projection Method for Variable- Density Flows* , 1992 .

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

[18]  Yoshiaki Oka,et al.  A particle method for calculating splashing of incompressible viscous fluid , 1995 .

[19]  Joseph Peter Morris.  A Study of the Stability Properties of SPH , 1995 .

[20]  J. Monaghan,et al.  SPH simulation of multi-phase flow , 1995 .

[21]  S. W. Attaway,et al.  On the feasibility of using Smoothed Particle Hydrodynamics for underwater explosion calculations , 1995 .

[22]  Giovanni Lapenta,et al.  Control of the number of particles in fluid and MHD particle in cell methods , 1995 .

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

[24]  Solving the Navier Stokes equations using an exact projection method on a 2-D staggered grid , 1997 .

[25]  S. Abdallah,et al.  Dilation-Free Solutions for the Incompressible Flow Equations on Nonstaggered Grids , 1997 .

[26]  S. Koshizuka,et al.  International Journal for Numerical Methods in Fluids Numerical Analysis of Breaking Waves Using the Moving Particle Semi-implicit Method , 2022 .

[27]  Wm. G. Hoover.  Isomorphism linking smooth particles and embedded atoms , 1998 .