Abstract: In this paper, based on non-Fouriers heat conduction equation, using the image method and the complex function method, the multiple scattering of the boundaries of the slab material is solved. The temperature field of the thermal wave scattering from two subsurface cylindrical defects in the semi-infinite slab materials is investigated, a general solution of the scattering field of thermal waves is obtained. The frontal surface temperature of the slab materials under different incident wave number, thermal diffusion length, buried depth and the thickness of the slab are analyzed and illustrated graphically. In addition, the thermal waves are excited at frontal surfaces of the slab materials by modulated ultra-short laser pulse with an adiabatic boundary condition for defects. The analytical methods and the numerical results are expected to provide references for the heat transfer analysis of engineering material, the thermal wave imaging, physical inverse problem and non-destructive evaluation of defects in materials.