On the stability of the unsmoothedFourier method for hyperbolic equations

Jonathan Goodman, Thomas Hou, Eitan Tadmor

Research output: Contribution to journalArticle

Abstract

It has been a long open question whether the pseudospectral Fourier methodwithout smoothing is stable for hyperbolic equations with variablecoefficients that change signs. In this work we answer this question with adetailed stability analysis of prototype cases of the Fourier method.We show that due to weighted (Formula presented.)-stability,the (Formula presented.)-degree Fourier solutionis algebraically stable in the sense that its (Formula presented.)amplification does not exceed (Formula presented.).Yet, the Fourier method is weakly(Formula presented.)- unstablein the sense that it does experience such (Formula presented.)amplification. The exact mechanism of thisweak instability is due the aliasing phenomenon, which isresponsible for an (Formula presented.) amplification of the Fourier modes atthe boundaries of the computed spectrum.Two practical conclusions emerge from our discussion. First,the Fourier method is required to have sufficiently many modes in order toresolve the underlying phenomenon. Otherwise, the lack ofresolution will excite the weak instability which willpropagate from the slowly decaying high modes to the lower ones.Second -- independent of whether smoothing was used or not,the small scale information contained in the highestmodes of the Fourier solution will bedestroyed by their (Formula presented.) amplification. Happily, with enoughresolution nothing worse can happen.

Original languageEnglish (US)
Pages (from-to)93-129
Number of pages37
JournalNumerische Mathematik
Volume67
Issue number1
DOIs
StatePublished - 1994

Fingerprint

Hyperbolic Equations
Amplification
Fourier Method
Smoothing
Aliasing
Sign Change
Stability Analysis
Exceed
Prototype

Keywords

  • Mathematics Subject Classification (1991): 65M12

ASJC Scopus subject areas

  • Applied Mathematics
  • Computational Mathematics

Cite this

On the stability of the unsmoothedFourier method for hyperbolic equations. / Goodman, Jonathan; Hou, Thomas; Tadmor, Eitan.

In: Numerische Mathematik, Vol. 67, No. 1, 1994, p. 93-129.

Research output: Contribution to journalArticle

Goodman, Jonathan ; Hou, Thomas ; Tadmor, Eitan. / On the stability of the unsmoothedFourier method for hyperbolic equations. In: Numerische Mathematik. 1994 ; Vol. 67, No. 1. pp. 93-129.
@article{f070de59b5644fceb6f8ee8097f2ac76,
title = "On the stability of the unsmoothedFourier method for hyperbolic equations",
abstract = "It has been a long open question whether the pseudospectral Fourier methodwithout smoothing is stable for hyperbolic equations with variablecoefficients that change signs. In this work we answer this question with adetailed stability analysis of prototype cases of the Fourier method.We show that due to weighted (Formula presented.)-stability,the (Formula presented.)-degree Fourier solutionis algebraically stable in the sense that its (Formula presented.)amplification does not exceed (Formula presented.).Yet, the Fourier method is weakly(Formula presented.)- unstablein the sense that it does experience such (Formula presented.)amplification. The exact mechanism of thisweak instability is due the aliasing phenomenon, which isresponsible for an (Formula presented.) amplification of the Fourier modes atthe boundaries of the computed spectrum.Two practical conclusions emerge from our discussion. First,the Fourier method is required to have sufficiently many modes in order toresolve the underlying phenomenon. Otherwise, the lack ofresolution will excite the weak instability which willpropagate from the slowly decaying high modes to the lower ones.Second -- independent of whether smoothing was used or not,the small scale information contained in the highestmodes of the Fourier solution will bedestroyed by their (Formula presented.) amplification. Happily, with enoughresolution nothing worse can happen.",
keywords = "Mathematics Subject Classification (1991): 65M12",
author = "Jonathan Goodman and Thomas Hou and Eitan Tadmor",
year = "1994",
doi = "10.1007/s002110050019",
language = "English (US)",
volume = "67",
pages = "93--129",
journal = "Numerische Mathematik",
issn = "0029-599X",
publisher = "Springer New York",
number = "1",

}

TY - JOUR

T1 - On the stability of the unsmoothedFourier method for hyperbolic equations

AU - Goodman, Jonathan

AU - Hou, Thomas

AU - Tadmor, Eitan

PY - 1994

Y1 - 1994

N2 - It has been a long open question whether the pseudospectral Fourier methodwithout smoothing is stable for hyperbolic equations with variablecoefficients that change signs. In this work we answer this question with adetailed stability analysis of prototype cases of the Fourier method.We show that due to weighted (Formula presented.)-stability,the (Formula presented.)-degree Fourier solutionis algebraically stable in the sense that its (Formula presented.)amplification does not exceed (Formula presented.).Yet, the Fourier method is weakly(Formula presented.)- unstablein the sense that it does experience such (Formula presented.)amplification. The exact mechanism of thisweak instability is due the aliasing phenomenon, which isresponsible for an (Formula presented.) amplification of the Fourier modes atthe boundaries of the computed spectrum.Two practical conclusions emerge from our discussion. First,the Fourier method is required to have sufficiently many modes in order toresolve the underlying phenomenon. Otherwise, the lack ofresolution will excite the weak instability which willpropagate from the slowly decaying high modes to the lower ones.Second -- independent of whether smoothing was used or not,the small scale information contained in the highestmodes of the Fourier solution will bedestroyed by their (Formula presented.) amplification. Happily, with enoughresolution nothing worse can happen.

AB - It has been a long open question whether the pseudospectral Fourier methodwithout smoothing is stable for hyperbolic equations with variablecoefficients that change signs. In this work we answer this question with adetailed stability analysis of prototype cases of the Fourier method.We show that due to weighted (Formula presented.)-stability,the (Formula presented.)-degree Fourier solutionis algebraically stable in the sense that its (Formula presented.)amplification does not exceed (Formula presented.).Yet, the Fourier method is weakly(Formula presented.)- unstablein the sense that it does experience such (Formula presented.)amplification. The exact mechanism of thisweak instability is due the aliasing phenomenon, which isresponsible for an (Formula presented.) amplification of the Fourier modes atthe boundaries of the computed spectrum.Two practical conclusions emerge from our discussion. First,the Fourier method is required to have sufficiently many modes in order toresolve the underlying phenomenon. Otherwise, the lack ofresolution will excite the weak instability which willpropagate from the slowly decaying high modes to the lower ones.Second -- independent of whether smoothing was used or not,the small scale information contained in the highestmodes of the Fourier solution will bedestroyed by their (Formula presented.) amplification. Happily, with enoughresolution nothing worse can happen.

KW - Mathematics Subject Classification (1991): 65M12

UR - http://www.scopus.com/inward/record.url?scp=21344494565&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=21344494565&partnerID=8YFLogxK

U2 - 10.1007/s002110050019

DO - 10.1007/s002110050019

M3 - Article

VL - 67

SP - 93

EP - 129

JO - Numerische Mathematik

JF - Numerische Mathematik

SN - 0029-599X

IS - 1

ER -