Cures for the shock instability: Development of a shock-stable Roe scheme.

*(English)*Zbl 1062.76538Summary: This paper deals with the development of an improved Roe scheme that is free from the shock instability and still preserves the accuracy and efficiency of the original Roe’s Flux Difference Splitting (FDS). Roe’s FDS is known to possess good accuracy but to suffer from the shock instability, such as the carbuncle phenomenon. As the first step towards a shock-stable scheme, Roe’s FDS is compared with the HLLE scheme to identify the source of the shock instability. Through a linear perturbation analysis on the odd–even decoupling problem, damping characteristic is examined and Mach number-based functions \(f\) and \(g\) are introduced to balance damping and feeding rates, which leads to a shock-stable Roe scheme. In order to satisfy the conservation of total enthalpy, which is crucial in predicting surface heat transfer rate in high-speed steady flows, an analysis of dissipation mechanism in the energy equation is carried out to find out the error source and to make the proposed scheme preserve total enthalpy. By modifying the maximum-minimum wave speed, the problem of expansion shock and numerical instability in the expansion region is also remedied without sacrificing the exact capturing of contact discontinuity. Various numerical tests concerned with the shock instability are performed to validate the robustness of the proposed scheme. Then, viscous flow test cases ranging from transonic to hypersonic regime are calculated to demonstrate the accuracy, robustness, and other essential features of the proposed scheme.

##### MSC:

76M20 | Finite difference methods applied to problems in fluid mechanics |

76L05 | Shock waves and blast waves in fluid mechanics |

65M06 | Finite difference methods for initial value and initial-boundary value problems involving PDEs |

##### Keywords:

Shock instability; Contact discontinuity; Expansion shock; Total enthalpy; preservation; Roe’s FDS; HLLE; RoeM scheme##### Software:

AUSMPW+
PDF
BibTeX
XML
Cite

\textit{S.-s. Kim} et al., J. Comput. Phys. 185, No. 2, 342--374 (2003; Zbl 1062.76538)

Full Text:
DOI

##### References:

[1] | Godunov, S.K., A difference scheme for the numerical computation of discontinuous solutions of hydrodynamic equations, Math. sbornik, 47, 271, (1959) · Zbl 0171.46204 |

[2] | Roe, P.L., Approximate Riemann solvers, parameter vectors, and difference schemes, J. comput. phys., 43, 357, (1981) · Zbl 0474.65066 |

[3] | Einfeldt, B.; Munz, C.D.; Roe, P.L.; Sjögreen, B., On Godunov-type methods near low densities, J. comput. phys., 92, 273, (1991) · Zbl 0709.76102 |

[4] | Osher, S.; Chakravarthy, S., Upwind scheme and boundary conditions with applications to Euler equations in generalized geometries, J. comput. phys., 50, 447, (1983) · Zbl 0518.76060 |

[5] | K.M. Peery, S.T. Imlay, Blunt-body flow simulations, AIAA Paper 88-2904, 1988 |

[6] | Quirk, J.J., A contribution to the great Riemann solver debate, Int. J. numer. meth. fluids, 18, 555, (1994) · Zbl 0794.76061 |

[7] | Wada, Y.; Liou, M.S., An accurate and robust flux splitting scheme for shock and contact discontinuities, SIAM J. sci. comput., 18, 3, 633, (1997) · Zbl 0879.76064 |

[8] | Lin, H.C., Dissipation addition to flux-difference splitting, J. comput. phys., 117, 20, (1995) · Zbl 0836.76061 |

[9] | Liou, M.S., Mass flux schemes and connection to shock instability, J. comput. phys., 160, 623, (2000) · Zbl 0967.76062 |

[10] | K. Xu, Gas-kinetic schemes for unsteady compressible flow simulations, VKI LS 1998-03 (1998) 1 |

[11] | Sanders, R.; Morano, E.; Druguet, M.C., Multidimensional dissipation for upwind schemes: stability and applications to gas dynamics, J. comput. phys., 145, 511, (1998) · Zbl 0924.76076 |

[12] | Gressier, J.; Moschetta, J.M., Robustness versus accuracy in shock-wave computations, Int. J. numer. meth. fluids., 33, 313, (2000) · Zbl 0980.76072 |

[13] | Robinet, J.C.; Gressier, J.; Casalis, G.; Moschetta, J.M., Shock wave instability and the carbuncle phenomenon: same intrinsic origin?, J. fluid mech., 417, 237, (2000) · Zbl 0961.76030 |

[14] | Pandolfi, M.; D’Ambrosio, D., Numerical instabilities in upwind methods: analysis and cures for the “carbuncle” phenomenon, J. comput. phys., 166, 271, (2001) · Zbl 0990.76051 |

[15] | Kim, K.H.; Rho, O.H., An improvement of AUSM schemes by introducing the pressure-based weight functions, Computers & fluids, 27, 3, 311, (1998) · Zbl 0964.76064 |

[16] | Kim, K.H.; Kim, C.; Rho, O.H., Methods for the accurate computations f hypersonic flows I. AUSMPW+scheme, J. comput. phys., 174, 38, (2001) · Zbl 1106.76421 |

[17] | K.H. Kim, C. Kim, Uncertainty in the Euler equations and numerical instability, prepared for publication 2002 |

[18] | Jameson, A., Analysis and design of numerical schemes for gas dynamics 2. artificial diffusion and discrete shock structure, Int. J. comput. fluid dynamics, 5, 1, (1995) |

[19] | D. Hänel, On the accuracy of upwind schemes in the solutions of the Navier-Stokes equations, AIAA Paper 87-1105-CP, 1987 |

[20] | Woodward, P.; Collela, P., The numerical simulation of two-dimensional fluid flow with strong shocks, J. comput. phys., 54, 115, (1984) |

[21] | R.J. Hakkinen, I. Greber, L. Trilling, S.S. Abarnel, The interaction of an oblique shock wave with a laminar boundary layer, Tr 57-1, Fluid Dynamics Research Group, MIT, May 1957 · Zbl 0079.40806 |

[22] | P.H. Cook, M.A. McDonald, M.C.P. Firmin, AIRFOIL RAE 2822 pressure distributions, boundary layer and wake measurements, AGARD Advisory Report 138, 1979 |

[23] | Menter, F.R., Two-equation eddy-viscosity turbulence models for engineering applications, Aiaa j., 32, 1598, (1994) |

[24] | M.S. Holden, A.R. Weiting, J.R. Moselle, C. Glass, Studies of aerothermal loads generated in regions of shock/shock interaction in hypersonic flow, AIAA Paper 88-0047, Jan 1988 |

[25] | V. Schmitt, F. Charpin, Pressure distribution on the ONERA M6 Wing at Transonic Mach Numbers, AGARD-AR-138, 1979 |

This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.