Combined Microstructure and Heat Transfer Modeling of Carbon Nanotube Thermal Interface Materials1

2016 ◽  
Vol 138 (4) ◽  
Author(s):  
Sridhar Sadasivam ◽  
Stephen L. Hodson ◽  
Matthew R. Maschmann ◽  
Timothy S. Fisher

A microstructure-sensitive thermomechanical simulation framework is developed to predict the mechanical and heat transfer properties of vertically aligned CNT (VACNT) arrays used as thermal interface materials (TIMs). The model addresses the gap between atomistic thermal transport simulations of individual CNTs (carbon nanotubes) and experimental measurements of thermal resistance of CNT arrays at mesoscopic length scales. Energy minimization is performed using a bead–spring coarse-grain model to obtain the microstructure of the CNT array as a function of the applied load. The microstructures obtained from the coarse-grain simulations are used as inputs to a finite volume solver that solves one-dimensional and three-dimensional Fourier heat conduction in the CNTs and filler matrix, respectively. Predictions from the finite volume solver are fitted to experimental data on the total thermal resistance of CNT arrays to obtain an individual CNT thermal conductivity of 12 W m−1 K−1 and CNT–substrate contact conductance of 7 × 107 W m−2 K−1. The results also indicate that the thermal resistance of the CNT array shows a weak dependence on the CNT–CNT contact resistance. Embedding the CNT array in wax is found to reduce the total thermal resistance of the array by almost 50%, and the pressure dependence of thermal resistance nearly vanishes when a matrix material is introduced. Detailed microstructural information such as the topology of CNT–substrate contacts and the pressure dependence of CNT–opposing substrate contact area are also reported.

Author(s):  
Koji Nishi ◽  
Tomoyuki Hatakeyama ◽  
Shinji Nakagawa ◽  
Masaru Ishizuka

The thermal network method has a long history with thermal design of electronic equipment. In particular, a one-dimensional thermal network is useful to know the temperature and heat transfer rate along each heat transfer path. It also saves computation time and/or computation resources to obtain target temperature. However, unlike three-dimensional thermal simulation with fine pitch grids and a three-dimensional thermal network with sufficient numbers of nodes, a traditional one-dimensional thermal network cannot predict the temperature of a microprocessor silicon die hot spot with sufficient accuracy in a three-dimensional domain analysis. Therefore, this paper introduces a one-dimensional thermal network with average temperature nodes. Thermal resistance values need to be obtained to calculate target temperature in a thermal network. For this purpose, thermal resistance calculation methodology with simplified boundary conditions, which calculates thermal resistance values from an analytical solution, is also introduced in this paper. The effectiveness of the methodology is explored with a simple model of the microprocessor system. The calculated result by the methodology is compared to a three-dimensional heat conduction simulation result. It is found that the introduced technique matches the three-dimensional heat conduction simulation result well.


