Time Evolution Features of Entropy Generation Rate in Turbulent Rayleigh-Bénard Convection with Mixed Insulating and Conducting Boundary Conditions

Time evolution features of kinetic and thermal entropy generation rates in turbulent Rayleigh-Bénard (RB) convection with mixed insulating and conducting boundary conditions at Ra = 109 are numerically investigated using the lattice Boltzmann method. The state of flow gradually develops from laminar flow to full turbulent thermal convection motion, and further evolves from full turbulent thermal convection to dissipation flow in the process of turbulent energy transfer. It was seen that the viscous, thermal, and total entropy generation rates gradually increase in wide range of t/τ < 32 with temporal evolution. However, the viscous, thermal, and total entropy generation rates evidently decrease at time t/τ = 64 compared to that of early time. The probability density function distributions, spatial-temporal features of the viscous, thermal, and total entropy generation rates in the closed system provide significant physical insight into the process of the energy injection, the kinetic energy, the kinetic energy transfer, the thermal energy transfer, the viscous dissipated flow and thermal dissipation.


Introduction
The Rayleigh-Bénard (RB) convection is one of most classical natural convections [1][2][3][4][5], which widely occur in a range of natural and industrial applications [1,2], such as in the Earth's core and mantle, atmosphere, oceans and stars, nuclear reactors, crystallization processes, solar heating devices and so on. The RB convection has been extensively investigated by several experiments in the last few decades [6,7], mostly in slender cells of aspect ratio smaller than or equal to unity in order to reach the largest possible Rayleigh numbers (Ra) or to reveal the detailed characteristic mechanisms of turbulent viscous dynamic and heat transport near the walls or central domains [3,4]. The detailed dynamical and statistical insights of the included turbulent transport and their coherent structures also have been increasingly studied in detail by direct numerical simulations (DNS) [8,9].
It is well known that when the DNS of involving no parametrization of subgrid-scale is carried out, all the dynamically important scales are resolved to faithfully represent the flow. Bailon et al. [10] and Schell et al. [11] reported the derived resolution criteria, thus starting the pioneering work [10,11]; subsequent refinements of this criterion were studied [12,13] and the fine boundary layer dynamics resolution were the main focus [14]. It is only recently that the focus of DNS investigations was presented by the bulk of the research in a convection cell with detailed discussions of the scaling statistics and properties of the dissipation fields [15].
In order to understand the global flow and heat transport loss mechanisms in detail, we will present the detailed statistical characteristic mechanisms of the velocity and temperature gradient fields, in particular the related viscous and thermal local entropy generation rates. The viscous and thermal local entropy generation rates as a criterion are used to provide insight into the local viscous and thermal flow loss in the flow field [16][17][18][19][20][21][22][23][24][25][26][27][28][29][30]. The viscous and thermal components of the local entropy generation rate are derived in the two-dimensional Cartesian space [30]. Their expressions are as follows, respectively [20,28]. .
where µ is the dynamics viscosity of fluid, k denotes the thermal conductivity, θ is the temperature, u represents the x-direction velocity and v is the y-direction velocity, respectively. The total entropy generation rate is the summation of the viscous and thermal entropy generation rates, its expression is as follows [22]: The Bejan number (Be) is regarded as an effective approach to judge the importance of heat transfer irreversibility in the domain [23]. Rejane et al. proposed the contribution of heat transfer entropy generation on over all entropy generation by using the Be [27]. The Be is defined by the following equation [29]: The range of Be is from 0 to 1. When Be is equal to 0, the irreversibility is dominated by fluid friction. Correspondingly, the irreversibility is dominated by heat transfer when the Be is equal to 1. The irreversibility due to heat transfer dominates in the flow when the Be is greater than 1/2. Correspondingly, Be < 1/2 implies that the irreversibilities due to the viscous effects dominate the processes. Meanwhile, it is also noted that the heat transfer and fluid friction entropy generation are equal in Be = 0.5 [27].
A wide variety of thermal plumes caused by the buoyancy in turbulent RB convection with mixed insulating and conducting boundary conditions play dominant role in the heat transfer. Once time evolution of the heat transfer has been described and understood in classical turbulent RB convection, the time evolution features of kinetic and thermal entropy generation rates in turbulent RB convection with mixed insulating and conducting boundary conditions still be further expanded. How does the mixed insulating and conducting boundary conditions reduce or improve the time evolution features of entropy generation rate? The mixed insulating and conducting boundary conditions considerably affect the time evolution characteristics of thermal plumes and entropy generation rate.
Based on above discussions, our work mainly focuses on the effect of the mixed insulating and conducting boundary conditions on the time evolution features of thermal plumes, the viscous, thermal and total entropy generation rates. The physical insight features of kinetic and thermal entropy generation rates with time evolution are discussed in in turbulent RB convection with the mixed insulating and conducting boundary conditions, which mainly tried to understand the dynamics of fluid. The remainder of this paper is divided into the following parts. In Section 2, the thermal fluid dynamics equations and numerical method will be briefly depicted. In Section 3, the detailed results of numerical simulation and some discussions are presented. Finally, some conclusions are addressed.

