A shock-stable numerical scheme accurate for contact discontinuities: applications to 3D compressible flows (Q6144089)

In this paper, the authors studied the mechanism for three-dimensional shock instability by means of linearized stability analysis and the dissipation-controlling approach which introduces the required transverse dissipation in the numerical chock layer and is adopted to enhance the Roe scheme's shock stability. The corresponding innovative healing method base on combination of a novel Roe-type flux and the second-order MUSCL reconstruction is proposed. The deficiency of unphysical expansion shocks is remedied by incorporating the adjacent cell's information in calculation of numerical signal wave speeds. In order to improve the accuracy for contact discontinuities and shear waves, the authors use an algebraic method combining with the Tangent of Hyperbola for Interface Capturing scheme and the Boundary Variation Diminishing algorithm to further minimize density difference in the numerical diffusion term. A series of numerical experiments fully demonstrates an excellent robustness against shock instability and high accuracy for contact discontinuities and shear waves.