A High-Order Well-Balanced Positivity-Preserving Moving Mesh DG Method for the Shallow Water Equations With Non-Flat Bottom Topography

Abstract

A rezoning-type adaptive moving mesh discontinuous Galerkin method is proposed for the numerical solution of the shallow water equations with non-flat bottom topography. The well-balance property is crucial to the simulation of perturbation waves over the lake-at-rest steady state such as waves on a lake or tsunami waves in the deep ocean. To ensure the well-balance and positivity-preserving properties, strategies are discussed in the use of slope limiting, positivity-preservation limiting, and data transferring between meshes. Particularly, it is suggested that a DG-interpolation scheme be used for the interpolation of both the flow variables and bottom topography from the old mesh to the new one and after each application of the positivity-preservation limiting on the water depth, a high-order correction be made to the approximation of the bottom topography according to the modifications in the water depth. Moreover, mesh adaptation based on the equilibrium variable and water depth is shown to give more desirable results than that based on the commonly used entropy function. Numerical examples in one and two spatial dimensions are presented to demonstrate the well-balance and positivity-preserving properties of the method and its ability to capture small perturbations of the lake-at-rest steady state.