Convection Diffusion Equation and Numerical Method
In this section, the convection diffusion equation of thermal fluid, and numerical method of solving them are represented, respectively.

Convection Diffusion Equation of Thermal Fluid
The convection diffusion equation of thermal fluid is the classical Oberbeck-Boussinesq equations [2,3], and their expressions can be given as follows [6,8].
in which ρ is the fluid density, u represents the macroscopic velocity, ν is the kinematic viscosity, p is the pressure of fluid, κ denotes the diffusivity, β is the thermal diffusion coefficient, g is the force of gravity, ∆θ represents the difference of temperature and θ denotes macroscopic temperature of fluid. A large number of numerical methods are widely used to solve the classical Oberbeck-Boussinesq equations [31][32][33][34][35][36]. The finite element methods [36], finite difference method [34] and the finite volume method [35] are traditional macroscopic methods for Computational Fluid Dynamics (CFD) calculation. The lattice Boltzmann method (LBM) is a computational fluid dynamics method based on mesoscopic simulation scale [37][38][39][40][41]. Compared with other traditional CFD calculation methods, this method has mesoscopic model characteristics between a micromolecular dynamics model and a macrocontinuous model. LBM also has the advantage of a simple description of fluid interaction, and is easier to set a complex boundary, reach a parallel calculation, and implement a program, etc. [42]. LBM has been widely considered as an effective method to describe fluid motion and deal with engineering problems [43,44]. In the subsection, double distribution LBM for will be introduced.

Numerical Method for Convection Diffusion Equation of Thermal Fluid
In the present paper, the double distributions of LBM are implemented to study the convection diffusion equation of thermal fluid, respectively [38,39]. A lattice Boltzmann equation is implemented to simulate the fluid flow field. Its expression is as follows [42]: where f i (x, t) denotes the density distribution functions at (x, t), c i represents the discrete velocity. F i is the discrete force term, f eq i (x, t) is the equilibrium function of density distribution, and τ ν denotes the relaxation time. The equilibrium function for the density is given as: where w i denotes the weight coefficient [42]. The kinematic viscosity ν is computed by the following equation The lattice Boltzmann equation for the temperature field is given by the following equation where g i (x, t) is the temperature distribution function at (x, t), τ θ denotes the relaxation times for temperature evolution equation. g eq i is the equilibrium function for temperature distribution. Its expression is given as [40]: The diffusivity κ is as follows: The density, macroscopic velocity, and temperature are as follows: The Chapman-Enskog expansions of Equations (8) and (11) are used to derive the classical Oberbeck-Boussinesq equations [42]. A macroscopic length scale (x 1 = εx) and two macroscopic time scales (t 1 = εt, t 2 = εt) are implemented in the Chapman-Enskog expansion. Two time scales ∂t = ε∂ t1 + ε 2 ∂ t2 and one spatial scale ∂ x = ε∂ α are used for the Frisch, Hasslacher, and Pomeau (FHP) model [38]. The classical Oberbeck-Boussinesq equations can be derived by executing the streaming step and using the above Chapman-Enskog expansion [42]. Figure 1 shows the computational model of geometrical schematic. As shown in Figure 1, the inhomogeneities heat plates are restricted only in the bottom condition (y = 0), and are made of alternating regions of either the isothermal boundary condition, θ = θ down , where the discrete black region denotes heat source, or adiabatic boundary condition, ∂ y θ = 0. The upper boundary keeps a constant temperature, θ = θ up . In this physical model, the ratio of width of dividing height, ξ = H/L, and another two nondimensionless parameters are implemented to define the geometrical configuration of the discrete heat source; the ratio of single heat source is defined as λ = l/L, and the total ratio of discrete heat source area, η = nl/L, in which n denotes the heat source number and l represents the single heat-source length. When η is equal to 1, the model becomes the classical RB convection. The above several nondimensionless parameters are implemented to obtain a better understanding of heat transfer transport in turbulent Rayleigh-Bénard convection with mixed insulating and conducting boundary conditions. The Rayleigh number is one of the most important dimensionless parameter in RB convection. Its expression is as follows: The Nusselt number is one of the most important dimensionless parameter in RB convection to reflect the performance of heat transfer system. It is obtained by: The Rayleigh number is one of the most important dimensionless parameter in RB convection. Its expression is as follows: Entropy 2020, 22, 672

