DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Classical 1/3 scaling of convection holds up to Ra = 10 15

Abstract

The global transport of heat and momentum in turbulent convection is constrained by thin thermal and viscous boundary layers at the heated and cooled boundaries of the system. This bottleneck is thought to be lifted once the boundary layers themselves become fully turbulent at very high values of the Rayleigh number Ra —the dimensionless parameter that describes the vigor of convective turbulence. Laboratory experiments in cylindrical cells for Ra1 0 12 have reported different outcomes on the putative heat transport law. Here we show, by direct numerical simulations of three-dimensional turbulent Rayleigh–Bénard convection flows in a slender cylindrical cell of aspect ratio 1/10 , that the Nusselt number—the dimensionless measure of heat transport—follows the classical power law of Nu= ( 0.0525 ± 0.006 ) × R a 0.331 ± 0.002 up to Ra=1 0 15 . Intermittent fluctuations in the wall stress, a blueprint of turbulence in the vicinity of the boundaries, manifest at all Ra considered here, increasing with increasing Ra , and suggest that an abrupt transition of the boundary layer to turbulence does not take place.

Authors:
 [1]; ORCiD logo [2]; ORCiD logo [3];  [4]
  1. Tandon School of Engineering, New York University, New York, NY 11201,
  2. Department of Physics, Occidental College, Los Angeles, CA 90041,
  3. Tandon School of Engineering, New York University, New York, NY 11201,, Department of Mechanical Engineering, Technische Universität Ilmenau, D-98684 Ilmenau, Germany,
  4. Tandon School of Engineering, New York University, New York, NY 11201,, Courant Institute of Mathematical Sciences, New York University, New York, NY 10012,, Department of Physics, New York University, New York, NY 10012,, Center for Space Science, New York University Abu Dhabi, Abu Dhabi 129188, United Arab Emirates
Publication Date:
Research Org.:
Argonne National Laboratory (ANL), Argonne, IL (United States)
Sponsoring Org.:
USDOE
OSTI Identifier:
1606457
Alternate Identifier(s):
OSTI ID: 1625066
Grant/Contract Number:  
AC02_06CH11357; AC02-06CH11357
Resource Type:
Published Article
Journal Name:
Proceedings of the National Academy of Sciences of the United States of America
Additional Journal Information:
Journal Name: Proceedings of the National Academy of Sciences of the United States of America Journal Volume: 117 Journal Issue: 14; Journal ID: ISSN 0027-8424
Publisher:
Proceedings of the National Academy of Sciences
Country of Publication:
United States
Language:
English
Subject:
Science & Technology - Other Topics

Citation Formats

