The computational model that is able to estimate the temperature distribution inside a solid specimen during the film boiling phase of immersion quenching (water entry) process has been presented in paper, and it is based on the prescribed initial temperatures of the solid specimen and the liquid quenchant. In addition, the turbulence effects have to be considered using the assumed turbulence kinetic energy value, i.e., the “frozen turbulence” approach, that remains constant thorough the simulation. The studied material is nickel alloy, Inconel 600, for which the extensive experimental data is available. The work has been carried out using ANSYS Fluent computational fluid dynamics software and the methods for solution of Stefan problem by Eulerian two fluid VOF model.. A satisfactory agreement between the experimental and the calculated data has been achieved, yielding thereby the computationally obtained data that fits to a great extent to the prescribed error band of ± 10 % during the estimated film boiling phase of the immersion quenching process itself.