Thermal and Fluid Dynamic Performance Comparison of Three Nanoﬂuids in Microchannels Using Analytical and Computational Models

: The ﬂuid dynamic and thermal performance of three nanoﬂuids containing aluminum oxide, copper oxide, and silicon dioxide nanoparticles dispersed in 60:40 ethylene glycol and water base ﬂuid as a coolant in a microchannel heatsink are compared here by two methods. The ﬁrst is a simple analytical analysis, which is acceptable for very low nanoparticle volumetric concentration (1–2%). The second method is a rigorous three-dimensional ﬁnite volume conjugate heat transfer and ﬂuid dynamic model based upon a constant heat ﬂux boundary condition, which is applicable for cooling electronic chips. The ﬂuids’ thermophysical properties employed in the modeling are based on empirically derived, temperature dependent correlations from the literature. The analytical and computational results for pressure drop and Nusselt number were in good agreement with the nanoﬂuids showing a maximum di ﬀ erence of 4.1% and 2.9%, respectively. Computations cover the practical range of Reynolds number from 20 to 200 in the laminar regime. Based on equal Reynolds number, all of the nanoﬂuids examined generate a higher convective heat transfer coe ﬃ cient in the microchannel than the base ﬂuid, while copper oxide provided the most signiﬁcant increase by 21%. Based on the analyses performed for this study, nanoﬂuids can enhance the cooling performance of the heatsink by requiring a lower pumping power to maintain the same maximum wall temperature. Aluminum oxide and copper oxide nanoﬂuids of 2% concentration reduce the pumping power by 23% and 22%, respectively, while maintaining the same maximum wall temperature as the base ﬂuid.


Introduction
High performance computer chips are cooled by joining them to heatsinks within which microchannels. The high heat flux produced by chips is carried away by the coolant flowing through the channels. The same concept is employed in cold plates, which are used to cool the electronic equipment in aircraft, spacecraft, or the space station, where a large amount of heat is generated in a confined space. Examples of such microchannel cooling have been described in Bergman et al. [1]. Microchannel heat exchangers (MCHX) are compact and increase the heat transfer surface area significantly compared to traditional heat exchangers. Therefore, they hold a great deal of promise in spacecraft. While working on developing innovative heat exchanger to cool chips, Tuckerman and Pease [2] introduced microchannel heatsink. They proposed that the microscale passages can be etched directly within the chip in the non-circuit areas. Microchannel heatsinks are a noteworthy improvement in heat exchanger design over the years to increase thermal performance. While the heat exchanger design has made substantial progress, the heat transfer fluid has remained unchanged.
Depending on the application, the heat transfer liquids widely used are water in warm environments or ethylene glycol and water (EG/W) or propylene glycol and water (PG/W) mixtures at sub-freezing temperatures. National Aeronautics and Space Administration (NASA) spacecrafts operating under subzero condition require glycol as the coolant. Although EG/W has superior thermal properties than PG/W, due to the toxicity of EG/W, NASA prefers the PG/W mixture. Among these three liquids, water has the highest thermal conductivity of 0.60 W/m·K, EG/W 0.36 W/m·K, and PG/W 0.30 W/m·K at room temperature, 20 • C. Comparing these values to aluminum oxide (Al 2 O 3 ) nanomaterials, which has a thermal conductivity of 36.0 W/m·K. By dispersing nanometer size solid particles (≤ 100 nm) into a fluid, the thermal conductivity of the nanofluid is higher than the base fluid. Past research [3][4][5][6] has shown a low volumetric concentration of nanofluids up to 2% can enhance the convective heat transfer coefficient. It has been demonstrated by Ray et al. [7], that the increased thermal performance of nanofluids results in the reduction of the pumping power and surface area (thus size and weight) of a compact heat exchanger, such as an automotive radiator by as much as 30% or 7%, respectively.
In recent years, computational analyses of nanofluid flow in microchannels have appeared in the literature. The book edited by Minkowycz et al. [6] contains extensive research articles on a wide array of topics on nanofluids. Klienstreuer et al. [8] presented results on heat transfer and pressure drop of water based Al 2 O 3 nanofluid in a trapezoidal microchannel using an effective viscosity model proposed by Masoumi et al. [9]. Another recent book edited by Bianco et al. [10] presented extensive results on heat transfer and fluid flow with nanofluids, including some analysis of their performance in microchannels. The treatise by Kandlikar et al. [11] is exclusively devoted to fluid dynamics and heat transfer research in minichannels and microchannels.
As microchannel heat exchangers hold a great deal of promise for the future, the objective of the present study has been to compare the thermal and fluid dynamic performance of aluminum oxide (Al 2 O 3 ), copper oxide (CuO), and silicon dioxide (SiO 2 ) dispersed 60:40 ethylene glycol and water (EG/W) by mass, circulating in a commonly accepted microchannel geometry. Two approaches have been adopted in this paper: (i) a computational fluid dynamic (CFD) modeling using the well-known software ANSYS Fluent [12] Version 18.1 and (ii) a comprehensive analytical modeling following the analytical correlations presented in the authoritative book on microchannel by Kandlikar et al. [11].

