A weakly compressible SPH method with WENO reconstruction

Abstract

In this work, we improve a weakly-compressible SPH method which is based on our previous work of Zhang et al. (2017) [30] by employing a weighted essentially non-oscillatory (WENO) reconstruction. The key idea is to construct along each interacting particle pair a 4-point stencil and to apply a WENO reconstruction for determining the particle interaction with a low-dissipation Riemann solver. Several numerical examples on Taylor-Green vortex flow, dam break and non-linear liquid sloshing demonstrate that the method preserves the capability of producing smooth and accurate pressure fields of the original method and now achieves also very small numerical dissipation.