Next Article in Journal
Faraday Waves-Based Integrated Ultrasonic Micro-Droplet Generator and Applications
Next Article in Special Issue
1600 Parallel Microchamber Microfluidic Device for Fast Sample Array Preparation Using the Immiscibility of Two Liquids
Previous Article in Journal
Ultrasonic Spray-Coating of Large-Scale TiO2 Compact Layer for Efficient Flexible Perovskite Solar Cells
Previous Article in Special Issue
A Microfluidic Chip for Liquid Metal Droplet Generation and Sorting
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Performance Analysis of Double-Layer Microchannel Heat Sinks under Non-Uniform Heating Conditions with Random Hotspots

Department of Mechanical Engineering, Inha University 100, Inha-Ro, Nam-Gu, Incheon 22212, Korea
*
Author to whom correspondence should be addressed.
Micromachines 2017, 8(2), 54; https://doi.org/10.3390/mi8020054
Submission received: 24 December 2016 / Revised: 20 January 2017 / Accepted: 7 February 2017 / Published: 14 February 2017

Abstract

:
Performance analysis of double-layer microchannel heat sinks was performed under non-uniform heating conditions having randomly distributed hotspots. Two parallel-channel (parallel-flow and counter-flow) and one cross-channel (transverse-flow) designs of double-layer heat sink were evaluated with three sets of heating schemes. Each set of heating scheme consisted of eleven randomly distributed hotspots generated by Latin hypercube sampling. The heat flux, area, and location of the hotspots were selected as the design parameters. Conjugate heat transfer analysis of the heat sinks was performed by solving three-dimensional Navier–Stokes and energy equations. Water with temperature-dependent properties was selected as the coolant. The thermal resistance, pressure drop, maximum temperature rise, and temperature variation among hotspots were evaluated for all the heat sinks. The transverse-flow microchannel heat sink exhibited the lowest thermal resistance, temperature rise and temperature variation among the hotspots throughout the specified range of flow rate. The lowest pressure drop was exhibited by the counter-flow heat sink.

Graphical Abstract

1. Introduction

A microprocessor is composed of various blocks producing different amounts of heat, and the logic-block being the highest contributor creates a localized high heat flux region known as “hotspot” [1]. Design complexity and dynamic power dissipation in a single core processor caused the design architecture to shift towards multicore [2]. In a multicore processor, most of the heat is concentrated at the cores creating multiple hotspots. High temperatures at the hotspots create large temperature gradients which induce thermal stresses creating circuit imbalances that could affect the efficiency and service life of a microprocessor [3]. With the advancement in manufacturing techniques [4], a microprocessor with 5.4 billion transistors was recently developed at IBM [5]. The thermal design powers (TDPs) of central processing units (CPUs) and graphics processing units (GPUs) have reached 150 W and 250 W, respectively [6]. Air-cooling techniques have reached their limits [7], and thus there is an eminent need to explore more effective cooling solutions that can keep up with the heat dissipation requirement of the advanced microprocessors.
Tuckerman and Pease [8] showed promising results by using a microchannel heat sink with water as a coolant. Following their work, numerous studies have been presented to explore and enhance the thermo-hydraulic performance of microchannel heat sinks. Liu and Garimella [9] performed numerical and experimental study on fluid flow in microchannels with wide ranges of hydraulic diameters and Reynolds numbers. They concluded that the conventional Navier–Stokes theory can be employed using commercial software packages to understand the flow in microchannels. Several active [10,11] and passive [11,12,13,14] heat transfer augmentation techniques have been explored to enhance the heat transfer capacity. Passive techniques (as considered in the present work) are the preferred choice over active techniques as they have simple design and also do not require external power therefore eliminating additional complexity and costs [11].
Double-layer microchannel heat sink gained attention for its higher heat transfer capability over single-layer microchannel heat sink [15,16]. A better temperature uniformity is offered by the double-layer counter-flow arrangement as compared to the double-layer parallel-flow arrangement [17]. Performance enhancements of double-layer microchannel heat sinks were reported by using converging channels [18,19] and passive microstructures [12,13,14,20], and by optimizing channel shapes [21,22,23]. Along with double-layer parallel-flow and counter-flow heat sinks, Ansari and Kim [24] analyzed the performance of a double-layer transverse-flow heat sink design with various flow configurations. They reported that the transverse-flow heat sink showed better temperature uniformity and lower thermal resistance as compared to double-layer parallel-flow and counter-flow heat sinks. Use of temperature-dependent fluid properties for numerical analysis is recommended by many researchers unless the temperature difference is very small as the constant fluid property method is crude approximation of the real phenomenon [25,26].
In a single core processor, the maximum clock frequency has saturated around 4 GHz due to increased design complexity [2]. Since 2004, the processor architecture has shifted towards multicore to keep up with the increasing requirement for computational capability [2,27]. The high heat flux concentrated at the cores as compared to the remaining chip area [28,29,30,31], creates large temperature gradient across the chip which leads to circuit imbalances and eventually degradation in the microchip performance [32]. Using spatially-resolved imaging of microprocessor power (SIMP) technique, Hamann et al. [33] analyzed the power distribution of a fully working dual core microprocessor. They presented that the flux intensity as well as the actual position of a hotspot in a dual core microprocessor changes quite drastically depending upon the level of utilization.
Hotspots in a microchip are currently being identified as a prevalent challenge in the design of an efficient heat sink. If heat sinks are designed considering the maximum heat flux at the hotspots, it would lead to over-cooling of the remaining background area of the microchip, while designing it considering average heat flux generated by microchip would lead to under-cooling of high heat flux regions. Abdoli et al. [34] numerically investigated performance of a micro pin-fin array heat sink with a single hotspot at the center and presented the effects of the pin-fin shape and the height on the thermo-hydraulic performance using six different pin-fin shapes. Effect of manifold design on the performance of a microchannel heat sink with non-uniform heating was presented by Liu et al. [35]. Various cooling techniques have been proposed for hotspot management using minichannel heat sink and thermoelectric cooler [36], hierarchical manifold microchannel heat sink [37], and micro heat sink with non-uniform pin-fin distribution [38]. Feng et al. [39] performed an analysis of a pin-fin heat sink under non-uniform impingement heating, and presented an analogy model for the fast calculation of thermal resistance with reasonable accuracy as compared to numerical calculation.
With maximum clock frequency of a single core processor at its limits [2], almost all future CPUs and GPUs would have multiple cores, creating multiple high heat flux zones. In addition, the intensity and the location of high heat flux varies with the processors utilization [33]. Thermal analysis of heat sinks considering random non-uniform heating schemes, would closely imitate the actual heating conditions of future multicore processors, and the results would be more reliable and practical. In the present study, performance of a previously proposed flow configuration of double-layer transverse-flow microchannel heat sink [24], was evaluated under random non-uniform heating condition and compared with that of double-layer parallel-flow and counter-flow heat sinks. Three different heating schemes with eleven random hotspots in each scheme were selected as heating conditions. Heat flux, area, and location were selected as the hotspots parameters. Latin hypercube sampling (LHS) [40,41] method was used to generate the hotspots in each set. Full heat sink was selected as the computational domain and a three-dimensional (3-D) conjugate heat transfer analysis was performed using Navier–Stokes and energy equations with temperature-dependent fluid properties. The performances of the heat sinks were evaluated in terms of maximum temperature rise at the hotspots, temperature variation among the hotspots, pressure drop and the total thermal resistance.

2. Heat Sink Models and Heating Conditions

