Gaussian FSBP operators: Comparison and application to numerical methods for hyperbolic conservation laws
arXiv:2607.13224
The paper compares open and closed function-space summation-by-parts (FSBP) operators built with generalized Gaussian quadrature and introduces exact extrapolation operators for open node distributions, demonstrating improved efficiency and accuracy in numerical solutions of hyperbolic conservation laws.
Abstract
Function-space summation-by-parts (FSBP) operators enable conservative and energy-stable numerical methods for hyperbolic conservation laws based on general, non-polynomial approximation spaces. Recent works show that using generalized Gaussian quadrature significantly reduces the number of grid points required compared to existing constructions that have mostly focused on equidistant grids. In this paper, we compare open and closed FSBP operators constructed with generalized Gaussian quadratures and apply them to numerically solve hyperbolic conservation laws. Furthermore, to support open node distributions, we extend the FSBP framework by introducing function-space exact extrapolation operators and operationalize them in numerical schemes for solving hyperbolic conservation laws. Our numerical experiments include the one-dimensional linear advection, non-viscous Burgers, and compressible Euler equations of gas dynamics. We observe that applying FSBP operators in numerical schemes can improve efficiency and accuracy. Notably, we demonstrate these advantages in more challenging time-dependent settings compared to other recent works on Gaussian FSBP operators.