We adopted the finite volume method to solve the governing equations, while the SIMPLE [30] algorithm was employed to couple the velocity field and pressure. A structured grid was generated for the supersonic separator while a finer grid scheme in the boundary layer is employed in Laval nozzle and supersonic channel. The grid independence was tested before we carried out