Generalised and efficient wall boundary condition treatment in GPU-accelerated smoothed particle hydrodynamics

Abstract: This paper presents a generalised and efficient wall boundary treatment in the smoothed particle hydrodynamics (SPH) method for 3-D complex and arbitrary geometries with single- and multi-phase flows to be executed on graphics processing units (GPUs). Using a force balance between the wall and fluid particles with a novel penalty method, a pressure boundary condition is applied on the wall dummy particles which effectively prevents non-physical particle penetration into the wall boundaries also in highly violent impacts and multi-phase flows with high density ratios. A new density reinitialisation scheme is also presented to enhance the accuracy. The proposed method is very simple in comparison with previous wall boundary formulations on GPUs that enforces no additional memory caching and thus is ideally suited for heterogeneous architectures of GPUs. The method is validated in various test cases involving violent single- and multi-phase flows in arbitrary geometries and demonstrates very good robustness, accuracy and performance. The new wall boundary condition treatment is able to improve the high accuracy of its previous version \citep{ADAMI2012wall} also in complex 3-D and multi-phase problems, while it is efficiently executable on GPUs with single precision floating points arithmetic which makes it suitable for a wide range of GPUs, including consumer graphic cards. Therefore, the method is a reliable solution for the long-lasting challenge of the wall boundary condition in the SPH method for a broad range of natural and industrial applications.