A Boundary Integral SPH Formulation: Consistency and Applications to ISPH and WCSPH

Abstract

One of the historical problems appearing in SPH formulations is the inconsistencies coming from the inappropriate implementation of boundary conditions. In this work, this problem has been investigated; instead of using typical methodologies such as extended domains with ghost or dummy particles where severe inconsistencies are found, we included the boundary terms that naturally appear in the formulation. First, we proved that in the 1D smoothed continuum formulation, the inclusion of boundary integrals allows for a consistent O(h) formulation close to the boundaries. Second, we showed that the corresponding discrete version converges to a certain solution when the discretization SPH parameters tend to zero. Typical tests with the first and second derivative operators confirm that this boundary condition implementation works consistently. The 2D Poisson problem, typically used in ISPH, was also studied, obtaining consistent results. For the sake of completeness, two practical applications, namely, the duct flow and a sloshing tank, were studied with the results showing a rather good agreement with former experiments and previous results.