Effects of Movable-Baffle on Heat Transfer and Entropy Generation in a Cavity Saturated by CNT Suspensions : Three-Dimensional Modeling

Convective heat transfer and entropy generation in a 3D closed cavity, equipped with adiabatic-driven baffle and filled with CNT (carbon nanotube)-water nanofluid, are numerically investigated for a range of Rayleigh numbers from 103 to 105. This research is conducted for three configurations; fixed baffle (V = 0), rotating baffle clockwise (V+) and rotating baffle counterclockwise (V−) and a range of CNT concentrations from 0 to 15%. Governing equations are formulated using potential vector vorticity formulation in its three-dimensional form, then solved by the finite volume method. The effects of motion direction of the inserted driven baffle and CNT concentration on heat transfer and entropy generation are studied. It was observed that for low Rayleigh numbers, the motion of the driven baffle enhances heat transfer regardless of its direction and the CNT concentration effect is negligible. However, with an increasing Rayleigh number, adding driven baffle increases the heat transfer only when it moves in the direction of the decreasing temperature gradient; elsewhere, convective heat transfer cannot be enhanced due to flow blockage at the corners of the baffle.


Introduction
Mixed convection heat transfer can be found in many industrial and engineering applications such as electronic component cooling, the food drying process, nuclear reactors, chemical processing equipment and so on.Several heat transfer enhancement techniques have been used to make efficient thermal devices by improving their thermal performance.For convective heat transfer, the increase of the effective surface area is one of the ways to ameliorate the heat transfer rate but it leads to an increase of the manufacturing cost.Removing this problem can be achievable by using a water-based nanofluid and inserts to generate the swirl in the bulk of the fluids and disturb the actual boundary layer, resulting in an increase of the effective surface area, residence time and consequently the convective heat transfer coefficient in the existing system.Nanofluids are stable suspension of the nanoparticles in a base fluid such as water, ethylene glycol and industrial oils.Nanoparticles can be metallic such as Copper and Silver, metal oxide such as CuO and Al 2 O 3 and nonmetallic such as carbon nanotubes (CNT) and SiC [1].Nanofluids have been applied in many energy systems, especially renewable energy-based systems [2][3][4].
A sizable amount of research interests has been addressed in the last decades on natural and mixed convection in 2D nanofluid-filled cavities.The effect of different nanofluid on the entropy generation in 2D porous cavity was studied numerically by Ashorynejad and Hoseinpour [5].In addition to the type of nanofluids, authors investigated the porosity and volume fraction of nanoparticles.The main result was that entropy generation is increasing with porosity and is decreasing with the nanoparticle volume fraction for all type of nanofluids.Cianfrini et al. [6] carried out a numerical study on the effect of the length and the position of the heater on nanofluid laminar natural convection inside a partially heated square cavity.Based on the obtained results and in addition to the flow structure and temperature field, authors presented correlations for the Nusselt number.A numerical study was conducted by Wang et al. [7] on 2D natural convection of a square enclosure with cooled right vertical wall.The left vertical wall was sinusoidally oscillated at the constant average temperature of the cavity with adiabatic horizontal walls.The cavity was filled by nanofluid, and the effect of the nanofluid volume fraction on heat transfer was the paramount objective of this research.The conclusions drawn show that the heat transfer and the oscillating behaviors were influenced by nanoparticles.Sheikhzadeh et al. [8] analyzed numerically a 2D steady laminar natural convection heat transfer and flow patterns of Cu-water nanofluid in a cavity with partially active vertical walls, while the horizontal walls were assumed to be thermally insulated.Authors concluded that heat transfer increases with increasing both the Rayleigh number and nanoparticle volume fraction.A magnetohydrodynamics (MHD) mixed convection of nanofluid in a 2D lid-driven square cavity that contains a rotating cylinder was simulated numerically by Selimefendigil and Öztop [9].The main objective was the investigation of the effects of the dimensionless parameters such as the Hartmann number, Richardson number, rotational speed of the cylinder, and the concentration of the nanoparticles.It was found that an increasing Richardson number leads to an increment of heat transfer.However, by increasing the value of the Hartmann number, heat transfer is reduced.Rotation of the cylinder has a considerable effect on the heat transfer enhancement.A 2D mixed convection nanofluid flow in an enclosure with moving top wall and wavy bottom wall was numerically investigated by Abu-Nada and Chamkha [10].It was proved that at the considered Ri, the nanoparticles have a significant effect on the heat transfer enhancement as well as the wavy bottom wall geometry ratios.
As mixed convection in industrial/chemical applications is often turbulent, recent development in the field of entropy generation and irreversibility in buoyancy-driven convection was achieved.Zonta et al. [11] performed DNS (Direct Numerical Simulations) to evaluate NOB (Non-Oberbeck-Boussinesq) effects on momentum and heat transfer mechanisms for stably-stratified turbulence in water flows.They have demonstrated that, for all situations in which NOB effects become significant, the Reynolds analogy is not justified and accurate physical modeling is required to improve parameterization of mixing.In another work, Zonta et al. [12] have analyzed the behavior of entropy fluctuations in turbulent thermal convection; they have performed DNS of a turbulent Rayleigh-Bénard flow inside a vertically confined fluid layer and they have followed the dynamics of pointwise Lagrangian tracers to measure the local quantities of the flow.They have shown that entropy production can be evaluated by looking at the work done by buoyancy on fluid particles.
The above presented literature survey shows that the majority of numerical investigations are restricted to the two-dimensional geometry and only very limited work has been done on 3D lid-driven cavity problems [13][14][15].Furthermore, no attention has been paid to mixed convection in three-dimensional closed cavities with inserted driven baffle.This paper aims to undertake a Entropy 2017, 19, 200 3 of 17 comprehensive analysis of the flow patterns and temperature distributions of mixed convection heat transfer in a 3D enclosure filled with CNT-water nanofluid and equipped with adiabatic-driven baffle.

