8,820 research outputs found
Summation-By-Parts Operators and High-Order Quadrature
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 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
norm for functions, and multi-dimensional SBP discretizations
accurately mimic the divergence theorem.Comment: 18 pages, 3 figure
Review of Summation-by-parts schemes for initial-boundary-value problems
High-order finite difference methods are efficient, easy to program, scales
well in multiple dimensions and can be modified locally for various reasons
(such as shock treatment for example). The main drawback have been the
complicated and sometimes even mysterious stability treatment at boundaries and
interfaces required for a stable scheme. The research on summation-by-parts
operators and weak boundary conditions during the last 20 years have removed
this drawback and now reached a mature state. It is now possible to construct
stable and high order accurate multi-block finite difference schemes in a
systematic building-block-like manner. In this paper we will review this
development, point out the main contributions and speculate about the next
lines of research in this area
High-order numerical methods for 2D parabolic problems in single and composite domains
In this work, we discuss and compare three methods for the numerical
approximation of constant- and variable-coefficient diffusion equations in both
single and composite domains with possible discontinuity in the solution/flux
at interfaces, considering (i) the Cut Finite Element Method; (ii) the
Difference Potentials Method; and (iii) the summation-by-parts Finite
Difference Method. First we give a brief introduction for each of the three
methods. Next, we propose benchmark problems, and consider numerical tests-with
respect to accuracy and convergence-for linear parabolic problems on a single
domain, and continue with similar tests for linear parabolic problems on a
composite domain (with the interface defined either explicitly or implicitly).
Lastly, a comparative discussion of the methods and numerical results will be
given.Comment: 45 pages, 12 figures, in revision for Journal of Scientific Computin
- …