Next Article in Journal
Perfusion System for Modification of Luminal Contents of Human Intestinal Organoids and Realtime Imaging Analysis of Microbial Populations
Previous Article in Journal
Dynamic Modeling and Experimental Validation of a Water Hydraulic Soft Manipulator Based on an Improved Newton—Euler Iterative Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation and Optimization of a Cross-Rib Micro-Channel Heat Sink

1
Institute of Microelectronics of Chinese Academy of Sciences, Beijing 100029, China
2
University of Chinese Academy of Sciences, Beijing 100029, China
*
Author to whom correspondence should be addressed.
Micromachines 2022, 13(1), 132; https://doi.org/10.3390/mi13010132
Submission received: 10 December 2021 / Revised: 11 January 2022 / Accepted: 12 January 2022 / Published: 14 January 2022

Abstract

:
This article presents a novel cross-rib micro-channel (MC-CR) heat sink to make fluid self-rotate. For a thermal test chip (TTC) with 100 w/cm2, the cross-ribs micro-channel were compared with the rectangular (MC-R) and horizontal rib micro-channel (MC-HR) heat sinks. The results show that, with the cross-rib micro-channel, the junction temperature of the thermal test chip was 336.49 K, and the pressure drop was 22 kPa. Compared with the rectangular and horizontal ribs heat sink, the cross-rib micro-channel had improvements of 28.6% and 14.3% in cooling capability, but the pressure drop increased by 10.7-fold and 5.5-fold, respectively. Then, the effects of the aspect ratio (λ) of micro-channel in different flow rates were studied. It was found that the aspect ratio and cooling performance were non-linear. To reduce the pressure drop, the inclination (α) and spacing (S) of the cross-ribs were optimized. When α = 30°, S = 0.1 mm, and λ = 4, the pressure drop was reduced from 22 kPa to 4.5 kPa. In addition, the heat dissipation performance of the rectangular, staggered fin (MC-SF), staggered rib (MC-SR) and cross-rib micro-channels were analyzed in the condition of the same pressure drop, MC-CR still has superior heat dissipation performance.

1. Introduction

With the rapid development of integrated circuits in terms of high integration and miniaturization, the thermal design power and heat flux of electronics are increasing. It is essential to dissipate the heat accumulated inside electronics [1]. According to statistics, 55% of electronic failures are caused by heat [2,3]. In the past few decades, heat dissipation has mainly been through natural and forced convection [4,5]. Among them, natural convection does not rely on external force; fluid flows because of the temperature difference inside fluid, and the heat generated by the microelectronics is taken away and released. This method is low in cost and is only suitable for devices with low power density and no strict temperature requirements, because the heat dissipation capacity of this method is very limited [6]. Common forced convection mainly includes air-cooling and water-cooling. Air-cooling accelerates the removal of heat, which increases the velocity of the air through the fan. Additionally, water-cooling is realized by installing a liquid cooling plate on the heat source, cooperating with a heat exchanger and a pump to dissipate heat by means of fluid flow. Compared with water-cooling, the cost of air-cooling is low and the performance is stable, although the cooling capacity is limited; therefore, it is widely used in scenarios with low power density and no specific wind speed requirements. However, the heat transfer coefficient of water is 1000~15,000 W/m2·k in the case of forced convection, which is more than 100-fold that of air. In addition, water-cooling is quieter.
Nowadays, micro-channel heat sinks are a mainstream liquid cooling technology, with the concept first proposed by Tuckerman and Pease [7]. The technology is mainly designed to continuously take away the heat generated by microelectronics through the circulation of fluid. With satisfactory heat transfer capability, compact structure, and relatively light weight, micro-channel heat sinks can be regarded as a suitable and feasible solution to meet the heat dissipation requirements [8,9]. The key factors affecting the heat dissipation performance of the micro-channel are the type of fluid and the structure of the micro-channel. There has been much research on the shape of the micro-channel. Wang et al. [10] compared the flow characteristics and heat transfer performance of triangular, trapezoidal and rectangular micro-channels through numerical simulations. The result showed that rectangular micro-channels have the lowest thermal resistance and the best heat dissipation performance, followed by trapezoidal micro-channels and triangular micro-channels. In addition, rectangular micro-channels have the lowest thermal resistance in the aspect ratio range of 8.904–11.442. Lin et al. [11] proposed a novel wavy micro-channel heat sink that changes the wavelength and amplitude along the flow path. When the wavelength of the wave unit decreases or the amplitude increases, the wavy micro-channel heat sink has a lower thermal resistance and exhibits better overall performance. Park et al. [12] optimized the critical dimensions of wave micro-channel heat sinks with grooves, such as the width, depth and spacing of the grooves, reducing the thermal resistance and friction coefficients by 1.55% and 3%, respectively. Huang et al. [13] conducted a numerical study on the heat transfer characteristics of rectangular parallel grooves, rectangular staggered grooves, and trapezoidal staggered groove micro-channel heat sinks. The results show that under conditions of constant power, the heat transfer performance of the trapezoidal staggered groove micro-channel heat sink is the best.
In addition, changing the flow state in the micro-channel has also attracted a large number of researchers. Ghani et al. [14] studied the fluid flow and heat transfer characteristics (MC-SCRR) of a novel sinusoidal cavity rectangular rib micro-channel heat sink. The performance of rectangular rib micro-channels (MC-RR) and sinusoidal cavity micro-channels (MC-SC) was compared and analyzed. The results show that the thermal performance of MC-SCRR is better than that of MC-RR and MC-SC. Chai et al. [15,16,17,18,19,20] studied the effect of ribs in different arrangements on the heat dissipation capacity of the micro-channel. The study found that both aligned ribs and offset ribs can increase the heat transfer area and enhance fluid disturbance. Xie et al. [21] investigated the influence of the inclined angle of the pin fins. When inserting pin fins with an inclined angle of 30° in the MCHS, intense secondary flow was generated to facilitate disturbance flow. Deng et al. [22] and Carlo et al. [23] studied the performance of double-layered micro-channel heat sinks (DL-MCHS) with different cross sections and an unequal number of micro-channels in two layers. The result shows that DL-MCHS significantly reduced wall temperature and thermal resistance, as well as pressure drop and pumping power compared with single-layered micro-channel heat sinks (SL-MCHS). In addition, when the number of micro-channels in the top layer is smaller than that in the bottom layer, there may be lower thermal resistance and better temperature uniformity. Jing et al. [24] investigated the thermal characteristics of an interleaved double-layer micro-channel heat sink with upper and lower micro-channels offset along the width direction. It was found that when the rib thickness is small, the thermal resistance increases with the increase in the offset value, and when the rib thickness is large enough, the thermal resistance decreases first and then increases with the increase in the offset value. Therefore, it will significantly enhance heat transfer but also increase pressure drop, especially aligned ribs. In summary, increasing the heat exchange area, reducing convective thermal resistance, and changing the flow state of the fluid can all improve the heat dissipation performance of the micro-channel.
This paper proposes a novel micro-channel that makes the fluid rotate. The heat dissipation performance of MC-CR was evaluated through numerical analysis and simulation; then, the influence of the aspect ratio of the micro-channel on the heat dissipation performance was analyzed. Then, the pressure drop was reduced by optimizing the inclination angle, spacing of the cross-rib, and aspect ratio of the micro-channels, which provided a reference for the optimization of the MC-CR heat sink. Finally, compared with MC-R, MC-SF, MC-SR, and MC-CR at the same pressure drop, MC-CR still had superior heat dissipation performance.