of 16
The Nusselt number is one of the most important dimensionless parameter in RB convection to reflect the performance of heat transfer system. It is obtained by: (16) where ∆θ is the temperature difference between the top boundary and the bottom boundary, H denotes the channel height, u y is the y-velocity, and · represents the average value of entire domain. The boundary conditions play a key role in the computational stability and accuracy. The periodic boundary condition and approach of nonequilibrium extrapolation are carried out in this paper. Their ideas will be introduced, respectively. The idea of the periodic boundary condition approach is as follows [42]: in which the vector L represents length of the flow pattern. The approach of nonequilibrium extrapolation is as follows [42]: (20) in which the nonequilibrium contribution can be derived from the fluid node x f next to x b along the boundary normal vector. During propagation, the unknown incoming populations can be obtained by leaving the domain at the opposite side. As illustrated in Figure 1, the inhomogeneities heat plates are implemented in all numerical simulations. For Ra = 10 9 , 4000 × 2000 lattices in two-dimensional space are implemented to study the temperature fields, viscous, thermal and total entropy generation rates. The parameter λ is equal to 1/9, η is equal to 5/9, the nonequilibrium extrapolation is applied at the top and bottom boundary conditions, the periodic boundary condition is used at left and right boundaries, and the Prandtl number (Pr = Entropy 2020, 20, x 6 of 17 number (Pr = ѵ ) is equal to 1. The dimensionless temperature of discrete heat source equals to 300 in Figure 1, and the dimensionless initial temperature of the fluid is 299.

Results and Discussions
The analysis of temperature field, flow streamlines, and various entropy generation rates will be discussed with spatial-temporal evolution in this section, respectively. According to the development phenomenon of thermal plumes at times t/τ = 8 and t/τ = 16, the thermal convective motion of the whole field is still in the initial stage of turbulent development. It was seen that with time evolution, a large-scale thermal plume ascends, strikes on the top plate, and in-volutes several thermal plumes to both sides in the left half of the system, and two large-scale thermal plumes descend at t/τ = 32. These thermal plumes interact and strike on the top and bottom plates with time evolution, a number of small-scale thermal plumes appear at t/τ = 32, which demonstrates that the physical system of thermal convection gradually evolves from the large-scale to small-scale thermal plumes with time evolution [3]. In the process of energy cascade of turbulent thermal convection, the energy of the first large vortex comes from the thermal buoyancy of the outside world, which produces the second small vortex. After the small vortex loses its stability, it produces a smaller vortex process [5]. At t/τ = 64, some smallest plumes can be coagulated into the big plumes, several big plumes reappear with time evolution again. The above phenomenon of temperature distributions with time evolution is consistent with the previous studies [11].

Analysis of Flow and Temperature Field
/κ) is equal to 1. The dimensionless temperature of discrete heat source equals to 300 in Figure 1, and the dimensionless initial temperature of the fluid is 299.

