DualSPHysics+: An enhanced DualSPHysics with improvements in accuracy, energy conservation and resolution of the continuity equation

To further reduce both the instantaneous and accumulated errors in velocity divergence, a Hyperbolic/Parabolic Divergence Cleaning (HPDC) scheme is incorporated in addition to the VEM scheme. The implementations of the introduced schemes on both CPU and GPU-based versions of the DualSPHysics+ code along with details on the compilation, running and computational performance are presented. Validations in terms of accuracy, energy conservation and convergence of DualSPHysics+ are shown via several relevant benchmarks.

It is demonstrated that a better velocity divergence error cleaning in both instantaneous and accumulated errors can be achieved by the combination of VEM and HPDC. Meanwhile, the excessive energy dissipation by the artificial viscosity is shown to be suppressed by adopting the Riemann stabilization term. Enhanced resolution of the continuity equation along with improved energy conservation of DualSPHysics+ advance the SPH-based simulation of incompressible free-surface fluid flows.