Properties of Materials
Four types of materials were used in the present analysis. The base fluid (60:40 EG/W), is the dispersing medium and the nanoparticles, Al 2 O 3 , CuO, and SiO 2 .

Base Fluid Properties
A 60:40 EG/W mixture by mass was adopted. This mixture provides a low freezing point of about −54 • C [13]. The lower temperature is favorable in spacecraft applications or sub-arctic climates. EG/W thermophysical properties were obtained from the American Society of Heating, Refrigerating, and Air-Conditioning Engineers (ASHRAE) Fundamentals Handbook [13]. A correlation was developed for each thermophysical property as a function of temperature, as shown in Table 1. Density, specific heat, and thermal conductivity correlations were modeled as a nondimensional form following Yaws [14] model covering the temperature range of 238 K ≤ T ≤ 398 K. The viscosity correlation follows the log-quadratic empirical fit recommend by White [15] for liquids. The viscosity correlation was improved by splitting the temperature range into two segments 238 K ≤ T ≤ 273 K and 273 K ≤ T ≤ 398 K to achieve an error of less than 0.9%. The correlation in the range of 273 K ≤ T ≤ 398 K is only presented here due to the temperature range within which computations in the present paper have been performed. The full set of correlations are available in Ray et al. [7]. All the thermophysical correlations show a coefficient of determination R 2 ≈ 1. The density and specific heat correlations have an absolute error of less than 0.1%, while viscosity and thermal conductivity are 0.91% and 0.11%, respectively. Compared to other important transport properties of fluids, viscosity and thermal conductivity have a strong dependence on temperature. Therefore, they would require more complex empirical formulas to reduce the curve-fit error. However, that will compromise the simplicity of the present curve fit formulas. For practical engineering design, it is acceptable to have correlations with an error of less than 1%.
The subscript "0" refers to the fluid property at the standard reference temperature of 273.15 K (T 0 ).

Nanofluid Properties
EG/W nanofluids have the most comprehensive thermophysical properties measurement. General correlations were developed by various researchers for density, specific heat, thermal conductivity, and viscosity. For the present study, EG/W nanofluids with Al 2 O 3 , CuO, and SiO 2 will be examined. The thermophysical properties of the nanoparticles are summarized in Table 2.
The majority of the nanofluids studied in this paper were purchased from Alfa Aesar [16]. The thermophysical properties measurements of these EG/W based nanofluids followed the strict nanofluid preparation and characterization procedures outlined in [5,17,18].  (1)) to predict the density of nanofluids. Vajjha et al. [21] conducted experimental measurements of density for three EG/W nanofluids (aluminum oxide, antimony-tin oxide, and zinc oxide). Their measured values were in good agreement with Equation (1).
The "nf ", "p", and "bf " subscripts denote nanofluid, particle and base fluid, respectively. Specific Heat Xuan and Roetzl [22] presented theoretical Equation (2) for determining the specific heat of nanofluids. The present study uses Equation (2) to determine the specific heat of copper oxide and silicon dioxide dispersed in 60:40 ethylene glycol and water mixture.
While for aluminum oxide, it uses the correlation (Equation (3)) that was developed by Vajjha and Das [17]. The curve-fit coefficients A, B, and C for Al 2 O 3 are summarized in Table 3. The percentage deviation in Table 3 represents the difference between the experimental data and curve-fit Equation (2).

