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.