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 this paper. 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 extensive experimental data are 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 fit to a great extent the prescribed error band of ±10% during the estimated film boiling phase of the immersion quenching process itself. It was, however, found that the temperature calculated in the center of a specimen fits this error band until reaching somewhere t < 6 s due to low presumed turbulence level in the domain. In addition, the explosion of the vapor phase after the body reaches the free surface of the quenchant has also been successfully tracked using the numerical simulation model proposed herein. The major novelty of the present research lies in the fact that a moving boundary problem has been successfully resolved in conjunction with, to a great extent, basic-principle-based heat and mass transfer in a turbulent flow conjugate heat transfer (CHT) numerical simulation using moderate computational resources.