The stability of numerical boundary treatments for compact high-order finite-difference schemes

Mark H. Carpenter, David Gottlieb, Saul Abarbanel

Research output: Contribution to journalArticlepeer-review

355 Scopus citations

Abstract

The stability characteristics of various compact fourth- and sixth-order spatial operators are assessed with the theory of Gustafsson, Kreiss, and Sundström (G-K-S) for the semidiscrete initial boundary value problem. These results are generalized to the fully discrete case with a recently developed theory of Kreiss and Wu. In all cases, favorable comparisons are obtained between G-K-S theory, eigenvalue determination, and numerical simulation. The conventional definition of stability then is sharpened to include only those spatial discretizations that are asymptotically stable (bounded, left half-plane eigenvalues). Many of the higher order schemes that are G-K-S stable are not asymptotically stable. A series of compact fourth- and sixth-order schemes that are both asymptotically and G-K-S stable for the scalar case are then developed.

Original languageEnglish
Pages (from-to)272-295
Number of pages24
JournalJournal of Computational Physics
Volume108
Issue number2
DOIs
StatePublished - 1993
Externally publishedYes

Fingerprint

Dive into the research topics of 'The stability of numerical boundary treatments for compact high-order finite-difference schemes'. Together they form a unique fingerprint.

Cite this