SPH particle boundary forces for arbitrary boundaries

This paper is concerned with approximating arbitrarily shaped boundaries in SPH simulations. We model the boundaries by means of boundary particles which exert forces on a fluid. We show that, when these forces are chosen correctly, and the boundary particle spacing is a factor of 2 (or more) less than the fluid particle spacing, the total boundary force on a fluid SPH particle is perpendicular to boundaries with negligible error. Furthermore, the variation in the force as a fluid particle moves, while keeping a fixed distance from the boundary, is also negligible. The method works equally well for convex or concave boundaries. The new boundary forces simplify SPH algorithms and are superior to other methods for simulating complicated boundaries. We apply the new method to (a) the rise of a cylinder contained in a curved basin, (b) the spin down of a fluid in a cylinder, and (c) the oscillation of a cylinder inside a larger fixed cylinder. The results of the simulations are in good agreement with those obtained using other methods, but with the advantage that they are very simple to implement.