Statistics of Heat Transfer in Two-Dimensional Turbulent Rayleigh-Bénard Convection at Various Prandtl Number.

Statistics of heat transfer in two-dimensional (2D) turbulent Rayleigh-Bénard (RB) convection for Pr=6,20,100 and 106 are investigated using the lattice Boltzmann method (LBM). Our results reveal that the large scale circulation is gradually broken up into small scale structures plumes with the increase of Pr, the large scale circulation disappears with increasing Pr, and a great deal of smaller thermal plumes vertically rise and fall from the bottom to top walls. It is further indicated that vertical motion of various plumes gradually plays main role with increasing Pr. In addition, our analysis also shows that the thermal dissipation is distributed mainly in the position of high temperature gradient, the thermal dissipation rate εθ already increasingly plays a dominant position in the thermal transport, εu can have no effect with increase of Pr. The kinematic viscosity dissipation rate and the thermal dissipation rate gradually decrease with increasing Pr. The energy spectrum significantly decreases with the increase of Pr. A scope of linear scaling arises in the second order velocity structure functions, the temperature structure function and mixed structure function(temperature-velocity). The value of linear scaling and the 2nd-order velocity decrease with increasing Pr, which is qualitatively consistent with the theoretical predictions.

1. Introduction

Thermal convection generally occurs in natural world and industrial field. Hartmann et al. 2001 [1] argued that for weather predictions the flow of atmosphere and thermal convection flow are not only related with smaller length scales and time scales, but also is closely greater scales for climate forecast. Marshall et al. (1999) [2] have investigated a key enforcing mechanical properties of ocean circulation in the ocean. Cardin and Olson (1994) [3] have studied that the enforced convection arises in the earth outer core. In general, the effect of rotation, the changing phases, complex boundary conditions and nonlinear dynamic can play a main role in many thermal convections. Lohse and Xia (2010) [4] and Chilla and Schumacher (2012) [5] respectively introduced research processes for turbulent RB convection. Turbulent RB convection is the most common phenomenon in natural convection [6]. Scheel et al. [7] presented the direct numerical calculations in a cylinder, which mainly focuses on the optimal scales of turbulent convective, especially, the analysis of the kinetic energy rates and thermal dissipation rates in the whole cell [7]. Hu [8] argued that the critical (Rayleigh number) decreases for the beginning of convection, and the aspect ratio effects the critical value with increasing density inversion parameter in RB convection. Zhang [9] have studied the topic of the buoyancy ratio enhancing the fluid stability, and the critical also gradually enhances with increasing value of the buoyant force ratio in thermal convection. Turbulent RB convection [10,11,12,13,14,15,16,17,18,19] has been studied numerically, experimentally and analytically. Shishkina [15] reported that massive thermal plumes rising in turbulent RB convection are analysed on account of 3D numerical calculations, which indicates that the thermal boundary layers are crucial to the temperature fields and velocity fields.

In RB convection, the (a measure of the buoyant force), and the Prandtl number (the kinematic viscosity divides the thermal diffusivity) are imporpant dimensionless numbers. Krishnamurti [20,21] performed extensive convection experiments on mercury (), air (), water (), freon (), and silicon oil (). Busse and Whitehead [22] also reported the jagged instability and cross-roll instability of silicon oil in experiment, which indicated that the results of two-dimensional (2D) convection are closely similar to that of 3D convection for large convection [10,11].

This survey to investigate the bifurcation and chaos in convection with large and different using a low-dimensional model containing is exploited in this paper. The main purpose of our paper is to offer ‘reference’ cases for comparison with a host of various flows. In general, the laboratory experiments are performed by using fluids with up to behave. However, a similar flow for infinite does not exist. Although a model of a plume is not investigated in the Earth’s mantle, the flow possessing overwhelmingly large is close to the mantle investigations since the mantle possesses a of . Our results mainly appraise the variety of physical mechanism phenomenon about the kinematic viscosity dissipation, thermal dissipation, energy spectra, temperature spectra and the 2nd structure function with increasing at the same . Our studies mainly focus on the characteristics of the above physical mechanism phenomena with increasing , the characteristics of flow and scale properties in turbulent RB system at large Prandtl number. In addition, the Bolgiano-Obukhov-like (BO59) scaling is used to explore the profound insight of the velocity fluctuation and the temperature fluctuation with the increase of in turbulent RB convection.

