Thermal-Hydraulic performance in a microchannel heat sink Thermal-Hydraulic performance in a microchannel heat sink equipped with longitudinal vortex generators (LVGs) and nanofluid equipped with longitudinal vortex generators (LVGs) and nanofluid

: In this study, the numerical conjugate heat transfer and hydraulic performance of nanoﬂuids ﬂow in a rectangular microchannel heat sink (RMCHS) with longitudinal vortex generators (LVGs) was investigated at di ﬀ erent Reynolds numbers (200–1200). Three-dimensional simulations are performed on a microchannel heated by a constant temperature with ﬁve di ﬀ erent conﬁgurations with di ﬀ erent angles of attack for the LVGs under laminar ﬂow conditions. The study uses ﬁve di ﬀ erent nanoﬂuid combinations of Al 2 O 3 or CuO, containing low volume fractions in the range of 0.5% to 3.0% with various nanoparticle sizes that are dispersed in pure water, PAO (Polyalphaoleﬁn) or ethylene glycol. The results show that for Reynolds number ranging from 100 to 1100, Al 2 O 3 –water has the best performance compared with CuO nanoﬂuid with Nusselt number values between 7.67 and 14.7, with an associated increase in Fanning friction factor by values of 0.0219–0.095. For the case of di ﬀ erent base ﬂuids, the results show that CuO–PAO has the best performance with Nusselt number values between 9.57 and 15.88, with an associated increase in Fanning friction factor by 0.022–0.096. The overall performance of all conﬁgurations of microchannels equipped with LVGs and nanoﬂuid showed higher values than the ones without LVG and water as a working ﬂuid


