• Thermal Management of Electronic Assemblies and Enclosures

Applications of CFD tools in PCB, Board-level and Chassis Level Thermal Simulations

Electronics Cooling


Table of Contents: THERMAL CHARACTERIZATION OF IC PACKAGES | MOSFET | The Gerber file | ICEPAK: Basic Solver Setting | Heat Transfers in Fins | Thermo-Electric Cooling | Thermal Simulations in Data Centres | EMN-EMP Files | FPGA | Efficiency and Power Losses in MOSFETS | Joint Electron Device Engineering Council (JEDEC) standards | Heat Generation Rates in Electronic Parts | IDF and IDX: ECAD - EMN/EMP | Thermal Resistances


chip Cooling

There are vast range of applications of thermal simulations in electronics industry dealing with Printed Circuit Board (PCB), IC Chips, Bipolar junction transistor (BJT), Diodes, MOSFET, FPGA, Power Transistors, Bipolar and Darlington Transistors, IGBT, Capacitors, Power Supply, Power Switches, Voltage Regulators, Controllers, Inductors, Solid State Relays, Bulk Acoustic Wave (BAW) Oscillator, Ultra-Low Jiter Oscillator, Flash PROM, Flash Mux, DDR4 SDRAM, DDR Clock Buffer, Digital Step Attenuator, Digital Isolator, Power Sequencer, Phase Locked Loop, Transceiver chip, Control relays, Op-Amp, Analog Multiplier, Silicon-Controlled Rectifiers (SCRs) ... In addition to board and chassis level simulation, thermal management of data centres include specialized cooling methods like heat sinks, TEC (Thermo-Electric Cooling), heat pipes, PCM (Phase Change Material), heat spreaders are also presented in this page. A thermal network model for MOSFET-on-PCB is shown below. The applications include Synchronous Step Down Switcher, Synchronous buck converter, Inverting Buffers / Drivers, DC/DC Controllers, Inverting Schmitt Triggers...

chip Thermal Network

The thermal network is a reasonably accurance and extermely fast approach to evaluate casing temperature of heat generating devices. TNsolver is a free open-source code written in GNU Octave to solver such thermal networks. For application of Machine Leaning in predicting thermal aspects of such devices, refer the article "Approximating the Steady-State Temperature of 3D Electronic Systems with Convolutional Neural Networks" by Monika Stipsitz and Hèlios Sanchis-Alepuz.

Equivalence of Radiative Heat Transfer Rate at Low Temperature and Small Temperature Differences

Convective heat transfer rate: qCNV'' = h x [T - TREF].

Radiative heat flux rate: qRAD'' = ε x σ x [TW14 - TW24] = ε x σ x [TW1 - TW2] x [TW1 + TW2] x [TW12 + TW22]

Comparing the two equations: hRAD ≡ ε x σ x [TW1 + TW2] x [TW12 + TW22]

  For ε = 1.0, TW1 = 100 [°C] and TW2 = 40 [°C], hRAD = 9.23 [W/m2.K]

  For ε = 1.0, TW1 = 80 [°C] and TW2 = 40 [°C], hRAD = 8.41 [W/m2.K]

  For ε = 1.0, TW1 = 60 [°C] and TW2 = 50 [°C], hRAD = 8.01 [W/m2.K]

As you can see, the convective heat transfer coefficient in natural convection with air ranges between 5 to 15 [W/m2.K]. Hence, the contribution of radiation in natural convection cases (such as electronic cooling) is approximatly equal to the convective heat transfer rate. In other words, convection and radiation contribute almost equal in natural convection heat transfer cases dealing with low temperature differences.

In general, the heat dissipation through the bodies of devices become difficult when the volume density of heat generation rate exceeds 108 [W/m3]. Such devices are often connected to the PCB with thermal vias having thickness ~ 0.5 [mm].


Types of Simulations

  • PCB or Package Level: Detail modeling of copper traces, thermal vias and dielectric material inside the PCB panel.
  • Board Level: PCB modeled as simplified orthotropic thermal conductivity.
  • Chassis Level: Similar to board-level simulation, bigger in computational domain.
  • Rack Level such as in Data Centres: Most of the heat sources are modelled as lumped blocks.

Electronics Cooling


Four ways PCB are modeled in ICEPAK

"PCB is developed to mechanically support and electrically connect electronic components through conductive tracks and pads. Components are typically soldered onto the PCB to be mechanically and electrically connected to it. Many active (for example, operational amplifiers and batteries) and passive components (such as inductors, resistors, and capacitors) are mounted on the PCBs."
  1. Hollow PCB: useful only in concept stage
  2. Compact PCB: all layers lumped in a single isotropic material (orthotropic thermal conductivities in all 3 directions)
  3. Detailed PCB: Each layer of PCB modeled with 'UNIFORM' thermal conductivity of that layer, calculated as volume weighted average of copper and substrate (FR4)
  4. ECAD PCB: ECAD files with detail traces and vias.
    • Actual geometry of the traces and vias are not meshed along with the rest of the CFD model
    • Each of the metal and dielectric layer is modeled separately
    • Effect of the traces and vias is modeled in each of the CFD mesh cell by computing orthotropic thermal conductivity for the imported traces and vias.

Inputs Required for Thermal Simulations:

The inputs from PCB Design (Hardware Team) to the thermal simulation team should be provided in the format described below.

S. No.Reference DesignatorQtyVendor Part NumberOutput Power [W]Efficiency [%]Heat Generation Rate [W]Mounting type on PCB
01U11LTM1235.075%1.25Vias and pin with airgap
02U25MJ12010.280%0.04Full Contact
.....................TIM
.....................Soldered Pins with airgap
Power Budget of the PCB Board: Energy balance helps in checking for correctness of the inputs - "worst-case estimates" sometimes have too much margin (extra heat generation rate) to keep the temperatures low.
Total input power to the PCB Board [W]:20.0
Total ouput power to the PCB Board [W]:16.0
Heat generation in the copper traces inside PCB Board [W]:0.50
Voltage regulators take power from the overall energy and make it into usable energy for circuit blocks. E.g. CPU cores are driven off low-voltage maybe 1.2 [V] or lower. The IOs on CPU may be driven from a 3.3 [V] regulator. Other circuit blocks on PCB may need 5 [V], and remainings something else. Power budget needs to be arrived by supply current and supply voltage or output current and output voltage for each device on the PCB.

For linear regulators, total current is sum of all currents, regardless of voltage. For switchers, efficiency is required to calculate required current at given voltage. Power distribution network (PDN) ensures stable power supply to all electronic components. PDN consists of traces, vias, planes, Voltage Regulators, and decoupling capacitors. PDN distributes power from the primary source throughout the PCB board to ensure voltage supply to various components. A good source to understand basics of PCB design is "tessolve.com/blogs/power-distribution-network-in-pcb-design-ensuring-stable-power-delivery".


Heat Generation and Heat Transfer Path

The units of power and heat, both are [W] and it creates confusion while interpreting the information among power-electronics or hardware designers and mechanical thermal engineers. The heat generation rate of a device is difference between power input and power output and it is important to note that power delivered by a device to a load is not power dissipated in the device as heat. When no specific efficiency curves are available in a data sheet for the application, an assumption of the efficiency is to be considered to calculate the input power. Typically, this value can range between 70% to 90%. in case CDR (Direct Current Resistance) is known from the data sheet, it is easy to calculate heat generation rate in the device. For many devices, the power dissipation consists of two basic components - the unloaded power dissipation inherent to the device and the load power dissipation which is a function of the device loading. For example, the loading of a logic device can significantly effect the power dissipation. Most of the logic loads are capacitive, leading to more of dynamic power dissipation.

The effect of radiation heat transfer is very important in natural convection, as it can contribute to approximately 25% of the total heat dissipation. Unless the components are facing hotter surface nearby such as enclosure, the radiative heat transfer must be accounted for.

Integrated circuit chip

Many devices are mounted to the PCB though the connector pins and the case is not in full contact with the PCB. This mounting arrangement must be taken into account to model the heat transfer path from case to the PCB. One option is to compare the contact area of the pins and PCB with that of the case and PCB and apply appropriate contact resistance.

Construction Features on IC Packages

Die: It designates the piece of semiconductor on which all the active circuits lie. Dies are made of Silicon (110 – 150 W/m.K) or Gallium Arsenide (GaAs, 45 – 60 W/m.K) is used in special applications such as microwaves. The circuitry is present within a thin layer on one side only, known as active surface. The concept of "junction temperature" applies to the top surface of the die. The heat flux in the die and die-attach is very high, the junction temperature is very sensitive to the thermal conductivity and the thickness of the die and die-attach.

