About this Journal Submit a Manuscript Table of Contents
Active and Passive Electronic Components
Volume 2011 (2011), Article ID 823654, 8 pages
Research Article

Application of Thermal Network Model to Transient Thermal Analysis of Power Electronic Package Substrate

1Department of Mechanics System Engineering, Toyama Prefectural University, Imizu 939-0398, Japan
2COSEL CO. LTD., 1-6-43 Kamiakae-machi, Toyama 930-0816, Japan

Received 12 April 2011; Revised 4 August 2011; Accepted 6 September 2011

Academic Editor: P. Markondeya Raj

Copyright © 2011 Masaru Ishizuka et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


In recent years, there is a growing demand to have smaller and lighter electronic circuits which have greater complexity, multifunctionality, and reliability. High-density multichip packaging technology has been used in order to meet these requirements. The higher the density scale is, the larger the power dissipation per unit area becomes. Therefore, in the designing process, it has become very important to carry out the thermal analysis. However, the heat transport model in multichip modules is very complex, and its treatment is tedious and time consuming. This paper describes an application of the thermal network method to the transient thermal analysis of multichip modules and proposes a simple model for the thermal analysis of multichip modules as a preliminary thermal design tool. On the basis of the result of transient thermal analysis, the validity of the thermal network method and the simple thermal analysis model is confirmed.

1. Introduction

Laptop personal computers and cellular phones have become smaller and thinner in recent years, while their functions have become more sophisticated. As a result, the heat dissipation density and dissipation power per volume of these devices have significantly increased. Cooling design is the key to the prevention of heat problems. At the same time, there are growing demands for the acceleration of product development and reduction of costs.

Computational fluid dynamics (CFDs) codes [1, 2] have proved their high potential as a tool of thermal design of electronic equipment. However, as the product development cycle is shortened, the CFD-based thermal design needs a new format that allows the packaging designer fast and versatile searches for better design options. The most serious factor that slows the CFD-based design is the geometric complexity created by packing various components in a tight space of the system box.

Recently, in order to reduce the CFD workload, several approaches for the thermal design of electronic equipment have been reported. Nakayama et al. [35] proposed a methodology coined build-up approach (BUA) and have applied it to the board- and box-level thermal analyses. Minichiello and Belady [6] proposed a thermal design methodology and applied it to the design of a multi-processor enterprise server, the RP8400. Their methodology combines the well-known analytical and experimental thermal design tools, including heat transfer correlations, flow network modeling, and CFD techniques and experimental measurements.

In a broader perspective, the issue of our interest is concerned with model reduction methodologies, and on this subject, several reports have been presented. Assumption of velocity profiles in computational elements reduces the nonlinear momentum equation to a linear equation, and such linearized models have been proposed by Moosmann et al. [7], D’Amore et al. [8], and Rudnyi et al. [9]. Molina and Clemente [10] reported on an automated technique for reducing actual systems to lumped-parameter thermal models.

Ishizuka and Hayama [11] have also proposed a thermal analysis scheme, in which CFD analysis and thermal network method are coupled. In this scheme, the essential factors affecting thermal fluid phenomena in electronic equipment are first clarified using CFD analysis, then, using the CFD results obtained under suitable boundary conditions, main heat flow passages are identified. At the end of the model-reduction process, a thermal network model is developed. The thermal network model facilitates the design parameter survey for choosing a set of optimum parameters without resorting to time-consuming CFD analysis. Ishizuka et al. have proposed some simplified models such as heat diffusion model in [12], thermal convection model in [13], phase changing process model in [14, 15], and 3D model in [15] to simulate various thermal behaviors in electronic devices. From those results, it is found that if we would understand the heat phenomenon in a system, some temperature values for the system could be predictable by using a simplified model for a thermal network method.

The present study focuses on the derivation of a simpler network model to predict the transient temperature distribution of the package substrate covered by the cap, where the package consists of high-power transistor chips, other IC chips, resistor chips, and capacitor chips. The analysis was carried out to simulate the results of the transient temperature rise on the surface of a substrate, taking into consideration natural convection and radiation, in the heat transfer mode for a cap-attached circuit substrate in addition to the thermal conduction. The alumina substrate has a circuit composed of various chips centering around four power transistor chips (hereinafter called power transistors). The heat transfer modes for packages are very complicated. The validity of the proposed thermal network model is demonstrated on the basis of a comparison between calculated values and measured values.