Physical Problem and Governing Equations
The main cubic geometry parameters are sketched in Figure 1.As shown, the left side of the closed space has a higher temperature than the right side while other walls are assumed to be insulated.The adiabatic baffle inserted on the top wall is allowed to rotate clockwise or counterclockwise with uniform velocity.The enclosure is filled with CNT/water nanofluid and the mixed convection flow is a result of movement of the adiabatic-driven baffle and the temperature gradient.

Physical Problem and Governing Equations
The main cubic geometry parameters are sketched in Figure 1.As shown, the left side of the closed space has a higher temperature than the right side while other walls are assumed to be insulated.The adiabatic baffle inserted on the top wall is allowed to rotate clockwise or counterclockwise with uniform velocity.The enclosure is filled with CNT/water nanofluid and the mixed convection flow is a result of movement of the adiabatic-driven baffle and the temperature gradient.The nanofluid is assumed to be incompressible and the flow is considered laminar.The base fluid and nanoparticles are assumed to be at thermal equilibrium.Table 1 presents the thermophysical properties of the base fluid (water) and the nanoparticles.

Governing Equations
As a numerical method, we had recourse to the potential vector-vorticity formalism which allows, in a 3D configuration, the elimination of the pressure, which is a delicate term to treat [16,17].To eliminate this term, one applies the rotational to the equation of momentum.The potential vector and the vorticity are, respectively, defined by the two following relations: The governing equations are written in non-dimensional forms as: The nanofluid is assumed to be incompressible and the flow is considered laminar.The base fluid and nanoparticles are assumed to be at thermal equilibrium.Table 1 presents the thermo-physical properties of the base fluid (water) and the nanoparticles.

Property
Water CNT

Governing Equations
As a numerical method, we had recourse to the potential vector-vorticity formalism ( → ψ − → ω) which allows, in a 3D configuration, the elimination of the pressure, which is a delicate term to treat [16,17].To eliminate this term, one applies the rotational to the equation of momentum.The potential vector and the vorticity are, respectively, defined by the two following relations: The governing equations are written in non-dimensional forms as: ∂T ∂t The time (t ), velocity ( → V ), vector potential ( → ψ ) and vorticity ( → ω ), are put respectively in their dimensional forms by l 2 /α, α/l, α and l 2 /α.The non-dimensional temperature is defined by: Prandtl and Rayleigh numbers are defined respectively as: In addition, the Richardson number is defined as: The effective density of the nanofluid is given by Kahveci [18] as: The heat capacitance of the nanofluid is expressed by Kahveci [16] as: The effective thermal conductivity of the nanofluid is approximated by Xue [19] as: The effective dynamic viscosity of the nanofluid is given by the Brinkman model [20] as: 2.5  (11)

Boundary Conditions
The boundary conditions for the present problem are given as follows: • Temperature: T = 0 at x = 0, and T = 1 x = 1; ∂T ∂n = 0 at all other walls (i.e., adiabatic).

•
Velocity : • Re • Pr on the vertical sides of the moving fin; • Re • Pr on the bottom horizontal side of the moving fin; The generated entropy (S gen ) is written by Kolsi et al. [17] as: The dimensionless local generated entropy (N s ) is written as: N s_th and N s_ f r are respectively the entropy generations due to heat transfer and viscous effect with: and where ϕ S = α l∆T 2 T 0 is the irreversibility coefficient.
The total dimensionless generated entropy (S tot ) is written as: The Bejan number (Be) is defined as: The local Nusselt number at the hot wall (Nu) is defined as follows: Entropy 2017, 19, 200 6 of 17 while, the average Nusselt number (Nu av ), on the isothermal walls is expressed by: The mathematical model described above is written using FORTRAN PowerStation 4.0 The control volume method is used to discretize Equations (Equations ( 2)-( 4) and ( 12)).The central-difference scheme is used for treating convective terms while the fully implicit procedure is used to discretize the temporal derivatives.The resulting nonlinear algebraic equations are solved using the successive relaxation-iterating scheme.The grids are considered uniform in all directions with clustering nodes on boundaries.