The Die is often attached to the substrate or the die pad by an adhesive known as the "die attach" which is made of an epoxy based compound having thickness 0.025 - 0.050 [mm] and thermal conductivity in the range 1 - 2 [W/m.K]. The semiconductor material of a die has a conductivity approaching that of a metal. If the active layer is assumed to generate constant power per unit area (an assumption that may not always be valid), the die will be practically isothermal. In practice, applications exist for which the heat flux varies significantly across the die, in which case a temperature gradient may exist on its active surface. For most packages, the thermal resistance offered by the die is small in comparison with that offered by the rest of the package. Referecnce: FloTHERM Pack User's Guide

Due to low thickness of the "die attach", it causes very little heat spreading. At the same time, it offers significant thermal resistance in out-of- the-plane direction due to its poor thermal conductivity value. Si is very widely used because of its cheap price, easiness in processing and fairly high thermal conductivity. SiC is an alternative material [though very costly] that has high thermal conductivity of 370 [W/m-K].

Die Pad or Die Flag: The Die is placed in insulated boxed often plastic packages, on a thin metal plate (made of copper and larger than the die) known as the die flag or die pad. It helps both in the manufacturing and thermal (heat dissipation) function. The metallic die flag acts as an effective heat spreader due to very high thermal conductivity of copper which can reduce the thermal resistance of a package by up to 15%.

Componenets of a Die

Bond wires, made of Gold or Aluminum having diameter of the order of 0.025 [mm], are characteristics of wire-bonded packages. The number of bond wires in an IC package are of same order as the number of external leads/pins. Excerpts from "FloTHERM Pack User's Guide": In most ceramic packages, a negligible portion of the heat from the Die flows to the substrate through the bond wires. However, in plastic packages this may not be the case. Bond wires play a significant role in the heat transfer within peripheral leaded packages such as the PQFP. In area-array plastic packages such as the PBGA, bond wires can be important, especially for a 2-layer substrate.

Componenets of an Integrated Circuit - IC Chip

Encapsulant: Also known as 'Overmold', it is an epoxy based compound with a thermal conductivity 0.6 - 0.7 [W/m.K].

Construction Detail of an Integrated Circuit - IC Chip

Thermal Resistance an Integrated Circuit - IC Chip

LED - Integrated Circuit - IC Chip


Reference: FloTHERM PACK User Guide

CBGA = Ceramic Ball Grid Array, PBGA = Plastic Ball Grid Array, TBGA: Tape Ball Grid Array, CPGA: Ceramic Pin Grid Array

Flip Chip CBGA

Bond Wire TBGA

Thermal conductivities: Alumina - 21 [W/m-K], Kovar - 17 [W/m-K], Typical encapsulant and elastomer material - 0.4 ~ 0.6 [W/m-K], Polyimide: 0.2 [W/m-K], 37Pb/63Sn solder: 50 [W/m-K]

Ceramic Pin Grid Arrays


THERMAL CHARACTERIZATION OF IC PACKAGES
The construction of an IC (Integrated Circuit) chip is not homogeneous in 3D space. The actual heat generating segment is made-up of silicon (with thermal conductivity of 96 [W/m-K]) which will be surrounded by an Encapusulant before being attached to a Printed Circuit Board. Hence, a thermal resistance of an IC package is specified by its manufacturers which is the measure of ability of the IC package to transfer heat generated by the IC (die) to the printed circuit board or the ambient.

two Resistor Model

  1. The junction temperature is represented with a single temperature node.
  2. The thermal flow path is represented by a thermal resistance network.
  3. Junction Node represents the thermal source of the chip.
  4. Case Node represents the upper surface of the package.
  5. Board Node represents the circuit board temperature at a position of 1 [mm] from the edge of the device.
  6. Two thermal resistors are connected between "the junction and the case" and between "the junction and the board".

Two-resistor model is a relatively simple model achieved by dividing a package vertically at the junction and is good for single function devices. However, this model has the least precise as compared to multi-resistor network and detailed thermal models. This model (like any thermal resistor model) does not support "Transient Analyses". The methods of thermal resistance measurement are described for θJB in JEDEC Standard JESD51-8 and for θJC in JESD51-14.

In general, thermal resistance Θ or RTH = ΔT / P where P is power dissipation from the chip package in [W] and ΔT = TJ - TA.

Sample from datasheet of a DD-DC Step-down Power Supply from Linear Technology:

DC-DC Step-Down Module Regulator

TJ is the junction temperature and TA is the ambient temperature. In other words, the thermal resistance of an IC package is defined as the amount of heat generated or a rise in temperature when 1 [W] of power is dissipated in the IC. E.g. if an IC package generate 0.50 [W] with Junction-to-Ambient thermal resistance of 100 [K/W], the temperature of junction will rise by 0.5 [W] × 100 [K/W] = 50 [K] over the ambient. The construction of IC packages decides the thermal resistances: few types are Dual Inline Packages (DIP), Ceramic through-hole package, Plastic through-hole package, Plastic leaded chip carrier (PLCC), Quad Flat Package (QFP), Plastic dual construction surface mounte package, Thin small outline package (TSOP)...
  • ΘJA: also designated as RΘJA, it is the thermal resistance from Junction to Ambient, measured as [K/W]. Ambient is a generic term to thermal 'ground'. This value depends on the package design, PCB, and heat transfer modes (airflow, radiation). As per "Semiconductor and IC Package Thermal Metrics" by Texas Instruments, it is a measure of the thermal performance of an IC package mounted on a specific test coupon.

    The standard test conditions are specified by Joint Electron Device Engineering Council (JEDEC) such as JESD15-3: Two-Resistor Compact Thermal Model Guideline, JESD51: Methodology for the Thermal Measurement of Component Packages (Single Semiconductor Device), JESD51-9 "Test Boards for Area Array Surface Mount Package Thermal Measurements" and JESD51-12: Guidelines for Reporting and Using Package Thermal Information.

    As described in webinar "Understanding Datasheet Thermal Parameters and IC Junction Temperatures" by Monolithic Power Systems - ΘJA is valid only for its defined PCB and it is not a constant which can be used on all PCBs. ΘJA allows comparison of different packages on a common PCB.

    As per datasheet for LTM8023: "θJA is the natural convection junction-to-ambient air thermal resistance measured in a one cubic foot sealed enclosure. This environment is sometimes referred to as still air although natural convection causes the air to move. This value is determined with the part mounted to a JESD 51-9 defined test board, which does not reflect an actual application or viable operating condition."

    Another use of RJA or ΘJA, is to calculate parameter called derating factor when the power dissipation values are unspecified in the datasheet. Using the formula TJ = TA + P × ΘJA, the permissible heat generation rate can be estimate for various ambient temperatures say between 25 [°C] to 75 [°C].

  • ΘJC: also designated as RΘJC, it is the thermal resistance from Junction to Case where case is a specified point on the outside surface of the package. It depends on the thermal conductivity of package materials (the lead frame, mold compound, die attach adhesive) and on the specific package design (thickness of the die body and die attach, exposed pads and internal thermal vias). ΘJC considers only the resistance of heat flow paths to the surface of the package and hence represents only the conductive heat transfer path thermal resistance. EIA/JESD51-1 defines RΘJC as "the thermal resistance from the operating portion of a semiconductor device to outside surface of the package (case) closest to the chip mounting area when that same surface is properly heat sunk so as to minimize temperature variation across that surface".

    As described in datasheet for LTM8023: "θJCbottom is the junction-to-board thermal resistance with all of the component power dissipation flowing through the bottom of the package. In the typical μModule regulator, the bulk of the heat flows out the bottom of the package, but there is always heat flow out into the ambient environment. As a result, this thermal resistance value may be useful for comparing packages but the test conditions don’t generally match the user’s application. θJCtop is determined with nearly all of the component power dissipation flowing through the top of the package. As the electrical connections of the typical μModule regulator are on the bottom of the package, it is rare for an application to operate such that most of the heat flows from the junction to the top of the part. As in the case of θJCbottom, this value may be useful for comparing packages but the test conditions don’t generally match the user’s application."

  • ΘCA: it is the thermal resistance from Case to Ambient. Thus: ΘJA = ΘJC + ΘCA
  • ΨJB: junction-to-board thermal-characterization parameter and ΨJT: junction-to-top thermal-characterization parameters are the "characterization parameters" that measure temperature change between the junction temperature and the temperatures of the board and top of the package respectively. These are not thermal resistances in true sense though they are useful to estimate the junction temperature when the temperature on top of the package or the board and the power dissipation are known.

    As described in datasheet for LTM8023: "θJB is the junction-to-board thermal resistance where almost all of the heat flows through the bottom of the μModule regulator and into the board, and is really the sum of the θJCbottom and the thermal resistance of the bottom of the part through the solder joints and through a portion of the board. The board temperature is measured a specified distance from the package, using a two sided, two layer board. This board is described in JESD 51-9."

