We have developed a cylindrical mixing layer model of a stellar jet including cooling effect in order to understand an optical emission mechanism along collimated high velocity stellar jets associated with young stellar objects. The cylindrical results have been calculated to be the same as the 2D ones presented by Canto & Raga(1991) because the entrainment efficiency in our cylindrical model has been obtained to be the same value as the 2D model has given. We have discussed the morphological and physical characteristics of the mixing layers by the cooling effect. As the jet Mach number increases, the initial temperature of the mixing layer goes high because the kinetic energy of the jet partly converts to the thermal energy of the mixing layer. The initial cooling of the mixing layer is very severe, changing its outer boundary radius. A subsequent change becomes adiabatic. The number of the Mach disks in the stellar jet and the total radiative luminosity of the mixing layer, based on our cylindrical calculation, have quite agreed with the observations.