We study droplet impact problems in a 3D cylindrical or equivalent 2D Cartesian geometry. Such structures do have an approximate experimental realization, and they are often simulated as a testbed for computational methods. We focus on droplet impact on a smooth homogeneous surface as well as head-on collision of two droplets. We perform an energy-budget analysis and introduce a correlation, which predicts the maximum spreading radius as a function of Reynolds number and Weber number. We show how the dissipation term in this analysis can be decomposed into boundary-layer dissipation in the droplet lamella (where applicable) and head loss. We use existing results in the literature (simulations and experiments) as well as our own simulation results to validate the correlation. Dissipation by head loss is a key term in the analysis: only by modeling it accurately, one can obtain good agreement between the simulations and the theory.

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.