A Staggered Semi-implicit Discontinuous Galerkin Scheme with a Posteriori Subcell Finite Volume Limiter for the Euler Equations of Gasdynamics

Abstract

In this paper we propose a novel semi-implicit Discontinuous Galerkin (DG) finite element scheme on staggered meshes with a posteriori subcell finite volume limiting for the one and two dimensional Euler equations of compressible gasdynamics. We therefore extend the strategy adopted by Dumbser and Casulli (Appl Math Comput 272:479–497, 2016), where the Euler equations have been solved solved using a semi-implicit finite volume scheme based on the flux-vector splitting method recently proposed by Toro and Vazquez-Cendon (Comput Fluids 70:1–12, 2012). In our scheme, the nonlinear convective terms are discretized explicitly, while the pressure terms are discretized implicitly. As a consequence, the time step is restricted only by a mild CFL condition based on the fluid velocity, which makes this method particularly suitable for simulations in the low Mach number regime. However, the conservative formulation of the scheme, together with the novel subcell finite volume limiter allows also the numerical simulation of high Mach number flows with shock waves. Inserting the discrete momentum equation into the discrete total energy conservation law yields a mildly nonlinear system with the scalar pressure as the only unknown. Due to the use of staggered meshes, the resulting pressure system has the most compact stencil possible and can be efficiently solved with modern iterative methods. In order to deal with shock waves or steep gradients, the new semi-implicit DG scheme proposed in this paper includes an a posteriori subcell finite volume limiting technique. This strategy was first proposed by Dumbser et al. (J Comput Phys 278:47–75, 2014) for explicit DG schemes on collocated grids and is based on the a posteriori MOOD algorithm of Clain, Loubere and Diot. Recently, this methodology was also extended to semi-implicit DG schemes on staggered meshes for the shallow water equations in Ioriatti and Dumbser (Appl Numer Math 135:443–480, 2019). Within the MOOD approach, an unlimited DG scheme first produces a so-called candidate solution for the next time level $$t^{n+1}$$. Later on, the control volumes with a non-admissible candidate solution are identified by using physical and numerical detection criteria, such as the positivity of the solution, the absence of floating point errors and the satisfaction of a relaxed discrete maximum principle (DMP). Then, in the detected troubled cells a more robust first order semi-implicit finite volume (FV) method is applied on a sub-grid composed of $$2P + 1$$ subcells, where P denotes the polynomial degree used in the DG scheme. For that purpose, the nonlinear convective terms are recomputed in the troubled cells using an explicit finite volume scheme on the subcell level. Also the linear system for the pressure needs to be assembled and solved again, but where now a low order semi-implicit finite volume scheme is used on the sub-cell level in all troubled DG elements, instead of the original high order DG method. Finally, the higher order DG polynomials are reconstructed from the piecewise constant subcell finite volume averages and the scheme proceeds to the next time step. In this paper we present, discuss and test this novel family of methods and simulate a set of classical numerical benchmark problems of compressible gasdynamics. Great attention is dedicated to 1D and 2D Riemann problems and we also show that for these test cases the scheme responds appropriately in the presence of shock waves and does not produce non-physical spurious numerical oscillations.