The numerical manifold method (NMM), combined with the precise time integration method (PTIM), is proposed for thermal shock fracture analysis. The temperature and displacement discontinuity across crack faces is naturally portrayed attributing to the cover systems in the NMM. The crack tip singularities are characterized through the use of asymptotic bases in the approximations. The discrete equations for transient thermal analysis are firstly solved with the PTIM and then the thermoelastic study is performed. With the interaction integral, the stress intensity factors are computed. Several examples are tested and the nice consistency between the present and existing results is found.