A discussion on numerical shock stability of unstructured finite volume method: Riemann solvers and limiters

Numerical shock instability is a complexity which may occur in supersonic simulations. Riemann solver is usually the crucial factor that affects both the computation accuracy and numerical shock stability. In this paper, several classical Riemann solvers are discussed, and the intrinsic mechanism of shock instability is especially concerned. It can be found that the momentum perturbation traversing shock wave is a major reason that invokes instability. Furthermore, slope limiters used to depress oscillation across shock wave is also a key factor for computation stability. Several slope limiters can cause significant numerical errors near shock waves, and make the computation fail to converge. Extra dissipation of Riemann solvers and slope limiters can be helpful to eliminate instability, but reduces the computation accuracy. Therefore, to properly introduce numerical dissipation is critical for numerical computations. Here, pressure based shock indicator is used to show the position of shock wave and tunes the numerical dissipation. Overall, the presented methods are showing satisfactory results in both the accuracy and stability.