Entropy Generation and Heat Transfer Performances of Al 2 O 3 -water Nanofluid Transitional Flow in Rectangular Channels with Dimples and Protrusions

Nanofluid has great potentials in heat transfer enhancement and entropy generation decrease as an effective cooling medium. Effects of Al 2 O 3-water nanofluid flow on entropy generation and heat transfer performance in a rectangular conventional channel are numerically investigated in this study. Four different volume fractions are considered and the boundary condition with a constant heat flux is adopted. The flow Reynolds number covers laminar flow, transitional flow and turbulent flow. The influences of the flow regime and nanofluid volume fraction are examined. Furthermore, dimples and protrusions are employed, and the impacts on heat transfer characteristic and entropy generation are acquired. It is found that the average heat transfer entropy generation rate descends and the average friction entropy generation rate rises with an increasing nanofluid volume fraction. The effect of nanofluid on average heat transfer entropy generation rate declines when Reynolds number ascends, which is inverse for average friction entropy generation rate. The average wall temperature and temperature uniformity both drop accompanied with increasing pumping power with the growth in nanofluid volume fraction. The employment of dimples and protrusions significantly decreases the average entropy generation rate and improve the heat transfer performance. The effect of dimple-case shows great difference with that of protrusion-case.


Introduction
As the remarkable development of the electronic equipment technology and consequent continuous growth of the device output power, the requirement of superior cooling technique has been dramatically drawing the researchers' attention.Due to limited cooling capacity of the conventional working substance, it is unable to meet the highly demanding cooling requirement.In the past few years, the nanofluids cooling technology has emerged since the higher thermal conductivities compared with the base fluids, and it has great potential in heat transfer enhancement as a promising cooling medium.
Choi [1] defined the dilute suspensions in which nano-sized particles are dispersed in traditional liquids as "nanofluids" in which the nanoparticles cover metals, metal-oxides, polymers, silica or even carbon nanotubes and the base fluids generally include water, oil, or ethylene glycol.Saidur et al. [2] discussed the various applications of nanofluids and listed the challenges of nanofluids in thermal engineering systems applications.A great number of investigations have been conducted on nanofluids thermophysical properties.Haddad et al. [3] numerically and experimentally investigated the natural convection of nanofluids in various types of cavities.Khanafer and Vafai [4] examined published models for calculating nanofluids thermophysical properties, mainly including thermal conductivity, viscosity, specific heat, and density.Mahbubul et al. [5] accomplished a comparison of different effective models for the viscosity of nanofluids.Fan and Wang [6] provided a review on heat conduction of nanofluids in which they focused on thermal conductivity of nanofluids.Brinkman [7] obtained an expression for predicting viscosities of solutions and suspensions within finite concentrations, and the influence of the additional one solute-molecule was taken into consideration.The Brinkman model was obtained as the most common relation that has been used to calculate the viscosity in the entropy generation problems.Corcione [8] carried out a theoretical investigation on the heat transfer characteristics of buoyancy-driven nanofluids in rectangular enclosures with vertical walls heated discriminately and presented an empirical correlation by regression analysis of experimental data.Duangthongsuk and Wongwises [9] experimentally examined the viscosity and thermal conductivity of nanofluids where TiO 2 nanoparticles are dispersed in water, based on which the expressions for predicting the viscosity of TiO 2 -water nanofluids were presented.The Maxwell model [10] is one of the earliest approaches to calculate the thermal conductivity of solid-liquid mixtures.Bruggeman [11] presented a thermal conductivity model considering interactions among the spherical particles, which can be used for high quantities of particle loadings of nanoparticles with in spherical shape.Hamilton and Crosser [12] carried out an investigation on the effect of included nanoparticle shape, composition, and pure component on the thermal conductivity of various two-component mixtures composed by continuous and discontinuous phases.They acquired a modified Maxwell model by taking the shape of nanoparticles into consideration.Nan et al. [13] proposed a methodology on the evaluation of the effective thermal conductivity of the arbitrary particulate composite.Based on the effective medium approach, the interfacial thermal resistance is considered, and the essential concept of Kapitza thermal contact resistance is also taken into account.A good agreement was found in the comparison of the results predicted with the existing models and experimental results.Taking the effect of the nanolayer between the liquid molecules and solid surface into consideration, the Maxwell model for predicting thermal conductivities has been modified by Yu and Choi [14] and a renovated Maxwell model was obtained.Then, they [15] proposed a renovated Hamilton-Crosser model that can be employed for nonspherical particles.With the renovated model, a correct prediction of the magnitude of the thermal conductivity of nanotube-in-oil nanofluids was acquired.Koo and Kleinstreuer [16] presented a model for calculating the thermal conductivity of nanofluids (K-K model) and they took account of the influences of the particle size, volume fraction and temperature dependence, and effects of the properties of base fluid and nanoparticles are also seriously considered.Feng and Kleinstreuer [17] proposed a more advanced theory that was free of any matching function or coefficient (F-K model), based on an analogy between the random nanoparticle fluctuations generated by the Brownian motion effect and turbulent fluctuations.Considering the extended irreversible thermodynamics, Machrafi and Lebon [18] presented a model for predicting nanofluids thermal conductivities, and the effect of coupled heat transfer mechanisms is taken into great consideration, including the interfacial layering between the base fluid and nanoparticles, Brownian motion and particles agglomeration.Besides, an assessment of the effect on the thermal conductivity of each mechanism was given.
Besides, a large number of studies on convective heat transfer enhancement by nanofluids have also been accomplished.Seyf and Feizbakhshi [19] carried out a numerical investigation on the application of nanofluids in Micro-Pin-Fin Heat Sinks.De-ionized water was selected to be the base fluid and the nanoparticles employed were CuO and Al 2 O 3 nanoparticles with different mean diameters.Selvakumar and Suresh [20] conducted a study on the convective heat transfer of CuO-water nanofluids and a thin-channeled copper water block was adopted.Diao et al. [21] achieved a research on the heat transfer characteristic of a microchannel surface at different pressures, and the effect of Al 2 O 3 -R141b nanofluids with various particle concentrations was considered.Wen and Ding [22] conducted a study on the convective heat transfer of γ-Al 2 O 3 -deionized water nanofluids in a copper tube under laminar flow.Suresh et al. [23] experimentally examined the convective heat transfer and friction characteristics of distilled water and CuO-water nanofluids in flat and dimpled tube under laminar flow, and a constant heat flux boundary condition was adopted.Later on, they conducted a similar study under turbulent flow [24].Vakili et al. [25] completed a research on the heat transfer behavior of TiO 2 nanofluid with different nanoparticle concentrations in a vertical pipe at different constant heat fluxes.Xuan and Li [26] built an experimental system to investigate the convective heat transfer and flow characteristics of nanofluids in a tube, acquiring a new convective heat transfer correlation for correlating experimental data of nanofluids.Gavili et al. [27] conducted a study on the heat transfer and flow characteristics of Al 2 O 3 -water nanofluid in a two-sided lid-driven rectangular cavity with walls heated discriminately.Mohammed et al. [28] achieved an investigation on the heat transfer and flow characteristics of nanofluids in a channel with ribs and grooves under turbulent flow.Nanoparticles, Al 2 O 3 , CuO, SiO 2 , and ZnO, and base fluids, water, glycerin, and engine oil, are all examined, respectively.Besides, the effects of the volume fraction and nanoparticle diameter are further studied.
Furthermore, except for enhancing the convection heat transfer as much as possible, it is also of significant importance to reduce the entropy generation to the greatest extent.Entropy generation determines the level of irreversibility accumulation during the flow and heat transfer process.Consequently, entropy generation is usually used for evaluating the performance of engineering devices.Oztop and Al-Salem [29] reported a review on the entropy generation in natural and mixed convection heat transfer in energy systems.Singh et al. [30] theoretically investigated the effect on entropy generation of Al 2 O 3 -water nanofluid in channels in three different sizes containing microchannels, minichannels, and conventional channels.The results under laminar flow and turbulent flow are both considered, and a prediction of the entropy generation rate was provided through an order of magnitude approach.Li and Kleinstreuer [31] numerically studied the entropy generation of CuO-water nanofluid in trapezoidal microchannels.The volume fractions lower than 4% and steady laminar developing flow are employed.Moghaddami et al. [32] presented an estimation of the entropy generation of Al 2 O 3 -water and Al 2 O 3 -EG nanofluids in a circular tube under both laminar and turbulent flows using a constant heat flux boundary condition.Leong et al. [33] studied the entropy generation of TiO 2 -water and Al 2 O 3 -water nanofluids flows in a circular tube with a wall at constant temperature.Mahian et al. [34] conducted an investigation on the entropy generation due to nanofluids between two isoflux rotating cylinders and two types of nanofluids, Al 2 O 3 -EG and TiO 2 -water, were considered.Shahi et al. [35] conducted a numerical investigation on the entropy generation of Cu-water nanofluids in a square cavity with four different designs, where a heat source was arranged.Mahmoudi et al. [36] considered the entropy generation of Cu-water nanofluid natural convection in a cavity, in which three walls were set as adiabatic, and a constant temperature was employed on the right wall.Khorasanizadeh et al. [37] carried out an investigation on the entropy generation of Cu-water nanofluids in a square cavity where the top wall moved with an invariable velocity and the bottom wall was set as adiabatic.Khorasanizadeh et al. [38] also finished a research on the entropy generation inside a cavity where a baffle is arranged on the bottom hot surface.Esmaeilpour and Abdollahzadeh [39] reported a study on the entropy generation of Cu-water nanofluids natural convection in a cavity with wavy walls.Cho et al. [40] conducted an investigation on natural convection of water-based nanofluids in an enclosure with wavy walls, and three different types of nanoparticles, Cu, Al 2 O 3 and TiO 2 , are adopted, respectively.Boghrati et al. [41] achieved a numerical investigation on the entropy generation of the Al 2 O 3 -water and carbon nanotubes-water nanofluids flows through two horizontal parallel plates with a rectangular barrier.Sarkar et al. [42] conducted a study on the entropy generation of water-based Al 2 O 3 and Cu nanofluids mixed convection flowing past a square barrier in the middle of two parallel plates.Leong et al. [43] carried out an investigation on the heat transfer performance of Cu-water nanofluids in three types of shell-and-tube heat exchangers.
Matin et al. [44] accomplished a numerical research on MHD mixed convection flow using SiO 2 -water nanofluids over a non-linear stretching sheet.Selimefendigil et al. [45] carried out an investigation on entropy generation of nanofluids natural convection in entrapped trapezoidal cavities taking the influence of magnetic field into consideration.
According to Singh et al. [30], from the view point of decreasing entropy generation, the Al 2 O 3 -water nanofluids behaves better in conventional channels and minichannels under laminar flow, and similar behavior can be found in minichannels and microchannels under turbulent flow, which indicates the significant influences of nanofluids flow regime and channel size on entropy generation.Furthermore, dimples and protrusions are two excellent types of passive flow control turbulator, which are widely used in the current heat transfer augment technique and have great heat transfer enhancement capacity without excessive additional flow resistance penalty.In the present paper, the heat transfer performance and flow characteristics of Al 2 O 3 -water nanofluid at four different volume fractions in a rectangular conventional channel with dimples and protrusions are numerically investigated.The different flow regimes are considered including laminar flow, turbulent flow and transitional flow.The effects of the nanofluid volume fraction and dimples/protrusions on entropy generation and heat transfer performance are obtained.Both heat transfer entropy generation and friction entropy generation are taken into consideration.Besides, the effects on average wall temperature, wall temperature uniformity and pumping power are also analyzed to provide assessments of the thermal performance.

