Communications on Applied Mathematics and Computation ›› 2024, Vol. 6 ›› Issue (3): 1860-1898.doi: 10.1007/s42967-023-00321-6

• ORIGINAL PAPERS • Previous Articles     Next Articles

Monolithic Convex Limiting for Legendre-Gauss-Lobatto Discontinuous Galerkin Spectral-Element Methods

Andrés M. Rueda-Ramírez1, Benjamin Bolm1, Dmitri Kuzmin2, Gregor J. Gassner1,3   

  1. 1 Department of Mathematics and Computer Science, University of Cologne, Weyertal 86-90, Cologne 50931, Germany;
    2 Institute of Applied Mathematics (LS III), TU Dortmund University, Vogelpothsweg 87, Dortmund 44227, Germany;
    3 Center for Data and Simulation Science, University of Cologne, Weyertal 86-90, Cologne 50931, Germany
  • Received:2023-05-01 Revised:2023-11-15 Accepted:2023-09-17 Published:2024-12-20
  • Contact: Andrés M. Rueda-Ramírez,aruedara@uni-koeln.de;Benjamin Bolm,bbolm@smail.uni-koeln.de;Dmitri Kuzmin,kuzmin@math.uni-dortmund.de;Gregor J. Gassner,ggassner@uni-koeln.de E-mail:aruedara@uni-koeln.de;bbolm@smail.uni-koeln.de;kuzmin@math.uni-dortmund.de;ggassner@uni-koeln.de
  • Supported by:
    The authors thank Dr. Hennes Hajduk for performing comparative studies with his code and giving a deeper insight into the subcell DG-MCL schemes he developed in [22] for high-order Bernstein finite elements. Moreover, the authors would like to thank Prof. Dr. Hendrik Ranocha for his support and advice during the implementation of FCT methods in Trixi.jl [57, 62, 63]. Gregor Gassner and Andrés M. RuedaRamírez acknowledge funding through the Klaus Tschira Stiftung via the project “HiFiLab”. Gregor Gassner further acknowledges funding by the German Research Foundation (DFG) under the grant number DFGFOR5409. Dmitri Kuzmin acknowledges DFG support under the grant number KU 1530/23-3. We furthermore thank the Regional Computing Center of the University of Cologne (RRZK) for providing computing time on the High Performance Computing (HPC) system ODIN, as well as for technical support.

Abstract: We extend the monolithic convex limiting (MCL) methodology to nodal discontinuous Galerkin spectral-element methods (DGSEMS). The use of Legendre-Gauss-Lobatto (LGL) quadrature endows collocated DGSEM space discretizations of nonlinear hyperbolic problems with properties that greatly simplify the design of invariant domain-preserving high-resolution schemes. Compared to many other continuous and discontinuous Galerkin method variants, a particular advantage of the LGL spectral operator is the availability of a natural decomposition into a compatible subcell flux discretization. Representing a highorder spatial semi-discretization in terms of intermediate states, we perform flux limiting in a manner that keeps these states and the results of Runge-Kutta stages in convex invariant domains. In addition, local bounds may be imposed on scalar quantities of interest. In contrast to limiting approaches based on predictor-corrector algorithms, our MCL procedure for LGL-DGSEM yields nonlinear flux approximations that are independent of the time-step size and can be further modified to enforce entropy stability. To demonstrate the robustness of MCL/DGSEM schemes for the compressible Euler equations, we run simulations for challenging setups featuring strong shocks, steep density gradients, and vortex dominated flows.

Key words: Structure-preserving schemes, Subcell flux limiting, Monolithic convex limiting (MCL), Discontinuous Galerkin spectral-element methods (DGSEMS), Legendre-Gauss-Lobatto (LGL) nodes

CLC Number: