On the order reduction of entropy stable DGSEM for the compressible Euler equations

On the order reduction of entropy stable DGSEM for the compressible Euler equations

Florian J. Hindenlang and Gregor J. Gassner Florian J. Hindenlang Max Planck Institute for Plasma Physics, Boltzmannstraße 2, D-85748 Garching, Germany, 22email: florian.hindenlang@ipp.mpg.deGregor J. GassnerDepartment for Mathematics and Computer Science; Center for Data and Simulation Science, University of Cologne, Cologne, Germany 44email: ggassner@math.uni-koeln.de
Abstract

Is the experimental order of convergence lower when using the entropy stable DGSEM-LGL variant? Recently, a debate on the question of the convergence behavior of the entropy stable nodal collocation discontinuous Galerkin spectral element method (DGSEM) with Legendre-Gauss-Lobatto nodes has emerged. Whereas it is well documented that the entropy conservative variant with no additional interface dissipation shows an odd-even behavior when testing its experimental convergence order, the results in the literature are less clear regarding the entropy stable version of the DGSEM-LGL, where explicit Riemann solver type dissipation is added at the element interfaces. We contribute to the ongoing discussion and present numerical experiments for the compressible Euler equations, where we investigate the effect of the choice of the numerical surface flux function. In our experiments, it turns out that the choice of the numerical surface flux has an impact on the convergence order. Penalty type numerical fluxes with high dissipation in all waves, such as the LLF and the HLL flux, appear to affect the convergence order negatively for odd polynomial degrees , in contrast to the entropy conserving variant, where even polynomial degrees are negatively affected. This behavior is more pronounced in low Mach number settings. In contrast, for numerical surface fluxes with less dissipative behavior in the contact wave such as e.g. Roe’s flux, the HLLC flux and the entropy conservative flux augmented with 5-wave matrix dissipation, optimal convergence rate of independent of the Mach number is observed.

Keywords:
entropy stability, discontinuous Galerkin spectral element method, summation-by-parts, order reduction, compressible Euler equations
\theoremstyle

plain \theoremstyleremark

1 Introduction

Discontinuous Galerkin spectral element collocation method (DGSEM) with either Legendre-Gauss or Legendre-Gauss-Lobatto (LGL) nodes (see e.g. [13]) are among the most efficient variants in the class of element based high order methods, such as e.g. discontinuous Galerkin, flux reconstruction, or summation-by-parts (SBP) finite differences. In particular, the LGL variant, starting in [8], turned out to be similar to a SBP finite difference approximation with simultaneous-approximate-term technique (SAT). This relationship allowed to construct conservative skew-symmetric approximations, e.g. [8, 10, 18], and later enabled DGSEM-LGL approximations that are discretely entropy stable, e.g. [2, 14, 15, 5, 11, 16, 1, 7], and/or kinetic energy preserving [12]. These novel variants of nodal split form DG methods feature drastically increased non-linear robustness towards aliasing induced instabilities and favourable properties regarding the simulation of unresolved turbulence, e.g. [6, 19].

In addition to the very robust dissipative entropy stable versions, it is also possible to construct virtually dissipation free variants by choosing appropriate element interface numerical fluxes. These entropy conserving variants all show an odd-even behavior when experimentally testing the order of convergence, e.g. [8, 18], where the observed convergence order for even polynomial degrees is and for odd is . Lately, a discussion emerged in the community, with interesting debates during the recent ICOSAHOM conference in London, where researchers reported non-optimal convergence behavior of the entropy stable DGSEM-LGL even with dissipative numerical surface fluxes, e.g. [5].

This paper contributes to this discussion and presents results of an experimental convergence order study for the compressible Euler equations with (i) the standard DGSEM with either Gauss and LGL nodes, (ii) the entropy stable DGSEM with LGL nodes. For these nodal schemes, we test the convergence order with different numerical surface fluxes and report the results depending on the Mach number of the test case. The remainder of the paper is organized as follows: in the next section we describe the numerical model for our numerical experiments, in Sec. 3 we present our observed experimental convergence orders for different configurations and draw our conclusion in Sec. 4.

2 Numerical Model

We consider the compressible Euler equations defined in the domain

(1)

The state vector contains the conservative variables and the advective flux components are

(2)

Here, are the mass density, fluid velocities, pressure and total energy. We close the system with the ideal gas assumption, which relates the total energy and pressure

(3)

where denotes the adiabatic coefficient.

For our discretization, we subdivide the domain into non-overlapping hexahedral elements. For each element, we define a transfinite mapping to a unit reference space and use this mapping to transform the equations (1) from physical to reference space. A weak form is created by taking the inner product of the transformed equation with a test function. We use integration-by-parts for the flux term and approximate the resulting weak form as follows: the conservative variables are approximated by a polynomial in reference space with degree N, interpolated at the Gauss or LGL nodes. The volume fluxes are replaced by a standard interpolation of the non-linear flux function at the same Gauss/LGL nodes (standard DGSEM-Gauss or DGSEM-LGL), see e.g. [13]. For the LGL variant, we are also able to introduce the split form volume integral based on entropy conserving and kinetic energy preserving numerical volume fluxes (Split-DGSEM), e.g. [12] and  [17], resulting in either the entropy conserving or entropy stable DGSEM variants, depending on the choice of numerical surface flux.

3 Convergence results

In this section, we compare the convergence of the standard DGSEM and the entropy conservative and entropy stable discretization for different choices of the numerical flux and polynomial degrees .

