On the Convergence of SPH Method for Scalar Conservation Laws with Boundary Conditions

Abstract

This paper is the third of a series where the convergence analysis of SPH method for multidimensional conservation laws is analyzed. In this paper, two original numerical models for the treatment of boundary conditions are elaborated. To take into account nonlinear effects in agreement with Bardos, LeRoux and Nedelec boundary conditions ([1], [14]), the state at the boundary is computed by solving appropriate Riemann problems. The first numerical model is developed around the idea of boundary forces in surrounding walls, recently initiated in [33] by Monaghan in his simulation of gravity currents. The second one extends the well-known approach of ghost particles for plane boundaries to the case of general curved boundaries. The convergence analysis in L p loc (p < ∞) is achieved thanks to the uniqueness result of measure-valued solutions recently established in [3] for L ∞ initial and boundary data.