On non-Oberbeck-Boussinesq effects in three-dimensional Rayleigh-Bénard convection in glycerol

Susanne Horn, Olga Shishkina, Claus Wagner

Research output: Contribution to journalArticle

28 Citations (Scopus)
9 Downloads (Pure)

Abstract


Rayleigh–Bénard convection in glycerol (Prandtl number ) in a cylindrical cell with an aspect ratio of was studied by means of three-dimensional direct numerical simulations (DNS). For that purpose, we implemented temperature-dependent material properties into our DNS code, by prescribing polynomial functions up to seventh order for the viscosity, the heat conductivity and the density. We performed simulations with the common Oberbeck–Boussinesq (OB) approximation and with non-Oberbeck–Boussinesq (NOB) effects within a range of Rayleigh numbers of . For the highest temperature differences, , the viscosity at the top is times higher than at the bottom, while the differences of the other material properties are less than . We analysed the temperature and velocity profiles and the thermal and viscous boundary-layer thicknesses. NOB effects generally lead to a breakdown of the top–bottom symmetry, typical for OB Rayleigh–Bénard convection. Under NOB conditions, the temperature in the centre of the cell increases with increasing and can be up to higher than under OB conditions. The comparison of our findings with several theoretical and empirical models showed that two-dimensional boundary-layer models overestimate the actual , while models based on the temperature or velocity scales predict very well with a standard deviation of . Furthermore, the obtained temperature profiles bend closer towards the cold top plate and further away from the hot bottom plate. The situation for the velocity profiles is reversed: they bend farther away from the top plate and closer towards to the bottom plate. The top boundary layers are always thicker than the bottom ones. Their ratio is up to 2.5 for the thermal and up to 4.5 for the viscous boundary layers. In addition, the Reynolds number and the Nusselt number were investigated: is higher and is lower under NOB conditions. The Nusselt number is influenced in a nonlinear way by NOB effects, stronger than was suggested by the two-dimensional simulations. The actual scaling of with in the NOB case is and is in excellent agreement with the experimental data.
Original languageEnglish
Pages (from-to)175-202
Number of pages28
JournalJournal of Fluid Mechanics
Volume724
Early online date29 Apr 2013
DOIs
Publication statusPublished - 10 Jun 2013
Externally publishedYes

Fingerprint

glycerols
Glycerol
convection
Boundary layers
Nusselt number
direct numerical simulation
temperature profiles
boundary layers
velocity distribution
two dimensional boundary layer
Direct numerical simulation
viscosity
Temperature
boundary layer thickness
temperature scales
Materials properties
Rayleigh number
Prandtl number
cells
Viscosity

Cite this

On non-Oberbeck-Boussinesq effects in three-dimensional Rayleigh-Bénard convection in glycerol. / Horn, Susanne; Shishkina, Olga; Wagner, Claus .

In: Journal of Fluid Mechanics, Vol. 724, 10.06.2013, p. 175-202.

Research output: Contribution to journalArticle

@article{6fd08b30e2d84e18b5352484926f4fc4,
title = "On non-Oberbeck-Boussinesq effects in three-dimensional Rayleigh-B{\'e}nard convection in glycerol",
abstract = "Rayleigh–B{\'e}nard convection in glycerol (Prandtl number ) in a cylindrical cell with an aspect ratio of was studied by means of three-dimensional direct numerical simulations (DNS). For that purpose, we implemented temperature-dependent material properties into our DNS code, by prescribing polynomial functions up to seventh order for the viscosity, the heat conductivity and the density. We performed simulations with the common Oberbeck–Boussinesq (OB) approximation and with non-Oberbeck–Boussinesq (NOB) effects within a range of Rayleigh numbers of . For the highest temperature differences, , the viscosity at the top is times higher than at the bottom, while the differences of the other material properties are less than . We analysed the temperature and velocity profiles and the thermal and viscous boundary-layer thicknesses. NOB effects generally lead to a breakdown of the top–bottom symmetry, typical for OB Rayleigh–B{\'e}nard convection. Under NOB conditions, the temperature in the centre of the cell increases with increasing and can be up to higher than under OB conditions. The comparison of our findings with several theoretical and empirical models showed that two-dimensional boundary-layer models overestimate the actual , while models based on the temperature or velocity scales predict very well with a standard deviation of . Furthermore, the obtained temperature profiles bend closer towards the cold top plate and further away from the hot bottom plate. The situation for the velocity profiles is reversed: they bend farther away from the top plate and closer towards to the bottom plate. The top boundary layers are always thicker than the bottom ones. Their ratio is up to 2.5 for the thermal and up to 4.5 for the viscous boundary layers. In addition, the Reynolds number and the Nusselt number were investigated: is higher and is lower under NOB conditions. The Nusselt number is influenced in a nonlinear way by NOB effects, stronger than was suggested by the two-dimensional simulations. The actual scaling of with in the NOB case is and is in excellent agreement with the experimental data.",
author = "Susanne Horn and Olga Shishkina and Claus Wagner",
year = "2013",
month = "6",
day = "10",
doi = "10.1017/jfm.2013.151",
language = "English",
volume = "724",
pages = "175--202",
journal = "Journal of Fluid Mechanics",
issn = "0022-1120",
publisher = "Cambridge University Press",

}

