numerical technique for multiple shock capturing in steady, quasi one-dimensional flows

A numerical technique is given to capture multiple shocks in steady, quasi one-dimensional flows by solving the Euler equations from a sequence of implicit/explicit solutions for the Riemann variables. A supersonic wind tunnel with a variable area diffuser is analyzed with the results compared to exact solutions. Examples are given with both one and two standing shocks. The technique given is an extension of Moretti's scheme for a single discontinuity in a De Laval nozzle. It is shown that this efficient technique is easily adaptable and is equally accurate for multiple discontinuities as it is for a single discontinuity.