Thermal Conductivity
Hamilton and Crosser [23] presented a thermal conductivity model for a solid-liquid mixture. Koo and Kleinstreuer [3] expanded the thermal conductivity model by including Brownian motion for nanofluids, as shown in Equation (4). Following Koo and Kleinstreuer [3] model, Vajjha and Das [5] and Sahoo et al. [24] have a developed correlation (Equation (5)) and the β term (Table 4) for nanoparticles dispersed in 60:40 EG/W mixture. Vajjha and Das [5] and Sahoo et al. [24] utilized steady-state measuring techniques to determine the thermal conductivity of aluminum oxide, copper oxide, zinc oxide, and silicon dioxide nanoparticles dispersed in EG/W. Using Equations (4) and (5), an average deviation for Al 2 O 3 , CuO, and SiO 2 EG/W nanofluids was 0.23%, 5.74%, and 1.97%, respectively. The curve-fit relations for β proposed by Vajjha and Das [5] for Al 2 O 3 , CuO, and Sahoo et al. [24] for SiO 2 are summarized in Table 4. Viscosity Namburu et al. [25,26] and Sahoo et al. [27] measured the viscosities of aluminum oxide, copper oxide, and silicon dioxide dispersed in EG/W with a temperature ranging from 273 to 363 K. Utilizing the data from the previous researchers, Vajjha et al. [20] developed a nondimensional correlation, Equation (6) for those three nanofluids. The curve-ft coefficients A and B for three nanofluids proposed by Vajjha et al. [20] are summarized in Table 5. The maximum deviation between the experimental data and the curve-fit Equation (6) is reported by Vajjha et al. [20] to be ±12%.

Geometries of the Heat Exchanger
The computational fluid dynamic and analytical models were conducted using the geometries of a microchannel presented by Kandlikar et al. [11]. A simple diagram of the microchannel heatsink is shown in Figure 1. The dimensions and material properties are noted in Table 6.
K and specified concentration ranges.

Type of Particle
Namburu et al. [25,26] and Sahoo et al. [27] measured the viscosities of aluminum oxide, copper oxide, and silicon dioxide dispersed in EG/W with a temperature ranging from 273 to 363 K. Utilizing the data from the previous researchers, Vajjha et al. [20] developed a nondimensional correlation, Equation (6) for those three nanofluids. The curve-ft coefficients A and B for three nanofluids proposed by Vajjha et al. [20] are summarized in Table 5. The maximum deviation between the experimental data and the curve-fit Equation (6) is reported by Vajjha et al. [20] to be ±12%.

Geometries of the Heat Exchanger
The computational fluid dynamic and analytical models were conducted using the geometries of a microchannel presented by Kandlikar et al. [11]. A simple diagram of the microchannel heatsink is shown in Figure 1. The dimensions and material properties are noted in Table 6.   Additional surface information is needed before performing an analysis. The equations below show how these surface related characteristics were calculated as presented by Kandlikar et al. [11].
Surface Area: Channel Cross sectional area: Heat Transfer Surface Area: Aspect ratio: Hydraulic diameter: Number of Channels:

