﻿ Two-Dimensional Numerical Investigation on Applicability of 45<sup>。</sup>Heat Spreading Angle

Journal of Electronics Cooling and Thermal Control
Vol.4 No.1(2014), Article ID:43494,11 pages DOI:10.4236/jectc.2014.41001

Two-Dimensional Numerical Investigation on Applicability of 45˚ Heat Spreading Angle

Yasushi Koito, Shoryu Okamoto, Toshio Tomimura

Department of Advanced Mechanical Systems, Kumamoto University, Kumamoto, Japan

Email: koito@gpo.kumamoto-u.ac.jp  Received 9 December 2013; revised 9 January 2014; accepted 16 January 2014

Abstract

Keywords

45˚ Heat Spreading Angle; Electronic Board; Electronics Cooling; Thermal Design; Heat Conduction 1. Introduction

A heat spreading angle is one of the concepts to simplify the heat transfer calculation and has been used for thermal design of electronic devices. This angle is applied for the board, e.g. printed circuit board, placed between a heat dissipating element and a relatively large heat sink, and then assumes a heat spreading area inside the board. While high-performance CFD simulators provide detailed numerical information, such a simplified calculation is also effective in a practical situation because the estimation can be made quickly without using high-performance computers. In addition, the rough design can be made easily by changing design parameters such as dimensions, physical properties and operating conditions, etc.

According to Vermeersch and Mey  , the early paper about this spreading angle goes back to the publication by Balents et al. in 1969. Balents et al.  used the constant value of 45˚ as the heat spreading angle and conducted an approximate heat transfer analysis concerning a heat dissipating element mounted directly on a ceramic substrate which was bonded to a metal baseplate. After a while, the 45˚ heat spreading model of Balents et al. was used by Cook et al.  in 1976, and thereafter many attempts have been made to modify the value or extend the theory of the heat spreading angle. Nevertheless, the above-mentioned 45˚ heat spreading angle is still familiar among thermal designers.

The applicability of the 45˚ heat spreading angle was described by several researchers. Masana  showed the comparison between the 45˚ heat spreading approach and the Fourier solution for a substrate (thickness: w) on which a square (dimensions: l ´ l) heat dissipating element was mounted. It was described that the common assumption of 45˚ for the spreading angle got close to the Fourier solution for large w/l. Guenin  also compared the 45˚ heat spreading analysis to the finite element analysis, and then concluded that the 45˚ heat spreading approximation provided the reasonable accuracy for boards in which the size of the heat source area was greater than the board thickness. Malhammer  described briefly that if a 20% error could be accepted, the 45 degree rule was always relevant for small heat sources. Lasance  reported that the 45˚ heat spreading angle was acceptable in a single layer board for h/k < 1 and d > 2 mm, where h was the heat transfer coefficient at the cooled section, k and d were the thermal conductivity and the thickness of the board, respectively. Furthermore, concerning the limitation, Lasance  and Ha and Graham  mentioned that the 45˚ heat spreading approach could not be used for multi-layer boards.

The above descriptions are guides in conducting the heat transfer calculation using the 45˚ heat spreading angle. However, they seem to be incoherent; accordingly, the comprehensive investigation is needed for the rational thermal design using the 45˚ heat spreading angle. In the present study, therefore, the extensive numerical investigation is conducted on the applicability of the 45˚ heat spreading angle. For simplicity, the two-dimensional mathematical model is considered and the calculation is conducted in a cylindrical coordinate system. Based on the numerical results, the validity of the 45˚ heat spreading angle and its relation with design parameters are discussed.

2. Analytical Method

2.1. Mathematical Modeling

Figure 1 shows the analytical system. The heat transfer characteristics inside the board are analyzed in a cylindrical (r − z) coordinate system. This board is essentially desk-shaped, and has the thickness of zb. The center (radius: rh) of the top is heated by a heat source while the bottom (radius: rc) is entirely cooled by a heat sink. The 45˚ heat spreading angle is shown in Figure 2.

The governing equation is given by , (1)

and the following boundary conditions are applied for the heated section (0 ≤ r ≤ rh, z = 0) and the cooled section (0 ≤ r ≤ rc, z = zb), respectively: (2) (3) Figure 1. Analytical system. Figure 2. 45˚ heat spreading angle.

