Simultaneous approximation terms and functional accuracy for diffusion problems discretized with multidimensional summation-by-parts operators

Abstract

Several types of simultaneous approximation term (SAT) for diffusion problems discretized with diagonal-norm multidimensional summation-by-parts (SBP) operators are analyzed based on a common framework. Conditions under which the SBP-SAT discretizations are consistent, conservative, adjoint consistent, and energy stable are presented. For SATs leading to primal and adjoint consistent discretizations, the error in output functionals is shown to be of order $h^{2p}$ when a degree $p$ multidimensional SBP operator is used to discretize the spatial derivatives. SAT penalty coefficients corresponding to various discontinuous Galerkin fluxes developed for elliptic partial differential equations are identified. We demonstrate that the original method of Bassi and Rebay, the modified method of Bassi and Rebay, and the symmetric interior penalty method are equivalent when implemented with SBP diagonal-E operators that have diagonal norm matrix, e.g., the Legendre-Gauss-Lobatto SBP operator in one space dimension. Similarly, the local discontinuous Galerkin and the compact discontinuous Galerkin schemes are equivalent for this family of operators. The analysis remains valid on curvilinear grids if a degree $\le p+1$ bijective polynomial mapping from the reference to physical elements is used. Numerical experiments with the two-dimensional Poisson problem support the theoretical results.

arXiv comment44 pages, 11 figures

References in corpus (7)

• Review of Summation-by-parts schemes for initial-boundary-value problems
• On discretely entropy conservative and entropy stable discontinuous Galerkin methods
• Programming Languages for Scientific Computing
• Summation-By-Parts Operators and High-Order Quadrature
• Summation-by-parts operators for correction procedure via reconstruction
• A dual consistent finite difference method with narrow stencil second derivative operators
• A unifying algebraic framework for discontinuous Galerkin and flux reconstruction methods based on the summation-by-parts property

Cited by in corpus (1)

• Stability and Functional Superconvergence of Narrow-Stencil Second-Derivative Generalized Summation-By-Parts Discretizations