Thermal and Fluid Dynamic Calculations
A MATLAB [28] code was developed incorporating the analytical technique outlined by Kandlikar et al. [11] to determine the thermal and fluid dynamic performance of the microchannel heat exchanger. The equations necessary for this analytical scheme are summarized below.
Energy Balance Equation: Using Equation (14), one can determine the exit temperature of the fluid (T 2 ) for a given inlet temperature (T 1 ) and mass or volumetric flow rate, velocity, or Reynolds number.
Han [29] provided entrance length correlations for various aspect ratios. Unfortunately, our given aspect ratio was not available; therefore, interpolating his results was necessary. The hydrodynamic entry length Equation (15) was interpolated for our given aspect ratio of 1/7 from the data presented by Han [29].
Hydrodynamic Entrance Length: Using Equation (15), the predicted hydrodynamic entrance length for the base fluid and nanofluids are the same for equal Reynolds number, if nanofluids are treated as a single-phase fluid. Most researchers thus far have found that dilute concentrations, 1-2 %, of nanofluids can be treated as a single-phase fluid to make computations easier and achieve results of reasonable accuracy. For a higher volumetric concentration of nanoparticles, a two-phase flow analysis will be essential.
Kandlikar et al. [11] presented a correlation for the thermal entrance length for rectangular microchannels from the data presented by Phillip [30].
Thermal Entrance Length:  Figure 2 shows the increase in thermal entry length at equal Reynolds number due to the dependence on the Prandtl number, which is dependent upon the concentration and the material of the nanoparticle.
nanofluids are the same for equal Reynolds number, if nanofluids are treated as a single-phase fluid. Most researchers thus far have found that dilute concentrations, 1-2 %, of nanofluids can be treated as a single-phase fluid to make computations easier and achieve results of reasonable accuracy. For a higher volumetric concentration of nanoparticles, a two-phase flow analysis will be essential.
Kandlikar et al. [11] presented a correlation for the thermal entrance length for rectangular microchannels from the data presented by Phillip [30].
Thermal Entrance Length: Figure 2 shows the increase in thermal entry length at equal Reynolds number due to the dependence on the Prandtl number, which is dependent upon the concentration and the material of the nanoparticle. According to Garimella and Singhal [31], it is reasonable to use fully developed hydrodynamic conditions in heat transfer analysis for microchannels. This is valid for fluids with a low Prandtl Number such as water (Pr = 7.03 @ 293K), but for EG/W (Pr = 44.37 @ 293K) and nanofluids (Al2O3 -2% Pr = 48.30 @ 293K) the analysis would require accounting for developing fluid flow. Kandlikar et al. [11] presented Nusselt Number correlations for four-sided heating in fully developed flow as a function of aspect ratio, as presented by Equation (17). For the microchannel analyzed here with an aspect ratio (αc = 0.143), the fully developed Nusselt Number four-sided heating is 6.23. According to Garimella and Singhal [31], it is reasonable to use fully developed hydrodynamic conditions in heat transfer analysis for microchannels. This is valid for fluids with a low Prandtl Number such as water (Pr = 7.03 @ 293K), but for EG/W (Pr = 44.37 @ 293K) and nanofluids (Al 2 O 3 -2% Pr = 48.30 @ 293K) the analysis would require accounting for developing fluid flow. Kandlikar et al. [11] presented Nusselt Number correlations for four-sided heating in fully developed flow as a function of aspect ratio, as presented by Equation (17). For the microchannel analyzed here with an aspect ratio (α c = 0.143), the fully developed Nusselt Number four-sided heating is 6.23.
Kandlikar et al. [11] presented correlations for the local Nusselt Number in the developing region for four-sided heating. Unfortunately, our chosen aspect ratio was not available and correlation was developed (Equation (18)) by extrapolating the data for aspect ratios of 0.25 and 1/3.
An average Nusselt Number is used to account for the developing flow effects, as shown in Equation (19). x * = x D h RePr (20) x 1 * = 0 (21) After determining the Nusselt number, the average convective heat transfer coefficient can be determined using Equation (24).
Assuming constant and uniform heat transfer coefficient, we can determine the effective heat flux into the heat transfer fluid using Equation (25).
The maximum surface temperature inside the heatsink can be determined using Equation (26) which is an essential thermal performance parameter to examine.
where h 2 represents the heat transfer coefficient at the outlet. The pressure drop in the microchannel heat exchanger can be determined using Equation (27).
Shah and London [32] developed a polynomial correlation, Equation (28), for determining fRe as a function of aspect ratio.
Steinke and Kandlikar [33] developed a fifth-order polynomial curve-fit (Equation (29)) for incremental pressure drop number, K(∞), using the trapezoidal data presented by Kakac et al. [34] taking a rectangular channel as a subset of a trapezoid. The correlation has an accuracy of 0.04% compared to the data [34].
Pumping Power: Friction power per unit surface area [35]: An earlier investigation by Vajjha et al. [18] had shown that the dilute concentration of nanofluids of about 1-2% yielded the best benefit. The higher concentration makes the fluid more viscous, requiring more pumping power, whereas a dilute concentration does not impose a heavy pumping power penalty while giving a significant enhancement in heat transfer. Therefore, concentrations of 1 and 2% are preferred here to compare the thermal and fluid dynamic performances of different fluids.
Using Equation (31), a comparison can be made between the base fluids and three nanofluids for 1 and 2% concentration, as shown in Figure 3. From Figure 3, it is observed the nanofluids perform better than the EG/W base fluid. Aluminum and copper oxide nanofluids of 1% concentration perform better than their 2% concentration. As the concentration increases, the performance gain diminishes.
An earlier investigation by Vajjha et al. [18] had shown that the dilute concentration of nanofluids of about 1-2% yielded the best benefit. The higher concentration makes the fluid more viscous, requiring more pumping power, whereas a dilute concentration does not impose a heavy pumping power penalty while giving a significant enhancement in heat transfer. Therefore, concentrations of 1 and 2% are preferred here to compare the thermal and fluid dynamic performances of different fluids.
Using Equation (31), a comparison can be made between the base fluids and three nanofluids for 1 and 2% concentration, as shown in Figure 3. From Figure 3, it is observed the nanofluids perform better than the EG/W base fluid. Aluminum and copper oxide nanofluids of 1% concentration perform better than their 2% concentration. As the concentration increases, the performance gain diminishes.

