Supersonic coaxial jets are investigated numerically by using the axisymmetric, Wavier-Stokes equations which are solved using a fully implicit finite volume method. Three different kinds of coaxial nozzles are employed to understand the flow physics involved in the supersonic coaxial jets. Two convergent-divergent supersonic nozzles are designed to have the same Mach number 2.0, and used to compare the coaxial jet flows with those discharging from one constant-area nozzle. The impingement angle of the annular jets are varied. The primary pressure ratio is changed in the range from 2.0 to 10.0 and the assistant jet ratio from 1.0 to 3.0. The results obtained show that the fluctuations of the total pressure and Mach number along the jet axis are much higher in the constant-area nozzle than those in the convergent-divergent nozzles, and the constant-area nozzle lead to higher total pressure losses, compared with the convergent-divergent nozzles. The assistant jets from the annular nozzle affect the coaxial jet flows within the distance less than about ten times the nozzle throat diameter, but beyond it the coaxial jet is conical with self-similar velocity profiles. Increasing both the primary jet pressure ratio and the assistant jet pressure ratio produces a longer coaxial jet core.