Numerical study of PPE source term errors in the incompressible SPH models

Abstract

In a recent work by Gui et al. , an incompressible SPH model was presented that employs a mixed pressure Poisson equation (PPE) source term combining both the density‐invariant and velocity divergence‐free formulations. The present work intends to apply the model to a wider range of fluid impact situations in order to quantify the numerical errors associated with different formulations of the PPE source term in incompressible SPH (ISPH) models. The good agreement achieved between the model predictions and the documented data is taken as a further demonstration that the mixed source term formulation can accurately predict the fluid impact pressures and forces, both in the magnitude and in the spatial and temporal patterns. Furthermore, an in‐depth numerical analysis using either the pure density‐invariant or velocity divergence‐free formulation has revealed that the pure density‐invariant formulation can lead to relatively large divergence errors while the velocity divergence‐free formulation may cause relatively large density errors. As compared with these two approaches, the mixed source term formulation performs much better having the minimum total errors in all test cases. Although some recent studies found that the weakly compressible SPH models perform somewhat better than the incompressible SPH models in certain fluid impact problems, we have shown that this could be largely caused by the particular formulation of PPE source term in the previous ISPH models and a better formulation of the source term can significantly improve the accuracy of ISPH models. Copyright © 2014 John Wiley & Sons, Ltd.

On-Bottom Stability Analysis of Cylinders under Tsunami-Like Solitary Waves

“…More complex models as compared to Morison-type ones and dealing with an improved description of flow-cylinder interaction processes have been also developed [20,21]. More recently, Aristodemo et al [22] performed a laboratory study on non-breaking solitary wave forces with respect to a horizontal cylinder placed at half-water depth that was supported by numerical simulations based on the smoothed particle hydrodynamics technique (e.g., [23][24][25][26][27][28][29]). In this case, the effect of the free surface was negligible, i.e., there was no scattering, and that related to the bottom was weak.…”

Improved Incompressible Material Point Method Based on Particle Density Correction

“…A better scheme in which a hybrid PPE source including both the velocity divergence term and relative density variance term is proposed by Asai et al [2012]. Numerical results show that the scheme can significantly improve the distribution of particles and velocity field accuracy [Gui et al (2013[Gui et al ( , 2015]. This hybrid idea has also been applied in the MPS method [Tanaka and Masunaga (2010); Khayyer and Gotoh (2011)].…”

A Symmetric Particle-Based Simulation Scheme towards Large Scale Diffuse Fluids

“…This method was first proposed by Cummins et al [23] in 1999, projecting the velocity field onto a divergence-free space by solving a pressure Poisson equation. To improve the ISPH modeling capacity, Gui et al [24] combined the standard density-invariant and the velocity divergence-free formulations in a weighted average form and proposed a mixed source term model. This model can obtain more accurate impact pressure and force as compared with the results obtained by using either the density-invariant or the velocity divergence-free ISPH model.…”