Results and Discussions
The analysis of temperature field, flow streamlines, and various entropy generation rates will be discussed with spatial-temporal evolution in this section, respectively. Figure 2 describes the isotherms' temperature distributions with time evolution at t/τ = 8, t/τ = 16, t/τ = 32, and t/τ = 64. Here, τ (τ = H/βg∆θ) is the characteristic time of the computing system. As described in Figure 2, it can be seen that a few thermal plumes ascend in the region of dimensionless bottom boundary (0.5), two big thermal plumes descend in the region of dimensionless top boundary (1.5) at time t/τ = 8, and the large-scale thermal plumes descend in the region of dimensionless top boundary (0.5) and ascend in the region of dimensionless bottom boundary (1.5) at time t/τ = 16. According to the development phenomenon of thermal plumes at times t/τ = 8 and t/τ = 16, the thermal convective motion of the whole field is still in the initial stage of turbulent development. It was seen that with time evolution, a large-scale thermal plume ascends, strikes on the top plate, and in-volutes several thermal plumes to both sides in the left half of the system, and two large-scale thermal plumes descend at t/τ = 32. These thermal plumes interact and strike on the top and bottom plates with time evolution, a number of small-scale thermal plumes appear at t/τ = 32, which demonstrates that the physical system of thermal convection gradually evolves from the large-scale to small-scale thermal plumes with time evolution [3]. In the process of energy cascade of turbulent thermal convection, the energy Entropy 2020, 22, 672 6 of 16 of the first large vortex comes from the thermal buoyancy of the outside world, which produces the second small vortex. After the small vortex loses its stability, it produces a smaller vortex process [5]. At t/τ = 64, some smallest plumes can be coagulated into the big plumes, several big plumes reappear with time evolution again. The above phenomenon of temperature distributions with time evolution is consistent with the previous studies [11].

Analysis of Flow and Temperature Field
was seen that with time evolution, a large-scale thermal plume ascends, strikes on the top plate, and in-volutes several thermal plumes to both sides in the left half of the system, and two large-scale thermal plumes descend at t/τ = 32. These thermal plumes interact and strike on the top and bottom plates with time evolution, a number of small-scale thermal plumes appear at t/τ = 32, which demonstrates that the physical system of thermal convection gradually evolves from the large-scale to small-scale thermal plumes with time evolution [3]. In the process of energy cascade of turbulent thermal convection, the energy of the first large vortex comes from the thermal buoyancy of the outside world, which produces the second small vortex. After the small vortex loses its stability, it produces a smaller vortex process [5]. At t/τ = 64, some smallest plumes can be coagulated into the big plumes, several big plumes reappear with time evolution again. The above phenomenon of temperature distributions with time evolution is consistent with the previous studies [11].  To further demonstrate the above thermal convection flow phenomenon of the whole field, the streamlines of the thermal convection flow at four same time evolution steps are shown in Figure 3. As illustrated in Figure 3, one can clearly see that two large vortexes occur in the central region of the whole field due to the injection of energy at the early characteristic time; two small vortexes appear at dimensionless bottom boundary (0.5) and dimensionless top boundary (1.5) at time t/τ = 8 respectively. In addition, two small vortexes ascend in the region of the dimensionless bottom boundary (0.5), two small vortexes descend in the region of dimensionless top boundary (1.5). At time t/τ = 16, two large vortexes in central region of the whole field become unstable, and more small vortexes appear in the dimensionless top and bottom boundaries (0.5 and 1.5). It was seen that at time t/τ = 32, the early large vortexes evolve into a large number of small scale vortexes, and a large number of small scale vortexes generate due to energy transfer process in the whole field. However, many small vortexes disappear in main flow field, and two relatively big vortexes reappear t/τ = 64. The above phenomenon demonstrates that several large vortexes interact and develop to a large number of small vortexes and a few small vortexes dissipate and big vortexes reappear with temporal evolution, which qualitatively depicts that the state of flow gradually develops from laminar flow to full turbulent thermal convection motion, and further evolve from full turbulent thermal convection to dissipation flow in the process of turbulent energy transfer. above phenomenon demonstrates that several large vortexes interact and develop to a large number of small vortexes and a few small vortexes dissipate and big vortexes reappear with temporal evolution, which qualitatively depicts that the state of flow gradually develops from laminar flow to full turbulent thermal convection motion, and further evolve from full turbulent thermal convection to dissipation flow in the process of turbulent energy transfer.

Analysis of Entropy Generation Rate
The isotherms temperature distributions and streamlines with time evolution are represented in the above section and several analyses of the entropy generation rate will be discussed in the following section. Figure 4 describes the viscous entropy generation rate at times t/τ = 8, t/τ = 16, t/τ = 32, and = t/τ = 64. As shown in Figure 4, it is clearly seen that the high viscous entropy generation rate mainly appears in the intersectional region between the main flow and the top and bottom boundaries and in the intersectional region between big vortexes at times t/τ = 8 and t/τ = 16. It was seen that the high viscous entropy generation rate mainly appears in the high shear region between main flow region and vortex, the low viscous entropy generation rate occurs near the central region