The Thermal Resistance is calculated in a condition where nearly all of the component power dissipation flows through either the top or the bottom of the package. In calculating the Thermal Characterization Parameter, power is the total power dissipation in the chip which can flow out of the chip through any thermal path, not just from the top or the bottom of the package.

Junction to case resistance

The heat transfer path and various thermal resistance as described above are explained in the schematic from Aavid Thermalloy Inc which is shown below.

heat Transfer Path MOSFETS

From the definitions and explanation above, it can be deduced that "Junction-to-ambient thermal resistance" should always be higher than "Junction-to-case thermal resistance" and "Junction-to-board thermal resistance". However, this is not the case: refer the screenshot from device LMK6X, Texas Instruments' Bulk-Acoustic Wave (BAW)

Rth-LMK6X

Example Calculation - 1: A transistor with power rating of 10 [W] and internal thermal resistance of 1.5 [K/W] has a case temperature of 60 [°C]. What is the actual value of junction temperature?

Here:
  • Junction-to-case resistance: ΘJC = 1.5 [K/W]
  • Heat generation rate ≡ Power dissipation rate: P = 10 [W]
  • Case temperature: TC = TJ - P × ΘJC. Thus: TJ = 60 [°C] + 10 [W] × 1.5 [K/W] = 75 [°C]

Example Calculation - 2: A chip with 5 [W] rating has a maximum junction temperature of 120 [°C] and an internal resistance of 0.5 [K/W] at an ambient of 40 [°C] with aluminium oxide wafers. What is the maximum permissible thermal resistance of the heatsink?

  • Junction-to-case resistance: ΘJC = 2.5 [K/W]
  • Heat generation rate ≡ Power dissipation rate: P = 5 [W]
  • Junction temperature: TJ = 120 [°C] As a "margin of safety", the permissible limit of junction temperature should be reduced by 20-30 [°C] from the value specified by the manufacturers.
  • Ambient temperature: TA = 40 [°C]. In case of natural convection (passive cooling) arrangment, due to the rise in temperature caused by heating of air between adjacent fins of the heatsink, the effective value of TA should be increased by a margin of 10-20 [°C].
  • Junction temperature: TJ = TA + P × [ΘHS + ΘJC]. Thus: 120 [°C] = 40 [°C] + 5 [W] × (ΘHS + 0.5 [K/W]). Hence, ΘHS = 15.5 [K/W]
Note that the above calculation method does not apply when multiple chips are connected to same heat sink. Using the catalogue of heat sink manufacturers, the heat sink [reference: Fischer Elektronik] can be selected as shown below.

Heat Sink Thermal Resistance Curve


Excerpts from "Semiconductor and IC Package Thermal Metrics" by Texas Instruments: RΘJA is a variable function of not just the package, but of many other system level characteristics such as the design and layout of the PCB on which the part is mounted. In effect, the test board is a heat sink that is soldered to the leads of the device. Changing the design or configuration of the test board changes the efficiency of the heat sink and therefore the measured RθJA. In fact, in still-air JEDEC-defined RΘJA measurements, almost 70 ~ 95% of the power generated by the chip is dissipated from the test board, not from the surfaces of the package. Because a system board rarely approximates the test coupon used to determine RΘJA, application of ΘJA using

[TJ = TA + RΘJA × P]

results in extremely erroneous values.

The document further elaborates: "In light of the fact that RΘJA is not a characteristic of the package by itself but of the package, PCB and other environmental factors, it is best used as a comparison of package thermal performance between different companies. For example, if TI reports an RΘJA of 40 [°C/W] or 40 [K/W] for a package compared to a competitor's value of 45 [K/W], the TI part will likely run 10% cooler in an application than the competitor's part." As a "margin of safety", the permissible limit of junction temperature should be reduced by 20-30 [°C] from the value specified by the manufacturers.


Excerpts from "Thermal Design Basics" by Analog Devices: In ICs, one temperature reference point is always the device junction, taken to mean the hottest spot inside the chip operating within a given package. The other relevant reference point will be either TC - the case of the device, or TA - that of the surrounding air. This then leads in turn to the above mentioned individual thermal resistances ΘJC and ΘJA. Taking the most simple case first, ΘJA is the thermal resistance of a given device measured between its junction and the ambient air. This thermal resistance is most often used with small, relatively low power ICs such as op-amps, which often dissipate 1 [W] or less. Generally, ΘJA figures typical of op-amps and other small devices are on the order of 90-100 [°C/W] for a plastic 8-pin DIP package, as well as the better SOIC packages.

Excerpts from a datasheet: "The LTM8023MP is guaranteed to meet specifications over the full –55 °C to 125 °C temperature range. Note that the maximum internal temperature is determined by specific operating conditions in conjunction with board layout, the rated package thermal resistance and other environmental factors."

LTM8023 Thermal Network

S. No.Field ParameterJESD51Application / Thermal SimulationRemark
01PCB1.6 [mm] thick, 4.0 x 4.5[in], 2-4 layers6-15 layers, high thermal conductivity1s0p, 2s2p (s: signal layer, p: power layer)
02Enclosure30 x 30 x 30 [cc]Smaller or biggerEnclosure shape and size irregular
03Board positionHorizontalCan be horizontal, vertical or inclinedMinimize convection to air
04Package mountingUpUp, Down, Back-to-backPackage can be mounted on both sides of PCB
05Measurement pointCentre of top caseMaximum in component volumeMaximum in volume is taken to make a safer estimate
06Heat dissipationTop face of package case onlyTop face, side faces and through pins to PCB: volumetric heat sourceRefer the image below for test set-up
07Calculation of RθJC(TJ - TC) / PTJ = TC + RθJC x PCASE where PCASE = heat transfer through case onlyRθJC remains constant
08Recommended methodRθJC is for 'theoretical' comparison of two designsUse ΨJT which is better representation of real world: TJ = TC + ΨJT x PRθJC measurement 'forces' all heat to casing

JEDEC51 Set-up


JEDEC51-8 Set-up


FPGA: Field Programmable Gate Array

As per Intel documents for Agilex 7 FPGA: "Thermal parameters do not include the traditional junction-to-case thermal resistance (θJC) and junction-to-board thermal resistance (θJB) values, due to its multi-chip package construction." The system-level thermal analysis of this product requires the use of its compact thermal model (CTM) in a computational fluid dynamic (CFD) tool. The CTMs are simplified mechanical models of the packages with modified thermal properties so they can predict an accurate case temperature with uniform power distribution for each die. The results of the CFD analysis are valid only to evaluate the TCASE which is temperature at the top center of the IHS - (Integrated Heat Spreader, case of an FPGA) of the package. TDP: Thermal Design Power, the power dissipated in a die that is used for thermal analysis purposes.

Virtex UltraScale and Virtex UltraScale+ are, ASIC-class architecture, FPGA product family from AMD. Excerpts from UltraScale+ Device Packaging and Pinouts Product Specification User Guide (UG575): "Unlike features in an ASIC or a microprocessor, the combination of FPGA features used in a user application is not known to the component supplier. Therefore, it remains a challenge for AMD to predict the power requirements of a given FPGA when it leaves the factory. Accurate estimates are obtained when the board design takes shape. For this purpose, AMD offers and supports a suite of integrated device power analysis tools to help users quickly and accurately estimate their design power requirements." The document further contains thermal resistance data with following notes: "The data includes junction-to-ambient in still air, junction-to-case, and junction-to-board data based on standard JEDEC four-layer measurements. The data in Table 10-1 is for device/package comparison purposes only. Attempts to recreate this data are only valid using the transient 2-phase measurement techniques outlined in JESD51-14. This data is not to be used in place of thermal simulation. Instead, refer to the thermal models provided for each device." In the chapter "Thermal Management Strategy", the document elaborates: "These resistances are measured using a prescribed JEDEC standard that might not necessarily reflect your actual board conditions and environment. The quoted θJA and θJC numbers are environmentally dependent, and JEDEC has traditionally recommended that these be used with that awareness."

AMD Flip-Chip BGA

