Rotating Rayleigh-Bénard convection
How does rotation influence heat transport? Our simulations and experiments show how rotation changes the flow organization and can enhance heat transfer through Ekman pumping. The enhancement depends on the Rayleigh and Prandtl numbers, rotation rate and confinement. At fixed buoyancy forcing, the Rossby number Ro is inversely proportional to the rotation rate.
Later developments
- Rotation and confinement (2022): the sampled cylindrical simulations contain separate heat-transport maxima associated with different vortex states; the sampled map does not establish universal optima.
- Optimal rotation at large Rayleigh numbers (2023): the inverse Rossby number giving the sampled maximum heat transport first increases and then decreases with Rayleigh number, while the maximum enhancement weakens at higher Rayleigh numbers in the studied cases.
Foundational results: Ekman pumping
During my PhD we investigated the influence of rotation by using direct numerical simulations and laboratory experiments in the Turbulence and Vortex Dynamics group in Eindhoven in collaboration with Herman Clercx. In the simulations we solve the incompressible Navier-Stokes and temperature equations with rotation in the Boussinesq approximation in a cylinder. Figure 2 shows that the Nusselt number can increase by up to 20% for the water cases studied. At the same imposed temperature difference, this is also a 20% increase in dimensional heat flux. The enhancement is due to Ekman pumping, which draws hot fluid out of the thermal boundary layer; see figure 1.

Figure 3 shows that the heat transfer enhancement as function of the Prandtl number has an optimum. This suggests that the effect of Ekman pumping is determined by two different mechanisms. We find that in the low Prandtl number regime the effect of Ekman pumping is limited by the larger thermal diffusivity of the fluid. Therefore the heat that is sucked into the vortices spreads out in the middle of the cell and this limits the effect of the process. In the high Prandtl number regime the kinetic boundary layer becomes thicker with respect to the thermal boundary layer. Because the vortices are formed at the edge of the kinetic boundary layer the temperature of fluid that is sucked into the vortices decreases with increasing Prandtl and this limits the effect of Ekman pumping in the high Prandtl number regime. A visualisation of these vortices close to the plates is shown in figure 4.


Data for rotating spherical convection
Simulation parameters and heat-transfer results for the spherical-convection study are available from Zenodo as version DOI 10.5281/zenodo.5034407. The record contains Rayleigh and Rossby numbers, grid resolution, rotational symmetry, and global and regional Nusselt numbers for simulations at Pr = 1 and radius ratio 0.6. The associated article is Wang et al., Diffusion-Free Scaling in Rotating Spherical Rayleigh-Bénard Convection. Zenodo's concept DOI 10.5281/zenodo.4606724 represents all versions of the dataset.
Computing project report
The Gauss Centre report Heat Transport Enhancement in Rotating Rayleigh–Bénard Convection (November 2022) describes related simulations on SuperMUC-NG and cites studies co-authored by Richard Stevens. Detlef Lohse is listed as principal investigator.
References
-
K.L. Chong, Y. Yang, S.-D. Huang, J.-Q. Zhong, R.J.A.M. Stevens, R. Verzicco, D. Lohse, K.-Q. Xia,
Confined Rayleigh-Bénard, Rotating Rayleigh-Bénard, and Double Diffusive Convection: A unifying view on turbulent transport enhancement through coherent structure manipulation,
Phys. Rev. Lett., 119, 064501 (2017) -
R.J.A.M. Stevens, H.J.H. Clercx, D. Lohse,
Heat transport and flow structure in rotating Rayleigh-Bénard convection,
European Journal of Mechanics B/Fluids 40, 41-49 (2013).
Read the finding -
R.J.A.M. Stevens, H.J.H. Clercx, D. Lohse,
Breakdown of the large-scale circulation in Γ = 1/2 rotating Rayleigh-Bénard flow,
Phys. Rev. E 86, 056311 (2012).
Read the finding -
R.P.J. Kunnen, R.J.A.M. Stevens, J. Overkamp, C. Sun, G.J.F. van Heijst, H.J.H. Clercx,
The role of Stewartson and Ekman layers in turbulent rotating Rayleigh-Bénard convection,
J. Fluid Mech. 688, 422-442 (2011).
Read the finding -
R.J.A.M. Stevens, J. Overkamp, D. Lohse, H.J.H. Clercx,
Effect of aspect ratio on vortex distribution and heat transfer in rotating Rayleigh-Bénard convection,
Phys. Rev. E 84, 056313 (2011).
Read the finding -
S. Weiss, R.J.A.M. Stevens, J.Q. Zhong, H.J.H. Clercx, D. Lohse, G. Ahlers,
Finite-size effects lead to supercritical bifurcations in turbulent rotating Rayleigh-Bénard convection,
Phys. Rev. Lett. 105, 224501 (2010),
Featured in the FOM News, 29 November 2010.
Read the finding -
R.J.A.M. Stevens, H.J.H. Clercx, D. Lohse,
Optimal Prandtl number for heat transfer in rotating Rayleigh-Bénard convection,
New J. Phys. 12, 075005 (2010),
Selected as one of New J. of Phys. best papers in 2010. Read the finding -
R.J.A.M. Stevens, H.J.H. Clercx, D. Lohse,
Boundary layers in rotating weakly turbulent Rayleigh-Bénard convection,
Phys. Fluids 22, 085103 (2010),
Read the finding -
J.Q. Zhong, R.J.A.M. Stevens, H.J.H. Clercx, R. Verzicco, D. Lohse, G. Ahlers,
Prandtl-, Rayleigh-, and Rossby-Number dependence of heat transport in turbulent rotating Rayleigh-Bénard convection,
Phys. Rev. Lett. 102, 044502 (2009),
Selected for the cover illustration of that issue.
Read the finding -
R.J.A.M. Stevens, J.Q. Zhong, H.J.H. Clercx, G. Ahlers, D. Lohse,
Transitions between turbulent states in rotating Rayleigh-Bénard convection,
Phys. Rev. Lett. 103, 024503 (2009),
Featured in the FOM News (in Dutch) FOM News, 14 July 2009
Read the finding