The phase change by the heat transfer accompanies the sharp discontinuity in thermal properties, the complex physical phenomenon and chaotic change in the shape of phase interface. Thus, the numerical simulation of the two-phase flow and heat transfer with phase change has been slow in progress. In this study, the new phase change model is proposed to solve the two-phase flow with the phase change based on the finite volume method with the Volume-of-fluid method. Present phase change model estimate the mass and heat transfer on the reconstructed phase interface using the volume fraction with the information of heat flux from the gas and liquid phase at the phase interface. The proposed model results in a more realistic description of the mass and heat transfer with a discontinuity on the thermal properties at the phase interface. Simulations of one-dimensional Stefan condensation and two-dimensional infinite corner solidification that has the analytic solution demonstrate the application of the present phase change model. The numerical results indicate that the present phase change model accurately predicts the shape of phase interface, the interface displacements and the reasonable temperature distributions compared with several different phase change models. A planar simulation of film boiling and bubble formation is performed in water at near critical pressure on a horizontal surface. And the transient flow condensation of refrigerants in a microchannel was numerically studied. The computational results with the proposed phase change model provided an improved understanding of film boiling and flow condensation yielding quantitative information on unsteady two-phase flow patterns and on the spatially and temporally varying local heat flux along the flow direction.