Computation of Weak Steady Shock Reflections by means of an Unstructured Shock-Fitting Solver

Weak steady Mach reflections are numerically simulated using unstructured grids by means of either “shock-fitting” or “shock-capturing” techniques. It is shown that shock-fitting allows using coarser meshes than those required by shock-capturing, since the latter needs mesh refinement in the direction normal to the discontinuities, which is not needed using the former approach. The shock-fitted solution is also free from the numerical disturbances that arise along the captured discontinuities and pollute the captured solution in the smooth flowfield region. Finally, the shock-fitting solutions show the presence of a small region next to the Mach stem and the reflected shock downstream of the triple point, characterized by very high gradients.