Abstract. First order errors downstream of shocks have been detected in computations with higher order shock capturing schemes in one and two dimensions. Based on a matched asymptotic expansion analysis we show how to modify the artificial viscosity and raise the order of accuracy. [20]). The small scale behavior is of significance in the applications above, and the phenomena involved cannot be captured without an appropriate resolution of this fine scale. The efficiency [10] of high order finite difference methods compared to low order methods can, of course, also be used to reduce the computational cost for a given level of accuracy.
KeyHowever, recently it has been reported (see [1], [2], [5]) that solutions of conservation laws obtained by formally high order accurate schemes degenerate to first order downstream of a shock layer. The effect is seen only when (i) the characteristics come out of the shock region and (ii) the solution is nonconstant. Examples in one space dimension where both these conditions are satisfied include steady state calculations for systems with a source term and time dependent calculations for systems with nonconstant initial data.The downstream degeneration of accuracy has also been observed in calculations of acoustic waves in the presence of shocks. This case can be modeled by a scalar, linear equation where the wave speed changes value in a thin region without changing sign (see [2] and references therein). The degeneracy in accuracy is troublesome, even though the first order terms for reasonable mesh-sizes seem to be small in many cases.In [4], an explanation of the degeneracy in accuracy on solutions of conservation laws were given. Steady state solutions to systems with source terms were analyzed by using matched asymptotic expansions for a corresponding viscous equation, the so-called modified equation; see [13]. It was assumed that an inner and an outer solution exist. The inner solution is valid in the shock region and the outer solution elsewhere. The two solutions are matched in a so-called matching zone. The inner solution, the outer solution, and the shock position are expanded in powers of the small viscosity coefficient ε. With ε = O(h) in the shock region, the outer solution contains a term of O(h) downstream of the shock. Upstream this term is zero.