We choose the test case of a two-dimensional density wave, with a constant pressure and transported with a constant velocity, which was proposed for one-dimensional convergence tests in [4]. The density evolves as

(4)

with a prescribed velocity . The pressure is chosen as with , so that the sound speed ranges between . Thus, by changing the velocity, we change the Mach number of the flow . Three Mach numbers are chosen: with , with and with . The experimental order of convergence (EOC) is computed with the error of the density at .

The convergence study is performed with the open source, three-dimensional curvilinear split-form DG framework FLUXO (www.github.com/project-fluxo). As the test case is two-dimensional, we use fully periodic cartesian meshes of the domain with an equal number of elements in x- and y-directions and always element in z-direction. Note that in the convergence tables refers to the coarsest mesh level, which is elements for () and elements for ().

All simulation results are obtained with an explicit five stage, fourth order accurate low storage Runge-Kutta scheme [3], where a stable time step is computed according to the adjustable coefficient CFL, the local maximum wave speed, and the relative grid size, e.g. [9]. We made sure that the time integrator did not influence the spatial convergence order, by adjusting the CFL number accordingly.

Volume disc. Mach mesh
+Surface flux    level EOC EOC EOC EOC
DGSEM-Gauss + HLL 1.87e-04 3.34 8.57e-06 4.02 1.03e-05 5.02 6.76e-07 6.07
2.27e-05 3.04 5.35e-07 4.00 3.30e-07 4.96 1.07e-08 5.99
2.82e-06 3.00 3.34e-08 4.00 1.02e-08 5.01 1.66e-10 6.01
3.53e-07 3.00 2.09e-09 4.00 3.22e-10 4.99 2.60e-12 5.99
2.60e-04 2.55 5.92e-06 4.55 1.15e-05 4.39 6.74e-07 6.84
3.74e-05 2.80 3.27e-07 4.18 4.94e-07 4.54 7.23e-09 6.54
4.88e-06 2.94 1.95e-08 4.07 1.73e-08 4.84 9.98e-11 6.18
6.18e-07 2.98 1.20e-09 4.03 5.60e-10 4.95 1.52e-12 6.04
4.87e-04 1.70 4.36e-06 5.05 1.57e-05 2.56 9.80e-07 6.66
1.08e-04 2.18 1.10e-07 5.31 9.90e-07 3.99 4.47e-09 7.77
1.95e-05 2.46 5.84e-09 4.23 5.06e-08 4.29 3.96e-11 6.82
3.06e-06 2.67 2.34e-10 4.64 2.26e-09 4.49 4.57e-13 6.44
DGSEM-Gauss + Roe 1.87e-04 3.34 8.57e-06 4.02 1.03e-05 5.02 6.76e-07 6.07
2.27e-05 3.04 5.35e-07 4.00 3.30e-07 4.96 1.07e-08 5.99
2.82e-06 3.00 3.34e-08 4.00 1.02e-08 5.01 1.66e-10 6.01
3.53e-07 3.00 2.09e-09 4.00 3.22e-10 4.99 2.60e-12 5.99
1.82e-04 3.07 8.76e-06 4.04 1.11e-05 5.07 6.95e-07 6.08
2.26e-05 3.00 5.40e-07 4.02 3.44e-07 5.00 1.07e-08 6.01
2.82e-06 3.00 3.35e-08 4.01 1.05e-08 5.04 1.67e-10 6.01
3.53e-07 3.00 2.09e-09 4.00 3.25e-10 5.01 2.61e-12 6.00
2.14e-04 2.65 1.04e-05 3.78 1.16e-05 4.53 7.76e-07 5.78
2.22e-05 3.26 5.49e-07 4.25 3.49e-07 5.05 1.08e-08 6.17
2.82e-06 2.98 3.76e-08 3.87 1.04e-08 5.07 1.70e-10 5.99
3.53e-07 3.00 2.07e-09 4.19 3.38e-10 4.94 2.64e-12 6.01
Table 1: Experimental order of convergence of error to the exact density (4), using the standard DGSEM-Gauss with HLL and Roe fluxes. Full order is marked with () and an order reduction with .
Volume disc. Mach mesh
+Surface flux    level EOC EOC EOC EOC
DGSEM-LGL + HLL 1.22e-03 3.38 3.85e-05 4.03 4.34e-05 4.93 2.70e-06 5.93
1.26e-04 3.27 2.41e-06 4.00 1.39e-06 4.96 4.33e-08 5.96
1.48e-05 3.10 1.51e-07 4.00 4.34e-08 5.01 6.67e-10 6.02
1.81e-06 3.03 9.42e-09 4.00 1.36e-09 4.99 1.05e-11 5.98
1.04e-03 2.30 3.44e-05 4.44 4.11e-05 4.89 2.89e-06 6.33
1.58e-04 2.72 1.88e-06 4.20 1.77e-06 4.54 3.85e-08 6.23
2.11e-05 2.91 1.16e-07 4.02 6.25e-08 4.82 5.37e-10 6.16
2.69e-06 2.97 7.18e-09 4.01 2.05e-09 4.93 8.36e-12 6.00
1.20e-03 1.96 4.86e-05 3.51 5.16e-05 3.91 3.85e-06 5.31
2.72e-04 2.14 1.85e-06 4.71 2.99e-06 4.11 4.25e-08 6.50
5.57e-05 2.29 1.21e-07 3.94 1.56e-07 4.26 5.45e-10 6.28
1.01e-05 2.47 5.95e-09 4.34 7.19e-09 4.44 7.11e-12 6.26
DGSEM-LGL + Roe 1.22e-03 3.38 3.85e-05 4.03 4.34e-05 4.93 2.70e-06 5.93
1.26e-04 3.27 2.41e-06 4.00 1.39e-06 4.96 4.33e-08 5.96
1.48e-05 3.10 1.51e-07 4.00 4.34e-08 5.01 6.67e-10 6.02
1.81e-06 3.03 9.42e-09 4.00 1.36e-09 4.99 1.05e-11 5.98
9.17e-04 2.86 3.96e-05 3.94 4.41e-05 4.94 2.76e-06 6.02
1.15e-04 2.99 2.41e-06 4.04 1.47e-06 4.90 4.44e-08 5.96
1.44e-05 3.00 1.51e-07 4.00 4.38e-08 5.07 6.84e-10 6.02
1.80e-06 3.00 9.42e-09 4.00 1.35e-09 5.02 1.08e-11 5.99
9.26e-04 2.35 4.63e-05 3.45 4.26e-05 4.27 2.97e-06 5.37
1.19e-04 2.96 2.40e-06 4.27 1.59e-06 4.74 4.48e-08 6.05
1.43e-05 3.06 1.57e-07 3.93 4.34e-08 5.19 6.83e-10 6.04
1.80e-06 2.99 9.37e-09 4.07 1.49e-09 4.86 1.09e-11 5.97
Table 2: Experimental order of convergence of error to the exact density (4), using DGSEM-GL with HLL and Roe fluxes. Full order is marked with () and an order reduction with .

