Transition to turbulence scaling in Rayleigh-Bénard convection

Jörg Schumacher, Ambrish Pandey, Victor Yakhot, Katepalli Sreenivasan

Research output: Contribution to journalArticle

Abstract

If a fluid flow is driven by a weak Gaussian random force, the nonlinearity in the Navier-Stokes equations is negligibly small and the resulting velocity field obeys Gaussian statistics. Nonlinear effects become important as the driving becomes stronger and a transition occurs to turbulence with anomalous scaling of velocity increments and derivatives. This process has been described by Yakhot and Donzis [Phys. Rev. Lett. 119, 044501 (2017)PRLTAO0031-900710.1103/PhysRevLett.119.044501] for homogeneous and isotropic turbulence. In more realistic flows driven by complex physical phenomena, such as instabilities and nonlocal forces, the initial state itself, and the transition to turbulence from that initial state, is much more complex. In this paper, we discuss the Reynolds-number dependence of moments of the kinetic energy dissipation rate of orders 2 and 3 obtained in the bulk of thermal convection in the Rayleigh-Bénard system. The data are obtained from three-dimensional spectral element direct numerical simulations in a cell with square cross section and aspect ratio 25 by Pandey et al. [Nat. Commun. 9, 2118 (2018)2041-172310.1038/s41467-018-04478-0]. Different Reynolds numbers 1 Re 1000 which are based on the thickness of the bulk region and the corresponding root-mean-square velocity are obtained by varying the Prandtl number Pr from 0.005 to 100 at a fixed Rayleigh number Ra=105. A few specific features of the data agree with the theory. The normalized moments of the kinetic energy dissipation rate En show a nonmonotonic dependence for small Reynolds numbers before obeying the algebraic scaling prediction for the turbulent state. Implications and reasons for this behavior are discussed.

Original languageEnglish (US)
Article number033120
JournalPhysical Review E
Volume98
Issue number3
DOIs
StatePublished - Sep 28 2018

Fingerprint

Transition to Turbulence
Rayleigh
Convection
Reynolds number
convection
turbulence
Scaling
Energy Dissipation
Kinetic energy
scaling
energy dissipation
kinetic energy
Moment
moments
Anomalous Scaling
Spectral Elements
Thermal Convection
homogeneous turbulence
isotropic turbulence
Rayleigh number

ASJC Scopus subject areas

  • Statistical and Nonlinear Physics
  • Statistics and Probability
  • Condensed Matter Physics

Cite this

Transition to turbulence scaling in Rayleigh-Bénard convection. / Schumacher, Jörg; Pandey, Ambrish; Yakhot, Victor; Sreenivasan, Katepalli.

In: Physical Review E, Vol. 98, No. 3, 033120, 28.09.2018.

Research output: Contribution to journalArticle

Schumacher, Jörg ; Pandey, Ambrish ; Yakhot, Victor ; Sreenivasan, Katepalli. / Transition to turbulence scaling in Rayleigh-Bénard convection. In: Physical Review E. 2018 ; Vol. 98, No. 3.
@article{cc7c00c27f104c0aa736e1603085cf8a,
title = "Transition to turbulence scaling in Rayleigh-B{\'e}nard convection",
abstract = "If a fluid flow is driven by a weak Gaussian random force, the nonlinearity in the Navier-Stokes equations is negligibly small and the resulting velocity field obeys Gaussian statistics. Nonlinear effects become important as the driving becomes stronger and a transition occurs to turbulence with anomalous scaling of velocity increments and derivatives. This process has been described by Yakhot and Donzis [Phys. Rev. Lett. 119, 044501 (2017)PRLTAO0031-900710.1103/PhysRevLett.119.044501] for homogeneous and isotropic turbulence. In more realistic flows driven by complex physical phenomena, such as instabilities and nonlocal forces, the initial state itself, and the transition to turbulence from that initial state, is much more complex. In this paper, we discuss the Reynolds-number dependence of moments of the kinetic energy dissipation rate of orders 2 and 3 obtained in the bulk of thermal convection in the Rayleigh-B{\'e}nard system. The data are obtained from three-dimensional spectral element direct numerical simulations in a cell with square cross section and aspect ratio 25 by Pandey et al. [Nat. Commun. 9, 2118 (2018)2041-172310.1038/s41467-018-04478-0]. Different Reynolds numbers 1 Re 1000 which are based on the thickness of the bulk region and the corresponding root-mean-square velocity are obtained by varying the Prandtl number Pr from 0.005 to 100 at a fixed Rayleigh number Ra=105. A few specific features of the data agree with the theory. The normalized moments of the kinetic energy dissipation rate En show a nonmonotonic dependence for small Reynolds numbers before obeying the algebraic scaling prediction for the turbulent state. Implications and reasons for this behavior are discussed.",
author = "J{\"o}rg Schumacher and Ambrish Pandey and Victor Yakhot and Katepalli Sreenivasan",
year = "2018",
month = "9",
day = "28",
doi = "10.1103/PhysRevE.98.033120",
language = "English (US)",
volume = "98",
journal = "Physical Review E - Statistical Physics, Plasmas, Fluids, and Related Interdisciplinary Topics",
issn = "1063-651X",
publisher = "American Physical Society",
number = "3",

}