Introduction
In the past few years, due to the fast development of micro/nanosystems in different engineering applications, like micro-electronics cooling applications, medical instruments, etc., the need to develop high performance compact thermal devices with high efficiency, minimum cost, lightweight and with the smallest size possible is necessary. Pioneered by Tuckerman and Pease [1], the concept of a rectangular microchannel heat sink made of silicon was first introduced. Since then active research on its thermal-hydraulic performance has been focused on by several investigators as reviewed by Adham et al. [2], Agostani et al. [3] and Morini. [4]. In general, liquid coolants have higher heat transfer coefficients as compared with gaseous coolants; that is why liquid coolants are most common to utilize in microchannel heat sinks [5]. Xu et al. [6] found that experimental results in the microchannel are compatible with the numerical results predicted by the conventional Navier-Stokes equation. Quand and Mudawar [7] found that the heat transfer properties of microchannel heat sinks can be predicted by conventional Navier-Stokes and energy equations. Harms et al. [8] found that with decrease of the width of the channel and increase of the depth of the channel the heat transfer increased for a range of Reynolds numbers between 173 and 12,900. Their study also showed that the critical Reynolds number was 1500. Tuckerman [9] reported that laminar flow is the best for heat removal through microchannels due to the development of a thin thermal boundary layer. But, in order to get higher performance, this barrier must be overcome with some surface geometry modification. In the past, thermal and hydraulic effects on the channels of trigonometric and reconstructed geometrical shapes were studied with circular cavities with rectangular ribs, triangular, and trapezoidal circles [10,11]. The significant improvement in heat transfer is due to the disruption of the boundary layer by ribs, which improves the mixing of the layer, leading to reduced heat transfer resistance [12,13].
In 1969 for the first time, Johnson and Joubert [14] investigated vortex generators' (VGs) effect on heat transfer performance. The generated pressure difference along the two faces of a vortex generator creates a flow-separation in the side edges and induced transverse, longitudinal and horseshoe vortices [14,15], which increase the heat transfer rate. When the angle of attack of the VG is small, the vortices generated is basically longitudinal. When the VG is perpendicular to the mainstream direction, the vortices generated will be transverse [16,17]. Yadav et al. [18] studied micro sinks with a cylindrical longitudinal vortex generator (LVG) and reported a significant effect on the heat transfer process. It was noticed that in comparison with transverse vortices, the longitudinal vortices are more efficient in increasing the heat transfer [12]. Fiebig et al. [19] in his experimental investigation found that channels equipped with VGs in the laminar regime showed remarkable improvement in local heat transfer by a factor of three relative to the one without VGs. They found that the critical Reynolds number could be reduced by a factor of 10 and more. Sohankar and Davidson [20] studied the effects of the size, tilt and position of VGs, using a pair of slant blocks in the channel, for a Re number between 400 and 1500. For the 30 • slants, the flow was unstable between a Re number of 1000 and 1500. They reported that with the increase in slant from 10 • to 30 • , the Nusselt number showed a sharp increase. Mohammed et al. [21] numerically studied the thermal and hydraulic properties of turbulent nanofluids in ribbed channels. They studied nine different shapes of the ribs, of which four Al 2 O 3 , CuO, SiO 2 and ZnO nanoparticles were distributed in a different set of base fluids (water, glycerol and motor oil). Reynolds numbers from 5000 to 20,000 increased the heat transfer. The results show that the polygonal rectangular groove has the best overall performance in all forms of polygonal sides, and the SiO 2 nanofluid has the best overall performance in all types of nanofluids.
Islam et al. [22] experimentally studied heat transfer performance in rib-roughened narrow millichannels for high heat flux system removal using water. The results showed that in the range of Reynolds number between 8000 and 30,000, heat transfer rates increased 2-2.5 times in comparison with a smooth channel, while the frictional pressure drop increased around 2.5 times. Leu et al. [23] and Wu and Tao [24] carried out experimental and numerical investigations on LVGs using water and air as working fluids, respectively. They found that, with various winglet attack angels, the attack angle of 45 • provides the best heat transfer performance, while the average Nusselt number was significantly enhanced with the increase of the attack angle. Liu et al. [25] experimentally employed a rectangular vortex generator with varying the number of VGs pairs and angles of attack in the rectangular microchannels to improve the heat transfer for Re numbers in the range between 170 and 1200. They reported that the transition zone was observed for Re numbers in the range from 600 to 720, while the best performance was with three pairs of LVG combination with angles of 150 • , 30 • and 150 • . They also found that in comparison to channels without LVGs, the range of critical Reynolds number decreases by adding LVGs.
Based on the work of Liu et al. [25], Chen et al. [26] used the LVG range of Re values between 350 and 1500. They found that the maximum heat transfer occurred at 45 • . An LVG with a height one-fourth of the length. Recently, Salleh et al. [27] experimentally investigated the heat transfer augmentation in a plain fin-and-tube heat exchanger (FTHE) with a trapezoidal winglet vortex generator in a Reynolds number range from 500 to 2500, as well as different angles of attack. They reported that the FTW (flat trapezoidal winglet) in CFU (common flow up) arrangement has the best overall performance across the FTHE channel but with the penalty of increasing pressure drop. Zli et al. [28] numerically studied, in two-dimensions, a flapping vortex generator employed for heat sinks to improve heat transfer. They used a vortex generator made of a thin flexible plate that was attached to the inner wall of the radiator channel at an oblique angle. They found that Young's vortex generator with a modulus of 1 MPa had the best performance among all three other options in the study and that it could improve heat removal efficiency by 140% at the same speed and the same total pump power. When the total extraction power is the same, the heat dissipation rate can be increased by 87%, and the average Nusselt average is increased by 200% compared with the ordinary channel at the same Reynolds number. Ebrahimi et al. [29] numerically studied the thermal and hydraulic properties of rectangular microchannel with LVG. He reported that with the use of LVG the heat transfer enhancement gets higher, but with a higher pressure drop.
Xie et al. [30] used Al 2 O 3 -nanoparticles in water to numerically study the entropy generation and heat transfer performance in traditional rectangular channels with dimples and protrusions. Four different-volume fractions and Reynolds numbers were used to cover laminar, transitional and turbulent flows, respectively. They found that the rate of average thermal entropy generation decreases, while the rate of average frictional entropy generation increases by increasing nanofluid volume fraction. Ebrahimi et al. [31] carried out the research and numerically studied the thermal and hydraulic properties of rectangular microchannel heat sinks with LVG using Al 2 O 3 or CuO nanoparticles dispersed in water. They reported that with the increase of the nanoparticle concentrations of nanofluids the heat transfer enhancement gets higher, but with a higher pressure drop.
Through the years, different working fluids have been used in microchannels such as water and oil, as well as in organic chemicals such as ethylene glycol. Still, due to the limited thermal conductivity of the conventional cooling substance, it is unable to keep with the requirement of rapid development in various applications. Nanofluids cooling technology can offer a solution because their higher thermal conductivities in comparison with conventional fluids. In 1995 Choi [32] proposed a theoretical model to improve the thermal properties of traditional liquids by using the two-component mixture model developed by Hamilton and Crosser [33] for the effective thermal conductivity, in which a dilute suspension, nano-sized particle dispersed in a base fluid, and thus called nanofluids. The nanoparticles in general include metals, metal-oxides, polymers, silica or even carbon nanotubes, while the base fluids cover water, oil or ethylene glycol. Brinkman [34] provided the most common expression used for the prediction of the viscosity of solutions and suspensions within specific concentrations and considered the effect of the additional solute molecule. Corcione [35] provided an empirical correlation through theoretically studying the heat transfer characteristics of buoyancy-driven nanofluids in a rectangular channel with heated vertical walls. Maxwell [36] proposed the first model to calculate the thermal conductivity of solid-liquid mixtures. Bruggeman [37] extended Maxwell's work by providing a thermal conductivity model to consider the interactions between the spherical particles. Yu and Choi [38] proposed an improved Hamilton-Crosser model that can correctly predict the thermal conductivity magnitude of nanotube-in-oil nanofluids for nonspherical particles. Machrafi and Lebon [39] provided a model for predicting the thermal conductivity of nanofluids, considering the intermediate layer between base liquids, nanoparticles, Brownian motion and particle complexes. Assessments were also provided to study the effect of each mechanism on thermal conductivity. PAO (Polyalphaolefin) oil-based nanoparticles seems to have great potential in enhancing the thermal conductivity. Shaikh et al. [40] found that by using nanofluids with nanocarbon fibers and exfoliated graphite dispersed in PAO oil, the thermal conductivity is improved by more than 100%. Lee et al. [41] experimentally studied the suspension of 4.0% volume 35 nm CuO particles in ethylene glycol and reported that the thermal conductivity was increased by 20%. Suresh et al. [42] experimentally studied the thermal and hydraulic properties of carbon nanofluids in laminar flat and dimpled tube distilled water. The results show that the use of nanofluid as the working fluid in a helically dimpled tube can improve the heat transfer performance, but the friction coefficient will increase slightly. Kalteh et al. [43] carried out a numerical study on heat transfer inside a microchannel with laminar flow and water-Cu nanofluid as a working fluid. They reported that with increasing the particle volume fraction and Reynolds number the heat transfer performance gets higher at the expense of a higher pressure drop. Gavili [44] performed a numerical study on heat transfer in a compact heat exchanger with VGs. He studied the effect of the angle of attack for VGs, and the position of the winglet on heat transfer for different Reynolds numbers. The results show that the position of the winglet has a huge impact on the enhancement of heat transfer performance of compact heat exchangers with the penalty of increasing pressure loss. Kalteh et al. [45] extended their previous work by numerically and experimentally studying the thermal-hydraulic performance of a microchannel heat sink with the use of water-alumina nanofluid as a working fluid. They reported that the temperature and velocity differences between the phases could be negligible.
As can be seen from the above literature review, the use of nanofluid together with LVG for different cooling applications in microchannels is still in its early stages. Thus, more studies are needed in order to get a better understanding of heat transfer and flow features of such a combination. This article aims to obtain a thermal-hydraulic analysis by applying a three-dimensional simulation in rectangular microchannels equipped with different LVGs configurations and different nanofluids as working fluids using validated and efficient thermal conductivity and viscosity models. That includes the effects of different LVG configurations, nanoparticle sizes, bulk temperature, the Brownian motion of particles with various nanofluids and various Reynolds numbers under laminar flow condition. To the best of the authors' knowledge, this is the first detailed analysis of the thermal-hydrodynamic performance of the combination of LVGs and nanofluids with different base fluids. A different sort of base fluid in microchannels is not yet addressed in the literature and may introduce a new point-of-view for novel methods of enhancing heat transfer in microchannels combined with LVGs.

