Truncation error in mesh-free particle methods
Quinlan, Nathan J. Basa, Mihai Lastiwka, Martin
Quinlan, Nathan J.
Basa, Mihai
Lastiwka, Martin
Loading...
Identifiers
http://hdl.handle.net/10379/1170
https://doi.org/10.13025/19145
https://doi.org/10.13025/19145
Repository DOI
Publication Date
2006-06-25
Type
Article
Downloads
Citation
Quinlan, N. J., Basa, M., & Lastiwka, M. (2006). Truncation error in mesh-free particle methods. International Journal for Numerical Methods in Engineering, 66(13), 2064-2085.
Abstract
A truncation error analysis has been developed for the approximation of spatial derivatives in smoothed particle hydrodynamics (SPH) and related first-order consistent methods such as the first-order form of the reproducing kernel particle method. Error is shown to depend on both the smoothing length h and the ratio of particle spacing to smoothing length, x/h. For uniformly spaced particles in one dimension, analysis shows that as h is reduced while maintaining constant x/h, error decays as h2 until a limiting discretization error is reached, which is independent of h. If x/h is reduced while maintaining constant h (i.e. if the number of neighbours per particle is increased), error decreases at a rate which depends on the kernel function's smoothness. When particles are distributed non-uniformly, error can grow as h is reduced with constant x/h. First-order consistent methods are shown to remove this divergent behaviour. Numerical experiments confirm the theoretical analysis for one dimension, and indicate that the main results are also true in three dimensions. This investigation highlights the complexity of error behaviour in SPH, and shows that the roles of both h and x/h must be considered when choosing particle distributions and smoothing lengths.
Funder
Publisher
Publisher DOI
Rights
Attribution-NonCommercial-NoDerivs 3.0 Ireland