Summation-by-parts operators and high-order quadrature (Q455831)

This paper highlights that a quadrature rule is particularly well suited for high-order summation-by-parts (SBP) finite-difference methods. The high-order quadrature in question is based on the weight matrix that forms a part of the definition of SBP operators. This result is somewhat surprising, because the accuracy of the quadrature induced by the weight matrix is not explicitly part of the SBP definition. The relationship between SBP operators and the quadrature has not been discussed previously in the literature. In this paper, this relationship is presented and its importance is demonstrated. The authors derive conditions on the quadrature weights for the class of trapezoid rules with end corrections. They use these conditions to establish the accuracy of SBP-based quadrature. They consider the impact of coordinate transformations on diagonal-norm SBP quadrature and show that the quadrature remains accurate on curvilinear multi-dimensional domains. Finally, the theoretical results are verified with several numerical examples.