Physical and Computational Model
Three-dimensional simulations were performed on five different configurations of a microchannel with rectangular winglet-type LVGs to examine the impact of using nanofluids on the conjugate heat transfer characteristics and flow field. Figure 1 shows the dimensions of the physical model with the relevant parameters of the microchannels while the geometry parameters for the different configurations are presented in detail in Table 1. All parameters in Table 1 are related to H, (H = 1), which is the height of the microchannel. LVGs are used in CFD (common flow down) for an arrangement with 30 • and CFU (common flow up) for arrangement with 150 • . The fluid flow was described by the Cartesian coordinate system in which the Z-axis direction is the stream-wise direction. The simulation consists of two parts, where the first part contains a microchannel with two pairs of LVGs, which are located in the heated region. The second part is a microchannel, including the adiabatic walls with the same weight of the original one but without VGs and its length is where it is used to develop a fully developed flow for the outlet while the velocity boundary condition on the inlet was set as a fully developed velocity profile. Because of symmetric arrangements of the microchannel and LVGs, only the hatched region that is shown in Figure 1 was included in the numerical simulations in order to decrease the computational time and costs.

Mathematical Models, Boundary Conditions and Governing Equations
The analysis presented in this study was performed on microchannel equipped with LVGs made of silicon, a conductive material, and based on the following assumptions: 1. The flow is considered as three dimensional, incompressible, hydraulically fully developed, steady and laminar due to the low velocity of the fluid and small pitch of the winglet. 2. The effects of radiation, gravity, body forces, viscous dissipation [1] and surface roughness are neglected [46]. 3. Working fluids were selected from different nanofluids with low volume-fraction and different sizes of Al2O3 or CuO nanoparticles dispersed in pure water, POA and ethylene glycol as base fluids. 4. A single-phase model was used to effectively describe the heat transfer behavior of nanofluids and validated based on low concentration of nanoparticles with diameters smaller than 100 nm [47]. 5. LVG was generated based on the quasi-steady phenomenon as reported by Ferrouillat et al. [48].   Because of symmetric arrangements of the microchannel and LVGs, only the hatched region that is shown in Figure 1 was included in the numerical simulations in order to decrease the computational time and costs.

Mathematical Models, Boundary Conditions and Governing Equations
The analysis presented in this study was performed on microchannel equipped with LVGs made of silicon, a conductive material, and based on the following assumptions: The flow is considered as three dimensional, incompressible, hydraulically fully developed, steady and laminar due to the low velocity of the fluid and small pitch of the winglet.

3.
Working fluids were selected from different nanofluids with low volume-fraction and different sizes of Al 2 O 3 or CuO nanoparticles dispersed in pure water, POA and ethylene glycol as base fluids.