Analysis of Entropy Generation Rate
The isotherms temperature distributions and streamlines with time evolution are represented in the above section and several analyses of the entropy generation rate will be discussed in the following section. Figure 4 describes the viscous entropy generation rate at times t/τ = 8, t/τ = 16, t/τ = 32, and t/τ = 64. As shown in Figure 4, it is clearly seen that the high viscous entropy generation rate mainly appears in the intersectional region between the main flow and the top and bottom boundaries and in the intersectional region between big vortexes at times t/τ = 8 and t/τ = 16. It was seen that the high viscous entropy generation rate mainly appears in the high shear region between main flow region and vortex, the low viscous entropy generation rate occurs near the central region of various vortex, which indicates that the viscous flow loss mainly occurs in the high shear region. Meanwhile, at time step t/τ = 32, the viscous entropy generation rate evidently increases with temporal evolution. Nevertheless, the viscous entropy generation rate evidently decreases at time t/τ = 64 compared to that of t/τ = 32, which indicates that the whole mainstream field has already entered the state of turbulent dissipation. Figure 5 illustrates the thermal entropy generation rate at times t/τ = 8, t/τ = 16, t/τ = 32, and t/τ = 64. Plotted in Figure 5, it is obviously observed that at times t/τ = 8 and t/τ = 16, the high distribution value of thermal entropy generation rate mainly dominates in the high gradient fields of temperature, especially near the top and bottom boundaries compared with the corresponding temperature fields in Figure 2. The low distribution value of thermal entropy generation rate mainly occurs in the homogenetic temperature fields. It is seen that with spatial-temporal evolution, the high distribution value of thermal entropy generation rate gradually increases due to the interaction and strike of these thermal plumes at time t/τ = 32, which indicates that the order degree of thermal movement gradually tends to be disordered in the whole closed system. However, the plume scale of thermal entropy generation rate gradually decreases at time t/τ = 64 compared to that of t/τ = 32, which further demonstrates that a great deal of large scale turbulent structures interact and develop into a large of various vortex, which indicates that the viscous flow loss mainly occurs in the high shear region. Meanwhile, at time step t/τ = 32, the viscous entropy generation rate evidently increases with temporal evolution. Nevertheless, the viscous entropy generation rate evidently decreases at time t/τ = 64 compared to that of t/τ = 32, which indicates that the whole mainstream field has already entered the state of turbulent dissipation.  Figure 5 illustrates the thermal entropy generation rate at times t/τ = 8, t/τ = 16, t/τ = 32, and t/τ = 64. Plotted in Figure 5, it is obviously observed that at times t/τ = 8 and t/τ = 16, the high distribution value of thermal entropy generation rate mainly dominates in the high gradient fields of temperature, especially near the top and bottom boundaries compared with the corresponding temperature fields in Figure 2. The low distribution value of thermal entropy generation rate mainly occurs in the homogenetic temperature fields. It is seen that with spatial-temporal evolution, the high distribution value of thermal entropy generation rate gradually increases due to the interaction and strike of these thermal plumes at time t/τ = 32, which indicates that the order degree of thermal movement gradually tends to be disordered in the whole closed system. However, the plume scale of thermal entropy generation rate gradually decreases at time t/τ = 64 compared to that of t/τ = 32, which further demonstrates that a great deal of large scale turbulent structures interact and develop into a large number of small scale turbulent structures; the thermal dissipation also appears with time evolution in the closed system.  Figure 6 describes the total entropy generation rate at times t/τ = 8, t/τ = 16, t/τ = 32, and t/τ = 64. As described in Figure 6, it can be seen that the high distribution value of total entropy generation rate mainly dominates in the largest temperature velocity gradient compared with the corresponding temperature fields in Figure 2. The low total entropy generation rate mainly clusters in the region of the homogenetic temperature fields. The distribution size trend of total entropy generation rate is  Figure 6 describes the total entropy generation rate at times t/τ = 8, t/τ = 16, t/τ = 32, and t/τ = 64. As described in Figure 6, it can be seen that the high distribution value of total entropy generation rate mainly dominates in the largest temperature velocity gradient compared with the corresponding temperature fields in Figure 2. The low total entropy generation rate mainly clusters in the region of the homogenetic temperature fields. The distribution size trend of total entropy generation rate is well consistent with that of thermal entropy generation rate in the corresponding time step. In the spatial evolution, the shape of high entropy generation rate congeals into a large number of varied plumes, which indicates that the role of thermal entropy generation rate gradually improves with time evolution in the heat transfer irreversibility. It can be clearly seen that with time evolution, a great deal of large scale plumes interact and develop to a large number of small scale plumes in the closed system, and the value of total entropy generation rate increases, which indicates that the order degree of energy dissipation in the whole closed system gradually tends to be disordered and increase. The viscous, thermal and total entropy generation rates with evolution can promote the idea that the type of mixed bottom boundary condition and thermal configuration can be extensively applied in a wide variety of practical engineering applications, such as the solar thermal absorber plate or the electronic existing plates.