Thermal interface material is needed because even the largest heat sink and fan cannot effectively cool an UltraScale or UltraScale+ device unless there is good physical contact between the base of the heat sink and the top of the UltraScale or UltraScale+ device. The surfaces of both the heat sink and the UltraScale or UltraScale+ device silicon are not absolutely smooth. This surface roughness is observed when examined at a microscopic level. Because surface roughness reduces the effective contact area, attaching a heat sink without a thermal interface material is not sufficient due to inadequate surface contact.

AMD Bare Die Flip-Chip

AMD advises against direct use of the θJC parameters to determine the thermal performance of the device in your application. The calculation of these parameters are done in accordance with the JEDEC standard JESD51 where system parameters differ greatly from most applications. Instead, run thermal simulations of the system in worst-case environmental conditions using Delphi thermal models, which more accurately represent the device thermal performance under all boundary conditions.

The thermal validation of FPGA thus needs to be performed using Compat Thermal Models (CTM) such as DELPHI model shown below (reference: Delphi Compact Thermal Model Guideline).

DELPHI Model

It is important to keep in mind that the availability of the DELPHI compact model does not eliminate the need for understanding the application in which the package is to be used. In other words, it is the user’s responsibility to take into account the environment surrounding the package.
Efficiency and Power Losses in MOSFETS

Reference: MOSFET power losses and how they affect power-supply efficiency - by By George Lakkas, Product Marketing Manager, Power Management

MOSFETs have a finite switching time, therefore switching losses come from the dynamic voltages and currents the MOSFETs must handle during the time it takes to turn on or off. MOSFET switching losses are a function of load current and the switching frequency of power supply. There are 3 types of losses: conduction losses, switching losses and static (quiescent) losses.

Efficiency of MOSFET

Most of the power is in the MOSFET gate driver. Gatedrive losses are frequency dependent and are also a function of the gate capacitance of the MOSFETs. When turning the MOSFET on and off, the higher the switching frequency, the higher the gate-drive losses. This is another reason why efficiency goes down as the switching frequency goes up.


MOSFET

Discrete power devices are the workhorses for power management and conversion of the most systems which fuel the central processing units and digital signal process ICs. Typical discrete products include various Diodes, Bipolars, metal-oxide-semiconductor field effect transistors (MOSFET) and Insulated Gate Bipolar Transistor (IGBT)s. A MOSFET may look like a black box from outside as shown below. However, inside the insulating enclusure contains many heat dissipating components consisting of silicon, copper and soldier joints. Capacitors and metal-oxide-silicon field effect transistors (MOSFETs) are the main elements in the ECU, which generate heat during operation, due to the Joule effect. In addition, capacitors and MOSFETs also generate heat due to electro-chemical reactions, dielectric losses and entropy changes. "The insulated-gate bipolar transistor (IGBT) was developed to overcome the drawbacks of bipolar transistor of low current gain, by combination of MOSFET and bipolar transistor."

MOSFET package

Reference: "Enhance Reliability of Semiconductor Devices in Power Converters by Minh Hoang Nguyen and Sangshin Kwak" - IGBT and SiC MOSFET have a similar chip-level structure, except for an additional p+ layer above the collector in IGBT and an additional body diode part in SiC MOSFET.

IGBT SiC MOSFET


Reference: AN90003 - LFPAK MOSFET thermal design guide

MOSFET LFPACK Thermal Guide

  • Part number BUK7M15-60E (LFPAK33, 15 mΩ, 60 V) the maximum thermal resistance junction to mounting base is 2.43 [K/W], the total thermal resistance, junction-to-ambient is 59.4 [K/W] when using 65.4 °C as (ambient) reference point.
  • Part number BUK7S1R0-40H (LFPAK88, 1 mΩ, 40 V) the maximum thermal resistance junction to mounting base is 0.4 [K/W], the thermal resistance between mounting base and ambient is of much higher value (~ 30 K/W)

The list of Joint Electron Device Engineering Council (JEDEC) standards related to thermal characterization of IC chips are:

  • JESD51: Methodology for the Thermal Measurement of Component Packages (Single Semiconductor Device)
  • JESD51-1: Integrated Circuit Thermal Measurement Method—Electrical Test Method (Single Semiconductor Device)
  • JESD51-2: Integrated Circuit Thermal Test Method Environmental Conditions—Natural Convection (Still Air)
  • JESD51-3: Low Effective Thermal Conductivity Test Board for Leaded Surface Mount Packages
  • JESD51-4: Thermal Test Chip Guideline (Wire Bond Type Chip)
  • JESD51-5: Extension of Thermal Test Board Standards for Packages with Direct Thermal Attachment Mechanisms
  • JESD51-6: Integrated Circuit Thermal Test Method Environmental Conditions—Forced Convection (Moving Air)
  • JESD51-7: High Effective Thermal Conductivity Test Board for Leaded Surface Mount Packages
  • JESD51-8: Integrated Circuit Thermal Test Method Environmental Conditions—Junction-to-Board
  • JESD51-9: Test Boards for Area Array Surface Mount Package Thermal Measurements
  • JESD51-10: Test Boards for Through-Hole Perimeter Leaded Package Thermal Measurements
  • JEDEC51-12: Guidelines for Reporting and Using Electronic Package Thermal Information

Chip Construction Details


Some other examples are:

chip Package Architecture

chip Junction Temperature


Commercial Tools for PCB Design and thermal simulations

  • ANSYS ICEPAK: A customized GUI for pre- and post-processing but uses FLUENT as solver.
  • SIEMENS FloTHERM: similar in look and operation like ICEPAK with primitives, SmartPart and attributes.
  • Electrical Packages such as PSPICE and Cadence (new owner of 6SigmaET): These packages primarily meant for design of IC chips has thermal modeling capabilities to check 'what-if' scenarios. Other EDA (Electronics Design and Automation and sometimes also used for Exploratory Data Analysis) tools are Allegro, Board Station, Expedition and CR5000.
  • ECXML format for Thermal Simulations is equivalent to STEP format for CAD tools. The data from ICEPAK to FloTHERM can be transferred in ECXML format which is a vendor neutral format.
  • DELPHI: Compact Thermal Model - the standard to create a thermal model for further use in thermal simulations

Heat Generation Rates in Electronic Parts

  • Power Density, PD: The power density of any device is the ratio of output power and the physical space it occupies. Thus, the unit of power density is [W/m3] or [W/cc].
  • Efficiency, η: This is the ratio of output power and the input power, η = POUT / PIN. The loss of energy (or power) gets converted into heat. Hence, when power density is increased, the volume and surface area of MOSFET and any such device will decrease. To keep temperature rise of the system within permissible limit, the amount of heat dissipation has to be decreased which in turn require that the efficiency of the device must be increased. The efficiency data is available in datasheets such as the one shown below for a DC-DC Step Down Power Supply (LTM8023):

    LTM8023 Power Supply

    Note that the calculated value as per method described below for heat gneration (Power Loss) does not match with the graph provided in the datasheet.
  • Heat generation rate = PIN - POUT = (1 - η).POUT = (1/η - 1).PIN and the efficiency values of are available in datasheets published by the manufacturers. For example, if the output of a device is 5 [A] at 3 [V] and efficiency is 75%, the heat generated in the device would be (1.0-0.75).(3 × 5) = 3.75 [W].
  • Heat dissipation: q = PIN - POUT = h x A x [TS - TAMB] where h is the convective heat transfer coefficient of the surface exposed to cooling media, A is the skin (fluid-wetted) surface area, TS is the surface temperature of the device casing and TAMB is the bulk mean temperature of cooling media.
  • Efficiency: η = POUT / PIN = POUT / [POUT + q]
  • η = 1 / [q/(V x PD) + 1] where V = volume of the device
  • η = 1 / [h x A x ΔT/(V x PD) + 1] where ΔT = [TS - TAMB]
  • Power Density of a MOSFET

MOSFET Cooling Arrangement

The conventional single side cooling 5x6mm MOSFET called SOP Advance (Figure 1) is covered with an insulating mold. Therefore, the heat dissipation path of the package is mainly its bottom drain side cooling plate.

The DSOP Advance package has a top source side cooling plate in addition to the bottom drain side plate as shown in Figure 2. These cooling plates contribute in reducing the thermal resistance (Rth).


Reference: TOSHIBA Technical Application Note

At thermal equilibrium, the maximum power dissipation PDMAX of a power MOSFET can be expressed in terms of the ambient temperature TA, the maximum channel temperature TCHMAX of the MOSFET and the channel-to-ambient thermal resistance RTHCH-A

