Metallic cylindrical materials are undertaken to different heat treatment processes in industry. The effects of velocity and subcooling during the investigation of quenching of steel with water are of great interest. In addition, variation of surface and center temperatures of steel and water temperature are common research efforts. In this study, temperature profiles in cylindrical ANSI 1025 steel bar under quenching boundary conditions are investigated theoretically using numerical methods to simulate heat transfer process in iron and steel production. While making the calculations for cooling process, variable physical properties of ANSI 1025 steel is taken into consideration. The results indicate that the software written in FORTRAN computer programming language can be used to obtain the cooling curves during the quenching process.