where λ is the thermal conductivity of the board, q the heat flux from the heat source, α the heat transfer coefficient at the cooled section and Tf the cooling fluid temperature. The following adiabatic boundary condition is applied at the surfaces except for the heated and the cooled section: , (4)

where n is the coordinate normal to the boundary surface.

The governing equation is discretized by the control volume method  and solved numerically with the boundary conditions. The analytical solution of Equations (1)-(4) having infinite series and Bessel functions was also shown by Lee et al.  . However, because the analytical solution is complicated and time-consuming, the numerical procedure is adopted here to obtain the results quickly. The numerical calculation is conducted by changing the parameters as shown in Table 1 at zb = 4.0 mm, q = 1.0 W/cm2 and Tf = 20˚C.

2.2. Dimensionless Parameters

The dimensionless temperature, θ, and the dimensionless coordinates, R, Z, N, are defined as follows: . (5)

Accordingly, the dimensionless governing equation and the dimensionless boundary conditions for the heated section (0 ≤ R ≤ Rh, Z = 0), the cooled section (0 ≤ R ≤ Rc, Z = 1) and the other adiabatic section are expressed as follows: , (6) (7) (8) (9)

where Rh, Rc and Bi are the following dimensionless radii and the Biot number: . (10)

It should be noted that the heat transfer behavior inside the present model is characterized by Rh, Rc, and Bi.

Table 1. Numerical conditions.

3. Results and Discussion

3.1. Temperature Distribution

The representative numerical results are shown in Figures 3 and 4. The effect of λ is shown at α = 100 W/(m2×K) and α = 5000 W/(m2∙K). Figure 3 is the temperature contour inside the board, and Figure 4 the temperature distribution at the cooled section (the top surface of the board). It should be noted that in Figure 4, the scales of the vertical axis for α = 100 W/(m2×K) and α = 5000 W/(m2×K) are different. In addition, if the conductive thermal resistance inside the board is neglected, the temperature of the board becomes uniform and then the uniform temperature, , is obtained by the following equation: . (11)

The value of is also shown in Figure 4. When α = 100 W/(m2×K) and λ = 0.4 W/(m×K), it is observed that the contour lines inside the board are crowded near the heated section. Accordingly, the temperature at the cooled section changes conspicuously near the heated section and then approaches the cooling temperature at Tf = 20˚C. In this case, because the conductive thermal resistance inside the board is relatively large, the heat flow is not spread out but concentrated near the heated section. When α = 100 W/(m2×K) and λ = 400 W/(m×K), on the other hand, it is observed that except for the region near the heated section, the contour lines are almost perpendicular to the top/bottom surfaces and then the temperature at the cooled section is very close to the uniform temperature at . In this case, the heat is spread out over the board. It is also confirmed that the temperature distribution becomes flatter and approaches with the increase in λ.

When α = 5000 W/(m2×K), it is observed that the temperature inside the board is entirely lower than that in the case of α = 100 W/(m2×K). Furthermore, the temperature at the cooled section is very close to the uniform temperature at even in the case of λ = 0.4 W/(m×K). This is attributed to the decrease in the convective thermal resistance at the cooled section, which enhances the heat transfer in z direction.

The temperature distribution at the cooled section is moreover shown in Figure 5 in the dimensionless form. The numerical results are shown by changing all dimensionless parameters of Rh, Rc and Bi. It should be noted that the scales of the vertical axis for Rh = 1, Rh = 2, Rh = 4 are different. In addition, the value of Bi is also indicated in Figure 3 for reference. Based on Bi, the heat transfer behavior inside the board is characterized as follows. When Bi is very small, namely Bi = 0.001, because the conductive thermal resistance inside the board is very small compared to the convective thermal resistance at the cooled section, the temperature is almost uniform irrespective of Rh and Rc. When Bi = 0.1, the temperature distribution is conspicuous due to the increase in the conductive thermal resistance inside the board. Furthermore, when Bi = 10, because of the decrease in the convective thermal resistance at the cooled section, the temperature is entirely lower than that in the case of Bi = 0.1.

Because the heated and the cooled area are prescribed by Rh and Rc respectively, it is expected that the temperature entirely increases with the increase in Rh or the decrease in Rc. However, due to the concentration of heat flow near the heated section, it is found that the temperature distribution is hardly affected by Rc, namely, the board width when Bi is large (Bi = 10).

