Recovery by equilibrium in patches (REP) technique is adopted for the non-linear transient heat conduction adaptive FEM analysis of casting system, which takes account of thermal conductivity of the interface. In order to test the present program, ANSYS and superconvergent patch recovery (SPR) methods are also employed for analysing the same problem respectively. The results of numerical analysis are compared with the experiment results, which show that REP method is cost effective.