Numerical analysis is performed for the unsteady axisymmetric two dimensional phase change problem of freezing of water around a vertical tube. Heat conduction in the tube wall and solid phase, natural convection in liquid phase and volume expansion caused by density difference between solid and liquid phases are included in the numerical analysis. Existing correlation is used for estimating density-temperature relation of water, and the effect of volume expansion is reflected as fluid velocity at the interface and the free surface. As pure water has maximum density at 4.deg. C, it is found that there exists an initial temperature at which the flow direction reverses near the interface and by this effect the slope of interface becomes reversed depending on the initial temperature of water. By considering natural convection and solid-liquid density difference in the calculation, their effects on phase change process are studied and the effects of various parameters are also studied quantitatively.