2. Geometrical Model and Method

In this study, the application scenario of the heat sink was an FCBGA thermal chip with a size of 25 × 25 × 0.75 mm. In order to realize the miniaturization design of the package, the heat sink was integrated inside the lid. A layer of Tim with a thickness of 100 μm was coated between the thermal test chip (TTC) and the heat sink for connection. As shown in Figure 1b, the size of the MC-CR heat sink is the same as that of the TTC, and there are 25 channels in total. Figure 1a shows a cross-sectional view of the assembly. Figure 1c shows the internal structure of MC-CR. The two crossed squares with an inclination angle of 30°, which were called cross-ribs, were added on the inner wall of the straight micro-channel. The cross-rib makes the fluid self-rotate along it, which promotes the mixing of the upper and lower parts of the fluid, and maximizes the fluid’s ability to absorb heat. Additionally, it also enhances the disturbance of the fluid, breaks the boundary layer, and significantly improves the heat dissipation. Figure 2 shows the assembly process. First, the aluminum plate was machined with a milling cutter to mill out the grooves, then the cross-ribs were manually assembled, and finally, they were glued to the manifold. However, the processing difficulty is inversely proportional to the inclination angle of the cross-ribs. Therefore, the inclination angle should be as large as possible on the premise of ensuring the heat dissipation performance.

2.1. Geometrical Model

This paper studied the influence of the micro-channel structure on heat dissipation performance, and used ANSYS fluent to solve the steady-state solution of the micro-channel heat dissipation performance. To ensure the uniqueness of variables, the width, length, depth of the heat source chip, Tim, base, and micro-channel remained the same. Figure 3 shows the basic unit model of the three micro-channels and their corresponding cross-sectional diagrams. Figure 3e shows a horizontal rib micro-channel (MC-HR) heat sink. The horizontal ribs on both sides were placed in the center of the wall and were in contact with each other, and its liquid–solid contact area was equal to MC-CR. The parameters of the positions marked in Figure 3 are tabulated in Table 1, and the ratio of the width of the fluid to the channel in the three micro-channels was 1/2 (Wc = 2 Wm).

2.2. Numerical Model and Parameter Definition