Comparison of Performances
Three analyses were conducted:

Comparison of Performances
Three analyses were conducted:

Operational Parameters
A MATLAB code was developed using the analysis scheme presented in Figure 4. The analysis will be conducted using three nanofluids (Al 2 O 3 , CuO, and SiO 2 ) of 1 and 2% concentrations and the base fluid, EG/W. The Reynolds number will be varied from 20 to 200. A uniform heat flux of 10 6 W/m 2 is applied to the bottom of the microchannel heatsink. The fluid will have an inlet temperature of 308.15 K. All the testing parameters are summarized in Table 7. performance. 3. Constant Pumping Power: While maintaining a constant pumping power based on the base fluid, compare the surface temperature as the cooling performance of nanofluids.

Operational Parameters
A MATLAB code was developed using the analysis scheme presented in Figure 4.  The analysis will be conducted using three nanofluids (Al2O3, CuO, and SiO2) of 1 and 2% concentrations and the base fluid, EG/W. The Reynolds number will be varied from 20 to 200. A uniform heat flux of 10 6 W/m 2 is applied to the bottom of the microchannel heatsink. The fluid will have an inlet temperature of 308.15 K. All the testing parameters are summarized in Table 7. q" 10 6 Initial temperature (K) T0 308.15

Computational Fluid Dynamic Analysis
A three-dimensional conjugate heat transfer model using Fluent Release 18.1 [12] was developed to validate the analytical scheme by comparing key parameters such as pressure drop and Nusselt Number. With the assumption of uniform heat flux, fluid distribution and neglecting the effects of the outer channels, one can analyze the system by examining half of a single channel, as shown in

Computational Fluid Dynamic Analysis
A three-dimensional conjugate heat transfer model using Fluent Release 18.1 [12] was developed to validate the analytical scheme by comparing key parameters such as pressure drop and Nusselt Number. With the assumption of uniform heat flux, fluid distribution and neglecting the effects of the outer channels, one can analyze the system by examining half of a single channel, as shown in Figure 5, by applying appropriate boundary conditions of symmetry. The cross-sectional coordinates are x and y, while the axial coordinate is z. The origin is at the bottom left corner.  The testing parameters for the numerical model are inlet velocity and fluid. The velocity is varied from 0.25 to 3.5 m/s for all three nanofluids (Al 2 O 3 , CuO, and SiO 2 ) at 1 and 2% concentration as well as the base fluid, EG/W. For this model, nanofluids are assumed to be a homogeneous engineered fluid with enhanced thermophysical properties. Past research has shown that the single-phase computational model is quite accurate for dilute nanoparticle concentrations of 1 to 2% but is acceptable up to 4%. Using a single-phase fluid model, the nanofluids are assumed to have thermophysical properties determined by correlations described in Section 2.1.2.

Governing Equations
The governing Equations (32)-(34) are for steady-state, incompressible, and laminar fluid flow in the microchannel while neglecting the gravitational force, natural convection, and the heat dissipation caused by viscosity. Equation (35) accounts for the heat conduction through the solid walls.

Numerical Methods
The system of governing Equations (32)- (35) were solved by the control volume approach using ANSYS Fluent [12]. For the solver setting, the standard pressure-based and steady-state were chosen. Using the laminar and energy model, the solution method utilized the Semi-Implicit Method for Pressure Linked Equations (SIMPLE) scheme with spatial discretization for gradient, pressure, and momentum being least squares cell based, second order, and second order upwind, respectively.
The models were initialized using the hybrid method, after which the residuals for continuity and velocities were closely monitored. Convergence was achieved for the simulation when the residuals were less than 10 −6 .

Boundary Conditions
The model has a uniform axial velocity and temperature applied at the inlet (z = 0 µm). At the base of the microchannel heatsink (y = 0 µm), a uniform heat flux is applied. The no-slip condition (V = 0 m/s) is applied to the channel wall surfaces. Symmetry boundary conditions were applied at y = 0 µm and y = s c /2 + a c .

