Designing an efficient solution strategy for fluid flows. 1. A stable high order finite difference scheme and sharp shock resolution for the Euler equations
- Stanford Univ., CA (United States)
We derive high-order finite difference schemes for the compressible Euler (and Navier-Stokes equations) that satisfy a semidiscrete energy estimate and present and efficient strategy for the treatment of discontinuities that leads to sharp shock resolution. The formulation of the semidiscrete energy estimates 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. 29 refs., 7 figs., 2 tabs.
- OSTI ID:
- 530639
- Journal Information:
- Journal of Computational Physics, Vol. 129, Issue 2; Other Information: PBD: Dec 1996
- Country of Publication:
- United States
- Language:
- English
Similar Records
Three-dimensional adaptive grid-embedding Euler technique
The discrete maximum principle for finite element approximations of anisotropic diffusion problems on arbitrary meshes