We assumed that the flow and heat transfer were in steady state, the coolant was an incompressible and Newtonian fluid, and the heat conducted from the chip was only dissipated by coolant convection; thus, heat loss to the ambient was ignored. With these assumptions, the governing equations are as follows:
Continuity equation [25]:
u 1 x + u 2 y + u 3 z = 0
Momentum equations
ρ ( u 1 u 1 x + u 2 u 1 y + u 3 u 1 z ) = p x + u 2 u 1 ρ ( u 1 u 2 x + u 2 u 2 y + u 3 u 2 z ) = p y + u 2 u 2 ρ ( u 1 u 3 x + u 2 u 3 y + u 3 u 3 z ) = p z + u 2 u 3
Energy equation for the coolant
ρ c p l ( u 1 T x + u 2 T y + u 3 T z ) = λ l 2 T
where u1, u2 and u3 are the velocity in the x, y and z axis directions, respectively. ρ is the density of water.
For turbulence flow, the turbulence kinetic energy, k, and its rate of dissipation, ε, are obtained from the following transport equations:
t ( ρ k ) + x i ( ρ k v i ) = x j [ ( μ + μ t σ k ) k x j ] + G k + G b ρ ε Y M + S k t ( ρ ε ) + x i ( ρ ε v i ) = x j [ ( μ + μ t σ ε ) ε x j ] + C 1 ε ε k ( G k + C 3 ε G b ) C 2 ε ρ ε 2 k + S ε
In these equations, Gk represents the generation of turbulence kinetic energy due to the mean velocity gradients, Gb is the generation of turbulence kinetic energy due to buoyancy, YM represents the contribution of the fluctuating dilatation in compressible turbulence to the overall dissipation rate. C, C, and C are constants. σk and σε are the turbulent Prandtl numbers for k and ε, respectively.
G k = ρ v i v j ¯ v j x i
G b = β g i μ t T P r t x i
Y M = 2 ρ ε M t 2
where β is the coefficient of thermal expansion, gi is the component of the gravitational vector in the i direction, and Mt is the turbulent Mach number.
In general, the Nusselt number, thermal resistance and convective heat transfer coefficient are used to evaluate the heat dissipation performance of micro-channels.
The Reynolds number (Re) is a dimensionless number used to characterize fluid flow and it’s defined as follows
R e = ρ v d μ
D = 4 A P = 2 W c H c W c + H c
where v is the flow velocity, D is the hydraulic diameter of the micro-channel, μ is the kinematic viscosity coefficient of the fluid, A is the cross-sectional area of the flow, and P is the wet perimeter.
The aspect ratio (λ) is defined as follows:
λ = H c / W c
Water is a Newtonian fluid; therefore, the calculation equation of the pressure drop in the tube is for a turbulent Newtonian fluid, and the friction factor is defined as follows [26]:
f = 2 D Δ P v 2 L ρ = P o R e
where L is the characteristic length of the tube and ΔP is the pressure drop, which can be calculated with Equation (12)
Δ P = P i n P o u t
In laminar flow, the Poiseuille number equation in friction factor form is defined as follows [27]:
P o = 24 ( 1 1.13553 λ + 1.9467 λ 2 1.7012 λ 3 + 0.9564 λ 4 0.2537 λ 5 )
The Nusselt number (Nu) is a dimensionless number that expresses the intensity of convective heat transfer, and is defined as follows [28]:
N u = h a v e L / k p
where kp is the thermal conductivity of the plate material and have is the average convective heat transfer coefficient, which is defined as follows [12]:
h a v e = q / A c ( T s T i o )
where q is the heat flux, Ac is the area of convective heat transfer between water and solid, Ts is the average temperature at the solid, and Tio is the average of the inlet and outlet temperature.
The equivalent thermal resistance (Req) is defined as [29]:
R eq = R tot R c h i p R t i m
R = L p k A p
R tot = T 2 T 1 P tot = T 2 T 1 q A p
where T2 is the average temperature of the heat source, T1 is the water inlet temperature, Ptot is the total power of the heat source, Rtot is the total thermal resistance, Lp is the thickness of the plate, and Ap is the cross-sectional area of the plate perpendicular to the direction of heat flow.

2.3. Meshing Independence and Validation

Eliminating the influence of meshing on the results ensures the reliability of the simulation results. Taking the MC-R heat sink as an example, the average temperature of five meshing methods is compared. Table 2 lists the minimum element size, the total number of grids, and the average temperature of the heat source for the five meshing methods. It can be seen that with the refinement of the grid, the simulation results are more accurate and more consistent with the theoretical values. However, the difference between the results of Method Ⅳ and Method Ⅴ is 0.87%. Therefore, in order to ensure accurate simulation results and higher calculation efficiency, the subsequent simulations all used Method Ⅳ (the minimum unit size was 0.04 mm and the number of grids is 187,500). Figure 4 shows the mesh model and residual curve of Method Ⅳ.
To verify the reliability of numerical results, the numerical results of the present method were compared with experimental data [30]. As shown in Figure 5, the numerical results agreed well with the experimental data. Additionally, the maximum error was less than 5% in different Reynolds numbers, which is within the acceptable range.

2.4. Setup and Boundary Conditions

At a reference temperature of 300 K, the kinematic viscosity of water is 8.937 × 10−5 m2/s [31]. The flow velocity was relatively uniform in the MC-R heat sink, and the Re was calculated as 2237.9 (Equations (8) and (9)); therefore, the laminar flow model was used to solve it. However, in the heat sink with special-shaped structures of MC-CR and MC-HR heat sinks, the fluid in the channels no longer maintained an orderly stratified flow, so turbulence models were used to solve it. In order to simplify the simulation setup, it was assumed that the microelectronic chip heated up uniformly. The details of the boundary conditions are as follows,
  • The inlet adopted velocity-inlet flow rate was 2 m/s;
  • The heat flux of the TTC was 100 w/cm2;
  • The material of fluid was water (0.6 W/m·K), the material of TTC was Si (148 W/m·K), the material of micro-channel was Al (202.2 W/m·K), and the thermal conductivity of Tim was 5 W/m·K;
  • In the adopted pressure-outlet, the pressure was set to atmospheric pressure.

3. Result and Discussion

3.1. Effect of Ribs