2. Experiment

2.1. Package Structure

Various kinds of chips such as IC, resistor, and capacitor, centered around four power transistors (square of 4.5 mm in size, 0.3 mm in thickness, each having power dissipation of 8 W), are mounted and wire-bonded to form an electronic circuit on a 1.6 mm-thick and 6 cm-square substrate formed from 92% alumina ceramic. The vicinity of the power transistors is illustrated in Figure 1. The power transistors, a heat source in this circuit, are arranged at 4 mm intervals.

Figure 1: (a) Power-transistor layout. (b) Power transistors (photo).

The power dissipation value contributed from other chips may be neglected. Interconnections are provided within the substrate, forming a multilayer printed circuit board. This chip group is covered with a 5 cm-square, 2.8 mm-high and 0.25 mm-thick kovar cap, in which nitrogen gas (N2 gas) is sealed. Figure 2 shows a cross-sectional view of the substrate. The power transistor chips were mounted on the molybdenum heat sink blocks, which were aligned on the copper heat sink blocks with electroconductive epoxy. The other IC chips, resistor chips and capacitor chips were also mounted on the substrate top conductor metal die pads with electroconductive epoxy. They were then wire-bonded with 30 μm diameter gold wire.

Figure 2: Cross-sectional view of the substrate.
2.2. Temperature Measurement

The module to be measured was held horizontally 10 cm above the surface of a desk and supported at its edges by three thin glass epoxy sticks (3 mm in diameter) in a clean room with free airflow. Temperature was controlled to 298 K. A point equally distant from the four power transistors was considered as the center (radius ), and two points were located on the back of the substrate () and one on the front ( mm). A copper-constantan thermocouple (0.1 mm in diameter) for use in making the measurements was bonded by good conductive epoxy material to each of the measuring points. As a measuring procedure, temperature rise values after the application of electric power from a DC power source to the power transistors were recorded by means of an infrared camera. The experiment was stopped when temperature at measuring point on the back of the substrate exceeded 410 K.

3. Analysis

3.1. Heat Transfer Model for Multichip Modules

The paths through which package heat transfer takes place was considered as follows. The substrate bottom area and substrate top area, with the exception of the kovar cap area, include heat transport by thermal radiation and natural convection. Inside the kovar cap, heat transport by heat conduction of the nitrogen gas and thermal radiation was considered. At the kovar cap top surface area, heat transport occurs by thermal radiation and natural convection. Furthermore, inside the substrate and the cap, heat transport occurs only by thermal conduction. These aspects are illustrated in Figure 3.

Figure 3: Heat transfer form for package.

Figure 3 indicates the paths by which heat transfer from the module takes place, where the following major heat transfer factors were considered:(1)thermal conduction inside the substrate,(2)natural convection and thermal radiation from the top and bottom surfaces of the substrate to the atmosphere,(3)thermal conduction through N2 gas and thermal radiation from the substrate surfaces to the cap,(4)thermal conduction from the substrate to the cap through the welded area between the substrate and the cap (not shown in Figure 3).

3.2. Thermal Network Model

The existence of the cap is considered to have an influence on the temperature distribution; therefore, the effect of the cap should be taken into consideration. However, since considering the cap results in a complicated three-dimensional analytical model, there is no advantage in adopting the thermal network method. Therefore, the analytical model is formed on the basis of the following assumptions to reduce the model’s dimensions from three to one.(1)The temperature of the center of the cap surface is the same as that of the substrate area where the cap is welded.(2)The 4 power transistors are regarded as one nodal point.(3)The heat dissipated by the power transistors flows out in the radial direction.

These assumptions will make the analytical model one dimensional (axial symmetry).

The substrate was divided, as shown in Figure 4, centering on the power transistor (). It is assumed that there is no temperature difference between individual divided partitions and that the partitions are numbered , respectively, from the center. In this study, was taken as 6. Representing the heat resistance and the thermal capacitance for each partition by and, respectively, the thermal network model for the heat transfer is formed, where heat flux is generated by the power transistors, as shown in Figure 5. The total of the heat flux was assumed to flow out toward the region 1 in the substrate. Since the area of the region 0 including the power transistor chips were smaller than those of the other regions, the heat transfer for the central region 0 was neglected. In Figure 5, in the case of the first partition where the temperature is represented by , heat flow flows out toward external region 2 through heat resistance due to the thermal conduction. Some of the heat is stored in the thermal capacitance , or dissipated directly or indirectly through cap, out of the substrate by natural convection and thermal radiation through thermal resistance values such as and . It is assumed that value is infinite. Note that although the actual module forms a regular square and the power transistors are arranged quadrilaterally, these may be represented approximately by circles, having diameters equal to the length of one side of the respective rectangle, as shown in Figure 4. On the basis of the previously mentioned assumptions, the center of the cap was linked with the welded portion of the substrate without any thermal resistance.

