A 2D Hydraulic Simulation Model Including Dynamic Piping and Overtopping Dambreach

15 Pags.- 12 Figs.- 3 Tabls. © 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license.Numerical simulation of unsteady free surface flows using depth averaged equations that consider the presence of initial discontinuities has been often reported for situations dealing with dam break flow. The usual approach is to assume a sudden removal of the gate at the dam location. Additionally, in order to prevent any kind of dam risk in earthen dams, it is very interesting to include the possibility of a progressive dam breach leading to dam overtopping or dam piping so that predictive hydraulic models benefit the global analysis of the water flow. On the other hand, when considering a realistic large domain with complex topography, a fine spatial discretization is mandatory. Hence, the number of grid cells is usually very large and, therefore, it is necessary to use parallelization techniques for the calculation, with the use of Graphic Processing Units (GPU) being one of the most efficient, due to the leveraging of thousands of processors within a single device. The aim of the present work is to describe an efficient GPU-based 2D shallow water flow solver (RiverFlow2D-GPU) supplied with the formulation of internal boundary conditions to represent dynamic dam failure processes. The results obtained indicate that it is able to develop a transient flow analysis taking into account several scenarios. The efficiency of the model is proven in two complex domains, leading to >76× faster simulations compared with the traditional CPU computation.This work was partially funded by the University of Zaragoza under the research project JIUZ2022-IAR-03: High Performance Computing Tools for Smart Agriculture and ForestryPeer reviewe