TY - JOUR

T1 - On non-Oberbeck-Boussinesq effects in three-dimensional Rayleigh-Bénard convection in glycerol

AU - Horn, Susanne

AU - Shishkina, Olga

AU - Wagner, Claus

PY - 2013/6/10

Y1 - 2013/6/10

N2 - Rayleigh–Bénard convection in glycerol (Prandtl number ) in a cylindrical cell with an aspect ratio of was studied by means of three-dimensional direct numerical simulations (DNS). For that purpose, we implemented temperature-dependent material properties into our DNS code, by prescribing polynomial functions up to seventh order for the viscosity, the heat conductivity and the density. We performed simulations with the common Oberbeck–Boussinesq (OB) approximation and with non-Oberbeck–Boussinesq (NOB) effects within a range of Rayleigh numbers of . For the highest temperature differences, , the viscosity at the top is times higher than at the bottom, while the differences of the other material properties are less than . We analysed the temperature and velocity profiles and the thermal and viscous boundary-layer thicknesses. NOB effects generally lead to a breakdown of the top–bottom symmetry, typical for OB Rayleigh–Bénard convection. Under NOB conditions, the temperature in the centre of the cell increases with increasing and can be up to higher than under OB conditions. The comparison of our findings with several theoretical and empirical models showed that two-dimensional boundary-layer models overestimate the actual , while models based on the temperature or velocity scales predict very well with a standard deviation of . Furthermore, the obtained temperature profiles bend closer towards the cold top plate and further away from the hot bottom plate. The situation for the velocity profiles is reversed: they bend farther away from the top plate and closer towards to the bottom plate. The top boundary layers are always thicker than the bottom ones. Their ratio is up to 2.5 for the thermal and up to 4.5 for the viscous boundary layers. In addition, the Reynolds number and the Nusselt number were investigated: is higher and is lower under NOB conditions. The Nusselt number is influenced in a nonlinear way by NOB effects, stronger than was suggested by the two-dimensional simulations. The actual scaling of with in the NOB case is and is in excellent agreement with the experimental data.

AB - Rayleigh–Bénard convection in glycerol (Prandtl number ) in a cylindrical cell with an aspect ratio of was studied by means of three-dimensional direct numerical simulations (DNS). For that purpose, we implemented temperature-dependent material properties into our DNS code, by prescribing polynomial functions up to seventh order for the viscosity, the heat conductivity and the density. We performed simulations with the common Oberbeck–Boussinesq (OB) approximation and with non-Oberbeck–Boussinesq (NOB) effects within a range of Rayleigh numbers of . For the highest temperature differences, , the viscosity at the top is times higher than at the bottom, while the differences of the other material properties are less than . We analysed the temperature and velocity profiles and the thermal and viscous boundary-layer thicknesses. NOB effects generally lead to a breakdown of the top–bottom symmetry, typical for OB Rayleigh–Bénard convection. Under NOB conditions, the temperature in the centre of the cell increases with increasing and can be up to higher than under OB conditions. The comparison of our findings with several theoretical and empirical models showed that two-dimensional boundary-layer models overestimate the actual , while models based on the temperature or velocity scales predict very well with a standard deviation of . Furthermore, the obtained temperature profiles bend closer towards the cold top plate and further away from the hot bottom plate. The situation for the velocity profiles is reversed: they bend farther away from the top plate and closer towards to the bottom plate. The top boundary layers are always thicker than the bottom ones. Their ratio is up to 2.5 for the thermal and up to 4.5 for the viscous boundary layers. In addition, the Reynolds number and the Nusselt number were investigated: is higher and is lower under NOB conditions. The Nusselt number is influenced in a nonlinear way by NOB effects, stronger than was suggested by the two-dimensional simulations. The actual scaling of with in the NOB case is and is in excellent agreement with the experimental data.

U2 - 10.1017/jfm.2013.151

DO - 10.1017/jfm.2013.151

M3 - Article

VL - 724

SP - 175

EP - 202

JO - Journal of Fluid Mechanics

JF - Journal of Fluid Mechanics

SN - 0022-1120

ER -