Validation of the Code and Grid Dependency
The numerical code is validated with the 2D work of Sun et al. [21] (Figure 2), that considered the case of a differentially heated cavity with a driven top-wall and conductive triangular fin in a different location.The tests are done by comparing the flow field and temperature distribution and it is found that the written homemade code shows good agreement with literature.For nanofluids' buoyancy-induced flow, the code was validated by comparing with the study of Jahanshahi et al. [22] (Figure 3).Jahanshahi et al. [22] studied the effect of adding SiO 2 nanoparticles in water as base fluid on heat transfer and flow structure inside a differentially heated square cavity.It can be concluded from this comparison that there is good agreement.
The sensitivity analysis has been performed with various grids.The tests have been performed for the spatial meshes of 61 3 , 71 3 , 81 3 and 91 3 .The average Nusselt number on the hot wall is considered as a sensitive parameter.The results of the analysis with Ra = 10 5 , V− and ϕ = 5% are presented in Table 2.The increase in percentage of Nu avg for the grid 91 3 is only 0.149%.Hence considering the computational economy and accuracy, the time step (10 −4 ) and spatial mesh size of 81 3 is selected for all the simulations.The solution is considered acceptable when the following convergence criterion is satisfied for each step of time as 1,2,3 Entropy 2017, 19, 200 6 of 17 while, the average Nusselt number ( av Nu ), on the isothermal walls is expressed by: The mathematical model described above is written using FORTRAN PowerStation 4.0 The control volume method is used to discretize Equations (Equations ( 2)-( 4) and ( 12)).The centraldifference scheme is used for treating convective terms while the fully implicit procedure is used to discretize the temporal derivatives.The resulting nonlinear algebraic equations are solved using the successive relaxation-iterating scheme.The grids are considered uniform in all directions with clustering nodes on boundaries.

Validation of the Code and Grid Dependency
The numerical code is validated with the 2D work of Sun et al. [21] (Figure 2), that considered the case of a differentially heated cavity with a driven top-wall and conductive triangular fin in a different location.The tests are done by comparing the flow field and temperature distribution and it is found that the written homemade code shows good agreement with literature.For nanofluids' buoyancy-induced flow, the code was validated by comparing with the study of Jahanshahi et al. [22] (Figure 3).Jahanshahi et al. [22] studied the effect of adding SiO2 nanoparticles in water as base fluid on heat transfer and flow structure inside a differentially heated square cavity.It can be concluded from this comparison that there is good agreement.
The sensitivity analysis has been performed with various grids.The tests have been performed for the spatial meshes of 61 3 , 71 3 , 81 3 and 91 3 .The average Nusselt number on the hot wall is considered as a sensitive parameter.The results of the analysis with Ra = 10 5 , V− and φ = 5% are presented in Table 2.The increase in percentage of Nuavg for the grid 91 3 is only 0.149%.Hence considering the computational economy and accuracy, the time step (10 −4 ) and spatial mesh size of 81 3 is selected for all the simulations.The solution is considered acceptable when the following convergence criterion is satisfied for each step of time as     It should be noted that the results are present for high concentrations, as much as 15%, to make the variations in the heat transfer rate more visible.It should be assumed that the nanofluids remain Newtonian, so the governing equations are the same.
Figure 4 illustrates the flow behavior for Ra = 10 5 , for (V−), (V = 0) and (V+) cases and different values of carbon nanotube (CNT) concentrations (φ = 0 and φ = 0.15).It is clear that the flow is fully disrupted in the case (V−); in fact, the flow is bypassed by the baffle that moves in the opposite direction of the spontaneous movement of the fluid inside the differentially heated cavity.It becomes obvious, therefore, that the flow becomes multicellular and characterized by the existence of secondary eddies near the top corners of the driven baffle which are known as downstream and upstream secondary eddies.Especially for higher CNT concentration, lower secondary eddies are also revealed near the bottom wall.However, when the baffle is fixed, the flow is stratified horizontally near the bottom wall and vertically near the active walls and characterized by a major clockwise vortex occupying most of the enclosure.The vortex body comprises two cores which are located in the lower part of the cavity compared to the case (V+) where they are located on either side of the driven baffle.These results are more tangible from Figure 5 representing the velocity vector projection in the XY−plane.Table 2. Grid sensitivity analysis; V−, Ra = 10 5 and ϕ = 0.05.