Figure 4: Substrate partition.
Figure 5: Thermal network model (Model 1).

In addition, the cap thermal capacity was neglected, because it is infinitesimally small compared to that of substrate. In addition, it is assumed, as described previously, that cap temperature is equal to substrate temperature in partition 5.

3.3. Evaluation of Thermal Resistance

The thermal conduction partition model in the substrate is illustrated in Figure 6. It was assumed that a circular heat source, whose radius is  cm, exists at the center of a circular multilayer ceramic substrate, whose thickness is cm, and the heat diffuses in the radial direction. If the temperature at a point on radius is and the temperature at a point of radius is , the heat flow rate along the radial direction is given by where , is the thermal conductivity of the substrate.

Figure 6: Thermal conduction partition model in the substrate.

Therefore, the thermal resistance, in region covered by radius and radius is given by

The thermal resistance for nitrogen gas inside the kovar cap is given by where is thermal conductivity, is the substrate-to-cap distance, and represents the thermal conduction area.

If there exists a temperature difference () between the air and the substrate, the heat flow rate which passes through the area of the substrate is given by

In general, the heat transfer coefficient, is related to the heat transfer phenomenon. If the flow is laminar and the heat transfer occurs through natural convection, the natural convection heat coefficient in the case of two-dimensional flat plate is given considering that Gr·Pr value (Grashof and Prandtl number) in this study was the order of 105 as follows where is the temperature rise with reference to ambient temperature ;  K is selected as a reference temperature and the value of is determined by the orientation location of the flat plate, shown in Figure 7. Facing up it is equal to 0.54, while facing down it is 0.27. is the characteristic length which is determined by the flow path of the cooling air on the hot surface of the material. In this divided doughnut-shaped plate, having inner and outer radius of and , respectively, the assumed characteristic length is given by the ratio of area/perimeter of the doughnut-region as follows: Then, thermal resistance is expressed as follows, where represents the thermal radiation area: In the exchange of thermal radiation energy between two materials which are not black bodies, the general thermal radiation equation is given by where is Stefan-Boltzmann’s constant which is  W/m2 K4, is the thermal radiation area, and and are absolute temperature, respectively. Furthermore, is the shape factor which is defined as where and are emissivities for the materials. If one side is the atmosphere, becomes infinity; therefore, becomes .

Figure 7: Value in natural heat transfer equation (5).

From (8), thermal resistance is expressed as follows:

3.4. Calculation of Thermal Capacitance

The calculation of thermal capacitance of each partition, for example, is for the th partition, is expressed as follows, by defining the substrate density as , specific heat at constant pressure as , and substrate thickness as :

3.5. Equation Formulation

The thermal network method is based on the similarity between the diffusion equation for thermal engineering and that for electric engineering, and it is found that the analogous quantities are voltage-temperature and current-heat flow.

 Therefore, the transient analog method is formed on the basis of a number of lumped tee sections, as shown in Figure 8, with subsequent employment of the node analysis method. The following equations, which can be written as examples, are based on heat flow conservation for a nodal point and for heat paths based on thermal resistance definition where represents the heat flow value through thermal resistance and represents the temperature value.

Figure 8: Heat flow at tee section.
3.6. Calculation Procedure

The above equations are nonlinear in nature, due to (7) and (10), and iterative calculations are required. However, to simplify the method of calculating parameters, temperature , required for calculating (7) and (10), is approximated by a value calculated with instead of , assuming that temperature change is very small during the period when is very small; that is, . To solve these equations, the Runge-Kutta secondary approximate method is used. The calculation flow chart is given below.(1)Set the applied power .(2)Give the initial temperature value for the initial time .(3)Calculate the thermal resistance values from (7) and (10) using the above temperature value and thermal capacitance .(4)Solve the linearized equation system to yield the temperature for the new time level using the Runge-Kutta method.(5)Advance the time level by .(6)The newly predicted temperatures are assigned to the initial temperature.(7)The procedure is repeated from step (5) for the subsequent time level until .