Quantitative Analysis of Entropy Generation Rate with Time Evolution
The probability density function (PDF) is used to reveal the distribution aggregation situation of physics variable. Wei et al. [24] argued that the PDFs of Su, Sθ and S with increase of Prandtl number, the tails of high entropy generation rates can fit well into the curve of the log-normal coordinate and the departure and the distribution of log-normality, gradually becoming more robust with the decrease of Prandtl number. In this paper, an exponential expression is implemented for PDF. Its exponential expression is as follows [24]:

Quantitative Analysis of Entropy Generation Rate with Time Evolution
The probability density function (PDF) is used to reveal the distribution aggregation situation of physics variable. Wei et al. [24] argued that the PDFs of S u , S θ and S with increase of Prandtl number, the tails of high entropy generation rates can fit well into the curve of the log-normal coordinate and the departure and the distribution of log-normality, gradually becoming more robust with the decrease of Prandtl number. In this paper, an exponential expression is implemented for PDF. Its exponential expression is as follows [24]: in which m, α and C represent the fitted parameters, and Y = X − X mp with X = S u /(S u ) rms , S θ /(S θ ) rms , S/(S) rms and X mp being the abscissa of the most probable amplitude. The best fit of Equation (21)  S θ and S in the whole region. Figure 7 describes the PDFs' distributions of · S u at four times t/τ = 8, t/τ = 16, t/τ = 32 and t/τ = 64. As described in Figure 7, we can see that the high value of · S u decreases in a range of · S u > 10 with time evolution. This is mainly due to the fact that the flow characteristic velocity of the large-scale flow in early characteristic time at t/τ= 8 is relatively large, the large-scale flow is broken into more small-scale flows, the viscosity entropy rate decreases in high value with time evolution, and the viscosity entropy generation rate of the small-scale flow is smaller than that of the large-scale flow.  Figure 9 illustrates the PDF distributions of total entropy generation rate 0 S ⋅ at four times t/τ = 8, t/τ = 16, t/τ = 32, and t/τ = 64.
As illustrated in Figure 9, it can be seen that the high values of 0 S ⋅ keep almost the same with θ S ⋅ , which indicates that the thermal entropy generation rate has a dominant position in the total entropy generation rate with time evolution.     Figure 9 illustrates the PDF distributions of total entropy generation rate · S 0 at four times t/τ = 8, t/τ = 16, t/τ = 32, and t/τ = 64.
As illustrated in Figure 9, it can be seen that the high values of · S 0 keep almost the same with · S θ , which indicates that the thermal entropy generation rate has a dominant position in the total entropy generation rate with time evolution.
To further reveal the distribution differences of S u , S θ and S with time evolution, the average value of  Figure 10 shows the time evolution of average viscous entropy generation rate from t/τ = 0 to t/τ = 100 in the whole field. Plotted in Figure 10, it is clearly observed that the average value of _ · S u alternately increases, three peaks successively appear from the time step of t/τ = 0 to t/τ = 32 with time evolution. One strong peak appears at the time step of t/τ = 32, however, the average value of _ · S u gradually decreases from the time step of t/τ = 32 to t/ τ = 64, and the average value of _ · S u gradually increases in a range of t/τ > 64. This is mainly due to the fact that the largest length-scales eddy is produced owing to the injection of energy at an early characteristic time; the decrease of flow eddies and the geometric eddy size is associated with the characteristic time scales (t/τ < 32). However, with time evolution, the large-scale flow is broken into more small-scale flows in a range of t/τ from 32 to 64; the viscosity entropy rate decreases in high value with time evolution, and the viscosity entropy generation rate of the small-scale flow is smaller than that of the large-scale flow. In a range of t/τ > 64, some of the smallest eddies can be distorted in this distortion process, which further indicates that the kinetic energy may be dissipated from the dissipation of the smallest eddies owing to the effect of viscous flow.     Figure 10 shows the time evolution of average viscous entropy generation rate from t/τ = 0 to t/τ = 100 in the whole field. Plotted in Figure 10, it is clearly observed that the average value of   Figure 11 illustrates the time evolution of average thermal entropy generation rate from t/τ = 0 to t/τ = 100 in the whole field. As illustrated in Figure 11, one can clearly see that at first _ · S θ is very large due to the extremely thin boundary layer. As time goes by, the boundary layer thickness rapidly increases to the normal level, and _ · S θ decreases rapidly. After the initial period, the average value of the temperature generation rate _ · S θ alternately increases, several peaks periodically appear from the time step of t/τ = 0 to t/τ = 32 with time evolution. One strong peak appears at the time step of t/τ = 32, however, the average value of _ · S θ periodically decreases from the time step of t/τ = 32 to t/τ = 64.
The average value of _ · S θ periodically and lightly increases in a range of t/τ > 64. This is mainly due to the fact that the large scale plumes produce owing to the injection of energy at the early characteristic time, the decrease of thermal plumes size is associated with the characteristic time-scales (t/τ < 32). Nevertheless, with time evolution, the large-scale plumes are broken into more small-scale plumes in a range of t/τ from 32 to 64, the thermal entropy rate decreases in high value with time evolution, and the thermal entropy generation rate of the small-scale plumes is smaller than that of the large-scale plumes. In a range of t/τ > 64, some of the smallest plumes can be coagulated into the big plumes.
from the time step of t/τ = 0 to t/τ = 32 with time evolution. One strong peak appears at the time step of t/τ = 32,however, the average value of _ u S ⋅ gradually decreases from the time step of t/τ = 32 to t/τ = 64, and the average value of _ u S ⋅ gradually increases in a range of t/τ > 64. This is mainly due to the fact that the largest length-scales eddy is produced owing to the injection of energy at an early characteristic time; the decrease of flow eddies and the geometric eddy size is associated with the characteristic time scales (t/τ < 32). However, with time evolution, the large-scale flow is broken into more small-scale flows in a range of t/τ from 32 to 64; the viscosity entropy rate decreases in high value with time evolution, and the viscosity entropy generation rate of the small-scale flow is smaller than that of the large-scale flow. In a range of t/τ > 64, some of the smallest eddies can be distorted in this distortion process, which further indicates that the kinetic energy may be dissipated from the dissipation of the smallest eddies owing to the effect of viscous flow.  Figure 11 illustrates the time evolution of average thermal entropy generation rate from t/τ = 0 to t/τ = 100 in the whole field. As illustrated in Figure 11, one can clearly see that at first The average value of _ θ S ⋅ periodically and lightly increases in a range of t/τ > 64. This is mainly due to the fact that the large scale plumes produce owing to the injection of energy at the early characteristic time, the decrease of thermal plumes size is associated with the characteristic timescales (t/τ < 32). Nevertheless, with time evolution, the large-scale plumes are broken into more smallscale plumes in a range of t/τ from 32 to 64, the thermal entropy rate decreases in high value with time evolution, and the thermal entropy generation rate of the small-scale plumes is smaller than that of the large-scale plumes. In a range of t/τ > 64, some of the smallest plumes can be coagulated into the big plumes. Figure 11. Time evolution of average thermal entropy generation rate in the whole two-dimensional field. Figure 12 shows the time evolution of average total entropy generation rate from t/τ = 0 to t/τ = 100 in the whole field. As shown in Figure 12, one can clearly see that at first indicates that the thermal entropy generation rate plays a dominated role in the heat transfer irreversibility-the viscous entropy generation can be neglected time evolution. The above phenomenon is well consistent with the importance of heat transfer irreversibility in the previous Figure 11. Time evolution of average thermal entropy generation rate in the whole two-dimensional field. Figure 12 shows the time evolution of average total entropy generation rate from t/τ = 0 to t/τ = 100 in the whole field. As shown in Figure 12, one can clearly see that at first _ · S 0 is very large due to the extremely thin boundary layer. With time evolution, it is clearly seen that the boundary layer thickness rapidly increases to the normal level, and _ · S 0 decreases rapidly. Plotted in Figure 12, it can be seen that the high values of · S 0 remain almost the same with · S θ with time evolution, which indicates that the thermal entropy generation rate plays a dominated role in the heat transfer irreversibility-the viscous entropy generation can be neglected time evolution. The above phenomenon is well consistent with the importance of heat transfer irreversibility in the previous studies [23][24][25]. Wei et al. [24] studied the effect of changing the Prandtl number on the entropy generation rate in two-dimensional RB convection, and argued that the thermal entropy generation rate has a dominant role in the heat transfer irreversibility-the viscous entropy generation can be neglected with the increasing Prandtl