In the same conditions, the heat dissipation performance of the micro-channel heat sink was characterized by the temperature of the heat source and pressure drop. It can be seen from Figure 6a that the temperatures of the TTC in the MC-CR and MC-HR were 14.64 K and 7.3 K lower than that in the MC-R, respectively. Therefore, the performance improvement in MC-HR was 50% that of MC-CR. The MC-CR and MC-HR had the same liquid–solid contact area, but the fluid was still in laminar flow. Therefore, the effects of flow disturbance and contact area on the heat dissipation performance of MC-CR are equal.
As shown in Figure 7a,c,e, the maximum temperature difference of MC-CR is 0.12 K, and that of MC-R is 2.46 K. Therefore, MC-CR has excellent cooling capacity and heat dissipation uniformity. Figure 7b,d,f shows the sidewall temperature distribution of MC-CR and MC-R heat sinks, respectively. It can be seen that in the direction of flow, the thermal boundary layer of MC-R gradually thickens (Equations (19) and (20)), which will greatly reduce the heat conduction at the contact interface; the tendency of the thermal boundary layer of MC-CR to thicken is broken due to the disturbance of the fluid. At the same time, the cross-ribs make the fluid flow upward after absorbing heat, and the cold fluid on it will be replenished in time, so that the upper and lower fluids can fully absorb the heat. For MC-HR, although the horizontal ribs did not make the boundary layer thinner, they played a balanced role, which reduced the temperature difference of the TTC to 0.53 K.
The thickness of the boundary layer (δ(x)) is defined as follows:
δ 1 ( x ) x / Re x
δ 2 ( x ) x /   5 Re x
where δ1(x) and δ2(x) are the boundary layer thickness of the laminar and turbulent, respectively; x is the distance of fluid flow, and Rex is the Reynolds number at x.
Figure 6a illustrates that the pressure drop of MC-CR is 10.7-fold larger than that of MC-R, and the MC-HR is 1.94-fold larger than that of MC-R. It can be seen from Figure 8a–f that the pressure is mainly concentrated where the rib is perpendicular to the flow direction. Therefore, the cross-ribs severely block the fluid flow and forcibly change the flow direction, which is one of the reasons for increasing the pressure drop. This also shows that the improvement in MC-CR’s heat dissipation performance comes at the expense of pressure drop. However, this is in contrast to the concept of energy saving, so this article provides some ideas to reduce the pressure drop in Section 4.
Figure 9 shows the streamlines of the fluids and velocity contours in the MC-R, MC-CR, and MC-HR heat sinks. It can be seen that the fluid flow in the MC-R has obvious stratification. The velocity of the center is the largest. The position near the wall is the velocity boundary layer in which the flow velocity is almost 0, due to the viscous force of the fluid, which cannot be ignored. The thermal conductivity of water is 0.6 W/m·K; therefore, the boundary layer would hinder the heat conduction. In the MC-CR, the cross-ribs structure makes the fluid self-rotate, which greatly disrupts the tendency of flow rate stratification and increases the flow rate and disturbance. Both of these methods increase the have and Nu, and thin the boundary layer. As shown in Figure 6b, the MC-CR has the largest Nu and the smallest Req. In the MC-HR, although the horizontal ribs increase the heat transfer area and have (Figure 6c), they do not improve the flow rate stratification (Figure 9c). Above and below the horizontal ribs, the flow velocity boundary layer still exists.

3.2. Effect of Aspect Ratio λ on Heat Dissipation Performance

In order to further study the relationship between the heat dissipation performance and the v, Wc and Hc, this study designed eight MC-CRs with different aspect ratios (Table 3), and the heat dissipation performance of eight MC-CRs with different aspect ratios was analyzed at flow rates of 1 m/s, 2 m/s and 3 m/s.
Figure 10 shows the Tave, have, Nu, Req and ΔP at different flow rates. As the flow rate increased, the improvements in heat dissipation performance of the microchannel tended to be saturated, but the pressure drop greatly increased. In addition, the thickness of the ribs also changed with the aspect ratio for MC-R. Therefore, it needed to be compared in groups. In the first, second, third, and fifth experiments, Hc remained constant, Wc decreased, L decreased, the λ increased, and the fluid–solid contact area decreased. This phenomenon weakened the turbulence of the fluid and increased the heat transfer coefficient and Nusselt number, which also indicates that the heat convection was weakened; therefore, the heat dissipation performance became worse. In addition, D decreased, the pressure drop should have increased, but for MC-CR, the pressure drop became smaller with the resistance to the fluid dropping due to the thinner thickness of the cross-ribs. In the fourth, sixth, seventh, and eighth experiments, Wc remained constant, Hc increased, L increased, λ increased, and the fluid–solid contact area increased. This phenomenon enhanced the turbulence of the fluid and increased the heat transfer coefficient and Nusselt number. Therefore, the heat dissipation performance was improved. In addition, D increased, so the pressure drop decreased. In summary, when λ was 4 and the flow rate was 2 m/s, the MC-CR had better heat dissipation performance and a lower pressure drop.

4. Parameter Optimization

This section focuses on optimizing the above-mentioned MC-CR. Even if MC-CR has better heat dissipation performance, its pressure drop must be reduced as much as possible. Therefore, the inclination angle, α, and spacing of the cross-ribs and the aspect ratio, λ, of the micro-channel were optimized.
When the flow velocity, tube length and fluid density are constant, reducing the inclination angle of the rifling structure and increasing the spacing will decrease the pressure drop (Equation (11)) because of the reduced friction coefficient; the effect of the aspect ratio on pressure drop was analyzed in the previous section, and when λ was 4, the pressure drop was the smallest.

