Enhanced spectral viscosity approximations for conservation laws

Anne Gelb, Eitan Tadmor

Research output: Contribution to journalArticle

46 Citations (Scopus)

Abstract

In this paper we construct, analyze and implement a new procedure for the spectral approximations of nonlinear conservation laws. It is well known that using spectral methods for nonlinear conservation laws will result in the formation of the Gibbs phenomenon once spontaneous shock discontinuities appear in the solution. These spurious oscillations will in turn lead to loss of resolution and render the standard spectral approximations unstable. The Spectral Viscosity (SV-) method (Tadmor, 1989) was developed to stabilize the spectral method by adding a spectrally small amount of high-frequencies diffusion carried out in the dual space. The resulting SV-approximation is stable without sacrificing spectral accuracy. The SV-method recovers a spectrally accurate approximation to the projection of the entropy solution; the exact projection, however, is at best a first order approximation to the exact solution as a result of the formation of the shock discontinuities. The issue of spectral resolution is addressed by post-processing the SV-solution to remove the spurious oscillations at the discontinuities, as well as increase the first-order - O(1/N) accuracy away from the shock discontinuities. Successful post-processing methods have been developed to eliminate the Gibbs phenomenon and recover spectral accuracy for the SV-approximation. However, such reconstruction methods require a priori knowledge of the locations of the shock discontinuities. Therefore, the detection of these discontinuities is essential to obtain an overall spectrally accurate solution. To this end, we employ the recently constructed enhanced edge detectors based on appropriate concentration factors (Gelb and Tadmor, 1999). Once the edges of these discontinuities are identified, we can utilize a post-processing reconstruction method, and show that the post-processed SV-solution recovers the exact entropy solution with remarkably high-resolution. We apply our new numerical method, the Enhanced SV-method, to two numerical examples, the scalar periodic Burgers' equation and the one-dimensional system of Euler equations of gas dynamics. Both approximations exhibit high accuracy and resolution to the exact entropy solution.

Original languageEnglish (US)
Pages (from-to)3-21
Number of pages19
JournalApplied Numerical Mathematics
Volume33
Issue number1
DOIs
StatePublished - May 2000

Fingerprint

Viscosity Approximation
Spectral Approximation
Conservation Laws
Discontinuity
Conservation
Spectral Methods
Viscosity
Viscosity Method
Entropy Solution
Shock
Post-processing
Gibbs Phenomenon
Entropy
Viscosity Solutions
High Resolution
Approximation
Processing
Projection
Oscillation
First-order

ASJC Scopus subject areas

  • Applied Mathematics
  • Computational Mathematics
  • Modeling and Simulation

Cite this

Enhanced spectral viscosity approximations for conservation laws. / Gelb, Anne; Tadmor, Eitan.

In: Applied Numerical Mathematics, Vol. 33, No. 1, 05.2000, p. 3-21.

Research output: Contribution to journalArticle

Gelb, Anne ; Tadmor, Eitan. / Enhanced spectral viscosity approximations for conservation laws. In: Applied Numerical Mathematics. 2000 ; Vol. 33, No. 1. pp. 3-21.
@article{c753d0f0d03842379ea7799efdcb994d,
title = "Enhanced spectral viscosity approximations for conservation laws",
abstract = "In this paper we construct, analyze and implement a new procedure for the spectral approximations of nonlinear conservation laws. It is well known that using spectral methods for nonlinear conservation laws will result in the formation of the Gibbs phenomenon once spontaneous shock discontinuities appear in the solution. These spurious oscillations will in turn lead to loss of resolution and render the standard spectral approximations unstable. The Spectral Viscosity (SV-) method (Tadmor, 1989) was developed to stabilize the spectral method by adding a spectrally small amount of high-frequencies diffusion carried out in the dual space. The resulting SV-approximation is stable without sacrificing spectral accuracy. The SV-method recovers a spectrally accurate approximation to the projection of the entropy solution; the exact projection, however, is at best a first order approximation to the exact solution as a result of the formation of the shock discontinuities. The issue of spectral resolution is addressed by post-processing the SV-solution to remove the spurious oscillations at the discontinuities, as well as increase the first-order - O(1/N) accuracy away from the shock discontinuities. Successful post-processing methods have been developed to eliminate the Gibbs phenomenon and recover spectral accuracy for the SV-approximation. However, such reconstruction methods require a priori knowledge of the locations of the shock discontinuities. Therefore, the detection of these discontinuities is essential to obtain an overall spectrally accurate solution. To this end, we employ the recently constructed enhanced edge detectors based on appropriate concentration factors (Gelb and Tadmor, 1999). Once the edges of these discontinuities are identified, we can utilize a post-processing reconstruction method, and show that the post-processed SV-solution recovers the exact entropy solution with remarkably high-resolution. We apply our new numerical method, the Enhanced SV-method, to two numerical examples, the scalar periodic Burgers' equation and the one-dimensional system of Euler equations of gas dynamics. Both approximations exhibit high accuracy and resolution to the exact entropy solution.",
author = "Anne Gelb and Eitan Tadmor",
year = "2000",
month = "5",
doi = "10.1016/S0168-9274(99)00067-7",
language = "English (US)",
volume = "33",
pages = "3--21",
journal = "Applied Numerical Mathematics",
issn = "0168-9274",
publisher = "Elsevier",
number = "1",

}