Physical Model and Numerical Method
In this study, the heat transfer performance and entropy generation characteristic of the Al 2 O 3 -water nanofluid flow in a conventional rectangular channel with staggered dimples or protrusions are numerically studied.Figure 1 shows the schematic diagram of the flow domain in this paper.The x, y, z represent the streamwise direction, spanwise direction and normal direction, respectively.The whole calculation domain consists of three parts, and the investigated part that is the middle heated section of the domain has a length of L = 245 mm.An inlet extension and an outlet extension, both of which are long enough, are set at the inlet and outlet, respectively, to guarantee a fully developed flow and avoid the outlet effect.The width W and height H of the flow channel are 125 mm and 10 mm, respectively.The staggered dimples or protrusions are arranged on one side of the channel.The detailed geometric parameters of dimples and protrusions are given in Figure 2. As two of the most important parameters of dimple or protrusion structure, the print diameter D = 30 mm while the ratio of depth-to-print diameter δ/D is 0.2.All dimples or protrusions are uniformly distributed on the heated surface, and the opposite side is also heated with a constant flux of 50,000 W/m 2 while the other walls are adiabatic.All the walls are set to be nonslip boundary in the computation.The flow Reynolds number Re ranges from 1000 to 40,000.Fully developed velocity boundary is employed on inlet, and the outlet pressure is set as the atmospheric pressure.The detailed information about the turbulator structures and nanofluid volume fractions of all cases are shown in Table 1.
Entropy 2016, 18, 148 4 of 19 exchangers.Matin et al. [44] accomplished a numerical research on MHD mixed convection flow using SiO2-water nanofluids over a non-linear stretching sheet.Selimefendigil et al. [45] carried out an investigation on entropy generation of nanofluids natural convection in entrapped trapezoidal cavities taking the influence of magnetic field into consideration.According to Singh et al. [30], from the view point of decreasing entropy generation, the Al2O3-water nanofluids behaves better in conventional channels and minichannels under laminar flow, and similar behavior can be found in minichannels and microchannels under turbulent flow, which indicates the significant influences of nanofluids flow regime and channel size on entropy generation.Furthermore, dimples and protrusions are two excellent types of passive flow control turbulator, which are widely used in the current heat transfer augment technique and have great heat transfer enhancement capacity without excessive additional flow resistance penalty.In the present paper, the heat transfer performance and flow characteristics of Al2O3-water nanofluid at four different volume fractions in a rectangular conventional channel with dimples and protrusions are numerically investigated.The different flow regimes are considered including laminar flow, turbulent flow and transitional flow.The effects of the nanofluid volume fraction and dimples/protrusions on entropy generation and heat transfer performance are obtained.Both heat transfer entropy generation and friction entropy generation are taken into consideration.Besides, the effects on average wall temperature, wall temperature uniformity and pumping power are also analyzed to provide assessments of the thermal performance.

