pith. sign in

arxiv: 2604.11670 · v2 · submitted 2026-04-13 · ⚛️ physics.flu-dyn

Influence of plume activity on thermal convection in a rectangular cell

Pith reviewed 2026-05-10 15:18 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn
keywords Rayleigh-Benard convectionplume activitythermal boundary layersviscous dissipationRayleigh numberrectangular celllarge-scale circulationturbulent convection
0
0 comments X

The pith

In regions of high plume activity within a rectangular convection cell, temperature fluctuations and dissipation rates decay more slowly with rising Rayleigh number than in low-activity zones.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper uses three-dimensional simulations of turbulent Rayleigh-Benard convection in a rectangular box to isolate how local plume activity shapes scaling behaviors. With aspect ratios that produce a stable pair of counter-rotating rolls, the setup creates distinct ejection, impact, and shear regions that remain independent of sidewalls. In the high-activity ejection zones, temperature fluctuations together with normalized thermal and viscous dissipation rates fall off less steeply as the Rayleigh number increases from 10^5 to 10^10. Boundary layers also thin faster with distance from ejection sites and with increasing Rayleigh number there than elsewhere. These local differences leave the global heat-transport laws essentially unchanged from those seen in other cell geometries.

Core claim

In this rectangular cell with width 0.8H and length 2.4H, the stable large-scale circulation fixes plume-ejection regions away from the sidewalls. There, temperature fluctuations and the normalized thermal and viscous dissipation rates decay more slowly with Rayleigh number than in impact or shear regions. Both viscous and thermal boundary layers thin rapidly with distance from the ejection zone; their local thicknesses also decline more rapidly with Rayleigh number in the ejection zone than in the other zones. Global heat transport nevertheless follows the same laws reported for cells of low to moderate aspect ratio.

What carries the argument

The stable unidirectional large-scale circulation of two counter-rotating rolls that fixes plume-ejection and shear-dominated regions with sufficient fetch for boundary-layer development.

If this is right

  • Local boundary-layer thicknesses decrease more rapidly with Rayleigh number inside the ejection region than in impact or shear zones.
  • Global heat-transport scaling remains the same as in cylindrical or other low-aspect-ratio cells despite the local variations.
  • The normalized thermal and viscous dissipation rates fall off more slowly wherever plumes move incessantly.
  • Both thermal and viscous boundary layers thin rapidly once the flow leaves the immediate ejection area.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The rectangular geometry could be used to test whether deliberately altering plume frequency in one sub-region changes the local dissipation scaling without affecting the cell-wide heat transport.
  • Averaging over regions of differing plume activity may be required to recover the familiar global scaling laws, suggesting that spatial inhomogeneity is hidden inside many existing measurements.
  • At still higher Rayleigh numbers the slower local decay in active zones might eventually alter the global exponents once the boundary layers become sufficiently thin everywhere.

Load-bearing premise

The plume-ejection regions remain independent of the sidewalls because the cell length provides enough space for velocity and thermal boundary layers to develop along the flow direction.

What would settle it

If local temperature-fluctuation amplitudes or normalized dissipation rates in the identified ejection regions were found to decay at the same rate with Rayleigh number as those in the impact and shear regions, the claimed influence of plume activity would not hold.

Figures

Figures reproduced from arXiv: 2604.11670 by Ambrish Pandey, J\"org Schumacher, Katepalli R. Sreenivasan, Matteo Parsani.