4.
A single-phase model was used to effectively describe the heat transfer behavior of nanofluids and validated based on low concentration of nanoparticles with diameters smaller than 100 nm [47].
The coolant considered to be Newtonian and the thermophysical properties are dependent on the temperature, volume fraction and size of nanoparticles.
Based on the previous assumptions, the equations for steady-state flow of continuity, momentum and energy can be expressed as follows [49]: The energy equation for the solid zone is expressed as: As depicted in Figure 1, the boundary conditions are: • At the inlet: the velocity is fully developed.
• No-slip boundary condition is applied at the top, bottom and side walls as • Bottom and side walls are considered as adiabatic walls: • At the top surface, a uniform temperature is applied as • The conjugate heat transfer between LVG surfaces (solid) and nanofluid is applied as where n is a normal vector on the LVGs surface pointed out of the boundary. Table 2 shows the thermo-physical properties of pure water, POA, ethylene glycol, silicon, Al 2 O 3 and CuO nanoparticles. The nanofluids' effective thermophysical properties are calculated using the correlations illustrated below.
By the use of the mixture model approach the density of the nanofluid can be calculated as follows [50]: By taking the energy balance between the base fluid and nanoparticles the specific heat of nanofluid can be determined as Both dynamic viscosity and the effective thermal conductivity of nanofluids are determined by using a model developed by Koo and Kleinstreuer [57], which takes into consideration the effects of particle size, Brownian motion and particle volume fraction with different types of base fluid and nanoparticle combinations. This model is valid for a temperature range between 300 K and 325 K, while the volume fraction is up to 4% [16].
where k static can be optioned from Maxwell's model and k Brownian is developed based on both kinetic theory and micromixing method as shown in the relationship below [58]: where R f is Kapitza resistance, which equals to 4 × 10 −8 km 2 /W. By using the same approach, the effective viscosity of nanofluids can be predicted as follows: where g function is a function that considered the effects of particle diameter, temperature and volume fraction. For various based fluids and nanoparticles, the relation for Al 2 O 3 -water nanofluids and CuO-water nanofluids based on the experimental data with regression values of 98% and 96%, respectively, are shown below [59]: where a, b, c, d, e, m, h, I, j and k are coefficients related to the type of based fluid and nanoparticle.

Numerical Procedures and Parameter Definitions
Commercial software for computational fluid dynamics ANSYS Fluent Software Release 16.0 was used to perform the simulations. Using ANSYS Meshing the meshes of the non-uniform grid for the computational domains. The grid generation of the model is shown in details in Figures 2 and 3. The COUPLED algorithm was used for the pressure-velocity coupling based on the following reasons: First, by using it the approximations associated were eliminated in comparison to the segregated solution method, such as SIMPLE and SIMPLEC. Secondly, it reduces the errors that caused and related to physical model nonlinearities, the initial conditions and meshes that have been stretched and skewed [60]. The solution converges when the residual values are less than 10 −6 for the momentum equation, 10 −6 for the continuity equation and 10 −8 for the energy equation.
where a, b, c, d, e, m, h, I, j and k are coefficients related to the type of based fluid and nanoparticle.

Numerical Procedures and Parameter Definitions
Commercial software for computational fluid dynamics ANSYS Fluent Software Release 16.0 was used to perform the simulations. Using ANSYS Meshing the meshes of the non-uniform grid for the computational domains. The grid generation of the model is shown in details in Figures 2-3. The COUPLED algorithm was used for the pressure-velocity coupling based on the following reasons: First, by using it the approximations associated were eliminated in comparison to the segregated solution method, such as SIMPLE and SIMPLEC. Secondly, it reduces the errors that caused and related to physical model nonlinearities, the initial conditions and meshes that have been stretched and skewed [60]. The solution converges when the residual values are less than 10 −6 for the momentum equation, 10 −6 for the continuity equation and 10 −8 for the energy equation.   The following parameters are defined to make a framework for the numerical results of this study.
• Reynolds number: The following parameters are defined to make a framework for the numerical results of this study.

Grid Independency and Model Validation
In order to keep the simulation accurate with reasonable computational costs, a grid independence test was performed on four grids with different sizes as shown in Table 3. The tests were performed on A 1 configuration with β1, β2 angle of 30 • each using pure water as a working fluid under Re = 800. It can be seen from the result shown in Table 3 that there is a small difference between (490198) grid and (931236) grid in the computed results which was 0.038 and 0.034, respectively, hence the (490198) grid was selected for further simulations. The results were validated based on the experimental results presented by Liu et al. [25]. They used in their study deionized water was a working fluid with Reynolds number in the range of 170 to 1200. The bottom wall of the microchannel was maintained at a constant temperature of 323 K. Figures 4 and 5 illustrate a comparison between Nu and f values from the present numerical simulation and the experimental results of Liu et al. [25]. It can be seen from these figures that the maximum and minimum deviations were 11.5% and 3.8%, respectively, for the Nu values while the maximum and minimum deviations were 13.34% and 2.7%, respectively, for the values of friction factor. Liu et al. [25] reported that the average uncertainties are caused by the equipment that has been used in determining the Nu number for the microchannel with and without LVGs is 19.6% and 12.8%, respectively.

