Abstract
A large Courant–Friedrichs–Lewy (CFL) algorithm is presented for the explicit, finite volume solution of hyperbolic systems of conservation laws, with a focus on the shallow water equations. The Riemann problems used in the flux computation are determined using averaging kernels that extend over several computational cells. The usual CFL stability constraint is replaced with a constraint involving the kernel support size. This makes the method unconditionally stable with respect to the size of the computational cells, allowing the computational mesh to be refined locally to an arbitrary degree without altering solution stability. The practical implementation of the method is detailed for the shallow water equations with topographical source term. Computational examples report applications of the method to the linear advection, Burgers and shallow water equations. In the case of sharp bottom discontinuities, the need for improved, well-balanced discretisations of the geometric source term is acknowledged.