David C. Del Rey Fernández

Daniel Bach, Andrés M. Rueda-Ramírez, Eric Sonnendrücker et al.

We demonstrate that we can carry over the strategy of Finite Element Exterior Calculus (FEEC) to Summation-by-Parts (SBP) Finite Difference (FD) methods to achieve divergence- and curl-free discretizations. This is not obvious at first sight, as for SBP-FD no basis functions are known, but only values and derivatives at points. The key is a remarkable analytic relationship that enables us to construct compatible operators using integral and nodal degrees of freedom. Pre-existing SBP-FD matrix operators can then be used to obtain nodal values from the integral degrees of freedom to derive a scheme with the desired properties.

Lucas Friedrich, Gero Schnücke, Andrew R. Winters et al.

This work examines the development of an entropy conservative (for smooth solutions) or entropy stable (for discontinuous solutions) space-time discontinuous Galerkin (DG) method for systems of non-linear hyperbolic conservation laws. The resulting numerical scheme is fully discrete and provides a bound on the mathematical entropy at any time according to its initial condition and boundary conditions. The crux of the method is that discrete derivative approximations in space and time are summation-by-parts (SBP) operators. This allows the discrete method to mimic results from the continuous entropy analysis and ensures that the complete numerical scheme obeys the second law of thermodynamics. Importantly, the novel method described herein does not assume any exactness of quadrature in the variational forms that naturally arise in the context of DG methods. Typically, the development of entropy stable schemes is done on the semi-discrete level ignoring the temporal dependence. In this work we demonstrate that creating an entropy stable DG method in time is similar to the spatial discrete entropy analysis, but there are important (and subtle) differences. Therefore, we highlight the temporal entropy analysis throughout this work. For the compressible Euler equations, the preservation of kinetic energy is of interest besides entropy stability. The construction of kinetic energy preserving (KEP) schemes is, again, typically done on the semi-discrete level similar to the construction of entropy stable schemes. We present a generalization of the KEP condition from Jameson to the space-time framework and provide the temporal components for both entropy stability and kinetic energy preservation. The properties of the space-time DG method derived herein is validated through numerical tests for the compressible Euler equations.

Lucas Friedrich, Andrew R. Winters, David C. Del Rey Fernández et al.

This work presents an entropy stable discontinuous Galerkin (DG) spectral element approximation for systems of non-linear conservation laws with general geometric (h) and polynomial order (p) non-conforming rectangular meshes. The crux of the proofs presented is that the nodal DG method is constructed with the collocated Legendre-Gauss-Lobatto nodes. This choice ensures that the derivative/mass matrix pair is a summation-by-parts (SBP) operator such that entropy stability proofs from the continuous analysis are discretely mimicked. Special attention is given to the coupling between nonconforming elements as we demonstrate that the standard mortar approach for DG methods does not guarantee entropy stability for non-linear problems, which can lead to instabilities. As such, we describe a precise procedure and modify the mortar method to guarantee entropy stability for general non-linear hyperbolic systems on h/p non-conforming meshes. We verify the high-order accuracy and the entropy conservation/stability of fully non-conforming approximation with numerical examples.

Anita Gjesteland, Sigrun Ortleb, Salim Elghawi et al.

We develop an unconditionally energy-stable tensor-product space-time discretization framework for the solution of a linear kinetic transport equation in one space dimension. The kinetic equation is a simplified model of radiative transfer formulated as a hyperbolic balance law in diffusive scaling for a particle distribution function of the independent variables space, time and velocity. Our numerical discretization is based on the well-known technique of micro-macro decomposition which results in a system of balance laws for equilibrium and non-equilibrium quantities and facilitates preservation of the asymptotic limit for vanishing scaling parameters at the discrete level. We prove fully discrete stability and asymptotic preservation for general spatial and temporal discretizations having the summation-by-parts property. A new provably energy-stable Dirichlet boundary treatment for the micro-macro decomposed system is developed based on the introduction of simultaneous approximation terms. Numerical results show convergence for smooth problems and demonstrate energy stability of the proposed boundary treatment.

David C. Del Rey Fernández, Jason E. Hicken, David W. Zingg

This paper is concerned with the accurate, conservative, and stable imposition of boundary conditions and inter-element coupling for multi-dimensional summation-by-parts (SBP) finite-difference operators. More precisely, the focus is on diagonal-norm SBP operators that are not based on tensor products and are applicable to unstructured grids composed of arbitrary elements. We show how penalty terms --- simultaneous approximation terms (SATs) --- can be adapted to discretizations based on multi-dimensional SBP operators to enforce boundary and interface conditions. A general SAT framework is presented that leads to conservative and stable discretizations of the variable-coefficient advection equation. This framework includes the case where there are no nodes on the boundary of the SBP element at which to apply penalties directly. This is an important generalization, because elements analogous to Legendre-Gauss collocation, \ie without boundary nodes, typically have higher accuracy for the same number of degrees of freedom. Symmetric and upwind examples of the general SAT framework are created using a decomposition of the symmetric part of an SBP operator; these particular SATs enable the pointwise imposition of boundary and inter-element conditions. We illustrate the proposed SATs using triangular-element SBP operators with and without nodes that lie on the boundary. The accuracy, conservation, and stability properties of the resulting SBP-SAT discretizations are verified using linear advection problems with spatially varying divergence-free velocity fields.

Jason E. Hicken, David C. Del Rey Fernández, David W. Zingg

Summation-by-parts (SBP) finite-difference discretizations share many attractive properties with Galerkin finite-element methods (FEMs), including time stability and superconvergent functionals; however, unlike FEMs, SBP operators are not completely determined by a basis, so the potential exists to tailor SBP operators to meet different objectives. To date, application of high-order SBP discretizations to multiple dimensions has been limited to tensor product domains. This paper presents a definition for multi-dimensional SBP finite-difference operators that is a natural extension of one-dimensional SBP operators. Theoretical implications of the definition are investigated for the special case of a diagonal norm (mass) matrix. In particular, a diagonal-norm SBP operator exists on a given domain if and only if there is a cubature rule with positive weights on that domain and the polynomial-basis matrix has full rank when evaluated at the cubature nodes. Appropriate simultaneous-approximation terms are developed to impose boundary conditions weakly, and the resulting discretizations are shown to be time stable. Concrete examples of multi-dimensional SBP operators are constructed for the triangle and tetrahedron; similarities and differences with spectral-element and spectral-difference methods are discussed. An assembly process is described that builds diagonal-norm SBP operators on a global domain from element-level operators. Numerical results of linear advection on a doubly periodic domain demonstrate the accuracy and time stability of the simplex operators.

David C. Del Rey Fernández, David W. Zingg

The comprehensive generalization of summation-by-parts of Del Rey Fernández et al.\ (J. Comput. Phys., 266, 2014) is extended to approximations of second derivatives with variable coefficients. This enables the construction of second-derivative operators with one or more of the following characteristics: i) non-repeating interior stencil, ii) nonuniform nodal distributions, and iii) exclusion of one or both boundary nodes. Definitions are proposed that give rise to generalized SBP operators that result in consistent, conservative, and stable discretizations of PDEs with or without mixed derivatives. It is proven that such operators can be constructed using a correction to the application of the first-derivative operator twice that is the same as used for the constant-coefficient operator. Moreover, for operators with a repeating interior stencil, a decomposition is proposed that makes the application of such operators particularly simple. A number of novel operators are constructed, including operators on pseudo-spectral nodal distributions and operators that have a repeating interior stencil, but unequal nodal spacing near boundaries. The various operators are compared to the application of the first-derivative operator twice in the context of the linear convection-diffusion equation with constant and variable coefficients.