Finite element analyses for three dimensional problems in solidification have been made. For handling the release of latent heat of solidification, the temperature recovery method was applied. In analyses, the accuracy of the computed results depended upon the element size and the time increment used in simulations. An interface heat transfer coefficient between the mold and the casting was determined numerically as a function of time and used in calculations. The predicted temperatures were in good agreement with the available solutions in the literature.