3.1 Standard DGSEM

The convergence of the standard DGSEM with Gauss-Legendre nodes (DGSEM-Gauss) and with Legendre-Gauss-Lobatto (DGSEM-LGL) is shown in Table 1 and Table 2, for the three Mach numbers and two choices of the numerical flux, namely the HLL and the Roe flux. The results of the local Lax-Friedrichs (LLF) and the HLLC fluxes are reported in the Appendix, as HLL is similar to LLF, and HLLC is the same as Roe, see Table 4 and Table 5.

For the HLL flux and the low Mach number , we observe an odd-even behavior with an order reduction for even polynomial degrees . Also for , the convergence for even degrees is slightly affected, whereas for the high Mach number, all fluxes converge with full order. Comparing the errors of the finest mesh for HLL and Roe for the low Mach number, HLL is less accurate for and more accurate for .

All numerical fluxes are approximate Riemann solvers, but the LLF and HLL only use the maximum wave speeds, whereas the HLLC and Roe also take the contact wave into account, and therefore keep the full order of the scheme for all Mach numbers for this test case.

Volume disc. Mach mesh
+Surface flux    level EOC EOC EOC EOC
Split-DGSEM + ECKEP 1.62e-03 4.23 8.45e-05 2.77 5.80e-05 5.88 4.63e-06 5.63
1.30e-04 3.64 7.14e-06 3.56 1.60e-06 5.18 8.38e-08 5.79
1.05e-05 3.62 7.90e-07 3.18 4.56e-08 5.13 3.59e-09 4.54
1.69e-06 2.64 9.58e-08 3.04 1.23e-09 5.21 1.18e-10 4.93
1.41e-03 3.89 9.45e-05 2.34 7.71e-05 5.86 3.49e-06 6.00
1.25e-04 3.49 1.26e-05 2.90 1.98e-06 5.28 5.00e-08 6.12
1.48e-05 3.09 1.60e-06 2.98 4.06e-08 5.61 1.22e-09 5.36
1.32e-06 3.48 2.01e-07 3.00 1.13e-09 5.17 3.85e-11 4.98
1.13e-03 2.26 8.03e-05 2.98 5.95e-05 4.14 4.11e-06 5.14
1.13e-04 3.32 1.02e-05 2.97 1.85e-06 5.01 1.33e-07 4.95
1.12e-05 3.34 1.29e-06 2.99 4.02e-08 5.52 4.21e-09 4.98
1.98e-06 2.50 1.61e-07 3.00 1.43e-09 4.81 1.32e-10 5.00
Split-DGSEM + HLL 1.23e-03 3.36 3.88e-05 4.06 4.49e-05 4.94 3.01e-06 5.97
1.27e-04 3.27 2.42e-06 4.00 1.43e-06 4.97 4.84e-08 5.96
1.48e-05 3.10 1.51e-07 4.00 4.44e-08 5.01 7.48e-10 6.02
1.81e-06 3.03 9.46e-09 4.00 1.40e-09 4.99 1.18e-11 5.98
1.04e-03 2.30 3.48e-05 4.44 4.16e-05 4.90 3.48e-06 6.10
1.58e-04 2.72 1.88e-06 4.21 1.79e-06 4.54 4.66e-08 6.23
2.11e-05 2.90 1.16e-07 4.02 6.39e-08 4.81 6.08e-10 6.26
2.69e-06 2.97 7.20e-09 4.01 2.11e-09 4.92 9.28e-12 6.03
1.20e-03 1.96 4.90e-05 3.51 5.23e-05 3.92 4.47e-06 5.16
2.72e-04 2.14 1.87e-06 4.71 3.02e-06 4.12 6.14e-08 6.19
5.58e-05 2.29 1.21e-07 3.95 1.57e-07 4.26 6.18e-10 6.64
1.01e-05 2.47 5.98e-09 4.34 7.28e-09 4.43 8.20e-12 6.23
Split-DGSEM + ECKEP-Roe 1.23e-03 3.36 3.88e-05 4.06 4.49e-05 4.94 3.01e-06 5.97
1.27e-04 3.27 2.42e-06 4.00 1.43e-06 4.97 4.84e-08 5.96
1.48e-05 3.10 1.51e-07 4.00 4.44e-08 5.01 7.48e-10 6.02
1.81e-06 3.03 9.46e-09 4.00 1.40e-09 4.99 1.18e-11 5.98
9.18e-04 2.86 3.98e-05 3.94 4.51e-05 4.92 3.08e-06 5.95
1.15e-04 2.99 2.42e-06 4.04 1.51e-06 4.90 4.97e-08 5.95
1.44e-05 3.00 1.51e-07 4.00 4.49e-08 5.07 7.66e-10 6.02
1.80e-06 3.00 9.46e-09 4.00 1.38e-09 5.02 1.21e-11 5.99
9.26e-04 2.35 4.65e-05 3.45 4.37e-05 4.26 3.26e-06 5.34
1.19e-04 2.96 2.40e-06 4.27 1.63e-06 4.74 5.03e-08 6.02
1.43e-05 3.06 1.58e-07 3.93 4.45e-08 5.19 7.66e-10 6.04
1.80e-06 2.99 9.40e-09 4.07 1.53e-09 4.86 1.22e-11 5.97
Table 3: Experimental order of convergence of error to the exact density (4), using entropy conservative ECKEP flux and entropy stable HLL and ECKEP-Roe fluxes. Full order is marked with () and an order reduction with .

