Numerical simulation was applied to investigate the Taylor vortex flow inside the concentric cylinders with a constant radial temperature gradient. The reliability of numerical simulation method was verified by the experimental results of PIV. The radial velocity and temperature distribution in plain and 12-slit model at different axial locations were compared, and the heat flux distributions along the inner cylinder wall at different work conditions were obtained. In the plain model, the average surface heat flux of inner cylinder increased with the inner cylinder rotation speed. In slit model, the slit wall significantly changed the distribution of flow field and temperature in the annulus gap, and the radial flow was strengthen obviously, which promoted the heat transfer process at the same working condition.