4.1. The Inclination Angle (α) of the Cross-Ribs

Figure 8a–f shows that the pressure drop was concentrated where the rib was perpendicular to the flow direction. Therefore, keeping the other parameters consistent, this section analyzes the pressure drop and the average temperature of the TTC when the inclination angles, α, of the cross-ribs are 45°, 40°, 35°, 30°, and 25°.
The resistance coefficient of the fluid in the MC-CR heat sinks is derived from the resistance of the cross-ribs structure and the wall friction. Equation (11) illustrates that the resistance coefficient is proportional to the pressure loss. Figure 11 shows that as the inclination angle, α, decreases, the heat dissipation performance becomes worse, and the pressure drop increases. The inclination angle, α, was reduced from 45° to 25°, and the average temperature rose by 0.3% in total; therefore, the change in inclination angle, α, has little effect on the heat dissipation performance. However, when the inclination angle, α, was reduced from 45° to 25°, the pressure drop was reduced by 97.7% in total, which is significant.

4.2. The Spacing (S) of the Cross-Rib

Through the preliminary simulation, it was found that fluid rotation can also be guaranteed in the condition where cross-ribs structure does not make contact. Therefore, in the condition that other parameters remain constant, this section analyzes the pressure drop and the average temperature of the heat source when the cross-rib spacing, S, is 0 mm, 0.05 mm, 0.1 mm, 0.15 mm, and 0.2 mm.
It can be seen from the streamline (Figure 9) that larger spacing in the cross-rib structure allows part of the fluid to flow directly through the middle, which reduces the area where the fluid is blocked, therefore reducing the pressure drop. Figure 12 shows that as the spacing was increased from 0 mm to 0.2 mm, the pressure drop was reduced by 86.3% in total. However, at the same time, the increase in spacing also slowed down the rotation of the fluid, causing an average temperature increase of 0.84 K.
Due to process limitations, as α decreases, it becomes more difficult for the milling cutter to mill the groove. Therefore, without affecting the heat dissipation performance, in the condition where α was 30°, the S was 0.1 mm, and λ was 4, the pressure was reduced from 22 kPa to 4.5 kPa.

4.3. Effect of ribs in Same Pressure Drop

At the same inlet flow rate, comparison between MC-CR, MC-R, MC-SF and MC-HR showed that MC-CR had better heat dissipation performance. However, that superior performance comes with the expense of high pressure drop in the structure. Therefore, this section compares the heat dissipation performance of MC-CR, MC-SF and MC-SR heat sinks with the same pressure drop (4.5 kPa). Figure 13a shows the unit model of the MC-SF, which diverted the water continuously to reduce the thickness of the boundary layer. As shown in Figure 13b, the MC-SR is a staggered rib micro-channel, as proposed by Chai et al. [15].
It can be seen from Figure 14a that the MC-CR still has obvious advantages at the same pressure drop. Compared with MC-R, the cooling performance of MC-CR increased by 14.25%, whereas MC-SF and MC-SR increased by 2.16% and 3.74%, respectively. Figure 14b shows that the Nu of MC-R is larger than that of MC-SF and MC-SR, but its heat dissipation performance is worse. The reason is that the heat dissipation performance is determined by both heat conduction and convection, and the boundary layer with poor thermal conductivity of MC-R makes heat conduction weaker. One of the advantages of MC-CR is that the boundary layer is well broken; thus, the heat dissipation performance has been greatly improved.

5. Conclusions

In this paper, a new cross-rib micro-channel heat sink has been proposed. The heat dissipation performance of the MC-CR was studied through numerical analyses, and the key parameters were optimized. The conclusions are as follows:
(1)
Compared with MC-R, the heat dissipation performance of MC-CR increased by 28.6%, and the pressure drop increased by 10.7-fold. Comparing with MC-HR, which has the same liquid–solid contact area, the heat dissipation performance increased by 14.3%. Therefore, the effects of fluid disturbance and contact area on the heat dissipation performance of MC-CR are equal;
(2)
Increasing the aspect ratio will increase the liquid–solid contact area, which is conducive to improving the heat dissipation performance. However, for MC-R, the thickness of the ribs will also change with the aspect ratio. Therefore, in the condition with constant Hc, the heat dissipation performance is inversely proportional to the aspect ratio; when the Wc is constant, the heat dissipation performance is proportional to the aspect ratio;
(3)
As the inclination angle, α, decreases, the average temperature rises by 5.7% overall, although the overall pressure drop is reduced by 97%. Meanwhile, the heat dissipation performance deteriorates by 5.6%, and the pressure drop is reduced by 86.3% as the spacing, S, increases. The heat dissipation performance is increased by 6.6%, and the pressure drop is reduced by 78.7% with the spacing, S, increases. In the condition where α = 30°, S = 0.1 mm, and λ = 4, the pressure drop is 4.5 kPa;
(4)
With the same pressure drop, the cooling performance of MC-CR increased by 14.25%, and the convection heat transfer coefficient increased by 2.27-fold compared with MC-R. The cooling performance of MC-CR increased by 10.51%, and the convection heat transfer coefficient increased by 2.3-fold compared with MC-SR.