3.2 Entropy conservative and entropy stable DGSEM

Now, we investigate the order reduction of the entropy conservative and entropy stable discretizations. Here, the standard DGSEM volume integral is replaced by split-form formulation (Split-DGSEM) using a two-point entropy conservative and kinetic energy preserving flux ECKEP. If we choose the ECKEP flux at the surface, we get an entropy-conserving scheme. For entropy stability, we can use the LLF or HLL flux directly at the surface, or use the ECKEP flux and add a dissipation term, which must still satisfy the entropy inequality condition. In Winters et al. [17], such dissipation terms are carefully derived, using either only the maximum wave speed (LLF-type) or incorporating all waves (Roe-type), which we will refer to as ECKEP-LLF and ECKEP-Roe fluxes.

In Table 3, we summarize the convergence of the dissipation-free ECKEP flux, the HLL and ECKEP-Roe flux. The results for LLF and ECKEP-LF fluxes are found in the Appendix in Table 6, as they have the same convergence and error levels as the HLL flux. As expected, the dissipation-free surface flux (ECKEP) produces an order reduction for all Mach numbers for , and for full order not kept in the last refinement step.

If we simply use the HLL flux, we have an entropy stable scheme, but an order reduction for can be observed for the low Mach number flow, analogously to the standard DGSEM-LGL scheme. Interestingly, the odd-even behavior switches between entropy conserving and entropy stable fluxes.

The ECKEP-Roe entropy stable flux accounts for all waves of the Riemann problem and adjusts the dissipation for each wave accordingly, which gives full order convergence for all Mach numbers.

4 Conclusions

In this work, we report the convergence of standard DGSEM Gauss and Gauss-Lobatto schemes to entropy conservative (EC) and entropy stable (ES) DGSEM schemes for the Euler equations, as there have been findings of order reduction for EC and ES schemes. We choose a simple density transport test case on a periodic domain and investigate the influence of the Mach number of the transport velocity.

The EC scheme is dissipation free and an order reduction is observed by the convergence study presented here, confirming many similar observations found in literature. We also confirm that the ES scheme can have an order reduction for low Mach numbers, but only if the entropy stable numerical flux relies on simple approximate Riemann solvers such as local Lax-Friedrichs or HLL. If all waves are accounted for in the dissipation term of the entropy stable flux as presented in [17], the full order is observed for all Mach numbers. IN addition, we reproduce the same behavior for the standard DGSEM Gauss and Gauss-Lobatto schemes, where the LLF and HLL fluxes suffer from order reduction at low Mach number, and HLLC and Roe fluxes have full order for all Mach numbers.

We want to emphasize that the present convergence study should be seen merely as an observation, confirming that the numerical flux can have strong influence on the convergence order for both the standard DGSEM and the entropy stable DGSEM. Also, we stress that in our tests the order reduction is related to the form of the dissipation term in the numerical surface flux and is not related to the insufficient integration precision of the LGL-quadrature.

In the authors’ experience, a convergence study using a manufactured solution technique can be misleading, as full convergence order is found independent of the choice of numerical flux. Hence, the introduction of a source term to balance the prescribed solution overcomes possible deficiencies of the surface fluxes, showing the limit of the manufactured solution technique in this context. In the Appendix, the convergence results of a manufactured solution are reported.

Acknowledgements

Gregor Gassner thanks the European Research Council for funding through the ERC Starting Grant “An Exascale aware and Un-crashable Space-Time-Adaptive Discontinuous Spectral Element Solver for Non-Linear Conservation Laws” (Extreme), ERC grant agreement no. 714487. Florian Hindenlang thanks Eric Sonnendrücker and the Max-Planck Institute for Plasma Physics in Garching for their constant support. We would also like all participants of the ICOSAHOM 2018 for the valuable discussions on the topic of entropy stable schemes, which motivated this work.

