Supersonic turbine is widely used in the turbo pump of modern rocket. A preliminary design method for supersonic turbine has been developed considering the coupling effects of turbine and nozzle. Numerical simulation has been proceeded to validate the feasibility of the design method. As the strong shockwave reflected on the mixing plane, additional numerical simulated error would be produced by the mixing plane model in the steady CFD. So unsteady CFD is employed to investigate the aerodynamic performance of the turbine and flow field in passage. Results showed that the preliminary design method developed in this paper is suitable for designing supersonic turbine. This periodical variation of complex shockwave system influences the development of secondary flow, wake and shock-boundary layer interaction, which obviously affect the secondary loss in vane passage. The periodical variation also influences the strength of reflecting shockwave, which affects the profile loss in vane passage. Besides, high circumferential velocity at vane outlet and short blade lead to high radial pressure gradient, which makes the low kinetic energy fluid moves towards hub region and produces additional loss.