2021-04-13T13:57:38.000-0400

Other Hardware

The major cause of system breakdown used to be overheating of critical components. But failure analysis shows that in today’s system designs, this is not the only problem.

Attachments: | Siemens SW Using Simcenter to test the thermal performance of electronics Solution Guide.pdf (7 MB) |

The major cause of system breakdown used to be overheating of critical components. At the hottest points in some of these systems, semiconductor junctions can reach 150°C or more, getting very close to the limits of their operation. These high temperatures will modify and eventually destroy the circuit operation unless excess heat is removed from the chip. But failure analysis shows that in today’s system designs, this is not the only problem.

Typical component breakdown can also be caused by repeated thermal transients. Heating and cooling induces shear stresses at the material interfaces in the package structure (die attach, solder joint), which result in delamination, tear off, etc. The lack of heat removal through the resulting diminished contact area can cause thermal runaway.

For a long time, packaged electronic components were represented by a single thermal resistance in the datasheet. Power devices were usually encapsulated in packages with a dedicated cooling surface, called the “case.” In discrete semiconductors (diodes, transistors), the hottest portion of the device was the pn-junction. As a first estimation, the engineer accepted that the temperature rise of the junction above the surface on which the package was mounted was the supplied thermal resistance multiplied by the applied power.

In conventional electronics as well as in solid-state lighting, the junction temperature (TJ) is a primary

quantity that influences system reliability and lifetime. The junction temperature of an LED is a performance indicator of the thermal design, and many properties of the light output of an LED depend on the absolute junction temperature.