Grid size
Nu avg Percentage Increase Incremental Increase 61 3  2.512 --71 3  2.645 5.302 -81 3  2.686 6.921 1.537 91 3  2.69 7.081 0.149 It should be noted that the results are present for high concentrations, as much as 15%, to make the variations in the heat transfer rate more visible.It should be assumed that the nanofluids remain Newtonian, so the governing equations are the same.
Figure 4 illustrates the flow behavior for Ra = 10 5 , for (V−), (V = 0) and (V+) cases and different values of carbon nanotube (CNT) concentrations (ϕ = 0 and ϕ = 0.15).It is clear that the flow is fully disrupted in the case (V−); in fact, the flow is bypassed by the baffle that moves in the opposite direction of the spontaneous movement of the fluid inside the differentially heated cavity.It becomes obvious, therefore, that the flow becomes multicellular and characterized by the existence of secondary eddies near the top corners of the driven baffle which are known as downstream and upstream secondary eddies.Especially for higher CNT concentration, lower secondary eddies are also revealed near the bottom wall.However, when the baffle is fixed, the flow is stratified horizontally near the bottom wall and vertically near the active walls and characterized by a major clockwise vortex occupying most of the enclosure.The vortex body comprises two cores which are located in the lower part of the cavity compared to the case (V+) where they are located on either side of the driven baffle.These results are more tangible from Figure 5 representing the velocity vector projection in the XY-plane.Figure 6 shows the iso-surfaces of temperature for Ra = 10 5 , and different values of carbon nanotube (CNT) concentrations (φ = 0 and φ = 0.15).As it can be seen from this figure, the structure of the iso-surfaces of temperature is extremely affected by the motion of the driven baffle and its direction.However, when the baffle is fixed, the flow is vertically stratified at the lower part of the hot wall and the upper part of the cold wall while it is laminated obliquely at the central portion of the cavity.It is noticed that the distortions at the core of the cavity are more pronounced for the "V+" case due to the cooperation between the buoyancy forces and the displacement of the baffle (the reverse is true for the case "V−").The effect of adding nanoparticles is also found to increase this distortion for all cases due to heat transfer enhancement.Figure 6 shows the iso-surfaces of temperature for Ra = 10 5 , and different values of carbon nanotube (CNT) concentrations (ϕ = 0 and ϕ = 0.15).As it can be seen from this figure, the structure of the iso-surfaces of temperature is extremely affected by the motion of the driven baffle and its direction.However, when the baffle is fixed, the flow is vertically stratified at the lower part of the hot wall and the upper part of the cold wall while it is laminated obliquely at the central portion of the cavity.It is noticed that the distortions at the core of the cavity are more pronounced for the "V+" case due to the cooperation between the buoyancy forces and the displacement of the baffle (the reverse is true for the case "V−").The effect of adding nanoparticles is also found to increase this distortion for all cases due to heat transfer enhancement.V− V = 0 V+   V− V = 0 V+   Figure 7 depicts isotherms at the z = 0.5 plane to see the effects of the nanofluid volume fraction at different baffle velocities.However, the fluid flows when the lid moves in a different direction and the heat transfer mode changes to convection even at the same value of Rayleigh number which is taken as Ra = 10 5 for all cases.It is noticed that the effects of the nanofluid volume fraction for moving baffle are more pronounced than that for a fixed baffle.The structure of the isotherms is more disturbed for the case "V−" because of the opposition of the convection motion and that of the baffle and the case of "V+" is found to accelerate the flow in the direction of decreasing the temperature gradients, especially for the highest nanofluid volume fraction.
baffle are more pronounced than that for a fixed baffle.The structure of the isotherms is more disturbed for the case "V−" because of the opposition of the convection motion and that of the baffle and the case of "V+" is found to accelerate the flow in the direction of decreasing the temperature gradients, especially for the highest nanofluid volume fraction.
V− V = 0 V+ The average Nusselt number at the hot wall versus the nanoparticle volume fraction is given in Figure 8 for different Rayleigh numbers and situations of the moving baffle.For all cases, the average Nusselt number is a linear function of the nanoparticle volume fraction.It should be noted that, for a lower Rayleigh number (Ra = 10 3 ), the average Nusselt number is almost constant with variation of the nanofluid volume fraction.For Ra = 10 4 , the average Nusselt number is increased almost linearly with the increasing nanoparticle volume fraction for both cases of "V+" and V = 0 due to the increasing thermal conductivity of the nanofluids.However, it decreases for the case of "V−" due to the resistance to the spontaneous flow from the hot wall to the cold one caused by the rotational motion of the baffle in the opposite direction.Moreover, it is observed that, by increasing the Rayleigh number, the mean Nusselt number continues to decrease for the case of "V−".Obviously, the highest heat transfer occurred for the case of "V+" and for the highest value of Rayleigh number.It is indeed an interesting result that the rotation direction of the driven baffle can be a control parameter of the heat transfer rate.The average Nusselt number at the hot wall versus the nanoparticle volume fraction is given in Figure 8 for different Rayleigh numbers and situations of the moving baffle.For all cases, the average Nusselt number is a linear function of the nanoparticle volume fraction.It should be noted that, for a lower Rayleigh number (Ra = 10 3 ), the average Nusselt number is almost constant with variation of the nanofluid volume fraction.For Ra = 10 4 , the average Nusselt number is increased almost linearly with the increasing nanoparticle volume fraction for both cases of "V+" and V = 0 due to the increasing thermal conductivity of the nanofluids.However, it decreases for the case of "V−" due to the resistance to the spontaneous flow from the hot wall to the cold one caused by the rotational motion of the baffle in the opposite direction.Moreover, it is observed that, by increasing the Rayleigh number, the mean Nusselt number continues to decrease for the case of "V−".Obviously, the highest heat transfer occurred for the case of "V+" and for the highest value of Rayleigh number.It is indeed an interesting result that the rotation direction of the driven baffle can be a control parameter of the heat transfer rate.
baffle are more pronounced than that for a fixed baffle.The structure of the isotherms is more disturbed for the case "V−" because of the opposition of the convection motion and that of the baffle and the case of "V+" is found to accelerate the flow in the direction of decreasing the temperature gradients, especially for the highest nanofluid volume fraction.
V− V = 0 V+ The average Nusselt number at the hot wall versus the nanoparticle volume fraction is given in Figure 8 for different Rayleigh numbers and situations of the moving baffle.For all cases, the average Nusselt number is a linear function of the nanoparticle volume fraction.It should be noted that, for a lower Rayleigh number (Ra = 10 3 ), the average Nusselt number is almost constant with variation of the nanofluid volume fraction.For Ra = 10 4 , the average Nusselt number is increased almost linearly with the increasing nanoparticle volume fraction for both cases of "V+" and V = 0 due to the increasing thermal conductivity of the nanofluids.However, it decreases for the case of "V−" due to the resistance to the spontaneous flow from the hot wall to the cold one caused by the rotational motion of the baffle in the opposite direction.Moreover, it is observed that, by increasing the Rayleigh number, the mean Nusselt number continues to decrease for the case of "V−".Obviously, the highest heat transfer occurred for the case of "V+" and for the highest value of Rayleigh number.It is indeed an interesting result that the rotation direction of the driven baffle can be a control parameter of the heat transfer rate.Entropy generation is calculated via the obtained values of temperature and velocities.Thus, local entropy generation due to heat transfer is presented in Figure 9 for different cases and nanofluid volume fractions for Ra = 10 5 .As shown from the figure, the edge of the baffle becomes less effective on entropy generation for V = 0.The accumulation of the entropy contours is at the lower part of the hot wall and the upper part of the cold wall.Moreover, the increase of the nanofluid volume fraction is accompanied by an increase of the developed thermal entropy.However, when the baffle is driven (case V+), the accumulation of the entropy contours spreads over all the active walls and the effect of the nanofluid volume fraction is more perceptible for the case of "V−" for which contours invade the major part of the cavity.It means that flow entrainment becomes a very effective parameter on entropy generation.Entropy generation is calculated via the obtained values of temperature and velocities.Thus, local entropy generation due to heat transfer is presented in Figure 9 for different cases and nanofluid volume fractions for Ra = 10 5 .As shown from the figure, the edge of the baffle becomes less effective on entropy generation for V = 0.The accumulation of the entropy contours is at the lower part of the hot wall and the upper part of the cold wall.Moreover, the increase of the nanofluid volume fraction is accompanied by an increase of the developed thermal entropy.However, when the baffle is driven (case V+), the accumulation of the entropy contours spreads over all the active walls and the effect of the nanofluid volume fraction is more perceptible for the case of "V−" for which contours invade the major part of the cavity.It means that flow entrainment becomes a very effective parameter on entropy generation.
V− V = 0 V+ Similarly, Figure 10 illustrates the local entropy generation due to friction at the z = 0.5 plane for different cases and Ra = 10 5 .For driven baffle and regardless of the motion direction, the edge of the baffle is effective on local entropy generation and contours are located only around the baffle.However, for V = 0, active walls and the bottom participate in entropy production.
V− V = 0 V+ Figure 11 displays the local total entropy generation for the same cases as in Figures 9 and 10 to make a comparison.Figure 12 shows the entropy generation due to heat transfer versus the nanoparticle volume fraction at different Rayleigh numbers and different cases of the baffle.For all cases, an increase of the nanofluid volume fraction is accompanied by an increase in the thermal entropy generation.The same observation can be drawn when increasing the Rayleigh number value except the case of "V−"for which the contrary is found.The highest value of entropy generation due to heat transfer is observed for the case "V+" and the highest value of the Rayleigh number.Similarly, Figure 10 illustrates the local entropy generation due to friction at the z = 0.5 plane for different cases and Ra = 10 5 .For driven baffle and regardless of the motion direction, the edge of the baffle is effective on local entropy generation and contours are located only around the baffle.However, for V = 0, active walls and the bottom participate in entropy production.
Entropy 2017, 19, 200 11 of 17 Entropy generation is calculated via the obtained values of temperature and velocities.Thus, local entropy generation due to heat transfer is presented in Figure 9 for different cases and nanofluid volume fractions for Ra = 10 5 .As shown from the figure, the edge of the baffle becomes less effective on entropy generation for V = 0.The accumulation of the entropy contours is at the lower part of the hot wall and the upper part of the cold wall.Moreover, the increase of the nanofluid volume fraction is accompanied by an increase of the developed thermal entropy.However, when the baffle is driven (case V+), the accumulation of the entropy contours spreads over all the active walls and the effect of the nanofluid volume fraction is more perceptible for the case of "V−" for which contours invade the major part of the cavity.It means that flow entrainment becomes a very effective parameter on entropy generation.
V− V = 0 V+ Similarly, Figure 10 illustrates the local entropy generation due to friction at the z = 0.5 plane for different cases and Ra = 10 5 .For driven baffle and regardless of the motion direction, the edge of the baffle is effective on local entropy generation and contours are located only around the baffle.However, for V = 0, active walls and the bottom participate in entropy production.
V− V = 0 V+ Figure 11 displays the local total entropy generation for the same cases as in Figures 9 and 10 to make a comparison.Figure 12 shows the entropy generation due to heat transfer versus the nanoparticle volume fraction at different Rayleigh numbers and different cases of the baffle.For all cases, an increase of the nanofluid volume fraction is accompanied by an increase in the thermal entropy generation.The same observation can be drawn when increasing the Rayleigh number value except the case of "V−"for which the contrary is found.The highest value of entropy generation due to heat transfer is observed for the case "V+" and the highest value of the Rayleigh number.make a comparison.Figure 12 shows the entropy generation due to heat transfer versus the nanoparticle volume fraction at different Rayleigh numbers and different cases of the baffle.For all cases, an increase of the nanofluid volume fraction is accompanied by an increase in the thermal entropy generation.The same observation can be drawn when increasing the Rayleigh number value except the case of "V−"for which the contrary is found.The highest value of entropy generation due to heat transfer is observed for the case "V+" and the highest value of the Rayleigh number.In the same way, Figure 13 illustrates the entropy generation due to fluid friction versus the nanoparticle volume fraction at different Rayleigh numbers and different cases of the baffle.This figure shows that the entropy generation is dominated by generation due to fluid friction which prevails over that which is due to the transfer of heat.When the baffle is fixed (V = 0), a small linear increase is obtained for a higher Rayleigh number while it is almost constant for the lower value of Rayleigh number.As seen from the figure, the addition of nanoparticles decreases the entropy generation due to friction and the highest value is obtained for the case of "V−" which favors the fluid friction against the driven baffle.Also, variation of overall entropy generation with the nanoparticle volume fraction is shown in Figure 14 for three studied situations of the baffle.The effect of the Rayleigh number is also illustrated in this figure.Regardless of the Rayleigh number, when the baffle is driven, a linear decrease of the overall entropy generation is observed with the increase of the nanofluid volume fraction.Also, the rotation direction has no significant effect on the total generated entropy.However, in all cases, the moving lid generates higher entropy than the fixed one.As mentioned above, the total entropy generation for different nanoparticle volume fractions for different cases and Rayleigh numbers is shown in Figure 14.As expected, the highest values are obtained for case "V−" and the highest value of Ra number.However, the value of Stot decreases with In the same way, Figure 13 illustrates the entropy generation due to fluid friction versus the nanoparticle volume fraction at different Rayleigh numbers and different cases of the baffle.This figure shows that the entropy generation is dominated by generation due to fluid friction which prevails over that which is due to the transfer of heat.When the baffle is fixed (V = 0), a small linear increase is obtained for a higher Rayleigh number while it is almost constant for the lower value of Rayleigh number.As seen from the figure, the addition of nanoparticles decreases the entropy generation due to friction and the highest value is obtained for the case of "V−" which favors the fluid friction against the driven baffle.Also, variation of overall entropy generation with the nanoparticle volume fraction is shown in Figure 14 for three studied situations of the baffle.The effect of the Rayleigh number is also illustrated in this figure.Regardless of the Rayleigh number, when the baffle is driven, a linear decrease of the overall entropy generation is observed with the increase of the nanofluid volume fraction.Also, the rotation direction has no significant effect on the total generated entropy.However, in all cases, the moving lid generates higher entropy than the fixed one.In the same way, Figure 13 illustrates the entropy generation due to fluid friction versus the nanoparticle volume fraction at different Rayleigh numbers and different cases of the baffle.This figure shows that the entropy generation is dominated by generation due to fluid friction which prevails over that which is due to the transfer of heat.When the baffle is fixed (V = 0), a small linear increase is obtained for a higher Rayleigh number while it is almost constant for the lower value of Rayleigh number.As seen from the figure, the addition of nanoparticles decreases the entropy generation due to friction and the highest value is obtained for the case of "V−" which favors the fluid friction against the driven baffle.Also, variation of overall entropy generation with the nanoparticle volume fraction is shown in Figure 14 for three studied situations of the baffle.The effect of the Rayleigh number is also illustrated in this figure.Regardless of the Rayleigh number, when the baffle is driven, a linear decrease of the overall entropy generation is observed with the increase of the nanofluid volume fraction.Also, the rotation direction has no significant effect on the total generated entropy.However, in all cases, the moving lid generates higher entropy than the fixed one.As mentioned above, the total entropy generation for different nanoparticle volume fractions for different cases and Rayleigh numbers is shown in Figure 14.As expected, the highest values are obtained for case "V−" and the highest value of Ra number.However, the value of Stot decreases with As mentioned above, the total entropy generation for different nanoparticle volume fractions for different cases and Rayleigh numbers is shown in Figure 14.As expected, the highest values are obtained for case "V−" and the highest value of Ra number.However, the value of S tot decreases with the increasing value of the nanoparticle volume fraction except the case of fixed baffle (V = 0) for which S tot value increases almost linearly for Ra = 10 5 .
Lastly, Figure 15 gives the Bejan number versus the nanoparticle volume fraction for different Rayleigh numbers and different cases of the baffle.As shown from the figure, the Bejan number increases almost linearly with the increasing nanoparticle volume fraction.The highest value of the Bejan number is obtained for the lowest Rayleigh number value (Ra = 10 3 ) and for fixed baffle (V = 0).The effect of the motion direction of the baffle on the Bejan number is significant only for higher values of Rayleigh number.The variation of Bejan has no significance for lower values of Rayleigh number due to the fact that forced convection dominates the natural convection.As defined above, the Bejan number is the ratio of entropy generation due to heat transfer to total entropy generation.Thus, it is almost 1 for the lowest value of Rayleigh number.It means that entropy generation due to fluid friction becomes insignificant in this case.Lastly, Figure 15 gives the Bejan number versus the nanoparticle volume fraction for different Rayleigh numbers and different cases of the baffle.As shown from the figure, the Bejan number increases almost linearly with the increasing nanoparticle volume fraction.The highest value of the Bejan number is obtained for the lowest Rayleigh number value (Ra = 10 3 ) and for fixed baffle (V = 0).The effect of the motion direction of the baffle on the Bejan number is significant only for higher values of Rayleigh number.The variation of Bejan has no significance for lower values of Rayleigh number due to the fact that forced convection dominates the natural convection.As defined above, the Bejan number is the ratio of entropy generation due to heat transfer to total entropy generation.Thus, it is almost 1 for the lowest value of Rayleigh number.It means that entropy generation due to fluid friction becomes insignificant in this case.The behavior of Nu av is shown in Figure 16 as a function of Ri for different nanoparticle volume fractions and different directions in which the baffle is moved.It can be observed that, under the value of Ri = 0.1, the heat transfer rate is nearly constant regardless of both the nanoparticle volume fraction and the baffle motion direction due to negligible buoyancy in the flow.The effect of motion direction is more pronounced when the Richardson number is in the order of unity and the flow is likely to be buoyancy-driven.The case of "V+" is found to enhance the heat transfer, however, the case of "V−", for which the motion direction of the baffles is opposed to that of the flow, induced entropy generation due to friction leading to a decrease in the heat transfer rate.The maximum in Nu av is reached for higher nanoparticle volume fractions with the case of "V+" of driven baffle.The behavior of Nuav is shown in Figure 16 as a function of Ri for different nanoparticle volume fractions and different directions in which the baffle is moved.It can be observed that, under the value of Ri = 0.1, the heat transfer rate is nearly constant regardless of both the nanoparticle volume fraction and the baffle motion direction due to negligible buoyancy in the flow.The effect of motion direction is more pronounced when the Richardson number is in the order of unity and the flow is likely to be buoyancy-driven.The case of "V+" is found to enhance the heat transfer, however, the case of "V−", for which the motion direction of the baffles is opposed to that of the flow, induced entropy generation due to friction leading to a decrease in the heat transfer rate.The maximum in Nuav is reached for higher nanoparticle volume fractions with the case of "V+" of driven baffle.< , the higher the nanoparticle volume fraction, the lower the total entropy generation.On the other hand, when Ri exceeds unity, there is an increase in the value of Stot in all cases due to an intensification of natural convection when buoyancy is dominant and there is insufficient kinetic energy to homogenize the flow.The maximum of total entropy generation occurs for the case of "V−" and the smaller nanoparticle volume fraction.