The lattice Boltzmann scheme is performed to simulate in all numerical examples. The lattice Boltzmann method (LBM) had possessed latent advantages in previous studies to multiphase, single, heat and mass transfer hydrodynamic problems [23,24,25,26,27,28,29,30,31,32,33,34,35,36]. The LBM is particularly propitious to attack complex boundary conditions. It is a attractive topic that the turbulent flows are simulated by the LBM [28]. The LBM on account of Direct Numerical Simulations (DNS), Large Eddy Simulations (LES) also possess great success in turbulence predictions [34]. In general, in turbulent RB convection, Rayleigh number () is a very important dimensionless parameter. When is less than , the flow state is laminar RB convection. When is between and , the flow state is soft turbulent RB convection. However, when is greater than , the flow state is full turbulent RB convection.

The main structure of the paper is introduced. Firstly, a thermal lattice Boltzmann model is introduced. Secondly, the numerical verifications are presented. Thirdly, the result analysis and discussions are provided. Finally, some conclusions are provided.

2. Mathematical Equation of Fluid and Lattice Boltzmann Method

The discrete kinetic model is based on LBM. A double-population approach using the lattice Boltzmann equation is proposed by Shan [35]. The lattice Boltzmann equation of computing dynamics and the lattice Boltzmann equation of computing advection equation are described by the following expressions, correspondingly [26,32]: where , denote the probability density functions at , belongs to a discrete velocity, is the mesoscopic buoyant body force, is the flow relaxation time and is the temperature relaxation time in the above equations, respectively. The local kinetic equilibrium is equilibrium function for flow, and is the equilibrium function for temperature [32]. where is the sound speed and denotes the weight coefficients.

The relation between the relaxation parameter and the kinematic viscosity is presented and the relaxation is also presented between the parameter and the fluid thermal diffusivity .

The mesoscopic buoyant body force with the Boussinesq approximation is formulated by the following equation in which g denotes the gravity acceleration, denotes the coefficient of thermal expansion, and denotes the y-component of correspondingly.

The relations are defined as coarse-grained (in velocity space) fields of the distribution functions for the macroscopical density, momentum, and temperature, which are given by

The Equations (1) and (2) are respectively expansed by a Chapman-Enskogm approach, which leads to the macroscopical equations of fluid. The inertial terms are reproduced by the streaming step on the left-hand side in the hydrodynamical equations, whereas the kinematic viscosity coefficient and the coefficient of thermal diffusion are respectively connected to the relaxation (towards equilibrium) properties in the right hand side.

The is an important variable in the natural convection. The definition of is in which is the kinematic viscosity, is the thermal diffusivity and is the coefficient of isobaric thermal expansion, is the temperature difference between the bottom boundary and top boundary in cavity, g is the gravity acceleration, and H is the height of cavity.

The perturbed is given in the initial stage of simulation. Once RB convection is taken shape, the heat transfer near wall is overwhelmingly strengthened. The Nusselt number denotes the enhancement of the heat transfer. The expression of is presented in our simulations. in which denotes the vertical velocity and is the spatial average in the whole flow domain.

3. Calculation Results and Discussions

In the present study, numerical calculations in 2D turbulent RB convection for and are carried out at . The computing grid is set to . The no-slip boundary condition is performed on top boundary and on the bottom boundary. And left and right boundary conditions are also performed by the no-slip boundary condition in all numerical calculations. The initial dimensionless temperature in bottom boundary is 1, and the initial dimensionless temperature of top boundary is 0, respectively. Meanwhile, a linear distribution of the dimensionless temperature from 0 to 1 is performed between top boundary and bottom boundary.

The budget relation of instantaneous kinetic energy is used to test and verify the accuracy of the numerical calculations. where is the total potential energy, denotes the total kinetic energy and denotes the dissipation rate of total kinetic energy. The ratio of the right-hand side to the left-hand side of Equation (14) as a relation of the normalized time is presented in Figure 1. One can see that the difference within only at all times is achieved at the energy balance equation, which is indicated that the accuracy of LBM is reliable.