TY - JOUR

T1 - Enhanced spectral viscosity approximations for conservation laws

AU - Gelb, Anne

AU - Tadmor, Eitan

PY - 2000/5

Y1 - 2000/5

N2 - In this paper we construct, analyze and implement a new procedure for the spectral approximations of nonlinear conservation laws. It is well known that using spectral methods for nonlinear conservation laws will result in the formation of the Gibbs phenomenon once spontaneous shock discontinuities appear in the solution. These spurious oscillations will in turn lead to loss of resolution and render the standard spectral approximations unstable. The Spectral Viscosity (SV-) method (Tadmor, 1989) was developed to stabilize the spectral method by adding a spectrally small amount of high-frequencies diffusion carried out in the dual space. The resulting SV-approximation is stable without sacrificing spectral accuracy. The SV-method recovers a spectrally accurate approximation to the projection of the entropy solution; the exact projection, however, is at best a first order approximation to the exact solution as a result of the formation of the shock discontinuities. The issue of spectral resolution is addressed by post-processing the SV-solution to remove the spurious oscillations at the discontinuities, as well as increase the first-order - O(1/N) accuracy away from the shock discontinuities. Successful post-processing methods have been developed to eliminate the Gibbs phenomenon and recover spectral accuracy for the SV-approximation. However, such reconstruction methods require a priori knowledge of the locations of the shock discontinuities. Therefore, the detection of these discontinuities is essential to obtain an overall spectrally accurate solution. To this end, we employ the recently constructed enhanced edge detectors based on appropriate concentration factors (Gelb and Tadmor, 1999). Once the edges of these discontinuities are identified, we can utilize a post-processing reconstruction method, and show that the post-processed SV-solution recovers the exact entropy solution with remarkably high-resolution. We apply our new numerical method, the Enhanced SV-method, to two numerical examples, the scalar periodic Burgers' equation and the one-dimensional system of Euler equations of gas dynamics. Both approximations exhibit high accuracy and resolution to the exact entropy solution.

AB - In this paper we construct, analyze and implement a new procedure for the spectral approximations of nonlinear conservation laws. It is well known that using spectral methods for nonlinear conservation laws will result in the formation of the Gibbs phenomenon once spontaneous shock discontinuities appear in the solution. These spurious oscillations will in turn lead to loss of resolution and render the standard spectral approximations unstable. The Spectral Viscosity (SV-) method (Tadmor, 1989) was developed to stabilize the spectral method by adding a spectrally small amount of high-frequencies diffusion carried out in the dual space. The resulting SV-approximation is stable without sacrificing spectral accuracy. The SV-method recovers a spectrally accurate approximation to the projection of the entropy solution; the exact projection, however, is at best a first order approximation to the exact solution as a result of the formation of the shock discontinuities. The issue of spectral resolution is addressed by post-processing the SV-solution to remove the spurious oscillations at the discontinuities, as well as increase the first-order - O(1/N) accuracy away from the shock discontinuities. Successful post-processing methods have been developed to eliminate the Gibbs phenomenon and recover spectral accuracy for the SV-approximation. However, such reconstruction methods require a priori knowledge of the locations of the shock discontinuities. Therefore, the detection of these discontinuities is essential to obtain an overall spectrally accurate solution. To this end, we employ the recently constructed enhanced edge detectors based on appropriate concentration factors (Gelb and Tadmor, 1999). Once the edges of these discontinuities are identified, we can utilize a post-processing reconstruction method, and show that the post-processed SV-solution recovers the exact entropy solution with remarkably high-resolution. We apply our new numerical method, the Enhanced SV-method, to two numerical examples, the scalar periodic Burgers' equation and the one-dimensional system of Euler equations of gas dynamics. Both approximations exhibit high accuracy and resolution to the exact entropy solution.

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

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

U2 - 10.1016/S0168-9274(99)00067-7

DO - 10.1016/S0168-9274(99)00067-7

M3 - Article

AN - SCOPUS:0034190412

VL - 33

SP - 3

EP - 21

JO - Applied Numerical Mathematics

JF - Applied Numerical Mathematics

SN - 0168-9274

IS - 1

ER -