The schematic configurations of two double-layer microchannel heat sinks (parallel-channel and cross-channel) analyzed under random non-uniform heating conditions are shown in Figure 1a,b. Each heat sink model was composed of total forty microchannels with twenty microchannels in each layer. In case of the parallel-channel design (Figure 1a), the channels in both the top and the bottom layers are aligned parallel to each other, whereas, in case of the cross-channel design (Figure 1b), the channels in the top layer are transversal (at right angles) to those in the bottom layer. To evaluate the performances of these two heat sink designs under non-uniform heating conditions, three flow configurations were selected. The two flow configurations from the parallel-channel design (PF and CF) and one flow configuration from the cross-channel design (TF) are shown in Figure 1d. In the PF configuration, the flow directions in both the bottom and the top layers are same. In the CF configuration, the flow direction in the top layer is opposite to that in the bottom layer. In the TF configuration, these two flow directions are transversal (at right angles); moreover, each layer is subdivided into quarters with opposite flow directions in the adjacent quarters. The TF flow configuration considered here showed the best performance among all the analyzed flow configurations of the parallel-channel and cross-channel designs in the previous study presented by the authors [24].
Both the heat sinks were designed with exactly same values of the geometric parameters as listed in Table 1. The analyzed heat sink designs can be manufactured using the process presented by Wei et al. [17]. Improper selection of inlet manifold could lead to non-uniformity in flow distribution resulting in high temperature rise [35] and increased pressure drop at high flow rates [35,42]. To uniformly distribute the coolant in microchannels, a multi-layer manifold design can be employed as presented by Wei et al. [17]. In this work, however, the inlet manifold design was not included in the computational domain due to the limitation of computational resources. Thus, to minimize the computational cost, the present research was concentrated to understand the effect of randomly distributed hotspots on the performance of the double-layer microchannel heat sinks.
Three hotspot schemes (HS1, HS2, and HS3) were considered with eleven hotspots in each scheme, where the hotspots were generated randomly by using a design-of-experiment technique, LHS [41]. The heat flux (qhs), area (Ahs), and location (x and z) of the hotspots were selected as the hotspot parameters. Table 2 shows the hotspot design parameters with their ranges. The hotspot design schemes, HS1, HS2, and HS3, are shown in Figure 2, and their hotspots are listed in Table 3, Table 4 and Table 5, respectively. In experimental case, the temperature measurements at the hotspots can be performed by utilizing infrared (IR) imaging technique using an IR-transparent heat sink as presented by Hamann et al. [33].

3. Numerical Analysis

A general purpose commercial CFD code, ANSYS CFX 15.0 [43] was used to perform three-dimensional conjugate heat transfer analyses in microchannel heat sinks. The code uses an implicitly coupled algebraic multigrid solver [44], which solves the governing equations for mass, momentum, and energy conservation. In the multigrid technique, the solver performs initial iterations on a fine mesh, then, on a coarser mesh which is created virtually, and finally the results are transferred back to the finer mesh to obtain accurate results. The coupled solver reduces the convergence time by solving an entire set of the equations simultaneously. The calculation of each equation was considered to have converged when the root-mean-square residual was less than 10−6.
Several assumptions were made to simplify the numerical calculation. Flow inside the microchannels was considered to be steady, incompressible, and laminar. The gravitational effect and radiation heat transfer were neglected. The steady-state governing equations for conjugate heat transfer in fluid and solid regions are describe as follows:
Continuity equation:
· ( ρ f V ¯ ) = 0
Momentum equations:
V ¯ · ( ρ f V ¯ ) = p + · ( μ f V ¯ )
Energy equation (liquid domain):
V ¯ · ( ρ f C p T f ) = · ( k f T f )
Energy equation (solid domain):
· ( k s T s ) = 0
To evaluate the performance of the heat sinks, total thermal resistance, maximum temperature rise at hotspots, temperature variation among hotspots, and total pressure drop, were calculated. The total thermal resistance of heat sinks under uniform as well as non-uniform heating conditions was defined as follows [8,39]:
R th , tot = T s , max   T f , in Q tot
where Ts,max is the maximum temperature at the bottom surface (solid) of heat sink, Tf,in is the inlet fluid temperature, and Qtot is the total heat applied at the bottom surface (solid) of the heat sink.
The Qtot for heat sink without hotspots was defined as:
Q tot , WHS = q b A b
where qb and Ab are the heat flux applied at the bottom surface and the total area of the bottom surface of the heat sink, respectively. For heat sink with hotspots, it was defined as:
Q tot , HS = q bg A bg + ( q hs A hs )
where qbg is the heat flux applied at the background area of the heat sink (area excluding the hotspots) and Abg is the background area of the heat sink. qhs and Ahs are the heat flux applied at each hotspot and the area of each hotspot, respectively. ( q hs A hs ) is the total heat applied at all (eleven) the hotspots.
The total pumping power required to pump the coolant through the channels was calculated as:
P tot = n ch u avg A c Δ p ch
where n ch is the total number of channels in both layers, u avg is the average inlet velocity of fluid, A c is the channels cross-sectional area, and Δ p ch is the pressure drop within a single channel.
The Reynolds number was defined as:
R e = ρ f u avg D h μ f
where ρf is the fluid (water) density, Dh is the channel hydraulic diameter, and μf is the dynamic viscosity of the fluid.
In the case of cross-channel design (TF configuration) (Figure 1b), selection of a periodic or symmetric computational domain is not possible, while it is possible for the parallel-channel design (PF and CF configurations) (Figure 1a). However, to maintain the consistency in result presentation, the complete heat sinks were selected as the computational domain for all the flow configurations (PF, CF, and TF). As boundary conditions, a uniform mass flow rate was assigned at all the channel inlets (thus, the flow rates in the bottom and top layers are same), and atmospheric pressure condition was assigned at the outlets. A uniform heat flux (100 W/cm2) was applied at the bottom surface of heat sinks without hotspots (qb) and the same was used as the background heat flux (area excluding hotspots) of heat sinks with hotspots (qbg).
Water represents superior thermophysical properties [16], and thus it was selected as the coolant. The inlet temperature of water ( T f , in ) and the heat sink (silicon) was set at 300 K. Silicon being the most widely used semiconductor material was selected as the heat sink material. The thermophysical properties of silicon, i.e., the density (ρs), specific heat (Cp,s), and the thermal conductivity (ks), were assigned as 2330 kg/m3, 712 J/kg·K, and 148 W/m·K, respectively [45].
The thermophysical properties of water were allowed to vary with temperature. Polynomial correlations, between the thermophysical properties of water and temperature in the property table [45], were developed and verified by using MATLAB [40] in this work. The correlations for density (ρw), specific heat (Cp,w), thermal conductivity(kw), and dynamic viscosity (μw) are described as follows:
ρ w ( T ) = 999.9 + 9.561 × 10 2 T 1.013 × 10 2 T 2 + 8.459 × 10 5 T 3 3.496   10 7 T 4
C p , w ( T ) = 4217 3.452 T + 1.155 × 10 1 T 2 1.862 × 10 3 T 3 + 1.538 × 10 5 T 4 4.850 × 10 8 T 5
k w ( T ) = 5.698 × 10 1 + 1.772 × 10 3 T 4.870 × 10 6 T 2 2.915 × 10 8 T 3 + 1.094 × 10 10 T 4
μ w ( T ) = 1.750 × 10 3 5.558 × 10 5 T + 1.172 × 10 6 T 2 1.579 × 10 8 T 3 + 1.169 × 10 10 T 4 3.535   × 10 13 T 5
where T is the temperature in degrees Celsius (°C).
The prediction accuracies of the above correlations were verified by comparing the predicted values obtained by the correlations with the values obtained from the property table [45], as shown in Figure 3. Excellent agreements between the predicted values from the correlations and the original data [45] can be observed in Figure 3. The maximum percentage residuals for the density, specific heat, thermal conductivity, and dynamic viscosity were 0.04%, 0.01%, 0.16%, and 0.71%, respectively. The goodness-of-fit parameters [40] to assess the prediction quality of the correlations, such as R-square ( R 2 ), adjusted R-square ( R adj 2 ), sum of squares due to error (SSE), and the root mean squared error (RMSE), are listed in Table 6. Values of R 2 and R adj 2 closer to unity indicate that the fit is accurate, whereas for SSE and RMSE, a value closer to zero implies a better fit.

