This paper presents a finite element method for solving non-linear heat conduction and thermal stress problems. The nonlinearity is caused by the variation of thermo-physical and mechanical properties of the material with respect to temperature. The distribution of temperature and stress in a cooling body can be followed from the onset of solidification. The method is capable of handling complex geometries and initial non-linear boundary conditions.