Finite volume model for two-dimensional shallow water flows on unstructured grids

A numerical model based upon a second-order upwind finite volume method on unstructured triangular grids is developed for solving shallow water equations. The HLL approximate Riemann solver is used for the computation of inviscid flux functions, which makes it possible to handle discontinuous solutions. A multidimensional slope-limiting technique is employed to achieve second-order spatial accuracy and to prevent spurious oscillations. To alleviate the problems associated with numerical instabilities due to small water depths near a wet/dry boundary, the friction source terms are treated in a fully implicit way. A third-order total variation diminishing Runge-Kutta method is used for the time integration of semidiscrete equations. The developed numerical model has been applied to several test cases as well as to real flows. Numerical tests prove the robustness and accuracy of the model.