Figure 1

Validation of the energy balance relation (14) for the LBM scheme.

3.1. Global Quantities of Turbulent RB Convection

In Figure 2, the temperature distribution with superimposed temperature fields at , (from a to d), and is displayed for several elective coarseness types, where is the flow characteristic of RB convection system, the blue regions denote cold fluid, and the red regions correspond to hot fluid. One can see that large-scale circulations are obtained in the center of the cavity at . A large number of thermal plumes are predominantly rising in corners of the cavity, and the cold plumes are falling in corners of the cavity, which can be main reason owing to a large scale circulation that is guided along one of diagonals of the cavity. It is demonstrated that large-scale circulations of the fluid in cavity are decomposed gradually with the increase of , which is consistent with the experiment for large [4]. It is also showed that vertical motion of various plumes gradually plays main role with increasing . Respectively, a large number of small-scale thermal plumes vertically rise from the bottom wall to top wall, a large number of small-scale cold plumes vertically fall from the top boundary to the bottom boundary, and a great number of small eddies appear at . The longitudinal velocity play a predominant role in cavity with increasing . It is further found that the average time of achieving steady-state increase overwhelmingly with increasing .

Figure 2

Temperature distributions with steady flow at , (a–d), and .

In the following section, Grossmann -Lohse theory for is used to investigate the dissipation rate of kinetic-energy rate and the thermal dissipation rate, which is the following equations [6]

The relations of Equations (15) and (16) are able to be derived from the Boussinesq equations [6].

The distributions of kinetic-energy dissipation rate at , (from a to d), and are shown in Figure 3. It is shown that the value of kinetic-energy dissipation rate decreases with the increase of in cavity. The kinetic-energy dissipation rate at plays a major role around cavity. The kinetic-energy dissipation rate is occupied gradually in the interior of cavity with increasing . The kinetic-energy dissipation rate becomes gradually dispersed at .

Figure 3

Kinetic energy dissipation rates , (a–d), and .

The distributions of thermal dissipation rate at , (from a to d) are shown in Figure 4, respectively. One can see that the value of thermal dissipation rate also decreases with the increase of in cavity. The thermal dissipation rate at plays a major role around cavity, similarly. The thermal dissipation rate gradually expands in the interior of cavity with increasing . A large number of thermal dissipation emerge in the whole domain at . From a detailed comparison of Figure 3 with Figure 4, we can see that the thermal dissipation is distributed mainly in the position of high temperature gradient. In other word, the intense thermal dissipation events mainly concentrate on the interfaces of hot and cold fluids and these interfaces between the hot fluid and the cold fluid develop to complex associate structures with large intricate.

Figure 4

Thermal energy dissipation rates at , (a–d), and .

Figure 5 and Figure 6 display the mean vertical dissipation rate profiles of kinetic-energy and the mean vertical thermal dissipation rate profiles for , , respectively, where represents a horizontal average. As observed in Figure 5 and Figure 6, the mean vertical dissipation rate profiles of kinetic-energy and the mean vertical thermal dissipation rate profiles significantly decrease with increasing . The strong dissipation events mainly arises near the top and bottom boundary in the turbulent range. It is further seen that the thermal dissipation rate already plays a key role in the thermal energy transport, thus may be neglected in virtue of inverse cascade of kinetic energy in thermal convection.

Figure 5

Mean vertical kinetic-energy dissipation rate profiles.

Figure 6

Mean vertical thermal dissipation rate profiles.

The contrast of the temperature on the midline (y = H/2) at , , and is displayed in Figure 7. Plotted in Figure 7, it is noted that the relative magnitude of temperature difference is increased overwhelmingly with increasing , and the variety of temperature is also enhanced drastically.

Figure 7

Instantaneous temperature distributions with superimposed temperature fields at = 6, 20, 100, , and .

The Bolgiano-Obukhov scaling is used to judge the dissipation scale in 2D turbulence thermal convection. The rigorous relations are obtained by assuming spatial homogeneity, within the Boussinesq approximations equation [4].

