An Analysis of Accuracy in One-Dimensional Smoothed Particle Hydrodynamics
17th AIAA Computational Fluid Dynamics Conference, 4622, AIAA, (2005)DOI: 10.2514/6.2005-4622
Abstract
A truncation error analysis has been developed for the approximation of spatial derivatives in one-dimensional Smoothed Particle Hydrodynamics (SPH). For this purpose, the SPH interpolation is understood as a two-step process of smoothing and discretisation. As smoothing length is reduced while maintaining a constant ratio of particle spacing Δx to smoothing length h, error decays as h2; however, there is a finite limiting discretisation error. If particle spacing Δx is reduced while holding a constant h, error decreases at a rate which depends on the kernel function's smoothness at its boundaries. When particles are distributed non-uniformly, error behaviour is complex, and discretisation error can diverge as h is reduced. A first-order consistent SPH method is shown to remove this behaviour. The findings of the theoretical analysis are confirmed by numerical experiments.