Designing an Efficient Solution Strategy for Fluid Flows: 1. A Stable High Order Finite Difference Scheme and Sharp Shock Resolution for the Euler Equations

We derive high-order finite difference schemes for the compressible Euler (and Navier–Stokes equations) that satisfy a semidiscrete energy estimate and present an efficient strategy for the treatment of discontinuities that leads to sharp shock resolution. The formulation of the semidiscrete energy...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Veröffentlicht in:Journal of Computational Physics 1996-12, Vol.129 (2), p.245-262
Hauptverfasser: Gerritsen, Margot, Olsson, Pelle
Format: Artikel
Sprache:eng
Schlagworte:
Online-Zugang:Volltext
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Beschreibung
Zusammenfassung:We derive high-order finite difference schemes for the compressible Euler (and Navier–Stokes equations) that satisfy a semidiscrete energy estimate and present an efficient strategy for the treatment of discontinuities that leads to sharp shock resolution. The formulation of the semidiscrete energy estimate is based on symmetrization of the equations, a canonical splitting of the flux derivative vector, and the use of difference operators that satisfy a discrete analogue to the integration-by-parts procedure used in the continuous energy estimate. For the Euler equations, the symmetrization is designed such as to preserve the homogeneity of the flux vectors. Around discontinuities or sharp gradients, refined grids are created on which the discrete equations are solved after adding artificial viscosity. The positioning of the subgrids and computation of the viscosity are aided by a detection algorithm which is based on a multiscale wavelet analysis of the pressure grid function. The wavelet theory provides easy-to-implement mathematical criteria to detect discontinuities, sharp gradients, and spurious oscillations quickly and efficiently. As the detection algorithm does not depend on the numerical method used, it is of general interest. The numerical method described and the detection algorithm are part of a general solution strategy for fluid flows, which is currently being developed by the authors and collaborators.
ISSN:0021-9991
1090-2716
DOI:10.1006/jcph.1996.0248