The time and volume -averaged kinetic dissipation rates versus the Rayleigh number , the time and volume-averaged thermal dissipation rates versus the and Nusselt number as a function of at = 6, 20, 100, and between numerical results and theoretical values are given in Figure 8, Figure 9 and Figure 10. Figure 8 represents the so-called volume and time-averaged kinetic dissipation rates versus the Rayleigh number . Figure 9 represents the thermal dissipation rates versus the in the so-called volume and time-averaged. The solid lines represent the theoretical relation as relations (17) and (18). The dispersed points are the results of the present LBM. It is mainly demonstrated that the numerical calculation results of the LBM are excellently consistent with the theoretical value for time-averaged kinetic dissipation rates and thermal dissipation rates versus the in Figure 8 and Figure 9.

Figure 8

Volume and time-averaged kinetic dissipation rates versus the Rayleigh number .

Figure 9

Volume and time-averaged thermal dissipation rates versus the Rayleigh number .

Figure 10

Nusselt number as a function of Rayleigh number at , and .

The relation of turbulent heat flux and the relation of the kinetic energy will be discussed, respectively. Figure 10 demonstrates the Nusselt number as a relation of Rayleigh number at , and . Figure 11 manifests the Reynolds number () as a relation of Rayleigh number, respectively. Nusselt is wall averaged Nusselt number in the whole computational domain. The Grossmann-Lohse theory had offered a good insight for the of comprehension of , (See Equations (19) and (20)) and even allowed multifarious predictions where the bulk turbulence plays an important role both the holistic kinetic dissipation and the system thermal dissipation [6].

Figure 11

Reynolds number as a function of Rayleigh number at , and .

The relations (19) and (20) respectively represent the steady state regime of RB turbulence, which was proposed for turbulent RB convection at overwhelmingly high by Kraichnan, and then developed by Grossmann and Lohse in their relations [6]. The ultimate state scaling was observed numerically [37] and experimentally [38,39], where the ultimate regime scaling is achieved when the solid boundaries are lacked. Plotted in Figure 10 and Figure 11, it is noted that a linear scaling can be presented for for both and for nearly four decades from to . It is validated that the energy spectrum significantly decreases with the increase of the . Our results further demonstrate that the numerical calculation results of the LBM are well consistent with the Grossmann-Lohse relations for and [6] from = to in Figure 10 and Figure 11.

3.2. Scaling of Energy Spectra, Fluxes and Spatial Intermittency

The scaling of energy spectra and fluxes will be discussed in the following section. For passive scalars, the value is in agreement with Obukhov-Corrsin scaling [6], and consistent with experimental findings of Kolmogorov scaling along the center line of turbulent convection [9]. Temperature power spectra have been investigated previously by some researchers [37,38]. and in which denotes a single dimensionless relation of a dimensionless argument. The kinetic spectra and temperature spectra in central point(Nx/2,Ny/2) at time scales are showed in Figure 12 and Figure 13, respectively. These results demonstrate that the energy spectrum is significantly decreased by the increase of the in Figure 12. For the dimensionless inertial range (), the numerical results are consistent with Kolmogorov’s 5/3 spectrum at . For , the numerical results well agree with Kolmogorov’s 5/3 spectrum at the dimensionless inertial range (). For the dimensionless inertial range (), the numerical results are consistent with Kolmogorov’s 5/3 spectrum at . For , excellent agreement between the numerical results and Kolmogorov’s 5/3 spectrum is presented at the dimensionless inertial range (). It is found that the dimensionless inertial range is enlarged with increasing , which is demonstrated that the energy dissipation rate scaling rises () and the thermal dissipation rate () increases gradually with increasing . It is also demonstrated that current results could be explained near at low and near at lightly higher frequencies, respectively, agree with the collapse well on top of each other in the Bolgiano-Obukhov-like scaling [4,6]. To investigate the spatial intermittency effects of turbulent RB convection for different Prandtl numbers, the longitudinal velocity (See Equation (19))and the structure function of temperature (See Equation (20))over horizontal separations are regarded. The 2nd-order structure functions of velocity fluctuations and the 2nd-order structure functions of temperature fluctuations are hoped to next expressions, which are the most-studied quantities in buoyancy-driven turbulence [4].

Figure 12

Kinetic energy spectra at different .

Figure 13

Temperature spectra at different .

