Robust HLLC Riemann solver with weighted average flux scheme for strong shock

Many researchers have reported failures of the approximate Riemann solvers in the presence of strong shock. This is believed to be due to perturbation transfer in the transverse direction of shock waves. We propose a simple and clear method to prevent such problems for the Harten-Lax-van Leer contact (HLLC) scheme. By defining a sensing function in the transverse direction of strong shock, the HLLC flux is switched to the Harten-Lax-van Leer (HLL) flux in that direction locally, and the magnitude of the additional dissipation is automatically determined using the HLL scheme. We combine the HLLC and HLL schemes in a single framework using a switching function. High-order accuracy is achieved using a weighted average flux (WAF) scheme, and a method for v-shear treatment is presented. The modified HLLC scheme is named HLLC-HLL. It is tested against a steady normal shock instability problem and Quirk's test problems, and spurious solutions in the strong shock regions are successfully controlled.

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

[2]  E. Toro.  Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

[3]  W. Rider,et al.  High-Resolution Methods for Incompressible and Low-Speed Flows , 2004 .

[4]  Eleuterio F. Toro,et al.  The weighted average flux method applied to the Euler equations , 1992, Philosophical Transactions of the Royal Society of London. Series A: Physical and Engineering Sciences.

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

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

[7]  Eleuterio F. Toro,et al.  Unsplit WAF-Type Schemes for Three Dimensional Hyperbolic Conservation Laws , 1998 .

[8]  Eleuterio F. Toro,et al.  Numerical Methods for Wave Propagation , 2011 .

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

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

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

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

[13]  In-Seuck Jeung,et al.  Realization of contact resolving approximate Riemann solvers for strong shock and expansion flows , 2009 .

[14]  J. J. Quirk,et al.  An adaptive grid algorithm for computational shock hydrodynamics , 1991 .

[15]  Derek M. Causon,et al.  On the Choice of Wavespeeds for the HLLC Riemann Solver , 1997, SIAM J. Sci. Comput..

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

[17]  Eleuterio F. Toro,et al.  AOn WAF-Type Schemes for Multidimensional Hyperbolic Conservation Laws , 1997 .

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

[19]  Keiichi Kitamura,et al.  Evaluation of Euler Fluxes for Hypersonic Flow Computations , 2009 .

[20]  Eleuterio F. Toro,et al.  A weighted average flux method for hyperbolic conservation laws , 1989, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.