In this study, a two-phase flow investigation on falling film evaporation phenomenon is described. To track liquid-vapor interface, the volume of fluid multiphase model based on the piecewise linear interface construction method in curvilinear coordinate is used. A finite volume method code based on the SIMPLE algorithm and non-orthogonal discretization grid system is used for solving the governing equations including continuity, energy, and Reynolds averaged Navier-Stokes equations with the k-ε turbulence model. Also, the energy and mass transfer during the phase change is taken into account. The effects of inlet mass flow rate and temperature difference between the tube temperature and the saturation temperature on the local and average heat transfer coefficient and the net vapor production are presented. Results show the average heat transfer coefficient increases/decreases by increasing in the inlet mass flow rate/the temperature difference at a constant temperature difference/inlet mass flow rate. Results also demonstrate at a constant value of DT= _ m f (MED-TVC design approach), as the inlet mass flow rate increases, the heat transfer coefficient may be decreased or increased. This increase or decrease depends on the values of DT= _ m f and ΔT.