Physical Model and Numerical Method
In this study, the heat transfer performance and entropy generation characteristic of the Al2O3-water nanofluid flow in a conventional rectangular channel with staggered dimples or protrusions are numerically studied.Figure 1 shows the schematic diagram of the flow domain in this paper.The x, y, z represent the streamwise direction, spanwise direction and normal direction, respectively.The whole calculation domain consists of three parts, and the investigated part that is the middle heated section of the domain has a length of L = 245 mm.An inlet extension and an outlet extension, both of which are long enough, are set at the inlet and outlet, respectively, to guarantee a fully developed flow and avoid the outlet effect.The width W and height H of the flow channel are 125 mm and 10 mm, respectively.The staggered dimples or protrusions are arranged on one side of the channel.The detailed geometric parameters of dimples and protrusions are given in Figure 2. As two of the most important parameters of dimple or protrusion structure, the print diameter D = 30 mm while the ratio of depth-to-print diameter δ/D is 0.2.All dimples or protrusions are uniformly distributed on the heated surface, and the opposite side is also heated with a constant flux of 50,000 W/m 2 while the other walls are adiabatic.All the walls are set to be nonslip boundary in the computation.The flow Reynolds number Re ranges from 1000 to 40,000.Fully developed velocity boundary is employed on inlet, and the outlet pressure is set as the atmospheric pressure.The detailed information about the turbulator structures and nanofluid volume fractions of all cases are shown in Table 1.The 3D Navier-Stokes (NS) equations in steady form are solved through the finite-volume based computational fluid dynamics solver CFX (version 15.0) to obtain the flow and heat transfer characteristic.The assumption that the flow is three-dimensional, steady, and incompressible is employed and the constant physical property is adopted.The resulting conservation equations of mass, momentum, and energy are given as follows (2 ) The NS equations are solved with implicit coupling method, the second order central differential scheme is employed in the diffusion term discretization and the advection discretization is achieved using high resolution scheme.The residues of continuity, energy and velocities are monitored to provide the judgment of the convergence of the computation, where the convergence criteria is set as 1 × 10 −6 .Besides, the temperature and velocity at the monitoring point are considered to judge the convergence of the computations in addition.
Considering the present Reynolds number range, which may cover laminar flow, transitional flow and turbulent flow, it is greatly appropriate to employ the shear stress transport (SST) turbulence model [46] coupled with Gamma-Theta transition model [47] as an effective approach to deal with the transitional flow in the present study.The transport equations of intermittency γ and transition momentum thickness Reynolds number Reθt are given as follows where Pγ1 and Eγ1 are the transition sources, and Pγ2 and Eγ2 are the destruction sources.μ and μt are the molecular viscosity coefficient and eddy viscosity coefficient.Pθt stands for the source term of the transition momentum thickness Reynolds number.The production terms of the turbulent kinetic energy downstream of the transition point is turned on by the employment of the intermittency.The empirical correlations are induced by the transition momentum thickness  The 3D Navier-Stokes (NS) equations in steady form are solved through the finite-volume based computational fluid dynamics solver CFX (version 15.0) to obtain the flow and heat transfer characteristic.The assumption that the flow is three-dimensional, steady, and incompressible is employed and the constant physical property is adopted.The resulting conservation equations of mass, momentum, and energy are given as follows ) The NS equations are solved with implicit coupling method, the second order central differential scheme is employed in the diffusion term discretization and the advection discretization is achieved using high resolution scheme.The residues of continuity, energy and velocities are monitored to provide the judgment of the convergence of the computation, where the convergence criteria is set as 1 ˆ10 ´6.Besides, the temperature and velocity at the monitoring point are considered to judge the convergence of the computations in addition.
Considering the present Reynolds number range, which may cover laminar flow, transitional flow and turbulent flow, it is greatly appropriate to employ the shear stress transport (SST) turbulence model [46] coupled with Gamma-Theta transition model [47] as an effective approach to deal with the transitional flow in the present study.The transport equations of intermittency γ and transition momentum thickness Reynolds number Re θt are given as follows where P γ1 and E γ1 are the transition sources, and P γ2 and E γ2 are the destruction sources.µ and µ t are the molecular viscosity coefficient and eddy viscosity coefficient.P θt stands for the source term of the transition momentum thickness Reynolds number.The production terms of the turbulent kinetic energy downstream of the transition point is turned on by the employment of the intermittency.The empirical correlations are induced by the transition momentum thickness Reynolds number; besides, the influences of the turbulence kinetic energy in the freestream and adverse pressure gradient are also captured.The transition model is realized based on the local variables, which can be easily achieved with the help of current CFD methods.Furthermore, the transition model is able to provide predictions of various transition processes using proper correlations based on experimental data, and the simulation results will be more reliable and accurate with the improvement of the correlation.The transition model coupled with SST turbulence model is adopted in this study.