Iyer, Kartik P., Scheel, Janet D., Schumacher, Jörg, and Sreenivasan, Katepalli R. Classical 1/3 scaling of convection holds up to Ra = 10 15. United States: N. p., 2020. Web. doi:10.1073/pnas.1922794117.
Iyer, Kartik P., Scheel, Janet D., Schumacher, Jörg, & Sreenivasan, Katepalli R. Classical 1/3 scaling of convection holds up to Ra = 10 15. United States. https://doi.org/10.1073/pnas.1922794117
Iyer, Kartik P., Scheel, Janet D., Schumacher, Jörg, and Sreenivasan, Katepalli R. Wed . "Classical 1/3 scaling of convection holds up to Ra = 10 15". United States. https://doi.org/10.1073/pnas.1922794117.
@article{osti_1606457,
title = {Classical 1/3 scaling of convection holds up to Ra = 10 15},
author = {Iyer, Kartik P. and Scheel, Janet D. and Schumacher, Jörg and Sreenivasan, Katepalli R.},
abstractNote = {The global transport of heat and momentum in turbulent convection is constrained by thin thermal and viscous boundary layers at the heated and cooled boundaries of the system. This bottleneck is thought to be lifted once the boundary layers themselves become fully turbulent at very high values of the Rayleigh numberRa—the dimensionless parameter that describes the vigor of convective turbulence. Laboratory experiments in cylindrical cells forRa≳1012have reported different outcomes on the putative heat transport law. Here we show, by direct numerical simulations of three-dimensional turbulent Rayleigh–Bénard convection flows in a slender cylindrical cell of aspect ratio1/10, that the Nusselt number—the dimensionless measure of heat transport—follows the classical power law ofNu=(0.0525±0.006)×Ra0.331±0.002up toRa=1015. Intermittent fluctuations in the wall stress, a blueprint of turbulence in the vicinity of the boundaries, manifest at allRaconsidered here, increasing with increasingRa, and suggest that an abrupt transition of the boundary layer to turbulence does not take place.},
doi = {10.1073/pnas.1922794117},
journal = {Proceedings of the National Academy of Sciences of the United States of America},
number = 14,
volume = 117,
place = {United States},
year = {Wed Mar 25 00:00:00 EDT 2020},
month = {Wed Mar 25 00:00:00 EDT 2020}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record
https://doi.org/10.1073/pnas.1922794117

Citation Metrics:
Cited by: 55 works
Citation information provided by
Web of Science

Figures / Tables:

Fig. 1 Fig. 1: Global scaling laws of turbulent trans- A B C port. (A) Turbulent heat transport law Nu(Ra) for two datasets on log–log coordinates. The present data (open circles) are for Pr = 1 and Γ=0.1 for 108 ≤ Ra ≤1015. Open squares representing DNS at Pr=0.7 and Γ=1 formore » 3 x 105 ≤ Ra ≤1010 are taken from ref. 23. Also added is the power-law fit of Nu= (0.0525 ± 0.006) x Ra0.331±0.002 which is obtained for 1010 ≤ Ra ≤1015. (B) Linear–log plot of the data of A, displayed in the form of compensated classical power law of 1/3 slope, with error bars computed from the standard deviation of the time series Nu(t), obtained at $\tilde{z}$= z/H=0 (bottom) and $\tilde{z}$=1 (top). The plot shows that the power-law exponent is smaller (0.29 ± 0.01) for Ra ≤ 109 and is 1/3 in the high-Ra range, with no tendency to a larger slope as would be required of the possible progression to the ultimate state. The dashed line corresponds to prefactor 0.0525, in close agreement with the classical prediction of 0.073 (6, 7). (C) Compensated power-law plot of the turbulent momentum transport law Re(Ra). A fit over the same Rayleigh number range as in B gives Re=(0.1555 ± 0.006) x Ra0.458 ± 0.006. Symbols mean the same in A–C.« less

Save / Share:

Works referenced in this record:

Efficiency of Heat Transfer in Turbulent Rayleigh-Bénard Convection
journal, July 2011


Convective Self-Aggregation in Numerical Simulations: A Review
journal, February 2017

  • Wing, Allison A.; Emanuel, Kerry; Holloway, Christopher E.
  • Surveys in Geophysics, Vol. 38, Issue 6
  • DOI: 10.1007/s10712-017-9408-4

Turbulent convection at very high Rayleigh numbers
journal, April 2000

  • Niemela, J. J.; Skrbek, L.; Sreenivasan, K. R.
  • Nature, Vol. 404, Issue 6780
  • DOI: 10.1038/35009036

Large granulation cells on the surface of the giant star π1 Gruis
journal, December 2017


Scaling of hard thermal turbulence in Rayleigh-Bénard convection
journal, July 1989

  • Castaing, Bernard; Gunaratne, Gemunu; Heslot, François
  • Journal of Fluid Mechanics, Vol. 204, Issue -1
  • DOI: 10.1017/S0022112089001643

Zhu et al. Reply:
journal, December 2019


Resolving the fine-scale structure in turbulent Rayleigh–Bénard convection
journal, November 2013


Irminger Sea deep convection injects oxygen and anthropogenic carbon to the ocean interior
journal, October 2016

  • Fröb, F.; Olsen, A.; Våge, K.
  • Nature Communications, Vol. 7, Issue 1
  • DOI: 10.1038/ncomms13244

New perspectives in turbulent Rayleigh-Bénard convection
journal, July 2012


Turbulent rotating convection confined in a slender cylinder: The sidewall circulation
journal, February 2020


Lithium–antimony–lead liquid metal battery for grid-level energy storage
journal, September 2014


An Overlapping Schwarz Method for Spectral Element Solution of the Incompressible Navier–Stokes Equations
journal, May 1997


Turbulent Thermal Convection at Arbitrary Prandtl Number
journal, January 1962


Transition to the Ultimate State of Turbulent Rayleigh-Bénard Convection
journal, January 2012


Wall-bounded turbulent flows at high Reynolds numbers: Recent advances and key issues
journal, June 2010

  • Marusic, I.; McKeon, B. J.; Monkewitz, P. A.
  • Physics of Fluids, Vol. 22, Issue 6
  • DOI: 10.1063/1.3453711

Influence of 2D and 3D convection–diffusion flow on tritium permeation in helium cooled solid breeder blanket units
journal, February 2006


High Rayleigh Number Convection
journal, January 1994


Confined turbulent convection
journal, April 2003


Absence of Evidence for the Ultimate Regime in Two-Dimensional Rayleigh-Bénard Convection
journal, December 2019


Turbulent Rayleigh–Bénard convection in gaseous and liquid He
journal, May 2001

  • Chavanne, X.; Chillà, F.; Chabaud, B.
  • Physics of Fluids, Vol. 13, Issue 5
  • DOI: 10.1063/1.1355683

Heat transfer and large scale dynamics in turbulent Rayleigh-Bénard convection
journal, April 2009

  • Ahlers, Guenter; Grossmann, Siegfried; Lohse, Detlef
  • Reviews of Modern Physics, Vol. 81, Issue 2
  • DOI: 10.1103/RevModPhys.81.503

Convection in Stars I. Basic Boussinesq Convection
journal, September 1971


Figures/Tables have been extracted from DOE-funded journal article accepted manuscripts.