Very simple, carbuncle-free, boundary-layer-resolving, rotated-hybrid Riemann solvers

In this paper, we propose new Euler flux functions for use in a finite-volume Euler/Navier-Stokes code, which are very simple, carbuncle-free, yet have an excellent boundary-layer-resolving capability, by combining two different Riemann solvers into one based on a rotated Riemann solver approach. We show that very economical Euler flux functions can be devised by combining the Roe solver (a full-wave solver) and the Rusanov/HLL solver (a fewer-wave solver), based on a rotated Riemann solver approach: a fewer-wave solver automatically applied in the direction normal to shocks to suppress carbuncles and a full-wave solver applied, again automatically, across shear layers to avoid an excessive amount of dissipation. The resulting flux functions can be implemented in a very simple and economical manner, in the form of the Roe solver with modified wave speeds, so that converting an existing Roe flux code into the new fluxes is an extremely simple task. They require only 7-14% extra CPU time and no problem-dependent tuning parameters. These new rotated fluxes are not only robust for shock-capturing, but also accurate for resolving shear layers. This is demonstrated by an extensive series of numerical experiments with standard finite-volume Euler and Navier-Stokes codes, including various shock instability problems and also an unstructured grid case.

[1]  P. Roe,et al.  On Godunov-type methods near low densities , 1991 .

[2]  Philip L. Roe,et al.  A Proposed Cure to the Carbuncle Phenomenon , 2009 .

[3]  Bram van Leer,et al.  Upwind-difference methods for aerodynamic problems governed by the Euler equations , 1985 .

[4]  V. Guinot.  Approximate Riemann Solvers , 2010 .

[5]  P. Janhunen,et al.  A Positive Conservative Method for Magnetohydrodynamics Based on HLL and Roe Methods , 2000 .

[6]  Jean-Marc Moschetta,et al.  Shock wave numerical structure and the carbuncle phenomenon , 2005 .

[7]  P. Lax,et al.  On Upstream Differencing and Godunov-Type Schemes for Hyperbolic Conservation Laws , 1983 .

[8]  S. Imlay,et al.  Blunt-body flow simulations , 1988 .

[9]  P. Roe.  Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

[10]  G. D. van Albada,et al.  A comparative study of computational methods in cosmic gas dynamics , 1982 .

[11]  J. Quirk.  A Contribution to the Great Riemann Solver Debate , 1994 .

[12]  Keiichi Kitamura,et al.  An Evaluation of Euler Fluxes for Hypersonic Flow Computations , 2007 .

[13]  Chongam Kim,et al.  Cures for the shock instability: development of a shock-stable Roe scheme , 2003 .

[14]  E. Toro,et al.  Restoration of the contact surface in the HLL-Riemann solver , 1994 .

[15]  R. LeVeque.  Approximate Riemann Solvers , 1992 .

[16]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[17]  S. F. Davis,et al.  A rotationally biased upwind difference scheme for the euler equations , 1984 .

[18]  C. Angelopoulos.  High resolution schemes for hyperbolic conservation laws , 1992 .

[19]  Richard Sanders,et al.  Regular ArticleMultidimensional Dissipation for Upwind Schemes: Stability and Applications to Gas Dynamics☆ , 1998 .

[20]  Domenic D'Ambrosio,et al.  Numerical Instablilities in Upwind Methods: Analysis and Cures for the “Carbuncle” Phenomenon , 2001 .

[21]  P. Woodward,et al.  The numerical simulation of two-dimensional fluid flow with strong shocks , 1984 .

[22]  Bernd Einfeld.  On Godunov-type methods for gas dynamics , 1988 .

[23]  Meng-Sing Liou,et al.  Mass Flux Schemes and Connection to Shock Instability , 2000 .

[24]  Meng-Sing Liou,et al.  A sequel to AUSM, Part II: AUSM+-up for all speeds , 2006, J. Comput. Phys..

[25]  Bram van Leer,et al.  Use of a rotated Riemann solver for the two-dimensional Euler equations , 1993 .

[26]  Philip L. Roe,et al.  On carbuncles and other excrescences , 2005 .

[27]  Chongam Kim,et al.  Methods for the accurate computations of hypersonic flows: II. shock-aligned grid technique , 2001 .

[28]  M. Liou.  A Sequel to AUSM , 1996 .

[29]  Richard Sanders,et al.  Multidimensional Dissipation for Upwind Schemes , 1998 .

[30]  Farzad Ismail,et al.  Toward a reliable prediction of shocks in hypersonic flow: Resolving carbuncles with entropy and vorticity control , 2006 .

[31]  Yu-Xin Ren,et al.  A robust shock-capturing scheme based on rotated Riemann solvers , 2003 .

[32]  Kun Xu,et al.  Dissipative mechanism in Godunov‐type schemes , 2001 .

[33]  Michael Dumbser,et al.  A matrix stability analysis of the carbuncle phenomenon , 2004 .

[34]  V. Rusanov,et al.  The calculation of the interaction of non-stationary shock waves and obstacles , 1962 .