Physical Property of Nanofluid
In the present study, the Al 2 O 3 -water nanofluids consisting of water as base fluid and Al 2 O 3 spherical nanoparticle with 30 nm diameter is employed.It is assumed that there is no motion slip between the continuous liquid and the discontinuous phase of the dispersed nanoparticles.Besides, it is also adopted that the local thermal equilibrium between base fluid and nanoparticles is achieved.Furthermore, the nanofluid is assumed to be in single phase [48] to evaluate the effective physical properties, including density, specific heat, dynamic viscosity and thermal conductivity.The effective models are given by: ‚ Density model ρ: ‚ Specific heat model C p [49]: ‚ Brinkman viscosity model µ [7], which is the most common relation to calculate the viscosity in the entropy generation problems: ‚ Bruggeman thermal conductivity model k [11] that can be applied to spherical particle with various concentrations: where In this paper, four different volume fractions are employed: 0%, 3%, 6% and 9%, to investigate the influence of volume fraction.The effective physical property models above are valid within the volume fraction range considered in this study.

Data Reduction
The flow Reynolds number Re in the present paper is defined by where U m,in stands for the average inlet velocity.The hydraulic diameter D h is defined as Entropy 2016, 18, 148 7 of 19 The pumping power P.P can be calculated as follows: where ∆P is the pressure drop between the inlet and outlet of the flow domain, and M represents the mass flow rate.The temperature uniformity ∆T s is given as follow, indicating the temperature distribution uniformity on the heated surface where T w,max and T w,min are the maximum and minimum temperature of the heated wall, respectively.The heated wall average temperature T w is defined by where A represents the area of the heated wall.Besides, the Nusselt number Nu and flow resistance coefficient f are given as follows, which are employed in the method validation and grid independence study section where h is the heat transfer coefficient and ∆P/L is the pressure gradient along the streamwise.
According to the relevant published literature, the local entropy generation rate S 3 gen in Cartesian coordinate systems can be obtained using the following relations [50] where the first term S 3 h and second term S 3 f on the right hand stand for the local entropy generation rates due to the heat transfer and flow friction, respectively.Through the numerical simulation, the local entropy generation rate can be obtained at every point in the flow domain.Furthermore, the total entropy generation can be obtained through the integral of the local entropy generation rate distribution over the whole domain.It can be found from Equation ( 17) that thermal conductivity and viscosity have significant influences on the entropy generation.The physical properties of nanofluids should be used in the equation above for the present study.Then, the average total entropy generation rate S rgen is given by where S rh and S rf are the average heat transfer entropy generation rate and average friction entropy generation rate, respectively.The dimensionless local entropy generation rate S 3 dgen is defined as where S

Method Validation and Grid Independence Study
In order to guarantee the validity of the method in the present study, a numerical method validation is completed through the comparison between the results in this study and other published results, as shown in Figure 3.It is shown that, considering the Re range in which Gnielinski [51] is valid, the present Nusselt numbers are in agreement with those based on Gnielinski when Re is no less than 3000.Besides, the Nusselt numbers at Re = 1000 and 3000 are in accordance with that of Shah [52].It is found from Figure 3b that the f values in the present study agree with those from Filonenko [53].In conclusion, the present method is deemed to be valid for predicting the flow and heat transfer in this study.Moreover, Gnielinski and Shah Equations are given as follows, as well as Filonenko Equation.

‚
Gnielinski: where c t " ´Pr f Pr w ¯0.01 , f is given as Filonenko Equation and L is the length of the domain considered.Pr f and Pr w are the Prandtl numbers based on the average fluid temperature and wall temperature, respectively.