3.7. Numerical Example

As the calculation condition, the power transistor region, the shapes of the cap, and the substrate are approximated as a circle of diameter 1.5 cm, 5 cm, and 6 cm, respectively, as shown as Figure 4. It is also assumed that the substrate is partitioned into six sections (number of partitions ), having radii of  cm (central partition),  cm,  cm,  cm,  cm,  cm, and  cm, respectively, where represents the radius for each circle. Substrate characteristics are as follows. Thickness:  mm, thermal conductivity:  W/(mK), density:  kg/m3 , and specific heat: J/kg K. Other factors are emissivity: , nitrogen thermal conductivity: W/mK, the cap height:  mm, and the cap thermal conductivity: W/mK. However, in this study, the property values of the cap were not decided using the assumptions. The calculations were carried out using a time increment  ms, where the conductivity of the substrate was substituted by the conductivity of alumina which was measured by the laser-flash method, because the volume of alumina occupies 98% or more of the total substrate volume.

4. Results

4.1. Comparison of Calculated Values and Measured Values

Figure 9 compares the measured and predicted values for the case of, , , and . is the temperature of the partition 0 as shown in Figure 4, near the center of the substrate bottom surface, is the temperature of the partition 3, near the midpoint between the edge of the substrate bottom surface and the center of the substrate bottom surface, is the temperature of the partition 5, near the edge of the substrate bottom surface, and is the temperature near the center of the kovar cap. For and , temperatures at two different points in the same region were averaged. Temperature differences between the two temperatures were within 3%. Agreement between measured and predicted values is very satisfactory. In order to prevent the damage of chip, the power supply switch was turned off in the case of temperature reaching to a value of 410 K.

Figure 9: Comparison between calculated values and measured values.

Figure 9 also shows result using the results of calculated temperature on the basis of another model in which thermal diffusion is caused only by the substrate thermal conduction (Model 2 as shown in Figure 10). It can be observed here that the temperature predicted by the model 2 is approximately 10% higher as compared to that obtained using the present model (Model 1). The trends of temperature rise for is similar to that in the case of the Model 2. The difference in temperatures other than was within the range of 5% to 8%. Here, since the heat sink effect due to the thermal capacity of the substrate is much large at only initial stage of the temperature rise, the results by model 2 considering only conduction was compared once by model 1. Of course, time passing, the difference becomes much larger. Of course, time passing, the difference between them becomes much larger.

Figure 10: Thermal network model (Model 2).

5. Discussion

In forming the heat transfer model for the present experiment condition, it was assumed that the thermal conduction inside the substrate plays a dominant role in contributing to the heat transfer of the package. The validity of this assumption was verified by comparing the results obtained by considering only thermal conduction and those obtained by considering natural convection and thermal radiation as well as thermal conduction, as shown in Figure 9. The temperature difference between the results obtained by using the heat transfer model and those obtained by using the model which considers only thermal conduction is at the most 10%. That is the reason that in the transient condition, the substrate capacitance works as a main heat sink, and in the steady state condition, the atmosphere does. However, the results obtained by using the proposed heat transfer model coincide closely with the results obtained by using the model which considers only thermal conduction. On the basis of the recent trends that the power dissipation density values for future electronic elements will become higher or they will be used in a steady-state condition, this 10% difference cannot be neglected. For the current power dissipation density (where power dissipation density value in this substrate is 32 watts for 36 cm2 area), this proposed model is considered to be useful as a simulation model for thermal design because of its simple and convenient form. It should also be noted that when a forced cooling method is applied instead of natural convection, which was adopted in this experiment, the effect of the proposed model should be higher if (5) were replaced by the forced convection heat transfer coefficient. The difference between results obtained by using the proposed model and by the other model where the natural convection and thermal radiation were neglected would be larger.

6. Conclusion

A thermal network model was applied to simulate a transient temperature rise in small-size multichip modules. The heat transfer process by natural convection and thermal radiation was also adopted in this model in addition to the thermal diffusion by thermal conduction. The calculated results using this model gave a very reasonable estimation of the measured substrate temperature values. Since substrate thermal designing becomes more important and great design accuracy is required, as chip power dissipation density becomes larger, this model is considered sufficiently practical as an engineering model for simulating a transient temperature rise in the substrate, because it is structurally simple and easy to use.