4. Results and Discussion

To create finite volumes in the computational domains for numerical calculations, unstructured tetrahedral and hexahedral elements were generated in the solid and fluid computational domains, respectively (with fine meshes near the walls of the fluid domain to resolve the boundary layer flow). A comprehensive grid sensitivity analysis was performed for consistent grid-independent results. Double-layer PF configuration with a single hotspot was selected for grid analysis at total flow rate of 0.006 kg/s. A wide range of grid nodes were tested; the number of nodes was varied from 3.9 million to 19.2 million. A grid with about 11 million nodes was determined to be optimum grid as no noticeable change in the results was observed with further increase in the grid density as shown in Figure 4 (13.1 million nodes showed 0.007% and 0.154% changes in the maximum temperature rise and the pressure drop, respectively, as compared to 11 million nodes and a grid with 8.7 million nodes showed 0.04% and 0.4% changes in the maximum temperature rise and the pressure drop, respectively, as compared to 11 million nodes).
The numerical results using temperature-dependent fluid properties were validated with the previous numerical and experimental results [17,46] for the temperature distribution and the pressure drop, as shown in Figure 5. The present numerical results show good agreements with the numerical and experimental results presented by Wei et al. [17] for the temperature variation at the base of a double-layer counter-flow heat sink with rectangular microchannels as shown in Figure 5a. Qu and Mudawar [46] presented experimental results for the pressure drop within a rectangular microchannel. Figure 5b shows good agreements between the present numerical and previous experimental data for the pressure drop in a range of Reynolds number.
In the present study, PF, CF, and TF heat sink designs were selected for the performance analysis under three randomly generated non-uniform heating schemes. The cross-channel heat sink design with the TF configuration considered in the present work showed the lowest thermal resistance as compared to various other cross-channel and parallel-channel designs in the previous work [24]. The maximum increase in temperature at the hotspots, the pressure drop variation, the thermal resistance and the temperature variation among all the hotspots were selected as the performance parameters to be evaluated. A total flow rate range of 0.006 to 0.014 kg/s (Re range from 400 to 933) was used to analyze the performance variation with flow rate. The minimum value of flow rate was selected so that the maximum temperature rise at the hotspots became below the recommended temperature range for reliable performance of a microprocessor [7].
The local temperature contours at y/Ly = 0 (bottom) and y/Ly = 0.5 (middle) for the three heat sinks without hotspots (PF-WHS, CF-WHS, and TF-WHS) at total flow rate of 0.01 kg/s are shown in Figure 6. In the case of PF heat sink, the hottest zone is observed near the outlet region for both the bottom (y/Ly = 0) and the top layers (y/Ly = 0.5) (Figure 6a). Due to the opposite flow arrangement in the top layer of the CF heat sink, the hottest zone is almost at the center in the top layer, however, in the bottom layer, the hottest zone is still closer to the outlet with a slight shift towards the center (Figure 6b). The hottest zone in the TF heat sink is concentrated around the central region for both the bottom and the top layer as shown in Figure 6c. The TF heat sink showed lower thermal resistance and better temperature uniformity as compared to the PF and CF heat sinks in the previous work [24]. The better performance of the TF heat sink was attributed to the ubiquitous heat transfer present all over the heat sink as a result of quarterly opposite flow in the same layer and the transversal flow arrangement between the bottom and the top layers. The detailed analysis was presented by Ansari and Kim [24]. The wavy pattern observed in the temperature contours was due to the difference in conduction heat transfer in the heat sink walls and convection heat transfer in the channels. The total thermal resistance and maximum temperature increase at the bottom surface of the heat sinks without hotspots are listed in Table 7. The TF heat sink showed a 7.5% and 0.1% lower thermal resistance (and also temperature increase at the bottom surface) as compared to the PF and CF heat sinks, respectively. A variation in total pressure drop was observed among the PF, CF, and TF heat sinks even though all the design parameters were exactly the same [24]. This variation in the pressure drop was due to the consideration of variable fluid properties. The CF heat sink showed the lowest total pressure drop among the three heat sinks due to lowest average fluid temperature in both the layers [24].
The local temperature contours for the three heat sink designs (PF, CF, and TF) of the HS1 scheme (PF-HS1, CF-HS1, and TF-HS1), at total flow rate of 0.01 kg/s, are shown in Figure 7. The temperature contours are plotted at the bottom (y/Ly = 0) and middle height of the heat sink (y/Ly = 0.5). The comparison among the maximum temperatures at the hotspots of the PF, CF, and TF with the HS1 scheme is shown in Figure 8. The PF heat sink has a large temperature gradient in the direction of flow (Figure 6a), which causes overcooling of the hotspots (S1, S2, …) close to the inlet (x = 0) and undercooling of the hotspots (S11, S10, …) near the outlet (x = Lx) creating large temperature difference between hotspots (Figure 7a). Due to the opposite flow arrangement in the top layer, the CF heat sink shows lower temperatures at the hotspots, S10 and S11, as compared to the PF heat sink (Figure 8). The TF heat sink shows relatively uniform temperature distribution all over the heat sink (Figure 6c). Due to the uniform temperature distribution in the TF heat sink, the hotspots near x = 0 are not overcooled (exhibits higher temperatures than the PF and CF heat sinks) and the hotspots near x = Lx are not undercooled (exhibits lower temperatures than the PF and CF heat sinks) (Figure 8). The maximum temperatures in the PF and CF heat sinks are observed at the hotspot S10, whereas the maximum temperature in the TF heat sink is found at the hotspot S4, as shown in Figure 7 and Figure 8. Table 8 shows the maximum temperature rise at the hotspots, temperature variation between the hotspots, and the total thermal resistance for the HS1 scheme. The PF heat sink exhibits the highest temperature variation of 21.959 K among the hotspots. The CF heat sink exhibits a temperature variation of 17.692 K, which is 24.1% lower than that of the PF heat sink. Due to ubiquitous temperature distribution, the TF heat sink shows the lowest temperature variation among the three heat sinks (14.980 K), which is 46.6% and 18.1% lower than those of the PF and CF heat sinks, respectively. Along with significantly lower temperature variation among hotspots, the TF heat sink exhibits 9.3% and 6.2% lower thermal resistance (or maximum temperature rise) as compared to those of the PF and CF heat sinks, respectively.
Figure 9 shows the local temperature contours for the three heat sink designs (PF, CF, and TF) of the HS2 scheme (PF-HS2, CF-HS2, and TF-HS2), at total flow rate of 0.01 kg/s. The temperature contours are plotted at the bottom (y/Ly = 0) and middle height of the heat sink (y/Ly = 0.5). The maximum temperature rise at the hotspots, the temperature variation between the hotspots, and the total thermal resistance for the HS2 scheme are shown in Table 9. All the heat sinks with the HS2 scheme show the results similar to those for the HS1 scheme. The PF heat sink exhibits lower temperatures at the hotspots closer to the inlet due to overcooling, and due to undercooling the hotspots closer to outlet exhibits higher temperatures, producing 23.278 K temperature variation among the hotspots (Table 9). The temperature variation among hotspots in the CF heat sink is 17.278 K, which is 34.7% lower than that of the PF heat sink (Table 9). The TF heat sink exhibits the lowest temperature variation of 14.232 K, which is 63.6% and 21.4% lower than those of the PF and CF heat sinks, respectively (Table 9). Thermal resistance (or maximum temperature rise) is also lower for the TF heat sink, which is 13.3% and 6.3% lower than those of the PF and CF heat sinks, respectively (Table 9). The maximum temperatures in the PF and CF heat sinks occur at the hotspot S11 whereas the maximum temperature in the TF heat sink is found at the hotspot S8, as shown in Figure 9 and Figure 10.
The local temperature contours for all the heat sink designs (PF, CF, and TF) of the HS3 scheme (PF-HS3, CF-HS3, and TF-HS3), at total flow rate of 0.01 kg/s are shown in Figure 11. The temperature contours are plotted at the bottom (y/Ly = 0) and middle height of the heat sink (y/Ly = 0.5). The results are consistent with those for the HS1 and HS2 schemes. The TF heat sink exhibits the lowest temperature variation of 17.106 K, which is 48.0% and 16.1% lower than those of the PF (25.316 K) and CF (19.854 K) heat sinks, respectively (Table 10). The thermal resistance (or maximum temperature rise) of the TF heat sink is 11.7% and 6.4% lower than those of the PF and CF heat sinks, respectively (Table 10). The maximum temperatures at the hotspots are shown in Figure 12, which exhibits the highest temperatures at the hotspot S11 for the PF and CF heat sinks and at the hotspot S9 for the TF heat sink.
Figure 13 shows the temperature variations among the hotspots for different total flow rate (in a range from 0.006 to 0.014 kg/s) for all the heat sink designs with the HS1, HS2, and HS3 schemes. The TF heat sink exhibits the lowest temperature variations for all the three hotspot schemes (HS1, HS2, and HS3) throughout the flow rate range, while the PF heat sink exhibits the highest temperature variations (Figure 13). The PF heat sink exhibits the highest decrease in the temperature variation with the increase in the total flow rate, whereas the CF and TF heat sinks show similar trends (Figure 13).
Although all the heat sinks (PF, CF, and TF) were designed using exactly same microchannels, a slight variation in the pressure drop is observed in Table 11. The variation in the pressure drop was caused by the temperature-dependent fluid properties. With the increase in fluid temperature, the pressure drop reduces [47,48]. Table 11 shows the pressure drop and average fluid temperature variation in each layer of heat sink at the total flow rate of 0.01 kg/s. The pressure drops in the bottom layer (Δpbl) are lower than that of the top layer (Δptl) in all the cases, because the average fluid temperatures in the bottom layer are higher than those of the top layer (Table 11). The total pressure drops in all the heat sinks (PF, CF, and TF)—with all the hotspot schemes (HS1, HS2, and HS3) at the total flow rate of 0.01 kg/s are shown in Figure 14. The lowest total pressure drop is exhibited by the CF heat sink for all the hotspots schemes, because of the highest average fluid temperature in each layer as compared to those of the PF and TF heat sinks (Figure 14 and Table 11). On the other hand, the PF heat sink exhibits the lowest fluid temperature in each layer resulting in the highest total pressure drop. The relative pressure drop variation among the three heat sink designs is only about 1% for all the hotspot schemes at the total flow rate of 0.01 kg/s. The total pumping power is minimum for the CF-HS2, and maximum for the PF-HS1, but the total variation is less than 1 mW as shown in Table 11.

