An algorithmic exploration of the existence of high-order summation by parts operators with diagonal norm

DOI10.1007/S10915-016-0211-ZMaRDI QIDQ2014030zbMATHOpenAlexFDO

Authors Nathan Albin, Joshua Klarmann

Publication date 10 August 2017

Published in Journal of Scientific Computing (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1403.5750

zbMATH Keywords

algorithmnumerical experimenthigh-order finite difference methodssummation by partsdiagonal energy norm

Mathematics Subject Classification ID

Finite difference methods for initial value and initial-boundary value problems involving PDEs (65M06) Initial-boundary value problems for first-order hyperbolic equations (35L04)

Abstract: This paper explores a common class of diagonal-norm summation by parts (SBP) operators found in the literature, which can be parameterized by an integer triple $(s,t,r)$ representing the interior order of accuracy ($2s)$, the boundary order of accuracy ($t$), and the dimension of the boundary closure ($r$). There is no simple formula for determining whether or not an SBP operator exists for a given triple of parameters. Instead, one must check that certain compatibility conditions are met: namely that a particular linear system of equations has a positive solution. Partly because of the complexity involved, not much is known about diagonal-norm SBP operators with $2s>10$. By utilizing a new algorithm for answering the question "Does an SBP operator exist for the parameters $(s,t,r)$?", it is possible to explore the existence of SBP operators with high order accuracy, and previously unknown SBP operators with interior order of accuracy as large as $2s=30$ are found. Additionally, a method for optimizing the spectral radius of the SBP derivative is introduced, and the effectiveness of this method is explored through numerical experiment.

Recommendations

• Optimization of multidimensional diagonal-norm summation-by-parts operators on simplices
• Summation-by-parts operators and high-order quadrature
• Corner-corrected diagonal-norm summation-by-parts operators for the first derivative with increased order of accuracy
• On the order of Accuracy of Finite Difference Operators on Diagonal Norm Based Summation-by-Parts Form
• Review of summation-by-parts operators with simultaneous approximation terms for the numerical solution of partial differential equations

Cites work

• scientific article; zbMATH DE number 3553701 (Why is no real title available?)
• scientific article; zbMATH DE number 2107836 (Why is no real title available?)
• scientific article; zbMATH DE number 6733711 (Why is no real title available?)
• A fourth order accurate finite difference scheme for the elastic wave equation in second order formulation
• A generalized framework for nodal first derivative summation-by-parts operators
• A solution to the stability issues with block norm summation by parts operators
• A stable high-order finite difference scheme for the compressible Navier-Stokes equations, far-field boundary conditions
• A stable high-order finite difference scheme for the compressible Navier-Stokes equations: No-slip wall boundary conditions
• On coordinate transformations for summation-by-parts operators
• Optimal diagonal-norm SBP operators
• Optimized high-order derivative and dissipation operators satisfying summation by parts, and applications in three-dimensional multi-block evolutions
• Review of summation-by-parts schemes for initial-boundary-value problems
• Stable Difference Approximations for the Elastic Wave Equation in Second Order Formulation
• Summation by parts for finite difference approximations for \(d/dx\)
• Summation-by-parts operators and high-order quadrature
• Time-stable boundary conditions for finite-difference schemes solving hyperbolic systems: Methodology and application to high-order compact schemes

Cited in

(7)

• Boundary optimized diagonal-norm SBP operators
• Optimization of multidimensional diagonal-norm summation-by-parts operators on simplices
• High-fidelity sound propagation in a varying 3D atmosphere
• Summation-by-parts operators for general function spaces: the second derivative
• Compatible diagonal-norm staggered and upwind SBP operators
• Simulation of acoustic and flexural-gravity waves in ice-covered oceans
• Corner-corrected diagonal-norm summation-by-parts operators for the first derivative with increased order of accuracy

Describes a project that uses

Uses Software

• SymPy
• CVXOPT

This page was built for publication: An algorithmic exploration of the existence of high-order summation by parts operators with diagonal norm

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