Equivalent thermal network of a MOSFET

  • θi: Internal thermal resistance (channel-to-package)
  • θb: External thermal resistance (package-to-ambient air), varies with the material and shape of the case and is significantly larger than θi, θc, θs and θf.
  • θs: Thermal resistance of insulation shield
  • θc: Contact thermal resistance (at the interface with a thermal fin)
  • θf: Thermal resistance of the heat sink
Since no heat sink is generally used for power MOSFETs in a thermally conductive package RTHCH-A can be calculated as RTHCH-A = θi + θb. Since θb >> [θc + θs + θf]:

Simplified thermal network MOSFET

Reference: Numerical Study on Heat Transfer Characteristics of the 36V Electronic Control Unit System for an Electric Bicycle by Gihan Ekanayake, Mahesh Suresh Patil, Jae-Hyeong Seo and Moo-Yeon Lee.
ComponentMaterialThermal Conductivity [W/m-K]DensitySpecific Heat Capacity [J/kg-K]
CapacitorAluminium Oxide303890880
MOSFET Silicon80 ~ 150*2390712
Thermal GreaseSiC0.422500750
Remarkably, the thermal conductivity (in W/mK) defined in various papers and handbooks varies widely for pure Si.

MOSFET Thermal Network


Reference - [AN90003: LFPAK MOSFET thermal design guide] "The PCB material FR4 has a maximum operating temperature of around 130 [°C], depending on manufacturer and chemistry, this is much lower than the limit specified for the junction of a silicon die (175 °C)".

A PCB trace is a thin line of conducting copper placed on a non-conductive base material usually called FR4 that carries the signal and power to the whole circuit. A copper trace has a specific width called trace width, and a particular thickness. The thickness of PCB is specified in ounce/ft2 = 28.35 gm/ft2 = 0.3052 [kg/m2] = 0.3042/8900 = 3.43 x 10-5 [m] = 0.0343 [mm]. For typical PCBs, the most common copper thickness is specified as 35 [μm] which is equivalent to 1.0 [oz/ft2]. With electrical resistivity of 1.72 x 10-8 [Ω.m] at 25 [°C] for copper, the electrical resistance per inch width per meter length of copper trace is 0.02 [Ω]. The electrical resistance per inch width per inch length of copper trace would be 0.50 [mΩ]. The resistance value shall increase in direct proportion of the length and inverse proportion of the width and/or thickness (mm or oz/ft2).

  For example, electrical resistance per inch width per inch length of copper trace having thickness of 2.0 [oz/ft2] would be 0.50 / 2 = 0.25 [mΩ]. Electrical resistance per 0.1 inch width per inch length of copper trace having thickness of 2.0 [oz/ft2] would be 0.50 / 2 / 0.1 = 2.5 [mΩ]. Electrical resistance per 0.1 inch width per meter length of copper trace having thickness of 1.0 [oz/ft2] would be 0.50 * 39.37 / 0.1 = 0.197 [Ω].

 With temperature coefficient of electrical resistance 0.0039 [1/K], electrical resistance per 0.1 inch width per meter length of copper trace having thickness of 1.0 [oz/ft2] at 100 [°C] would be 0.197 x [1 + 0.0039 x 75] = 0.255 [Ω]

Trace Heating: In the electrical wires the cooling rate initially increase with thickness of insulation, reaches a peak and then starts decreasing. Similarly, the temperature of traces are dependent on the board thickness, traces on top of a thin board get hotter than one on a thicker board. This is because a thicker board has higher cross-section for the heat to conduct through. Beyond a thickness of board, there is more material under the trace than the trace can efficiently utilize - the path to surface on which convection occurs gets longer - and hence there is no improvement in temperature or in fact the temperature of traces may increase further. PCB trace temperatures are very sensitive to thermal conductivity though the in-plane thermal conductivity has higher influence than through-plane value.

IDF and IDX: ECAD - EMN/EMP

Some key features of such simulations are usage of Intermediate Data Format (IDF) and Incremental Data eXchange (IDX) files that have been exported from an ECAD package. These files contain informations of traces in PCB. In ANSYS ICEPAK, while importing traces the default materials are Cu-pure for metal and FR4 for dielectric. PCB construction is a layered design along the thickness direction and hence the thermal conductivity is necessarily orthotropic. Here, the conductivity value along the thickness direction - known as through-the-plane conductivity is far less than in-the-plane conductivity values. .emn file supports importing the data into a part that would represent the PCB or panel without any component placement. The .emn file has a section for the placement info with the header that can be edited in a text editor to delete it. This section starts with ".PLACEMENT" should you need to remove the placement data from the file. The IDF files are text files and can be easily manipulated by scripting tools such a Python and PERL. Section 3.10 of the IDF 3.0 specification which defines how drilled holes are handled. In summary:
  1. EMN File Contains:
    • PCB Outline
    • Component Location
    • Component Orientation
    • Hole Information
    • Keep In and Keep Out Regions
  2. EMP File Contains:
    • ECAD outlines for every component
    • ECAD component height information

There are many programs used to generate design of electronic items such as PCB and chips. Cadence, EasyEDA, Eagle, Altium Designer, Siemens Expedition Enterprise to name few. However, there is a need to have interface between ECAD and MCAD. This is accomplished by EMN and EMP files. The EMN file is related to the board and the EMP file is the library list of the components on the board. These files can generate a 3D view of the PCB Assembly. IDF files default to the following extensions: *.emn - neutral file of the board outline and component placement and *.emp - profile file that contains component outlines.. A sample EMN file commented with explanation of the information can be found here.

While reading .EMN file in other MCAD programs such as ANSYS Discovery or Creo, the model tree may contain only the names of size designators and reference designator or actual part number may be missing. A reference designators are combination of letters and numbers assigned to PCB components, such as resistors, capacitors, and other electronic elements which provide a standardized way to identify and reference components on the board.This may only give a good representation of the board but one cannot do further processing such as assigning heat generation rates, material properties... To get Mechanical Models with different part names, ECAD Design Software should be set to output package types as the ecad_name and (a) either an internal corporate part number or (b) a vendor part number as the ecad_alt_name. Do not set up ECAD IDF output to use identical ecad and ecad_alt names. AR = Amplifier, C = Capacitor, D = Diode, F = Fuse, FB = Ferrite Bead, J = Connector / Jack Connector, K = Relay, L = Inductor, LED = Light Emitting Diode, M = Motor, P = Plug, PS = Power supply, Q = Transistor, R = Resistor, S = Switch, T or XMER = Transformer, TP = Test Point, TR = Transistor or transducer, U = Integrated Circuit. For example, R1 might refer to the first resistor on the board, C2 could be the second capacitor, and U3 might represent the third integrated circuit. More information at resources.altium.com/p/altium-designer-helps-you-track-reference-designators-your-pcb such as the image below.

Refererence designators

Reference: simplifiedsolutionsinc.com/images/Steps-to-create-3D-PCBs-in-ProE.pdf: In Pro/E (predecessor to Creo), there was an option to create an ecad_hint.map file which linked components in CAD Tool to 3D Pro/Engineer Model. When an IDF (emn) File is imported into Pro/Engineer, Pro/E cross-referenced the ecad_name and ecad_alt_name from the emn file against the ecad_hint.map. If a match was found, Pro/Engineer replaced "on the fly" geometry with a real 3D Pro/Engineer part or assembly. More information can be found at support.ptc.com/.../Map_File_Standard_Conventions.html

Sample ecad_hint.map file. Each section begins with the purpose, followed by '->'. Each section ends with 'end'. # is the comment character. Wildcard (*) is valid for 'all'. Object and value fields are separated by a space, spaces are permitted in value strings if the string is surrounded by quotation marks.

mcad_in_ignore ->
ecad_name "resistor"
ecad_alt_name "res_5"
ecad_type "part"
ref_des "*"
end
As per Alitum documentation under "Mechanical Data Import-Export Support": "In the STEP file, each component is identified by its designator. If the MCAD designer needs to import multiple boards into a single MCAD file there is likely to be designator clashes, to avoid this include a Component Suffix."

PWA = Printed Wiring Assemblies. Excerpts from "Intermediate Data Format Specification, Version 3.0": Structure of the Intermediate Data Format - The Intermediate Data Format consists of three files: the Board File, the Library File and the Panel File.

The Board File: It contains a description of a single PWA, including the board shape, layout restrictions and component placement.

The Library File: It contains descriptions of components used by one or more PWAs. *.emp is the extension of the library file.

