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

DOI10.1007/S10915-016-0211-ZzbMATH Open1370.65039arXiv1403.5750OpenAlexW3102712297MaRDI QIDQ2014030FDOQ2014030

Joshua Klarmann, Nathan Albin

Publication date: 10 August 2017

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

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.

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)

Cites Work

• Title not available (Why is that?)
• Title not available (Why is that?)
• Title not available (Why is that?)
• Time-stable boundary conditions for finite-difference schemes solving hyperbolic systems: Methodology and application to high-order compact schemes
• Summation by parts for finite difference approximations for \(d/dx\)
• A stable high-order finite difference scheme for the compressible Navier-Stokes equations, far-field boundary conditions
• 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
• A generalized framework for nodal first derivative summation-by-parts operators
• 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
• A solution to the stability issues with block norm summation by parts operators
• A fourth order accurate finite difference scheme for the elastic wave equation in second order formulation
• Optimal diagonal-norm SBP operators
• Summation-by-parts operators and high-order quadrature
• Stable Difference Approximations for the Elastic Wave Equation in Second Order Formulation

Cited In (5)

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

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)