Towards consistence and convergence of conservative SPH approximations

Abstract

Typical conservative smoothed particle hydrodynamics (SPH) approximations of the gradient of a scalar field introduce two errors: one (smoothing error) is due to smoothing of the gradient by an integration associated with a kernel function; the other (integration error) is due to approximating the integration by summation over all particles within the kernel support. When particles are not on a uniform grid, the integration error leads to violation of zero-order consistency, i.e. the inability to reproduce a constant field. In this paper we confirm that partition of unity is the condition under which the conservative SPH approximation achieves both consistence and convergence. We show that this condition can be achieved by relaxing a particle distribution under a constant pressure field and invariant particle volume. The resulting particle distribution is very similar to that is typical for liquid molecules. We further show that with two different typical kernel functions the SPH approximation, upon satisfying the partition of unity property, is able to achieve very high-order of the integration error, which previously could be shown only with particles on a uniform grid. The background pressure used in a weakly compressible SPH simulation implies a self-relaxation mechanism, which explains that convergence with respect to increasing particle numbers could be obtained in SPH simulations, although not predicted by previous numerical analysis. Furthermore, by relating the integration error to the background pressure, we explain why the previously proposed transport-velocity formulation of SPH (S. Adami et al. (2013) [1]) is able to achieve unprecedented accuracy and stability.