Theoretical analysis of SPH in simulating free-surface viscous flows

A theoretical analysis on the performance, close to a free surface, of the most used SPH formulations for Newtonian viscous terms is carried out in this paper. After an introduction of the SPH formalism, the SPH expressions for the viscous term in the momentum equation are analyzed in their continuous form. Using a Taylor expansion, a reformulation of those expressions is undertaken which allows to characterize the behavior of the viscous term close to the free surface. Under specific flow conditions, we show that the viscous term close to the free surface is singular when the spatial resolution is increased. This problem is in essence related to the incompleteness of the kernel function close to the free surface and appears for all the formulations considered. In order to assess the impact of such singular behavior, an analysis of the global energy dissipation is carried out, which shows that such a free-surface singularity vanishes when the integral quantities are considered. Notwithstanding that, not all the SPH viscous formulas allow the correct evaluation of the energy dissipation rate and, consequently, they may lead to an inaccurate modelling of viscous free-surface flows.