Author Contributions

Conceptualization, H.C. and C.C.; validation, H.C., C.C., G.S. and C.Y.; resources, B.J. and R.L.; writing—original draft preparation, H.C.; writing—review and editing, H.C. and F.H.; project administration, Y.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by [the National Natural Science Foundation of China] grant number [62104249].

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Laguna, G.; Vilarrubí, M.; Ibañez, M.; Betancourt, Y.; Illa, J.; Azarkish, H.; Amnache, A.; Collin, L.M.; Coudrain, P.; Fréchette, L.; et al. Numerical parametric study of a hotspot-targeted microfluidic cooling array for microelectronics. Appl. Therm. Eng. 2018, 144, 71–80. [Google Scholar] [CrossRef]
  2. Bailey, C. Thermal management technologies for electronic packaging: Current capabilities and future challenges for modelling tools. In Proceedings of the 10th Electronics Packaging Technology Conference, Singapore, 9–12 December 2008; pp. 527–532. [Google Scholar]
  3. Krishnan, S.; Garimella, S.V.; Chrysler, G.M.; Mahajan, R.V. Towards a Thermal Moore’s Law. IEEE Trans. Adv. Packag. 2007, 30, 462–474. [Google Scholar] [CrossRef] [Green Version]
  4. Xie, G.; Sundén, B.; Wang, L.; Utriainen, E. Parametric study on heat transfer enhancement and pressure drop of an internal blade tip-wall with pin-fin arrays. Heat Mass Transf. 2011, 47, 45–57. [Google Scholar] [CrossRef]
  5. Guan, N.; Liu, Z.-G.; Zhang, C.-W. Numerical investigation on heat transfer of liquid flow at low Reynolds number in micro-cylinder-groups. Heat Mass Transf. 2012, 48, 1141–1153. [Google Scholar] [CrossRef]
  6. Sauciuc, L.; Chrysler, G.; Mahajan, R.; Szleper, M. Air-cooling Extension-Performance Limits for Processor cooling applicatioins. In Proceedings of the 19th Annual IEEE Semiconductor Thermal Measurement and Management Symposium, San Jose, CA, USA, 11–13 March 2003; pp. 74–81. [Google Scholar]
  7. Tuckerman, D.F.; Pease, R.F.W. High-Performance heat sinking for VLSI. IEEE Electron Device Lett. 1981, 2, 126–129. [Google Scholar] [CrossRef]
  8. Rajalingam, A.; Chakraborty, S. Effect of micro-structures in a microchannel heat sink—A comprehensive study. Int. J. Heat Mass Transf. 2020, 154, 119617. [Google Scholar]
  9. Zhang, B.-C.; Li, Q.-L.; Wang, Y.; Zhang, J.-Q.; Song, J.; Zhuang, F.-C. Experimental investigation of nitrogen flow boiling heat transfer in a single mini-channel. J. Zhejiang Univ. A 2020, 21, 147–166. [Google Scholar] [CrossRef]
  10. Wang, H.; Chen, Z.; Gao, J. Influence of geometric parameters on flow and heat transfer performance of micro-channel heat sinks. Appl. Therm. Eng. 2016, 107, 870–879. [Google Scholar] [CrossRef]
  11. Lin, L.; Zhao, J.; Lu, G.; Wang, X.-D.; Yan, W.-M. Heat transfer enhancement in microchannel heat sink by wavy channel with changing wavelength/amplitude. Int. J. Therm. Sci. 2017, 118, 423–434. [Google Scholar] [CrossRef]
  12. Park, M.-C.; Ma, S.-B.; Kim, K.-Y. Optimization of a Wavy Microchannel Heat Sink with Grooves. Processes 2021, 9, 373. [Google Scholar] [CrossRef]
  13. Hung, T.C.; Huang, Y.X.; Yan, W.M. Thermal performance of porous microchannel heat sink: Effects of enlarging channel outlet. Int. Commun. Heat Mass Tranf. 2013, 48, 86–92. [Google Scholar] [CrossRef]
  14. Ghani, I.; Sidik, N.A.C.; Mamat, R.; Najafi, G.; Ken, T.L.; Asako, Y.; Japar, W.M.A.A. Heat transfer enhancement in microchannel heat sink using hybrid technique of ribs and secondary channels. Int. J. Heat Mass Transf. 2017, 114, 640–655. [Google Scholar] [CrossRef]
  15. Chai, L.; Xia, G.D.; Wang, H.S. Parametric study on thermal and hydraulic characteristics of laminar flow in microchannel heat sink with fan-shaped ribs on sidewalls e Part 3: Performance evaluation. Int. J. Heat Mass Transf. 2016, 97, 1091–1101. [Google Scholar] [CrossRef]
  16. Chai, L.; Xia, G.D.; Wang, H.S. Parametric study on thermal and hydraulic characteristics of laminar flow in microchannel heat sink with fan-shaped ribs on sidewalls e Part 1: Heat transfer. Int. J. Heat Mass Transf. 2016, 97, 1069–1080. [Google Scholar] [CrossRef]
  17. Chai, L.; Xia, G.D.; Wang, H.S. Parametric study on thermal and hydraulic characteristics of laminar flow in microchannel heat sink with fan-shaped ribs on sidewalls e Part 2: Pressure drop. Int. J. Heat Mass Transf. 2016, 97, 1081–1090. [Google Scholar] [CrossRef]
  18. Chai, L.; Xia, G.D.; Wang, H.S. Numerical study of laminar flow and heat transfer in microchannel heat sink with offset ribs on sidewalls. Appl. Therm. Eng. 2016, 92, 32–41. [Google Scholar] [CrossRef]
  19. Chai, L.; Wang, L.; Bai, X. Thermohydraulic performance of microchannel heat sinks with triangular ribs on sidewalls e Part 1: Local fluid flow and heat transfer characteristics. Int. J. Heat Mass Transf. 2018, 127, 1124–1137. [Google Scholar] [CrossRef]
  20. Chai, L.; Wang, L.; Bai, X. Thermohydraulic performance of microchannel heat sinks with triangular ribs on sidewalls e Part 2: Average fluid flow and heat transfer characteristics. Int. J. Heat Mass Transf. 2019, 128, 634–648. [Google Scholar] [CrossRef]
  21. Xie, G.-F.; Zhao, L.; Dong, Y.-Y.; Li, Y.-G.; Zhang, S.-L.; Yang, C. Hydraulic and Thermal Performance of Microchannel Heat Sink Inserted with Pin Fins. Micromachines 2021, 12, 245. [Google Scholar] [CrossRef]
  22. Deng, D.; Pi, G.; Zhang, W.; Wang, P.; Fu, T. Numerical Study of Double-Layered Microchannel Heat Sinks with Different Cross-Sectional Shapes. Entropy 2019, 21, 16. [Google Scholar] [CrossRef] [Green Version]
  23. Nonino, C.; Savino, S. Temperature Uniformity in Cross-Flow Double-Layered Microchannel Heat Sinks. Fluids 2020, 5, 143. [Google Scholar] [CrossRef]
  24. Jing, D.; He, L. Thermal Characteristics of Staggered Double-Layer Microchannel Heat Sink. Entropy 2018, 20, 537. [Google Scholar] [CrossRef] [Green Version]
  25. Zhang, F.; Wu, B.; Du, B. Heat transfer optimization based on finned microchannel heat sink. Int. J. Therm. Sci. 2021, 172, 107357. [Google Scholar] [CrossRef]
  26. Sui, Y.; Lee, P.; Teo, C. An experimental study of flow friction and heat transfer in wavy microchannels with rectangular cross section. Int. J. Therm. Sci. 2011, 50, 2473–2482. [Google Scholar] [CrossRef]
  27. Khan, M.Z.U.; Uddin, E.; Akbar, B.; Akram, N.; Naqvi, A.A.; Sajid, M.; Ali, Z.; Younis, M.Y.; García Márquez, F.P. Investigation of Heat Transfer and Pressure Drop in Microchannel Heat Sink Using Al2O3 and ZrO2 Nanofluids. Nanomaterials 2020, 10, 1796. [Google Scholar] [CrossRef] [PubMed]
  28. Wang, S.-L.; Zhu, J.-F.; An, D.; Zhang, B.-X.; Chen, L.-Y.; Yang, Y.-R.; Zheng, S.-F.; Wang, X.-D. Heat transfer enhancement of symmetric and parallel wavy microchannel heat sinks with secondary branch design. Int. J. Therm. Sci. 2021, 171, 107229. [Google Scholar] [CrossRef]
  29. Ansari, D.; Husain, A.; Kim, K.-Y. Multiobjective Optimization of a Grooved Micro-Channel Heat Sink. IEEE Trans. Components Packag. Technol. 2010, 33, 767–776. [Google Scholar] [CrossRef]
  30. 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]
  31. Kestin, J.; Sokolov, M.; Wakeham, W.A. Viscosity of liquid water in the range −8 °C to 150 °C. J. Phys. Chem. Ref. Data 1978, 7, 941–948. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (a) Schematic diagram of the micro-channel and thermal chip; (b) dimensional schematic diagram of the MC-CR heat sink; (c) a basic unit in the MC-CR heat sink.
