Influence of plume activity on thermal convection in a rectangular cell
Pith reviewed 2026-05-10 15:18 UTC · model grok-4.3
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.
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
- 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
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.
Referee Report
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)
- [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.
- [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)
- [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.
- [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
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
-
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
-
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
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
free parameters (3)
- Rayleigh number =
10^5 to 10^10
- Prandtl number =
1
- Aspect ratios =
2.4 and 0.8
axioms (2)
- standard math Incompressible Navier-Stokes equations with Boussinesq buoyancy term govern the flow
- domain assumption No-slip velocity and fixed-temperature boundary conditions on all walls
Reference graph
Works this paper leans on
-
[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]
" 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]
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
work page 2009
-
[4]
Belmonte, A. , Tilgner, A. & Libchaber, A. 1993 Boundary layer length scales in thermal turbulence . Phys. Rev. Lett. 70 , 4067--4070
work page 1993
-
[5]
Belmonte, A. , Tilgner, A. & Libchaber, A. 1994 Temperature and velocity boundary layers in turbulent convection . Phys. Rev. E 50 , 269--279
work page 1994
-
[6]
Bhattacharya, S. , Samtaney, R. & Verma, M. K. 2019 Scaling and spatial intermittency of thermal dissipation in turbulent convection . Phys. Fluids 31 (7), 075104
work page 2019
-
[7]
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
work page 2021
-
[8]
Busse, F. H. 1978 Non-linear properties of thermal convection . Rep. Prog. Phys. 41 (12), 1929--1967
work page 1978
-
[9]
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
work page 1989
-
[10]
Chandra, M. & Verma, M. K. 2013 Flow reversals in turbulent convection via vortex reconnections . Phys. Rev. Lett. 110 , 114503
work page 2013
-
[11]
Chong, K. L. & Xia, K.-Q. 2016 Exploring the severely confined regime in R ayleigh– B \'enard convection . J. Fluid Mech. 805 , R4
work page 2016
-
[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
work page 2012
-
[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
work page 1997
- [14]
-
[15]
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
work page 2007
-
[16]
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
work page 2019
-
[17]
Howard, L. N. 1972 Bounds on flow quantities . Annu. Rev. Fluid Mech. 4 (1), 473--494
work page 1972
-
[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
work page 2016
-
[19]
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
work page 2024
-
[20]
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
work page 2002
-
[21]
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
work page 2021
-
[22]
Lohse, D. & Shishkina, O. 2024 Ultimate R ayleigh- B \'enard turbulence . Rev. Mod. Phys. 96 , 035001
work page 2024
-
[23]
Lui, S.-L. & Xia, K.-Q. 1998 Spatial structure of the thermal boundary layer in turbulent convection . Phys. Rev. E 57 , 5494--5503
work page 1998
-
[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
work page 2011
-
[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
work page 2000
-
[26]
Niemela , J. J. , Skrbek , L. , Sreenivasan , K. R. & Donnelly , R. J. 2001 The wind in confined thermal convection . J. Fluid Mech. 449 , 169
work page 2001
-
[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
work page 2003
-
[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
work page 2021
-
[29]
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
work page 2022
-
[30]
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
work page 2022
-
[31]
Pandey, A. , Scheel, J. D. & Schumacher, J. 2018 Turbulent superstructures in R ayleigh- B \'enard convection . Nat. Commun. 9 , 2118
work page 2018
-
[32]
Pandey, A. & Sreenivasan, K. R. 2024 Turbulent convection in rotating slender cells . J. Fluid Mech. 999 , A28
work page 2024
-
[33]
Pandey, A. & Sreenivasan, K. R. 2025 Transient and steady convection in two dimensions . J. Fluid Mech. 1015 , A42
work page 2025
-
[34]
Pandey, A. & Verma, M. K. 2016 Scaling of large-scale quantities in R ayleigh- B \'e nard convection . Phys. Fluids 28 (9), 095105
work page 2016
-
[35]
Qiu, X.-L. & Xia, K.-Q. 1998 a\/ Spatial structure of the viscous boundary layer in turbulent convection . Phys. Rev. E 58 , 5816--5820
work page 1998
-
[36]
Qiu, X.-L. & Xia, K.-Q. 1998 b\/ Viscous boundary layers at the sidewall of a convection cell . Phys. Rev. E 58 , 486--491
work page 1998
-
[37]
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
work page 2008
-
[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
work page 2024
-
[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
work page 2025
-
[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
work page 2013
-
[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
work page 2012
-
[42]
Schumacher, J. , Bandaru, V. , Pandey, A. & Scheel, J. D. 2016 Transitional boundary layers in low- P randtl-number convection . Phys. Rev. Fluids 1 , 084402
work page 2016
-
[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
work page 2025
-
[44]
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
work page 2012
-
[45]
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
work page 2010
-
[46]
Shraiman, B. I. & Siggia, E. D. 1990 Heat transport in high- R ayleigh-number convection . Phys. Rev. A 42 , 3650--3653
work page 1990
-
[47]
Siggia, E. D. 1994 High R ayleigh number convection . Annu. Rev. Fluid Mech. 26 (1), 137--168
work page 1994
-
[48]
Sreenivasan, K. R. , Bershadskii, A. & Niemela, J. J. 2002 Mean wind and its reversal in thermal convection . Phys. Rev. E 65 , 056306
work page 2002
-
[49]
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
work page 2010
-
[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
work page 2015
-
[51]
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
work page 2012
-
[52]
Wu, X.-Z. & Libchaber, A. 1992 Scaling relations in thermal turbulence: The aspect-ratio dependence . Phys. Rev. A 45 , 842--845
work page 1992
- [53]
-
[54]
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
work page 2024
-
[55]
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
work page 2018
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.