In order to evaluate the 45˚ heat spreading characteristics, the index, η, is defined here by the ratio of the heat transfer rate inside the 45˚ heat spreading area to the total heat transfer rate through the board. η is calculated by, Figure 3. Temperature distribution inside board (rh = 8 mm, rc = 40 mm). Figure 4. Temperature distribution at cooled section (rh = 8 mm, rc = 40 mm). Figure 5. Dimensionless temperature distribution at cooled section. . (12)

Figure 6 shows the relation between η and λ changing α and rc as parameters. As reference, the values of η = 0.80 and η = 0.90 are indicated by lines in this figure. As expected, it is observed that η increases with the decrease in λ because of the concentration of heat flow near the heated section; η increases with α because of the enhancement of heat transfer in z direction. It is also confirmed that the high value of η such as η = 0.80 and η = 0.90 is obtained for small λ and large α, which corresponds to relatively large Bi.

The numerical calculation is moreover conducted by using the following isothermal boundary conditions: . (13)

This boundary condition is applied at the cooled section and then the numerical result of η is also shown in Figure 6. It is found that irrespective of the numerical conditions, the high value of η (η » 0.91) is obtained under Equation (13), implying that the 45˚ heat spreading angle is applicable when the isothermal boundary condition is used at the cooled section.

By changing all dimensionless parameters of Rh, Rc and Bi, η is rearranged as shown in Figure 7 in the dimensionless form. It should be noted that η = 1 is obtained irrespective of Bi when Rh = 4 and Rc = 5, because Figure 6. 45˚ heat spreading characteristics (rh = 8 mm).

the cooled section is completely included within the 45˚ heat spreading area. This trivial case is not shown in this figure. The values of η = 0.80 and η = 0.90 are also indicated by lines in this figure. As shown, η is well characterized by Rh, Rc and Bi. It is observed that irrespective of Rh and Rc, η increases with Bi and then relatively large increase in η is obtained between Bi = 0.01 and Bi = 10. Furthermore, it is also found that the value of η is not affected by Rc when Bi is large. This is due to the concentration of heat flow near the heated section.

From Figure 7, the values of Bi when η = 0.80 and η = 0.90 are obtained, and then these values denoted by Bi|η=0.80 and Bi|η=0.90 respectively are shown in Figure 8. Since η is not affected by Rc for large Bi as mentioned above, Bi|η=0.80 and Bi|η=0.90 depend on only Rh. Moreover, the relations between Bi|η=0.80, Bi|η=0.90 and Rh are expressed as follows: (14) (15)

Equations (14) and (15) are also shown in Figure 8. It is confirmed that the 45˚ heat spreading angle is appli- Figure 7. 45˚ heat spreading characteristics in dimensionless form. Figure 8. Biot number index to use 45˚ heat spreading angle.

cable when Bi is relatively large as expected, and then the specific numerical values of Bi are obtained by Equations (14) and (15) when 20% and 10% errors are accepted, respectively.

3.3. Thermal Resistance

When the 45˚ heat spreading angle is used, the thermal resistance of the board, R45, is simply obtained as . (16)

The derivation of Equation (16) is shown in Appendix. From the numerical results, on the other hand, the thermal resistance of the board, Rb, is also obtained by subtracting the convective thermal resistance at the cooled section from the total thermal resistance. Rb is calculated by, , (17)

where Th,ave is the average temperature at the heated section.

By using the ratio of R45 to Rb, the comparison between these thermal resistances is shown in Figure 9 changing all dimensionless parameters of Rh, Rc and Bi. Irrespective of the numerical conditions, it is found that the value of R45/Rb is less than 1, implying that R45 is smaller than Rb. This is owing to the fact that the effect of the conductive thermal resistance inside the board in r direction is neglected in Equation (16). It is also observed that the change in R45/Rb with Bi and Rc is similar to that in η; while in the range of R45/Rb > 0.8, the effect of Rh on R45/Rb is smaller than the case of η. Therefore, Bi = 10 is the index to use Equation (16) when 20% error is accepted. Figure 9. Thermal resistance of board: comparison between 45˚ heat spreading approach and numerical results.

4. Conclusions

The numerical analysis is conducted on the applicability of the 45˚ heat spreading angle inside the board. In the present calculation range, the main findings can be summarized as follows:

1) The heat transfer behavior inside the board is well characterized by the three dimensionless parameters: the two dimensionless widths of the heated and the cooled section, and the Biot number.

2) However, when the Biot number is large, the heat transfer behavior is hardly affected by the dimensionless width of the cooled section due to the concentration of heat flow near the heated section.

3) The validity of the 45˚ heat spreading angle is confirmed when the isothermal boundary condition is used at the cooled section.

4) The 45˚ heat spreading angle is applicable when the Biot number is large. The equations are proposed to calculate the Biot number index to use the 45˚ heat spreading angle.

5) The thermal resistance of the board obtained by the 45˚ heat spreading angle can be used with 20% error when the Biot number is 10.

References

1. Vermeersch, B. and Mey, G.D. (2008) A Fixed-Angle Dynamic Heat Spreading Model for (an) Isotropic Rear-Cooled Substrates. ASME Journal of Heat Transfer, 130, 121301. http://dx.doi.org/10.1115/1.2976557
2. Balents, L., Gold, R.D., Kaiser, A.W. and Peterson, W.R. (1969) Design Considerations for Power Hybrid Circuits. Proceedings of the International Hybrid Microelectronics Symposium, Dallas, 29-30 October 1969, 323-344.
3. Cook, K.B., Kerns, D.V., Nagle, H.T., Slagh, T.D. and Ruwe, V.W. (1976) Computer-Aided Thermal Analysis of a Hybrid Multistage Active Bandpass Filter/Amplifier. IEEE Transactions on Parts, Hybrids, and Packaging, 12, 344- 350. http://dx.doi.org/10.1109/TPHP.1976.1135151
4. Masana, F.N. (1996) A Closed Form Solution of Junction to Substrate Thermal Resistance in Semiconductor Chips. IEEE Transactions on Components, Packaging, and Manufacturing Technology, Part A, 19, 539-545. http://dx.doi.org/10.1109/95.554935
5. Guenin, B. (2003) The 45˚ Heat Spreading Angle—An Urban Legend? Electronics Cooling. http://www.electronics-cooling.com/2003/11/the-45-heat-spreading-angle-an-urban-legend/
7. Lasance, C.J.M. (2010) How to Estimate Heat Spreading Effects in Practice. ASME Journal of Electronic Packaging, 132, 031004. http://dx.doi.org/10.1115/1.4001856
8. Ha, M. and Graham, S. (2012) Development of a Thermal Resistance Model for Chip-On-Board Packaging of High Power LED Arrays. Microelectronics Reliability, 52, 836-844. http://dx.doi.org/10.1016/j.microrel.2012.02.005
9. Patankar, S.V. (1980) Numerical Heat Transfer and Fluid Flow. Hemisphere Pub. Corp., Washington DC.
10. Lee, S., Song, S., Au, V. and Moran, K.P. (1995) Constriction/Spreading Resistance Model for Electronics Packaging. Proceedings of the 4th ASME/JSME Thermal Engineering Joint Conference, Hawaii, 19-24 March 1995, 4, 199-206.

Appendix

The thermal resistance, R, is defined by , (18)

where Q(z) is the heat transfer rate inside the board at the position of z. When the 45˚ heat spreading angle is used, Q(z) is expressed as (19)

Substituting Equation (19) into Equation (18) and then integrating from z = 0 to z = zb, R45 is obtained as shown in Equation (16).

Nomenclature

Bi :        Biot number (-);

N  :        dimensionless coordinate normal to surface (-);

n   :        coordinate normal to surface (m);

Q  :        heat transfer rate (W);

q   :        heat flux (W/cm2, W/m2);

R   :        thermal resistance (K/W); dimensionless radial coordinate (-); dimensionless radius (-);

T   :        temperature (˚C);

Tf    :        cooling temperature (˚C); :        uniform temperature (˚C);

Z   :        dimensionless axial coordinate (-);

z    :        axial coordinate (mm, m); thickness (mm, m).

Greek Symbols

a   :        heat transfer coefficient (W/(m2×K));

h   :        index defined by Equation (12) (-);

l   :        thermal conductivity (W/(m×K));

q   :        dimensionless temperature (-).

Subscripts

45 :        45˚ heat spreading angle;

ave:       average;

b   :        board;

c   :        cooled section;

h   :        heated section.