Method Validation and Grid Independence Study
In order to guarantee the validity of the method in the present study, a numerical method validation is completed through the comparison between the results in this study and other published results, as shown in Figure 3.It is shown that, considering the Re range in which Gnielinski [51] is valid, the present Nusselt numbers are in agreement with those based on Gnielinski when Re is no less than 3000.Besides, the Nusselt numbers at Re = 1000 and 3000 are in accordance with that of Shah [52].It is found from Figure 3b that the f values in the present study agree with those from Filonenko [53].In conclusion, the present method is deemed to be valid for predicting the flow and heat transfer in this study.Moreover, Gnielinski and Shah Equations are given as follows, as well as Filonenko Equation.


Gnielinski: where , f is given as Filonenko Equation and L is the length of the domain considered.Prf and Prw are the Prandtl numbers based on the average fluid temperature and wall temperature, respectively. Filonenko:  Shah: All hexahedral mesh is adopted in this study and the y + in all cases is less than 1 to guarantee the computation precision.The O-type block mesh is employed near the dimples and protrusions in order to improve the local mesh quality.In order to achieve a balance between the calculation precision and computational resource, a grid independence study is carried out to obtain the reasonable grid for computation.As is shown in Table 2, the relative discrepancies of Nu and f are All hexahedral mesh is adopted in this study and the y + in all cases is less than 1 to guarantee the computation precision.The O-type block mesh is employed near the dimples and protrusions in order to improve the local mesh quality.In order to achieve a balance between the calculation precision and computational resource, a grid independence study is carried out to obtain the reasonable grid for computation.As is shown in Table 2, the relative discrepancies of Nu and f are 0.566% and 0.603% when the mesh changes from Mesh 3 to Mesh 4, respectively.Therefore, the proposed mesh is Mesh 3.

Results and Discussion
The Reynolds number in this paper ranges from 1000 to 40,000, which covers laminar flow, transitional flow and turbulent flow.After checking the turbulence intermittency of all cases, it is found that when Re reaches 10,000 or higher, the flow has successfully turned into fully developed turbulent flow.Moreover, the variation tendencies of entropy generation and other relevant parameters in the cases with Re at 10,000 or higher show no aberrations or unpredictable changes.In other words, the variation tendency of the considered parameters when Re is no less than 10,000 is in correspondence with that when Re ranges 1000 from 10,000.Thus, for a more distinct understanding of the results and discussion, this section provides the illustration and explanation of the cases with Re from 1000 to 10,000.

Flow Regimes and Flow Structures
Considering the focus of this paper-the entropy generation and heat transfer characteristics near the transitional flow, it is necessary to give an estimate about the flow regimes at different Reynolds number based on the simulation results.Taking the cases with nanofluids volume fraction at 0% as examples, the turbulence intermittency γ distributions contour and streamlines on the middle plane vertical to the spanwise are shown in Figure 4, where the legend levels are mandatorily set from 0 to 1 and the mainstream flows from the left to right.When γ equals to 0, the flow regime is considered to be laminar flow, and it is supposed to be turbulent flow when γ is 1.The flow with γ between 0 and 1 is defined as the transitional flow.It is shown that when Re = 1000, the mainstream is principally considered to be under laminar flow, and with the growth in Re, the turbulivity gets more intensified especially at the center of dimples and back porch of protrusions.When Re = 5000, the turbulent flow almost occupies the major part of the domain and the laminar flow only exists at the flow boundary layers neat the wall, leading to a laminar-layer.With the further increase in Re, the thickness of the laminar-layer keeps continuously declining.It can be concluded that when Re reaches 10,000 or even higher, the mainstream has completely turned into turbulent flow, and the Reynolds number range in this study covers laminar flow, transitional flow and turbulent flow.The cases with Reynolds numbers higher than 10,000 are additionally accomplished to study the influence of higher Re, which may better meet the practical industrial requirement.
With regard to the flow structures, a large flow separation region occurs inside the dimple when Re = 1000 which occupies the major section of the dimple, and the core of the vortex gradually moves towards the dimple center with the rise in Re.In protrusion cases, there is a large separation region at the protrusion trailing edge and the reattachment turns up at the leading edge of the adjacent protrusion.When Re ascends, the flow separation region size slightly increases.There exists a dense-streamline region near the protrusion top, indicating a local high velocity region.

Effect of Nanofluids on Entropy Generation
The variation tendencies of average heat transfer entropy generation rate Srh and average friction entropy generation rate Srf with Re, as well as the average total entropy generation rate Srgen, are shown in Figure 5a,c,e.Furthermore, the entropy generation differences of Srh, Srf and Srgen due to the nanofluid with Re are given in Figure 5b,d,f, respectively.It can be found in Figure 5a that Srh values in all cases observably decrease with the increase in Re, and the slopes get increasingly lower when Re rises, indicating the weaker influence to Srh of Re when flow velocity ascends.Furthermore, from the viewpoint of the effect of nanofluids, it can apparently be seen that the employment of nanofluid markedly reduces the Srh values in all Reynolds numbers in both dimple-case and protrusion-case.The Srh value keeps decreasing at a given Reynolds number with an increasing volume fraction.It should be noted that the decrease of the Srh value owing to nanofluids declines with the increase in Re.Taking the Srh value differences between the cases with φ = 0% and φ = 9% as examples (Figure 5b), the decrease extent dramatically descends with the growth in Re, from about 550 when Re = 1000 to about 100 when Re = 10,000 for dimple-case, which keeps dropping when Re further grows.The results above can be explained by that due to the high thermal conductivity of nanofluid, the physical property has great effects on the heat transfer performance when the mainstream keeps in laminar flow with low Re, which signally enhances the local heat transfer and decrease the local temperature gradient, leading to a lower heat transfer entropy generation rate.With the increase in Re, the mainstream gradually turns into turbulent flow, and the turbulence acts as the predominant impetus in heat transfer enhancement.The effect of the nanofluid thermal conductivity is not as significant as that when Re is low.As a result, the effect of nanofluids relatively weakens.