5. Conclusions

In this study, non-uniform heating with random hotspots was introduced for the analysis of double-layer microchannel heat sinks. Eleven hotspots were randomly generated using LHS to constitute three sets of hotspot schemes for three double-layer heat sinks: two heat sinks from the double-layer parallel-channel design (PF and CF) and one heat sink from the cross-channel design (TF). The entire heat sink was comprised of two layers with 20 microchannels in each layer. All the microchannels had same dimensions. The heat flux, area, and location (x and z) of the hotspots were selected as the design parameters. Three-dimensional conjugate heat transfer analysis was performed using 3-D Navier–Stokes equations. The thermophysical properties of water (i.e., the density, specific heat, thermal conductivity, and dynamic viscosity) were allowed to change with temperature. Polynomial correlations between the thermophysical properties of water and temperature in the property table were developed and verified in this work. The entire heat sink was selected as computational domain. The total flow rate for analysis was varied in a range from 0.006 to 0.014 kg/s (Re range from 400 to 933). The numerical model was validated using experimental data for pressure drop (for rectangular channel single-layer heat sink) and average temperature variation at the base (for rectangular channel double-layer counter-flow heat sink). Among the heat sinks without hotspots, the lowest thermal resistance was exhibited by the TF heat sink. The TF heat sink design also exhibited the lowest thermal resistance, maximum temperature rise among the hotspots, and minimum temperature variation among the hotspots for all the three hotspot heating schemes (HS1, HS2, and HS3) at all the flow rates. The PF heat sink exhibited the largest change in the temperature variation with the change of the flow rate as compared to the CF and TF heat sinks for all the three heating schemes. The lowest pressure drop was exhibited by the CF heat sink (with and without hotspots) due to the highest average fluid temperature in both the layers. The maximum relative variation in pressure drop among the heat sinks under three heating schemes was only about 1% at total flow rate of 0.01 kg/s. From the overall estimation for all the three random heating schemes, it can be concluded that the TF heat sink can be the best choice for cooling of a microchip with multiple hotspots as it exhibits the lowest thermal resistance (also the lowest maximum temperature rise at the hotspots) and the lowest temperature variation among the hotspots. The random hotspot analysis presented in this work is expected to be useful to assess the performance of heat sinks with realistic heating conditions of advance multicore processors.

Acknowledgments

This work was supported by INHA UNIVERSITY Research Grant.

Author Contributions

Danish Ansari conceived the idea, performed the numerical calculations, analyzed the results, and drafted the manuscript; Kwang-Yong Kim impeccably revised and improved the manuscript for the final submission.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Abtotal area of the heat sink base (m2)
Abgarea of background base (excluding hotspot area) (m2)
Ahsarea of single hotspot (m2)
Acmicrochannel cross-sectional area (m2)
BLbottom layer
CFcounter flow
Cpspecific heat (J/kg·K)
Dhhydraulic diameter of the channel (m)
hconvconvective heat transfer coefficient (W/m2·K)
Hcmicrochannel height (m)
HS1, HS2, HS3three random hotspot layout schemes
kthermal conductivity (W/m·K)
LBlower bound
Lslength of hotspot (m)
Lx, Ly, Lzlength, total height, and total width of the heat sink, respectively (m)
m tot total mass flow rate (kg/s)
n ch total numbers of channels
Δpchpressure drop in one channel (Pa)
Ptottotal pumping power required for a heat sink (W)
PFparallel flow
q hs , q b , q bg heat fluxes applied at the hotspot, base, and background base, respectively (W/m2)
Qtottotal heat applied at the base (W)
ReReynolds number
Rth,tottotal thermal resistance (K/W)
Ttemperature (K)
TLtop layer
TFtransverse flow
uavgaverage fluid velocity (m/s2)
UBupper bound
V ¯ velocity vector
Wb1thickness of base at the bottom (m)
Wb2thickness of base between the bottom and top layer (m)
Wcmicrochannel width (m)
Wcpthickness of the cover plate (m)
Wswidth of a hotspot (m)
Wwwidth of microchannel walls (m)
X, Y, Zorthogonal coordinate system

