Lee Model is widely used to model mass transfer in two-phase flow. The phase change model involves a time relaxation coefficient, which has a significant impact on numerical simulation. Two benchmark problems related to two-phase flow have been numerically simulated using the volume of fluid (VOF) approach with a customized solver in OpenFOAM 9. The influence of the time relaxation coefficient for 1-D Stefan and 2-D film boiling problems is studied, and the same is discussed here. The results indicate that for different fluids, the time relaxation parameter yields accurate results within certain ranges; however, outside these ranges, the precision of the numerical simulation diminishes significantly. While analyzing various fluids, it was found that good numerical simulations are achieved with a time relaxation coefficient ranging between \(10^3\) and \(10^7\) . For this range of \(r_e\) , the error is less than \(1.5\%\) for the Stefan problem and within \(15\%\) with correlations for 2-D film boiling. The time relaxation coefficient is influenced by the phase initiating the transition, which shows variations in bubble growth during the simulation of 2-D film boiling with different values of \(r_e\) . This variation is also mirrored in the periodicity of the average Nusselt number.