Riemann solver which is capable of solving contact discontinuities accurately but with
robustness for strong shock. It is well known that the original HLL, which has incomplete
wave structures, is too dissipative to capture contact discontinuities accurately. On the other
side, contact-capturing approximate Riemann solvers such as HLL with Contact (HLLC)
usually suffer from spurious solutions, also called carbuncle phenomenon, for strong shock …