Table of Contents

Thermal Boundary Layers in Turbulent Rayleigh-Bénard Convection with Rough and Smooth Surfaces

Experiments by du Puits


Description

In this data set, we provide wall-normal profiles of the mean temperature \(\overline{T}(z)\) (\(z\) is the coordinate perpendicular to the wall) and its fluctuations \(T'(z)\) in turbulent Rayleigh-Bénard convection with rough and smooth walls. The data has been obtained in the large-scale Rayleigh-Bénard experiment “Barrel of Ilmenau (BOI)”. We used very tiny micro-thermistors of 150 \(\mu\)m diameter and 350 \(\mu\)m length to measure the temperature in and outside the thermal boundary layer at the heated bottom plate with, up to now, unsurpassed spatial and temporal resolution. Our measurements covered Rayleigh numbers between \(Ra=5.4\times 10^9\) and \(Ra<9.6\times 10^{11}\) and two aspect ratios \(\Gamma=1.1\) and \(\Gamma=2.9\). The Prandtl number was fixed at \(Pr=0.7\). The definition of these numbers is as follows:

\[ Ra = \frac{\beta g \Delta T H^3}{\nu \kappa} \qquad\qquad Pr = \frac{\nu}{\kappa} \qquad\qquad \Gamma=\frac{D}{H}. \]

\(\Delta T\) is the temperature difference between the heated bottom and the cooled top plate. \(H\) and \(D\) are their distance and the diameter of the test cell. \(\beta\), \(\nu\), and \(\kappa\) are the thermal expansion coefficient, the kinematic viscosity and the thermal diffusivity of the fluid. \(g\) is the gravitational constant of the Earth. The mean temperatures \(\overline{T}\) and the fluctuations were calculated from 90 min time series with a sampling rate of 50 samples/s. A more detailed description of the experimental facility, the measurement technique and a table with a selection of the parameters will be provided in the subsequent section. Moreover, we refer to a recently published paper in Physical Review Fluids (du Puits, 2024).

Experimental set-up, measurement technique and main parameters

All measurements have been undertaken in the large-scale RB experiment “Barrel of Ilmenau” (figure 1) (see also Cases 90, 91 and 94). The experiment consists of an adiabatic container of 8.0 m in height and 7.1 m in diameter which is filled with air (\(Pr=0.7\)). At the bottom of the container, a heating plate composed of an electrical underfloor heating and a water-flown plate of Aluminum on top heats the air. The cooling plate at the top of the container consists of 16 segments with water circulation inside as well. The segments are mounted on a solid steel construction and are separately leveled perpendicular to the vector of gravity. The entire construction with a weight of about 6 tons is mounted on a crane and can be lifted up and down. The sandwich technology of the heating and the cooling plates guarantees a uniformity of the temperature at their surfaces with a deviation of the local temperature less than 1 per cent (1.5 per cent for the cooling plate) of the total temperature drop \(\Delta T\) between the plates. Separate PID-controllers keep the temperature of both plates in a band of ±0.02 K around the value adjusted. The sidewall of the container is shielded by an active compensation heating system. Eighteen electrical heating elements are arranged between an inner and an outer isolation of 16 and 12 cm thickness, respectively. Their temperatures are controlled to be equal to the temperatures at the inner surface of the wall with the consequence of virtually no heat flux across the sidewall. More detailed information on the experimental facility can be found in du Puits et al (2013).

We placed 336 cuboid roughness elements of 30 mm by 30 mm cross section and 12 mm height at the surface of the heating plate. They are made of aluminum and form a chessboard-like pattern with a period of 2d = 60 mm (figure 1(b)). The cuboids are mounted at a circular aluminum disk, 950 mm in diameter and 3 mm thick, using adhesive pads with high thermal conductivity. This “rough” disk was placed at the center of the heating plate adding a 0.15-mm-thick foil of silicon rubber in the gap between the surface of the plate and the disk. The foil significantly improves the thermal contact and ensures that the temperature drop between the plate and the disk remains as small as possible. For the measurements at the smooth surface, we removed the disk temporarily.

Experimental set-up sketchFig. 1: (a) Principal setup of the large-scale convection facility ``Barrel of Ilmenau''. (b) Detailed sketch of a fraction of the rough surface with the temperature sensors. The sketch shows the elements in the correct proportion except the temperature sensors, which are up-scaled by a factor of 2 for better visibility.

We used very tiny thermistors to measure the temperature field near the rough surface. They have an ellipsoidal shape. Their diameter and their length amount to 150 \(\mu\)m and 350 \(\mu\)m, respectively. They are fixed at a holder which can be moved in steps of 10 \(\mu\)m up and down. The temperature at the top of one cuboid and in one valley in between the cuboids was measured using two of them in an arrangement shown in figure 1(b). The resistance of the thermistors was measured using a PC based data acquisition system. The data rate amounts to 50 samples/s and is sufficiently high to capture even the fastest temperature fluctuations. Before we used the thermistors in the BOI, we calibrated them against a certified PT-100 in a calibration chamber. The remaining uncertainty of ±20 mK was further reduced by an in-situ calibration in the measurement chamber.

Our experiments spanned two and a half decades in Rayleigh number \(Ra=5.4\times 10^9\) to \(Ra<9.6\times 10^{11}\) and two aspect ratios \(\Gamma=1.1\) and \(\Gamma=2.9\). To vary the Rayleigh number within a single aspect ratio, we changed the temperature difference \(\Delta T\). The temperature of the heating and the cooling plates was always set symmetrically above and below a temperature of 30.0oC, except for the measurements at \(\Gamma=1.1\), \(Ra=9.6\times 10^{11}\) and \(\Gamma=2.9\), \(Ra=5.9\times 10^{10}\). For those measurements, the theoretical center temperature amounted to 42.5oC. The main parameters of all measurements are listed in the table below, where the Nusselt number \(Nu\) results from direct measurements of the local wall heat flux at the smooth plate, and the boundary layer thickness \(\delta_{th}=H/(2Nu)\) is given as a reference in comparison to the roughness height \(h\).

\(\Gamma\) \(Ra\) \(\Delta T\) [K] \(T_b\) [oC] \(H\) [m] \(Nu\) \(\delta_{th}\) [mm]
1.1 \(5.8\times 10^{10}\) 2.7 30.4 6.3 164.2 19.18
1.1 \(8.6\times 10^{10}\) 4.0 30.5 6.3 256.8 12.26
1.1 \(1.5\times 10^{11}\) 7.0 30.6 6.3 313.9 10.04
1.1 \(2.2\times 10^{11}\) 10.0 30.8 6.3 361.6 8.71
1.1 \(4.3\times 10^{11}\) 20.0 31.2 6.3 429.2 7.34
1.1 \(6.3\times 10^{11}\) 30.0 31.6 6.3 471.5 6.68
1.1 \(8.0\times 10^{11}\) 38.0 31.8 6.3 542.8 5.80
1.1 \(9.6\times 10^{11}\) 55.0 45.0 6.3 559.0 5.64
2.9 \(5.4\times 10^{9}\) 4.0 30.3 2.5 71.7 17.44
2.9 \(9.4\times 10^{9}\) 7.0 30.5 2.5 90.4 13.82
2.9 \(1.3\times 10^{10}\) 10.0 30.6 2.5 106.4 11.75
2.9 \(2.7\times 10^{10}\) 20.0 31.1 2.5 127.1 9.83
2.9 \(4.0\times 10^{10}\) 30.0 31.5 2.5 143.4 8.72
2.9 \(5.1\times 10^{10}\) 38.0 30.5 2.5 160.6 7.78
2.9 \(5.9\times 10^{10}\) 54.7 43.7 2.5 165.0 7.58

In addition, we mounted heat flux plates to measure the local heat flux released and absorbed from the surfaces of the hot bottom and the cold top plates, respectively. The sensors, developed and provided by the company Phymeas GbR are little discs with a well known heat conductivity. They have a total diameter of 33 mm, however, only the inner 25 mm area is being active for the measurement. We installed five such sensors each in the bottom and the top plates of our experiment. One is located at the center of each plate and the other four are equally spaced at a radius of \(R=0.45\) m. They are flush-mounted in a notch in the surface of the bottom (top) plate to prevent measurement errors caused by artificial turbulence at the edges of the sensors.

Available data

The data set incorporates wall-normal profiles of the normalized mean temperature \(\overline{\Theta}(z/H)\) and its fluctuations \(\Theta'(z/H)\):

\[ \overline{\Theta}=\frac{\overline{T}(z)-\overline{T}_b}{\overline{T}_{hp}-\overline{T}_b} \qquad\qquad \Theta'=\sqrt{\frac{1}{N}\sum_{i=1}^{N}(\Theta_i-\overline{\Theta})^2} \]

\(z/H\) is the coordinate perpendicular to the wall normalized by the plate distance.

Sample plots of the data are available.

The data files can be downloaded as compressed archives, or individually from the table below.

Data files contain profiles of mean and rms normalized temperature obtained above the top of a cuboid, above the valley in between the cuboids, and from a smooth surface as a reference, for each aspect ratio and Rayleigh number. The Excel file contains the same data, with a sheet corresponding to each specified Rayleigh number and aspect ratio. The underlying time series can be requested from the author.

\(\Gamma\) \(Ra\) File
1.1 \(5.8\times 10^{10}\) rbc-g11-ra58e10.dat
1.1 \(8.6\times 10^{10}\) rbc-g11-ra86e10.dat
1.1 \(1.5\times 10^{11}\) rbc-g11-ra15e11.dat
1.1 \(2.2\times 10^{11}\) rbc-g11-ra22e11.dat
1.1 \(4.3\times 10^{11}\) rbc-g11-ra43e11.dat
1.1 \(6.3\times 10^{11}\) rbc-g11-ra63e11.dat
1.1 \(8.0\times 10^{11}\) rbc-g11-ra80e11.dat
1.1 \(9.6\times 10^{11}\) rbc-g11-ra96e11.dat
2.9 \(5.4\times 10^{9}\) rbc-g29-ra54e09.dat
2.9 \(9.4\times 10^{9}\) rbc-g29-ra94e09.dat
2.9 \(1.3\times 10^{10}\) rbc-g29-ra13e10.dat
2.9 \(2.7\times 10^{10}\) rbc-g29-ra27e10.dat
2.9 \(4.0\times 10^{10}\) rbc-g29-ra40e10.dat
2.9 \(5.1\times 10^{10}\) rbc-g29-ra51e10.dat
2.9 \(5.9\times 10^{10}\) rbc-g29-ra59e10.dat
Excel file with all data rbc-temperature-data.xlsx

References

  1. du Puits, R., Resagk, C., Thess, A. (2013), Thermal boundary layers in turbulent Rayleigh–Bénard convection at aspect ratios between 1 and 9, New J. Physics, Vol. 15, p. 013040.

Acknowledgment

The authors wish to acknowledge the funding by the Deutsche Forschungsgemeinschaft (DFG,German Research Foundation) under Grant No. 495678007. Moreover, they wish to acknowledge Sabine Scherge for her technical assistance running the experiments.


Indexed data:

case095 (dbcase, confined_flow)
case095
titleThermal boundary layers in turbulent Rayleigh-Benard convection with rough and smooth surfaces
authordu Puits
year2024
typeEXP
flow_tag3d, scalar, buoyant