Results and Discussions
The numerical computations were conducted for five different geometrical configurations on the thermo-hydraulic performance of microchannels with VGs and compared to the smooth microchannels using an Al 2 O 3 -water nanofluid as a working fluid. In this case, the effects of the nanofluids with different nanoparticle types (Al 2 O 3 and CuO) and different base fluid (water, PAO and ethanol glycol) with a nanoparticles volume fraction of 2% and various nanoparticles diameters (29 nm and 28.6 nm) were tested using A 1 configuration and compared to pure water. The effects of different geometry configurations and Reynolds numbers on the temperature and velocity distributions, Nusselt number, fanning friction coefficient and JF performance factor for different Reynolds numbers are interpreted and presented in this section.

The Effect of Different Geometry Configurations
The temperature distribution on different cross-sections along the stream-wise direction for A 1 configurations at Re = 800 is shown in Figure 6. It is observed that the interactions between the working fluid and hot surfaces and longitudinal vortex generator cause a change in the temperature distribution. By increasing the Reynolds number, the thermal boundary layer on the microchannel walls and the area near the sides of longitudinal vortex generator gets thinner causing more change in temperature distribution. The higher temperature gradients that led to a higher heat transfer rate were caused by longitudinal and horseshoe vortices generated by the longitudinal vortex generator and the decrease in thermal boundary-layer thickness. It is also observed that by increasing the Reynolds numbers the temperature difference between the working fluid and hot surfaces gets higher, which causes a higher heat transfer rate. This can be explained that the forced convection becomes more dominant over the diffusive of heat transfer by increasing the Reynolds numbers, and vice versa [5,26]. Additionally, the area of the vortex generator is too low than the overall heat transfer surface and, therefore, the main reason to enhance the heat transfer is due to the increase in the convective heat transfer coefficient by flow distortion (secondary flow).  Figure 7 shows the flow pattern behind the VG at Re = 1100. The temperature distribution for different configurations at different Reynolds numbers with Al 2 O 3 -water as a working fluid is shown in Figure 8. It shows clearly that for all configurations the vortices were generated in the area behind the longitudinal vortex generator. Two different zones of high temperature occurred along the lateral side of the longitudinal vortex generator. The vortices generated by the vortex generator start to cause the two high-temperature zones to merge together and gradually starts to spread in a region close to the heated wall. The longitudinal vortex generator causes an up-wash and a down-wash in the area behind it. The up-wash and a down-wash flows have dual effects: (i) The downwash flow causes the hotter fluid near the heated wall moves at the central region, which leads to a decrease in the thickness of the thermal boundary layer. (ii) The up-wash flow causes the cooler fluid to be transferred from the central region to the heated wall, which causes an increase in the thickness of the thermal boundary layer; that explains why the temperature gradient gets higher near the walls [27].  Figure 7 shows the flow pattern behind the VG at Re = 1100. The temperature distribution for different configurations at different Reynolds numbers with Al2O3-water as a working fluid is shown in Figure 8. It shows clearly that for all configurations the vortices were generated in the area behind the longitudinal vortex generator. Two different zones of high temperature occurred along the lateral side of the longitudinal vortex generator. The vortices generated by the vortex generator start to cause the two high-temperature zones to merge together and gradually starts to spread in a region close to the heated wall. The longitudinal vortex generator causes an up-wash and a down-wash in the area behind it. The up-wash and a down-wash flows have dual effects: (i) The downwash flow causes the hotter fluid near the heated wall moves at the central region, which leads to a decrease in the thickness of the thermal boundary layer. (ii) The up-wash flow causes the cooler fluid to be transferred from the central region to the heated wall, which causes an increase in the thickness of the thermal boundary layer; that explains why the temperature gradient gets higher near the walls [27].   Figure 9 illustrates the velocity distribution for A1 configuration at various Reynolds numbers with Al2O3-water as a nanofluid. It is noticed that a zone of high-velocity forms in the area between the longitudinal vortex generator and near the walls of the microchannel. The recirculation regions are also formed in the area behind the longitudinal vortex generator, which increases with the increase in Reynolds number. Furthermore, it is observed that the mixing of fluid becomes more intense with an increase in Reynolds number, and because of the shear layer generated behind the longitudinal vortex generator, an interaction between the vortices occurs; that causes the vortices to become more intense with the increase in Reynolds number. At higher Reynolds numbers, the flow covers longer distances going through a wavy path downstream of the longitudinal vortex generator; that is one of the reasons behind the enhancement in heat transfer by using a longitudinal vortex generator.
The Nusselt number and Fanning friction factor variations for different configurations for various Reynolds numbers using Al2O3-water as a working fluid are presented in Figure 10. It can be obviously seen that the Nusselt number for different configurations of the microchannels with a longitudinal vortex generator is greater than the plain microchannels in all ranges of the Reynolds  Figure 9 illustrates the velocity distribution for A 1 configuration at various Reynolds numbers with Al 2 O 3 -water as a nanofluid. It is noticed that a zone of high-velocity forms in the area between the longitudinal vortex generator and near the walls of the microchannel. The recirculation regions are also formed in the area behind the longitudinal vortex generator, which increases with the increase in Reynolds number. Furthermore, it is observed that the mixing of fluid becomes more intense with an increase in Reynolds number, and because of the shear layer generated behind the longitudinal vortex generator, an interaction between the vortices occurs; that causes the vortices to become more intense with the increase in Reynolds number. At higher Reynolds numbers, the flow covers longer distances going through a wavy path downstream of the longitudinal vortex generator; that is one of the reasons behind the enhancement in heat transfer by using a longitudinal vortex generator. caused by the presence of secondary flow and the interaction that happened between the vortices and the microchannel walls. It is also found that A1 configuration has the lowest friction factor among all other configurations while other configurations have almost the same friction factor in the considered Reynolds numbers range of this study. In general, the longitudinal vortex generator causes the frictional resistance of the channel walls to increase, due to more pressure drag brought by it [27].  The Nusselt number and Fanning friction factor variations for different configurations for various Reynolds numbers using Al 2 O 3 -water as a working fluid are presented in Figure 10. It can be obviously seen that the Nusselt number for different configurations of the microchannels with a longitudinal vortex generator is greater than the plain microchannels in all ranges of the Reynolds numbers. It is found that the Nusselt number improved by 0.9%-28.1%. It is also noticed that when Reynolds number increases the Nusselt number increases as well. This can be explained as follows: By increasing the Reynolds number, the vortices become more powerful, which augment the fluid mixing and boost the swirl. It is observed that when the Re number reaches approximately 800, a rapid change in Nusselt number for A 1 starts to occur. It can be explained that the longitudinal vortex generator with an acute angle of attack on the first row induces more powerful vortices in comparison to the ones with the obtuse angle of attack. As the Reynolds number increases, the vortices generated in the microchannel become stronger, which causes the decrease in thermal boundary thickness and increase in flow disturbance which improves the swirl in the microchannel. The JF factor gives an evaluation of the overall thermal-hydraulic performance, which indicates the best performance of heat transfer related to the pressure drop. Figure 11 shows the JF values for different configurations for various Reynolds numbers by using CuO-Water as a working fluid. The results indicate that the A1 configuration has the highest JF factor values while other configurations have almost the same the JF factor values in Reynolds numbers range considered for this study.  As shown in Figure 10, microchannels with longitudinal vortex generators have better thermal performance, but it is associated with a pressure drop along the channel. The pressure drop may be caused by the presence of secondary flow and the interaction that happened between the vortices and the microchannel walls. It is also found that A 1 configuration has the lowest friction factor among all other configurations while other configurations have almost the same friction factor in the considered Reynolds numbers range of this study. In general, the longitudinal vortex generator causes the frictional resistance of the channel walls to increase, due to more pressure drag brought by it [27].
The JF factor gives an evaluation of the overall thermal-hydraulic performance, which indicates the best performance of heat transfer related to the pressure drop. Figure 11 shows the JF values for different configurations for various Reynolds numbers by using CuO-Water as a working fluid. The results indicate that the A 1 configuration has the highest JF factor values while other configurations have almost the same the JF factor values in Reynolds numbers range considered for this study.
The JF factor gives an evaluation of the overall thermal-hydraulic performance, which indicates the best performance of heat transfer related to the pressure drop. Figure 11 shows the JF values for different configurations for various Reynolds numbers by using CuO-Water as a working fluid. The results indicate that the A1 configuration has the highest JF factor values while other configurations have almost the same the JF factor values in Reynolds numbers range considered for this study. Figure 11. JF factor versus Reynolds number for various microchannel configurations. Figure 11. JF factor versus Reynolds number for various microchannel configurations.

