Coupled 1D2D models emerged as an efficient solution for a two-dimensional (2D) representation of the floodplain combined with a fast one-dimensional (1D) schematization of the main channel. At the same time, high-performance computing (HPC) has appeared as an efficient tool for model acceleration. In this work, a previously validated 1D2D Central Processing Unit (CPU) model is combined with an HPC technique for fast and accurate flood simulation. Due to the speed of 1D schemes, a hybrid CPU/GPU model that runs the 1D main channel on CPU and accelerates the 2D floodplain with a Graphics Processing Unit (GPU) is presented. Since the data transfer between sub-domains and devices (CPU/GPU) may be the main potential drawback of this architecture, the test cases are selected to carry out a careful time analysis. The results reveal the speed-up dependency on the 2D mesh, the event to be solved and the 1D discretization of the main channel. Additionally, special attention must be paid to the time step size computation shared between sub-models. In spite of the use of a hybrid CPU/GPU implementation, high speed-ups are accomplished in some cases.

The vast majority of reservoirs, although built for irrigation and water supply purposes, are also used as regulation tools during floods in river basins. Thus, the selection of the most suitable model when facing the simulation of a flood wave in a combination of river reach and reservoir is not direct and frequently some analysis of the proper system of equations and the number of solved flow velocity components is needed. In this work, a stretch of the Ebro River (Spain), which is the biggest river in Spain, is simulated solving the Shallow Water Equations (SWE). The simulation model covers the area of river between the city of Zaragoza and the Mequinenza dam. The domain encompasses 721.92 km${}^2$ with 221 km of river bed, of which the last 75 km belong to the Mequinenza reservoir. The results obtained from a one-dimensional (1D) model are validated comparing with those provided by a two-dimensional (2D) model based on the same numerical scheme and with measurements. The 1D modelling loses the detail of the floodplain, but nevertheless the computational consumption is much lower compared to the 2D model with a permissible loss of accuracy. Additionally, the particular nature of this reservoir might turn the 1D model into a more suitable option. An alternative technique is applied in order to model the reservoir globally by means of a volume balance (0D) model, coupled to the 1D model of the river (1D-0D model). The results obtained are similar to those provided by the full 1D model with an improvement on computational time. Finally, an automatic regulation is implemented by means of a Proportional-Integral-Derivative (PID) algorithm and tested in both the full 1D model and the 1D-0D model. The results show that the coupled model behaves correctly even when controlled by the automatic algorithm.

The computational simulation of rivers is a useful tool that can be applied in a wide range of situations from providing real time alerts to the design of future mitigation plans. However, for all the applications, there are two important requirements when modeling river behavior: accuracy and reasonable computational times. This target has led to recent developments in numerical models based on the full two-dimensional (2D) shallow water equations (SWE). This work presents a GPU accelerated 2D SW model for the simulation of flood events in real time. It is based on a well-balanced explicit first-order finite volume scheme able to run over dry beds without the numerical instabilities that are likely to occur when used in complex topography. The model is applied to reproduce a real event in the reach of the Ebro River (Spain) with a downstream reservoir, in which a study of the most appropriate boundary condition (BC) for modeling of the dam is assessed (time-dependent level condition and weir condition). The whole creation of the model is detailed in terms of mesh optimization and validation. The simulation results are compared with field data over the flood duration (up to 20 days), allowing an analysis of the performance and time saved by different GPU devices and with the different BCs. The high values of fit between observed and simulated results, as well as the computational times achieved, are encouraging to propose the use of the model as a forecasting system.