High-order LES of turbulent heat transfer in a rotor–stator cavity

International audienceThe present work examines the turbulent flow in an enclosed rotor–stator system subjected to heat transfer effects. Besides their fundamental importance as three-dimensional prototype flows, such flows arise in many industrial applications but also in many geophysical and astrophysical settings. Large eddy simulations(LES) are here performed using a spectral vanishing viscosity technique. The LES results have already been favorably compared to velocity measurements in the isothermal case (Séverac et al. 2007) for a large range of Reynolds numbers in an annular cavity of large aspect ratio and weak curvature parameter. The purpose of this paper is to extend these previous results in the non-isothermal case using the Boussinesq approximation to take into account the buoyancy effects. Thus, the effects of thermal convection have been examined for a turbulent flow of air in the same rotor–stator system for Rayleigh numbers up to 100 millions. These LES results provide accurate, instantaneous quantities which are of interest in understanding the physics of turbulent flows and heat transfers in an interdisk cavity. Even at high Rayleigh numbers, the structure of the iso-values of the instantaneous normal temperature gradient at the disk surfaces resembles the one of the iso-values of the tangential velocity with large spiral arms along the rotor and more thin structures along the stator. The averaged results show small effects of density variation on the mean and turbulent fields. The turbulent Prandtl number is a decreasing function of the distance to the wall with 1.4 close to the disks and about 0.3 in the outer layers. The local Nusselt number is found to be proportional to the local Reynolds number to the power 0.7. The evolution of the averaged Bolgiano length scale with the Rayleigh number indicates that temperature fluctuations may have a large influence on the dynamics only at the largest scales of the system for Ra larger than 10 millions, since the averaged Bolgiano length scale remains lower than the thermal boundary layer thicknesses