A weakly compressible SPH method based on a low-dissipation Riemann solver

We present a low-dissipation weakly-compressible SPH method for modeling free-surface flows exhibiting violent events such as impact and breaking. The key idea is to modify a Riemann solver which determines the interaction between particles by a simple limiter to decrease the intrinsic numerical dissipation. The modified Riemann solver is also extended for imposing wall boundary conditions. Numerical tests show that the method resolves free-surface flows accurately and produces smooth, accurate pressure fields. The method is compatible with the hydrostatic solution and exhibits considerably less numerical damping of the mechanical energy than previous methods.