Summation-by-parts operators and high-order quadrature

From MaRDI portal
(Redirected from Publication:455831)



Abstract: Summation-by-parts (SBP) operators are finite-difference operators that mimic integration by parts. This property can be useful in constructing energy-stable discretizations of partial differential vequations. SBP operators are defined by a weight matrix and a difference operator, with the latter designed to approximate d/dx to a specified order of accuracy. The accuracy of the weight matrix as a quadrature rule is not explicitly part of the SBP definition. We show that SBP weight matrices are related to trapezoid rules with end corrections whose accuracy matches the corresponding difference operator at internal nodes. The accuracy of SBP quadrature extends to curvilinear domains provided the Jacobian is approximated with the same SBP operator used for the quadrature. This quadrature has significant implications for SBP-based discretizations; for example, the discrete norm accurately approximates the L2 norm for functions, and multi-dimensional SBP discretizations accurately mimic the divergence theorem.


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.




Cited in
(67)








This page was built for publication: Summation-by-parts operators and high-order quadrature

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q455831)