This paper deals with the problem of conjugate heat transfer of a rotary disk which is uniformly heated. Computational fluid dynamics simulations are performed for different ranges of rotational Reynolds number (20,000-50,000) and heat flux (100, 200, and 400 kW/m 3) on this surface. Two numerical approaches including multiple reference frame and sliding mesh are employed to simulate the interaction between rotary and stationary domains. It is shown that both methods yield satisfactory results, and an excellent agreement with the results obtained from the correlation available in the literature is also found. Furthermore, the effects of involved parameters on average and local convective heat transfer are investigated. Results show that the temperature profiles on the disk surface remain rather unchanged for all ranges of those parameters under study here. It is also shown that local Nusselt number decreases more steeply for higher ratios of fluid conductivity to solid conductivity.