A numerical analysis program is developed by FDM scheme for the prediction of microstructural transformation during heat treatment of steels. In this study, multi-phase model was used fur description of diffusional austenite transformations in low-alloy hypoeutectoid steels during cooling after austenitization. A fundamental property of the model consisting of coupled differential equations is that by taking into account the rate of austenite grain growth, it permits the prediction of the progress of ferrite, pearlite, and bainite transformations simultaneously during quenching and estimate the amount of martensite also by using K-M eq. In order to simulate the microstructural evolution during tempering process, another Avrami-type eq. was adopted and method for vickers hardness prediction was also proposed. To verify the developed program, the calculated results are compared with experimental ones of casting product. Based on these results, newly designed heat treatment process is proposed and it was proved to be effective for industry.