In present study, volume of fluid method in OpenFOAM open source CFD package will be extended to consider phase change phenomena due to condensation process. Both phases (liquid – vapor) are incompressible and immiscible. Vapor phase is assumed in saturated temperature. Interface between two phases are tracked with color function volume of fluid (CF-VOF) method. ُSurface Tension is taken accounted by Continuous Surface Force (CSF) model and mass transfer occurs along interface is considered by Lee mass transfer model. Pressure-Velocity coupling will be solved with PISO algorithm in the collocated grid. This solver is validated with Stefan problem. In one dimensional Stefan problem, the desistance of interface motion from cold wall is compared by the analytical solution. Then condensate laminar liquid film flow over vertical plate is simulated in the presence of gravity. Numerical result shows calculated film thickness from numerical simulation is thinner than analytical solution. Also, it shows Nusselt number is a function of vapor specific heat which neglected in existing correlations, therefore analytical solution and experimental correlation should modified to consider this effect on the Nusselt Number.