Conclusions
A numerical analysis of mixed convection heat transfer in a 3D enclosure filled with CNT-water nanofluid and equipped with adiabatic-driven baffle was performed using computational fluid dynamics (CFD).Some conclusions can be drawn as follows: 1. Driven baffle direction plays an important role in enhancing heat transfer and can be the main parameter to control the desired rate.2. Adding a low concentration of nanoparticles to pure fluid enhances the heat transfer, especially at high Rayleigh number values.Figure 17 displays changes of the total entropy generation in terms of Ri.As it can be seen, for both motion directions there is only a slight variation in the value of S tot which is nearly constant until the Richardson number reaches unity for a fixed nanoparticle volume fraction value.Moreover, for Ri < 1, the higher the nanoparticle volume fraction, the lower the total entropy generation.On the other hand, when Ri exceeds unity, there is an increase in the value of S tot in all cases due to an intensification of natural convection when buoyancy is dominant and there is insufficient kinetic energy to homogenize the flow.The maximum of total entropy generation occurs for the case of "V−" and the smaller nanoparticle volume fraction.The behavior of Nuav is shown in Figure 16 as a function of Ri for different nanoparticle volume fractions and different directions in which the baffle is moved.It can be observed that, under the value of Ri = 0.1, the heat transfer rate is nearly constant regardless of both the nanoparticle volume fraction and the baffle motion direction due to negligible buoyancy in the flow.The effect of motion direction is more pronounced when the Richardson number is in the order of unity and the flow is likely to be buoyancy-driven.The case of "V+" is found to enhance the heat transfer, however, the case of "V−", for which the motion direction of the baffles is opposed to that of the flow, induced entropy generation due to friction leading to a decrease in the heat transfer rate.The maximum in Nuav is reached for higher nanoparticle volume fractions with the case of "V+" of driven baffle.< , the higher the nanoparticle volume fraction, the lower the total entropy generation.On the other hand, when Ri exceeds unity, there is an increase in the value of Stot in all cases due to an intensification of natural convection when buoyancy is dominant and there is insufficient kinetic energy to homogenize the flow.The maximum of total entropy generation occurs for the case of "V−" and the smaller nanoparticle volume fraction.