The Panel File: It contains a description of a manufacturing panel including the panel shape, layout restrictions and the placement of boards and components on the panel.


Data is organized by sections in these files. Each section begins with a keyword indicating the type of data the section contains and a matching keyword at the end of the section. All data between the section keyword and its corresponding ending keyword pertains to that section. Sections cannot be nested. Unless otherwise noted, sections within a file can be in any order. Data within the sections is represented by one or more records consisting of one or more fields. Each line in a file is a separate record: fields within a record are separated by one or more blanks. Records within a section and fields within a record must be in a specific order. Records are free format which means that the fields they contain can be any length, and each field can begin in any column as long as the order of fields is maintained.

The Panel File is an optional file, similar to the Board File, that contains the physical description of a manufacturing step-and-repeat panel and the locations of boards and components on that panel. The Panel File references one or more PWAs described in separate Board Files. Any component placed on the panel itself is referenced in a Library File.

Note: "The comment character is the pound or hash sign (#). A comment must be a separate line (record) and the comment character must be in column 1. Comments should be located between, but not within sections of the IDF files." The Header section must be the first section in the file, the second section must be the Outline section, and the last section must be the Placement section. All other sections may be in any order. Exporting IDF files in Atlium wil1 generate two files - one containing information about the physical size and shape of PCB and positions of components, the other containing information about each component including name, size, and shape. These are typically referred to as the board and library files, respectively. Different CAD packages use different file extensions for the board and library files. The board file and library file extensions of generated files have following pairs: .brd and -pro, .brd and .lib, emn and emp, . bdf and.ldf, .idb and .idl, .idf and .lib. Note that EMP is for describing parts that equip the board. Thus, for a PCB (without parts on it), no EMP file will be generated.

The Gerber file

It is a connector and bridge between designers, engineers and PCB manufacturers. It needs to go through every manufacturing process and the factory can clearly define the customer's needs. According to UCAMCO (the company that currently owns the rights to Gerber File format): "Gerber file format" is a standard for PCB design data storage or transfer. Gerber file describes and communicates the constituents of a PCB image like the number of copper layers, solder masks and many others such attributes. Gerber files also act as input files to PCB printing devices like photo-plotters and Automated Optical Inspection (AOI) machines to print or compare circuit board images for different gadgets. Gerber files may also include metadata (data about other constituting data within a file) like solder mask, legend/silk and number of copper layers among other relevant printing information.

Following files comprise the full list of GERBER files which are usually zipped and shared to thermal simulation engineer or PCB Manufacturer.

Gerber File TypeExtension
Top side (copper) Layer.GTL
Bottom side (copper) Layer.GBL
Top Overlay.GTO
Bottom Overlay.GBO
Top Paste Mask.GTP
Bottom Paste Mask.GBP
Top Solder Mask.GTS
Bottom Solder Mask.GBS
Keep-Out Layer.GKO
Drill Drawing.GD1
Drill Guide.GG1
Internal Plane Layer 1, 2 ... 16.GP1, .GP2 ... .GP16
Ref: resources.pcb.cadence.com/blog/what-is-documentation-in-pcb-manufacturing - "Manufacturing Data File Formats: The manufacturing process for PCBs have relied on Gerbers, which is an open ASCII vector format file of 2D binary images. However, there are also other file formats that that you can use such as ODB++ and IPC-2581; keep in mind that it depends on which file format your manufacturer accepts. Manufacturing data are files that will be used to create the etched base laminate on your PCB, including copper layers, solder mask, and other elements."

"Electronic Schematics: Schematics communicate specific information about how electronics in a design should be connected to each other. Components are labeled with their electrical characteristics, such as capacitance or resistance, with the complete circuit being illustrated across the PCB. The schematic is an organized view of the electrical circuit and provides necessary information for manufacturing."


Typical *.emn File Entry
#------------------------------------------------------------------------------
#The comment character is the pound sign (#). A comment must be a separate line
#(record) and the comment character must be in column 1. Comments should be 
#located between, but not within sections of the IDF files.
#------------------------------------------------------------------------------
.HEADER
BOARD_FILE 3.0 "Sample File Generator" 2020/07/01.16:02:44 1
sample_board THOU
.END_HEADER
#Unit MM or THOU = milli-inch
#------------------------------------------------------------------------------
#Section keyord: .BOARD_OUTLINE or .PANEL_OUTLINE
#MCAD - Outline is owned by the Mechanical system and should not be modified in
#the Electrical system
#ECAD - Outline is owned by the Electrical system and should not be modified in 
#the Mechanical system
#UNOWNED - Outline can be modified in either system
#62.0 is thickness in milli-inch
#------------------------------------------------------------------------------
.BOARD_OUTLINE MCAD
62.0
0 5030.5 -120.0  0.0
0 5187.5 -120.0  0.0
0 5187.5  130.0  0.0
0 5155.0  130.0  0.0
0 5155.0  550.0 -180.0
...
.END_BOARD_OUTLINE
#------------------------------------------------------------------------------
.ELECTRICAL
EthnetBrd 135792468 THOU 59.0
0 -92.0  63.0 0
0 -92.0 -63.0 0
0  92.0 -63.0 0
0  92.0  63.0 0
0 -92.0  63.0 0
.END_ELECTRICAL
#------------------------------------------------------------------------------
EthnetBrd 135792468 E7
20.5 32.5 0.0 90.0 BOTTOM PLACED
  • EthnetBrd is the ECAD name
  • 135792468 is the ecad_alt_name (unique part number)
  • E7 is the reference designator
  • 20.5, 32.5 0.0 represent the location with respect to the origin of PCB (0,0) point X, Y, Z directions
  • 90.0 represents the orientation of the component
  • BOTTOM specifies the side of the PCB
  • PLACED specifies component will be added to 3D PCB

These softwares are meant for electronics industry mainly and hence contains many built-in objects to expedite the simulation process. They can be summarized as follows:

Electronic Cooling Software

  • Primitive - Fundamental geometric entities in FloTHERM and ICEPAK: Cuboids, Prisms and Flow Resistances
  • SmartPart - Object parametrically created out of Primitives: e.g. Enclosure, fan, PCB, cylinders, volume or surface heat source, heat sinks (described by base dimension, number of fins, fin width and fin height), perforated plate (fully designated by hole size and arrangement, pitch, free area ratio)
  • Assembly - A group of Primitives, SmartParts and Sub-Assemblies
  • Attribute - A property that can be attached to Primitives and SmartParts (e.g. material properties)
Primitives defined in ICEPAK

ICEPAK Primitives


Reading Mechanical CAD (MCAD) data
  • The MCAD data can be read either in their native format such as ProE, Solidworks or CATIA files (parts and assemblies) or neutral formats such as IGES, STEP or PARASOLID.
  • After initial defeaturing (removal of chamfers and fillets) and simplifications (removal of small holes, branding logos, part identifiers), the MCAD geometry needs to be converted into ICEPAK / FloTHERM entities.
  • Simplification is permissible to the to the extent where the geometry could be created manually using primitives and smart-parts in ICEPAK / FloTHERM!
  • The conversion into ICEPAK / FloTHERM entities is a process of replacing a detailed geometry say perforated plate with a plate of same size without perforation and specifying the perforation details as attributes to the ICEPAK / FloTHERM entities.

Mesh generation process and recommendations
  • For naturally convected cases, the computational domain needs to be made bigger than the chassis or the board.
  • For the space above the chassis extend the domain 2~3 × height of chassis. Note that the buoyancy may not be strong enough to establish a flow from lower face to upper face and reverse flow can be observed on both of these faces. In order to reduce the reverse flows, the domain boundaries may need to be kept closer to the heat source.
  • For space below chassis extend the domain equal to the height of chassis.
  • For the remaining sides (front, rear, left and right) of the chassis, extend the domain 0.5 ~ 1.0 × depth and width of the chassis.
  • In case on natural convection, the heat transfer through convection is low and hence radiation also contributed significantly. Up to 50% of heat transfer is by radiation and remaining 50% by natural convection.This can be justified by the fact that the HTC value for natural convection in air varies between 5~10 [W/m2.K] whereas equivalent HTC for radiative heat transfer = ε × σ × (TWALL4 - TAMB4) / (TWALL - TAMB) = 4.9 [W/m2.K] for ε = 0.5, TWALL = 100 [°C] and TAMB = 50 [°C].
  • The components are mostly represented as volume with sharp corners such as cuboid and rectangles. A cut-cell method also known as trimmer mesh or Cartesian mesh or snappyHexMesh (in OpenFOAM) is used in ICEPAK and FloTHERM.
  • For printed circuit boards as in any surface with a significant amount of heat flux, 3 cells in the first millimeter above the board (air volume) and 3 cells in the first millimeter below the board (solid volume) is recommended.
  • Maintain aspect ratio of cells < 100.
  • For fins, maintain > 5 cells cross a channel for better heat transfer and pressure drop prediction, use two or more cells across the thickness in the solid.