Single thermal resistance values such as the junction-to-case (R

The thermal transient measurement method was developed to address these issues, providing a better solution than using thermal sensors. Today, accurate thermal characterization is needed to be able to create an optimal design for electronic systems. Thermal characterization techniques help to provide better results when developing models for simulation. Using the combination of thermal characterization, thermal transient testing and 3D thermal modeling saves time and costs by reducing the need for iterations of physical prototypes, re-spins during production and recalls of defective products from the field.

At its basic level, thermal transient testing means applying a steady, low power level to the device and then instantly switching to a higher level to observe a heating transient. Similarly, switching from a higher power level to a lower one, the cooling of the device can be monitored. These transients can be completely captured until the steady state is reached. For many more details on the math, theory and techniques behind the development of thermal transient testing techniques from this base point to what is commonly used today, see the book Thermal Management for LED Applications

The development of thermal transient measurements for model creation has culminated in the Cauer-ladder network model. This model is useful for associating circuit elements with physical regions. The behavior described by the model is the basis for identifying the heat-flow path using “structure functions.” Researchers at the University of Budapest developed structure-function analysis methods to fulfill the need for greater interpretation inside packages

Structure functions transform the thermal transient measurement result into a thermal resistance versus thermal capacitance profile, which provides detailed thermal information of each layer the heat passes through, from junction to ambient. This enables identification of the physical characteristics of layers such as die attach, base, package, heatsink and even cooling devices such as a fan.

With this method, engineers are now able to identify the die-attach failures of single-die packages

Some portions of an LED package are highly stable (for example, the chip, submount, heat slug) However, the thermal interface material (TIM) layers, which are applied to fill the small gaps between attachment surfaces, can show large variation even in samples from the same manufacturing lot. Testing the TIM itself yields no hint on the actual thermal resistance achieved in the fabricated TIM layer. These unavoidable variations in the product can be best studied using structure functions.

Figure 1 shows the typical thermal interfaces in an LED application. Their quality can be measured as the length along the thermal resistance axis as shown in Figure 2. There could be many reasons why interfacial thermal resistance changes: curing/soldering temperature change, variation in the thickness of the TIM layer, aging, or the quality is changed on purpose.

In the setup shown in figure 3, the ∆T

As the sample thickness is changed with a precise, dedicated mechanical system, the total measured thermal impedance of the complete test setup will change. The structure functions in Figure 4 show that the change of the total junction-to-cold plate resistance of the test setup was caused solely by the changes in material sample thickness. The thermal behavior of the power diode can be seen to be unchanged, as are the interfacial resistances on either side of the sample under test.

The λ thermal conductivity of the sample can be calculated as follows:

where A is the cross-sectional area of the heat flow path across the sample, ∆L is the bondline thickness change, and ∆R

According to the equation, the thermal conductivity of the TIM sample tested is directly proportional to the slope of the R

There are some advantages of this TIM testing method compared to other techniques. For example, measurement uncertainty is reduced by the differential approach. The test method is a quasi in-situ technique because the test fixture resembles the real-life application conditions of TIM materials. Last but not least, every measurement includes an inherent self-test of the measuring system. Based on the obtained structure functions, the structural integrity of the test fixture can always be checked. This approach is implemented in the Simcenter DYNTIM™ software

A simulation model is only as good as the available input data, that is, the device geometry and material properties. This is always an issue when creating detailed models for a computational fluid dynamics (CFD) simulation tool, even for semiconductor vendors who in principle must be aware of at least the device geometry. But many times, material parameters and effective volumes or area raise questions. As noted, a possible source of uncertainty in detailed thermal simulation models is the interface thermal resistance, both at TIM1 (die attach) and TIM2 (for example, thermal grease).

The idea behind using structure functions for detailed model calibration/validation is that if both the geometry/material properties in the simulation model and the boundary conditions correspond to the real-life situation, then there should be no differences between the measured and simulated thermal-impedance curves.

Any small difference caused by geometry or material mismatch should, therefore, be visible in the corresponding structure functions. The following case study shows how a power semiconductor device package model could be fine-tuned with the help of structure functions

Figures 6 through 8 show the major iteration steps taken to adjust the simulation model. In the early iteration steps, the die size and the area of the active (dissipating) chip-surface region were matched. With this modification, the first section of the heat-flow–path model was corrected (Figure 6). The match between the simulated thermal impedance and the measured thermal impedance was perfect up to a cumulative thermal resistance value of approximately 2.5 K/W.

After the die-attach layer’s properties were also modified (interface resistance set to the proper value by adjusting the thermal conductivity of the TIM1 material), the matching of the structure functions was perfect up to roughly 4.4 K/W (Figure 7).

The last step in the model calibration was to properly set the thermal resistance of the applied TIM2 layer (Figure 8). Thus, the model calibration was completed. The remaining difference was attributed to the modeling of the cold plate that was used.

This procedure can be done automatically using the CFD simulation tool, Simcenter Flotherm™ software command center, which uses data from the Simcenter T3STER transient thermal tester.

Just two examples were discussed on how structure functions can be used to analyze what is happening

thermally inside a semiconductor package, or any complex electronic system for that matter. Structure functions also can be used for other applications such as to obtain thermal measurements in multichip packages, to conduct in-situ testing of TIMs in a package under different environmental conditions, to provide data for reliability analysis using temperature and power cycling, and to test AC-driving LEDs.

These capabilities are available today in testing equipment by the Simcenter T3STER thermal transient analysis system, Simcenter TERALED™ software for analyzing LEDs, Simcenter DYNTIM for testing TIMs and the Simcenter POWERTESTER™ 1500A software for power cycling and testing packages in the lab or on the factory floor.

1.Thermal Management for LED Applications, Clemens J.M. Lasance and Andras Poppe, Eds., Springer, New York, 2014.

2. Lasance C J M, den Hertog D., Stehouwer P., “Creation and evaluation of compact models for thermal characterisation using dedicated optimisation software,” Proceedings of the 15th IEEE Semiconductor Thermal Measurement and Management Symposium (SEMITHERM’ 99), San Diego, USA, 9-11 March 1999, pp. 189-200.

3. Székely V., Bien T.V., “Fine structure of heat flow path in semiconductor devices: a measurement and identification method,” Solid-State Electronics 31(9):1363-1368 (1988).

4. Farkas G., Poppe A., Kollár E., Stehouwer P., “Dynamic Compact Models of Cooling Mounts for Fast Board Level Design,” Proceedings of the 19th IEEE Semiconductor Thermal Measurement and Management Symposium (SEMI-THERM’03). San Jose, 11-13 March 2003. pp. 255-262 (2003).

5. Szabó P., Poppe A., Rencz M., “Studies on the possibilities of in-line die attach characterization of semiconductor devices,” Proceedings of the 9th Electronics Packaging Technology Conference (EPTC’07), 10-12 December 2007, Singapore, pp. 779-784 (ISBN: 978-1-4244-1324-9).

6. Szabó P., Rencz M., Farkas G., Poppe, A. “Short time die attach characterization of LEDs for in-line testing application,” Proceedings of the 8th Electronics Packaging Technology Conference (EPTC’06): Volume One. Singapore, 6-8 December 2008, pp. 360-366 (ISBN: 1-4244-0664-1).

7. ASTM Standard D 5470, 2012, Standard Test Method for Thermal Transmission Properties of Thin Thermally Conductive Solid Electrical Insulation Materials. ASTM International, West Conshohocken, Pennsylvania, USA, 2012, DOI: 10.1520/D5470-12; www.astm.org.

8. http://www.mentor.com/products/mechanical/micred/dyntim/

9. Vass Várnai A., Bornoff R., Sárkány Z., Ress S., Rencz M., “Measurement Based Compact Thermal Model Creation – Accurate Approach to Neglect Inaccurate TIM Conductivity Data,” Proceedings of the 13th Electronics Packaging Technology Conference (EPCT’11), 7-9 December 2011, Singapore, pp. 67-72.

Structure functions are obtained by direct mathematical transformations from the measured or simulated thermal transient response functions of the system. They are a “graphical” representation of the RC model of a thermal system. In the case of essentially 1D heat-flow (such as longitudinal flow in a rod, or radial spreading in homogeneous material layers, or even in cylindrical or spherical spreading), structure functions can be considered as direct models of the thermal system.

Cauer-equivalent models of thermal impedances also can be generated from the time-constant spectra (for more details and theory behind this, see Thermal Management for LED Applications

The practical problem with such a network model is that 150-200 individual thermal resistance and thermal capacitance values cannot be interpreted. But by introducing two simple definitions, the Cauer-equivalent model can then be represented graphically. The cumulative thermal resistance is defined as the thermal resistance between the n-th element of the model network and the heat source (driving point)

and the cumulative thermal capacitance is

where Ri and Ci denote the element values of the i-th stage of the Cauer-type model network.

Starting from the driving point (the junction heat source), the sum of the partial thermal resistance and

partial thermal capacitance values for all subsequent heat-flow-path sections can be calculated.

Interpreting the cumulative thermal capacitance as a function of the cumulative thermal resistance, the

so-called cumulative structure function (CSF) is obtained:

It can then be proved that the derivative of C

The origin of the function corresponds to the junction. Because all thermal capacitance values are positive, it follows that the cumulative structure function should be a monotonically increasing function. The heat-conduction path ends in the ambient, which has an infinite heat sinking capacity; therefore, the cumulative thermal capacitance must tend to infinity. This means that the cumulative structure function should end with a singularity (at the location corresponding to the ambient). As a further consequence, the distance between the singularity and the origin is the junction-to-ambient thermal resistance (R

If there is heat flow through a small portion of material, two effects can be observed. The first effect is a temperature drop between two (isothermal) surfaces of the material (assuming an adiabatic condition at the other four faces of the cuboid).

If the material has thermal conductivity (λ) and power (P) flows through the surfaces a and b, they will have T

where the expression in the bracket on the right side is the R

The second effect arises because the same material slice can store thermal energy. If there is heat flow into the material, then in a short dt = t

The cumulative structure function is a graphic tool that has proven to be useful for analyzing physical structure. In low-gradient sections, it shows that a small amount of material with low capacitance is causing a large change in thermal resistance. These regions have low thermal conductivity and/or small cross-sectional area. Steep sections correspond to material regions of high thermal conductivity or large cross-sectional area. Sudden breaks of the slope belong to material or geometry changes. In such a way, thermal resistance and capacitance values, geometrical dimensions, heat-transfer coefficients and material parameters can be indirectly interpreted from cumulative structure functions.

In a real environment, the heat flow can have various shapes – longitudinal along a beam, radial in a board, or conical in a heatsink below a smaller package. In most cases, a “proper” slicing of the material can be made on the isothermal surfaces, perpendicular to the direction of the flow. These slices must be narrow, but not always of a very small cross-sectional area (Figure A1), and the surfaces are usually not planes.

It is sometimes easier to identify the interface between the sections using the derivative of the cumulative curve, known as the differential structure function (DSF). Here, peaks correspond to regions of high thermal conductivity, such as the chip or a heatsink, and valleys show regions of low thermal conductivity, such as the die attach or air. Interface surfaces are represented as inflexion points between peaks and valleys.

The differential structure function, derived from the previous two equations, also yields information on the cross-sectional area along the heat conduction path:

For example, consider the homogeneous rod with thermal boundary conditions shown in Figure A2. This rod can be considered as a series of infinitesimally small material sections. Consequently, the network model of its thermal impedance would also be a series connection of the single RC stages. With this slicing along the heat conduction path, a ladder of lateral thermal resistances is created between two thermal nodes and thermal capacitances between a node and the ambient.

Assuming homogeneity, the ratio of the elementary thermal capacitances and thermal resistances in the network model shown in Figure A2 would be constant. As shown in Figure A3, this means that the cumulative structure function of the rod would be a straight line – its slope is determined by the C

Assume that in a given section in the middle of the rod, the Cth/Rth ratio is halved. As Figure A4 shows, this results in a shallower middle section in the cumulative structure function (with the slope halved), and a corresponding dip in the differential structure function (being half the height of the sections either side of it).

Therefore, the differential structure function is a powerful tool for detecting slight changes in the cumulative structure-function curve, aiding the interpretation of the results.

Questions? Post a reply to this article or email Al Prosuk (albert.prosuk@siemens.com) for more information.

Simcenter T3STER: Non-destructive semiconductor thermal testing