Subscripts

avgaverage value
bbase of heat sink
blbottom layer
ffluid
HShotspot
ininlet
maxmaximum value
ssolid
tltop layer
tottotal
wwater
WHSwithout hotspot

Greek Symbols

µdynamic viscosity of fluid (kg/m·s)
ρdensity of fluid (kg/m3)

References

  1. Tiwari, V.; Singh, D.; Rajgopal, S.; Mehta, G.; Patel, R.; Baez, F. Reducing power in high-performance microprocessors. In Proceedings of the Design Automation Conference, San Francisco, CA, USA, 15–19 June 1998; pp. 732–737.
  2. Parkhurst, J.; Darringer, J.; Grundmann, B. From single core to multi-core: Preparing for a new exponential. In Proceedings of the IEEE/ACM International Conference on Computer-Aided Design (ICCAD ’06), San Jose, CA, USA, 5–9 November 2006; pp. 67–72.
  3. Semenov, O.; Vassighi, A.; Sachdev, M. Impact of self-heating effect on long-term reliability and performance degradation in CMOS circuits. IEEE Trans. Device Mater. Reliab. 2006, 6, 17–27. [Google Scholar] [CrossRef]
  4. Traub, M.C.; Longsine, W.; Truskett, V.N. Advances in nanolithography/imprint lithography. Ann. Rev. Chem. Biomol. Eng. 2016, 7, 583–604. [Google Scholar] [CrossRef] [PubMed]
  5. Merolla, P.A.; Arthur, J.V.; Alvarez-Icaza, R.; Cassidy, A.S.; Sawada, J.; Akopyan, F.; Jackson, B.L.; Imam, N.; Guo, C.; Nakamura, Y.; et al. A million spiking-neuron integrated circuit with a scalable communication network and interface. Science 2014, 345, 668–673. [Google Scholar] [CrossRef] [PubMed]
  6. Suszko, A.; El-Genk, M.S. Thermally anisotropic composite heat spreaders for enhanced thermal management of high-performance microprocessors. Int. J. Therm. Sci. 2016, 100, 213–228. [Google Scholar]
  7. Sauciuc, L.; Chrysler, G.; Mahajan, R.; Szleper, M. Air-cooling extension-performance limits for processor cooling applications. In Proceedings of the Ninteenth Annual IEEE Semiconductor Thermal Measurement and Management Symposium, San Jose, CA, USA, 11–13 March 2003; pp. 74–81.
  8. Tuckerman, D.B.; Pease, R.F.W. High-performance heat sinking for VLSI. IEEE Electron Device Lett. 1981, 2, 126–129. [Google Scholar] [CrossRef]
  9. Liu, D.; Garimella, S.V. Investigation of liquid flow in microchannels. J. Thermophys. Heat Transf. 2004, 18, 65–72. [Google Scholar] [CrossRef]
  10. Al-Rjoub, M.F.; Roy, A.K.; Ganguli, S.; Banerjee, R.K. Assessment of an active-cooling micro-channel heat sink device, using electro-osmotic flow. Int. J. Heat Mass Transf. 2011, 54, 4560–4569. [Google Scholar] [CrossRef]
  11. Steinke, M.E.; Kandlikar, S.G. Review of single-phase heat transfer enhancement techniques for application in microchannels, minichannels and microdevices. Int. J. Heat Technol. 2004, 22, 3–11. [Google Scholar]
  12. Ansari, D.; Husain, A.; Kim, K.-Y. Multiobjective optimization of a grooved micro-channel heat sink. IEEE Trans. Compon. Packag. Technol. 2010, 33, 767–776. [Google Scholar] [CrossRef]
  13. Husain, A.; Kim, K.-Y. Microchannel heat sink with designed roughness: analysis and optimization. J. Thermophys. Heat Transf. 2008, 22, 342–351. [Google Scholar] [CrossRef]
  14. Moon, M.-A.; Kim, K.-Y. Exergetic analysis for optimization of a rotating equilateral triangular cooling channel with staggered square ribs. Int. J. Fluid Mach. Syst. 2016, 9, 229–236. [Google Scholar]
  15. Vafai, K.; Zhu, L. Analysis of two-layered micro-channel heat sink concept in electronic cooling. Int. J. Heat Mass Transf. 1999, 42, 2287–2297. [Google Scholar] [CrossRef]
  16. Hung, T.-C.; Yan, W.-M.; Li, W.-P. Analysis of heat transfer characteristics of double-layered microchannel heat sink. Int. J. Heat Mass Transf. 2012, 55, 3090–3099. [Google Scholar] [CrossRef]
  17. Wei, X.; Joshi, Y.; Patterson, M.K. Experimental and numerical study of a stacked microchannel heat sink for liquid cooling of microelectronic devices. J. Heat Transf. 2007, 129, 1432. [Google Scholar] [CrossRef]
  18. Wong, K.-C.; Ang, M.-L. Thermal hydraulic performance of a double-layer microchannel heat sink with channel contraction. Int. Commun. Heat Mass Transf. 2016. [Google Scholar] [CrossRef]
  19. Osanloo, B.; Mohammadi-Ahmar, A.; Solati, A.; Baghani, M. Performance enhancement of the double-layered micro-channel heat sink by use of tapered channels. Appl. Therm. Eng. 2016, 102, 1345–1354. [Google Scholar] [CrossRef]
  20. Cheng, Y.J. Numerical simulation of stacked microchannel heat sink with mixing-enhanced passive structure. Int. Commun. Heat Mass Transf. 2007, 34, 295–303. [Google Scholar] [CrossRef]
  21. Kulkarni, K.; Afzal, A.; Kim, K.-Y. Multi-objective optimization of a double-layered microchannel heat sink with temperature-dependent fluid properties. Appl. Therm. Eng. 2016, 99, 262–272. [Google Scholar] [CrossRef]
  22. Leng, C.; Wang, X.-D.; Wang, T.-H.; Yan, W.-M. Optimization of thermal resistance and bottom wall temperature uniformity for double-layered microchannel heat sink. Energy Convers. Manag. 2015, 93, 141–150. [Google Scholar] [CrossRef]
  23. Hung, T.-C.; Yan, W.-M.; Wang, X.-D.; Huang, Y.-X. Optimal design of geometric parameters of double-layered microchannel heat sinks. Int. J. Heat Mass Transf. 2012, 55, 3262–3272. [Google Scholar] [CrossRef]
  24. Ansari, D.; Kim, K.-Y. Double-layer microchannel heat sinks with transverse-flow configurations. J. Electron. Packag. 2016, 138, 031005. [Google Scholar] [CrossRef]
  25. Herwig, H.; Wickern, G. The effect of variable properties on laminar boundary layer flow. Wärme Stoffübertragung 1986, 20, 47–57. [Google Scholar] [CrossRef]
  26. Li, Z.; Huai, X.; Tao, Y.; Chen, H. Effects of thermal property variations on the liquid flow and heat transfer in microchannel heat sinks. Appl. Therm. Eng. 2007, 27, 2803–2814. [Google Scholar] [CrossRef]
  27. Esmaeilzadeh, H.; Blem, E.; St. Amant, R.; Sankaralingam, K.; Burger, D. Power limitations and dark silicon challenge the future of multicore. ACM Trans. Comput. Syst. 2012, 30, 11. [Google Scholar] [CrossRef]
  28. Wojciechowski, B.; Nikodem, M.; Weber, P.; Zagrabski, M.; Kepa, K. Hardware microprocessor thermal emulation using synthetic heat sources and temperature sensors in FPGA. In Proceedings of the 2014 20th International Workshop on Thermal Investigations of ICs and Systems (THERMINIC), London, UK, 24–26 September 2014; pp. 1–5.
  29. Rotem, E.; Hermerding, J.; Cohen, A.; Cain, H. Temperature measurement in the Intel(R) CoreTM Duo Processor. arXiv 2007. [Google Scholar]
  30. Xu, G. Thermal modeling of multi-core processors. In Proceedings of the Tenth Intersociety Conference on Thermal and Thermomechanical Phenomena in Electronics Systems (ITHERM’06), San Diego, CA, USA, May 30–June 2006; pp. 96–100.
  31. Tavakkoli, F.; Ebrahimi, S.; Wang, S.; Vafai, K. Thermophysical and geometrical effects on the thermal performance and optimization of a three-dimensional integrated circuit. J. Heat Transf. 2016, 138, 082101. [Google Scholar] [CrossRef]
  32. Joint Electron Device Engineering Council. Failure Mechanisms and Models for Semiconductor Devices; JEDEC: South Arlington, VA, USA, 2010. [Google Scholar]
  33. Hamann, H.F.; Weger, A.; Lacey, J.A.; Hu, Z.; Bose, P.; Cohen, E.; Wakil, J. Hotspot-limited microprocessors: Direct temperature and power distribution measurements. IEEE J. Solid-State Circ. 2007, 42, 56–65. [Google Scholar] [CrossRef]
  34. Abdoli, A.; Jimenez, G.; Dulikravich, G.S. Thermo-fluid analysis of micro pin-fin array cooling configurations for high heat fluxes with a hot spot. Int. J. Therm. Sci. 2015, 90, 290–297. [Google Scholar] [CrossRef]
  35. Liu, C.-K.; Yang, S.-J.; Chao, Y.-L.; Liou, K.Y.; Wang, C.-C. Effect of non-uniform heating on the performance of the microchannel heat sinks. Int. Commun. Heat Mass Transf. 2013, 43, 57–62. [Google Scholar] [CrossRef]
  36. Hao, X.; Peng, B.; Xie, G.; Chen, Y. Efficient on-chip hotspot removal combined solution of thermoelectric cooler and mini-channel heat sink. Appl. Therm. Eng. 2016, 100, 170–178. [Google Scholar] [CrossRef]
  37. Drummond, K.P.; Weibel, J.A.; Garimella, S.V.; Back, D.; Janes, D.B.; Sinanis, M.D.; Peroulis, D. Evaporative intrachip hotspot cooling with a hierarchical manifold microchannel heat sink array. In Proceedings of the 2016 15th IEEE Intersociety Conference on Thermal and Thermomechanical Phenomena in Electronic Systems (ITherm), Las Vegas, NV, USA, 31 May–3 June 2016; pp. 307–315.
  38. Green, C.E.; Kottke, P.A.; Sarvey, T.E.; Fedorov, A.G.; Joshi, Y.; Bakir, M.S. Performance and integration implications of addressing localized hotspots through two approaches: Clustering of micro pin-fins and dedicated microgap coolers. In Proceedings of the ASME 2015 International Technical Conference and Exhibition on Packaging and Integration of Electronic and Photonic Microsystems collocated with the ASME 2015 13th International Conference on Nanochannels, Microchannels, and Minichannels, San Francisco, CA, USA, 6–9 July 2015.
  39. Feng, S.S.; Kim, T.; Lu, T.J. Thermal resistance analysis of pin-fin heat sinks under nonuniform impingement heating. J. Thermophys. Heat Transf. 2011, 25, 119–129. [Google Scholar] [CrossRef]
  40. The MathWorks Inc. MATLAB User’s Guide; The MathWorks Inc.: Natick, MA, USA, 2014. [Google Scholar]
  41. McKay, M.D.; Beckman, R.J.; Conover, W.J. A Comparison of three methods for selecting values of input variables in the analysis of output from a computer code. Technometrics 1979, 21, 239–245. [Google Scholar] [CrossRef]
  42. Harms, T.M.; Kazmierczak, M.J.; Gerner, F.M. Developing convective heat transfer in deep rectangular microchannels. Int. J. Heat Fluid Flow 1999, 20, 149–157. [Google Scholar] [CrossRef]
  43. ANSYS. CFX 15.0. In Solver Theory Guide; ANSYS Inc.: Canonsburg, PA, USA, 2013. [Google Scholar]
  44. Raw, M. Robustness of coupled algebraic multigrid for the Navier–Stokes equations. In proceedings of the 34th Aerospace Sciences Meeting and Exhibit, Aerospace Sciences Meetings, Reno, NV, USA, 15–18 January 1996.
  45. Bergman, T.L.; Lavine, A.S.; Incropera, F.P.; Dewitt, D.P. Fundamentals of Heat and Mass Transfer, 7th ed.; John Wiley & Sons: Hoboken, NJ, USA, 2011. [Google Scholar]
  46. Qu, W.; Mudawar, I. Experimental and numerical study of pressure drop and heat transfer in a single-phase micro-channel heat sink. Int. J. Heat Mass Transf. 2002, 45, 2549–2565. [Google Scholar] [CrossRef]
  47. Toh, K.C.; Chen, X.Y.; Chai, J.C. Numerical computation of fluid flow and heat transfer in microchannels. Int. J. Heat Mass Transf. 2002, 45, 5133–5141. [Google Scholar] [CrossRef]
  48. Incropera, F.P. Liquid Cooling of Electronic Devices by Single-Phase Convection; Wiley: New York, NY, USA, 1999. [Google Scholar]
