Heat transfer in rough-wall turbulent thermal convection in the ultimate regime
Pith reviewed 2026-05-25 08:52 UTC · model grok-4.3
The pith
Rough-wall heat transfer follows Nu ~ Ra^{0.42} scaling and never reaches the ultimate regime at finite Ra.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Using logarithmic temperature profiles with a roughness-induced shift in analyses of forced convection over rough walls predicts an effective scaling Nu ∼ Ra^{0.42} corresponding to Ch ∼ Re^{-0.16}. Unlike the skin-friction coefficient which becomes independent of Re due to pressure drag, the absence of an analog to pressure drag in the temperature advection equation produces this different scaling. Consequently the asymptotic ultimate regime where Nu ∼ Ra^{1/2} will never be reached for heat transfer at finite Ra.
What carries the argument
Logarithmic temperature profiles with a roughness-induced shift, applied to extend prior ultimate-regime analyses to rough-wall data and derive the effective scaling.
If this is right
- Heat transfer scales differently from momentum transfer in the fully rough regime.
- The effective exponent for Nu versus Ra is 0.42 rather than approaching 0.5.
- The ultimate regime is not attained for heat transfer even as Ra increases.
- Stanton number decreases as Re to the power -0.16.
Where Pith is reading between the lines
- The analysis implies that increasing Ra further will not change the 0.42 exponent toward the ultimate value.
- Similar logarithmic assumptions might be tested in other scalar transport problems with rough boundaries.
- Experiments could check whether temperature profiles stay logarithmic at the highest accessible Ra.
Load-bearing premise
The temperature profiles stay logarithmic with a fixed roughness-induced shift even in the ultimate regime.
What would settle it
A simulation or experiment at sufficiently high Ra showing either non-logarithmic temperature profiles or a Nu-Ra exponent approaching 0.5 would contradict the prediction.
Figures
read the original abstract
Heat and momentum transfer in wall-bounded turbulent flow, coupled with the effects of wall-roughness, is one of the outstanding questions in turbulence research. In the standard Rayleigh-B\'enard problem for natural thermal convection, it is notoriously difficult to reach the so-called ultimate regime in which the near-wall boundary layers are turbulent. Following the analyses proposed by Kraichnan [Phys. Fluids vol 5., pp. 1374-1389 (1962)] and Grossmann & Lohse [Phys. Fluids vol. 23, pp. 045108 (2011)], we instead utilize recent direct numerical simulations of forced convection over a rough wall in a minimal channel [MacDonald, Hutchins & Chung, J. Fluid Mech. vol. 861, pp. 138--162 (2019)] to directly study these turbulent boundary layers. We focus on the heat transport (in dimensionless form, the Nusselt number $Nu$) or equivalently the heat transfer coefficient (the Stanton number $C_h$). Extending the analyses of Kraichnan and Grossmann & Lohse, we assume logarithmic temperature profiles with a roughness-induced shift to predict an effective scaling of $Nu \sim Ra^{0.42}$, where $Ra$ is the dimensionless temperature difference, corresponding to $C_h \sim Re^{-0.16}$, where $Re$ is the centerline Reynolds number. This is pronouncedly different from the skin-friction coefficient $C_f$, which in the fully rough turbulent regime is independent of $Re$, due to the dominant pressure drag. In rough-wall turbulence the absence of the analog to pressure drag in the temperature advection equation is the origin for the very different scaling properties of the heat transfer as compared to the momentum transfer. This analysis suggests that, unlike momentum transfer, the asymptotic ultimate regime, where $Nu\sim Ra^{1/2}$, will never be reached for heat transfer at finite $Ra$.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript extends Kraichnan and Grossmann-Lohse analyses of ultimate-regime thermal convection by using forced-convection rough-wall DNS data to assume logarithmic temperature profiles with a fixed roughness-induced shift. This yields the effective scaling Nu ∼ Ra^{0.42} (equivalently C_h ∼ Re^{-0.16}) and the conclusion that, unlike momentum transfer (where C_f becomes Re-independent due to pressure drag), the asymptotic ultimate regime Nu ∼ Ra^{1/2} is never reached for heat transfer at finite Ra because the temperature equation lacks an analog to pressure drag.
Significance. If the assumed logarithmic form with constant shift persists, the work supplies a concrete, assumption-driven extrapolation that explains the distinct scaling behaviors of heat versus momentum transfer in rough-wall turbulence. The approach leverages existing DNS to address a regime that remains inaccessible in natural-convection simulations, and the explicit contrast with pressure-drag effects is a useful conceptual contribution.
major comments (2)
- [Abstract and analysis section] Abstract and analysis section: the exponent 0.42 is obtained by direct substitution of the logarithmic temperature profile (with roughness shift taken from the authors' prior DNS) into the Grossmann-Lohse framework; the numerical value is therefore fixed by the profile assumption rather than by an independent constraint from the natural-convection equations.
- [Analysis and conclusion] Analysis and conclusion: the claim that Nu ∼ Ra^{1/2} is unreachable at finite Ra rests entirely on the premise that the temperature profiles remain logarithmic with a Ra-independent roughness shift even when the boundary layers become fully turbulent. No verification, sensitivity test, or argument is supplied showing why the shift remains constant once the near-wall region is in the ultimate regime.
minor comments (1)
- The manuscript would benefit from an explicit statement of the Reynolds-number range of the underlying DNS data and the precise functional form used for the roughness shift when inserted into the scaling relations.
Simulated Author's Rebuttal
We thank the referee for their detailed review and valuable feedback on our manuscript. Below we provide point-by-point responses to the major comments. We have made revisions to clarify certain aspects of our analysis.
read point-by-point responses
-
Referee: [Abstract and analysis section] Abstract and analysis section: the exponent 0.42 is obtained by direct substitution of the logarithmic temperature profile (with roughness shift taken from the authors' prior DNS) into the Grossmann-Lohse framework; the numerical value is therefore fixed by the profile assumption rather than by an independent constraint from the natural-convection equations.
Authors: We concur that the specific exponent of 0.42 is determined by the substitution of the logarithmic temperature profile and the fixed roughness shift (obtained from our prior DNS) into the Grossmann-Lohse framework. This is the core of our method: to leverage high-fidelity forced-convection data to inform the profile assumptions in the natural-convection analysis. The manuscript does not claim an independent derivation from the natural-convection equations but rather an effective scaling based on this informed assumption. To address this comment, we have revised the abstract and the analysis section to more explicitly state that the exponent follows from the profile assumption. revision: yes
-
Referee: [Analysis and conclusion] Analysis and conclusion: the claim that Nu ∼ Ra^{1/2} is unreachable at finite Ra rests entirely on the premise that the temperature profiles remain logarithmic with a Ra-independent roughness shift even when the boundary layers become fully turbulent. No verification, sensitivity test, or argument is supplied showing why the shift remains constant once the near-wall region is in the ultimate regime.
Authors: The referee is correct that our conclusion depends on the persistence of the logarithmic profile with a constant roughness shift into the ultimate regime. This premise is based on the Reynolds-number independence of the shift observed in the forced-convection DNS over a range of Re in the fully rough regime. We argue that the analogy holds because the near-wall turbulent boundary layers in the ultimate regime of natural convection are expected to behave similarly to those in forced convection. However, we acknowledge the lack of direct verification or sensitivity tests, which is due to the inaccessibility of the ultimate regime in current natural-convection simulations. We have partially revised the conclusion to include a more detailed discussion of this assumption and its justification via the forced-convection analogy, while noting the limitation. revision: partial
- The absence of direct numerical verification or sensitivity analysis for the constancy of the roughness shift in the ultimate regime of Rayleigh-Bénard convection, as this regime is not yet accessible to simulation.
Circularity Check
No significant circularity; scaling follows explicitly from stated assumption
full rationale
The paper states its central premise as an explicit assumption ('we assume logarithmic temperature profiles with a roughness-induced shift') and then performs the algebraic extension of the Kraichnan/GL framework to obtain Nu ∼ Ra^{0.42}. This is a forward derivation from the chosen functional form rather than a self-referential loop in which the output is used to justify the input. Prior DNS data from the authors' own 2019 JFM paper supplies the roughness shift value but is external observational input, not a fitted parameter renamed as a prediction. No equation reduces the claimed exponent to the assumption by construction, and the 'never reaches 1/2' conclusion is presented as a logical consequence of the maintained log-law premise rather than an independent verification. The derivation chain therefore remains self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Temperature profiles are logarithmic with a roughness-induced additive shift even in the ultimate regime
Reference graph
Works this paper leans on
- [1]
-
[2]
D. Lohse and K.-Q. Xia, Small-scale properties of turbulent Rayleigh–B´ enard convection, Annu. Rev. Fluid Mech. 42, 335 (2010)
work page 2010
-
[3]
F. Chilla and J. Schumacher, New perspectives in turbulent Rayleigh–B´ enard convection, Eur. Phys. J. E 35, 58 (2012)
work page 2012
-
[4]
S. Grossmann and D. Lohse, Scaling in thermal convection: A unifying theory, J. Fluid. Mech. 407, 27 (2000)
work page 2000
-
[5]
S. Grossmann and D. Lohse, Thermal convection for large Prandtl number, Phys. Rev. Lett. 86, 3316 (2001)
work page 2001
-
[6]
R. J. A. M. Stevens, E. P. van der Poel, S. Grossmann, and D. Lohse, The unifying theory of scaling in thermal convection: the updated prefactors, J. Fluid Mech. 730, 295 (2013)
work page 2013
-
[7]
R. H. Kraichnan, Turbulent thermal convection at arbitrary Prandtl number, Phys. Fluids 5, 1374 (1962)
work page 1962
-
[8]
E. A. Spiegel, Convection in stars, Annu. Rev. Astron. Astrophys. 9, 323 (1971)
work page 1971
-
[9]
S. Grossmann and D. Lohse, Multiple scaling in the ultimate regime of thermal convection, Phys. Fluids 23, 045108 (2011)
work page 2011
-
[10]
S. Grossmann and D. Lohse, Logarithmic temperature profiles in the ultimate regime of ther- mal convection, Phys. Fluids 24, 125103 (2012)
work page 2012
-
[11]
P. R. Owen and W. R. Thomson, Heat transfer across rough surfaces, J. Fluid Mech. 15, 321 (1963)
work page 1963
-
[12]
L. N. Howard, Heat transport by turbulent convection, J. Fluid. Mech. 17, 405 (1963)
work page 1963
-
[13]
F. H. Busse, On Howard’s upper bound for heat transport in turbulent convection, J. Fluid Mech. 37, 457 (1969)
work page 1969
-
[14]
C. R. Doering and P. Constantin, Variational bounds on energy dissipation in incompressible flows: III. Convection, Phys. Rev. E 53, 5957 (1996)
work page 1996
-
[15]
X. He, D. Funfschilling, H. Nobach, E. Bodenschatz, and G. Ahlers, Transition to the ultimate state of turbulent Rayleigh–B´ enard convection, Phys. Rev. Lett.108, 024502 (2012). 14
work page 2012
-
[16]
X. He, D. Funfschilling, E. Bodenschatz, and G. Ahlers, Heat transport by turbulent Rayleigh– B´ enard convection forPr = 0.8 and 4× 1011 < Ra <2× 1014: ultimate-state transition for aspect ratio Γ = 1.00, New J. Phys. 14, 063030 (2012)
work page 2012
- [17]
-
[18]
X. Chavanne, F. Chilla, B. Castaing, B. Hebral, B. Chabaud, and J. Chaussy, Observation of the ultimate regime in Rayleigh–B´ enard convection, Phys. Rev. Lett.79, 3648 (1997)
work page 1997
-
[19]
X. Chavanne, F. Chilla, B. Chabaud, B. Castaing, and B. Hebral, Turbulent Rayleigh–B´ enard convection in gaseous and liquid He, Phys. Fluids 13, 1300 (2001)
work page 2001
-
[20]
P. E. Roche, G. Gauthier, R. Kaiser, and J. Salort, On the triggering of the ultimate regime of convection, New J. Phys. 12, 085014 (2010)
work page 2010
-
[21]
J. J. Niemela, L. Skrbek, K. R. Sreenivasan, and R. Donnelly, Turbulent convection at very high Rayleigh numbers, Nature 404, 837 (2000)
work page 2000
-
[22]
J. J. Niemela and K. R. Sreenivasan, Does confined turbulent convection ever attain the ‘asymptotic scaling’ with 1/2 power?, New J. Phys. 12, 115002 (2010)
work page 2010
- [23]
-
[24]
X. Zhu, V. Mathai, R. J. A. M. Stevens, R. Verzicco, and D. Lohse, Transition to the ultimate regime in two-dimensional Rayleigh–B´ enard convection, Phys. Rev. Lett.120, 144502 (2018)
work page 2018
-
[25]
D. Krug, X. Zhu, D. Chung, I. Marusic, R. Verzicco, and D. Lohse, Transition to ultimate Rayleigh–B´ enard turbulence revealed through extended self-similarity scaling analysis of the temperature structure functions, J. Fluid Mech. 851, R3 (2018)
work page 2018
-
[26]
Y. Shen, P. Tong, and K.-Q. Xia, Turbulent convection over rough surfaces, Phys. Rev. Lett. 76, 908 (1996)
work page 1996
-
[27]
S. Ciliberto and C. Laroche, Random roughness of boundary increases the turbulent convection scaling exponent, Phys. Rev. Lett. 82, 3998 (1999)
work page 1999
-
[28]
Y. B. Du and P. Tong, Turbulent thermal convection in a cell with ordered rough boundaries, J. Fluid Mech. 407, 57 (2000)
work page 2000
-
[29]
P. E. Roche, B. Castaing, B. Chabaud, and B. Hebral, Observation of the 1/2 power law in 15 Rayleigh–B´ enard convection, Phys. Rev. E63, 045303(R) (2001)
work page 2001
-
[30]
X. L. Qiu, K.-Q. Xia, and P. Tong, Experimental study of velocity boundary layer near a rough conducting surface in turbulent natural convection, J. Turb. 6, 1 (2005)
work page 2005
-
[31]
J. C. Tisserand, M. Creyssels, Y. Gasteuil, H. Pabiou, M. Gibert, B. Castaing, and F. Chilla, Comparison between rough and smooth plates within the same Rayleigh–B´ enard cell, Phys. Fluids 23, 015105 (2011)
work page 2011
- [32]
- [33]
-
[34]
Y.-C. Xie and K.-Q. Xia, Turbulent thermal convection over rough plates with varying rough- ness geometries, J. Fluid Mech. 825, 573 (2017)
work page 2017
-
[35]
M. J. Tummers and M. Steunebrink, Effect of surface roughness on heat transfer in Rayleigh- B´ enard convection, Intl J. Heat Mass Transfer139, 1056 (2019)
work page 2019
-
[36]
G. Stringano and R. Verzicco, Mean flow structure in thermal convection in a cylindrical cell of aspect-ratio one half, J. Fluid Mech. 548, 1 (2006)
work page 2006
-
[37]
O. Shishkina and C. Wagner, Modelling the influence of wall roughness on heat transfer in thermal convection, J. Fluid Mech. 686, 568 (2011)
work page 2011
-
[38]
S. Wagner and O. Shishkina, Heat flux enhancement by regular surface roughness in turbulent thermal convection, J. Fluid Mech. 763, 109 (2015)
work page 2015
-
[39]
X. Zhu, R. J. A. M. Stevens, R. Verzicco, and D. Lohse, Roughness-facilitated local 1/2 scaling does not imply the onset of the ultimate regime of thermal convection, Phys. Rev. Lett. 119, 154501 (2017)
work page 2017
-
[40]
A. J. Smits, B. J. McKeon, and I. Marusic, High-Reynolds Number Wall Turbulence, Annu. Rev. Fluid Mech. 43, 353 (2011)
work page 2011
-
[41]
A. J. Smits and I. Marusic, Wall-bounded turbulence, Phys. Today 66, 25 (2013)
work page 2013
-
[42]
S. Grossmann, D. Lohse, and C. Sun, High Reynolds number Taylor–Couette turbulence, Annu. Rev. Fluid Mech 48, 53 (2016)
work page 2016
-
[43]
F. H. Busse, Viewpoint: The twins of turbulence research, Physics 5, 4 (2012)
work page 2012
-
[44]
B. Eckhardt, S. Grossmann, and D. Lohse, Torque scaling in turbulent Taylor–Couette flow 16 between independently rotating cylinders, J. Fluid Mech. 581, 221 (2007)
work page 2007
- [45]
-
[46]
T. H. van den Berg, C. R. Doering, D. Lohse, and D. P. Lathrop, Smooth and rough boundaries in turbulent Taylor–Couette flow, Phys. Rev. E 68, 036307 (2003)
work page 2003
-
[47]
X. Zhu, R. A. Verschoof, D. Bakhuis, S. G. Huisman, R. Verzicco, C. Sun, and D. Lohse, Wall-roughness induces asymptotic ultimate turbulence, Nature Physics 14, 417 (2018)
work page 2018
-
[48]
Nikuradse, Laws of flow in rough pipes., Forsch
J. Nikuradse, Laws of flow in rough pipes., Forsch. Arb. Ing.-Wes. (1933), English translation published 1950 as NACA Tech. Memo. 1292
work page 1933
-
[49]
F. R. Hama, Boundary-layer characteristics for smooth and rough surfaces, Trans. Soc. Naval Arch. Mar. Engrs 62, 333 (1954)
work page 1954
-
[50]
Jim´ enez, Turbulent flows over rough walls, Annu
J. Jim´ enez, Turbulent flows over rough walls, Annu. Rev. Fluid Mech. 36, 173 (2004)
work page 2004
-
[51]
B. Eckhardt, S. Grossmann, and D. Lohse, Fluxes and energy dissipation in thermal convection and shear flows, Europhys. Lett. 78, 24001 (2007)
work page 2007
-
[52]
M. MacDonald, N. Hutchins, and D. Chung, Roughness effects in turbulent forced convection, J. Fluid Mech. 861, 138 (2019)
work page 2019
-
[53]
L. Chan, M. MacDonald, D. Chung, N. Hutchins, and A. Ooi, A systematic investigation of roughness height and wavelength in turbulent pipe flow in the transitionally rough regime, J. Fluid Mech. 771, 743 (2015)
work page 2015
- [54]
-
[55]
M. MacDonald, D. Chung, N. Hutchins, L. Chan, A. Ooi, and R. Garc´ ıa-Mayoral, The minimal-span channel for rough-wall turbulent flows, J. Fluid Mech. 816, 5 (2017)
work page 2017
-
[56]
B. A. Kader, Temperature and concentration profiles in fully turbulent boundary layers, Intl J. Heat Mass Transfer 24, 1541 (1981)
work page 1981
-
[57]
H. Kawamura, H. Abe, and Y. Matsuo, DNS of turbulent heat transfer in channel flow with respect to Reynolds and Prandtl number effects, Intl J. Heat Fluid Fl. 20, 196 (1999)
work page 1999
-
[58]
S. Pirozzoli, M. Bernardini, and P. Orlandi, Passive scalars in turbulent channel flow at high Reynolds number, J. Fluid Mech. 788, 614 (2016)
work page 2016
-
[59]
A. M. Yaglom, Similarity laws for constant-pressure and pressure-gradient turbulent wall flows, 17 Annu. Rev. Fluid Mech. 11, 505 (1979)
work page 1979
- [60]
-
[61]
S. Leonardi, P. Orlandi, and R. A. Antonia, Heat transfer in a turbulent channel flow with roughness, in Proc. 5th Intl Symp. on Turbulence and Shear Flow Phenomena (Munich, Ger- many, 2007) pp. 785–790
work page 2007
-
[62]
Schlichting, Experimental investigation of the problem of surface roughness, Ing.-Arch
H. Schlichting, Experimental investigation of the problem of surface roughness, Ing.-Arch. 7, 1 (1936), English translation published 1937 as NACA Tech. Memo. 823
work page 1936
-
[63]
M. P. Schultz and K. A. Flack, Turbulent boundary layers on a systematically varied rough wall, Phys. Fluids 21, 015104 (2009)
work page 2009
- [64]
-
[65]
J. P. Monty, E. Dogan, R. Hanson, A. J. Scardino, B. Ganapathisubramani, and N. Hutchins, An assessment of the ship drag penalty arising from light calcareous tubeworm fouling, Bio- fouling 32, 451 (2016)
work page 2016
-
[66]
W. M. Kays, M. E. Crawford, and B. Weigand, Convective Heat and Mass Transfer , 4th ed. (McGraw-Hill, 2005)
work page 2005
-
[67]
O. Shishkina, M. S. Emran, S. Grossmann, and D. Lohse, Scaling relations in large-Prandtl- number natural thermal convection, Phys. Rev. Fluids 2, 103502 (2017)
work page 2017
-
[68]
T. Cebeci and P. Bradshaw, Physical and computational aspects of convective heat transfer (Springer-Verlag, 1984)
work page 1984
- [69]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.