[en] Highlights: • A thermal LBM-LES model in body-fitted coordinate system is established. • Heat transfer around a circular cylinder is modeled in a wide Reynolds number range. • The effects of Reynolds number and Prandtl number on heat transfer are revealed. -- Abstract: In order to apply the lattice Boltzmann method (LBM) for modeling passive heat transfer at high Reynolds numbers, a number of models were proposed by introducing the large eddy simulation (LES) into the LBM framework to improve numerical stability. Our study finds that the generalized form of interpolation-supplemented LBM (GILBM), likewise, can locally modify the dimensionless relaxation time, thus enhancing the numerical stability at high Reynolds numbers. Given additional advantages of the GILBM in dealing with complicated geometries and improving computing accuracy, a thermal LBM-LES model in body-fitted coordinates is established in this paper. Numerical validation is performed by investigating the flow and heat transfer around a circular cylinder over a wide range of Reynolds numbers. The obtained results agree well with both experimental and numerical data in the previous work. Meanwhile, the effects of Reynolds number and Prandtl number on thermodynamic features of flow past a circular cylinder are revealed. It is found out that when the Reynolds number exceeds the critical value, the local Nusselt number fluctuates rapidly in a specific region of the rear cylinder surface affected by the Prandtl number. In the near-wake region, the temperature field exhibits significant dependence on the Prandtl number at moderate Reynolds numbers, while such effects turn to be slight at high Reynolds numbers.