A numerical study was conducted to characterize the thermal behavior of laminar circular liquid jets. The main goal of this work was to clarify the significant recovery effect of large Prandtl number liquid jets reported by the present authors. Radial distribution of local recovery factor was studied in numerical detail. The effects of many factors were investigated including Reynolds number, Prandtl number, nozzle-to-plate spacing, jet velocity, nozzle diameter and velocity profile at nozzle exit. Very large recovery factors over 20 were obtained with high Prandtl number fluid jets at medium velocity. The results exhibited good agreement with the experimental data by the present authors, verifying the great significance of recovery effect with large Prandtl number liquid jets.