TY - JOUR

T1 - A simple procedure to improve the pressure evaluation in hydrodynamic context using the SPH

AU - Molteni, Diego

AU - Colagrossi, Andrea

PY - 2009

Y1 - 2009

N2 - In literature, it is well know that the Smoothed Particle Hydrodynamics method can be affected by numerical noise on the pressure field when dealing with liquids. This can be highly dangerous when an SPH code is dynamically coupled with a structural solver. In this work a simple procedure is proposed to improve the computation of the pressure distribution in the dynamics of liquids. Such a procedure is based on the use of a density diffusion term in the equation for the mass conservation. This diffusion is a pure numerical effect, similar to the well known artificial viscosity originally proposed in SPH method to smooth out the shock discontinuities. As the artificial viscosity, the density diffusion used here goes to zero increasing the number of particles recovering consistency and convergence of the final numerical scheme adopted. Different artificial density diffusion formulas have been studied, paying attention to prevent unphysical changes of the flows. To show the improvements of the new scheme proposed here, a suitable set of examples, for which reference solutions or experimental data are available, has been tested.

AB - In literature, it is well know that the Smoothed Particle Hydrodynamics method can be affected by numerical noise on the pressure field when dealing with liquids. This can be highly dangerous when an SPH code is dynamically coupled with a structural solver. In this work a simple procedure is proposed to improve the computation of the pressure distribution in the dynamics of liquids. Such a procedure is based on the use of a density diffusion term in the equation for the mass conservation. This diffusion is a pure numerical effect, similar to the well known artificial viscosity originally proposed in SPH method to smooth out the shock discontinuities. As the artificial viscosity, the density diffusion used here goes to zero increasing the number of particles recovering consistency and convergence of the final numerical scheme adopted. Different artificial density diffusion formulas have been studied, paying attention to prevent unphysical changes of the flows. To show the improvements of the new scheme proposed here, a suitable set of examples, for which reference solutions or experimental data are available, has been tested.

KW - Convergence tests

KW - Fluid–structure impact problems

KW - Free surface flows

KW - SPH pressure evaluation

KW - Smoothed Particle Hydrodynamics

KW - Weak-compressibility

KW - Convergence tests

KW - Fluid–structure impact problems

KW - Free surface flows

KW - SPH pressure evaluation

KW - Smoothed Particle Hydrodynamics

KW - Weak-compressibility

UR - http://hdl.handle.net/10447/35057

M3 - Article

VL - 180

SP - 861

EP - 872

JO - Computer Physics Communications

JF - Computer Physics Communications

SN - 0010-4655

ER -