Construction of hybrid numerical flows that suppress the development of shock-wave instability

In this paper, we propose a new hybrid numerical flux for calculating the fluxes of the Euler part of the Navier–Stokes system of equations. This flux avoids the occurrence of instability and maintains high accuracy on shock waves and boundary layers. This flux is a combination of Godunov's numerical flow and Rusanov–Lax–Friedrichs numerical flow. Numerical simulation of supersonic flow around a cruise missile Tomahawk has been carried out.