A high-order well-balanced positivity-preserving moving mesh DG method for the shallow water equations with non-flat bottom topography

TL;DR

This paper introduces a high-order, well-balanced, positivity-preserving moving mesh discontinuous Galerkin method for accurately simulating shallow water equations with complex bottom topography, ensuring stability and precision in perturbation wave modeling.

Contribution

It develops a novel adaptive moving mesh DG method that maintains well-balance and positivity, with advanced data transfer and correction strategies for shallow water simulations.

Findings

Method effectively preserves lake-at-rest steady state.

Accurately captures small perturbation waves.

Demonstrates superior mesh adaptation based on equilibrium variables.

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.