Simulation of shallow water flows with shoaling areas and bottom discontinuities

A numerical method based on a second-order accurate Godunov-type scheme is described for solving the shallow water equations on unstructured triangular-quadrilateral meshes. The bottom surface is represented by a piecewise linear approximation with discontinuities, and a new approximate Riemann solver is used to treat the bottom jump. Flows with a dry sloping bottom are computed using a simplified method that admits negative depths and preserves the liquid mass and the equilibrium state. The accuracy and performance of the approach proposed for shallow water flow simulation are illustrated by computing one- and two-dimensional problems.