A novel approach for the solution of Stefan problem within the framework of the multi fluid model supplemented with Volume of Fluid (VOF) method, i.e. two-fluid VOF, is presented in this paper. The governing equation set is comprised of mass, momentum and energy conservation equations, written on a per phase basis and supplemented with closure models via the source terms. In our method, the heat and mass transfer is calculated from the heat transfer coefficient, which has a fictitious function and depends on the local cell size and the thermal conductivity, and the implementation is straightforward because of the usage of the local value instead of a global parameter. The interface sharpness is ensured by the application of the geometrical reconstruction scheme implemented in VOF. The model is verified for three types of computational meshes including triangular cells, and good agreement was obtained for the interface position and the temperature field. Although the developed method was validated only for Stefan problem, the application of the method to engineering problems is considered to be straightforward since it is implemented to a commercial CFD code only using a local value; especially in the field of naval hydrodynamics wherein the reduction of ship resistance using boiling flow can be computed efficiently since the method handles phase change processes using low resolution meshes.