Mesh Independence Study
A mesh independence study was conducted to achieve the least number of elements that can yield accurate computational results. Six different mesh sizes were chosen for discretizing the computational domain of the microchannel heatsink. They are Mesh-I, Mesh-II, Mesh-III, Mesh-IV, Mesh-V, and Mesh-VI. In each case, a converged solution of outlet velocity, outlet temperature, pressure drop, and average wall temperature was obtained. In Table 8, the difference between these parameters for each and the previous mesh are summarized. Mesh-I is the baseline with the least number of elements. For all six meshes, velocity, fluid temperature, and average wall temperature were in good agreement with each other. Pressure drop was the only parameter that slightly varied between each mesh with a maximum difference of 0.44%, which is quite low. Therefore, mesh independence was achieved for Mesh-I and the denser meshes. However, when the shear stress and the heat transfer coefficient in the rapidly diminishing entrance region were examined, they were not captured well by Mesh-I through Mesh-III due to the smaller number of elements, in comparison to the results from the refined Mesh-IV through Mesh-VI. From this mesh independence study, Mesh-IV was selected because of low computation time and used in all subsequent computational runs.

Analytical Results and Computational Model Validation
The computational results were validated for all fluids by comparing the pressure drop and Nusselt Number. For brevity, here we presented the comparison of parameters for EG/W and summarized the results for the nanofluids in Table 9. Figure 6 shows the pressure drop in the MCHX for EG/W. The numerical results from Fluent and the results from Equation (27) are in good agreement within a maximum difference of 4.4%. Table 9. Nanofluids results from a comparison between analytical and CFD schemes.  In Figure 7, a comparison between the numerical results and developing and fully developed Nusselt number for EG/W is shown. The developing flow average Nusselt Number for four-sided heating (Nuavg4, Equation (19)) matches well with the numerical results with a maximum difference of 3%.  Figure 7, a comparison between the numerical results and developing and fully developed Nusselt number for EG/W is shown. The developing flow average Nusselt Number for four-sided heating (Nu avg4 , Equation (19)) matches well with the numerical results with a maximum difference of 3%.  Figure 7, a comparison between the numerical results and developing and fully developed Nusselt number for EG/W is shown. The developing flow average Nusselt Number for four-sided heating (Nuavg4, Equation (19)) matches well with the numerical results with a maximum difference of 3%. The nanofluids had a similar agreement with the correlations with a maximum difference of 4.1% and 2.9% for pressure drop and Nusselt number, respectively, as shown in Table 9. This validates the ability of the analytical and computational scheme to predict the correct results for the thermal and fluid dynamic performance of the microchannel heatsink. The nanofluids had a similar agreement with the correlations with a maximum difference of 4.1% and 2.9% for pressure drop and Nusselt number, respectively, as shown in Table 9. This validates the ability of the analytical and computational scheme to predict the correct results for the thermal and fluid dynamic performance of the microchannel heatsink. Figure 8 shows that higher particle concentration and Reynolds number increases the convective heat transfer coefficient. The convective heat transfer coefficient was increased by 18% and 21% for Al 2 O 3 -2% and CuO -2% nanofluids, respectively, over the base fluid EG/W. The SiO 2 -2% nanofluid due to relatively lower thermal conductivity increased the heat transfer coefficient by 8% only.   Figure 8 shows that higher particle concentration and Reynolds number increases the convective heat transfer coefficient. The convective heat transfer coefficient was increased by 18% and 21% for Al2O3 -2% and CuO -2% nanofluids, respectively, over the base fluid EG/W. The SiO2 -2% nanofluid due to relatively lower thermal conductivity increased the heat transfer coefficient by 8% only. From Figure 9, we can see how nanofluids can reduce the maximum surface temperature of the heatsink when compared to the base fluid. Nanofluids show the greatest improvement at the lower Reynolds number with CuO -2% dropping the maximum surface temperature approximately by 13.5K. Al2O3 -2% also showed similar performance by lowering the maximum surface temperature by 10.4K. From Figure 9, we can see how nanofluids can reduce the maximum surface temperature of the heatsink when compared to the base fluid. Nanofluids show the greatest improvement at the lower Reynolds number with CuO -2% dropping the maximum surface temperature approximately by 13.5 K. Al 2 O 3 -2% also showed similar performance by lowering the maximum surface temperature by 10.4 K.

