The aim of the present study is to propose a convenient numerical calculation method combined with finite volume and temperature recovery methods, and to confirm the validity of the present numerical calculation method through a comparison with experiments and numerical calculation conducted by other researchers. In this paper, the details of the combined numerical calculation method have been first described. Secondary, the numerical results computed by the present proposed method has been compared with experimental and numerical calculation results obtained by other researchers in order to verify the effectiveness of the present combined numerical calculation method. As a result, present numerical calculation results for both the shape of ice formation and streamlines in unfrozen water were in good agreement with the experimental and numerical calculation ones previously obtained by other researchers.