The purpose of this research is to develope a FEM program for analyzing solidification processes of axisymmetric casting, considering phase changes and the contact between the metal and mold. Tempera- ture recovery method is employed fro considering the phase changes releasing the latent heat and the coin- cident node method is used for calculating the amount of heat transfer between the metal and mold. Tan- gent modulus algorithm is adopted for calculating flow stress and a gap element is employed for modeling the interface between the mold and metal in finding deformed shapes. In order to verify the developed program, axisymmetric aluminum and steel casting processes are simulated. Temperature distribution, phase front position, and shrinkage and porosity creation are compared with measurements, FIDAP results, and good agreements are examined.