Entropy Stable Discontinuous Galerkin Schemes on Moving Meshes with Summation-by-Parts Property for Hyperbolic Conservation Laws

TL;DR

This paper develops an entropy stable discontinuous Galerkin scheme on moving meshes for hyperbolic conservation laws, extending entropy analysis and flux functions to dynamic meshes, and validating with numerical experiments.

Contribution

It introduces a novel entropy analysis framework for moving mesh DG schemes, including entropy conservative fluxes and stability proofs, for hyperbolic PDEs.

Findings

Entropy conservative flux functions for shallow water and Euler equations are constructed.

The scheme is proven to be entropy conservative and entropy stable on moving meshes.

Numerical experiments confirm theoretical properties and stability.

Abstract

This work is focused on the entropy analysis of a semi-discrete nodal discontinuous Galerkin spectral element method (DGSEM) on moving meshes for hyperbolic conservation laws. The DGSEM is constructed with a local tensor-product Lagrange-polynomial basis computed from Legendre-Gauss-Lobatto (LGL) points. Furthermore, the collocation of interpolation and quadrature nodes is used in the spatial discretization. This approach leads to discrete derivative approximations in space that are summation-by-parts (SBP) operators. On a static mesh, the SBP property and suitable two-point flux functions, which satisfy the entropy condition from Tadmor, allow to mimic results from the continuous entropy analysis on the discrete level. In this paper, Tadmor's condition is extended to the moving mesh framework. Based on the moving mesh entropy condition, entropy conservative two-point flux functions for the homogeneous shallow water equations and the compressible Euler equations are constructed. Furthermore, it will be proven that the semi-discrete moving mesh DGSEM is an entropy conservative scheme when a two-point flux function, which satisfies the moving mesh entropy condition, is applied in the split form DG framework. This proof does not require any exactness of quadrature in the spatial integrals of the variational form. Nevertheless, entropy conservation is not sufficient to tame discontinuities in the numerical solution and thus the entropy conservative moving mesh DGSEM is modified by adding numerical dissipation matrices to the entropy conservative fluxes. Then, the method becomes entropy stable such that the discrete mathematical entropy is bounded at any time by its initial and boundary data when the boundary conditions are specified appropriately. The theoretical properties of the moving mesh DGSEM will be validated by numerical experiments for the compressible Euler equations.