Conclusions
In this paper, the time evolution features of entropy generation rate in turbulent Rayleigh-Bénard convection are investigated in mixed insulating and conducting boundary conditions. Several conclusions are given as follows.
The physical system of thermal convection gradually evolves from the large-scale to small-scale thermal plumes-some of the smallest plumes can be coagulated into the big plumes and several big plumes reappear in the time evolution. The state of flow gradually develops from laminar flow to turbulent thermal convection motion, and further evolves from turbulent thermal convection to dissipation flow in the process of turbulent energy transfer.
The viscous, thermal, and total entropy generation rates evidently increase in wide range of t/τ < 32 with temporal evolution. Nevertheless, the viscous, thermal, and total entropy generation rates evidently decreases at time t/τ = 64 compared to that of t/τ = 32.
The high value of u S ⋅ decreases in a range of 10 u S ⋅ > with time evolution. It is revealed that the flow characteristic velocity of the large-scale flow in early characteristic time at t/τ = 8 is relatively large: the large-scale flow is broken into more small-scale flows, the viscosity entropy rate decreases in high value with time evolution, and the viscosity entropy generation rate of the small-scale flow is smaller than that of the large-scale flow.
It was seen that the largest length-scale eddy produces owing to the injection of energy at the early characteristic time, the decrease of flow eddies and geometric eddy size are associated with the characteristic time-scales (t/τ < 32). However, the large-scale flow is broken into more small-scale flows in a range of t/τ from 32 to 64, the viscosity entropy rate decreases in high value with time evolution, and the viscosity entropy generation rate of the small-scale flow is smaller than that of the large-scale flow. In a range of t/τ > 64, some of the smallest eddies can be distorted in this distortion process. The average value of _ θ S ⋅ alternately increases from the time step of t/τ = 0 to t/τ = 32,however,

