Implicit multigrid Euler solutions with symmetric Total-Variation-Diminishing dissipation

A symmetric Total-Variation-Diminishing (TVD) formulation of the numerical dissipation terms has been incorporated into a diagonalized alternating direction implicit multigrid algorithm to solve the Euler equations of inviscid compressible flow. The new treatment of the dissipation makes is possible to capture both very strong and very weak shocks, virtually without oscillation for the steady flows of interest here. In addition, the TVD constraint fixes one of the two previously arbitrary constants in the formulation of the dissipation, and results in both converged solutions and convergence rates which are relatively insensitive to the choice of the remaining dissipation parameter.