Constant Reynolds Number
Processes 2020, 8,754 15 of 20 Figure 9. Effects of Reynolds number on maximum surface temperature.
From Figures 8 and 9, the nanofluids showed great thermal performance when compared to their base fluid (EG/W). This increase in thermal performance will come at the cost of increased pumping power, as shown in Figure 10. CuO -2% increased the pumping power requirements by as much as 180% (+1-16 W) for Re > 55 and for a Reynolds number of 20 the increase was 300% (+0.08 W). From Figures 8 and 9, the nanofluids showed great thermal performance when compared to their base fluid (EG/W). This increase in thermal performance will come at the cost of increased pumping power, as shown in Figure 10. CuO -2% increased the pumping power requirements by as much as 180% (+1-16 W) for Re > 55 and for a Reynolds number of 20 the increase was 300% (+0.08 W). From Figures 8 and 9, the nanofluids showed great thermal performance when compared to their base fluid (EG/W). This increase in thermal performance will come at the cost of increased pumping power, as shown in Figure 10. CuO -2% increased the pumping power requirements by as much as 180% (+1-16 W) for Re > 55 and for a Reynolds number of 20 the increase was 300% (+0.08 W).

Constant Maximum Surface Temperature
While maintaining a constant maximum surface temperature based on the base fluid calculations, determine the nanofluid fluid dynamic performance. With this analysis, we can determine how much pumping power is required for nanofluids while maintaining the same conditions of the base fluid.
In Figure 11, we can see that the nanofluids can meet the same constraints as the base fluid but with less pumping power requirements. The figure shows that nanofluids with concentrations of 1% and 2% can perform better (reduce pumping power) than the base fluid. The best performing nanofluids are Al 2 O 3 -1% and CuO -1% that reduced the pumping power by approximately 23.0 and 22.2%, respectively, while maintaining the maximum surface temperature at 315 K. This occurred with a Reynolds number of 153 and 144 for Al 2 O 3 and CuO, respectively. Here nanofluids exhibited better performance at the higher Reynolds number (≥100) versus the lower Reynolds number (≤ 100). In further analysis, nanofluids exhibited a continued decrease in pumping power to a Reynolds number of 300, beyond which the pumping power gradually increases slightly.

Constant Maximum Surface Temperature
While maintaining a constant maximum surface temperature based on the base fluid calculations, determine the nanofluid fluid dynamic performance. With this analysis, we can determine how much pumping power is required for nanofluids while maintaining the same conditions of the base fluid.
In Figure 11, we can see that the nanofluids can meet the same constraints as the base fluid but with less pumping power requirements. The figure shows that nanofluids with concentrations of 1% and 2% can perform better (reduce pumping power) than the base fluid. The best performing nanofluids are Al2O3 -1% and CuO -1% that reduced the pumping power by approximately 23.0 and 22.2%, respectively, while maintaining the maximum surface temperature at 315 K. This occurred with a Reynolds number of 153 and 144 for Al2O3 and CuO, respectively. Here nanofluids exhibited better performance at the higher Reynolds number (≥ 100) versus the lower Reynolds number (≤ 100). In further analysis, nanofluids exhibited a continued decrease in pumping power to a Reynolds number of 300, beyond which the pumping power gradually increases slightly. With a reduction in pumping power, the nanofluids would also reduce entropy generation, as shown in Figure 12. Using Equation (36), entropy generation per unit length [36] can be calculated. With a reduction in pumping power, the nanofluids would also reduce entropy generation, as shown in Figure 12. Using Equation (36), entropy generation per unit length [36] can be calculated.
The minimization of entropy generation is essential to enhance the thermal performance of heat exchangers. The lower is the entropy generation; the better is the heat exchanger. Entropy generation is proportional to the irreversibility in a thermodynamic process. By reducing entropy generation, we lower irreversibility and produce a better thermodynamic efficiency for the heat exchanger based on the second law of thermodynamics. Figure 12 shows 1% concentration nanofluids generated a less entropy than compared to the base fluid. Particularly, Al 2 O 3 -1% and CuO -1% nanofluids decreased entropy generation by 24% than the base fluid, while performing the same cooling. Thus, some nanofluids of low concentration can yield superior performance in heat exchangers.
. Figure 12. Entropy generation reduction while maintaining thermal performance.
The minimization of entropy generation is essential to enhance the thermal performance of heat exchangers. The lower is the entropy generation; the better is the heat exchanger. Entropy generation is proportional to the irreversibility in a thermodynamic process. By reducing entropy generation, we lower irreversibility and produce a better thermodynamic efficiency for the heat exchanger based on the second law of thermodynamics. Figure 12 shows 1% concentration nanofluids generated a less entropy than compared to the base fluid. Particularly, Al2O3 -1% and CuO -1% nanofluids decreased entropy generation by 24% than the base fluid, while performing the same cooling. Thus, some nanofluids of low concentration can yield superior performance in heat exchangers.