Conclusions
A numerical analysis of mixed convection heat transfer in a 3D enclosure filled with CNT-water nanofluid and equipped with adiabatic-driven baffle was performed using computational fluid dynamics (CFD).Some conclusions can be drawn as follows: 1. Driven baffle direction plays an important role in enhancing heat transfer and can be the main parameter to control the desired rate.2. Adding a low concentration of nanoparticles to pure fluid enhances the heat transfer, especially at high Rayleigh number values.

Conclusions
A numerical analysis of mixed convection heat transfer in a 3D enclosure filled with CNT-water nanofluid and equipped with adiabatic-driven baffle was performed using computational fluid dynamics (CFD).Some conclusions can be drawn as follows: 1.
Driven baffle direction plays an important role in enhancing heat transfer and can be the main parameter to control the desired rate.

2.
Adding a low concentration of nanoparticles to pure fluid enhances the heat transfer, especially at high Rayleigh number values.

3.
The maximum rate of heat transfer occurred for the case of "V+" and for the highest value of Rayleigh number.

4.
The highest values of total generated entropy are obtained for case "V−" and the highest value of Ra number.5.
Using moving baffles is not always an efficient technique for enhancing the convective heat transfer at low and moderate Rayleigh number values due to blocking flow at the corners of baffles.It can improve the heat transfer only at high Rayleigh number values where natural convection is the dominant mechanism.

