Gero Schnücke

Lucas Friedrich, Gero Schnücke, Andrew R. Winters et al.

This work examines the development of an entropy conservative (for smooth solutions) or entropy stable (for discontinuous solutions) space-time discontinuous Galerkin (DG) method for systems of non-linear hyperbolic conservation laws. The resulting numerical scheme is fully discrete and provides a bound on the mathematical entropy at any time according to its initial condition and boundary conditions. The crux of the method is that discrete derivative approximations in space and time are summation-by-parts (SBP) operators. This allows the discrete method to mimic results from the continuous entropy analysis and ensures that the complete numerical scheme obeys the second law of thermodynamics. Importantly, the novel method described herein does not assume any exactness of quadrature in the variational forms that naturally arise in the context of DG methods. Typically, the development of entropy stable schemes is done on the semi-discrete level ignoring the temporal dependence. In this work we demonstrate that creating an entropy stable DG method in time is similar to the spatial discrete entropy analysis, but there are important (and subtle) differences. Therefore, we highlight the temporal entropy analysis throughout this work. For the compressible Euler equations, the preservation of kinetic energy is of interest besides entropy stability. The construction of kinetic energy preserving (KEP) schemes is, again, typically done on the semi-discrete level similar to the construction of entropy stable schemes. We present a generalization of the KEP condition from Jameson to the space-time framework and provide the temporal components for both entropy stability and kinetic energy preservation. The properties of the space-time DG method derived herein is validated through numerical tests for the compressible Euler equations.

Pei Fu, Gero Schnücke, Yinhua Xia

In Klingenberg, Schnücke and Xia (Math. Comp. 86 (2017), 1203-1232) an arbitrary Lagrangian-Eulerian discontinuous Galerkin (ALE-DG) method to solve conservation laws has been developed and analyzed. In this paper, the ALE-DG method will be extended to several dimensions. The method will be designed for simplex meshes. This will ensure that the method satisfies the geometric conservation law, if the accuracy of the time integrator is not less than the value of the spatial dimension. For the semi-discrete method the L2-stability will be proven. Furthermore, an error estimate which provides the suboptimal (k+1/2) convergence with respect to the L-infinity-norm will be presented, when an arbitrary monotone flux is used and for each cell the approximating functions are given by polynomials of degree $k$. The two dimensional fully-discrete explicit method will be combined with the bound preserving limiter developed by Zhang, Xia and Shu in (J. Sci. Comput. 50 (2012), 29-62). This limiter does not affect the high order accuracy of a numerical method. Then, for the ALE-DG method revised by the limiter the validity of a discrete maximum principle will be proven. The numerical stability, robustness and accuracy of the method will be shown by a variety of two dimensional computational experiments on moving triangular meshes.

Gero Schnücke, Nico Krais, Thomas Bolemann et al.

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.