A Numerical Model for Two-Dimensional Flood Routing in Complex Terrains

In this paper, a new powerful numerical hydrodynamic in-house model is presented and tested. The model simulates flood routing in two dimensions. It is based on the solution of Shallow Water Equations using the Finite Difference Method according to the explicit McCormack numerical scheme which has shock capturing capability. The innovation of the proposed model lies in the modification of McCormack scheme by incorporating artificial viscosity through a diffusion factor in order to avoid oscillations as proposed by various researchers. Additionally, a threshold of water depth is introduced in order to distinguish the wet and dry cells of the computational domain. The model is capable of producing maps for the inundation extent, water depths and depth-averaged water velocities. Finally, the paper presents extensive testing of the model by comparison with analytical solution, experimental results and with the output of another software package in real world flood simulation studies.