Effect of Nanofluids on Entropy Generation
The variation tendencies of average heat transfer entropy generation rate S rh and average friction entropy generation rate S rf with Re, as well as the average total entropy generation rate S rgen , are shown in Figure 5a,c,e.Furthermore, the entropy generation differences of S rh , S rf and S rgen due to the nanofluid with Re are given in Figure 5b,d,f, respectively.It can be found in Figure 5a that S rh values in all cases observably decrease with the increase in Re, and the slopes get increasingly lower when Re rises, indicating the weaker influence to S rh of Re when flow velocity ascends.Furthermore, from the viewpoint of the effect of nanofluids, it can apparently be seen that the employment of nanofluid markedly reduces the S rh values in all Reynolds numbers in both dimple-case and protrusion-case.The S rh value keeps decreasing at a given Reynolds number with an increasing volume fraction.It should be noted that the decrease of the S rh value owing to nanofluids declines with the increase in Re.Taking the S rh value differences between the cases with φ = 0% and φ = 9% as examples (Figure 5b), the decrease extent dramatically descends with the growth in Re, from about 550 when Re = 1000 to about 100 when Re = 10,000 for dimple-case, which keeps dropping when Re further grows.The results above can be explained by that due to the high thermal conductivity of nanofluid, the physical property has great effects on the heat transfer performance when the mainstream keeps in laminar flow with low Re, which signally enhances the local heat transfer and decrease the local temperature gradient, leading to a lower heat transfer entropy generation rate.With the increase in Re, the mainstream gradually turns into turbulent flow, and the turbulence acts as the predominant impetus in heat transfer enhancement.The effect of the nanofluid thermal conductivity is not as significant as that when Re is low.As a result, the effect of nanofluids relatively weakens.On the contrary, Srf value exhibits an opposite variation tendency.It can be found from Figure 5c that the Srf values in all cases observably ascend with the increase in Re, and the Srf values rise faster when Re is higher, which suggests that the influence on Srf of Re gets increasingly significant with higher Re.From the view point of nanofluid effect, the employment of nanoparticles increases the Srf values in both dimple-case and protrusion-case at a given Re, and Srf ascends with an increasing volume fraction, which is mainly caused by the higher viscosity compared with the base fluid and the accompanying higher velocity gradients.Besides, the effect gets more significant with the growth in Re.As shown in Figure 5d, the Srf value increases due to nanofluids dramatically ascend when Re rises.It is caused by that with the increase in Re, the influence of the viscosity difference between the nanofluids with various volume fractions gradually gets more significant.Considering that the increase in Srf value is not comparable to the decrease in Srh value, as a consequence, the average total entropy generation rate Srgen expresses a variation tendency with high similarity with that of Srh as shown in Figure 5e.It is concluded that the Srgen values dramatically decrease when Re rises or nanofluid volume fraction increases, and the effect of nanofluid on Srgen gradually becomes inconspicuous with the growth in Re as shown in On the contrary, S rf value exhibits an opposite variation tendency.It can be found from Figure 5c that the S rf values in all cases observably ascend with the increase in Re, and the S rf values rise faster when Re is higher, which suggests that the influence on S rf of Re gets increasingly significant with higher Re.From the view point of nanofluid effect, the employment of nanoparticles increases the S rf values in both dimple-case and protrusion-case at a given Re, and S rf ascends with an increasing volume fraction, which is mainly caused by the higher viscosity compared with the base fluid and the accompanying higher velocity gradients.Besides, the effect gets more significant with the growth in Re.As shown in Figure 5d, the S rf value increases due to nanofluids dramatically ascend when Re rises.It is caused by that with the increase in Re, the influence of the viscosity difference between the nanofluids with various volume fractions gradually gets more significant.Considering that the increase in S rf value is not comparable to the decrease in S rh value, as a consequence, the average total entropy generation rate S rgen expresses a variation tendency with high similarity with that of S rh as shown in Figure 5e.It is concluded that the S rgen values dramatically decrease when Re rises or nanofluid volume fraction increases, and the effect of nanofluid on S rgen gradually becomes inconspicuous with the growth in Re as shown in Figure 5f.Furthermore, it should be paid great attention to that with the increase in Re, the difference between S rh value and S rf value gets lower, so does the influences of nanofluids on these two values.Consequently, it can be predicted that the in the flow and heat transfer with Re higher than that in this paper, called critical Re here, the S rf value will play a role as important as S rh value in the effect on average total entropy generation rate.The detailed study on the critical Re will be considered in the further investigation.