2001 ◽  
Vol 124 (1) ◽  
pp. 11-21 ◽  
Author(s):  
J. Cadafalch ◽  
C. D. Pe´rez-Segarra ◽  
R. Co`nsul ◽  
A. Oliva

This work presents a post-processing tool for the verification of steady-state fluid flow and heat transfer finite volume computations. It is based both on the generalized Richardson extrapolation and the Grid Convergence Index GCI. The observed order of accuracy and a error band where the grid independent solution is expected to be contained are estimated. The results corresponding to the following two and three-dimensional steady-state simulations are post-processed: a flow inside a cavity with moving top wall, an axisymmetric turbulent flow through a compressor valve, a premixed methane/air laminar flat flame on a perforated burner, and the heat transfer from an isothermal cylinder enclosed by a square duct. Discussion is carried out about the certainty of the estimators obtained with the post-processing procedure. They have been shown to be useful parameters in order to assess credibility and quality to the reported numerical solutions.


2015 ◽  
Vol 137 (3) ◽  
Author(s):  
Andrew J. McNamara ◽  
Yogendra Joshi ◽  
Zhuomin Zhang ◽  
Kyoung-sik Moon ◽  
Ziyin Lin ◽  
...  

Recently, much attention has been given to reducing the thermal resistance attributed to thermal interface materials (TIMs) in electronic devices, which contribute significantly to the overall package thermal resistance. Thermal transport measured experimentally through several vertically aligned carbon nanotube (VACNT) array TIMs anchored to copper and silicon substrates is considered. A steady-state infrared (IR) microscopy experimental setup was designed and utilized to measure the cross-plane total thermal resistance of VACNT TIMs. Overall thermal resistance for the anchored arrays ranged from 4 to 50 mm2 KW-1. These values are comparable to the best current TIMs used for microelectronic packaging. Furthermore, thermal stability after prolonged exposure to a high-temperature environment and thermal cycling tests shows limited deterioration for an array anchored using a silver-loaded thermal conductive adhesive (TCA).


Author(s):  
Sridhar Sadasivam ◽  
Stephen L. Hodson ◽  
Timothy S. Fisher

A growing interest has developed in the use of carbon nanotube (CNT) arrays as thermal interface materials (TIMs). However, theoretical modeling of CNT TIMs has largely been limited to semi-empirical methods without detailed consideration of array microstructure, primarily due to the inherent randomness of the microstructure and the computational complexity involved in full atomistic modeling of CNTs. In this work, we report combined thermo-mechanical simulation of CNT arrays with a coarse-grain approach for the mechanics modeling and a thermal network approach for the heat transfer modeling. Parametric studies on the effects of CNT height on the Young’s modulus and buckling load of CNT arrays are reported. The thermal network model is used to estimate the pressure dependence of diffusive and tip contact resistances of CNT arrays; the predictions are compared with thermal resistance measurements using the photoacoustic method. The resulting simulation framework enables a particularly rich and broad thermo-mechanical data set. Selected parametric variations are computed to assess the stress-strain behavior, effective conductivity within the CNT array, and aspects of the contact topologies of the CNT-substrate interface.


Author(s):  
E. Farsad ◽  
S. P. Abbasi ◽  
M. S. Zabihi

Performance of microchannel heatsink (MCHS) partially filled with foam is investigated numerically. The open cell copper foams have the porosity and pore density in the ranges of 60–90% and 60–100 PPI (pore per inch), respectively. The three-dimensional steady, laminar flow, and heat transfer governing equations are solved using finite volume method. The performance of microchannel heatsink is evaluated in terms of overall thermal resistance, pressure drop, and heat transfer coefficient and temperature distribution. It is found that the results of the surface temperature profile are in good agreement with numerical data. The results show the microchannel heatsink with insert foam appears to be good candidates as the next generation of cooling devices for high power electronic devices. The thermal resistance for all cases decreases with the decrease in porosity. The uniformity of temperature in this heatsink is enhanced compared the heatsink with no foam. The thermal resistance versus the pumping power is depicted, it is found that 80% is the optimal porosity for the foam at 60 PPI with a minimum thermal resistance 0.346 K/W. The results demonstrate the microchannel heatsink partially filled with foam is capable for removing heat generation 100 watt over an area of 9 × 10−6 m2 with the temperature of heat flux surface up to 59 °C.


2020 ◽  
Vol 14 (1) ◽  
pp. 6344-6361
Author(s):  
Pankaj Srivastava ◽  
Anupam Dewan

This paper presents the effects of microchannel shape with ribs and cavities on turbulent heat transfer. Three-dimensional conjugate heat transfer using the SST k-ω turbulence model has been investigated for four different microchannels, namely, rectangular, rectangular with ribs and cavities, convergent-divergent (CD) and convergent-divergent with Ribs and Cavities (CD-RC). The flow field, pressure and temperature distributions and friction factor are analyzed, and thermal resistance and average Nusselt number are compared. The thermal performance of the CD-RC microchannel is found to be better than that of other microchannels considered in terms of an average Nusselt number increased from 16% to 40%. Heat transfer increases due to a strong fluid mixing and periodic interruption of boundary-layer. It is observed that with an increase in Reynolds number (Re), the thermal resitance drops rapidly. The thermal resistance of the CD-RC microchannel is decreased by 30% than that of the rectangular microchannel for Re ranging from 2500 to 7000. However, such design of microchannel loses its heat transfer effectiveness due to a high pumping power at high values of Re.


2020 ◽  
Vol 10 (24) ◽  
pp. 9044
Author(s):  
Rafał Wyczółkowski ◽  
Marek Gała ◽  
Vazgen Bagdasaryan

During heat treatment of rectangular steel sections, a heated charge in the form of regularly arranged packages is placed in a furnace. The article presents a model of a complex heat transfer in such a package using the thermo-electric analogy. The model considers the following types of heat transfer: conduction in section walls, conduction and natural convection within gas, heat radiation between the walls of a section, as well as contact conduction between the adjacent sections. The results of our own experimental research were used for calculations of heat resistance applying to natural convection and contact conduction. We assumed that the material of sections was low-carbon steel and the gas was air. The result of the calculations of the presented model is total thermal resistance Rto. The calculations were performed for the temperature range 20–700 °C for four geometrical cases. Due to the variability of conditions for contact heat conduction, we assumed that total thermal resistance for a given charge is contained within a value range between Rto-min and Rto-max. We established that the value of Rto depends significantly on the section’s geometry. The larger the section sizes, the greater the changes of Rto. The minimal and maximal values of Rto for all packages were 0.0051 (m2·K)/W and 0.0238 (m2·K)/W, respectively. The correctness of model calculations was verified with the use of experimental data.


Sign in / Sign up

Export Citation Format

Share Document