The Effect of the Different Working Fluids
Based on the best configuration obtained from the previous section, different types of nanoparticles (Al 2 O 3 and CuO) dispersed in different types of base fluid (water, PAO and ethanol glycol) with the nanoparticles volume fraction of 2% and various nanoparticles diameters (29 nm and 28.6 nm) are applied to the A 1 configuration. This microchannel with a longitudinal vortex generator has better thermal performance, but it is associated with a pressure drop along the channel. The aim of this section is to determine the nanofluid with the best thermal-hydraulic performance. Figure 12 demonstrates the effect of changing nanoparticle material on the Nu number and Fanning friction factor for a range of Reynolds numbers and the particle diameter of 29 nm and 28.6 nm. It can be clearly seen that there is a directly proportional relation between Nu number and Re number for all working fluids in the present study. It can be explained that nanofluids have greater heat absorption and the use of metallic nanoparticles increases the thermal conductivity of the mixture, with a larger heat transfer surface area, and thus increases the collision rate of the nanoparticles that increases the nanoparticles' Brownian motion [45]. The increase in thermal conductivity leads to intensive energy transfer, which dominates the diffusion in the working fluid and causes greater bulk temperatures. It is found that CuO-water shows a greater increase in heat transfer compared to Al 2 O 3 -water at the same volume fraction and Reynolds numbers; that can be attributed to the fact that CuO nanoparticles have higher thermal conductivity with the comparison to Al 2 O 3 nanoparticles as well as by the effects of the alignment, structure and interactions of the nanoparticles [49]. conductivity leads to intensive energy transfer, which dominates the diffusion in the working fluid and causes greater bulk temperatures. It is found that CuO-water shows a greater increase in heat transfer compared to Al2O3-water at the same volume fraction and Reynolds numbers; that can be attributed to the fact that CuO nanoparticles have higher thermal conductivity with the comparison to Al2O3 nanoparticles as well as by the effects of the alignment, structure and interactions of the nanoparticles [49]. Based on Figure 12, the enhancement of heat transfer performance comes with the increase in Fanning friction factor; that can be explained as nanofluids have higher viscosity if compared with its base-fluid. Furthermore, at a higher volume fraction, the powerful hydrodynamic interactions between the nanoparticles cause an increase in nanofluids viscosity. Figure 12 depicts at the same Reynolds numbers and volume fraction, that CuO-water nanofluids have a higher Fanning friction Based on Figure 12, the enhancement of heat transfer performance comes with the increase in Fanning friction factor; that can be explained as nanofluids have higher viscosity if compared with its base-fluid. Furthermore, at a higher volume fraction, the powerful hydrodynamic interactions between the nanoparticles cause an increase in nanofluids viscosity. Figure 12 depicts at the same Reynolds numbers and volume fraction, that CuO-water nanofluids have a higher Fanning friction factor in comparison to Al 2 O 3 -water nanofluids. As the Reynolds number increases, the µ Brownian gets lower while µ static gets higher. If it is compared with pure water, this relative increase in nanofluid viscosity gets higher at lower Re number leading to further increase in Fanning friction factor at constant volume fraction; that causes the need for more pumping power to drive nanofluids inside the microchannels [35]. Figure 13 shows the effect of base fluid material on the Nu number and Fanning friction factor for a range of Reynolds numbers and the particle diameter of 29 nm. As shown in Figure 10, in the range of Re number between 200 and 600, it is revealed that the CuO-PAO nanofluid has a higher Nu number among another nanofluid in the scope of the study. But, in the range of Re between 600 and 1200, the Nusselt number for CuO-PAO starts to decrease while CuO-water maintains an increase in Nusselt number. With the increase of Re number the CuO-water takes the lead from CuO-PAO. In all ranges of Re number, CuO-water and EG have the lowest Nusselt numbers among all nanofluids used in this study. Regarding the Fanning friction factor, in the same range of Re number between 200 and 400, almost all base fluids have a higher Fanning friction factor compared to water. In the Re number range of study, CuO-PAO has the lowest Fanning friction factor while CuO-W:EG has the highest Fanning friction factor; that can be explained in that PAO has a low viscosity compared to ethylene glycol, which makes CuO-PAO have a lower Fanning friction factor than CuO-W:EG [40]. It is also shown that CuO-PAO has lower pressure drop than CuO-water, which can be explained that the PAO could have a lubrication effect that reduces the pressure drop along the microchannels. to water. In the Re number range of study, CuO-PAO has the lowest Fanning friction factor while CuO-W:EG has the highest Fanning friction factor; that can be explained in that PAO has a low viscosity compared to ethylene glycol, which makes CuO-PAO have a lower Fanning friction factor than CuO-W:EG [40]. It is also shown that CuO-PAO has lower pressure drop than CuO-water, which can be explained that the PAO could have a lubrication effect that reduces the pressure drop along the microchannels. As shown in Figure 14, the results indicate that Al2O3-water nanofluid has the highest JF factor compared with another nanofluid in the considered Reynolds numbers range; that indicates that Al2O3-water has the best performance in terms of high Nusselt number with a low friction factor among all nanofluids. With regard to the base fluid, the results revealed that CuO-PAO has the highest JF factor in the range of Re number between 200 and 1000, but in a high Reynolds number range between 100 and 1200 the JF factor decreases, which make the CuO-water nanofluid to have the highest JF factor as shown in Figure 15. In almost all the range of Reynolds number in the present As shown in Figure 14, the results indicate that Al 2 O 3 -water nanofluid has the highest JF factor compared with another nanofluid in the considered Reynolds numbers range; that indicates that Al 2 O 3 -water has the best performance in terms of high Nusselt number with a low friction factor among all nanofluids. With regard to the base fluid, the results revealed that CuO-PAO has the highest JF factor in the range of Re number between 200 and 1000, but in a high Reynolds number range between 100 and 1200 the JF factor decreases, which make the CuO-water nanofluid to have the highest JF factor as shown in Figure 15. In almost all the range of Reynolds number in the present study, CuO-water has an increase with respect to the increase in Reynolds number while CuO-W:EG keeps having the lowest JF factor values.   The use of nanofluids as working fluids in microchannels with longitudinal vortex generators cause an observable enhancement in overall thermal performance. Figure 16 represents the summary of thermal performance percentage of the A1 configuration microchannels with various nanofluids for different Reynolds numbers. According to Figure 16, the higher overall thermal performance is observed for the CuO-PAO, CuO-water and Al2O3-water nanofluid, in an ascending order. Furthermore, across the range of Reynolds numbers used in the present study CuO-W:EG   The use of nanofluids as working fluids in microchannels with longitudinal vortex generators cause an observable enhancement in overall thermal performance. Figure 16 represents the summary of thermal performance percentage of the A1 configuration microchannels with various nanofluids for different Reynolds numbers. According to Figure 16, the higher overall thermal performance is observed for the CuO-PAO, CuO-water and Al2O3-water nanofluid, in an ascending order. Furthermore, across the range of Reynolds numbers used in the present study CuO-W:EG The use of nanofluids as working fluids in microchannels with longitudinal vortex generators cause an observable enhancement in overall thermal performance. Figure 16 represents the summary of thermal performance percentage of the A1 configuration microchannels with various nanofluids for different Reynolds numbers. According to Figure 16, the higher overall thermal performance is observed for the CuO-PAO, CuO-water and Al 2 O 3 -water nanofluid, in an ascending order. Furthermore, across the range of Reynolds numbers used in the present study CuO-W:EG showed the lowest overall thermal performance. The CuO-PAO shows the best overall performance because it is due to the highest Nu number and lowest friction factor. It is worthy to mention that the oil in the POA reduced the pressure drop in the channel and the CuO particle increased the thermal properties of the oil, which increased the Nu number. According to previous result, the proposed CuO-PAO nanofluid has a promising solution to solve the problem of a high pressure drop in the microchannel since it also has one of the lowest Fanning frictions due to the lubrication effect of its base fluid. showed the lowest overall thermal performance. The CuO-PAO shows the best overall performance because it is due to the highest Nu number and lowest friction factor. It is worthy to mention that the oil in the POA reduced the pressure drop in the channel and the CuO particle increased the thermal properties of the oil, which increased the Nu number. According to previous result, the proposed CuO-PAO nanofluid has a promising solution to solve the problem of a high pressure drop in the microchannel since it also has one of the lowest Fanning frictions due to the lubrication effect of its base fluid. Figure 16. Thermal performance of the A1 configuration microchannels with various nanofluids for different Reynolds numbers.

