Abstract
Many modern discontinuous Galerkin (DG) methods for conservation laws make use of summation by parts operators and flux differencing to achieve kinetic energy preservation or entropy stability. While these techniques increase the robustness of DG methods significantly, they are also computationally more demanding than standard weak form nodal DG methods. We present several implementation techniques to improve the efficiency of flux differencing DG methods that use tensor product quadrilateral or hexahedral elements, in 2D or 3D, respectively. Focus is mostly given to CPUs and DG methods for the compressible Euler equations, although these techniques are generally also useful for other physical systems, including the compressible Navier-Stokes and magnetohydrodynamics equations. We present results using two open source codes, Trixi.jl written in Julia and FLUXO written in Fortran, to demonstrate that our proposed implementation techniques are applicable to different code bases and programming languages.
- [1] . 2020. Analysis of the SBP-SAT stabilization for finite element methods part I: Linear problems. J Sci. Comput. 85, 2 (2020), 1–29.
DOI: Google ScholarDigital Library - [2] . 2019. Encapsulated high order difference operators on curvilinear non-conforming grids. J. Comput. Phys. 385 (2019), 209–224.
DOI: Google ScholarDigital Library - [3] . 2020. Assessment of split form nodal discontinuous Galerkin schemes for the LES of a low pressure turbine profile. In Direct and Large Eddy Simulation XII (ERCOFTACSeries), Vol. 27. Springer Nature, Cham, 365–371.
DOI: Google ScholarCross Ref - [4] . 2017. Julia: A fresh approach to numerical computing. SIAM Rev. 59, 1 (2017), 65–98.
DOI: Google ScholarDigital Library - [5] . 2014. Entropy stable spectral collocation schemes for the Navier-Stokes equations: Discontinuous interfaces. SIAM J. Scient. Comput. 36, 5 (2014), B835–B867.
DOI: Google ScholarDigital Library - [6] . 2018. On discretely entropy conservative and entropy stable discontinuous Galerkin methods. J. Comput. Phys. 362 (2018), 346–374.
DOI: Google ScholarDigital Library - [7] . 2019. Skew-symmetric entropy stable modal discontinuous Galerkin formulations. J Sci. Comput. 81, 1 (2019), 459–485.
DOI: Google ScholarDigital Library - [8] . 2019. Efficient entropy stable Gauss collocation methods. SIAM J. Scient. Comput. 41, 5 (2019), A2938–A2966.
DOI: Google ScholarCross Ref - [9] . 2013. Kinetic energy preserving and entropy stable finite volume schemes for compressible Euler and Navier-Stokes equations. Commun. Comput. Phys. 14, 5 (2013), 1252–1286.
DOI: Google ScholarCross Ref - [10] . 2017. Entropy stable high order discontinuous Galerkin methods with suitable quadrature rules for hyperbolic conservation laws. J. Comput. Phys. 345 (2017), 427–461.
DOI: Google ScholarDigital Library - [11] . 2018. Entropy-stable summation-by-parts discretization of the Euler equations on general curved elements. J. Comput. Phys. 356 (2018), 410–438.
DOI: Google ScholarCross Ref - [12] . 2021. Roadmap to Julia BLAS and Linear Algebra. Retrieved from https://www.youtube.com/watch?v=KQ8nvlURX4MGoogle Scholar
- [13] . 2014. A generalized framework for nodal first derivative summation-by-parts operators. J. Comput. Phys. 266 (2014), 214–239.
DOI: Google ScholarDigital Library - [14] . 2014. Review of summation-by-parts operators with simultaneous approximation terms for the numerical solution of partial differential equations. Comput. Fluids 95 (2014), 171–196.
DOI: Google ScholarCross Ref - [15] . 2013. Discretely conservative finite-difference formulations for nonlinear conservation laws in split form: Theory and boundary conditions. J. Comput. Phys. 234 (2013), 353–375.
DOI: Google ScholarDigital Library - [16] . 2012. Arbitrarily high-order accurate entropy stable essentially nonoscillatory schemes for systems of conservation laws. SIAM J. Numer. Anal. 50, 2 (2012), 544–573.
DOI: Google ScholarDigital Library - [17] . 2017. On the use of kinetic energy preserving DG-schemes for large eddy simulation. J. Comput. Phys. 350 (2017), 782–795.
DOI: Google ScholarCross Ref - [18] . 2021. Instruction tables. (
08 2021). Retrieved from https://www.agner.org/optimize/instruction_tables.pdfGoogle Scholar - [19] . 2013. A skew-symmetric discontinuous Galerkin spectral element discretization and its relation to SBP-SAT finite difference methods. SIAM J. Scient. Comput. 35, 3 (2013), A1233–A1253.
DOI: Google ScholarCross Ref - [20] . 2014. A kinetic energy preserving nodal discontinuous Galerkin spectral element method. Int. J. Numer. Meth. Fluids 76, 1 (2014), 28–50.
DOI: Google ScholarCross Ref - [21] . 2013. On the accuracy of high-order discretizations for underresolved turbulence simulations. Theor. Comput. Fluid Dynam. 27, 3-4 (2013), 221–237.
DOI: Google ScholarCross Ref - [22] . 2011. A comparison of the dispersion and dissipation errors of Gauss and Gauss-Lobatto discontinuous Galerkin spectral element methods. SIAM J. Scient. Comput. 33, 5 (2011), 2560–2579.
DOI: Google ScholarDigital Library - [23] . 2021. A novel robust strategy for discontinuous Galerkin methods in computational fluid mechanics: Why? When? What? Where? Front. Phys. 8 (2021), 612.
DOI: Google ScholarCross Ref - [24] . 2016. Split form nodal discontinuous Galerkin schemes with summation-by-parts property for the compressible euler equations. J. Comput. Phys. 327 (2016), 39–66.
DOI: Google ScholarDigital Library - [25] . 2022. On the implementation of a robust and efficient finite element-based parallel solver for the compressible Navier-Stokes equations. Comput. Meth. Appl. Mechan. Eng. 389 (
02 2022), 114250.DOI: Google ScholarCross Ref - [26] . 1983. On upstream differencing and Godunov-type schemes for hyperbolic conservation laws. SIAM Rev. 25, 1 (1983), 35–61.
DOI: Google ScholarDigital Library - [27] . 2007. Nodal Discontinuous Galerkin Methods: Algorithms, Analysis, and Applications
(Texts in Applied Mathematics , Vol. 54). Springer Science & Business Media, New York.DOI: Google ScholarCross Ref - [28] . 2020. Entropy-stable, high-order summation-by-parts discretizations without interface penalties. J Sci. Comput. 82, 2 (2020), 50.
DOI: Google ScholarDigital Library - [29] . 2016. Multidimensional summation-by-parts operators: General theory and application to simplex elements. SIAM J. Scient. Comput. 38, 4 (2016), A1935–A1958.
DOI: Google ScholarDigital Library - [30] . 2007. A flux reconstruction approach to high-order schemes including discontinuous Galerkin methods. In 18th AIAA Computational Fluid Dynamics Conference. American Institute of Aeronautics and Astronautics.
DOI: Google ScholarCross Ref - [31] . 2009. Affordable, entropy-consistent euler flux functions II: Entropy production at shocks. J. Comput. Phys. 228, 15 (2009), 5410–5436.
DOI: Google ScholarDigital Library - [32] . 2008. Formulation of kinetic energy preserving conservative schemes for gas dynamics and direct numerical simulation of one-dimensional viscous compressible flow in a shock tube using entropy and kinetic energy preserving schemes. J Sci. Comput. 34, 2 (2008), 188–208.
DOI: Google ScholarDigital Library - [33] . 1994. Fourth Order 2N-Storage Runge-Kutta Schemes.
Technical Memorandum NASA-TM-109112. NASA, NASA Langley Research Center, Hampton, VA.Google Scholar - [34] . 2020. Assessing standard and kinetic energy conserving volume fluxes in discontinuous Galerkin formulations for marginally resolved Navier-Stokes flows. Comput. Fluids 205 (2020), 104557.
DOI: Google ScholarCross Ref - [35] . 2006. Metric identities and the discontinuous spectral element method on curvilinear meshes. J Sci. Comput. 26, 3 (2006), 301–327.
DOI: Google ScholarDigital Library - [36] . 2009. Implementing Spectral Methods for Partial Differential Equations: Algorithms for Scientists and Engineers. Springer Science & Business Media, New York.
DOI: Google ScholarCross Ref - [37] . 2010. On the quadrature and weak form choices in collocation type discontinuous Galerkin spectral element methods. J Sci. Comput. 44, 2 (2010), 136–155.
DOI: Google ScholarDigital Library - [38] . 2021. FLEXI: A high order discontinuous Galerkin framework for hyperbolic-parabolic conservation laws. Comput. Math. Applic. 81 (2021), 186–219.
DOI: Google ScholarCross Ref - [39] . 1974. Finite element and finite difference methods for hyperbolic partial differential equations. In Mathematical Aspects of Finite Elements in Partial Differential Equations, (Ed.). Academic Press, New York, 195–212.Google ScholarCross Ref
- [40] . 2002. Fully discrete, entropy conservative schemes of arbitrary order. SIAM J. Numer. Anal. 40, 5 (2002), 1968–1992.
DOI: Google ScholarDigital Library - [41] . 2021. Efficient parallel 3D computation of the compressible euler equations with an invariant-domain preserving second-order finite-element scheme. ACM Trans. Parallel Comput. 8, 3 (2021), 1–30.
DOI: Google ScholarDigital Library - [42] . 2001. Finite volume approximations and strict stability for hyperbolic problems. Appl. Numer. Math. 38, 3 (2001), 237–255.
DOI: Google ScholarDigital Library - [43] . 2003. Finite volume methods, unstructured meshes and strict stability for hyperbolic problems. Appl. Numer. Math. 45, 4 (2003), 453–473.
DOI: Google ScholarDigital Library - [44] . 2021. High-order accurate entropy-stable discontinuous collocated Galerkin methods with the summation-by-parts property for compressible CFD frameworks: Scalable SSDC algorithms and flow solver. J. Comput. Phys. 424 (2021), 109844.
DOI: Google ScholarCross Ref - [45] . 2017. DifferentialEquations.jl—A performant and feature-rich ecosystem for solving differential equations in Julia. J. Open Res. Softw. 5, 1 (2017), 15.
DOI: Google ScholarCross Ref - [46] . 2017. Shallow water equations: Split-form, entropy stable, well-balanced, and positivity preserving numerical methods. GEM—Int. J. Geomath. 8, 1 (
04 2017), 85–133.DOI: Google ScholarCross Ref - [47] . 2018. Comparison of some entropy conservative numerical fluxes for the euler equations. J Sci. Comput. 76, 1 (
07 2018), 216–242.DOI: arXiv:1701.02264 Google ScholarDigital Library - [48] . 2018. Generalised Summation-by-parts Operators and Entropy Stability of Numerical Methods for Hyperbolic Balance Laws. Ph.D. Dissertation. TU Braunschweig.Google Scholar
- [49] . 2020. Entropy conserving and kinetic energy preserving numerical methods for the euler equations using summation-by-parts operators. In Spectral and High Order Methods for Partial Differential Equations ICOSAHOM 2018 (Lecture Notes in Computational Science and Engineering), , , , , and (Eds.), Vol. 134. Springer, Cham, 525–535.
DOI: Google ScholarCross Ref - [50] . 2021. SummationByPartsOperators.jl: A Julia library of provably stable semidiscretization techniques with mimetic properties. J. Open Source Softw. 6, 64 (
08 2021), 3454.DOI: Google ScholarCross Ref - [51] . 2021. Optimized Runge-Kutta methods with automatic step size control for compressible computational fluid dynamics. Commun. Appl. Math. Comput. (
11 2021).DOI: Google ScholarCross Ref - [52] . 2021. Preventing pressure oscillations does not fix local linear stability issues of entropy-based split-form high-order schemes. Commun. Appl. Math. Comput. (
08 2021).DOI: Google ScholarCross Ref - [53] . 2021. A broad class of conservative numerical methods for dispersive wave equations. Commun. Comput. Phys. 29, 4 (
02 2021), 979–1029.DOI: Google ScholarCross Ref - [54] . 2016. Summation-by-parts operators for correction procedure via reconstruction. J. Comput. Phys. 311 (
04 2016), 299–328.DOI: Google ScholarDigital Library - [55] . 2021. Reproducibility repository for efficient implementation of modern entropy stable and kinetic energy preserving discontinuous Galerkin methods for conservation laws.
DOI: Google ScholarCross Ref - [56] . 2022. Adaptive numerical simulations with trixi.jl: A case study of Julia for scientific computing. Proc. JuliaCon Conf. 1, 1 (
01 2022), 77.DOI: arXiv:2108.06476 Google ScholarCross Ref - [57] . 2020. Final Report on HLST Project OPT-DG2.
Final Report . Max-Planck-Institut für Plasmaphysik, Boltzmannstraße 2, 85748 Garching, Germany.Google Scholar - [58] . 2022. Performance analysis of relaxation Runge-Kutta methods. Int. J. High Perform. Comput. Applic. (
05 2022).DOI: Google ScholarDigital Library - [59] . 2021. On the robustness and performance of entropy stable discontinuous collocation methods. J. Comput. Phys. 426 (
2 2021), 109891.DOI: Google ScholarCross Ref - [60] . 2021. An entropy stable nodal discontinuous Galerkin method for the resistive MHD equations. Part II: Subcell finite volume shock capturing. J. Comput. Phys. 444 (2021), 110580.
DOI: Google ScholarDigital Library - [61] . 2021. A purely hyperbolic discontinuous Galerkin approach for self-gravitating gas dynamics. J. Comput. Phys. 442 (
06 2021), 110467.DOI: Google ScholarCross Ref - [62] . 2020. Preventing spurious pressure oscillations in split convective form discretization for compressible flows. J. Comput. Phys. 427 (2020), 110060.
DOI: Google ScholarCross Ref - [63] . 1997. Essentially Non-oscillatory and Weighted Essentially Non-oscillatory Schemes for Hyperbolic Conservation Laws.
Final Report NASA/CR-97-206253. NASA, Institute for Computer Applications in Science and Engineering, NASA Langley Research Center, Hampton, VA.Google Scholar - [64] . 2018. High order entropy conservative central schemes for wide ranges of compressible gas dynamics and MHD flows. J. Comput. Phys. 364 (2018), 153–185.
DOI: Google ScholarDigital Library - [65] . 2017. Skew-symmetric splitting and stability of high order central schemes. In Journal of Physics: Conference Series, Vol. 837. IOP Publishing, 012019.
DOI: Google ScholarCross Ref - [66] . 1994. Summation by parts for finite difference approximations for \(d/dx\). J. Comput. Phys. 110, 1 (1994), 47–67.
DOI: Google ScholarDigital Library - [67] . 2014. Review of summation-by-parts schemes for initial-boundary-value problems. J. Comput. Phys. 268 (2014), 17–38.
DOI: Google ScholarCross Ref - [68] . 1987. The numerical viscosity of entropy stable schemes for systems of conservation laws. I. Math. Comput. 49, 179 (1987), 91–103.
DOI: Google ScholarCross Ref - [69] . 2003. Entropy stability theory for difference approximations of nonlinear conservation laws and related time-dependent problems. Acta Numer. 12 (2003), 451–512.
DOI: Google ScholarCross Ref - [70] . 1990. Table-driven implementation of the logarithm function in IEEE floating-point arithmetic. ACM Trans. Math. Softw. 16, 4 (1990), 378–400.
DOI: Google ScholarDigital Library - [71] . 1979. Geometric conservation law and its application to flow computations on moving grids. AIAA J. 17, 10 (1979), 1030–1037.
DOI: Google ScholarCross Ref - [72] . 2009. Riemann Solvers and Numerical Methods for Fluid Dynamics: A Practical Introduction. Springer, Berlin.
DOI: Google ScholarCross Ref - [73] . 2010. LIKWID: A lightweight performance-oriented tool suite for x86 multicore environments. In 39th International Conference on Parallel Processing Workshops. IEEE, 207–216.
DOI: Google ScholarDigital Library - [74] . 2009. Roofline: An insightful visual performance model for multicore architectures. Commun. ACM 52, 4 (2009), 65–76.
DOI: Google ScholarDigital Library - [75] . 2017. An entropy stable nodal discontinuous Galerkin method for the two dimensional shallow water equations on unstructured curvilinear meshes with discontinuous bathymetry. J. Comput. Phys. 340 (2017), 200–242.
DOI: Google ScholarDigital Library - [76] . 2015. A comparison of two entropy stable discontinuous Galerkin spectral element approximations for the shallow water equations with non-constant topography. J. Comput. Phys. 301 (2015), 357–376.
DOI: Google ScholarDigital Library - [77] . 2018. A comparative study on polynomial dealiasing and split form discontinuous Galerkin schemes for under-resolved turbulence computations. J. Comput. Phys. 372 (2018), 1–21.
DOI: Google ScholarDigital Library
Index Terms
- Efficient Implementation of Modern Entropy Stable and Kinetic Energy Preserving Discontinuous Galerkin Methods for Conservation Laws
Recommendations
An entropy stable spectral vanishing viscosity for discontinuous Galerkin schemes: Application to shock capturing and LES models
AbstractWe present a stable spectral vanishing viscosity for discontinuous Galerkin schemes, with applications to turbulent and supersonic flows. The idea behind the SVV is to spatially filter the dissipative fluxes, such that it concentrates ...
Highlights- Shock-capturing method for the DGSEM with SVV-filtered artificial viscosity.
- ...
On discretely entropy conservative and entropy stable discontinuous Galerkin methods
AbstractHigh order methods based on diagonal-norm summation by parts operators can be shown to satisfy a discrete conservation or dissipation of entropy for nonlinear systems of hyperbolic PDEs [1,2]. These methods can also be interpreted as ...
Entropy stable modal discontinuous Galerkin schemes and wall boundary conditions for the compressible Navier-Stokes equations
Highlights- We derive more general “modal” entropy stable discretizations of viscous terms in the compressible Navier-Stokes.
AbstractEntropy stable schemes ensure that physically meaningful numerical solutions also satisfy a semi-discrete entropy inequality under appropriate boundary conditions. In this work, we describe a discretization of viscous terms in the ...
Comments