In this study, we compare different diffuse and sharp interface schemes of direct-forcing immersed boundary — thermal lattice Boltzmann method (IB-TLBM) for non-Newtonian flow over a heated circular cylinder. Both effects of the discrete lattice and the body force on the momentum and energy equations are considered, by applying the split-forcing Lattice Boltzmann equations. A new technique based on predetermined parameters of direct forcing IB-TLBM is presented for computing the Nusselt number. The study covers both steady and unsteady regimes (20<Re<80) in the power-law index range of 0.6<n <1.4, encompassing both shear-thinning and shear-thickening non-Newtonian fluids. The numerical scheme, hydrodynamic approach and thermal parameters of different interface schemes are compared in both steady and unsteady cases. It is found that the sharp interface scheme is a suitable and possibly competitive method for thermal-IBM in terms of accuracy and computational cost.