Constant Pumping Power
In this analysis, the pumping power is held constant for all fluids. With this constraint, we can determine various thermal performance, such as maximum and entropy generation for the microchannel heatsink with nanofluids.
In Figure 13, it is observed that the nanofluids cannot reduce the maximum surface temperature while maintaining the equal pumping power requirements of the base fluid.

Constant Pumping Power
In this analysis, the pumping power is held constant for all fluids. With this constraint, we can determine various thermal performance, such as maximum and entropy generation for the microchannel heatsink with nanofluids.
In Figure 13, it is observed that the nanofluids cannot reduce the maximum surface temperature while maintaining the equal pumping power requirements of the base fluid.

Conclusions
From analytical and CFD modeling of nanofluids flow in a microchannel based on a detailed examination of the results, the following conclusions are drawn.
Considering the condition of equal friction power per unit surface area, nanofluids generated higher heat transfer coefficients than the base fluid, as shown in Figure 3.
The pressure drop and Nusselt number results for a range of Reynolds numbers agree quite well within 4.5% between the analytical model of Kandlikar et al. and our ANSYS model for dilute concentration of nanofluids.
From Figure 8, the CuO -2% exhibits the biggest improvement in the heat transfer coefficient

Conclusions
From analytical and CFD modeling of nanofluids flow in a microchannel based on a detailed examination of the results, the following conclusions are drawn.
Considering the condition of equal friction power per unit surface area, nanofluids generated higher heat transfer coefficients than the base fluid, as shown in Figure 3.
The pressure drop and Nusselt number results for a range of Reynolds numbers agree quite well within 4.5% between the analytical model of Kandlikar et al. and our ANSYS model for dilute concentration of nanofluids.
From Figure 8, the CuO -2% exhibits the biggest improvement in the heat transfer coefficient compared to the EG/W base fluid over the range 20 ≤ Re ≤ 200. The inside heat transfer coefficient for the CuO -2% 60/40 EG/W nanofluid exceeds that of the base fluid by 21% at Re = 200. The highest value is offered by CuO, followed by Al 2 O 3 and then SiO 2 nanofluid. An increase in nanofluid concentration increases the heat transfer coefficient.
From Figure 9, observe that considering equal inlet Reynolds number, all the nanofluids examined lowered the maximum temperature of the wall temperature of the heatsink. The wall temperature can be lowered from 349 K using the base fluid to 335 K at a Reynolds number of 20. In electronic cooling applications where microchannel heat exchangers are used, nanofluids hold promise for lowering operating temperatures of critical components and contribute to increased life and system reliability.
The higher heat transfer coefficient comes at a pumping power penalty. From Figure 10, one observes that the CuO -2% nanofluid requires the highest pumping power among all nanofluids analyzed here considering equal Reynolds number.
While maintaining equivalent thermal performance to the base fluid, nanofluids reduce the pump power requirements. From Figure 11, Al 2 O 3 -1% and CuO -1% reduced the pumping power by approximately 23.0 and 22.2%, respectively, while maintaining the maximum heatsink surface temperature at 315 K.
Entropy generation is proportional to the irreversibility in a thermodynamic process. By reducing entropy generation, the irreversibility is lowered and produces a better thermodynamic efficiency for the heat exchanger based on the second law of thermodynamics. Two of the nanofluids examined generated lower entropy compared to the base fluid, as shown in Figure 12. Al 2 O 3 -1% and CuO -1% nanofluids both decreased entropy generation by 24% when compared to the base fluid with maintaining the same maximum heatsink surface temperature as the base fluid.
The improvement in heat transfer performance in microchannels with nanofluids comes at the cost of higher pumping power. While maintaining an equal thermal performance to the base fluid, nanofluids can reduce the pumping power requirements and contribute to a more thermodynamically efficient heat exchange process.