The Volume of Fluid method(VOF) integrated in fluid dynamic program ANSYS FLUENT was adopted to calculate the process of single bubble growth in subcooled static water under two-dimensional condition. The phase change process could be simulated by adding extra mass source term and energy source term to the Navier-Stokes equations in the cells adjacent to vapor-liquid interface. The surface tension force, contact angle between bubble and wall and the evaporation of micro layer beneath the bubble had also been taken into consideration. The reliability of the calculation could be verified by comparing computational results with contrast experimental results. The numerical analysis could also obtain some physical parameters which are hard to be got from experiment directly, such as temperature fields and velocity fields.