Figure 2 .
Figure 2. Validation of the code with 2D solution from the literature for Ri = 10; (a) study of Sun et al. [21] (on the left) and (b) present study (on the right).Figure 2. Validation of the code with 2D solution from the literature for Ri = 10; (a) study of Sun et al. [21] (on the left) and (b) present study (on the right).

Figure 2 .
Figure 2. Validation of the code with 2D solution from the literature for Ri = 10; (a) study of Sun et al. [21] (on the left) and (b) present study (on the right).Figure 2. Validation of the code with 2D solution from the literature for Ri = 10; (a) study of Sun et al. [21] (on the left) and (b) present study (on the right).

Figure 7
Figure 7 depicts isotherms at the z = 0.5 plane to see the effects of the nanofluid volume fraction at different baffle velocities.However, the fluid flows when the lid moves in a different direction and the heat transfer mode changes to convection even at the same value of Rayleigh number which is taken as Ra = 10 5 for all cases.It is noticed that the effects of the nanofluid volume fraction for moving

Figure 7
Figure 7 depicts isotherms at the z = 0.5 plane to see the effects of the nanofluid volume fraction at different baffle velocities.However, the fluid flows when the lid moves in a different direction and the heat transfer mode changes to convection even at the same value of Rayleigh number which is taken as Ra = 10 5 for all cases.It is noticed that the effects of the nanofluid volume fraction for moving

