High Rayleigh number thermal convection

In high-Rayleigh-number convection, numerical resolution and cell geometry are central to reliable heat-transfer predictions. We showed excellent agreement between experiments and simulations when the resolution is sufficient [1]; see figure 1a. In later work we used DNS to study how the sidewall temperature boundary condition influences heat transport (Stevens, Lohse and Verzicco, 2014).

Heat-transfer comparisons

Experimental and numerical studies of Rayleigh-Bénard convection are complementary. Experiments provide direct heat-transfer measurements with controlled uncertainties, while direct numerical simulations provide the heat transfer and the complete modeled flow field simultaneously. Experimental data can therefore validate simulations, while simulations can help explain effects observed in experiments.

Compensated Nusselt number versus Rayleigh number comparing experiments and simulations at Pr=0.7 — Figure 1: Compensated Nusselt number vs the Rayleigh number for Pr=0.7. (a) Purple stars are experimental data from Niemela et al. (2000) and the green squares are from Chavanne et al. (2001). The DNS results from Verzicco & Camussi (2003) and Amati et al. (2005) are indicated in red and high resolution results by R.J.A.M. Stevens, R. Verzicco, D. Lohse, J. Fluid Mech. 643, 495-507 (2010) are indicated by the black dots. When the vertical errorbar is not visible the error is smaller than the dot size. The results of underresolved simulations of this study are indicated by the blue dots. (b) Comparison between simulations, experiments, and the 2013 updated-prefactor version of the Grossmann-Lohse theory published in the paper The unifying theory of scaling in thermal convection: The updated prefactors, J. Fluid Mech. 730 , 295-308 (2013). Good agreement is found up to Ra=10¹¹, but still unexplained differences are observed for higher Ra. Full-resolution image 1 · Full-resolution image 2

Comparison of simulations, experiments, and Grossmann-Lohse theory predictions for heat transport up to Ra=10^11 — Figure 1: Compensated Nusselt number vs the Rayleigh number for Pr=0.7. (a) Purple stars are experimental data from Niemela et al. (2000) and the green squares are from Chavanne et al. (2001). The DNS results from Verzicco & Camussi (2003) and Amati et al. (2005) are indicated in red and high resolution results by R.J.A.M. Stevens, R. Verzicco, D. Lohse, J. Fluid Mech. 643, 495-507 (2010) are indicated by the black dots. When the vertical errorbar is not visible the error is smaller than the dot size. The results of underresolved simulations of this study are indicated by the blue dots. (b) Comparison between simulations, experiments, and the 2013 updated-prefactor version of the Grossmann-Lohse theory published in the paper The unifying theory of scaling in thermal convection: The updated prefactors, J. Fluid Mech. 730 , 295-308 (2013). Good agreement is found up to Ra=10¹¹, but still unexplained differences are observed for higher Ra. Full-resolution image 1 · Full-resolution image 2

High-Rayleigh-number DNS (2011)

Movie 1 shows what was, when published in 2011, the highest-Rayleigh-number DNS of Rayleigh-Bénard convection [3]. The simulation has Ra = 2 x 10¹² and Pr = 0.7 in a cylindrical Γ = 0.5 cell, using a computational mesh of 2701 x 671 x 2501 grid points.

Movie 1: Movie of the temperature in a vertical cut for the simulation at Ra=2x10¹², Pr=0.7, and Γ=0.5. The dimensionless time is indicated in the top of the movie. This movie is rendered with a lower resolution for the website. Additional movies of this simulation are available on the PhD thesis supplemental material page.

Resolution and domain-size checks

References

R.J.A.M. Stevens, D. Lohse, R. Verzicco,
Sidewall effects in Rayleigh-Bénard convection,
J. Fluid Mech. 741, 1-27 (2014). Read the finding
G. Ahlers, E. Bodenschatz, D. Funfschilling, S. Grossmann, X. He, D. Lohse, R.J.A.M. Stevens, R. Verzicco,
Logarithmic temperature profiles in turbulent Rayleigh-Bénard convection,
Phys. Rev. Lett. 109, 114501 (2012),
Featured in FOM News, 4 September 2012. Read the finding
R.J.A.M. Stevens, D. Lohse, R. Verzicco,
Prandtl and Rayleigh number dependence of heat transport in high Rayleigh number thermal convection,
J. Fluid Mech. 688, 31-43 (2011). Read the finding
O. Shishkina, R.J.A.M Stevens, S. Grossmann, D. Lohse,
Boundary layer structure in turbulent thermal convection and its consequences for the required numerical resolution,
New J. Phys. 12, 075022 (2010). Read the finding
R.J.A.M. Stevens, R. Verzicco, D. Lohse,
Radial boundary layer structure and Nusselt number in Rayleigh-Bénard convection,
J. Fluid Mech. 643, 495-507 (2010). Read the finding