Research Papers

Assessment of Joule Heating and Temperature Distribution on Printed Circuit Boards Via Electrothermal Simulations

[+] Author and Article Information
M. Baris Dogruoz

Cisco Systems, Inc.,
San Jose, CA 95134
e-mail: dogruozb2@asme.org

Contributed by the Electronic and Photonic Packaging Division of ASME for publication in the JOURNAL OF ELECTRONIC PACKAGING. Manuscript received February 19, 2016; final manuscript received March 18, 2016; published online April 19, 2016. Assoc. Editor: Mehmet Arik.

J. Electron. Packag 138(2), 021004 (Apr 19, 2016) (8 pages) Paper No: EP-16-1035; doi: 10.1115/1.4033109 History: Received February 19, 2016; Revised March 18, 2016

A printed circuit board (PCB) comprises a solid piece of dielectric material with embedded layers of current carrying metal traces and vias. Geometric features of these metal traces and vias in modern PCBs are highly nonuniform and complicated such that the card level or system level numerical simulations by using the actual trace and via geometries are computationally expensive. The present study investigates the effects of Joule heating in current carrying traces on the temperature distribution of PCBs by conducting one-way and two-way direct current (DC) electric and computational fluid dynamics (CFD) simulations. DC electric field simulations are performed to determine the power map of trace layers which are modeled as planar heat generating sources by using the temperature-dependent electrical conductivity of the metal trace. The power distribution varies with the implemented size and power thresholds. Thermal conductivity map of the PCB is determined by using the electronic computer-aided design (ECAD) images of the individual layers. By using these planar source and thermal conductivity maps, CFD simulations are conducted to determine the resulting temperature distribution on the board. A methodology is developed and applied to a sample, complex PCB, and the generated results are compared with those of the previous studies and conventional models. The computational data show that the temperature distributions over the PCB and its mounted components experience large variations based on the implemented thermal conductivity mapping and the Joule heating modeling technique.

Copyright © 2016 by ASME
Your Session has timed out. Please sign back in to continue.


