A comprehensive study on the parameters setting in smoothed particle hydrodynamics (SPH) method applied to hydrodynamics problems

Smoothed particle hydrodynamics (SPH) is a meshfree, Lagrangian particle method which has advantages in handling solids with extremely large deformation. Like any other numerical methods, cares must be taken to ensure its desirable accuracy and stability through considering several correction techniques in calculation. The selection of values for parameters in those correction approaches is a key step in SPH simulation, which is always difficult for new beginners to deal well with effectively. This paper examines the common inconsistency and instability problems in SPH method and studies its computational efficiency when applied to hydrodynamics problems with material strength like soil column collapse. We analyzed in detail how the correction techniques mitigate these inconsistency and instability problems. Also, the numerical testing results associate with different values for the parameters used in the correction techniques are provided for better understanding the influence of these parameters and for finding out the desirable values. It is found that (1) the SPH method is easily subjected to an inconsistency problem in the boundary area due to the boundary deficiency, and it can be treated well by adopting “virtual particles” contributing to the particle summations. (2) The numerical oscillation in SPH simulation can be mitigated effectively by artificial viscosity with the suggested parameter values. (3) The tension cracking treatment, artificial viscosity and artificial stress work well in removing the tensile instability problem in SPH method. In addition, the nearest neighboring particle searching (NNPS) algorithm, spacing ratio, smoothing length and time step influence the efficiency and accuracy of SPH method significantly. It is shown that SPH method with suggested parameters values can produce a very good result compared with the experimental result.