Bibliography

  • [1] M. Bohm, A. R. Winters, G. J. Gassner, D. Derigs, F. Hindenlang, and J. Saur. An entropy stable nodal discontinuous Galerkin method for the resistive MHD equations. part I: Theory and numerical verification. Journal of Computational Physics, 2018.
  • [2] M. Carpenter, T. Fisher, E. Nielsen, and S. Frankel. Entropy stable spectral collocation schemes for the Navier–Stokes equations: Discontinuous interfaces. SIAM Journal on Scientific Computing, 36(5):B835–B867, 2014.
  • [3] M. Carpenter and C. Kennedy. Fourth-order 2N-storage Runge-Kutta schemes. Technical Report NASA TM 109111, 1994.
  • [4] J. Chan. On discretely entropy conservative and entropy stable discontinuous Galerkin methods. Journal of Computational Physics, 362:346 – 374, 2018.
  • [5] T. Chen and C.-W. Shu. Entropy stable high order discontinuous Galerkin methods with suitable quadrature rules for hyperbolic conservation laws. Journal of Computational Physics, 345:427 – 461, 2017.
  • [6] D. Flad and G. Gassner. On the use of kinetic energy preserving DG-schemes for large eddy simulation. Journal of Computational Physics, 350:782 – 795, 2017.
  • [7] L. Friedrich, A. R. Winters, D. C. Del Rey Fernández, G. J. Gassner, M. Parsani, and M. H. Carpenter. An entropy stable h / p non-conforming discontinuous Galerkin method with the summation-by-parts property. Journal of Scientific Computing, May 2018.
  • [8] G. Gassner. A skew-symmetric discontinuous Galerkin spectral element discretization and its relation to SBP-SAT finite difference methods. SIAM Journal on Scientific Computing, 35(3):A1233–A1253, 2013.
  • [9] G. Gassner and D. A. Kopriva. A comparison of the dispersion and dissipation errors of Gauss and Gauss–Lobatto discontinuous Galerkin spectral element methods. SIAM Journal on Scientific Computing, 33:2560–2579, 2011.
  • [10] G. J. Gassner. A kinetic energy preserving nodal discontinuous Galerkin spectral element method. International Journal for Numerical Methods in Fluids, 76(1):28–50, 2014.
  • [11] G. J. Gassner, A. R. Winters, F. J. Hindenlang, and D. A. Kopriva. The BR1 scheme is stable for the compressible Navier–Stokes equations. Journal of Scientific Computing, Apr 2018.
  • [12] G. J. Gassner, A. R. Winters, and D. A. Kopriva. Split form nodal discontinuous Galerkin schemes with summation-by-parts property for the compressible Euler equations. Journal Of Computational Physics, 327:39–66, 2016.
  • [13] D. A. Kopriva. Implementing Spectral Methods for Partial Differential Equations. Scientific Computation. Springer, May 2009.
  • [14] M. Parsani, M. H. Carpenter, and E. J. Nielsen. Entropy stable wall boundary conditions for the three-dimensional compressible Navier–Stokes equations. Journal of Computational Physics, 292:88 – 113, 2015.
  • [15] N. Wintermeyer, A. R. Winters, G. J. Gassner, and D. A. Kopriva. An entropy stable nodal discontinuous Galerkin method for the two dimensional shallow water equations on unstructured curvilinear meshes with discontinuous bathymetry. Journal of Computational Physics, submitted. arXiv:1509.07096 [math.NA], 2016.
  • [16] N. Wintermeyer, A. R. Winters, G. J. Gassner, and T. Warburton. An entropy stable discontinuous Galerkin method for the shallow water equations on curvilinear meshes with wet/dry fronts accelerated by GPUs. Journal of Computational Physics, 375:447 – 480, 2018.
  • [17] A. R. Winters, D. Derigs, G. J. Gassner, and S. Walch. A uniquely defined entropy stable matrix dissipation operator for high Mach number ideal MHD and compressible Euler simulations. Journal of Computational Physics, 332:274 – 289, 2017.
  • [18] A. R. Winters and G. J. Gassner. A comparison of two entropy stable discontinuous Galerkin spectral element approximations for the shallow water equations with non-constant topography. Journal Of Computational Physics, 301:357–376, 2015.
  • [19] A. R. Winters, R. C. Moura, G. Mengaldo, G. J. Gassner, S. Walch, J. Peiro, and S. J. Sherwin. A comparative study on polynomial dealiasing and split form discontinuous Galerkin schemes for under-resolved turbulence computations. Journal of Computational Physics, 372:1 – 21, 2018.

Appendix

Additional convergence results

In this section, we present additional convergence results of the density wave tst case for the DGSEM-Gauss and DGSEM-LGL with LLF and HLLC fluxes in Table 4 and Table 5, and also the entropy stable schemes with LLF and ECKEP-LLF fluxes in Table 6. The results for LLF-type fluxes behave like the HLL flux, and for the HLLC flux like the Roe-type fluxes presented in Table 3.