:Area, m2
:Thermal capacitance, J/K
:Specific heat, J/(kg K)
:Constant in the natural heat transfer equation
:Substrate thickness, m
:Shape factor
:Cap height, m
:Reference length, m
:Heat flow rate, W
:Thermal resistance, K/W
:Temperature, K
:Temperature rise, K.

Greek Symbols

:Heat transfer coefficient, W/m2 K
:Thermal conductivity, W/m K
:Density, kg/m3
:Stefan-Boltzmann constant = 5.67 × 10−8, W/(m2 K4)
:Time, s.


cd:Thermal conduction
cv:Natural convection
ra:Thermal radiation


  1. http://www.cradle.co.jp/en/index.htm.
  2. http://www.flomerics.com/.
  3. W. Nakayama, R. Matsuki, Y. Hacho, and K. Yajima, “A new role of CFD simulation in thermal design of compact electronic equipment: application of the build-up approach to thermal analysis of a benchmark model,” Journal of Electronic Packaging, Transactions of the ASME, vol. 126, no. 4, pp. 440–448, 2004. View at Publisher · View at Google Scholar · View at Scopus
  4. W. Nakayama, T. Nakajima, H. Koike, and R. Matsuki, “Heat conduction in printed circuit boards—part I; overview and the case of a JEDEC test board,” in Proceedings of IPACK, ASME InterPACK, Vancouver, BC, Canada, July 2007.
  5. W. Nakayama, T. Nakajima, H. Koike, and R. Matsuki, “Heat conduction in printed circuit boards—part I; part ii; small PCBs connected to large thermal mass at their edge,” in Proceedings of IPACK, ASME InterPACK, Vancouver, BC, Canada, July 2007.
  6. A. Minichiello and C. Belady, “Thermal design methodology for electronic systems,” in Proceedings of the Inter Society Conference on Thermal Phenomena, San Diego Calif, USA, June 2002.
  7. C. Moosmann, E. B. Rudnyi, A. Greiner, and J. G. Korvink, “Model order reduction for linear convective thermal flow,” in Proceedings of the 10th International Workshop on Thermal Investigations of ICs and Systems, Sophia Antipolis, France, October 2004.
  8. D. D'Amore, L. Codecasa, and P. Maffezzoni, “An Arnoldi based thermal network reduction method for electro-thermal analysis,” IEEE Transactions on Components and Packaging Technologies, vol. 26, no. 1, pp. 186–192, 2003. View at Publisher · View at Google Scholar
  9. E. B. Rudnyi, T. Bechtold, and J. G. Korvink, “Automatic generation of compact electrothermal models for semiconductor devices,” IEICE Transactions on Electronics, vol. 86C, pp. 459–465, 2003.
  10. M. Molina and C. Clemente, “Thermal model reduction: algorithms and validation techniques,” SAE Technical Papers 2006-01-2112, 2006.
  11. M. Ishizuka and S. Hayama, “Application of a semi-empirical approach to the thermal design of electronic equipment,” IMechE Journal of Power and Energy, vol. 214, no. 5, pp. 513–522, 2000. View at Scopus
  12. M. Ishizuka, S. Nakagawa, T. Hachiga, and H. Ishida, “Application of a thermal network method to non-steady thermal analysis of electronic equipment,” Far East Journal of Mathematical Sciences, vol. 28, no. 3, pp. 667–680, 2009.
  13. M. Ishizuka, S. Nakagawa, K. Koizumi, and Y. Kitamura, “Practical study of application of thermal network method to thermal design of electronic equipment (an example for thermal design of a compact self-ballasted fluorescent lamp),” Far East Journal of Applied Mathematics, vol. 26, no. 3, pp. 59–608, 2007.
  14. M. Ishizuka, “Operation time control of a high density packaging using a low melting point alloy,” IEEE Transactions on Components, Packaging Technologies, vol. 27, no. 2, pp. 239–243, 2004.
  15. S. Takakuwa, M. Ishizuka, and S. Nakagawa, “Application of a thermal network method to thermal (16)analysis of electronic devices using phase change materials,” in Proceedings of the IPACK, ASME InterPACK, San Francisco, Calif, USA, July 2009.