We use central differences to solve the time dependent Euler equations. The schemes are all advanced using a Runge-Kutta formula in time. Near shocks, a second difference is added as an artificial viscosity. This reduces the scheme to a first order upwind scheme at shocks. The switch that is used guarantees that the scheme is locally total variation diminishing (TVD). For steady state problems it is usually advantageous to relax this condition. Then small oscillations do not activate the switches and the convergence to a steady state is improved. To sharpen the shocks, different coefficients are needed for different equations and so a matrix valued dissipation is introduced and compared with the scalar viscosity. The connection between this artificial viscosity and flux limiters is shown. Any flux limiter can be used as the basis of a shock detector for an artificial viscosity. We compare the use of the van Leer, van Albada, mimmod, superbee, and the 'average' flux limiters for this central difference scheme. For time dependent problems, we need to use a small enough time step so that the CFL was less than one even though the scheme was linearly stable for larger time steps. Using a total variation bounded (TVB) Runge-Kutta scheme yields minor improvements in the accuracy.
[1]
A. Jameson,et al.
Numerical solution of the Euler equations by finite volume methods using Runge Kutta time stepping schemes
,
1981
.
[2]
H. Huynh,et al.
Accurate monotone cubic interpolation
,
1993
.
[3]
Eli Turkel,et al.
On Central-Difference and Upwind Schemes
,
1992
.
[4]
P. Sweby.
High Resolution Schemes Using Flux Limiters for Hyperbolic Conservation Laws
,
1984
.
[5]
Alain Lerat,et al.
Efficient solution of the steady Euler equations with a centered implicit method
,
1988
.
[6]
S. Osher,et al.
Efficient implementation of essentially non-oscillatory shock-capturing schemes,II
,
1989
.
[7]
Chi-Wang Shu.
Total-variation-diminishing time discretizations
,
1988
.