Special applications in Electronic Cooling
  • Thermo-Electric Cooling (TEC): This device is based on Peltier effect where thermal energy is absorbed at one dissimilar metal junction and discharged at the other junction when electric current flows within a closed circuit. It comprises of p-type and n-type semiconductors sandwiched between ceramic electrical insulators. TEC are solid state heat pumps for applications where cooling below ambient are required. The cold junction acts as 'evaporator' and hot junction as 'condenser' of a refrigeration cycle.

    TEC - Thermo-Electric Coolers

    • Seebeck Effect: ΔV = α × (TH - TC) where α [V/K] = differential Seebeck coefficient or (thermo electric power coefficient) between the two materials, positive when the direction of electric current is same as the direction of heat flow.
    • QH [W] = β [V] × I [A] where β is differential Peltier coefficient between given two materials
    • β < 0: Electric current and heat flow in opposite directions
    • β > 0: Electric current and heat flow in same directions
    • COP = coefficient of performance of the thermoelectric device = QC / J, which typically is between 0.4 and 0.7 for single stage applications.
  • Heat Pipe: Capillary effect and phase change (evaporation and condensation) are the phenomena which define operation of heat pipes. Due to phase change from liquid to vapour, the heat transfer coefficient for heat pipes is extremely high - of the order of 40,000 [W/m2.K]. Originally invented by NASA for space application, it has gained widespread application in electronics industry.

    Heat Pipe Section View

    Heat Pipe Appliction in CPU Cooling

  • Heat Spreader: This is an application similar to heat sink. The purpose of a heat spreader is to use material with a very high thermal conductivity such as graphite with k = 1400 [W/m-K] to make the heat flow in-plane over a larger area so that it can be further dissipated into ambient using heat sinks.
  • PCM - Phase Change Material: These are energy storage and release mechanism based on change of phase (typically solidification and melting). The material can be used to keep temperature fluctuations low in case of heating and cooling cycles. Sometimes, the PCM can also act as an insulator to heat dissipation based on thermal conductivity. Some key characteristics required for a energy storage and release type PCM are tabulated below.
    CharacterisitcsDesired valueRemark
    Melting point, TMPAs per temperature controlSelection of material will depend on temperature to be maintained
    Specific heat capacity, CpHighEnergy storage capacity ∝ Cp. Higher the Cp, lesser the mass required to store a given amount of energy.
    Density, ρHighEnergy storage capacity ∝ ρ and the volume required is also less as m = ρ * V
    Thermal Conductivity, kHigh for energy storage purposeLow value is required for insulation where heat is to be maintained near the source itself
    Coefficient of volume expansion, γLowThis governs flexibility or void space required in the storage container
    Chemical compatibilityNon-corrosiveShould not react with the container and other materials in case of leaks
    Thermal cycling (heating-cooling) stabilityNo degradationThe micro-structure and material properties should not degrade with heating-cooling cycles

TEC Performance Calculation

Field Variables

  1. i = Current [A]
  2. k = Thermal conductivity of TEC block [W/m-K]
  3. QH = heat rejection from hot surface (the 'condenser') [W]
  4. QC = heat absorption by cold surface (the 'evaporator') [W]
  5. R = Electric Resistance of TEC assembly [Ω]
  6. TC = cold side temperature (the 'evaporator') [K]
  7. TH = hot side temperature (the 'condenser') [K]
  8. P = power consumption by TEC device = (QH - QL) [W]
  9. S = seebeck co-efficient [V/K]
  10. ΔT = temperature difference between hot and cold junction [K]
  11. G = Ratio of cross-section area to the length (thickness) of the device [m]
  12. ρ = resistivity of TEC device [Ω-m]
  13. COP = Coefficient of Performance = QC / P

TEC Analogy

QC = -[S.i.TC − 0.5 × i2R − k.×G×(TH − TC)], negative sign is for heat flow into the device.

Or

QC = -[S.i.TC - 0.5 × i2ρ/G − k×G×(TH − TC)]

Similarly:

QH = [S.i.TH + 0.5 × i2R − k.×G×(TH − TC)]

Or

QC = [S.i.TH + 0.5 × i2ρ/G − k×G×(TH − TC)]

TEC Temperature Profile


Anisotropic or Orthotropic Thermal Conductivity: The printed-circuit boards are formed by many layers of copper wires known as traces and dielectric material (say FR4). They are so thin that they cannot be modeled (meshed and boundary conditions applied) separately. Hence, the thermal conductivity of board can be simplified using equivalent uniformed value along the thickness and in-the-plane direction. This simplification is done using series and parallel arrangement of thermal resistances analogous to electrical resistances.

Orthotropic conductivity based on lumped block assumptions. 10% Copper is a reasonable guess. In case one needs to use "Locally Varying Orthotropic" thermal conductivity, the detail layout of traces and FR4 layers need to be modelled.

  • kIN_PLANE = VFCu × kCu + (1-VFCu) × kDie
  • 1/kX_PLANE = VFCu / kCu + (1-VFCu) / kDie
  • Here, VF = Volume Fraction of copper in fractions (that is the value lies between 0 and 1).
CFD mesh vs. background mesh:

PCB Background Mesh in ICEPAK

PCB Local Thermal Conductivity in ICEPAK

A bit more detailed calculation of thermal conductivity is layer-by-layer estimation also known as "Discrete Layer Stack-up". For each layer:

kLAYER_i = kCu × ACu / APCB + kFRP × AFRP / APCB

Discrete layer stack-up in PCB

A detailed calculation of effective thermal conductivities in X-, Y- and Z-directions can be estimated by spliting the PCB in smaller segements and using smaller patches as described below. This method is useful for capturing conduction paths near heat generating components and heat sinks.

detailed Thermal Conductivity PCB

detailed Thermal Conductivity X, Y, Z directions


ICEPAK is a GUI for pre- and post-processing. It uses FLUENT as solver and in this process many files get created. Following is a list of files and its owner (ICEPAK or FLUENT?).

File Type Created by Used by Suffix / Filename Remark
ModelICEPAK ICEPAK model
Problem ICEPAK ICEPAK problem
Job ICEPAK ICEPAK job
Mesh inputICEPAK  meshergrid_inputInputs for the mesh generator.
Mesh outputmesherICEPAK grid_outputOutput from the mesh generator that is the mesh file
CaseICEPAK FLUENT .casContains all the information that is needed by ICEPAK to run the solver
DataFLUENT FLUENT .dat and .fdatFiles when it has finished calculating: *.dat and *.fdat. These data files can be used to restart the solver
ResidualFLUENT ICEPAK .resInformation about convergence monitors: Solve → Solution monitor or select Convergence plot in Post menu
ScriptICEPAK ICEPAK.SCRIPT or _sc.batRuns the solver executable and can also be used to run the solver in batch mode.
Solver inputICEPAK FLUENT .uns_inThe solver input file (projectname.uns_in) is read by the solver to start the calculation.
Solver outputFLUENT .uns_outInformation from solver that is displayed on screen during calculation - this file is written only on Linux systems
DiagnosticICEPAK .diagContains information about correspondence between object names in model file and object names in case file
OptimizationICEPAK optimizer.log, .dat, .tab, .post, .rptOptimization of field variables
PostprocessingFLUENTICEPAK .resdUsed by ICEPAK for post-processing. All solutions that exist for the current project are listed by solution ID.
Log ICEPAK ICEPAK .log
GeometryExternalICEPAK.igs, .stpCAD geometry - input to ICEPAK
PackagedICEPAK ICEPAK.tzrProject archive

Few keyboard short-cuts and special topics in ICEPAK
  • Move legend with Ctrl and the middle mouse button
  • Edit levels and set orientation with shift right-click on legend
  • Move the cut-plane in the domain with shift and the middle mouse button
  • Shift middle-click on CAD objects to graphically move the CAD geometry
  • Thermal Chokepoint: the dot product of heat flux and temperature gradient - it shows regions of high heat flux experiencing large thermal resistances
  • Thermal Cross: the cross product of heat flux and temperature gradient - it shows regions where large heat flux vectors not aligned with high thermal gradients

Basic Solver Setting
Note that ICEPAK uses FLUENT in background as solver.