Effect of Dimples/Protrusions on Entropy Generation
As to the comparison between the influences of dimple-case and protrusion-case, it can be seen from Figure 5a that the protrusion-case' S rh values are much lower than those of dimple-case, especially in low Re.The difference between the S rh values of protrusion-case and dimple-case gradually declines when Re rises, which is mainly caused by that in laminar flow, the protrusion-case's heat transfer enhancement capacity is markedly stronger than that of dimple-case, leading to the more uniform temperature field and lower temperature gradient.Consequently, the S rh value is much lower.With the growth in Re, the turbulence's effect on heat transfer becomes increasingly remarkable and the influences of dimples and protrusions are relatively reduced.With regard to S rf , it can be found in Figure 5c that the S rf values of protrusion-case are distinctly higher than those of dimple-case at given Reynolds numbers, of which the difference increasingly ascends with the rise in Re, which is considered to be the result of that when Re rises, the effect of protrusion-case's flow separation on the flow field grows and the large flow separation region brings about the local low velocity region and the accompanying increase in local friction entropy generation rate.
It should be noted that after checking the results of all cases, the S rh values and S rf values in dimple-case are markedly lower than S rh0 and S rf0 at given Reynolds numbers, respectively, suggesting that dimple-case has great benefit in decreasing entropy generation compared to smooth channel both in heat transfer entropy generation rate and friction entropy generation rate, which primarily owes to the heat transfer enhancement by dimples without significant additional flow resistance.Besides, the protrusion-case also provides notable advantage in diminishing heat transfer entropy generation, of which the mechanism is similar to that in dimple-case.Nevertheless, the S rf values in protrusion-case are considerably higher than S rf0 , which indicates that when compared with smooth channel, protrusion-case suffers from penalty of distinct friction entropy generation rate increase though it performs well in heat transfer entropy generation rate decline.
Furthermore, for a better understanding, the contours of S dh layer in protrusion-case is evidently lower than that in the dimple-case, which is caused by the acceleration of the mainstream near the protrusion leading edge and the local reattachment.As shown in Figure 6c, part of the speeded mainstream intensively impinges on the opposite smooth wall and then the developments of the flow and heat boundary layers are both strongly disturbed.The local heat transfer is enhanced and lower temperature gradients are achieved.Nevertheless, the disturbance effect above in the dimple-case not as significant as that in protrusion-case, bringing about a thicker high S

Heat Transfer Performance Analysis
Taking the cases with Re = 3000 as examples for the heat transfer performance analysis, Figure 7 provides the wall temperature distributions and limit streamlines on the heated walls in dimple-case and protrusion-case, and the cases with φ = 0% and φ = 9% are selected to illustrate the effect of nanofluids, where the mainstream flows from the left to right.In dimple-case, it is shown that two symmetrical large vortices occur inside the dimple, leading to two symmetrical high temperature regions in agreement with the vortices' locations.The reattachment on the dimple trailing edge significantly enhances the local heat transfer, resulting in a large low temperature region.With the employment of nanofluids, the whole temperature distribution pattern hardly alters.Nevertheless, it is quite obvious that the high temperature regions' area notably decreases and the low temperature regions' area markedly increases.Moreover, the entire temperature level prominently declines almost at all regions when nanofluids is adopted, indicating the great enhancement of heat transfer, especially on the dimple trailing edge.The wall temperature distribution on the opposite smooth wall is also given in Figure 7c,d.It is found that there is a notable decrease in the area of the high temperature regions with the employment of nanofluid.In protrusion-case, the low temperature regions are mainly located at the protrusion leading edges, and the high temperature regions are substantially detected near the flow separation regions between the two rows of protrusions.The acceleration and impinging enhance the local heat

Heat Transfer Performance Analysis
Taking the cases with Re = 3000 as examples for the heat transfer performance analysis, Figure 7 provides the wall temperature distributions and limit streamlines on the heated walls in dimple-case and protrusion-case, and the cases with φ = 0% and φ = 9% are selected to illustrate the effect of nanofluids, where the mainstream flows from the left to right.In dimple-case, it is shown that two symmetrical large vortices occur inside the dimple, leading to two symmetrical high temperature regions in agreement with the vortices' locations.The reattachment on the dimple trailing edge significantly enhances the local heat transfer, resulting in a large low temperature region.With the employment of nanofluids, the whole temperature distribution pattern hardly alters.Nevertheless, it is quite obvious that the high temperature regions' area notably decreases and the low temperature regions' area markedly increases.Moreover, the entire temperature level prominently declines almost at all regions when nanofluids is adopted, indicating the great enhancement of heat transfer, especially on the dimple trailing edge.The wall temperature distribution on the opposite smooth wall is also given in Figure 7c,d.It is found that there is a notable decrease in the area of the high temperature regions with the employment of nanofluid.In protrusion-case, the low temperature regions are mainly located at the protrusion leading edges, and the high temperature regions are substantially detected near the flow separation regions between the two rows of protrusions.The acceleration and impinging enhance the local heat transfer, which is deteriorated near the vortices.When nanofluids is employed, it is seen that the high temperature regions' area declines and the low temperature regions' area grows.The analogous effect also takes place on the opposite smooth wall, and the wall temperature level markedly drops.For a more comprehensive analysis, the variation tendencies of the average wall temperature T w , wall temperature uniformity ∆T s and pumping power P.P versus Re are provided in Figure 8.It can be found in Figure 8a that at a given Re, T w keeps decreasing with the increase in nanofluids volume fraction in both dimple-case and protrusion-case, which is primarily due to the higher effective thermal conductivity.In other words, higher volume fraction brings about more effective cooling.Besides, the effect of nanofluids gradually gets weak when Re rises.In industrial application, the high temperature rise results in the nonuniform temperature distribution on the equipment that may bring about the damaging thermal stress as the consequence of the difference in thermal expansion coefficient.Besides, the spatial temperature gradient also has a negative effect on the equipment's stable operation.It is shown in Figure 8b that with the growth in Re, the temperature uniformity ∆T s notably declines in both dimple-case and protrusion-case.Besides, the temperature uniformity also decrease with an increasing nanofluid volume fraction, indicating the improvement of temperature uniformity with the employment of nanofluid.It is found from Figure 8c that with the increase in volume fraction, the pumping power P.P ascends in both dimple-case and protrusion-case, which is mainly caused by the rise in the effective viscosity that results in a larger pressure drop.Moreover, the pumping power also grows with an increasing Re as expected.For a more comprehensive analysis, the variation tendencies of the average wall temperature Tw, wall temperature uniformity ΔTs and pumping power P.P versus Re are provided in Figure 8.It can be found in Figure 8a that at a given Re, Tw keeps decreasing with the increase in nanofluids volume fraction in both dimple-case and protrusion-case, which is primarily due to the higher effective thermal conductivity.In other words, higher volume fraction brings about more effective cooling.Besides, the effect of nanofluids gradually gets weak when Re rises.In industrial application, the high temperature rise results in the nonuniform temperature distribution on the equipment that may bring about the damaging thermal stress as the consequence of the difference in thermal expansion coefficient.Besides, the spatial temperature gradient also has a negative effect on the equipment's stable operation.It is shown in Figure 8b that with the growth in Re, the temperature uniformity ΔTs notably declines in both dimple-case and protrusion-case.Besides, the temperature uniformity also decrease with an increasing nanofluid volume fraction, indicating the improvement of temperature uniformity with the employment of nanofluid.It is found from Figure 8c that with the increase in volume fraction, the pumping power P.P ascends in both dimple-case and protrusion-case, which is mainly caused by the rise in the effective viscosity that results in a larger pressure drop.Moreover, the pumping power also grows with an increasing Re as expected.

