SPH for incompressible free-surface flows. Part I: Error analysis of the basic assumptions

• Published in: Computers & fluids Vol. 86; pp. 611 - 624
• Main Authors: Kiara, Areti, Hendrickson, Kelli, Yue, Dick K.P.
• Format: Journal Article
• Language: English
• Published: Elsevier Ltd 05.11.2013
• Subjects: Acoustics; Compressibility; Computational fluid dynamics; Error analysis; Fluid flow; Free-surface flows; High-frequency oscillations; Instability; Interpolation; Kernel interpolation; Kernels; Mathematical analysis; SPH; Weak compressibility assumption; High-frequency oscillations; Instability; Kernel interpolation; SPH; Weak compressibility assumption; Free-surface flows
• ISSN: 0045-7930; 1879-0747
• DOI: 10.1016/j.compfluid.2013.05.023

•We provide a systematic analysis of the assumptions fundamental to SPH.•We identify and quantify the high frequency oscillations in SPH dynamics.•We analyze consistency of kernel interpolation in non-uniform particle distribution.•We identify instability mechanisms due to weak incompressibility.•We demonstrate our findings in the illustrative case of a hydrostatic problem. We provide a quantitative error analysis of the fundamental assumptions of weak compressibility and kernel interpolation in SPH for incompressible free-surface flows. This provides a framework for understanding many of the existing SPH variations. Weak compressibility generates spurious high-frequency acoustic solutions. We quantify acoustic eigensolutions, demonstrating that these can dominate the dynamics (more so than the kinematics) in SPH. Kernel interpolation consistency is affected when particle distributions are non-uniform and near boundaries. We show that for non-uniform distributions obtained from the SPH governing equations, second order consistency is retained. Near the free surface, however, kernel interpolation implicitly imposes local spurious accelerations. These, in the presence of weak compressibility, are manifest in high-frequency acoustic dynamics of the order of gravity g throughout the domain. We identify instability mechanisms in the presence of density non-uniformities that result in unstable temporal growth of both depth decaying and depth oscillatory modes. The latter are also identified as spurious high-frequency SPH errors. We illustrate all of these findings in the simplest case of a hydrostatic problem. Finally, we discuss the performance of existing treatments to SPH in the context of the present findings.