Conclusions
This paper has presented a numerical investigation on the thermo-hydraulic performance of a three-dimensional rectangular microchannel with longitudinal vortex generators. The flow was assumed to be steady-state, single-phase, laminar flow. The study consisted of two parts: First, the effects of five different geometrical configurations on the thermo-hydraulic performance of the microchannels with VGs were studied and compared with the smooth microchannels, using the Al2O3-water nanofluid as a working fluid. Second, the effects of the nanofluids with different nanoparticle types (Al2O3 and CuO) and different base fluid (water, PAO and ethanol glycol) with a nanoparticle volume fraction of 2% and various nanoparticles diameters (29 nm and 28.6 nm) applied into the A1 configuration were studied and compared to pure water. The simulations were performed using a CFD code and the results were validated based on the experimental and numerical data leading to the following remarkable conclusions: 1. The assessment of thermal and hydraulic performance on various LVG configurations shows A1 to be the best configuration for LVG arrangement, and after that A2, A4 and A3 are in the list, respectively. The augmentation in Nusselt number was 0.9%-28.1% using the Al2O3-water nanofluid. However, it comes with the penalty of increasing the Fanning friction factor by 5.2%-28% for the Al2O3-water nanofluid with respect to the smooth microchannel.

Conclusions
This paper has presented a numerical investigation on the thermo-hydraulic performance of a three-dimensional rectangular microchannel with longitudinal vortex generators. The flow was assumed to be steady-state, single-phase, laminar flow. The study consisted of two parts: First, the effects of five different geometrical configurations on the thermo-hydraulic performance of the microchannels with VGs were studied and compared with the smooth microchannels, using the Al 2 O 3 -water nanofluid as a working fluid. Second, the effects of the nanofluids with different nanoparticle types (Al 2 O 3 and CuO) and different base fluid (water, PAO and ethanol glycol) with a nanoparticle volume fraction of 2% and various nanoparticles diameters (29 nm and 28.6 nm) applied into the A 1 configuration were studied and compared to pure water. The simulations were performed using a CFD code and the results were validated based on the experimental and numerical data leading to the following remarkable conclusions: 1.
The assessment of thermal and hydraulic performance on various LVG configurations shows A 1 to be the best configuration for LVG arrangement, and after that A 2 , A 4 and A 3 are in the list, respectively. The augmentation in Nusselt number was 0.9%-28.1% using the Al 2 O 3 -water nanofluid. However, it comes with the penalty of increasing the Fanning friction factor by 5.2%-28% for the Al 2 O 3 -water nanofluid with respect to the smooth microchannel.

2.
In case of different base fluids, CuO-PAO has the best performance. The Nusselt number values were 7.67-14.7 and 9.57-15.88, respectively, for Al 2 O 3 -water and CuO-PAO, with the penalty of increasing the Fanning friction factor by 5%-33.6% and 4.2%-26%, respectively, for Al 2 O 3 -water and CuO-PAO.
This device is mainly designed for chip cooling. So, it is recommended to use this technique under high flowrates in order to achieve a higher overall efficiency. Further investigations can be done on various other nano microchannel heat management applications, such as Bio-MEMS in nano drug delivery.