Conclusions
In this study, effect of Al2O3-water nanofluid flow on entropy generation and heat transfer performance in a rectangular conventional channel with staggered dimples and protrusions is numerically investigated.Some important conclusions from the results can be drawn as follows.

Conclusions
In this study, effect of Al 2 O 3 -water nanofluid flow on entropy generation and heat transfer performance in a rectangular conventional channel with staggered dimples and protrusions is numerically investigated.Some important conclusions from the results can be drawn as follows.

‚
When Re reaches 10,000 or even higher, the mainstream has completely turned into turbulent flow, and Re range in this study covers laminar flow, transitional flow and turbulent flow.

‚
The S rh keeps decreasing at a given Re with an increasing volume fraction, and the decrease of the S rh owing to nanofluids declines with the increase in Re.The employment of nanofluid increases S rf at a given Re, and S rf further ascends with an increasing volume fraction, in which the effect gets more significant with the growth in Re.The S rgen expresses a variation tendency with high similarity with that of S rh under the Re range in this study.

‚
The dimple-case has great benefit in decreasing entropy generation compared to smooth channel both in heat transfer entropy generation rate and friction entropy generation rate.The protrusion-case suffers from penalty of distinct friction entropy generation rate increase though it performs well in heat transfer entropy generation rate decline.The high temperature regions' area notably decreases and the low temperature regions' area markedly increases with an increasing volume fraction.The entire temperature level prominently declines almost at all regions in both the structured and smooth surfaces.

‚
With the increase in nanofluids volume fraction, T w keeps decreasing, ∆T s declines and P.P ascends in both dimple-case and protrusion-case at a given Re.

‚
Conclusively, it is preferred to employ nanofluids to obtain great heat transfer enhancement performance and significant entropy generation decrease under laminar flow.It is also recommended that protrusions be adopted under laminar flow for its superior heat transfer enhancement performance and remarkable entropy generation decrease.Besides, dimples are preferable under turbulent flow in order to avoid high pumping power and excessive friction entropy generation increase.
Author Contributions: Yonghui Xie and Lu Zheng performed the numerical simulations, finished the results analysis, and then wrote some sections of the manuscript.Di Zhang and Gongnan Xie prepared some other sections of the paper.All of the authors contributed equally for reviewing and revising the manuscript.All authors have read and approved the final manuscript.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript:

Figure 1 .
Figure 1.Schematic diagram and detailed information of the flow domain.

Figure 1 .
Figure 1.Schematic diagram and detailed information of the flow domain.

Figure 2 .
Figure 2. Detailed geometric parameters of dimples and protrusions.

Figure 2 .
Figure 2. Detailed geometric parameters of dimples and protrusions.

3 dh and S 3 df
are the dimensionless local heat transfer entropy generation rate and dimensionless local friction entropy generation rate, respectively.S rh0 and S rf0 are the average heat transfer entropy generation rate and average friction entropy generation rate of the smooth channel without dimples or protrusions.

Figure 3 .
Figure 3. Validation of the present method: (a) Nu vs. Re; and (b) f vs. Re.

Figure 3 .
Figure 3. Validation of the present method: (a) Nu vs. Re; and (b) f vs. Re.

Figure 5 .
Figure 5. Average entropy generation rate characteristic: (a) S rh vs. Re; (b) S rh decrease vs. Re; (c) S rf vs. Re; (d) S rf increase vs. Re; (e) S rgen vs. Re; and (f) S rgen decrease vs. Re.

3 dh and S 3 d
f are illustrated, here taking Re = 3000 as examples.Figure 6 provides the S 3 dh and S 3 d f distributions on the spanwise-vertical plane (abbreviated as Plane 1) and streamwise-vertical plane (abbreviated as Plane 2) in both dimple-case and protrusion-case when Re = 3000, where the mainstream flows from the left to right in Plane 1.It is shown that S 3 dh keeps at a quite low level in the major area of the domain, and the high S 3 dh region is mainly located at the boundary layer near the wall where high temperature gradients exist, resulting in a high S 3 dh layer.It can be found that the thickness of the high S 3

Figure 8 .
Figure 8.Average wall temperature, wall temperature uniformity and pumping power variation characteristics: (a) T w vs. Re; (b) ∆T s vs. Re; and (c) P.P vs. Re.

‚ The S 3 dh
keeps at a quite low level in the major area of the domain, and the high S 3 dh regions are mainly located near the boundary layer and the local vortices.There is a low S 3 d f region near dimple trailing edge, and high S 3 d f regions are detected in accordance with low velocity areas.‚

Table 1 .
Detailed information of the cases.

Table 1 .
Detailed information of the cases.