Summation-by-parts operators for general function spaces: Flux differencing, entropy stability, and accuracy

Abstract

Summation-by-parts (SBP) operators are a cornerstone of provably stable high-order methods for hyperbolic conservation laws. For nonlinear problems, the flux-differencing framework of Fisher and Carpenter combines such SBP operators with entropy-conservative two-point fluxes to obtain semi-discretizations that are conservative in the sense of Lax--Wendroff and discretely entropy-conservative or stable. So far, this framework has been developed exclusively for operators that are exact for polynomials. Function-space SBP (FSBP) operators lift this restriction by being exact on a general finite-dimensional function space that can be tailored to the problem at hand. In this work, we extend flux differencing to FSBP operators for general (potentially non-polynomial) function spaces. We prove that, whenever the operator has a diagonal norm and its function space contains the constants, the resulting semi-discretization is conservative in the sense of Lax--Wendroff, semi-discretely entropy-conservative within a single element, and entropy-stable across element interfaces. We also provide the first theoretical accuracy analysis of FSBP operators: We characterize their design order through an algebraic criterion on the underlying function space, prove that suitable non-polynomial spaces achieve the same design order as polynomial spaces of equal dimension, and show that the flux-differencing semi-discretization retains this order. Notably, the entropy-conservative two-point flux is reused unchanged; only the underlying function space differs from the polynomial case. Numerical experiments for the Burgers' equation and the compressible Euler equations in one and two dimensions confirm these properties.

A generating function identity for finite and symmetric multiple zeta values via Dougall's bilateral ₂H₂ summation
Takumi Maesaka, Masataka Ono
arXiv13 Oct 2026PreprintCC BY 4.0Full textAlgebra & number theoryarXiv
A summation rule for the four-coupling beta function of four-dimensional nonplanar Euclidean scalar φ⁴ theory under factorial coefficient bounds, and a domain on which the summed beta function is well defined
Puyu Zeng
arXiv1 Sep 2026PreprintCC BY 4.0Full textMathematical physics Statistical & nonlinear physicsarXiv
Stability of Serrin-Type Overdetermined Problems for the p-Laplacian and Related Operators
Yingfang Zhang, Wenming Zou
arXiv20 Oct 2026PreprintarXiv licenceFull textAnalysis & PDEsarXiv
Function spaces for weighted subcoercive operators
Tommaso Bruno, Mattia Calzi, Marco M. Peloso
arXiv16 Sep 2026PreprintarXiv licenceFull textAnalysis & PDEsarXiv
On the stability of proximal operators in Wasserstein spaces under different notions of convexity
Simone Di Marino, Sara Farinelli, Emanuele Naldi
arXiv24 Sep 2026PreprintarXiv licenceFull textAnalysis & PDEs Optimisation & controlarXiv

Cite

BibTeX RIS

BibTeX

@misc{haase2026summationbyparts,
  title = {{Summation-by-parts operators for general function spaces: Flux differencing, entropy stability, and accuracy}},
  author = {Henry M. Haase and Andrew R. Winters and Jan Glaubitz},
  howpublished = {arXiv},
  year = {2026},
  month = oct,
  doi = {10.48550/arxiv.2610.30272},
  eprint = {2610.30272},
  archivePrefix = {arXiv},
  url = {https://arxiv.org/abs/2610.30272},
}

RIS

TY  - UNPB
TI  - Summation-by-parts operators for general function spaces: Flux differencing, entropy stability, and accuracy
AU  - Henry M. Haase
AU  - Andrew R. Winters
AU  - Jan Glaubitz
PY  - 2026
DA  - 2026/10/21
PB  - arXiv
DO  - 10.48550/arxiv.2610.30272
UR  - https://arxiv.org/abs/2610.30272
AB  - Summation-by-parts (SBP) operators are a cornerstone of provably stable high-order methods for hyperbolic conservation laws. For nonlinear problems, the flux-differencing framework of Fisher and Carpenter combines such SBP operators with entropy-conservative two-point fluxes to obtain semi-discretizations that are conservative in the sense of Lax--Wendroff and discretely entropy-conservative or stable. So far, this framework has been developed exclusively for operators that are exact for polynomials. Function-space SBP (FSBP) operators lift this restriction by being exact on a general finite-dimensional function space that can be tailored to the problem at hand. In this work, we extend flux differencing to FSBP operators for general (potentially non-polynomial) function spaces. We prove that, whenever the operator has a diagonal norm and its function space contains the constants, the resulting semi-discretization is conservative in the sense of Lax--Wendroff, semi-discretely entropy-conservative within a single element, and entropy-stable across element interfaces. We also provide the first theoretical accuracy analysis of FSBP operators: We characterize their design order through an algebraic criterion on the underlying function space, prove that suitable non-polynomial spaces achieve the same design order as polynomial spaces of equal dimension, and show that the flux-differencing semi-discretization retains this order. Notably, the entropy-conservative two-point flux is reused unchanged; only the underlying function space differs from the polynomial case. Numerical experiments for the Burgers' equation and the compressible Euler equations in one and two dimensions confirm these properties.
ER  -

Abstract

Related

Cite