Adaptive physical-constraints-preserving unstaggered central schemes for shallow water equations on quadrilateral meshes

Abstract

A well-balanced and positivity-preserving adaptive unstaggered central scheme for two-dimensional shallow water equations with nonflat bottom topography on irregular quadrangles is presented. The irregular quadrilateral mesh adds to the difficulty of designing unstaggered central schemes. In particular, the integral of the source term needs to subtly be dealt with. A new method of discretizing the source term for the well-balanced property is proposed, which is one of the main contributions of this work. The spacial second-order accuracy is obtained by constructing piecewise bilinear functions. Another novelty is that we introduce a strong-stability-preserving \emph{Unstaggered-Runge-Kutta} method to improve the accuracy in time integration. Adaptive moving mesh strategies are introduced to couple with the current unstaggered central scheme. The well-balanced property is still valid. The positivity-preserving property can be proved when the cells shrink. We prove that the current adaptive unstaggered central scheme can obtain the stationary solution (``lake at rest" steady solutions) and guarantee the water depth to be nonnegative. Several classical problems of shallow water equations are shown to demonstrate the properties of the current numerical scheme.