A Cure for the Sonic Point Glitch

Among the various numerical schemes developed since the ’80s for the computation of the compressible Euler equations, the vast majority produce in certain cases spurious pressure glitches at sonic points. This flaw is particularly visible in the computation of transonic expansions and leads lo unphysical “expansion shocks” when the flow undergoes rapid change of direction.

The analysis of this flaw is presented, based on a series of numerical experiments. For Flux-Vector Splitting methods, it is suggested that it is not the order of differentiability of the numerical flux which is crucial but the way the pressure at an interface is calculated. A new way of evaluating the pressure at the interface is proposed, based upon kinetic theory, and is applied to most current available algorithms including Flux-Vector Splitting and Flux-Difference Splitting methods as well as recent hybrid schemes (AUSM, HUS).