A new insight into the consistency of smoothed particle hydrodynamics

Abstract

In this paper the problem of consistency of smoothed particle hydrodynamics (SPH) is solved. A novel error analysis is developed in $n$-dimensional space using the Poisson summation formula, which enables the treatment of the kernel and particle approximation errors in combined fashion. New consistency integral relations are derived for the particle approximation which correspond to the cosine Fourier transform of the classically known consistency conditions for the kernel approximation. The functional dependence of the error bounds on the SPH interpolation parameters, namely the smoothing length $h$ and the number of particles within the kernel support $\mathcal{N}$ is demonstrated explicitly from which consistency conditions are seen to follow naturally. As $\mathcal{N}\to \infty$, the particle approximation converges to the kernel approximation independently of $h$ provided that the particle mass scales with $h$ as $m\propto h^{\beta }$, with $\beta >n$. This implies that as $h\to 0$, the joint limit $m\to 0$, $\mathcal{N}\to \infty$, and $N\to \infty$ is necessary for complete convergence to the continuum, where $N$ is the total number of particles. The analysis also reveals the presence of a dominant error term of the form $(\ln \mathcal{N}{)}^n/\mathcal{N}$, which tends asymptotically to $1/\mathcal{N}$ when $\mathcal{N}\gg 1$, as it has long been conjectured based on the similarity between the SPH and the quasi-Monte Carlo estimates.

27 pages. v2: Matches published version

References in corpus (6)

• Improving convergence in smoothed particle hydrodynamics simulations without pairing instability
• Astrophysical Smooth Particle Hydrodynamics
• Resolving mixing in Smoothed Particle Hydrodynamics
• Numerical Convergence in Smoothed Particle Hydrodynamics
• Particle Methods in Astrophysical Fluid Dynamics
• Recent results in the systematic derivation and convergence of SPH

Cited by in corpus (5)

• Improving stability of moving particle semi-implicit method by source terms based on time-scale correction of particle-level impulses
• Impetus: Consistent SPH calculations of 3D spherical Bondi accretion onto a black hole
• REMIX SPH -- improving mixing in smoothed particle hydrodynamics simulations using a generalised, material-independent approach
• SPH modeling of natural convection around a heated horizontal cylinder: A comparison with experiments
• Multi-scale physics of cryogenic liquid helium-4: Inverse coarse-graining properties of smoothed particle hydrodynamics