Conclusions
In this paper, the time evolution features of entropy generation rate in turbulent Rayleigh-Bénard convection are investigated in mixed insulating and conducting boundary conditions. Several conclusions are given as follows.
The physical system of thermal convection gradually evolves from the large-scale to small-scale thermal plumes-some of the smallest plumes can be coagulated into the big plumes and several big plumes reappear in the time evolution. The state of flow gradually develops from laminar flow to turbulent thermal convection motion, and further evolves from turbulent thermal convection to dissipation flow in the process of turbulent energy transfer.
The viscous, thermal, and total entropy generation rates evidently increase in wide range of t/τ < 32 with temporal evolution. Nevertheless, the viscous, thermal, and total entropy generation rates evidently decreases at time t/τ = 64 compared to that of t/τ = 32.
The high value of · S u decreases in a range of · S u > 10 with time evolution. It is revealed that the flow characteristic velocity of the large-scale flow in early characteristic time at t/τ = 8 is relatively large: the large-scale flow is broken into more small-scale flows, the viscosity entropy rate decreases in high value with time evolution, and the viscosity entropy generation rate of the small-scale flow is smaller than that of the large-scale flow.
It was seen that the largest length-scale eddy produces owing to the injection of energy at the early characteristic time, the decrease of flow eddies and geometric eddy size are associated with the characteristic time-scales (t/τ < 32). However, the large-scale flow is broken into more small-scale flows in a range of t/τ from 32 to 64, the viscosity entropy rate decreases in high value with time evolution, and the viscosity entropy generation rate of the small-scale flow is smaller than that of the large-scale flow. In a range of t/τ > 64, some of the smallest eddies can be distorted in this distortion process. The average value of _ · S θ alternately increases from the time step of t/τ = 0 to t/τ = 32, however, the average value of _ · S θ periodically decreases from the time step of t/τ = 32 to t/τ = 64. Interestingly, the average value of _ · S θ periodically and lightly increases in a range of t/τ > 64. The above studies further demonstrate the process of the energy injection, the kinetic energy, the kinetic energy transfer, the thermal energy transfer, the viscous dissipated flow and thermal dissipation. In practical engineering, the type of mixed-bottom boundary condition and thermal configuration can be extensively applied in a wide variety of equipment, such as the solar thermal absorber plate or the electronic existing plates.