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

The 45 ̊ heat spreading angle is familiar among thermal designers. This angle has been used for thermal design of electronic devices, and provides a heat spreading area inside a board, e.g. printed circuit board, which is placed between a heat dissipating element and a relatively large heat sink. By using this angle, the heat transfer behavior can be estimated quickly without using high-performance computers. In addition, the rough design can be made easily by changing design parameters. This angle is effective in a practical situation; however, the discussion has not been made sufficiently on the applicability of the 45 ̊ heat spreading angle. In the present study, therefore, the extensive numerical investigation is conducted for the rational thermal design using the 45 ̊ heat spreading angle. The two-dimensional mathematical model of the board is considered; the center of the top is heated by a heat source while the bottom is entirely cooled by a heat sink. The temperature distribution is obtained by solving the heat conduction equation numerically with the boundary conditions. From the numerical results, the heat transfer behavior inside the board is shown and its relation with the design parameters is clarified. The heat transfer behavior inside the 45 ̊ heat spreading area is also evaluated. The applicability is moreover discussed on the thermal resistance of the board obtained by the 45 ̊ heat spreading angle. It is confirmed that the 45 ̊ heat spreading angle is applicable when the Biot number is large, and then the equations are proposed to calculate the Biot number index to use the 45 ̊ angle. Furthermore, the validity of the 45 ̊ heat spreading angle is also confirmed when the isothermal boundary condition is used at the cooled section of the board.


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 [1], the early paper about this spreading angle goes back to the publication by Balents et al. in 1969.Balents et al. [2] 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. [3] 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 [4] 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 [5] 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 [6] described briefly that if a 20% error could be accepted, the 45 degree rule was always relevant for small heat sources.Lasance [7] 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 [7] and Ha and Graham [8] 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.

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 z b .The center (radius: r h ) of the top is heated by a heat source while the bottom (radius: r c ) is entirely cooled by a heat sink.The 45˚ heat spreading angle is shown in Figure 2.
The governing equation is given by and the following boundary conditions are applied for the heated section (0 ≤ r ≤ r h , z = 0) and the cooled section (0 ≤ r ≤ r c , z = z b ), respectively:  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 T f the cooling fluid temperature.The following adiabatic boundary condition is applied at the surfaces except for the heated and the cooled section: where n is the coordinate normal to the boundary surface.
The governing equation is discretized by the control volume method [9] 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. [10].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 z b = 4.0 mm, q = 1.0 W/cm 2 and T f = 20˚C.

Dimensionless Parameters
The dimensionless temperature, θ, and the dimensionless coordinates, R, Z, N, are defined as follows: Accordingly, the dimensionless governing equation and the dimensionless boundary conditions for the heated section (0 ≤ R ≤ R h , Z = 0), the cooled section (0 ≤ R ≤ R c , Z = 1) and the other adiabatic section are expressed as follows: ( ) where R h , R c and Bi are the following dimensionless radii and the Biot number: , , It should be noted that the heat transfer behavior inside the present model is characterized by R h , R c , and Bi.

Temperature Distribution
The representative numerical results are shown in Figures 3 and 4. The effect of λ is shown at α = 100 W/(m 2 ⋅K) and α = 5000 W/(m 2 •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/(m 2 ⋅K) and α = 5000 W/(m 2 ⋅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, T , is obtained by the following equation: The value of T is also shown in Figure 4.When α = 100 W/(m 2 ⋅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 T f = 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/(m 2 ⋅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 24.0 C T =  .In this case, the heat is spread out over the board.It is also confirmed that the temperature distri- bution becomes flatter and approaches T with the increase in λ.
When α = 5000 W/(m 2 ⋅K), it is observed that the temperature inside the board is entirely lower than that in the case of α = 100 W/(m 2 ⋅K).Furthermore, the temperature at the cooled section is very close to the uniform temperature at 20.1 C T =  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 R h , R c and Bi.It should be noted that the scales of the vertical axis for R h = 1, R h = 2, R h = 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 R h and R c .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 R h and R c respectively, it is expected that the temperature entirely increases with the increase in R h or the decrease in R c .However, due to the concentration of heat flow near the heated section, it is found that the temperature distribution is hardly affected by R c , namely, the board width when Bi is large (Bi = 10).

45˚ Heat Spreading Characteristics
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 6 shows the relation between η and λ changing α and r c as parameters.As reference, the values of η = 0.80 and η = 0.90 are indicated by 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: 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 R h , R c 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 R h = 4 and R c = 5, because  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 R h , R c and Bi.It is observed that irrespective of R and R c , η 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 R c when Bi is large.This is due to the concentration of heat flow near the heated section.From ) Equations ( 14) and ( 15) are also shown in Figure 8.It is confirmed that the 45˚ heat spreading angle is appli-

Thermal Resistance
When the 45˚ heat spreading angle is used, the thermal resistance of the board, R 45 , is simply obtained as ( ) The derivation of Equation ( 16) is shown in Appendix.From the numerical results, on the other hand, the thermal resistance of the board, R b , is also obtained by subtracting the convective thermal resistance at the cooled section from the total thermal resistance.R b is calculated by, where T h,ave is the average temperature at the heated section.By using the ratio of R 45 to R b , the comparison between these thermal resistances is shown in Figure 9 changing all dimensionless parameters of R h , R c and Bi.Irrespective of the numerical conditions, it is found that the value of R 45 /R b is less than 1, implying that R 45 is smaller than R b .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 R 45 /R b with Bi and R c is similar to that in η; while in the range of R 45 /R b > 0.8, the effect of R h on R 45 /R b is smaller than the case of η.Therefore, Bi = 10 is the index to use Equation ( 16) when 20% error is accepted.

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.

Figure 7 , 8 .
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 Since η is not affected by R c for large Bi as mentioned above, Bi| η=0.80 and Bi| η=0.90 depend on only R h .Moreover, the relations between Bi| η=0.80 , Bi| η=0.90 and R h are expressed as follows:

Figure 8 .:
Figure 8. Biot number index to use 45˚ heat spreading angle.cablewhen 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.

Figure 9 .
Figure 9.Thermal resistance of board: comparison between 45˚ heat spreading approach and numerical results.