Figure 1
Figure 1. Figure 1: Instantaneous flow realisation for 𝑅𝑎 = 1010 in the closed rectangular cuboid with the red and blue colours representing, respectively, hot upwelling and cold downwelling structures. Fine thermal structures are revealed by temperature isosurfaces: hot and cold structures correspond to 𝑇 = 0.65 and 𝑇 = 0.35, respectively. Vertical velocity contours, shown in the vicinity of the sidewall at 𝑥 ≈ 𝐿𝑥 and on the… view at source ↗
Figure 2
Figure 2. Figure 2: An illustration depicting various regions of the flow in the present configuration, using a [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: (a) Global heat transport, 𝑁𝑢, increases with the thermal driving as 𝑅𝑎0.29 in the rectangular box. Present results agree excellently with those for 𝑃𝑟 = 1 from Blass et al. (2021) and van Reeuwijk et al. (2008) obtained respectively in Γ = 32 and Γ = 4 horizontally-periodic square cuboids as well as those for 𝑃𝑟 = 0.7 from Samuel et al. (2024). (b) Reynolds number based on 𝑢RMS increases almost as √ 𝑅𝑎. D… view at source ↗
Figure 4
Figure 4. Figure 4: (a) Fractions of the kinetic energy corresponding to velocity components as functions of [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Scaling relations in quiet (red) and active (blue) bulk regions. (a) Reynolds numbers based on the [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Mean thickness of the viscous layer at the horizontal plates, estimated using the slope method [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: (a) Normalised wall-parallel velocity profile [PITH_FULL_IMAGE:figures/full_fig_p011_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: (a) Thermal fluctuation thickness 𝛿𝜎 (symbols) decreases with 𝑅𝑎 with the best fit (solid curve) suggesting a 3.74𝑅𝑎−0.30 scaling for the entire range. The 𝛿𝜎 remains slightly smaller than the slope thickness 𝛿𝑇 that decreases as 3.4𝑅𝑎−0.29 (dashed line). The present data agrees well with the thermal fluctuation thickness observed for 𝑃𝑟 = 0.7 in a square cuboid of Γ = 4 (Samuel et al. 2024), shown as oran… view at source ↗
Figure 9
Figure 9. Figure 9: Structure of the viscous boundary layer at the bottom plate. (a) Local thickness estimated using the [PITH_FULL_IMAGE:figures/full_fig_p013_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Structure of the thermal boundary layer at the bottom plate. Both the local fluctuation thickness [PITH_FULL_IMAGE:figures/full_fig_p013_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Scaling of (a) local VBL thickness and (b) TBL thickness as functions of [PITH_FULL_IMAGE:figures/full_fig_p014_11.png] view at source ↗
read the original abstract

We present three-dimensional direct numerical simulations of turbulent Rayleigh-B\'enard convection in a closed rectangular box whose width $L_y$ and length $L_x$ are 0.8 and 2.4 times the height $H$, respectively. The Rayleigh number $Ra$ varies from $10^5$ to $10^{10}$, and the Prandtl number is unity. The advantages of the present configuration are: (a) A relatively stable unidirectional large-scale circulation, consisting of two counter-rotating rolls, fills the cell and fixes the thermal plume ejection- and shear-dominated regions, in contrast to those in closed cylindrical cells. (b) The regions of plume ejection are essentially independent of the sidewalls so that their autonomous existence can be studied. This is because there is some space, or "fetch", for the velocity and thermal boundary layers to develop along the length. (c) This geometry allows one to study the influence of locally thin and thick boundary layers (which follow larger or smaller plume activity) on the scaling of convection properties. In regions of larger plume activity (defined by an incessant movement of plumes), the temperature fluctuation as well as the normalised thermal and viscous dissipation rates decay more slowly with $Ra$ than in regions of lower activity. Both viscous and thermal boundary layers thin down rapidly with increasing distance from the plume ejection region. The local thicknesses of both boundary layers decline more rapidly with $Ra$ in the ejection region than in regions of impact and shear, where they are similar to each other. Despite these details, the global heat transport laws are practically the same as those in other configurations of low to moderate aspect ratios.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript reports three-dimensional direct numerical simulations of Rayleigh-Bénard convection in a rectangular cell of dimensions L_x/H=2.4 and L_y/H=0.8 at Pr=1, with Rayleigh numbers spanning 10^5 to 10^10. It claims that a stable two-roll large-scale circulation fixes distinct plume-ejection, impact, and shear regions; in high-plume-activity zones the temperature fluctuations and normalized thermal/viscous dissipation rates decay more slowly with Ra, while local viscous and thermal boundary-layer thicknesses thin more rapidly with Ra in ejection regions than in impact/shear regions (where the two thicknesses are comparable). Global heat transport is reported to follow the same scaling as in other low-to-moderate aspect-ratio cells.

Significance. If the local distinctions survive scrutiny, the work supplies a concrete geometric mechanism for why local scaling exponents can differ from global ones in turbulent convection, by tying them to plume activity and boundary-layer development along a controlled large-scale flow. The rectangular setup that stabilizes the circulation is a useful experimental design choice for isolating these effects.

major comments (2)
  1. [Abstract and §1] Abstract and §1 (Introduction): The central claim that 'the regions of plume ejection are essentially independent of the sidewalls' because 'there is some space, or fetch, for the velocity and thermal boundary layers to develop along the length' is load-bearing for all local scaling results. With L_x/H=2.4 and a stable two-roll LSC the effective streamwise distance from a sidewall to the interior of an ejection region is only ~1.2H. The manuscript must demonstrate, e.g., via streamwise profiles of mean velocity, temperature variance, or enstrophy, that sidewall-induced perturbations have decayed before the analyzed ejection zones; without such evidence the reported differences in Ra-scaling between regions remain vulnerable to contamination.
  2. [Results (local statistics)] Results (local statistics and boundary-layer sections): The slower Ra-decay of temperature fluctuations and normalized dissipation rates in high-plume-activity regions, together with the differing local boundary-layer scalings, are presented without error bars, statistical convergence tests, or grid-resolution studies. At Ra up to 10^10 such checks are required to confirm that the claimed distinctions in scaling exponents are not resolution artifacts or consequences of insufficient temporal averaging.
minor comments (2)
  1. [Abstract] Abstract: The statement that 'the local thicknesses of both boundary layers decline more rapidly with Ra in the ejection region than in regions of impact and shear, where they are similar to each other' would be clearer if the approximate scaling exponents or the relevant figure panels were referenced directly.
  2. [Throughout] Notation: Define the abbreviations for the three regions (ejection, impact, shear) at first use and ensure consistent symbols for local versus global quantities throughout the text and figures.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment below and agree that the suggested additions will strengthen the work.

read point-by-point responses
  1. Referee: [Abstract and §1] Abstract and §1 (Introduction): The central claim that 'the regions of plume ejection are essentially independent of the sidewalls' because 'there is some space, or fetch, for the velocity and thermal boundary layers to develop along the length' is load-bearing for all local scaling results. With L_x/H=2.4 and a stable two-roll LSC the effective streamwise distance from a sidewall to the interior of an ejection region is only ~1.2H. The manuscript must demonstrate, e.g., via streamwise profiles of mean velocity, temperature variance, or enstrophy, that sidewall-induced perturbations have decayed before the analyzed ejection zones; without such evidence the reported differences in Ra-scaling between regions remain vulnerable to contamination.

    Authors: We agree that explicit demonstration of the decay of sidewall-induced perturbations is necessary to support the independence of the plume-ejection regions. In the revised manuscript we will add streamwise profiles of mean velocity, temperature variance, and enstrophy (and, if space permits, enstrophy) extracted along the cell length. These profiles will show that perturbations originating at the sidewalls have decayed to negligible levels before reaching the interior of the ejection zones, thereby confirming that the reported local scaling differences are not contaminated by sidewall effects. revision: yes

  2. Referee: [Results (local statistics)] Results (local statistics and boundary-layer sections): The slower Ra-decay of temperature fluctuations and normalized dissipation rates in high-plume-activity regions, together with the differing local boundary-layer scalings, are presented without error bars, statistical convergence tests, or grid-resolution studies. At Ra up to 10^10 such checks are required to confirm that the claimed distinctions in scaling exponents are not resolution artifacts or consequences of insufficient temporal averaging.

    Authors: We acknowledge that quantitative error bars, convergence diagnostics, and resolution verification are essential at the highest Rayleigh numbers. In the revised version we will include error bars on all local scaling plots, report the duration and sensitivity of the temporal averaging windows used for each Ra, and add a dedicated subsection (or appendix) summarizing grid-resolution checks and Kolmogorov-scale resolution ratios. These additions will confirm that the observed differences in scaling exponents are statistically robust and not numerical artifacts. revision: yes

Circularity Check

0 steps flagged

No circularity: all quantities obtained by direct integration of the governing equations.

full rationale

The paper performs direct numerical simulations of the Boussinesq equations in a rectangular domain and extracts all reported statistics (temperature fluctuations, normalized dissipation rates, local boundary-layer thicknesses, and their Ra scalings) by spatial and temporal averaging of the computed fields. Plume-activity regions are identified directly from the instantaneous flow structures without any auxiliary fitting or redefinition of inputs as outputs. No self-citations, uniqueness theorems, or ansatzes are invoked to close the central claims; the geometry and aspect ratio are chosen explicitly to enable the reported comparison, but the comparison itself remains a numerical observation rather than a tautological reduction. The analysis is therefore self-contained as a controlled numerical experiment.

Axiom & Free-Parameter Ledger

3 free parameters · 2 axioms · 0 invented entities

The central claims rest on numerical solution of the incompressible Navier-Stokes equations under the Boussinesq approximation with three externally chosen parameters (Ra range, Pr = 1, and the two aspect ratios) that define the computational domain.

free parameters (3)
  • Rayleigh number = 10^5 to 10^10
    Varied parametrically from 10^5 to 10^10 to extract scaling trends
  • Prandtl number = 1
    Fixed at unity for all runs
  • Aspect ratios = 2.4 and 0.8
    Lx/H = 2.4 and Ly/H = 0.8 chosen to produce stable rolls and boundary-layer fetch
axioms (2)
  • standard math Incompressible Navier-Stokes equations with Boussinesq buoyancy term govern the flow
    Invoked implicitly as the basis for all direct numerical simulations of Rayleigh-Bénard convection
  • domain assumption No-slip velocity and fixed-temperature boundary conditions on all walls
    Standard closed-box setup required to produce the described large-scale circulation

pith-pipeline@v0.9.0 · 5616 in / 1621 out tokens · 53966 ms · 2026-05-10T15:18:54.583646+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

55 extracted references · 55 canonical work pages

  1. [1]

    , " * write output.state after.block = add.period write newline

    ENTRY address author booktitle chapter edition editor howpublished institution journal key month note number organization pages publisher school series title type volume year eprint label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 'mid.sentence ...

  2. [2]

    write newline

    " write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in capitalize " " * FUNCT...

  3. [3]

    , Grossmann, S

    Ahlers, G. , Grossmann, S. & Lohse, D. 2009 Heat transfer and large scale dynamics in turbulent R ayleigh- B \'enard convection . Rev. Mod. Phys. 81 , 503--537

  4. [4]

    , Tilgner, A

    Belmonte, A. , Tilgner, A. & Libchaber, A. 1993 Boundary layer length scales in thermal turbulence . Phys. Rev. Lett. 70 , 4067--4070

  5. [5]

    , Tilgner, A

    Belmonte, A. , Tilgner, A. & Libchaber, A. 1994 Temperature and velocity boundary layers in turbulent convection . Phys. Rev. E 50 , 269--279

  6. [6]

    , Samtaney, R

    Bhattacharya, S. , Samtaney, R. & Verma, M. K. 2019 Scaling and spatial intermittency of thermal dissipation in turbulent convection . Phys. Fluids 31 (7), 075104

  7. [7]

    , Verzicco, R

    Blass, A. , Verzicco, R. , Lohse, D. , Stevens, R. J. A. M. & Krug, D. 2021 Flow organisation in laterally unconfined R ayleigh– B \'enard turbulence . J. Fluid Mech. 906 , A26

  8. [8]

    Busse, F. H. 1978 Non-linear properties of thermal convection . Rep. Prog. Phys. 41 (12), 1929--1967

  9. [9]

    , Gunaratne , G

    Castaing , B. , Gunaratne , G. , Kadanoff , L. , Libchaber , A. & Heslot , F. 1989 Scaling of hard thermal turbulence in R ayleigh- B \'e nard convection . J. Fluid Mech. 204 , 1--30

  10. [10]

    & Verma, M

    Chandra, M. & Verma, M. K. 2013 Flow reversals in turbulent convection via vortex reconnections . Phys. Rev. Lett. 110 , 114503

  11. [11]

    Chong, K. L. & Xia, K.-Q. 2016 Exploring the severely confined regime in R ayleigh– B \'enard convection . J. Fluid Mech. 805 , R4

  12. [12]

    Emran , M. S. & Schumacher , J. 2012 Conditional statistics of thermal dissipation rate in turbulent R ayleigh- B \'e nard convection . Eur. Phys. J. E 35 , 12108

  13. [13]

    Fischer , P. F. 1997 An overlapping S chwarz method for spectral element solution of the incompressible N avier- S tokes equations . J. Comp. Phys. 133 (1), 84--101

  14. [14]

    , Bao, Y

    He, J.-C. , Bao, Y. & Chen, X. 2024 Turbulent boundary layers in thermal convection at moderately high R ayleigh numbers . Phys. Fluids 36 (2), 025140

  15. [15]

    , Tong , P

    He , X. , Tong , P. & Xia , K. Q. 2007 Measured thermal dissipation field in turbulent R ayleigh- B \'e nard convection . Phys. Rev. Lett. 98 , 144501

  16. [16]

    & Xia, K.-Q

    He, Y.-H. & Xia, K.-Q. 2019 Temperature fluctuation profiles in turbulent thermal convection: A logarithmic dependence versus a power-law dependence . Phys. Rev. Lett. 122 , 014503

  17. [17]

    Howard, L. N. 1972 Bounds on flow quantities . Annu. Rev. Fluid Mech. 4 (1), 473--494

  18. [18]

    Huang , S. D. & Xia , K. Q. 2016 Effects of geometric confinement in quasi-2-D turbulent R ayleigh- B \'e nard convection . J. Fluid Mech. 794 , 639--654

  19. [19]

    , Schindler, F

    Kim, N. , Schindler, F. , Vogt, T. & Eckert, S. 2024 Thermal boundary layer dynamics in low-Prandtl-number R ayleigh– B \'enard convection . J. Fluid Mech. 994 , A4

  20. [20]

    , Shang , X.-D

    Lam , S. , Shang , X.-D. , Zhou , S.-Q. & Xia , K.-Q. 2002 P randtl number dependence of the viscous boundary layer and the R eynolds numbers in R ayleigh- B \'e nard convection . Phys. Rev. E 65 , 066306

  21. [21]

    , Huang, S.-D

    Li, X.-M. , Huang, S.-D. , Ni, R. & Xia, K.-Q. 2021 Lagrangian velocity and acceleration measurements in plume-rich regions of turbulent R ayleigh– B \'enard convection . Phys. Rev. Fluids 6 , 053503

  22. [22]

    & Shishkina, O

    Lohse, D. & Shishkina, O. 2024 Ultimate R ayleigh- B \'enard turbulence . Rev. Mod. Phys. 96 , 035001

  23. [23]

    & Xia, K.-Q

    Lui, S.-L. & Xia, K.-Q. 1998 Spatial structure of the thermal boundary layer in turbulent convection . Phys. Rev. E 57 , 5494--5503

  24. [24]

    Mishra , P. K. , De , A. K. , Verma , M. K. & Eswaran , V. 2011 Dynamics of reorientations and reversals of large-scale flow in R ayleigh- B \'e nard convection . J. Fluid Mech. 668 , 480--499

  25. [25]

    Niemela , J. J. , Skrbek , L. , Sreenivasan , K. R. & Donnelly , R. J. 2000 Turbulent convection at very high R ayleigh numbers . Nature 404 , 837--840

  26. [26]

    Niemela , J. J. , Skrbek , L. , Sreenivasan , K. R. & Donnelly , R. J. 2001 The wind in confined thermal convection . J. Fluid Mech. 449 , 169

  27. [27]

    Niemela, J. J. & Sreenivasan, K. R. 2003 R ayleigh-number evolution of large-scale coherent motion in turbulent convection . Europhys. Lett. 62 (6), 829--833

  28. [28]

    2021 Thermal boundary layer structure in low- P randtl-number turbulent convection

    Pandey, A. 2021 Thermal boundary layer structure in low- P randtl-number turbulent convection . J. Fluid Mech. 910 , A13

  29. [29]

    , Krasnov, D

    Pandey, A. , Krasnov, D. , Schumacher, J. , Samtaney, R. & Sreenivasan, K. R. 2022 a\/ Similarities between characteristics of convective turbulence in confined and extended domains . Physica D 442 , 133537

  30. [30]

    , Krasnov, D

    Pandey, A. , Krasnov, D. , Sreenivasan, K. R. & Schumacher, J. 2022 b\/ Convective mesoscale turbulence at very low P randtl numbers . J. Fluid Mech. 948 , A23

  31. [31]

    , Scheel, J

    Pandey, A. , Scheel, J. D. & Schumacher, J. 2018 Turbulent superstructures in R ayleigh- B \'enard convection . Nat. Commun. 9 , 2118

  32. [32]

    & Sreenivasan, K

    Pandey, A. & Sreenivasan, K. R. 2024 Turbulent convection in rotating slender cells . J. Fluid Mech. 999 , A28

  33. [33]

    & Sreenivasan, K

    Pandey, A. & Sreenivasan, K. R. 2025 Transient and steady convection in two dimensions . J. Fluid Mech. 1015 , A42

  34. [34]

    & Verma, M

    Pandey, A. & Verma, M. K. 2016 Scaling of large-scale quantities in R ayleigh- B \'e nard convection . Phys. Fluids 28 (9), 095105

  35. [35]

    & Xia, K.-Q

    Qiu, X.-L. & Xia, K.-Q. 1998 a\/ Spatial structure of the viscous boundary layer in turbulent convection . Phys. Rev. E 58 , 5816--5820

  36. [36]

    & Xia, K.-Q

    Qiu, X.-L. & Xia, K.-Q. 1998 b\/ Viscous boundary layers at the sidewall of a convection cell . Phys. Rev. E 58 , 486--491

  37. [37]

    , Jonker, H

    van Reeuwijk, M. , Jonker, H. J. J. & Hanjali c \' c , K. 2008 Wind and boundary layers in R ayleigh- B \'enard convection. i. analysis and modeling . Phys. Rev. E 77 , 036311

  38. [38]

    Samuel, R. J. , Bode, M. , Scheel, J. D. , Sreenivasan, K. R. & Schumacher, J. 2024 No sustained mean velocity in the boundary region of plane thermal convection . J. Fluid Mech. 996 , A49

  39. [39]

    Samuel, R. J. & Schumacher, J. 2025 Plane-layer Rayleigh–B \' e nard convection up to Ra = 10^ 11 : near-wall fluctuations and role of initial conditions . Fluid Dyn. Res. 57 , 061401

  40. [40]

    Scheel , J. D. , Emran , M. S. & Schumacher , J. 2013 Resolving the fine-scale structure in turbulent R ayleigh- B \'e nard convection . New J. Phys. 15 , 113063

  41. [41]

    Scheel, J. D. , Kim, E. & White, K. R. 2012 Thermal and viscous boundary layers in turbulent R ayleigh– B \'enard convection . J. Fluid Mech. 711 , 281–305

  42. [42]

    , Bandaru, V

    Schumacher, J. , Bandaru, V. , Pandey, A. & Scheel, J. D. 2016 Transitional boundary layers in low- P randtl-number convection . Phys. Rev. Fluids 1 , 084402

  43. [43]

    Shevkar, P. P. , Samuel, R. J. , Zinchenko, G. , Bode, M. , Schumacher, J. & Sreenivasan, K. R. 2025 Hierarchical network of thermal plumes and their dynamics in turbulent R ayleigh- B \'enard convection . Proc. Natl. Acad. Sci. USA 122 (32), e2502972122

  44. [44]

    , Emran, M

    Shi, N. , Emran, M. S. & Schumacher, J. 2012 Boundary layer structure in turbulent R ayleigh– B \'e nard convection . J. Fluid Mech. 706 , 5–33

  45. [45]

    , Stevens , R

    Shishkina , O. , Stevens , R. , Grossmann , S. & Lohse , D. 2010 Boundary layer structure in turbulent thermal convection and its consequences for the required numerical resolution . New J. Phys. 12 , 075022

  46. [46]

    Shraiman, B. I. & Siggia, E. D. 1990 Heat transport in high- R ayleigh-number convection . Phys. Rev. A 42 , 3650--3653

  47. [47]

    Siggia, E. D. 1994 High R ayleigh number convection . Annu. Rev. Fluid Mech. 26 (1), 137--168

  48. [48]

    Sreenivasan, K. R. , Bershadskii, A. & Niemela, J. J. 2002 Mean wind and its reversal in thermal convection . Phys. Rev. E 65 , 056306

  49. [49]

    , Verzicco , R

    Stevens , R. , Verzicco , R. & Lohse , D. 2010 Radial boundary layer structure and N usselt number in R ayleigh- B \'e nard convection . J. Fluid Mech. 643 , 495--507

  50. [50]

    van der Poel , E. P. , Ostilla-M\' o nico , R. , Verzicco , R. , Grossmann , S. & Lohse , D. 2015 Logarithmic mean temperature profiles and their connection to plume emissions in turbulent R ayleigh- B \'e nard convection . Phys. Rev. Lett. 115 , 154501

  51. [51]

    , Shishkina, O

    Wagner, S. , Shishkina, O. & Wagner, C. 2012 Boundary layers and wind in cylindrical R ayleigh– B \'e nard cells . J. Fluid Mech. 697 , 336–366

  52. [52]

    & Libchaber, A

    Wu, X.-Z. & Libchaber, A. 1992 Scaling relations in thermal turbulence: The aspect-ratio dependence . Phys. Rev. A 45 , 842--845

  53. [53]

    , Zhou, Q

    Xi, H.-D. , Zhou, Q. & Xia, K.-Q. 2006 Azimuthal motion of the mean wind in turbulent thermal convection . Phys. Rev. E 73 , 056312

  54. [54]

    , Zhang, L

    Xu, F. , Zhang, L. & Xia, K.-Q. 2024 Experimental measurement of spatio-temporally resolved energy dissipation rate in turbulent R ayleigh– B \'enard convection . J. Fluid Mech. 984 , A8

  55. [55]

    , Mathai, V

    Zhu, X. , Mathai, V. , Stevens, R. J. A. M. , Verzicco, R. & Lohse, D. 2018 Transition to the ultimate regime in two-dimensional R ayleigh- B \'enard convection . Phys. Rev. Lett. 120 , 144502