On TVD difference schemes for the three-dimensional Euler equations in general co-ordinates

An improved treatment for the Harten–Yee and Chakravarthy–Osher TVD numerical flux functions in general co-ordinates is presented. The proposed formulation is demonstrated by a series of numerical experiments for three-dimensional flows around the ONERA-M6 wing. The numerical results indicate that it is important to use a suitable artificial compression parameter in order to obtain more accurate solutions around the leading edge of the wing. The two TVD numerical fluxes give excellent results: they capture the shock wave without numerical oscillations, they capture the rapid expansion around the leading edge sharply, they have self-adjusting mechanisms regarding numerical viscosity and they also have robustness.