Sadiku, M. , 2000, Numerical Techniques in Electromagnetics, CRC Press, Boca Raton.
Jin, J. , 2002, The Finite Element Method in Electromagnetics, 2nd ed., Wiley, New York.
Horowitz, M. , and Dutton, R. W. , 1983, “ Resistance Extraction From Mask Layout Data,” IEEE Trans. Comput. Aided Des. Integr. Circuit Syst., 2(3), pp. 145–150. [CrossRef]
Vithayathil, A. , Hu, X. , and White, J. , 2003, “ Substrate Resistance Extraction Using a Multi-Domain Surface Integral Formulation,” International Conference on Simulation of Semiconductor Processes and Devices (SISPAD 2003), Boston, MA, Sept. 3–5, SISPAD, pp. 323–326.
Wang, X. , Yu, W. , and Wang, Z. , 2006, “ Efficient Direct Boundary Element Method for Resistance Extraction of Substrate With Arbitrary Doping Profile,” IEEE Trans. Comput. Aided Des. Integr. Circuit Syst., 25(12), pp. 3035–3042. [CrossRef]
Yang, B. , and Murata, H. , 2007, “ A Finite Element Domain Decomposition Coupled Resistance Extraction Method With Virtual Terminal Insertion,” 50th Midwest Symposium on Circuit and Systems (MWSCAS 2007), Montreal, QC, Canada, Aug. 5–8, pp. 1425–1428.
Rajagopalan, S. , and Batterywala, S. , 2007, “ A 3-Dimensional FEM Based Resistance Extraction,” 20th International Conference on VLSI Design(VSLID), Bangalore, India, Jan. 6–10, pp. 565–570.
Oppeneer, M. , Sumant, P. , and Cangellaris, A. , 2008, “ Robust Iterative Finite Element Solver for Multi-Terminal Power Distribution Network Resistance Extraction,” IEEE Electrical Performance of Electronic Packaging (EPEP), San Jose, CA, Oct. 27–29, pp. 181–184.
Shakeri, K. , and Meindl, J. D. , 2005, “ Compact Physical IR Drop Models for Chip/Package Co-Design of Gigascale Integration (GSI),” IEEE Trans. Electron. Devices, 52(6), pp. 1087–1096. [CrossRef]
Kose, S. , and Freidman, E. G. , 2012, “ Efficient Algorithms for Fast IR Drop Analysis Exploiting Locality,” Integr. VLSI J., 45(2), pp. 149–161. [CrossRef]
Ozisik, M. N. , 1989, Boundary Value Problems of Heat Conduction, Dover Publications, Mineola, NY.
Ellison, G. N. , 1996, “ Thermal Analysis of Circuit Boards and Microelectronic Components Using an Analytic Solution to the Heat Conduction Equation,” 12th Annual IEEE Semiconductor Thermal Measurement and Management Symposium (SEMI-THERM XII), Austin, TX, Mar. 5–7, pp. 144–150.
Ellison, G. N. , 1992, “ Extensions of the Closed Form Method for Substrate Thermal Analyzers to Include Thermal Resistances From Source-to-Substrate and Source-to-Ambient,” IEEE Trans. Compon., Hybrids, Manuf. Technol., 15(5), pp. 658–666. [CrossRef]
Lemczyk, T. F. , Mack, B. L. , Culham, J. R. , and Yovanovich, M. M. , 1992, “ PCB Trace Thermal Analysis and Effective Conductivity,” ASME J. Electron. Packag., 114(4), pp. 413–419. [CrossRef]
Eveloy, V. , Lohan, J. , and Rodgers, P. , 2000, “ A Benchmark Study of Computational Fluid Dynamics Predictive Accuracy for Component-Printed Circuit Board Heat Transfer,” IEEE Trans. Compon. Packag. Technol., 23(3), pp. 568–577. [CrossRef]
Dogruoz, M. B. , and Nagulapally, M. K. , 2009, “ Effects of Trace Layers and Joule Heating on the Temperature Distribution of Printed Circuit Boards: A Computational Study,” ASME J. Therm. Sci. Eng. Appl., 1(2), p. 022003. [CrossRef]
Shankaran, G. V. , Dogruoz, M. B. , and deAraujo, D. , 2010, “ Orthotropic Thermal Conductivity and Joule Heating Effects on the Temperature Distribution of Printed Circuit Boards,” 12th Intersociety Conference on Thermal and Thermomechanical Phenomena in Electronics Systems (ITherm), Las Vegas, NV, June 2–5.
Dogruoz, M. B. , and Shankaran, G. V. , 2012, “ Spatial Variation of Temperature on Printed Circuit Boards: Effects of Anisotropic Thermal Conductivity and Joule Heating,” IEEE Trans. Compon., Packag., Manuf. Technol., 2(10), pp. 1649–1658. [CrossRef]
Stefani, G. G. , Goel, N. S. , and Jenks, D. B. , 1993, “ An Efficient Numerical Technique for Thermal Characterization of Printed Wiring Boards,” ASME J. Electron. Packag., 115(4), pp. 366–372. [CrossRef]
Agarwal, R. K. , Dasgupta, A. , Pecht, M. , and Barker, D. , 1991, “ Prediction of PWB/PCB Thermal Conductivity,” Int. J. Hybrid Microelectron., 14(3), pp. 83–95.
Graebner, J. E. , and Azar, K. , 1997, “ Thermal Conductivity Measurements in Printed Wiring Boards,” ASME J. Heat Transfer, 119(3), pp. 401–405. [CrossRef]
Sarvar, F. , Poole, N. J. , and Witting, P. A. , 1990, “ PCB Glass-Fibre Laminates: Thermal Conductivity Measurements and Their Effect on Simulation,” J. Electron. Mater., 19(12), pp. 1345–1350. [CrossRef]
Coppola, L. , Cottet, D. , and Wildner, F. , 2008, “ Investigation on Current Density Limits in Power Printed Circuit Boards,” 23rd IEEE Applied Power Electronics Conference and Exposition (APEC 2008), Austin, TX, Feb. 24–28, pp. 205–210.
Shankaran, G. V. , and Singh, R. K. , 2006, “ Selection of Appropriate Thermal Model For Printed Circuit Boards in CFD Analysis,” 10th Intersociety Conference on Thermal and Thermomechanical Phenomena in Electronics Systems (ITherm '06), San Diego, CA, May 30–June 2, pp. 234–242.
Galloway, J. , and Shidore, S. , 2004, “ Implementing Compact Thermal Models Under Non-Symmetric Trace Routing Conditions,” 20th Annual IEEE Semiconductor Thermal Measurement and Management Symposium (STHERM), San Jose, CA, Mar. 9–11, pp. 255–261.
Adam, J. , 2004, “ New Correlations Between Electrical Current and Temperature Rise in PCB Traces,” 20th Annual IEEE Semiconductor Thermal Measurement and Management Symposium (STHERM), San Jose, CA, Mar. 9–11, pp. 292–299.
IPC, 2003, “ Generic Standard on Printed Board Design,” Association Connecting Electronics Industries, Bannockburn, IL, Standard No. IPC-2221A.
IPC, 2009, “ Standard for Determining Current Carrying Capacity in Printed Circuit Board Design,” Association Connecting Electronics Industries, Bannockburn, IL, Standard No. IPC-2152.
Xie, J. , and Swaminathan, M. , 2011, “ Electrical-Thermal Co-Simulation of 3D Integrated Systems With Micro-Fluidic Cooling and Joule Heating Effects,” IEEE Trans. Compon., Packag., Manuf. Technol., 1(2), pp. 234–246. [CrossRef]
ANSYS, 2008, “ SIWave 4.0 User's Guide,” ANSYS, Canonsburg, PA.
ANSYS, 2015, “ Icepak 16.1 User's Guide,” ANSYS, Canonsburg, PA.
Cupertino, F. , and Ettorre, S. , 2013, “ Experimental Evaluation of Current Carrying Capacity of Printed Circuit Stator Coils,” 39th Annual Conference of the IEEE Industrial Electronics Society (IECON 2013), Vienna, Austria, Nov. 10–13, pp. 2810–2815.


Grahic Jump Location
Fig. 1

A complex PCB design

Grahic Jump Location
Fig. 2

Flowchart showing the computational methodology

Grahic Jump Location
Fig. 3

CFD model (H: high (top) side and L: low (bottom) side of the PCB)

Grahic Jump Location
Fig. 4

Model validation—symbols: experiments [32] and solid line: present computations

Grahic Jump Location
Fig. 5

Current density vectors (j) in the top trace layer: (a) PCB top surface and (b) zoom-in to the via array on the “NVDD” trace

Grahic Jump Location
Fig. 6

Two-way simulations (case 6) with 5 mW/0.127 mm (power/size) filters. Top—Joule heating power map (left) and temperature distribution (right) for the top layer. Bottom—Joule heating power map (left) and temperature distribution (right) for the bottom layer. Lines (plots on the left)—Bold: top surface components, medium: bottom surface components, light: sources in the layers. There are a total of 638 sources on the top layer, 204 sources on the bottom layer, and 1072 sources on the PCB.

Grahic Jump Location
Fig. 7

Thermal conductivity distribution in normal direction, κt,: top layer (top) and bottom layer (bottom)

Grahic Jump Location
Fig. 8

Detailed versus lumped thermal conductivity—individual component temperature values of cases 1 and 2, with qJ = 0

Grahic Jump Location
Fig. 9

Effects of Joule heating on individual component temperature values—cases 2 and 4

Grahic Jump Location
Fig. 10

Effects of Joule heating modeling on individual component temperature values—cases 4–6

Grahic Jump Location
Fig. 11

Lumped versus fully detailed model—individual component temperature values of cases 3 and 6

Grahic Jump Location
Fig. 12

Distribution of DR values for cases 2, 4, and 5




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In