Figure 1. (a) Schematic diagram of the micro-channel and thermal chip; (b) dimensional schematic diagram of the MC-CR heat sink; (c) a basic unit in the MC-CR heat sink.
Micromachines 13 00132 g001
Figure 2. Assembly process.
Figure 2. Assembly process.
Micromachines 13 00132 g002
Figure 3. The unit models (a,c,e) and cross-sections (b,d,f) in the micro-channel heat sink: (a,b) MC-R; (c,d) MC-CR; (e,f) MC-HR.
Figure 3. The unit models (a,c,e) and cross-sections (b,d,f) in the micro-channel heat sink: (a,b) MC-R; (c,d) MC-CR; (e,f) MC-HR.
Micromachines 13 00132 g003aMicromachines 13 00132 g003b
Figure 4. Meshed model and residuals.
Figure 4. Meshed model and residuals.
Micromachines 13 00132 g004
Figure 5. Comparison between numerical and experimental data.
Figure 5. Comparison between numerical and experimental data.
Micromachines 13 00132 g005
Figure 6. The performance of MC-R, MC-CR and MC-HR when the velocity was 2 m/s. (a) The average temperature and the pressure drop; (b) the Nu and the equivalent thermal resistance; (c) the average convective heat transfer coefficient.
Figure 6. The performance of MC-R, MC-CR and MC-HR when the velocity was 2 m/s. (a) The average temperature and the pressure drop; (b) the Nu and the equivalent thermal resistance; (c) the average convective heat transfer coefficient.
Micromachines 13 00132 g006
Figure 7. The temperature contours on the horizontal (a,c,e) and longitudinal (b,d,f) cross-sections: (a,b) MC-R; (c,d) MC-CR; (e,f) MC-HR.
Figure 7. The temperature contours on the horizontal (a,c,e) and longitudinal (b,d,f) cross-sections: (a,b) MC-R; (c,d) MC-CR; (e,f) MC-HR.
Micromachines 13 00132 g007
Figure 8. The pressure drop contours on the horizontal (a,c,e) and longitudinal (b,d,f) cross-sections: (a,b) MC-R; (c,d) MC-CR; (e,f) MC-HR.
Figure 8. The pressure drop contours on the horizontal (a,c,e) and longitudinal (b,d,f) cross-sections: (a,b) MC-R; (c,d) MC-CR; (e,f) MC-HR.
Micromachines 13 00132 g008
Figure 9. The streamline and velocity contours. (a) MC-R; (b) MC-CR; (c) MC-HR.
Figure 9. The streamline and velocity contours. (a) MC-R; (b) MC-CR; (c) MC-HR.
Micromachines 13 00132 g009
Figure 10. The heat dissipation performance of MC-CR with different aspect ratios. (a) The average temperature; (b) the average convective heat transfer coefficient; (c) Nusselt number; (d) the equivalent thermal resistance; (e) the pressure drop.
Figure 10. The heat dissipation performance of MC-CR with different aspect ratios. (a) The average temperature; (b) the average convective heat transfer coefficient; (c) Nusselt number; (d) the equivalent thermal resistance; (e) the pressure drop.
Micromachines 13 00132 g010aMicromachines 13 00132 g010b
Figure 11. The average temperature and pressure drop of the heat source in different inclination angles of the cross-ribs.
Figure 11. The average temperature and pressure drop of the heat source in different inclination angles of the cross-ribs.
Micromachines 13 00132 g011
Figure 12. The average temperature and pressure drop of the heat source with different spacing of the cross-ribs.
Figure 12. The average temperature and pressure drop of the heat source with different spacing of the cross-ribs.
Micromachines 13 00132 g012
Figure 13. The unit models of the micro-channel heat sinks: (a) MC-SF; (b) MC-SR.
Figure 13. The unit models of the micro-channel heat sinks: (a) MC-SF; (b) MC-SR.
Micromachines 13 00132 g013
Figure 14. The performance of MC-R, MC-SF, MC-SR, and MC-CR when the pressure drop is 4.5 kPa. (a) The average temperature and the equivalent thermal resistance; (b) the Nu and average convective heat transfer coefficient.
Figure 14. The performance of MC-R, MC-SF, MC-SR, and MC-CR when the pressure drop is 4.5 kPa. (a) The average temperature and the equivalent thermal resistance; (b) the Nu and average convective heat transfer coefficient.
Micromachines 13 00132 g014
Table 1. The size parameters of the marked position on the three views in Figure 3.
Table 1. The size parameters of the marked position on the three views in Figure 3.
ParameterDimensionParameterDimension
L5 mmH1.5 mm
W1 mmHch0.75 mm
Wc0.5 mmHc1 mm
Wm0.25 mmHt0.1 mm
W_HR0.2 mmL_HR1.5 mm
α30°S0
Table 2. Five meshing methods and simulation results.
Table 2. Five meshing methods and simulation results.
MethodIIIIIIIVV
Minimum unit size (mm)0.10.080.060.040.02
Grid number11,50023,75155,029187,5001,504,500
Temperature (K)360.771356.880353.179351.128351.576
Table 3. Detailed information about the aspect ratio of MC-CR.
Table 3. Detailed information about the aspect ratio of MC-CR.
NoHcWcλ(Hc/Wc)
110.751.333
210.6251.6
310.52
41.250.52.5
510.3752.666
61.50.53
71.750.53.5
820.54
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, H.; Chen, C.; Zhou, Y.; Yang, C.; Song, G.; Hou, F.; Jiao, B.; Liu, R. Evaluation and Optimization of a Cross-Rib Micro-Channel Heat Sink. Micromachines 2022, 13, 132. https://doi.org/10.3390/mi13010132

AMA Style

Chen H, Chen C, Zhou Y, Yang C, Song G, Hou F, Jiao B, Liu R. Evaluation and Optimization of a Cross-Rib Micro-Channel Heat Sink. Micromachines. 2022; 13(1):132. https://doi.org/10.3390/mi13010132

Chicago/Turabian Style

Chen, Haiying, Chuan Chen, Yunyan Zhou, Chenglin Yang, Gang Song, Fengze Hou, Binbin Jiao, and Ruiwen Liu. 2022. "Evaluation and Optimization of a Cross-Rib Micro-Channel Heat Sink" Micromachines 13, no. 1: 132. https://doi.org/10.3390/mi13010132

APA Style

Chen, H., Chen, C., Zhou, Y., Yang, C., Song, G., Hou, F., Jiao, B., & Liu, R. (2022). Evaluation and Optimization of a Cross-Rib Micro-Channel Heat Sink. Micromachines, 13(1), 132. https://doi.org/10.3390/mi13010132

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