Solver Settings in ICEPAK


Mesh Generation
The default option in ICEPAK is to generate a Cartesian mesh similar to snappyHexMesh in OpenFOAM, trimmer mesh in STAR-CCM+

Mesh Generation Settings in ICEPAK

Mesh Display Settings in ICEPAK

Mesh Quality Check in ICEPAK


Fan, Fins and Grilles
There are many in-built feature to model flow resistances and momentum sources / sinks.

Types of Fans in ICEPAK

Fan swirl: tangential component of flow at exit of a fan.

Modeling of Fan Swirl in ICEPAK

Heat Sinks:

Fin Types in ICEPAK

Fins in ICEPAK

Fin Geometry in ICEPAK

Grilles and Louvres:

Grilles in ICEPAK


Boundary Conditions
Wall boundary condition:

Wall B.C. in ICEPAK

PCB Gerber File Import

pcb Import Gerber in ICEPAK

PCB Stack Data

pcb stack data in ICEPAK

PCB Stack-up Data in Detail

pcb geometry detailed input in ICEPAK


Monitors and Runs Settings

Monitor Points in ICEPAK

Solver Run Settings in ICEPAK


Micro-Channels
Cooling systems in the field of microelectronics and microelectromechanical systems (MEMS) use small and microscale flow passages. The transport phenomena in microscale channels are very different as compared to conventional size channels or macroscale channels. There are several dimensionless numbers used to represent the feature of fluid flow in microscale channels, similar to Reynolds number in macro-chnnels. The distinction between macro- and microscale channels may also be classified according to these dimensionless numbers.
Heat Transfers from Extended Surfaces (Fins)
This section provides a way to calculate heat transfer rate and estimate typical range of values for designs. Two design criteria used to define performance of fins are "fin efficiency" and "fin effectiveness". Fin efficiency is ratio of "actual heat transferred through fins" to the "heat that would be transferred if the entire fin were as base temperature".

Fins Plate Type

Such fins are extruded or machined and needs to be fixed on the heat generating component using adhesive or screws. This creates a thermal contact resistance and needs to be accounted for temperature calculation of the mating surfaces. The parameter ΔT/[Q.A] can be used in fin selection for given heat dissipation and surface area (available space). Lower the number, better the heat sink design.

BFThickness of fins[mm]122255
LFHeight of fins[mm]202010202010
WFWidth of fins [depth perpendicular to screen][mm]250250250250100100
AFCross-section of fins = BF×WF[mm2]250500500500500500
PFPerimeter of fins = 2[BF+WF][mm]502504504504210210
kFThermal conductivity of fin[W/m-K]100100100100100100
hConvective Heat Transfer Coefficient[W/m2-K]101020101010
NFFin parameter[m-1]14.1710.0414.2010.046.486.48
RFThermal resistance of each fin[K/W]10.2310.059.9910.0523.9447.69
iFNumber of fins[nos]101010101010
RFThermal resistance of all fins[K/W]1.021.010.991.012.394.77
ΔTTemperature potential[K]606060606060
QHeat transfer rate (all fins)[W]58.759.760.159.725.112.6
AHeat transfer surface area[cm2]101051042
ΔT/[A.Q]Heat transfer rate[K/W/cm2]0.1020.1010.2000.1010.5992.384

Fins Pin Type

Following calculations is valid for pin-type and plate-type fins only.
Select the type of Fin:
Thermal conductivity of fin in [W/m-K]:
Specify length or height of fin [mm]:
Diameter (pin-type) or thickness (plate-type) of fin [mm]:
Width of plate-type fin [mm]:
Convective HTC on fin surface [W/m2.K]:
Reference temperature [°C]:
Base temperature [°C]:
Base thickness [mm]:
Tip convection [W/m2.K]:

The heat transfer rate of pin-type fin for natural convection in air is tabulated below. Note the impact of L/D ratio on heat transfer rate.

L[mm]100.0100.0100.0100.0100.0100.0100.0
D[mm]5.010.020.025.050.0100.0200.0
h[W/m2-K]10.010.010.010.010.010.010.0
ΔT[K]60.060.060.060.060.060.060.0
NF[m-1]8.96.34.54.02.82.01.40
A[cm2]15.731.462.878.5157.1314.2628.3
Q[W]0.75201.66803.53704.47609.1818.6037.45
ΔT/[A.Q][K/W/cm2]5.0791.1450.2700.17070.04160.01030.0025

The values for forced convection in air is tabulated below. 'A' is the heat transfer area of the fin and not the cross-section area. The parameter ΔT/[Q.A] can be used in fin selection for given heat dissipation and surface area (available space). For example, for a heat dissipation of 10 [W] using a fin of diameter 20 [mm] and height 100 [mm], the expected increase in temperature of the base of fin is 10 [W] x 62.8 [cm2] x 0.1431 = 89.9 [K].

L[mm]100.0100.0100.0100.0100.0100.0100.0
D[mm]5.010.020.025.050.0100.0200.0
h[W/m2-K]20.020.020.020.020.020.020.0
ΔT[K]60.060.060.060.060.060.060.0
NF[m]8.96.34.54.02.82.01.40
A[cm2]15.731.462.878.5157.1314.2628.3
Q[W]1.2703.0086.6738.53317.90536.72574.409
ΔT/[A.Q][K/W/cm2]3.00700.63500.14310.08950.02130.00520.0013

The chart below provides a comparison between two designs of circular fins. Here, L ≡ r2 - r1 and 'A' is based on section shaded in dark colour. Reference: A Textbook of Heat Transfer by Lienhard and Lienhard.

Circular disk fins

Triangular fins


More on LED

SMT High Power LED


LED Die Construction


Thermal Resistances

LED Analytical Solution

Refer to fischerelektronik.de/fileadmin/fischertemplates/download/Katalog/heatsinks.pdf for types of heat sinks and their thermal resistances. As per the document, the values indicated in the diagrams apply only for heatsinks with black anodised surface, mounted vertically and natural convection. Correction factors: natural surface: +10 to 15%, for horizontal mounting: +15 to 20%.

Thermal Simulations in Data Centres

The simulations of data centres can be groupled into sub-systems and multi-scale models. The super-system is the room itself followed by systems such as racks and plenum/ductings following by sub-systems like servers, server-mounted heat exchangers, appliances, mechanical systems (fans, blowers, dehumidifiers), cold-plates, and even Insulated-Gate Bipolar Transistor (IGBT). The simulation gets further complicated by change in weather conditions during the day and over a longer time periods.

Data Centre

Traditional cooling systems are air-based in which hot air from the servers is pulled to the back row of the racks into a "hot aisle" plenum and then pushed to an air handling unit where air-to-water heat exchangers remove the heat, recuperate the air and push the fresh, cold air back into the data centre. This is centralized way of handling the cooling media. Another approach is to mount the heat exchangers directly at the back row of the racks which are supplied with cold water to cool the hot air collected from the servers. For failure mode analysis, a transient simulation or modeling would be needed to check how long the system can be allowed to operate in safe limit due to thermal inertia of the coolant and solids.

Data Centre Cooling

The cooling and thermal management of data centres are similar to any other HVAC systems. Some key abbreviations used in HVAC industry for buildings are also used in data centre cooling technology. APD: Air Pressure Drop, CHWR: Chilled Water Return, CHWS: Chilled Water Supply, CRAC: Computer Room Air Conditioning, CRAH: Computer Room Air Handler, DCIM: Data Center Infrastructure Management, DCV: Demand Controlled Ventilation, EAHU Exhaust Air Handling Unit, FCA: Fan Coil Assembly, FCU: Fan Coil Unit, HAHX: Humid Air Heat Exchangers, HDD: Heating Degree Days, HEPA: High Efficiency Particulate Arresting, LAT / LWT: Leaving Air/Water Temperature, PTAC Packaged Terminal Air Conditioner, TEFC: Totally Enclosed Fan Cooled, VAV: Variable Air Volume, VOC: Volatile Organic Compounds. CRAC units are used for data centers with lower electrical loads and CRAH are used in case for higher electrical loads. CRAC used refrigerant, CRAH used chilled water. Service Level Agreement (SLA) for energy efficiency is an assurance or guarantee to deliver a service at lower than a specified Watt per unit of measure.
Contact us
Disclaimers and Policies

The content on CFDyna.com is being constantly refined and improvised with on-the-job experience, testing, and training. Examples might be simplified to improve insight into the physics and basic understanding. Linked pages, articles, references, and examples are constantly reviewed to reduce errors, but we cannot warrant full correctness of all content.