Figure 1. Schematic diagrams: (a) parallel-channel design (PF and CF); (b) cross-channel design (TF); (c) geometric parameters; and (d) flow conditions for top layer (TL) and bottom layer (BL).
Figure 1. Schematic diagrams: (a) parallel-channel design (PF and CF); (b) cross-channel design (TF); (c) geometric parameters; and (d) flow conditions for top layer (TL) and bottom layer (BL).
Micromachines 08 00054 g001
Figure 2. Hotspot layout schemes: (a) HS1; (b) HS2; and (c) HS3.
Figure 2. Hotspot layout schemes: (a) HS1; (b) HS2; and (c) HS3.
Micromachines 08 00054 g002
Figure 3. Verification of polynomial correlations between thermophysical properties of water and temperature with the property data from Bergman et al. [45]: (a) density ( ρ w ); (b) specific heat ( C p , w ); (c) thermal conductivity ( k w ); and (d) dynamic viscosity ( μ w ).
Figure 3. Verification of polynomial correlations between thermophysical properties of water and temperature with the property data from Bergman et al. [45]: (a) density ( ρ w ); (b) specific heat ( C p , w ); (c) thermal conductivity ( k w ); and (d) dynamic viscosity ( μ w ).
Micromachines 08 00054 g003
Figure 4. Grid independency test (PF design with single hotspot at the center at a total flow rate of 0.006 kg/s).
Figure 4. Grid independency test (PF design with single hotspot at the center at a total flow rate of 0.006 kg/s).
Micromachines 08 00054 g004
Figure 5. Validations of numerical results compared with: (a) experimental and numerical results obtained by Wei et al. [17] in a double-layer counter-flow heat sink for average temperature variation at the bottom surface in the flow direction; and (b) experimental results obtained by Qu and Mudawar [46] for the pressure drop in a rectangular microchannel with a heat flux of 200 W/cm2.
Figure 5. Validations of numerical results compared with: (a) experimental and numerical results obtained by Wei et al. [17] in a double-layer counter-flow heat sink for average temperature variation at the bottom surface in the flow direction; and (b) experimental results obtained by Qu and Mudawar [46] for the pressure drop in a rectangular microchannel with a heat flux of 200 W/cm2.
Micromachines 08 00054 g005
Figure 6. Temperature contours at the base (y/Ly = 0) and at the middle height (y/Ly = 0.5) of the heat sinks without hotspots at total flow rate of 0.01 kg/s: (a) PF-WHS; (b) CF-WHS; and (c) TF-WHS.
Figure 6. Temperature contours at the base (y/Ly = 0) and at the middle height (y/Ly = 0.5) of the heat sinks without hotspots at total flow rate of 0.01 kg/s: (a) PF-WHS; (b) CF-WHS; and (c) TF-WHS.
Micromachines 08 00054 g006
Figure 7. Temperature contours at the base (y/Ly = 0) and at the middle height (y/Ly = 0.5) of the heat sinks with the hotspot scheme HS1 at the total flow rate of 0.01 kg/s: (a) PF-HS1; (b) CF-HS1; and (c) TF-HS1.
Figure 7. Temperature contours at the base (y/Ly = 0) and at the middle height (y/Ly = 0.5) of the heat sinks with the hotspot scheme HS1 at the total flow rate of 0.01 kg/s: (a) PF-HS1; (b) CF-HS1; and (c) TF-HS1.
Micromachines 08 00054 g007
Figure 8. Maximum temperatures at the hotspots of PF, CF, and TF heat sinks with the hotspot scheme HS1 at the total flow rate of 0.01 kg/s.
Figure 8. Maximum temperatures at the hotspots of PF, CF, and TF heat sinks with the hotspot scheme HS1 at the total flow rate of 0.01 kg/s.
Micromachines 08 00054 g008
Figure 9. Temperature contours at the base (y/Ly = 0) and at the middle height (y/Ly = 0.5) of the heat sinks with the hotspot scheme HS2 at the total flow rate of 0.01 kg/s: (a) PF-HS2; (b) CF-HS2; and (c) TF-HS2.
Figure 9. Temperature contours at the base (y/Ly = 0) and at the middle height (y/Ly = 0.5) of the heat sinks with the hotspot scheme HS2 at the total flow rate of 0.01 kg/s: (a) PF-HS2; (b) CF-HS2; and (c) TF-HS2.
Micromachines 08 00054 g009
Figure 10. Maximum temperatures at the hotspots of PF, CF, and TF heat sinks with the hotspot scheme HS2 at the total flow rate of 0.01 kg/s.
Figure 10. Maximum temperatures at the hotspots of PF, CF, and TF heat sinks with the hotspot scheme HS2 at the total flow rate of 0.01 kg/s.
Micromachines 08 00054 g010
Figure 11. Temperature contours at the base (y/Ly = 0) and at the middle height (y/Ly = 0.5) of the heat sinks with the hotspot scheme HS3 at the total flow rate of 0.01 kg/s: (a) PF-HS3; (b) CF-HS3; and (c) TF-HS3.
Figure 11. Temperature contours at the base (y/Ly = 0) and at the middle height (y/Ly = 0.5) of the heat sinks with the hotspot scheme HS3 at the total flow rate of 0.01 kg/s: (a) PF-HS3; (b) CF-HS3; and (c) TF-HS3.
Micromachines 08 00054 g011
Figure 12. Maximum temperatures at the hotspots of PF, CF, and TF heat sinks with the hotspot scheme HS3 at the total flow rate of 0.01 kg/s.
Figure 12. Maximum temperatures at the hotspots of PF, CF, and TF heat sinks with the hotspot scheme HS3 at the total flow rate of 0.01 kg/s.
Micromachines 08 00054 g012
Figure 13. Temperature variation among the hotspots for different total flow rates for the PF, CF, and TF heat sinks (total flow rate varied from 0.006 to 0.014 kg/s): (a) HS1; (b) HS2; and (c) HS3.
Figure 13. Temperature variation among the hotspots for different total flow rates for the PF, CF, and TF heat sinks (total flow rate varied from 0.006 to 0.014 kg/s): (a) HS1; (b) HS2; and (c) HS3.
Micromachines 08 00054 g013
Figure 14. Total pressure drops in the PF, CF, and TF heat sink designs for the HS1, HS2, and HS3 hotspot schemes at the total flow rate of 0.01 kg/s.
Figure 14. Total pressure drops in the PF, CF, and TF heat sink designs for the HS1, HS2, and HS3 hotspot schemes at the total flow rate of 0.01 kg/s.
Micromachines 08 00054 g014
Table 1. Microchannel dimensions (mm).
Table 1. Microchannel dimensions (mm).
WcHcWwWb1Wb2WcpLxLyLz
0.250.500.250.200.200.20101.610
Table 2. Hotspot parameters and their ranges.
Table 2. Hotspot parameters and their ranges.
x (mm)z (mm)Ahs (mm2)qhs (W/cm2)
LB0.750.750.25200
UB9.259.252.25600
Table 3. Hotspots for the HS1 scheme.
Table 3. Hotspots for the HS1 scheme.
Hotspotx (mm)z (mm)Ahs (mm2)qhs (W/cm2)
S10.806.220.41374
S21.561.362.06225
S32.744.511.03466
S43.282.380.73558
S53.908.421.26256
S65.124.920.30582
S75.843.440.80311
S86.807.540.49359
S97.062.081.44403
S108.188.561.85435
S118.945.861.54258
Table 4. Hotspots for the HS2 scheme.
Table 4. Hotspots for the HS2 scheme.
Hotspotx (mm)z (mm)Ahs (mm2)qhs (W/cm2)
S11.344.620.63526
S21.586.781.53211
S32.302.002.18325
S43.224.221.35321
S54.129.161.90343
S65.240.800.51576
S75.468.160.39386
S86.502.501.25437
S97.105.480.34495
S107.787.461.08289
S118.823.320.93547
Table 5. Hotspots for the HS3 scheme.
Table 5. Hotspots for the HS3 scheme.
Hotspotx (mm)z (mm)Ahs (mm2)qhs (W/cm2)
S11.344.620.63594
S21.586.781.53300
S32.302.002.18201
S43.224.221.35379
S54.129.161.90237
S65.240.800.51488
S75.468.160.39272
S86.502.501.25393
S97.105.480.34503
S107.787.461.08312
S118.823.320.93452
Table 6. Goodness-of-fit parameters for the temperature-dependent property correlations.
Table 6. Goodness-of-fit parameters for the temperature-dependent property correlations.
Equations R 2 R a d j 2 SSERMSE
ρw0.99970.99971.12000.2567
Cp,w0.99950.99931.93400.3477
kw0.99990.99983.732 × 10−64.686 × 10−4
µw1.00001.00009.979 × 10−112.497 × 10−6
Table 7. Total thermal resistances and maximum temperature increases at the base of the PF, CF, and TF heat sinks without hotspot (WHS) at the total flow rate of 0.01 kg/s.
Table 7. Total thermal resistances and maximum temperature increases at the base of the PF, CF, and TF heat sinks without hotspot (WHS) at the total flow rate of 0.01 kg/s.
Temperature (K)Heat Sink Designs
PF-WHSCF-WHSTF-WHS
ΔTmax,b (K)20.82719.39019.374
7.5% ↓ PF-WHS
0.1% ↓ CF-WHS
Rth,tot (K/W)0.20830.19390.1937
7.5% ↓ PF-WHS
0.1% ↓ CF-WHS
Table 8. Total thermal resistances and maximum temperature increases at the hotspots of the PF, CF, and TF heat sinks for the hotspot scheme HS1 at the total flow rate of 0.01 kg/s.
Table 8. Total thermal resistances and maximum temperature increases at the hotspots of the PF, CF, and TF heat sinks for the hotspot scheme HS1 at the total flow rate of 0.01 kg/s.
Temperature (K)Heat Sink Designs
PF-HS1CF-HS1TF-HS1
ΔTmax,S119.36522.48822.837
ΔTmax,S221.83725.33925.657
ΔTmax,S331.55834.93737.583
ΔTmax,S434.42637.60637.817
ΔTmax,S524.28927.23728.807
ΔTmax,S631.42433.36333.107
ΔTmax,S730.26331.47030.284
ΔTmax,S830.18930.86530.119
ΔTmax,S937.48237.33536.440
ΔTmax,S1041.32440.18036.498
ΔTmax,S1131.61829.02229.544
ΔTmin,spot19.36522.48822.837
ΔTmax,spot41.32440.18037.817
9.3% ↓ PF-HS1
6.2% ↓ CF-HS1
Variation (range) ΔTmax,spotΔTmin,spot21.95917.69214.980
46.6% ↓ PF-HS1
18.1% ↓ CF-HS1
Rth,tot (K/W)0.31750.30870.2906
9.3% ↓ PF-HS1
6.2% ↓ CF-HS1
Table 9. Total thermal resistances and maximum temperature increases at the hotspots of the PF, CF, and TF heat sinks for the hotspot scheme HS2 at the total flow rate of 0.01 kg/s.
Table 9. Total thermal resistances and maximum temperature increases at the hotspots of the PF, CF, and TF heat sinks for the hotspot scheme HS2 at the total flow rate of 0.01 kg/s.
Temperature (K)Heat Sink Designs
PF-HS2CF-HS2TF-HS2
ΔTmax,S128.24931.75836.417
ΔTmax,S219.76223.08823.741
ΔTmax,S328.79632.42033.293
ΔTmax,S428.93832.06533.594
ΔTmax,S532.99335.64735.742
ΔTmax,S635.62437.55936.109
ΔTmax,S729.15630.78828.710
ΔTmax,S837.88338.55837.973
ΔTmax,S931.29431.22931.101
ΔTmax,S1030.34429.36028.405
ΔTmax,S1143.04040.36635.285
ΔTmin,spot19.76223.08823.741
ΔTmax,spot43.04040.36637.973
13.3% ↓ PF-HS2
6.3% ↓ CF-HS2
Variation (range) ΔTmax,spot–ΔTmin,spot23.27817.27814.232
63.6% ↓ PF-HS2
21.4% ↓ CF-HS2
Rth,tot (K/W)0.32550.30530.2872
13.3% ↓ PF-HS2
6.3% ↓ CF-HS2
Table 10. Total thermal resistances and maximum temperature increases at the hotspots of the PF, CF, and TF heat sinks for the hotspot scheme HS3 at the total flow rate of 0.01 kg/s.
Table 10. Total thermal resistances and maximum temperature increases at the hotspots of the PF, CF, and TF heat sinks for the hotspot scheme HS3 at the total flow rate of 0.01 kg/s.
Temperature (K)Heat Sink Designs
PF-HS3CF-HS3TF-HS3
ΔTmax,S126.99530.34930.778
ΔTmax,S224.71428.18331.222
ΔTmax,S318.31221.67021.936
ΔTmax,S432.71336.12438.604
ΔTmax,S525.65628.06128.363
ΔTmax,S629.15231.10830.364
ΔTmax,S729.56330.70430.572
ΔTmax,S837.71138.41535.066
ΔTmax,S939.08639.53039.042
ΔTmax,S1028.36727.04526.425
ΔTmax,S1143.62841.52436.871
ΔTmin,spot18.31221.67021.936
ΔTmax,spot43.62841.52439.042
11.7% ↓ PF-HS3
6.4% ↓ CF-HS3
Variation (range) ΔTmax,spot–ΔTmin,spot25.31619.85417.106
48.0% ↓ PF-HS3
16.1% ↓ CF-HS3
Rth,tot (K/W)0.30230.28770.2705
11.7% ↓ PF-HS3
6.4% ↓ CF-HS3
Table 11. Average fluid temperatures (in each layer) and pressure drops (in each layer and total) of the PF, CF, and TF heat sinks with the hotspot schemes HS1, HS2, and HS3 at the total flow rate of 0.01 kg/s.
Table 11. Average fluid temperatures (in each layer) and pressure drops (in each layer and total) of the PF, CF, and TF heat sinks with the hotspot schemes HS1, HS2, and HS3 at the total flow rate of 0.01 kg/s.
DesignTavg,f,bl (K)Tavg,f,tl (K)Δpbl (Pa)Δptl (Pa)Δptot = Δpbl + Δptl (Pa)Ptot (W) × 10−3
PF-HS1308.623305.973126,707132,092258,79964.894
CF-HS1309.520307.014125,429130,695256,12464.224
TF-HS1309.514306.802125,477130,816256,29364.226
PF-HS2308.879306.166126,302131,747258,04964.706
CF-HS2309.748307.123125,105130,448255,55364.080
TF-HS2309.635306.905125,243130,621255,86464.158
PF-HS3308.789306.090126,422131,861258,28364.765
CF-HS3309.693307.127125,151130,469255,62064.097
TF-HS3309.705306.998125,132130,550255,68264.113

Share and Cite

MDPI and ACS Style

Ansari, D.; Kim, K.-Y. Performance Analysis of Double-Layer Microchannel Heat Sinks under Non-Uniform Heating Conditions with Random Hotspots. Micromachines 2017, 8, 54. https://doi.org/10.3390/mi8020054

AMA Style

Ansari D, Kim K-Y. Performance Analysis of Double-Layer Microchannel Heat Sinks under Non-Uniform Heating Conditions with Random Hotspots. Micromachines. 2017; 8(2):54. https://doi.org/10.3390/mi8020054

Chicago/Turabian Style

Ansari, Danish, and Kwang-Yong Kim. 2017. "Performance Analysis of Double-Layer Microchannel Heat Sinks under Non-Uniform Heating Conditions with Random Hotspots" Micromachines 8, no. 2: 54. https://doi.org/10.3390/mi8020054

APA Style

Ansari, D., & Kim, K. -Y. (2017). Performance Analysis of Double-Layer Microchannel Heat Sinks under Non-Uniform Heating Conditions with Random Hotspots. Micromachines, 8(2), 54. https://doi.org/10.3390/mi8020054

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop