Entropy Stable Discontinuous Galerkin Research Articles

Entropy stable discontinuous Galerkin methods for the shallow water equations with subcell positivity preservation

AbstractHigh order schemes are known to be unstable in the presence of shock discontinuities or under‐resolved solution features, and have traditionally required additional filtering, limiting, or artificial viscosity to avoid solution blow up. Entropy stable schemes address this instability by ensuring that physically relevant solutions satisfy a semi‐discrete entropy inequality independently of discretization parameters. However, additional measures must be taken to ensure that solutions satisfy physical constraints such as positivity. In this work, we present a high order entropy stable discontinuous Galerkin (ESDG) method for the nonlinear shallow water equations (SWE) on two‐dimensional (2D) triangular meshes which preserves the positivity of the water heights. The scheme combines a low order positivity preserving method with a high order entropy stable method using convex limiting. This method is entropy stable and well‐balanced for fitted meshes with continuous bathymetry profiles.

A positivity preserving and oscillation-free entropy stable discontinuous Galerkin scheme for the reactive Euler equations

The reactive Euler equations are a basic model for fluid flows with chemical reactions. In this work, we construct a high order positivity preserving and oscillation-free entropy stable discontinuous Galerkin (DG) scheme for the reactive Euler equations. The main ingredients of the scheme include (i) entropy preserving and entropy stable fluxes to achieve entropy stability, (ii) artificial damping terms to restrain spurious oscillations near the shocks, and (iii) positivity preserving limiters to guarantee the positivity of solutions. These ingredients are compatible with each other so that our scheme simultaneously enjoys the properties of entropy stable, oscillation-free and positivity preserving. Another distinctive feature of our scheme is that it is entropy stable for both the thermodynamic and mathematical entropies. Numerical experiments validate the designed high convergence orders of the scheme and demonstrate its good performances for discontinuous problems.

An Entropy Stable Discontinuous Galerkin Method for the Two-Layer Shallow Water Equations on Curvilinear Meshes

We present an entropy stable nodal discontinuous Galerkin spectral element method (DGSEM) for the two-layer shallow water equations on two dimensional curvilinear meshes. We mimic the continuous entropy analysis on the semi-discrete level with the DGSEM constructed on Legendre–Gauss–Lobatto (LGL) nodes. The use of LGL nodes endows the collocated nodal DGSEM with the summation-by-parts property that is key in the discrete analysis. The approximation exploits an equivalent flux differencing formulation for the volume contributions, which generate an entropy conservative split-form of the governing equations. A specific combination of a numerical surface flux and discretization of the nonconservative terms is then applied to obtain a high-order path-conservative scheme that is entropy conservative. Furthermore, we find that this combination yields an analogous discretization for the pressure and nonconservative terms such that the numerical method is well-balanced for discontinuous bathymetry on curvilinear domains. Dissipation is added at the interfaces to create an entropy stable approximation that satisfies the second law of thermodynamics in the discrete case, while maintaining the well-balanced property. We conclude with verification of the theoretical findings through numerical tests and demonstrate results about convergence, entropy stability and well-balancedness of the scheme.

High order entropy stable discontinuous Galerkin spectral element methods through subcell limiting

Subcell limiting strategies for discontinuous Galerkin spectral element methods do not provably satisfy a semi-discrete cell entropy inequality. In this work, we introduce an extension to the subcell and monolithic convex limiting strategies [12,11,2] that satisfies the semi-discrete cell entropy inequality by formulating the limiting factors as solutions to an optimization problem. The optimization problem is efficiently solved using a deterministic greedy algorithm. We also discuss the extension of the proposed subcell limiting strategy to preserve general convex constraints. Numerical experiments confirm that the proposed limiting strategy preserves high-order accuracy for smooth solutions and satisfies the cell entropy inequality.

Entropy Stable Discontinuous Galerkin Schemes for Two-Fluid Relativistic Plasma Flow Equations

Entropy Stable Discontinuous Galerkin Schemes for Two-Fluid Relativistic Plasma Flow Equations

A positivity preserving strategy for entropy stable discontinuous Galerkin discretizations of the compressible Euler and Navier-Stokes equations

High-order entropy-stable discontinuous Galerkin methods for the compressible Euler and Navier-Stokes equations require the positivity of thermodynamic quantities in order to guarantee their well-posedness. In this work, we introduce a positivity limiting strategy for entropy-stable discontinuous Galerkin discretizations constructed by blending high order solutions with a low order positivity-preserving discretization. The proposed low order discretization is semi-discretely entropy stable, and the proposed limiting strategy is positivity preserving for the compressible Euler and Navier-Stokes equations. Numerical experiments confirm the high order accuracy and robustness of the proposed strategy.

Entropy stable discontinuous Galerkin methods for balance laws in non-conservative form: Applications to the Euler equations with gravity

In this work a non-conservative balance law formulation is considered that encompasses the rotating, compressible Euler equations for dry atmospheric flows. We develop a semi-discretely entropy stable discontinuous Galerkin method on curvilinear meshes using a generalization of flux differencing for numerical fluxes in fluctuation form. The method uses the skew-hybridized formulation of the element operators to ensure that, even in the presence of under-integration on curvilinear meshes, the resulting discretization is entropy stable. Several atmospheric flow test cases in one, two, and three dimensions confirm the theoretical entropy stability results as well as show the high-order accuracy and robustness of the method.

On the Entropy Projection and the Robustness of High Order Entropy Stable Discontinuous Galerkin Schemes for Under-Resolved Flows

High order entropy stable schemes provide improved robustness for computational simulations of fluid flows. However, additional stabilization and positivity preserving limiting can still be required for variable-density flows with under-resolved features. We demonstrate numerically that entropy stable Discontinuous Galerkin (DG) methods which incorporate an “entropy projection” are less likely to require additional limiting to retain positivity for certain types of flows. We conclude by investigating potential explanations for this observed improvement in robustness.

Entropy stable discontinuous Galerkin schemes for the special relativistic hydrodynamics equations

This article presents entropy stable discontinuous Galerkin numerical schemes for equations of special relativistic hydrodynamics with the ideal equation of state. The numerical schemes use the summation by parts (SBP) property of the Gauss-Lobatto quadrature rules. To achieve entropy stability for the scheme, we use two-point entropy conservative numerical flux inside the cells and a suitable entropy stable numerical flux at the cell interfaces. The resulting semi-discrete scheme is then shown to be entropy stable. Time discretization is performed using SSP Runge-Kutta methods. Several numerical test cases are presented to validate the accuracy and stability of the proposed schemes.

An entropy–stable p–adaptive nodal discontinuous Galerkin for the coupled Navier–Stokes/Cahn–Hilliard system

We develop a novel entropy–stable discontinuous Galerkin approximation of the incompressible Navier–Stokes/Cahn–Hilliard system for p–non–conforming elements. This work constitutes an evolution of the work presented by Manzanero et al. ((2020) [10]), as it extends the discrete analysis into supporting p–adaptation (p–refinement/coarsening). The scheme is based on the summation–by–parts simultaneous–approximation term property along with Gauss–Lobatto points and suitable numerical fluxes. The p–non–conforming elements are connected through the classic mortar method, the use of central fluxes for the inviscid terms, and the BR1 scheme with additional dissipation for the viscous fluxes. The scheme is proven to retain its properties of the original conforming scheme when transitioning to p–non–conforming elements and to mimic the continuous entropy analysis of the model. We focus on dynamic polynomial adaptation as the applications of interest are unsteady multiphase flows. In this work, we introduce a heuristic adaptation criterion that depends on the location of the interface between the different phases and utilises the convection velocity to predict the movement of the interface. The scheme is verified to be total phase conserving, entropy–stable and freestream preserving for curvilinear p–non–conforming meshes. We also present the results for a rising bubble simulation and we show that for the same accuracy we get a ×2 to ×6 reduction in the degrees of freedom and a 41% reduction in the computational time. We compare our results for the three–dimensional dam break test case against experimental and numerical data and we show that a ×4.3 to ×9.5 reduction of the degrees of freedom and a 51% reduction in the computational time can be achieved compared to the p–uniform solution.

An entropy–stable discontinuous Galerkin approximation of the Spalart–Allmaras turbulence model for the compressible Reynolds Averaged Navier–Stokes equations

We present an entropy–stable formulation for the compressible Reynolds Averaged Navier–Stokes (RANS) Discontinuous Galerkin (DG) equations and the Spalart–Allmaras one–equation closure. The model is designed to satisfy an entropy law, which includes free– and no–slip wall boundary conditions. We then construct a high–order DG approximation of the model that satisfies the summation–by–parts simultaneous–approximation–term (SBP-SAT) property. With the help of a discrete stability analysis, we construct two approximations: a kinetic energy preserving scheme based on Pirozzoli's two–point flux and a thermodynamic entropy conserving one based on Chandrashekar's split–form. The schemes are applicable to, and the stability proofs hold for, three–dimensional unstructured meshes with curvilinear hexahedral elements. We test the convergence of the schemes on a manufactured solution for increasing polynomial orders and mesh refinement levels, to then assess their numerical stability by propagating a flow from a random initial condition, and finally solve the flow around a two–dimensional flat plate and a NACA 0012 airfoil, comparing numerical results with those available in the literature. The proposed schemes are entropy–stable, and provide accurate solutions for the selected test cases.

High order entropy stable and positivity-preserving discontinuous Galerkin method for the nonlocal electron heat transport model

The nonlocal electron heat transport model in laser heated plasmas plays a crucial role in inertial confinement fusion (ICF), and it is important to solve it numerically in an accurate and robust way. In this paper, we first develop an one-dimensional high-order entropy stable discontinuous Galerkin method for the nonlocal electron heat transport model. We further design our DG scheme to have the positivity-preserving property by using the scaling positivity-preserving limiter, which is shown, by a computer-aided proof, to have no extra time step constraint than that required by L2 stability. Next, we extend our one-dimensional scheme to two-dimensions on rectangular meshes and tensor product polynomial spaces. Numerical examples are given to verify the high-order accuracy and positivity-preserving properties of our scheme. By comparing the local and nonlocal electron heat transport models, we also observe more physical phenomena such as the flux reduction and the preheat effect from the nonlocal model.

Entropy Stable Discontinuous Galerkin Finite Element Method with Multi-Dimensional Slope Limitation for Euler Equations

Abstract We present an entropy stable Discontinuous Galerkin (DG) finite element method to approximate systems of 2-dimensional symmetrizable conservation laws on unstructured grids. The scheme is constructed using a combination of entropy conservative fluxes and entropy-stable numerical dissipation operators. The method is designed to work on structured as well as on unstructured meshes. As solutions of hyperbolic conservation laws can develop discontinuities (shocks) in finite time, we include a multidimensional slope limitation step to suppress spurious oscillations in the vicinity of shocks. The numerical scheme has two steps: the first step is a finite element calculation which includes calculations of fluxes across the edges of the elements using 1-D entropy stable solver. The second step is a procedure of stabilization through a truly multi-dimensional slope limiter. We compared the Entropy Stable Scheme (ESS) versus Roe’s solvers associated with entropy corrections and Osher’s solver. The method is illustrated by computing solution of the two stationary problems: a regular shock reflection problem and a 2-D flow around a double ellipse at high Mach number.

A Positivity Preserving Strategy for Entropy Stable Discontinuous Galerkin Discretizations of the Compressible Euler and Navier-Stokes Equations

High-order entropy-stable discontinuous Galerkin methods for the compressible Euler and Navier-Stokes equations require the positivity of thermodynamic quantities in order to guarantee their well-posedness. In this work, we introduce a positivity limiting strategy for entropy-stable discontinuous Galerkin discretizations based on convex limiting. The key ingredient in the limiting procedure is a low order positivity-preserving discretization based on graph viscosity terms. The proposed limiting strategy is both positivity preserving and discretely entropy-stable for the compressible Euler and Navier-Stokes equations. Numerical experiments confirm the high order accuracy and robustness of the proposed strategy.

Entropy stable modal discontinuous Galerkin schemes and wall boundary conditions for the compressible Navier-Stokes equations

Entropy 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 compressible Navier-Stokes equations which enables a simple and explicit imposition of entropy stable no-slip and reflective (symmetry) wall boundary conditions for discontinuous Galerkin (DG) discretizations. Specifically, we derive methods for imposing adiabatic no-slip and reflective (symmetry) boundary conditions for modal entropy stable DG formulations which preserve a semi-discrete entropy inequality. Numerical results confirm the robustness and accuracy of the proposed approaches.

Entropy Stable Discontinuous Galerkin Methods for Nonlinear Conservation Laws on Networks and Multi-Dimensional Domains

We present a high-order entropy stable discontinuous Galerkin method for nonlinear conservation laws on both multi-dimensional domains and on networks constructed from one-dimensional domains. These methods utilize treatments of multi-dimensional interfaces and network junctions which retain entropy stability when coupling together entropy stable discretizations. Numerical experiments verify the stability of the proposed schemes, and comparisons with fully 2D implementations demonstrate the accuracy of each type of junction treatment.

Entropy stable discontinuous Galerkin methods for ten-moment Gaussian closure equations

In this article, we propose high order discontinuous Galerkin entropy stable schemes for ten-moment Gaussian closure equations, based on the suitable quadrature rules (see [1]). The key components of the proposed schemes are the use of an entropy conservative numerical flux [2] in each cell and an appropriate entropy stable numerical flux at the cell edges. These fluxes are then used in the entropy stable DG framework of [1] to obtain entropy stability of the semi-discrete schemes. We also extend these schemes to a source term that models plasma laser interaction. For the time discretization, we use strong stability preserving methods. The proposed schemes are then tested on several test cases to demonstrate stability, accuracy and robustness.

High-order entropy stable discontinuous Galerkin methods for the shallow water equations: Curved triangular meshes and GPU acceleration

We present a high-order entropy stable discontinuous Galerkin (ESDG) method for the two dimensional shallow water equations (SWE) on curved triangular meshes. The presented scheme preserves a semi-discrete entropy inequality and remains well-balanced for continuous bathymetry profiles. We provide numerical experiments which confirm the high-order accuracy and theoretical properties of the scheme, and compare the presented scheme to an entropy stable scheme based on simplicial summation-by-parts (SBP) finite difference operators. Finally, we report the computational performance of an implementation on Graphics Processing Units (GPUs) and provide comparisons to existing GPU-accelerated implementations of high-order DG methods on quadrilateral meshes.

A provably entropy stable subcell shock capturing approach for high order split form DG for the compressible Euler equations

The main result in this paper is a provably entropy stable shock capturing approach for the high order entropy stable Discontinuous Galerkin Spectral Element Method (DGSEM) based on a hybrid blending with a subcell low order variant. Since it is possible to rewrite a high order summation-by-parts (SBP) operator into an equivalent conservative finite volume form, we were able to design a low order scheme directly with the Legendre-Gauss-Lobatto (LGL) nodes that is compatible to the discrete entropy analysis used for the proof of the entropy stable DGSEM. Furthermore, we present a hybrid low order/high order discretisation where it is possible to seamlessly blend between the two approaches, while still being provably entropy stable. With tensor products and careful design of the low order scheme on curved elements, we are able to extend the approach to three spatial dimensions on unstructured curvilinear hexahedral meshes. We validate our theoretical findings and demonstrate convergence order for smooth problems, conservation of the primary quantities and discrete entropy stability for an arbitrary blending on curvilinear grids. In practical simulations, we connect the blending factor to a local troubled element indicator that provides the control of the amount of low order dissipation injected into the high order scheme. We modified an existing shock indicator, which is based on the modal polynomial representation, to our provably stable hybrid scheme. The aim is to reduce the impact of the parameters as good as possible. We describe our indicator in detail and demonstrate its robustness in combination with the hybrid scheme, as it is possible to compute all the different test cases without changing the indicator. The test cases include e.g. the double Mach reflection setup, forward and backward facing steps with shock Mach numbers up to 100. The proposed approach is relatively straight forward to implement in an existing entropy stable DGSEM code as only modifications local to an element are necessary.

Discontinuous Galerkin Method with an Entropic Slope Limiter for Euler Equations

The variational approach to obtaining equations of the entropy stable discontinuous Galerkin method is generalized. It is shown how the monotonicity property can be incorporated into this approach. As applied to Euler equations, the entropic slope limiter, a new effective approximate method for the studied approach, is designed. It guarantees the monotonicity of the numerical solution, as well as the nonnegativity of the pressure and entropy production for each finite element. This method is successfully tested on some well-known gas dynamics model problems.