Figure 11
Figure 11 displays the local total entropy generation for the same cases as in Figures 9 and 10to make a comparison.Figure12shows the entropy generation due to heat transfer versus the nanoparticle volume fraction at different Rayleigh numbers and different cases of the baffle.For all cases, an increase of the nanofluid volume fraction is accompanied by an increase in the thermal entropy generation.The same observation can be drawn when increasing the Rayleigh number value except the case of "V−"for which the contrary is found.The highest value of entropy generation due to heat transfer is observed for the case "V+" and the highest value of the Rayleigh number.
Figure 11 displays the local total entropy generation for the same cases as in Figures 9 and 10to make a comparison.Figure12shows the entropy generation due to heat transfer versus the nanoparticle volume fraction at different Rayleigh numbers and different cases of the baffle.For all cases, an increase of the nanofluid volume fraction is accompanied by an increase in the thermal entropy generation.The same observation can be drawn when increasing the Rayleigh number value except the case of "V−"for which the contrary is found.The highest value of entropy generation due to heat transfer is observed for the case "V+" and the highest value of the Rayleigh number.

Entropy 2017 ,
19, 200  13 of 17    the increasing value of the nanoparticle volume fraction except the case of fixed baffle (V = 0) for which Stot value increases almost linearly for Ra = 10 5 .

Figure 16 .
Figure 16.The Nusselt number as a function of Ri.

Figure 17
Figure17displays changes of the total entropy generation in terms of Ri.As it can be seen, for both motion directions there is only a slight variation in the value of Stot which is nearly constant until the Richardson number reaches unity for a fixed nanoparticle volume fraction value.Moreover, for Ri 1 < , the higher the nanoparticle volume fraction, the lower the total entropy generation.On the

Figure 17 .
Figure 17.Total entropy generation as a function of Ri.

15 Figure 16 .
Figure 16.The Nusselt number as a function of Ri.

Figure 16 .
Figure 16.The Nusselt number as a function of Ri.

Figure 17
Figure17displays changes of the total entropy generation in terms of Ri.As it can be seen, for both motion directions there is only a slight variation in the value of Stot which is nearly constant until the Richardson number reaches unity for a fixed nanoparticle volume fraction value.Moreover, for Ri 1 < , the higher the nanoparticle volume fraction, the lower the total entropy generation.On the

Figure 17 .
Figure 17.Total entropy generation as a function of Ri.

15 Figure 17 .
Figure 17.Total entropy generation as a function of Ri.

Table 1 .
Thermo-physical properties of water and carbon nanotube (CNT) nanoparticles at 25 °C.

Table 1 .
Thermo-physical properties of water and carbon nanotube (CNT) nanoparticles at 25 • C.