2016/05/14 by Leonardo Di G. Sigalotti, Sigalotti, Leonardo Di G., Jaime Klapp +7
Engineering · #Astrophysics of Galaxies (astro-ph.GA) #Computational Physics (physics.comp-ph) #Cosmology and Nongalactic Astrophysics (astro-ph.CO) #FOS: Mathematics #FOS: Physical sciences #Fluid Dynamics (physics.flu-dyn) #Fluid Dynamics Simulations and Interactions #Fluid Dynamics and Vibration Analysis #Numerical Analysis (math.NA) #Numerical methods in engineering
paper · pdf · doi:10.48550/arxiv.1605.05245
openalex publication_date 2016/05/14 · openalex created_date 2022/08/12 · openalex updated_date 2026/07/28
The problem of consistency of smoothed particle hydrodynamics (SPH) has demanded considerable attention in the past few years due to the ever increasing number of applications of the method in many areas of science and engineering. A loss of consistency leads to an inevitable loss of approximation accuracy. In this paper, we revisit the issue of SPH kernel and particle consistency and demonstrate that SPH has a limiting second-order convergence rate. Numerical experiments with suitably chosen test functions validate this conclusion. In particular, we find that when using the root mean square error as a model evaluation statistics, well-known corrective SPH schemes, which were thought to converge to second, or even higher order, are actually first-order accurate, or at best close to second order. We also find that observing the joint limit when N→∞, h→ 0, and n→∞, as was recently proposed by Zhu et al., where N is the total number of particles, h is the smoothing length, and n is the number of neighbor particles, standard SPH restores full C0 particle consistency for both the estimates of the function and its derivatives and becomes insensitive to particle disorder.