Modelling Dam Break Evolution over a Wet Bed with Smoothed Particle Hydrodynamics: A Parameter Study

Abstract

When investigating water flow in spillways and energy dissipation, it is important to know the behavior of the free surfaces. To capture the real dynamic behavior of the free surfaces is therefore crucial when performing simulations. Today, there is a lack in the possibility to model such phenomenon with traditional methods. Hence, this work focuses on a parameter study for one alternative simulation tool available, namely the meshfree, Lagrangian particle method Smoothed Particle Hydrodynamics (SPH). The parameter study includes the choice of equation-of-state (EOS), the artificial viscosity constants, using a dynamic versus a static smoothing length, SPH particle spatial resolution and the finite element method (FEM) mesh scaling of the boundaries. The two dimensional SPHERIC Benchmark test case of dam break evolution over a wet bed was used for comparison and validation. The numerical results generally showed a tendency of the wave front to be ahead of the experimental results, i.e. to have a greater wave front velocity. The choice of EOS, FEM mesh scaling as well as using a dynamic or a static smoothing length showed little or no significant effect on the outcome, though the SPH particle resolution and the choice of artificial viscosity constants had a major impact. A high particle resolution increased the number of flow features resolved for both choices of artificial viscosity constants, but at the expense of increasing the mean error. Furthermore, setting the artificial viscosity constants equal to unity for the coarser cases resulted in a highly viscous and unphysical solution, and thus the relation between the artificial viscosity constants and the particle resolution and its impact on the behavior of the fluid needed to be further investigated.