A numerical method for simulation of interfacial flows with phase change is presented. A volume-of-fluid (VOF) method based on piecewise linear interface construction (PLIC) is used to track interface. Spocial treatment is taken to deal with discontinuous velocity field due to phase change at liquid-vapor interface. The treatment is easy and effective. SIMPLE method is used for flow field and a CSF model for surface tension. The method is used to simulate three-dimensional film boiling on a horizontal plate with Cartesian coordinates and axisymmetric two-dimensional natural convection film boiling on a hemisphericallyended rod with body fitted coordinates. They are in good agreement with analytical correlations.