Volume disc. Mach mesh
+Surface flux    level EOC EOC EOC EOC
DGSEM-Gauss + LLF 2.42e-04 2.97 6.43e-06 4.41 1.05e-05 4.51 6.68e-07 6.60
3.24e-05 2.90 3.71e-07 4.11 4.32e-07 4.60 8.04e-09 6.38
4.15e-06 2.96 2.27e-08 4.03 1.47e-08 4.87 1.15e-10 6.13
5.22e-07 2.99 1.41e-09 4.01 4.73e-10 4.96 1.77e-12 6.02
3.13e-04 2.29 4.59e-06 4.84 1.18e-05 3.90 6.69e-07 7.40
5.30e-05 2.56 2.25e-07 4.35 6.08e-07 4.28 5.56e-09 6.91
7.43e-06 2.83 1.29e-08 4.12 2.47e-08 4.62 6.79e-11 6.35
9.61e-07 2.95 7.65e-10 4.07 8.53e-10 4.85 9.97e-13 6.09
4.95e-04 1.69 4.33e-06 5.07 1.58e-05 2.53 9.88e-07 6.68
1.12e-04 2.15 1.06e-07 5.35 1.01e-06 3.97 4.46e-09 7.79
2.06e-05 2.44 5.47e-09 4.28 5.23e-08 4.27 3.86e-11 6.85
3.29e-06 2.65 2.15e-10 4.67 2.38e-09 4.46 4.35e-13 6.47
DGSEM-Gauss + HLLC 1.87e-04 3.34 8.57e-06 4.02 1.03e-05 5.02 6.76e-07 6.07
2.27e-05 3.04 5.35e-07 4.00 3.30e-07 4.96 1.07e-08 5.99
2.82e-06 3.00 3.34e-08 4.00 1.02e-08 5.01 1.66e-10 6.01
3.53e-07 3.00 2.09e-09 4.00 3.22e-10 4.99 2.60e-12 5.99
1.82e-04 3.07 8.76e-06 4.04 1.11e-05 5.07 6.95e-07 6.08
2.26e-05 3.00 5.40e-07 4.02 3.44e-07 5.00 1.07e-08 6.01
2.82e-06 3.00 3.35e-08 4.01 1.05e-08 5.04 1.67e-10 6.01
3.53e-07 3.00 2.09e-09 4.00 3.25e-10 5.01 2.61e-12 6.00
2.14e-04 2.65 1.04e-05 3.78 1.16e-05 4.53 7.76e-07 5.78
2.22e-05 3.26 5.49e-07 4.25 3.49e-07 5.05 1.08e-08 6.17
2.82e-06 2.98 3.76e-08 3.87 1.04e-08 5.07 1.70e-10 5.99
3.53e-07 3.00 2.07e-09 4.19 3.38e-10 4.94 2.64e-12 6.01
Table 4: Experimental order of convergence of error to the exact density (4), using DGSEM-Gauss with LLF and HLLC fluxes. Full order is marked with () and an order reduction with .
Volume disc. Mach mesh
+Surface flux    level EOC EOC EOC EOC
DGSEM-LGL + LLF 1.23e-03 3.11 3.44e-05 4.36 4.02e-05 4.82 2.91e-06 6.14
1.51e-04 3.02 2.00e-06 4.10 1.61e-06 4.64 3.80e-08 6.26
1.89e-05 3.00 1.23e-07 4.02 5.55e-08 4.86 5.46e-10 6.12
2.36e-06 3.00 7.66e-09 4.01 1.79e-09 4.95 8.54e-12 6.00
1.10e-03 1.97 3.15e-05 4.83 4.02e-05 4.93 2.85e-06 6.73
2.03e-04 2.43 1.68e-06 4.23 1.98e-06 4.34 3.68e-08 6.28
2.99e-05 2.76 1.03e-07 4.02 8.35e-08 4.57 4.86e-10 6.24
3.94e-06 2.93 6.39e-09 4.01 2.96e-09 4.82 7.53e-12 6.01
1.21e-03 1.95 4.88e-05 3.51 5.20e-05 3.90 3.89e-06 5.31
2.77e-04 2.12 1.86e-06 4.72 3.05e-06 4.09 4.27e-08 6.51
5.78e-05 2.26 1.21e-07 3.94 1.61e-07 4.24 5.47e-10 6.29
1.06e-05 2.44 5.98e-09 4.34 7.57e-09 4.41 7.16e-12 6.26
DGSEM-LGL + HLLC 1.22e-03 3.38 3.85e-05 4.03 4.34e-05 4.93 2.70e-06 5.93
1.26e-04 3.27 2.41e-06 4.00 1.39e-06 4.96 4.33e-08 5.96
1.48e-05 3.10 1.51e-07 4.00 4.34e-08 5.01 6.67e-10 6.02
1.81e-06 3.03 9.42e-09 4.00 1.36e-09 4.99 1.05e-11 5.98
9.17e-04 2.86 3.96e-05 3.94 4.41e-05 4.94 2.76e-06 6.02
1.15e-04 2.99 2.41e-06 4.04 1.47e-06 4.90 4.44e-08 5.96
1.44e-05 3.00 1.51e-07 4.00 4.38e-08 5.07 6.84e-10 6.02
1.80e-06 3.00 9.42e-09 4.00 1.35e-09 5.02 1.08e-11 5.99
9.26e-04 2.35 4.63e-05 3.45 4.26e-05 4.27 2.97e-06 5.37
1.19e-04 2.96 2.40e-06 4.27 1.59e-06 4.74 4.48e-08 6.05
1.43e-05 3.06 1.57e-07 3.93 4.34e-08 5.19 6.83e-10 6.04
1.80e-06 2.99 9.37e-09 4.07 1.49e-09 4.86 1.09e-11 5.97
Table 5: Experimental order of convergence of error to the exact density (4), using DGSEM-GL with LLF and HLLC fluxes. Full order is marked with () and an order reduction with .
Volume disc. Mach mesh
+Surface flux    level EOC EOC EOC EOC
Split-DGSEM + LLF 1.24e-03 3.11 3.49e-05 4.37 4.16e-05 4.85 3.37e-06 6.05
1.52e-04 3.03 2.01e-06 4.12 1.64e-06 4.67 4.49e-08 6.23
1.89e-05 3.00 1.23e-07 4.03 5.67e-08 4.85 6.21e-10 6.18
2.36e-06 3.00 7.68e-09 4.01 1.84e-09 4.95 9.58e-12 6.02
1.09e-03 1.97 3.21e-05 4.82 4.07e-05 4.94 3.77e-06 6.34
2.03e-04 2.43 1.69e-06 4.25 1.99e-06 4.35 4.69e-08 6.33
2.99e-05 2.76 1.04e-07 4.03 8.46e-08 4.56 5.57e-10 6.39
3.94e-06 2.93 6.41e-09 4.02 3.03e-09 4.80 8.38e-12 6.05
1.21e-03 1.95 4.92e-05 3.50 5.26e-05 3.91 4.52e-06 5.15
2.77e-04 2.12 1.88e-06 4.71 3.07e-06 4.10 6.21e-08 6.19
5.79e-05 2.26 1.22e-07 3.95 1.63e-07 4.24 6.21e-10 6.64
1.06e-05 2.44 6.01e-09 4.34 7.66e-09 4.41 8.26e-12 6.23
Split-DGSEM + ECKEP-LLF 1.24e-03 3.10 3.49e-05 4.37 4.16e-05 4.84 3.38e-06 6.06
1.52e-04 3.03 2.01e-06 4.12 1.64e-06 4.67 4.49e-08 6.23
1.89e-05 3.01 1.23e-07 4.03 5.68e-08 4.85 6.21e-10 6.18
2.36e-06 3.00 7.68e-09 4.01 1.84e-09 4.95 9.58e-12 6.02
1.10e-03 1.97 3.20e-05 4.83 4.08e-05 4.94 3.77e-06 6.36
2.04e-04 2.43 1.69e-06 4.25 2.00e-06 4.35 4.69e-08 6.33
3.00e-05 2.76 1.04e-07 4.03 8.49e-08 4.56 5.57e-10 6.40
3.95e-06 2.93 6.41e-09 4.02 3.04e-09 4.81 8.38e-12 6.05
1.21e-03 1.95 4.93e-05 3.51 5.27e-05 3.91 4.53e-06 5.15
2.78e-04 2.12 1.88e-06 4.71 3.08e-06 4.10 6.21e-08 6.19
5.79e-05 2.26 1.22e-07 3.95 1.63e-07 4.24 6.21e-10 6.65
1.07e-05 2.44 6.01e-09 4.34 7.67e-09 4.41 8.26e-12 6.23
Table 6: Experimental order of convergence of error to the exact density (4), using entropy stable LLF and ECKEP-LLF flux. Full order is marked with () and an order reduction with .

