Browse Prior Art Database

ON [Equation ommitted] -INSTABILITY AND OSCTLtATION AT DISCONTINUITIES IN FINITE DIFFERENCE SCHEMES

IP.com Disclosure Number: IPCOM000128195D
Original Publication Date: 1984-Dec-31
Included in the Prior Art Database: 2005-Sep-15
Document File: 5 page(s) / 23K

Publishing Venue

Software Patent Institute

Related People

Lloyd N. Trefethen: AUTHOR [+3]

Abstract

It is known that even-order finite differ-ence models of hyperbolic partial differential equations generate spurious oscillations ar-ound discontinuities, and also that they are unstable in RP norms for p#2. This paper pre-sents an elementary argument involving disper-sion and dissipation which shows that these two phenomena are closely related, explains their physical basis, and reproduces the known estimates for the width of the region of os-cillations and for the strength of the insta-bility.

This text was extracted from a PDF file.
This is the abbreviated version, containing approximately 23% of the total text.

Page 1 of 5

THIS DOCUMENT IS AN APPROXIMATE REPRESENTATION OF THE ORIGINAL.

ON

(Equation Omitted)

-INSTABILITY AND OSCTLtATION AT DISCONTINUITIES IN FINITE DIFFERENCE SCHEMES

B y

Lloyd N. Trefethen*

April 1984

Report # 116

9' *Supported by an NSF Postdoctoral Fellowship and by the U.S. Department of Energy under contract DE-AG02-76-ER03077-V. ON LP-INSTABILITY AND OSCILLATION AT DISCONTINUITIES IN FINITE DIFFERENCE SCHEMES

Lloyd N. Trefethen Courant Institute of Mathematical Sciences New York University , New York, NY 10012

Abstract

It is known that even-order finite differ-ence models of hyperbolic partial differential equations generate spurious oscillations ar-ound discontinuities, and also that they are unstable in RP norms for p#2. This paper pre-sents an elementary argument involving disper-sion and dissipation which shows that these two phenomena are closely related, explains their physical basis, and reproduces the known estimates for the width of the region of os-cillations and for the strength of the insta-bility.

1. Introduction It is well known that finite difference formulas for hyperbolic partial differential equations often suffer from spurious oscilla-tions near discontinuities. As a model hyper-bolic equation it is customary to consider the simple linear first-order wave equation

(Equation Omitted)

Various precise results for the behavior of the region of oscillations generated by finite difference approximations to (1) have been ob-tained over the years by Apelkrans, Brenner, Chin, Hedstrom, Serdjukova, Thomee, and others [1,2,4,6,7]. In particular, Chin and Hedstrom
[4] have shown by saddle-point analysis that the numerical solutions to (1) behave approxi- mately like integrals of generalized Airy functions.

Another widely recognized fact, first proved by Thomee (1964, unpublished) and also mentioned on p. 100 of the book by Richtmyer and Morton (8], is that every finite differ-ence model of (I1 with even order of accuracy is unstable in the RP norm for every pe[1,-] with p#2. See also [5]. For example the leap frog, Lax-Wendroff, and Crank-Nicolson differ-ence formulas are all 0- unstable. The insta-bility is weak, but it may have undesirable consequences in extensions to

New York University Page 1 Dec 31, 1984

Page 2 of 5

ON [Equation ommitted] -INSTABILITY AND OSCTLtATION AT DISCONTINUITIES IN FINITE DIFFERENCE SCHEMES

nonlinear prob-lems, where both Rl and R have a natural sig-nificance. For this reason a number of mathe-maticians have studied 0-instability of fi-nite difference formulas during the past twen-ty-five years, including Brenner, Hedstrom, Serdjukova, Stetter, Strang, Thomee, and Wahl-bin. A wealth of results of this work are presented in the monograph [3]. Most of the proofs given there are based on techniques of Fourier multipliers, the specialization to constant coefficients of pseudodifferential operators. The purpose of this brief paper is to show that both Rp-instability and oscillations at discontinui...