It is found that tht value of linear scaling decreases with the increase of . in which △ denotes the difference between the top wall and the bottom wall, and H denotes height of cavity. These relations can be easily derived from the Boussinesq equations and the corresponding boundary conditions [6]. Secondly, the mixed velocity-temperature structure function is expected to following expression

Figure 14, Figure 15 and Figure 16 demonstrate respectively the second order structure function of velocity , structure function of temperature and the mixed structure function (velocity-temperature) at the whole computational domain computed at a late stage of the self-similar regime. As can be seen from Figure 14, as the increases, the 2nd-order velocity decreases, which is qualitatively consistent with the theoretical predictions. It is validated that the 2nd-order velocity structure functions display a series of linear scaling that extends to larger scales for different Prandtl numbers, i.e, for , for , for and for in Figure 14. It is validated that the value of linear scaling decreases with the increase of , which may be given rise to large-scale circulation disappearing, the secondary eddies arising, and the dissipation of energy playing a major role at in cavity, respectively. Noting that as the increases, the 2nd-order velocity decreases, which is qualitatively consistent with the theoretical predictions (See Equation (19)). It is also demonstrated that a range of linear scaling emerges in the 2nd-order structure functions of velocity fluctuations, the 2nd-order structure functions fluctuations of temperature and the 2nd-order structure functions of mixed structure function, which may be in virtue of the small-scale temperature structures, like the thermal spikes or plumes that possess strong corrections. Plotted in Figure 15 and Figure 16, the 2nd-order structure functions of temperature and the mixed structure function are also qualitative in agreement with those observed in Figure 14.

Figure 14

2nd-order velocity structure functions at different .

Figure 15

2nd-order temperature structure functions at different .

Figure 16

2nd-order mixed velocity-temperature structure function at different .

4. Conclusions

Numerical calculations of two-dimensional turbulent RB convection at large Prandtl number are investigated using LBM. Our results demonstrate that the variation characteristics of physical mechanism phenomenon about the kinematic viscosity dissipation, thermal dissipation, energy spectra, temperature spectra and the 2nd structure function with increasing at the same . The main conclusions are as follows:

Fisrtly, the large scale circulation in cavity is gradually broken up into small scale structures that thermal plumes rise and fall from the bottom to top walls with increasing , which agrees with the visualization of the experiment. The large-scale circulations of the fluid are decomposed gradually with the increase of . It is found that the vertical motion of various plumes gradually play a main role with increasing . Especially, a large number of small-scale thermal plumes vertically rise from the bottom wall to top wall, a large number of small-scale cold plumes vertically fall from the top wall to the bottom wall, and a great deal of small eddies appear at . Secondly, the kinematic viscosity dissipation rate and the thermal dissipation rate gradually decrease with increasing , which is qualitatively consistent with the theoretical predictions. The intense thermal dissipation events mainly concentrate on the interfaces of hot fluid and cold fluid. The thermal dissipation rate already plays a key role in the thermal energy transport, thus can have no effect in virtue of inverse cascade of kinetic energy in thermal convection. It is also demonstrated that the numerical calculation results of the LBM are excellently consistent with the theoretical value for time-averaged kinetic and thermal dissipation rates versus the . In addition, the energy spectrum significantly decreases with increasing , which indicates that the numerical results of the LBM agree with the Grossmann-Lohse theory for and [6] from to . Finally, the 2nd-order velocity decreases with increasing , which is qualitatively consistent with the theoretical value. It is also found that a scope of linear scaling appears in the second order structure functions of velocity fluctuations, the second order structure functions of temperature fluctuations and mixed structure function, which may be in virtue of the small-scale temperature structures, like the thermal plumes or spikes.

Due to natural parallelism of the LBM, the LBM may have potential prospects to simulate turbulent RB convection. Even though the results of two-dimensional cases are reported in this paper, similar scaling laws may appear (the Kolmogorov-Obukhov scenario) in 3D and it might be crucial to modify parametric assembles in numerical simulations. Accordingly, allowed us to evaluate the importance of the Kolmogorov-Obukhov scenario in three dimensions at large . Our results discussed in this paper provide insights into the flow dynamics of RB convection. These results will be useful in modelling convective flows in the atmospheres and interiors of stars and planets, as well as in the applications of engineer.