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.

This paper presents a Lagrangian model for particle transport driven by a two-dimensional (2D) shallow water model, assuming that the particles have negligible mass and volume, are located at the free surface, and without interactions between them. Particle motion is based on advection and turbulent diffusion, which is added using a random-walk model. The equations for particle advective transport are solved using the flow velocity provided by a 2D shallow water solver and an online first-order Euler method, an online fourth-order Runge–Kutta method and an offline fourth-order Runge–Kutta method. The primary objective of this work is to present the capabilities of the new Lagrangian particle transport (LPT) model, while also providing an analysis of the accuracy and computational efficiency of the numerical schemes and their implementation for particle transport. To verify the accuracy and computational cost, several test cases inspired by laboratory set-ups are simulated. In this analysis, the Euler online method provides the best compromise between accuracy and computational efficiency. Finally, a localized precipitation event in the Arnás catchment is simulated to test the model's capability to represent particle transport in overland flow over irregular topography.