An axisymmetric model has been developed to simulate conduction phase change during laser processing based on enthalpy method. The spatial and temporal profile of the laser beam and the temperature dependent properties of the material have been taken into account in the model. A Gaussian beam and a flat beam are used, respectively, to simulate the interaction with iron as target and the results for vapor-liquid and liquid-solid interface evolutions and temperature field development in the material are obtained.