TY - JOUR

T1 - Transition to turbulence scaling in Rayleigh-Bénard convection

AU - Schumacher, Jörg

AU - Pandey, Ambrish

AU - Yakhot, Victor

AU - Sreenivasan, Katepalli

PY - 2018/9/28

Y1 - 2018/9/28

N2 - If a fluid flow is driven by a weak Gaussian random force, the nonlinearity in the Navier-Stokes equations is negligibly small and the resulting velocity field obeys Gaussian statistics. Nonlinear effects become important as the driving becomes stronger and a transition occurs to turbulence with anomalous scaling of velocity increments and derivatives. This process has been described by Yakhot and Donzis [Phys. Rev. Lett. 119, 044501 (2017)PRLTAO0031-900710.1103/PhysRevLett.119.044501] for homogeneous and isotropic turbulence. In more realistic flows driven by complex physical phenomena, such as instabilities and nonlocal forces, the initial state itself, and the transition to turbulence from that initial state, is much more complex. In this paper, we discuss the Reynolds-number dependence of moments of the kinetic energy dissipation rate of orders 2 and 3 obtained in the bulk of thermal convection in the Rayleigh-Bénard system. The data are obtained from three-dimensional spectral element direct numerical simulations in a cell with square cross section and aspect ratio 25 by Pandey et al. [Nat. Commun. 9, 2118 (2018)2041-172310.1038/s41467-018-04478-0]. Different Reynolds numbers 1 Re 1000 which are based on the thickness of the bulk region and the corresponding root-mean-square velocity are obtained by varying the Prandtl number Pr from 0.005 to 100 at a fixed Rayleigh number Ra=105. A few specific features of the data agree with the theory. The normalized moments of the kinetic energy dissipation rate En show a nonmonotonic dependence for small Reynolds numbers before obeying the algebraic scaling prediction for the turbulent state. Implications and reasons for this behavior are discussed.

AB - If a fluid flow is driven by a weak Gaussian random force, the nonlinearity in the Navier-Stokes equations is negligibly small and the resulting velocity field obeys Gaussian statistics. Nonlinear effects become important as the driving becomes stronger and a transition occurs to turbulence with anomalous scaling of velocity increments and derivatives. This process has been described by Yakhot and Donzis [Phys. Rev. Lett. 119, 044501 (2017)PRLTAO0031-900710.1103/PhysRevLett.119.044501] for homogeneous and isotropic turbulence. In more realistic flows driven by complex physical phenomena, such as instabilities and nonlocal forces, the initial state itself, and the transition to turbulence from that initial state, is much more complex. In this paper, we discuss the Reynolds-number dependence of moments of the kinetic energy dissipation rate of orders 2 and 3 obtained in the bulk of thermal convection in the Rayleigh-Bénard system. The data are obtained from three-dimensional spectral element direct numerical simulations in a cell with square cross section and aspect ratio 25 by Pandey et al. [Nat. Commun. 9, 2118 (2018)2041-172310.1038/s41467-018-04478-0]. Different Reynolds numbers 1 Re 1000 which are based on the thickness of the bulk region and the corresponding root-mean-square velocity are obtained by varying the Prandtl number Pr from 0.005 to 100 at a fixed Rayleigh number Ra=105. A few specific features of the data agree with the theory. The normalized moments of the kinetic energy dissipation rate En show a nonmonotonic dependence for small Reynolds numbers before obeying the algebraic scaling prediction for the turbulent state. Implications and reasons for this behavior are discussed.

UR - http://www.scopus.com/inward/record.url?scp=85054546122&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85054546122&partnerID=8YFLogxK

U2 - 10.1103/PhysRevE.98.033120

DO - 10.1103/PhysRevE.98.033120

M3 - Article

AN - SCOPUS:85054546122

VL - 98

JO - Physical Review E - Statistical Physics, Plasmas, Fluids, and Related Interdisciplinary Topics

JF - Physical Review E - Statistical Physics, Plasmas, Fluids, and Related Interdisciplinary Topics

SN - 1063-651X

IS - 3

M1 - 033120

ER -