Manufactured solution with source term

Here, we run a convergence test with the method of manufactured solutions. To do so, we assume a two-dimensional solution of the form

(5)

Note that the average Mach number in the domain is . Inserting (5) into the Euler equations, and using the fact that spatial and time derivatives are , we get an additional residual

(6)

To solve the inhomogeneous problem, we subtract the residual from the approximate solution in each Runge-Kutta step. Moreover, we run the test case up to the final time .

In the convergence results for the standard DGSEM Gauss and Gauss-Lobatto, we see that the LLF flux still leads to an order reduction for , whereas full order is found for the HLL, HLLC and Roe fluxes, see Table 7 and Table 8.

In Table 9 the entropy conservative scheme shows again an order reduction for , and the LLF-Type dissipation too, for , and for this test case, all entropy stable schemes exhibit full order.

Volume disc. mesh
+Surface flux    level EOC EOC EOC EOC
DGSEM-Gauss + LLF 2.30e-03 2.20 4.54e-05 5.34 1.13e-04 6.35 4.52e-05 4.38
4.81e-04 2.25 1.99e-06 4.52 4.78e-06 4.56 2.37e-07 7.58
9.48e-05 2.34 1.02e-07 4.28 2.88e-07 4.05 1.53e-09 7.28
1.57e-05 2.60 6.25e-09 4.03 1.73e-08 4.05 1.57e-11 6.61
DGSEM-Gauss + HLL 1.24e-03 2.84 5.46e-05 4.35 1.32e-04 5.22 1.47e-05 6.30
1.17e-04 3.41 3.36e-06 4.02 2.89e-06 5.51 1.44e-07 6.67
1.41e-05 3.04 1.85e-07 4.18 7.14e-08 5.34 1.74e-09 6.37
1.76e-06 3.00 1.07e-08 4.11 2.15e-09 5.05 2.24e-11 6.28
DGSEM-Gauss + HLLC 1.24e-03 2.84 5.46e-05 4.35 1.32e-04 5.22 1.47e-05 6.30
1.17e-04 3.41 3.36e-06 4.02 2.89e-06 5.51 1.44e-07 6.67
1.41e-05 3.04 1.85e-07 4.18 7.14e-08 5.34 1.74e-09 6.37
1.76e-06 3.00 1.07e-08 4.11 2.15e-09 5.05 2.24e-11 6.28
DGSEM-Gauss + Roe 1.24e-03 2.84 5.46e-05 4.35 1.32e-04 5.22 1.47e-05 6.30
1.17e-04 3.41 3.36e-06 4.02 2.89e-06 5.51 1.44e-07 6.67
1.41e-05 3.04 1.85e-07 4.18 7.14e-08 5.34 1.74e-09 6.37
1.76e-06 3.00 1.07e-08 4.11 2.15e-09 5.05 2.24e-11 6.28
Table 7: Experimental order of convergence of error of density for the manufactured solution (5), using DGSEM-Gauss with LLF, HLL, HLLC and Roe fluxes. Full order is marked with () and an order reduction with .
Volume disc. mesh
+Surface flux    level EOC EOC EOC EOC
DGSEM-LGL + LLF 7.36e-03 2.85 3.15e-04 4.38 5.69e-04 5.78 9.33e-05 5.34
1.33e-03 2.47 1.43e-05 4.46 2.04e-05 4.80 9.27e-07 6.65
2.79e-04 2.25 7.99e-07 4.16 8.81e-07 4.54 9.34e-09 6.63
5.31e-05 2.39 4.72e-08 4.08 5.94e-08 3.89 1.37e-10 6.09
DGSEM-LGL + HLL 5.32e-03 3.16 2.52e-04 4.00 3.84e-04 5.36 4.33e-05 6.24
5.99e-04 3.15 1.38e-05 4.19 1.43e-05 4.75 4.58e-07 6.56
7.25e-05 3.05 7.69e-07 4.17 2.92e-07 5.61 7.08e-09 6.02
9.02e-06 3.01 4.74e-08 4.02 7.77e-09 5.23 1.10e-10 6.01
DGSEM-LGL + HLLC 5.32e-03 3.16 2.52e-04 4.00 3.84e-04 5.36 4.33e-05 6.24
5.99e-04 3.15 1.38e-05 4.19 1.43e-05 4.75 4.58e-07 6.56
7.25e-05 3.05 7.69e-07 4.17 2.92e-07 5.61 7.08e-09 6.02
9.02e-06 3.01 4.74e-08 4.02 7.77e-09 5.23 1.10e-10 6.01
DGSEM-LGL + Roe 5.32e-03 3.16 2.52e-04 4.00 3.84e-04 5.36 4.33e-05 6.24
5.99e-04 3.15 1.38e-05 4.19 1.43e-05 4.75 4.58e-07 6.56
7.25e-05 3.05 7.69e-07 4.17 2.92e-07 5.61 7.08e-09 6.02
9.02e-06 3.01 4.74e-08 4.02 7.77e-09 5.23 1.10e-10 6.01
Table 8: Experimental order of convergence of error of density for the manufactured solution (5), using DGSEM-LGL with LLF, HLL, HLLC and Roe fluxes. Full order is marked with () and an order reduction with .
Volume disc. mesh
+Surface flux    level EOC EOC EOC EOC
Split-DGSEM + ECKEP 1.31e-02 3.53 1.28e-03 2.88 5.62e-03 2.49 1.12e-03 4.05
1.30e-03 3.34 1.13e-04 3.49 3.12e-04 4.17 5.97e-06 7.56
1.23e-04 3.40 1.24e-05 3.20 2.00e-06 7.29 9.45e-08 5.98
1.76e-05 2.80 1.67e-06 2.89 3.41e-08 5.87 3.17e-09 4.90
Split-DGSEM + LLF 7.60e-03 3.13 3.70e-04 3.91 6.56e-04 4.44 1.28e-04 4.75
1.63e-03 2.22 1.90e-05 4.29 3.09e-05 4.41 2.01e-06 6.00
3.40e-04 2.26 9.89e-07 4.26 1.83e-06 4.08 1.80e-08 6.80
6.14e-05 2.47 6.41e-08 3.95 9.15e-08 4.32 2.17e-10 6.37
Split-DGSEM + ECKEP-LLF 7.66e-03 3.21 3.74e-04 3.98 6.79e-04 4.45 1.33e-04 4.75
1.64e-03 2.23 1.90e-05 4.30 3.08e-05 4.46 2.06e-06 6.01
3.41e-04 2.26 9.88e-07 4.26 1.83e-06 4.07 1.79e-08 6.85
6.16e-05 2.47 6.40e-08 3.95 9.17e-08 4.32 2.17e-10 6.37
Split-DGSEM + HLL 5.82e-03 3.10 3.01e-04 3.85 5.11e-04 4.35 7.06e-05 5.27
7.06e-04 3.04 2.04e-05 3.88 1.67e-05 4.94 1.08e-06 6.03
8.63e-05 3.03 1.16e-06 4.14 5.08e-07 5.04 1.67e-08 6.02
1.08e-05 3.00 7.20e-08 4.01 1.62e-08 4.97 2.64e-10 5.98
Split-DGSEM + ECKEP-Roe 5.81e-03 3.11 3.01e-04 3.85 5.11e-04 4.35 7.06e-05 5.27
7.06e-04 3.04 2.04e-05 3.88 1.67e-05 4.94 1.08e-06 6.03
8.63e-05 3.03 1.16e-06 4.14 5.08e-07 5.04 1.67e-08 6.02
1.08e-05 3.00 7.20e-08 4.01 1.62e-08 4.97 2.64e-10 5.98
Table 9: Experimental order of convergence of error of density for the manufactured solution (5), using entropy conservative and entropy stable schemes. Full order is marked with () and an order reduction with .
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
""
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
   
Add comment
Cancel
Loading ...
332289
This is a comment super asjknd jkasnjk adsnkj
Upvote
Downvote
""
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Submit
Cancel

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test
Test description