Multi-scale steady solution for Rayleigh–Bénard convection

We found a multi-scale steady solution of the Boussinesq equations for Rayleigh–Bénard convection in a three-dimensional periodic domain between horizontal plates with a constant temperature difference. This was realised using a homotopy from the wall-to-wall optimal transport solution provided by Motoki et al. (J. Fluid Mech., vol. 851, 2018, R4). A connected steady solution, which is a consequence of bifurcation from a thermal conduction state at Rayleigh number $Ra\sim 10^{3}$, is tracked up to $Ra\sim 10^{7}$ using a Newton–Krylov iteration. The three-dimensional exact coherent thermal convection exhibits a scaling of $Nu\sim Ra^{0.31}$ (where $Nu$ is the Nusselt number) as well as multi-scale thermal plume and vortex structures, which are quite similar to those in turbulent Rayleigh–Bénard convection. The mean temperature profiles and the root-mean-square of the temperature and velocity fluctuations are in good agreement with those of the turbulent states. Furthermore, the energy spectrum follows Kolmogorov's $-5/3$ scaling law with a consistent prefactor, and the energy transfer to small scales with a nearly constant flux in the wavenumber space is in accordance with the turbulent energy transfer.