CARBON NANOSTRUCTURES AS THERMAL INTERFACE MATERIALS:

PROCESSING AND PROPERTIES

Thesis

Submitted to

The School of Engineering of the

UNIVERSITY OF DAYTON

In Partial Fulfillment of the Requirements for

The Degree

Master of Science in Aerospace Engineering

By

Muhammad Omar Memon

Dayton, Ohio

May, 2011

CARBON NANOSTRUCTURES AS THERMAL INTERFACE MATERIALS:

PROCESSING AND PROPERTIES

Name: Memon, Muhammad Omar

APPROVED BY:

______Khalid Lafdi, Ph.D Kevin Hallinan, Ph.D Advisor Committee Chairman Committee Member Mechanical and Aerospace Engineering Mechanical and Aerospace Engineering

______Lawrance Flach, Ph.D Don Klosterman, Ph.D Advisor Committee Chairman Committee Member Chemical and Materials Engineering Chemical and Materials Engineering

______John G. Weber, Ph.D Tony Saliba, Ph.D Associate Dean Dean School of Engineering School of Engineering

ii

ABSTRACT

CARBON NANOSTRUCTURES AS THERMAL INTERFACE MATERIALS:

PROCESSING AND PROPERTIES

Name: Memon, Muhammad Omar University of Dayton

Adviser: Dr. Khalid Lafdi

The power density of electronic packages has substantially increased. The thermal interface resistance involves more than 50% of the total thermal resistance in current high-power packages. The portion of the thermal budget spent on interface resistance is growing because die-level power dissipation densities are projected to exceed 100 W/cm2 in near future. There is an urgent need for advanced thermal interface materials (TIMs) that would achieve order- of-magnitude improvement in performance.

Carbon nanotubes and nanofibers have received significant attention in the past because of its small diameter and high thermal conductivity. The present study is intended to overcome the shortcomings of commercially used thermal interface materials by introducing a compliant material which would conform to the mating surfaces and operate at higher temperatures.

iii

Thin film ―labeled buckypaper‖ of CNF based Materials was processed and optimized. An experimental setup was designed to test processed materials in terms of thermal impedance as a function of load and materials density, thickness and thermal conductivity. Results show that the thermal impedance decreased in conjunction with the increasing heat-treatment temperature of

CNFs. TIM using heat treated CNF showed a significant decrement of 54% in thermal impedance. Numerical simulations confirmed the validity of the experimental model. A parametric study was carried out which showed significant decrement in the thermal resistance with the decrease in TIM thickness.

A transient spike power was carried out using two conditions; uniform heat pulse of 24 Watts, and power spikes of 24-96 Watts. The results show that heat treated CNF was 12% more temperature resistant than direct contact with more than 50% enhancement in heat transport across it.

iv

ACKNOWLEDGEMENTS

I am very grateful to Allah (God) for blessing me throughout this project and my entire life. This accomplishment would not have been possible without

His help and desire.

It is my pleasure to thank those who made this research possible. I am heartily thankful to my advisor, Khalid Lafdi, for his guidance, encouragement and support from the initial to the final level of this research. I am grateful to him for being so kind and patient throughout this project. I would like to extend my sincere thanks to the committee members Dr. Kevin Hallinan, Dr. Lawrance

Flach and Dr. Don Klosterman, for evaluating this research.

I would like to thank my collaborators at the University of Dayton Research

Institute (UDRI), Dr. Lee, Alex and Matt Boehle for their tremendous support and guidance at each step of this project.

I owe my deepest gratitude to; my grandparents for their uncountable prayers for my health and success, my parents, for their incomparable love and support, my dear siblings Henna, Ali and Sana for being extremely patient and morally supportive, and last but not least, my immediate relatives for taking care of me at every single step throughout this endeavor. This work would not have been possible without the love and support of my entire family.

v

TABLE OF CONTENTS

ABSTRACT ...... iii

ACKNOWLEDGEMENTS ...... v

LIST OF FIGURES ...... viii

LIST OF TABLES ...... xii

Chapter I – Introduction ...... 1

Chapter II – Literature Review ...... 6

2.1 Introduction ...... 6

2.2 System Approach ...... 7

2.3 Materials Approach ...... 13

2.3.1 Thermal Interface Materials ...... 14

2.4 Thermal Resistance Measurement ...... 29

2.5 Proposed Methodology ...... 31

Chapter III –Material Fabrication and Characterization ...... 33

3.1 Materials Fabrication ...... 33

3.2 Validation of the Test Setup...... 52

vi

3.3 Experimental Characterization ...... 55

3.3.1 Thermal Conductivity calculation ...... 55

3.3.2 Thermal Impedance comparison...... 57

3.3.3 PVA/CNF Ratio ...... 60

3.3.3.1 PVA Effects ...... 66

3.3.3.2 CNF Effects ...... 68

3.4 Modeling ...... 70

3.5 Validation of the Model ...... 74

3.6 Numerical Simulation ...... 77

3.6.1 Parametric Study ...... 78

Chapter IV – Power Spike Analysis ...... 81

4.1 Experimental Setup ...... 86

4.2 Results...... 93

4.2.1 Uniform Power Pulse ...... 93

4.2.2 High Power Spikes ...... 93

4.3 Summary ...... 104

Chapter V – Conclusions ...... 106

Chapter VI – Recommendations ...... 111

BIBLIOGRAPHY ...... 113

APPENDIX ...... 130

vii

LIST OF FIGURES

Figure 2.1 Interface with both contact points and air gaps ...... 14

Figure 2.2 Nyquist plots of the thermal impedance of different interface materials

[77] ...... 23

Figure 2.3 Ratio of styrene to other hydrocarbon products as a function of time on stream [121] ...... 29

Figure 3.1 Bright-field image of pristine PS carbon nanofibers ...... 33

Figure 3.2 Models of various nanofiber configurations ...... 34

Figure 3.3 Carbon plane structure as a function of heat-treatment temperature

[134] ...... 35

Figure 3.4 Bright field micrograph of Dixie Cup carbon nanofiber structure ...... 36

Figure 3.5 High-resolution imaging of localized area of Dixie Cups carbon nanofiber structure ...... 37

Figure 3.6 Volume resistivity of as received carbon nanofibers at 108psi ...... 38

Figure 3.7 Fabrication process of buckypaper TIM ...... 39

Figure 3.8 a) Filtration process and b) Prepared sample ...... 40

Figure 3.9 Compliance buckypaper TIM ...... 40

Figure 3.10 Experimental Setup ...... 42

Figure 3.11 Control unit ...... 42

Figure 3.12 Schematic of the test setup with resistances ...... 43

viii

Figure 3.13 Global resistivity of the system ...... 45

Figure 3.14 Thermal resistances at each component ...... 46

Figure 3.15 Thermal linear fit as a function of thickness ...... 48

Figure 3.16 Calibration of the system ...... 49

Figure 3.17 Accumulation time of the setup ...... 50

Figure 3.18 Thermal impedance of Silver as a function of load ...... 52

Figure 3.19 Thickness calculations of samples by using crosshead position ..... 54

Figure 3.20 Thermal resistance of buckypaper TIMs as a function of thickness at various loads ...... 55

Figure 3.21 Plot showing increase in thermal conductivity of the sample as load increases ...... 56

Figure 3.22 Decrease in porosity of the sample as the load increases ...... 57

Figure 3.23 Comparison of thermal impedance of various heat treated CNF buckypapers with the dry contact ...... 59

Figure 3.24 Thermal resistance of BP1 as a function of thickness ...... 61

Figure 3.25 Thermal resistance of BP2 as a function of thickness ...... 61

Figure 3.26 Thermal resistance of BP3 as a function of thickness ...... 62

Figure 3.27 Thermal resistance of BP4 as a function of thickness ...... 62

Figure 3.28 Plot showing the thermal impedance of PVA/CNF ratio ...... 64

Figure 3.29 Thermal conductivity comparison of PVA to CNF ratio ...... 65

Figure 3.30 Heat transfer coefficient as a function of ratio of PVA to CNF ...... 65

Figure 3.31 Thermal impedance values as a function of PVA weight ...... 67

Figure 3.32 Thermal impedance of various TIMs as a function of load ...... 67

ix

Figure 3.33 Thermal impedance values as a function of CNF weight ...... 69

Figure 3.34 Thermal impedance of various Buckypaper TIMs as a function of load ...... 69

Figure 3.35 Schematic of the test model ...... 71

Figure 3.36 Validation of the test model ...... 75

Figure 3.37 Graph showing the thermal impedance comparison of Arctic Silver using various analysis techniques ...... 76

Figure 3.38 Contours of static temperature on the numerical model ...... 77

Figure 3.39 Static temperature distributions on the numerical model ...... 78

Figure 3.40 Thermal impedance of various TIMs as a function of thickness ...... 79

Figure 3.41 Thermal impedance of various TIMs as a function of load ...... 80

Figure 4.1 Original design of the experimental setup ...... 87

Figure 4.2 Schematic of the sample holder ...... 88

Figure 4.3 a) Experimental setup b) Enlarged view of sample holder ...... 90

Figure 4.4 Calibration of the six thermocouples...... 91

Figure 4.5 Direct contact analyses for uniform heat pulse ...... 95

Figure 4.6 Direct contact analyses for transient power spikes ...... 95

Figure 4.7 CNF-HHT buckypaper analyses for uniform heat pulse ...... 96

Figure 4.8 CNF-HHT buckypaper analyses for transient power spikes ...... 97

Figure 4.9 CNF-LHT buckypaper analyses for uniform heat pulse ...... 98

Figure 4.10 CNF-LHT buckypaper analyses for transient power spikes ...... 98

Figure 4.11 CNF-PS buckypaper analyses for uniform heat pulse ...... 99

Figure 4.12 CNF-PS buckypaper analyses for transient power spikes ...... 100

x

Figure 4.13 Temperature comparison of various heat treated CNF buckypapers and direct contact subjected to uniform heat pulse ...... 101

Figure 4.14 Temperature comparison of various heat treated CNF buckypapers and direct contact subjected to transient power spikes ...... 101

Figure 6.1 Adhesive forces of a single gecko foot ...... 112

Figure A.1 Cooling Plate ...... 130

Figure A.2 Aluminum Base Plate ...... 130

Figure A.3 Heater Plate ...... 131

Figure A.4 Plastic Insulator ...... 131

Figure A.5 1x1 inch Aluminum block ...... 132

Figure A.6 2x2 inch Aluminum block ...... 132

Figure A.7 Test Setup with 1x1 inch blocks installed ...... 133

Figure A.8 Test Setup with 2x2 inch blocks installed ...... 133

xi

LIST OF TABLES

Table 3.1 Description of tested carbon nanofibers ...... 37

Table 3.2 Configuration parameters of MTS Displacement Unit ...... 49

Table 3.3 Final error calculation of thermocouples ...... 51

Table 3.4 Validation of the setup ...... 52

Table 3.5 Thermal conductivity values at various loads ...... 56

Table 3.6 Thickness comparison of three types of CNF buckypapers ...... 58

Table 3.7 Thermal impedance comparison of direct contact and CNF TIMs ...... 58

Table 3.8 Thermal performance comparison of various TIMs ...... 59

Table 3.9 Thermal impedance ratio of PVA to CNF Buckypaper ...... 63

Table 3.10 Weight ratios of various buckypaper TIMs ...... 66

Table 3.11 Weight ratios of various buckypaper TIMs ...... 68

Table 3.12 Validation of the numerical model with experimental and published data for Arctic Silver 5 ...... 76

Table 4.1 Maximum temperature results for uniform pulse analyses ...... 102

Table 4.2 Temperature difference results for uniform pulse analyses ...... 102

Table 4.3 Maximum temperature results for power spikes analyses ...... 103

Table 4.4 Temperature difference results for power spikes analyses ...... 103

xii

Chapter I

Introduction

Thermal management is an important design consideration for number of engineered device components and packages. The need for integrating complex functions in a single circuit, together with a demand for thinner, light weight and most efficient products has been massively growing. High power and faster speeds are the major requirements of the global industry which result in excessive heating of the systems. The modern electronic systems are therefore characterized by increased density of circuits. As the circuit becomes denser, high power is required and hence high heat is generated throughout the system which should be dissipated. Rising amount of power has to be offset by efficient cooling in order for a system to be efficient. A major challenge in the field is, therefore, the ability to manage the heat without compromising on the performance of the system.

As component powers continue to grow, so do their cooling requirements.

One rule of thumb says that for every 10°C rise of the junction temperature, the failure rate doubles. Thus, there is an urgent need to remove heat from current electronic devices to the surrounding air stream. Most systems are designed to

1 minimize thermal resistance and maximize power dissipation. But with increased miniaturization of systems and increased circuit density, today‘s electronics are vulnerable and tend to generate substantial amount of heat. If the heat is not dissipated, the lifetime and reliability of the electronics will be at great risk. This is a problem that requires external thermal solutions, including heat sinks, fans, heat exchanger etc. Although these solutions are increasingly used to keep the temperature of the devices at a minimum, but materials also play a critical role.

For any system assembler, when two machined surfaces are in contact, there are only few contact points between surfaces. The insulating air gaps created by non contacting areas become a thermal barrier. Current technique to overcome this barrier is either by eliminating the air gaps through an increased surface wetability (using fluid based media) or to enhance the conduction process by using highly conductive material as a thermal interface material (TIM). TIMs are used to couple these devices and interfaces with lids and heat spreaders [1-4].

Yovanovich et al. [1] calculated that replacing air with grease can reduce the thermal resistance by a factor of five or so (depending on the surfaces and contact pressure).

Thermal interface materials are commonly in the form of paste () [5-7], polymers [8-9], flexible graphite [10-13], phase change materials

[14-15], low melting alloys [16-17] and nano-structured carbon materials [18-22].

The performance of a thermal interface material depends on its conformability, thermal conductivity and thickness [11]. An important property of any TIM is its thermal conductivity. Unfilled polymers have a thermal conductivity of about 0.1

2

W/m.K. All modern TIMs are composites containing fillers such as aluminum oxide, aluminum nitride, boron nitride, graphite and diamond powder, silver that increase thermal conductivity up to the 2 W/m-K. Unfortunately, high thermal conductivity alone is not enough to guarantee optimal system performance.

Thermal resistance obtained from a one-dimensional heat flow calculation should also be minimized. Other properties taken into account are interfacial thickness and applied forces. In addition to thermal performance, TIMs are selected based on their inherent thermal properties, ease of use in assembly and long-term stability (reliability).

Historically, greases typically yield an interface resistance of about 1 Kcm2/W.

They are used to eliminate microscopic air pockets. They are often characterized as ―messy‖ and difficult to apply due to their high viscosity. A more significant concern with application is reproducibly delivering the correct amount to achieve full coverage with the appropriately thin gap. The components require mechanical clamping pressure of about 300 kPa for optimal thermal performance. Elastomeric pads are the logical extension of greases but exhibit a higher thermal resistance range 1-3 Kcm2/W. They are supplied as die-cut performs in the precise shape needed for the application. Assembly is therefore very simple with these products but require high clamping pressures (~700 kPa) are needed to achieve an adequate interface. Thermal tapes were developed as a attachment method. They eliminate the need for external clamps but thermal performance is in the 1 – 4 Kcm2/W resistance range and is highly dependent on surface quality. Tapes also have very limited compliance and are

3 generally not suitable for current most electronic packages. However, phase change materials combine the thermal performance of grease with the convenience of an elastomer pad. Phase change materials are predominantly waxes that typically melt in the 50 – 80°C range. Essentially, they are effective conductors of heat both above and below the melt point. When operating above the melt point they are not effective as an adhesive and need mechanical support, therefore they are always used with a clamp applying pressure in the

300 kPa range. Performance levels are very close to grease, in the range 0.5 –

0.7 Kcm2/W. Solder is another TIM often overlooked. It represents the ultimate, a solid metal interface with resistance < 0.05 Kcm2/W. Despite the challenges of high temperature processing, solder is used as a thermal interface where no other viable option exists. The performance of a thermal interface material depends on its conformability, thermal conductivity and thickness. However, for the interfaces that are subjected to high temperature applications, have high roughness or not perfectly aligned, a compliant TIM must be needed. In contrast with the conventional TIMs, compliant materials conform to the mating surfaces, resist higher temperatures and do not leak out of the interface in complex applications.

Researchers are constantly seeking ―better‖ thermal interface products. This is usually expressed as a request for higher thermal conductivity and lower thermal impedance. But is that really possible? The current study attempts to use a nanotechnology approach to look for alternative solutions.

4

In the light of above mentioned issues in heat management, the present study is focused on minimizing thermal contact resistance using carbon nanofiber as thermal interface materials (TIM). The performance of these TIMs was explored in pulsed heat load conditions. Chapter 2 gives a comprehensive literature review on possible approaches to overcome thermal management issues. The thermal interface material fabrication process test results and numerical simulation of the test model in terms of a parametric study on thermal conductivity and thickness of the TIMs are described in chapter 3. Chapter 4 provides TIM behavior under transient spike powers. Finally, chapter 5 concludes the present study and its implications on the ongoing and future electronics applications.

5

Chapter II

Literature Review

2.1 Introduction

The need for integrating complex functions in a single circuit, together with a demand for thinner, light weight and most efficient products has been massively growing. The rapid increase in the study of controlling of matter at nano-scale is a driving force in consumer applications [1-2]. Heat dissipation has become one of the most critical problems that limit improvement in performance and reliability because no two surfaces can ever make a perfect contact due to machining limitations. Tiny air gaps always remain which provide poor transmission of heat through [3]. However within the contact points, at relatively high contact pressures, hot spots may occur at the contact interface. The rapid temperature rise caused by such intensive thermal energy dissipation may be sufficiently large to cause some damage to the system [4]. Improved functionality and efficient performance of the products require more heat that must be dissipated. Heat sinks like mechanical fans, coolers have traditionally been used to transfer heat to the ambient. This solution however has been a limiting factor in designing complex and most efficient systems.

6

The thermal contact resistance varies considerably depending on the interface geometry, the thermal and mechanical properties of the contacting materials and the interstitial fluid [23]. The contact peak criterion is critical to the contact peak characteristics which are required to evaluate the number and size of the contact points as well as air gaps. The thermal contact conductance is the sum of the conductance of several discrete spots that exist on the interface. The thermal contact resistance can hence be determined as a function of surface parameters, material properties and contact pressure [24]. Based on published studies, both system and material approaches were considered to overcome thermal management issues [6-12].

2.2 System Approach

The system approach seems very complex and limited because of material casting issues. For solids of high thermal conductivity, the contact resistance may be reduced by increasing the area of contact spots, accomplished either by increasing contact pressure which will ‗flatten‘ the peaks of the micro-roughness, and deflecting the mating surfaces to reduce any non- flatness, or reducing the roughness of the surfaces before the interface is formed by grinding the surfaces to remove non-flatness and buffing the surface to reduce micro-roughness [34]. Even by reducing the roughness through polishing the contact surfaces, it is impossible to obtain a perfect contact. Machining limitations and cost constraints are major aspects of this shortcoming. Adequate increase in the contact pressure is impractical due to the load requirements of the system.

7

Load constraints on electronic components and circuit boards make it unfeasible to use high contact pressure [35]. A finite distance below the heat source is required for the heat flux to become uniform. Similarly, the type and location of heat sinks also have significant effects on the overall performance of the device.

There have been a number of comprehensive reviews that explored various types of heat sinks designed for efficient heat dissipation [36-49]. Lelea [36] presented a geometric optimization of the micro-heat sink with straight circular micro-channels. The inlet cross-section of rectangular shape was positioned tangentially to the tube axis with four different geometries. The results of thermal performances of heat sink were presented in terms of the temperature distribution along the bottom wall at constant heat flux and temperature of 100

W/cm2 and 293 K respectively. The author observed that the lowest temperature was obtained for case 3 where the inlet channel width covered half of the tubes cross-section with minimum temperature of 313.2 K and ∆T of 10.17 K. The author reported that the results from tangential micro-heat sink showed better thermal performance than the conventional micro-heat sink with lateral inlet/outlet cross-section. Elshafei [37] performed experiments on natural convection heat transfer from circular pin fin heat sinks subject to the influence of its geometry, heat flux and orientation. The author found that the solid pin fin heat sink showed competitive performance for upward and sideward orientations. Higher heat transfer heat transfer coefficients were obtained in comparison to those of perforated/hollow pin fin ones in both arrangements. The augmentation factor was found to approximately 1.05-1.11. According to the author, the temperature

8 difference between the base plate and surrounding air, at the same heat input value, was found to be less for hollow/perforated pin fin heat sink than that for solid pin heat sink. Barba et al [38] evaluated the temperature distribution and thermal resistance of a polymeric heat sink with circular micro-channels. The bottom surface of the heat sink received a uniform heat flux, while the top surface was adiabatic. The experiments were carried out using two Nitrogen and Helium gases. The thermal conductivities and heat transfer coefficients were obtained to be 0.158 and 0.027 W/m-K and 3436 and 588 W/m2-K for helium and nitrogen respectively. The thermal resistance however was 5 ºC/W for Nitrogen and 3.7

ºC/W for Helium. Chiang et al. [39] used response surface methodology to identify the effects of design parameters of the pin-fin heat sink on the thermal performance. The height and diameter of pin-fin and the width of pitch between fins were considered as design parameters. The thermal resistance and pressure drop were calculated. Results show that the pin height and pin diameter were significant influential factors in minimizing the thermal resistance. Increase in these parameter values leads to the increase in the rate of removal heat capacity, hence thermal resistance value is decreased. The value of pressure drop was inversely proportional to the pin diameter and longitudinal and transverse pitch values. The authors also predicted the optimum values of thermal resistance and pressure drop by conducting confirmation experiments.

Khor et al [40] investigated the importance of the effects of thermal radiation and its view factor on the thermal performance of a straight-fin heat sink. Three different models were developed to explore these effects on the convection

9 coefficient as well as the fin performance of the heat sink. Results show that the average convection coefficient for the case neglecting thermal radiation was largest with 30% error, followed by that with thermal radiation including view factor, and that with thermal radiation excluding view factor as the lowest with

60% error. The fin effectiveness was overrated for the cases when thermal radiation and view factor were excluded with more than 40% error. The authors concluded that it is reasonable to exclude thermal radiation under certain circumstances but to consider thermal radiation without incorporating the view factor could be resulting in more critical and detrimental errors. Kim et al. [41] suggested closed form correlations for thermal optimization of vertical plate-fin heat sinks under natural convection in a fully-developed-flow regime. Analytical solutions for velocity and temperature distributions for high channel aspect ratios, high conductivity ratios, and low Rayleigh numbers were presented using the volume averaging approach. The analytical solutions explained the explicit correlation for optimal fin thickness and optimal channel width, which minimize thermal resistance for given height, width, and length of heat sink. The authors reported that these correlations show that the optimal fin thickness is a function of height, the solid conductivity, and the fluid conductivity only and is independent of the Rayleigh number, the viscosity of the fluid, and the length of the heat sink.

Huang et al [42] carried out experiments on natural convection heat transfer from square pin fin heat sinks subject to the influence of orientation. The authors tested a flat plate and seven square pin fin heat sinks with various arrangements under controlled environment. The results indicated that the upward and

10 sideward facing orientations were of comparable magnitude and showed competitive nature. Another significant factor highlighted by the authors was the porosity of the heat sink which had secondary effect on the performance of the pin fin. The optimal porosity of the heat sink was found to be around 83% for upward arrangement and 91% for the sideward arrangement. The augmentation factor (heat transfer of a heat sink relative to the base plate) was found to be around 1.1-2.5 for the upward arrangement and between 0.8–1.8 for the sideward arrangement. Jouhara and Axcell [43] described the thermal conditions within a heat sink with rectangular fins under forced convection cooling. The authors show a numerical study on key parameters of heat transfer that change with axial distance. The heat transfer coefficient and pin efficiency rapidly changed near leading edges of cooling fins. The authors reported that despite the rapid changes in these parameters, good engineering accuracy for heat sink performance is attainable using analytical methods which incorporate average values of heat transfer coefficient and fin efficiency. Lorenzini et al. [44] presented a modular heat sink consisting of elements with wavy fin profile which can be pressed together to construct the component. The temperature distribution and global thermal resistance were measured under steady state conditions. The experiments were carried out on both uniform and non uniform heat flux. The temperatures were measured by varying the location of the dissipating heaters, the power, and the number of spacers (one and five) between fan and sink. The thermal resistance is almost unaffected by the value of the ambient temperature [45]. Results show that the device performed better

11 when the heat source is distributed over the whole surface of the heat sink and at either side of it. Although the thermal resistances with one and five spaces were almost equal, 0.087 and 0.089 K/W respectively, the overall temperature distribution was better for single spacer compared to several spaces [44]. Ismail et al. [45] used four different types of heat sinks; Pentiums III and IV, AMD Athlon and Duron heat sinks; in order to analyze their performance. The simulation and experimental investigations were made at different Reynolds numbers. The experiments were carried out by using the air chamber with nozzle at different

Reynolds numbers. The authors reported that the total surface area and fin spacing significantly affects the heat sink performance. To obtain maximum performance i.e. minimum thermal impedance, the surface should be higher with suitable fin density. Although the surface area of all heat sinks was equal, the

Pentium IV heat sink showed highest performance in heat dissipation with thermal impedance of 0.4 ºC/W compared to other heat sinks. This was achieved due to larger fin spacing. Higher pressure drop across closer fins produced poor flows to pass through the fins and hence reduces the performance of other heat sinks. Wang [47] investigated the thermal performance of heat sinks with one and two pairs of embedded heat pipes. The embedded heat pipes transfer the total heat capacity from the heat source to the base plate and disperse heat into the ambient. The author found that that two and four heat pipes embedded in the base plate carry 36% and 48% of the total dissipated heat respectively. It was also reported that when the total heating power of the heat sink with two embedded heat pipes was 140 W, the total thermal resistance reached the

12 minimum value of 0.27 °C/W, while for the heat sink with four embedded heat pipes, when the total heating power was between 40 W and 240 W, the total thermal resistance was calculated to be 0.24 °C/W.

In the light of above cited literature, the performance of heat sinks largely depends on its geometry, the placement of the heat source and the amount of power dissipated. By the same token, heat sink‘s performance varies with smooth and wavy channels because of the difference in pitch.

2.3 Materials Approach

An important way to alleviate the overheating problem is to improve the thermal contact between the heat source and heat sink by using thermal interface materials (TIMs) on the interface of contacting surfaces [48-50]. Any interstitial substance that fills the gap between contacting surfaces, and whose thermal conductivity exceeds that of air, will decrease the contact resistance [51].

The characteristics of a TIM would be high thermal conductivity, minimal thickness, would not leak out of the interface, maintains performance indefinitely,

Non-toxic and manufacturing friendly (easy to apply and remove) [52]. A visual representation of interface micro-voids filled with air and TIM is shown in figure

2.1.

13

Figure 2.1 Interface with both contact points and air gaps

2.3.1 Thermal Interface Materials

Thermal interface materials (TIMs) are used to eliminate air gaps between contacting surfaces in order to maximize thermal transport. These materials are highly conductive and provide efficient heat path by conforming to the mating surfaces. Gwinn and Webb [52] presented a classification of high performance

TIMs used to minimize the thermal resistance between the CPU and the heat sink, and discussed their advantages and disadvantages.

Conventional Thermal Interface materials

The most commonly used method to reduce thermal contact resistance is to introduce a thin layer of oil or grease at the interface of contacting surfaces.

Mineral oil is one of the most common thermal fluids which is highly conformable and spreadable, but has low thermal conductivity. Thermal greases are composed of thermally conductive fillers dispersed in silicone or a hydrocarbon oil to form a paste. Thermal greases require no post-dispense processing and

14 have higher effective thermal conductivity compared to other classes of materials

[53]. Greases have been used very successfully in combination with various packaging form-factors and have shown excellent performance. Corrugated copper substrate with thermal grease is described by Tzeng [54] which has high lateral conductivity and distributes the heat in the lateral direction more evenly.

This would however be of value only if the heat flux were non-uniform in the interface region. Salerno et al. [55] measured thermal contact conductance of pressed contacts using Apiezon-N grease at the interface. The measurements were taken over the temperature range of 1.6-6.0 K and applied forces from 22 to

670 N. The authors reported that the addition of thermal grease between the contact surfaces resulted in an improvement of thermal conductance over uncoated surfaces to approximately a factor of 3. He [56] investigated thermal transport properties of three thermal greases (G9, G1 and S78) used as thermal interface materials in microelectronic packaging. The specific heat as a function of temperature, thermal conductivity and thermal diffusivity of these materials were measured. The author reported that among all three greases, G1 with more than 90 wt.% of highly conductive fillers, showed a thermal conductivity of

5 W/m K. Lallemand and Sartre [57] measured thermal contact resistance using several commercial greases and foils as a function of heat transfer rate and material thickness. The contact resistance decreased as the heat flux or the thickness of the material decreased. According to the authors, the highest dimensionless contact conductance factors (E) were achieved with greases (3 <

E < 6). Phase change material coated foils exhibit E-values ranging from 2.5 to

15

3.5. Graphite or metallic foils have E-values lower than 2 and for silicone foils E is significantly reduced to less than 1. Viswanath et al. [53] presented a detailed stress analysis on a number of interface materials including thermal greases. The authors found various phenomenon associated with thermal greases under stress like 1) power cycling (loss of material due to a relative motion between the die and the base of a heat sink. The authors indicate that the assembly between

0 and 100 ºC over 7500 cycles results in a four to six fold increase in thermal resistance compared to a negligible increase in resistance for a 0 to 80 ºC exposure over 2500 power cycles. 2) thermal bake (under high-temperature bake, the formulation chemistries utilized in typical thermal greases and cause separation of the polymer and filler matrix which results in poor wetability at the interfaces, 3) Mechanical Shock and Vibration: Data collected on mechanical shock and vibration of heat-sink masses between 200-250 grams indicate that the retention of the heat sink to the processor is critical to prevent mechanical damage to the die surface. Maguire et al. [5] undertook a comprehensive study understand the thermal conditions within high power, radio frequency signal amplifiers. Experiments were carried out to determine the interface thermal resistance of a thermal grease compound and comparison of results with other alternatives. The authors reported that at maximum power dissipation, interface temperature (which accounts for 23.4% of total temperature rise) was reduced by up to 20 ºC through the use of high conductivity thermal grease. The optimum interstitial thickness depends on the thermal conductivity of the grease, which can be improved by adding metallic particles to the silicone greases [58]. Chung

16

[59] explored thermal contact conductance values for a number of greases and the influence of various additives to those. The author found that PEG polymer based paste gives much higher contact conductance (11.0×104 W/m2 °C) than silicone paste (3.08×104 W/m2 °C), due to its relatively low viscosity. The addition of Li salt (1.5 wt.%) to PEG raised the conductance to 12.3×104 W/m2 °C. Further addition of water increased the conductance to 16.0×104 W/m2 °C and decreases the viscosity. Further addition of BN particles (18.0 vol.%) enhanced the conductance to 18.9×104 W/m2 °C. Sodium silicate based pastes were found to give higher thermal contact conductance than polymer based pastes and oils, due to their fluidity and the consequent greater conformability [7]. A gel-like TIM has been developed by Thermoset of Lord Chemical Products [60], named as gelease, that has thermal performance similar to high thermal conductivity greases, but does not leak out of the interface. Carbon black filler [9, 61], is relatively more conformable to the contacting surfaces because of its high thermal conductivity. One added advantage of carbon is its porosity which allows penetration of the organic material into a carbon particle and hence marks it superior over conventional interface material. In a recent study [62], Hu and

Chung used polyol-ester based carbon black pastes for graphite coating and/or penetration to increase the thermal contact conductance of flexible graphite between copper surfaces. The authors reported that the paste penetration by up to an effective thickness of 5 μm increased the conductance by up to 350%, 98% and 36% for thicknesses of 50, 130 and 300 μm, respectively. Whereas the paste coating up to 10 μm increases the conductance by up to 200%, 120% and 65%

17 for thicknesses of 50, 130 and 300 μm, respectively. Since the paste penetration effectively reduced the thermal resistivity compared to coating, the highest thermal contact conductance was found to be 1.4 × 105 W/m2 K by paste- penetrated flexible graphite of thickness 26 μm. Conversely, thermal greases are not manufacturing friendly as are messy, difficult to apply and remove, and can dry out with time resulting in increased thermal resistance [52]. Excess grease that flows out of the interface enables contamination and electrical shorts which limits its overall thermal performance.

In order to overcome the limitations of thermal greases, thermal pads were introduced which are composed of silicone gel in combination with a thermal medium. These pads have an inherent tack that aids placement during assembly, are more stable than phase change products, and have higher operating temperatures [63]. Due to this, thermal pads exhibit slightly better performance than greases. The Bergquist Inc. [64], one of the primary thermal interface material producers, has introduced a thermal pad named ―GAP-PAD‖. These pads can suit different application sectors hence the thermal conductivity could vary from 0.6 to 5 W/mK and thickness from 0.254 to 6.35 mm. Chomerics [65] produced a thermal pad named THERM-A-GAP which can offer thermal conductivity of up to 6.5 W/mK and thickness of 1 mm. Use of thermal pads could save time by speeding up the assembly process but they are too expensive to use.

Phase Change Materials (PCMs), because of its high storage capacity; have received considerable attention in recent years [66-68]. PCMs can be

18 classified as organic and inorganic compounds [69]. Organic PCMs can be further classified as paraffin and non-paraffin organics [75]. Abhat [69] reviewed various PCMs for low temperature latent heat storage in the range of 0-120 ºC.

The freezing and melting behaviors of these materials were investigated using various thermal analysis techniques. Himran et al. [71] reviewed a number of paraffin waxes for use as energy storage materials. Pal and Joshi [72] used paraffin wax, n-triacontane in a sealed heat sink to thermally control rising temperatures of the heat source. Clarksean et al. [73] performed a numerical analysis of a finned surface surrounded by a paraffin to investigate high heat flux levels. Liu and Chung [17] carried out a comparative study on phase change behavior of organic and inorganic PCMs (paraffin and microcrystalline respectively) with melting temperature close to room temperature. They found that paraffin wax was potentially good thermal interface material because of the negative super-cooling of −7°C, large heat of fusion (up to 142 J/g) and excellent thermal cycling stability. Paraffin also exhibited clear endothermic and exothermic DSC peaks which were not clear for microcrystalline. Incongruent melting and decomposition behavior was also seen in microcrystalline with high supercooling of 8°C or more and instable thermal cycling. Most PCMs with high energy storage density have an unacceptably low thermal conductivity. Therefore the performance of the PCMs can be enhanced by mixing the PCM with polymers and with particulate fillers of high thermal conductivity. Baker et al.

[138] used two nonmetallic PCMs (octadecane and pentacosane) for heat storage with metal fins to enhance the effective thermal conductivity. The results

19 show high energy storage rates in the package and in PCMs as well as peak operating temperatures. Vesligaj and Amon [140] used an epoxy polymer as a phase change material for passive thermal control during time-varying workloads on electronics. They indicated that the operational performance of the electronics improved when such a passive thermal storage device was used. Krishnan and

Garimella [143] carried out a transient analysis of the phase change process inside a rectangular enclosure for electronics cooling applications with pulsed power dissipation. The experimental investigation included the melting of a pure

PCM, n-eicosane, inside a rectangular aluminum container with multiple discrete heat sources mounted on one side. A copper heat spreader was used in conjunction with the heat sources. The influence of changes in the frequency of pulses, heat source location and aspect ratio of the containment volume on the thermal performance of the PCM unit was studied. The performance of the PCM system was quantified by studying maximum temperature observed in the container, cumulative heat energy absorbed heat loss through the container walls and instantaneous heat absorbed. It was reported by the authors that the heat source location and container aspect ratio played a very significant role in the performance of the PCM system with pulsed heat input. Krishnan and Garimella

[72] also performed a transient thermal analysis to investigate thermal control of power semiconductors using phase change materials, and to compare the performance of this approach to that of copper heat sinks. The authors concluded that a significant suppression of junction temperatures was achieved by the use of PCMs when compared with copper heat sinks. Pal and Joshi [73]

20 investigated transient transport characteristics of melting of PCMs from a uniformly dissipating source, with other walls adiabatic. It was observed that melting of PCM due to buoyancy driven natural convection, considerably dampened the effects of transient heat flux. It was shown that by estimating the heat transfer rate and molten volume fractions using the dimensionless numbers the amount of PCM and the surface area can be determined for thermal management of uniformly dissipating heat sources with various power dissipations. Alawadhi and Amon [74] carried out an experimental and numerical investigation of a TC unit for portable electronic devices. The experimental model was a prototype of a Technical Information Assistance, a wearable computer where the TC unit is embedded in the electronic device. The experimental model contained a TC unit, which is an aluminum enclosure, with aluminum foam impregnated with Eicosane PCM. The authors concluded that for the varying power operations, the TC unit highly reduced the heat source temperature fluctuations and its average temperature using PCM-foam system. Shaikh and

Lafdi [75] used paraffin wax (with carbon nanotube additives) for composite thermal control (TC) system for thermal protection of electronics against transient pulsed power heat loads. The samples were subjected to both uniform and varying power conditions. The PCM composites due to greater energy storage and heat dissipation rate resulted in lower values for maximum junction temperatures for both uniform and varying power loads. The results also showed improvement in overall thermal performance of TC unit for protection against pulsed heat loads. Nurmawati et al. [76] used PCMs as potential thermal

21 interface materials by using two Ni-plated Cu heat spreaders. Thermal resistances were measured at initial thicknesses of 0.20–0.22 mm at various pressures, temperatures and powers. Since thermal resistance is directly proportional to the thickness of the sample, the minimum resistance of 0.13 ºC/W was obtained at given thickness and 0.35 MPa. Gwinn and Webb [52] reported the thermal interface resistance values for several PCMs vary from 0.14 to

0.58 K cm2/W at 10 psi. The authors also indicated the improvement in PCM performance at higher contact pressures i.e. the thermal resistance decreased by

60% (from 0.14 to 0.058 K cm2/W) when the contact pressure was increased to

80 psi. In contrast to the above literature, there are a number of significant disadvantages of using PCMs. These fillers require moderate compression between surfaces which causes the TIM to flow. The poor thermal conductivity of

PCMs at high heat flux levels is one of the major limitations towards their overall thermal performance.

Selection of an interface material not only influences the overall thermal resistance (Rth=Zth(0)) but the entire shape of the Nyquist plot as well. Although the has almost the same total thermal resistance as compared to the case without any interface material, it is clearly observed from Figure 2.2 that the middle frequency parts of the two plots do not coincide. The middle arc of the

Nyquist plot is associated to the interface material. Hence a thermal resistance

Rth interface can be associated to each of the interface materials [77].

22

Figure 2.2 - Nyquist plots of the thermal impedance of different interface

materials [77]

Carbon Nanomaterial based Thermal Interface materials

Carbon nanotubes (CNTs) and Carbon nanofibers (CNFs) are new forms of carbon with outstanding electrical and thermal properties, high specific tensile strength and modulus [78-84]. These nanostructures are produced by many techniques including chemical vapor deposition (CVD) process [85-91]. Wang et al. [92] fabricated CNT-reinforced buckypapers by dispersing CNTs into water with the aid of epoxy resin. Multiple layers of these buckypapers with high CNT loading (up to 39 wt.%) were used to produce nanocomposites. The storage moduli of these nanocomposites were 429% higher compared to neat resin modulus. However, the major limitation of this material lies in its rigidness. Zaeri et al. [93] developed a mechanical model to determine the elastic properties of buckypaper composites. Various parameters were investigated including Young‘s modulus, shear modulus and porosity of the buckypaper composite. The Young‘s

23 modulus of buckypaper composite was 132.2 MPa compared to the predicted value of 148.6 MPa with the porosity of 67.5%. Hong and Tai [94] synthesized

SWCNTs and MWCNTs to form freestanding buckypapers. The thermal conductivities at room temperature of SWCNT and MWCNT buckypapers were found to be 2.24 W/mK, and 0.36 W/mK, respectively. Xu and Fisher [95] synthesized CNT arrays directly on silicon wafers using plasma-enhanced CVD to enhance thermal contact conductance of the system. The system was based on calorimeter principle with one-dimensional reference bars and carried out in a high-vacuum environment with radiation shielding. Results show that dry CNT arrays produce a minimum thermal interface resistance of 19.8 mm2 K/W, while the combination of a CNT array and a phase change material produces a minimum resistance of 5.2 mm2 K/W. Cross et al. [96] grew CNTs on Si substrate using CVD process followed by metallization using Ti/Au. The coated CNTs were bonded to metallized substrates at 220 ºC. Measurements were carried out by photoacoustic technique which showed that at no applied pressure, the thermal resistance for bonded CNT films of 30 μm in length was 1.7 mm2 KW−1 where as that of CNTs up to 130 μm in length was 10 mm2 KW−1. It was hence suggested by the authors [96] that the interface resistance at the two bonded interfaces in the CNT samples dominates the overall resistance of the structure. Cola et al.

[97] developed a contact resistance model to describe heat transfer across CNT array interfaces. The model revealed that the thermal resistances at single-CNT contacts dominate the thermal transport across CNT array interfaces. It was suggested that the total thermal resistance of CNT array interfaces can be

24 reduced by optimizing the array conformability and density such that the true contact area established in the interface is maximized. Wang et al. [98] fabricated a TIM by synthesizing aligned CNTs on both sides of a thin copper foil. The

Laser flash and hot disk methods were used to measure the thermal conductivity of the CNT-TIMs. Results showed that a thicker copper foil substrate or CNT layer led to a lower overall thermal resistance of 8.78 mm2 K/W. It was suggested that an enhancement in thermal conductivity of more than 290% could be obtained under an applied contact pressure of 0.01 MPa, as compared with two copper plates in direct contact. Shaikh et al. [1] fabricated a thin low density aligned CNT as TIM to explore the effect of thermal contact resistance in aluminum and graphite pieces. Compared to the direct contact, the thermal conductivity of CNT-TIM with aluminum surfaces was enhanced by 385% with an

80% decrease in thermal resistance. While with graphite surfaces, the thermal conductivity was increased by 362% with a 78% decrease in thermal resistance.

All these data show the superiority of carbon based TIMs over conventional materials. High fabrication cost is however a significant limitation of CNT based materials.

Wang in [98] synthesized aligned carbon nanotubes (CNT) as a thermal interface material on both sides of a thin copper foil to measure the thermal conductivity of CNT-TIMs. It was seen that aligned CNT layer led to a lower overall thermal resistance. An enhancement in thermal conductivity of nearly

300% was obtained under certain contact pressures. Results also show that a thicker aligned CNT film produced even lower thermal resistance but there is a

25 trade-off between thickness and thermal resistance which should be considered to get best results. The effect of interface on the overall conductivity of composite is very small for short nanotubes, whereas interface has a significant effect on the overall thermal conductivity of the composite for long nanotubes [99-102].

Hence the interface thickness has got small effect on the thermal performance of a material. It has been verified that the interface material between two contacting surfaces can have a major influence on the thermal impedance and not only on the thermal resistance. It was observed that the thermal impedance plots depend not only on the applied conditions but also on the interface material. Material thermal performance will be affected by surface finish, contact pressure, and application properties.

Conversely, carbon nanofibers (CNFs), due to their high tensile strength, modulus, and relatively low cost, are drawing significant attention for their potential applications in nano-scale polymer reinforcement [103-106]. They are synthesized from pyrolysis of hydrocarbons or carbon monoxide in the gaseous state, in the presence of a catalyst [107-109]. There are a number of studies that show the composites formulated using vapor grown carbon nanofibers

(VGCNFs) [67-72]. These CNFs distinguish themselves from other types of nano fibers, such as polyacylonitrile or mesophase pitch-based carbon fiber, in its method of production, physical properties and structure [110]. Mahanta et al.

[111] measured the thermal conductivity of vapor grown carbon nanofiber mats at different heat treatment temperatures and increasing volume fractions. The authors found that the in-plane thermal conductivity of the mats varied from 12

26

W/m-K to 157 W/m-K for volume fractions of 0.067 and 0.462 respectively.

Thermal conductivities through the thickness were measured to be 0.428 W/m-K and 0.711 W/m-K. The authors reported that thermal conductivities through the thickness increased by the order of magnitude when CNF mats were heat treated to temperatures above 3000ºC. Thermal conductivities of individual CNFs at temperatures of 1100 ºC and 3000 ºC were calculated to be close to 1400 W/m-K and 1600 W/m-K respectively. Patton et al. [112] measured thermal conductivities of CNF composites in epoxy of up to 0.8 W/m-K for 40 vol.%. This increment over the neat resin value of 0.26 W/m-K was attributed to the complexity of thermal energy transfer within fibers. Lafdi and Matzek [113] demonstrated a significant increase in thermal conductivity of CNF composite from 0.2 W/m-K for epoxy resin to 2.8 W/m-K for a 20 wt%. Zhou et al. [114] infused various percentages of CNFs to epoxy matrix through sonic cavitations at room temperature. It was seen that the storage modulus steadily increased with an increasing weight percent of the fiber. The addition of 3 wt% of CNF yielded a

65% increase in the storage modulus. A stable CNF layer was catalytically grown on Ni foam by decomposing ethylene [115]. The presence of two distinct carbon layers was found; an apparently dense layer ‗C-layer‘ at the carbon–Ni interface and a CNF layer on top of that. The nickel surface and the attached carbon layer had similar morphological features which provided strong adhesion of the C-layer to Ni which result in mechanical stability.

These properties can however be improved by tuning CNFs as a function of heat treatment temperatures. Tibbetts et al. [116] found that when CNFs are

27 heat treated at temperatures above 1500 °C, a significant rearrangement of the core morphology is occurred [117]. These CNFs have special structure with large inner and outer diameters, and a crystal imperfection. Besides the main product of CNFs, a small amount of amorphous carbon and metal particles can also be observed in as prepared materials. The performance of the carbon nanofibers

CNFs were related to their degree of graphitic order. The degree of graphitization is determined by the temperature and the preservation period. When the temperature is up to 1873 K, the crystallite increases resulting in more conductive fiber [118]. Endo et al. [119] reported that there was a sharp decrease in the interlayer spacing from 0.342 to 0.337 nm after heat treatment to temperatures above 2100 °C. Heat treatment at graphitizing temperatures (>

2400 °C) is an effective technique for purifying and annealing the graphitic structure [120]. Parrott and Zeitler in [121] performed a catalytic activity test on the three CNF samples which resulted in the formation of a variety of hydrocarbon products, with styrene and benzene accounting for in excess of

96% of these in all cases. The ratio of styrene to other hydrocarbon products is shown in Figure 2.3.

28

Figure 2.3 Ratio of styrene to other hydrocarbon products as a function of time

on stream [121]

2.4 Thermal Resistance Measurement

Thermal interface resistance and thermal conductivity measurements for thermally-enhanced interface materials are presented by a number of researchers in the published articles. Bolger [33] performed steady state thermal interface resistance measurements for single and multiple layers of epoxy tape adhesives with diamond, silver, aluminum and alumina fillers at high pressures between polished aluminum surfaces. Using these data, Bolger [33] calculated the effective value of thermal conductivity k for each of the adhesives and developed a correlation to predict the volume fraction of the filler required value of thermal conductivity. Mirmira et al. [34] performed thermal contact resistance measurements using a steady-state test for a variety of adhesive materials,

29 including epoxies, cements and silicone. Certain thickness of each material was considered and experiments were carried out as a function of overall interfacial resistance of the system. Kilik et al. [35] presented thermal conductivity data for a variety of copper and aluminum filled epoxy adhesives. The authors reported that the thermal conductivity of adhesives can be measured effectively using the unsteady method based on conduction into a high conductivity body through a low conductivity surface film of adhesive. The thermal conductivity of adhesive decreased with the increasing temperature however, it increased substantially upon the addition of copper and aluminum powders [35].

According to the above cited literature [33-35], these methods provide much faster results than a steady-state test but do not give any relevant data for thermal resistance between a thermal interface material (TIM) and the interfacing surface, which is as significant as the bulk conduction resistance. In this study however, a steady state method will be presented for determining the bulk thermal conductivity and thermal contact resistance across the interface using thermal impedance measurements. These tests will be performed with a hot and cold plate technique based on the calorimeter principle [122]. The TIM sheet was sandwiched between a hot and a cold plate and the heat flow is measured [123].

The thermal resistance is then calculated as the temperature difference between the hot and cold plate divided by the heat flow. The hot and cold plates are perfectly isothermal and they make a perfect contact with the interface material under test. Zhu and Wang [30] stated that the thermal interface resistance (TIR) can be reduced significantly by exerting the pressure or spreading an adhesive

30 between two contacting surfaces. Thermal conductivity is an important property of composites application for electronic packaging, thermal insulation, heat spreader, etc. [124]. Voller and Tirovic [125] presented some methods for decreasing the thermal resistance across a bolted joint, such as using thin aluminum gasket at the interface. Cheng and Madhusudana [126] investigated the contact heat transfer of fin–tube exchangers where zinc, tin, silver, and gold were chosen as plating materials. The minimum TIR was obtained when the tube was coated with tin. This indicates that, although the thermal conductivity is important, the softness of the plating material also plays an important role in enhancing the thermal conductance of the interface.

2.5 Proposed Methodology

To achieve high thermal conductivity, CNFs at various heat treatment temperatures (PR24-PS, PR24-LHT, and PR24-HHT) were used [133]. These

CNFs have special structure with large diameters (less than 100nm) and some crystal imperfections. The performance of the CNFs will be related to their degree of graphitizability. This study is focused on the fabrication and characterization of CNF buckypapers. The analysis is performed by varying the ratio of CNF and binder polyvinyl alcohol (PVA) and, heat treatment temperatures of CNFs. These altering parameters should lead to a buckypaper with specific density (i.e. thickness) and thermal conductivity.

There are numerous synthetic polymers that are used in fabrication of composite materials. Polymer based composites reinforced with a small percentage of strong fillers can significantly improve the mechanical, thermal,

31 and barrier properties of the pure polymer matrix [127]. Moreover, these improvements are achieved through conventional processing techniques without any detrimental effects on process ability, appearance, density, and ageing performance of the matrix. Many researchers synthesized micro-scale fillers with epoxy resin [128-130]. Zhou and Pervin [114] infused various percentages of

CNFs to epoxy matrix through sonic cavitation at room temperature. It was seen that the storage modulus steadily increased with an increasing fiber weight percent. The addition of 3 wt% of CNF yielded a 65% increase of the storage modulus. Another most common water-soluble polymer, Polyvinyl Alcohol (PVA), was used to make composite fibers which contained large fraction of nanotubes

[131]. Composite nanotube–PVA fibers exhibit an extremely high toughness, superior to that of any other materials [132-133]. They are electrically and thermally conductive and could therefore be the basis for lightweight, strong and conducting composites.

32

Chapter III

Material Fabrication and Characterization

3.1 Materials Fabrication

The as-received nanofibers consist of a variety of carbon configurations like helical, straight, nested, carbon blacks with narrow distribution of diameters and lengths (Figure 3.1).

Figure 3.1 Bright-field image of pristine PS carbon nanofibers.

33

The five major carbon constituents represented as carbon nanofibers are shown in Figure 3.2. Carbon blacks, sometimes recognized as soot, are nanometric spheres of agglomerated carbon. Helical carbon nanofibers, or nanocoils, consist of carbon nanofibers that have a configuration similar to that of a DNA strand. Straight carbon nanofibers exist as a series of coaxial carbon cylinders surrounding a central hollow tube. The bamboo species are similar to that of the straight carbon nanofibers, except that they are segmented along their length. Nested carbon nanofibers have an orientation similar to that of a set of stacked Dixie Cups with a hollow core and are also referred to as fishbone type carbon nanofibers.

Figure 3.2 Models of various nanofiber configurations

Carbon at low temperatures exhibits local molecular ordering. As a carbon sample is heat treated, the increase in temperature results in the aromatic

34 molecules stacking into a columnar structure. Further heat treatment causes these columns to join together forming a distorted, wavy structure. Heat treatment temperatures of 2500°C or higher compress the distorted graphene layers forming an aligned structure. Graphitic carbon will attain the minimum interlayer spacing in the graphite order between graphene layers (Figure 3.3).

Heat treatment of pristine carbon nanofibers to 3000°C results in a straightening of the graphene layers. The minimum interlayer spacing was attained for the

HHT nanofibers.

Figure 3.3 Carbon plane structure as a function of heat-treatment temperature [134]

As shown in a TEM micrograph (Figure 3.4) the layers within the Dixie

Cup carbon nanofiber coalesced following heat treatment resulting in localized ordering of the graphene layers and continuous planes. At this magnification the

35 inclination angle of each cup is apparent. The use of grey-scale allowed the walls of the nanofibers to appear dark due to the high electronic density in this region.

Figure 3.4 Bright field micrograph of Dixie Cup carbon nanofiber structure

High resolution imaging of the graphene layers (Figure 3.5) indicated the alignment of the layers without disclination defects. There is no change in the inclination angle to the central core axis. The edge of any pair of graphene layers is rounded to encapsulate the exposed edge of the carbon planes, allowing the exposed graphene planes to attain a level of maximum structural stability.

Each of the pristine carbon nanofibers was tested for electrical volume resistivity measurement. The description of each of the nanofibers is provided in

Table 3.1.

36

Figure 3.5 High-resolution imaging of localized area of Dixie Cups carbon nanofiber structure

Table 3.1 Description of tested carbon nanofibers

ASI Sample Description

AG Grade As-grown nanofiber (processing temperature 1100°C) PS Grade Pyrolytically stripped nanofiber (processing temperature 1100°C) LHT Grade Nanofiber after heat treatment temperature up to 1500°C HHT Grade Nanofiber after heat treatment temperature to 3000°C

The results of the volume resistivity testing are presented in Figure 3.6.

There is a significant drop in the volume resistivity between the AG and PS nanofibers and between the PS and LHT nanofibers. However, the change in

37 the temperature difference between the AG, PS, and LHT is only 400°C. This indicated an anomaly in the observed decrease in resistivity. Further heat treatment of the nanofibers does not significantly reduce the volume resistivity. If the volume resistivity were only a function of heat treatment temperature, the values should result in a relatively linear decrease with a sharp decrease after

2500°C. There is likely an additional factor in the determination of volume resistivity which resides either in the surface contact or deposition of a dielectric material.

108PSI

0.014

0.012

0.010

0.008

0.006

0.004

0.002

Vol. Resistivity (ohm-cm / %Vol) 0.000 As-Grown PR-19-PS PR-19-LHT PR-19-HHT Effect of Heat treatment Temperature

Figure 3.6 Volume resistivity of as received carbon nanofibers at 108psi

A mixture of carbon nanofibers (CNF) and water soluble binder polyvinyl alcohol (PVA) was subjected to ultrasonic dispersion process. A wide range of

38 parameters can be varied in this process to optimize the thermal properties of the thermal interface material (TIM). In the present study however, the parameters were kept constant to achieve best results in terms of thickness and conformability of the TIM. In the first step, 0.8 g of PVA was mixed with 400 ml of water. Secondly, 1.2 g of CNFs was added to the mixture before being subjected to ultrasonic dispersion process for four hours. The solution was then filtered using latex filter Millipore (0.7um). Figure 3.8 shows the filtration process and the filtered sample. The sample was then heat treated at 400 ºC for 15 minutes to evaporate residual PVA and water. The final sample obtained was a CNF compliant material (buckypaper) with less than 1% impurities (Figure 3.9). Figure

3.7 shows the step wise buckypaper fabrication process.

Figure 3.7 Fabrication process of buckypaper TIM

39

Figure 3.8 a) Filtration process b) Prepared sample

Figure 3.9 Compliance buckypaper TIM

The characterization tool used in this study is primarily based on the ASTM

D5470 test method for measuring thermal resistance. Major improvements have been made to the base design (Appendix A). It measures the thermal impedance and thermal conductivity under a range of clamping forces and thicknesses, for thermal interface materials (TIMs).

The construction of the experimental setup is based on a steady-state heat flow along the sample (TIM) placed between two aluminum blocks, one attached to the heating and another to the cooling plate, as shown in Figure 3.10. Each block is one square inch in size and has two thermocouples inserted into the heat

40 transport path. Power is applied to generate heat in the heating block. The material under test is placed between the upper block attached to the ‗hot surface‘ and lower block attached to the ‗cold surface‘. The lower block is designed utilizing a liquid cold plate which has a regulated supply of coolant of known temperature. Load is applied on top of the heating plate by the MTS

Electromechanical displacement unit. Data acquisition device reads the temperatures from four thermocouples when thermal equilibrium is established, the data is then transferred to the Excel spread sheet to perform the thermal interface resistance calculations.

The control unit (Figure 3.11) consists of a power controller, a temperature controller and a data acquisition device. A 120 V ~ 60 Hz dimmer is used as a power controller and VFL series PID as a temperature controller. NI-USB 9162 data acquisition device is used to read temperature reading from four thermocouples. Electrical power is applied to the controller via a power switch.

The MTS Electromechanical displacement unit was programmed and configured using Test Works 4.0. This includes pre-specimen, specimen and post-specimen programming of the crosshead movement. Appendix A shows the detailed parameters of each step in the program.

Prepared test sample is applied to the lower aluminum block. All tests are carried out at loads of 45, 90, 135 and 180 kg. The heater at the upper block was set to 80 ºC and the chilled water bath at the lower block was set to 10 ºC. The heat generated from the heater flows down thru the metering block, thru the TIM, thru the lower metering block and into the liquid cold plate controlled by the

41 chilled water bath. The free convection and radiation terms were neglected since the test rig is located in a controlled atmosphere room. At steady state, four temperature readings are taken with the help of the data acquisition. The temperature difference is taken to calculate the thermal impedance at each load.

The extension in the displacement cell measures the thickness of the TIM under test which is used to calculate the effective thermal conductivity of the TIM. The total heat loss from the system is estimated to be less than 3%.

Figure 3.10 Experimental Setup Figure 3.11 Control unit

42

Figure 3.12 Schematic of the test setup with resistances

Figure 3.12 represents total resistances of the two contacting surfaces as well as those at the interface. The determination of the heat flux through a thermal interface material is equal to the electrical power loss. Heat either transfers through the TIM contacts or the air gaps. The remaining air conductivity and the contact conductivity (between TIM and the aluminum block) could be added to achieve global thermal conductivity hglobal,Al as shown in Eq. 1.

43

Where;

2 hglobal,Al = Global thermal contact conductivity (W/m -K);

2 hAl-TIM = Thermal contact conductivity between Al and the TIM. (W/m -K);

2 hAl-Air = Thermal contact conductivity between Al and Air (W/m -K);

TIM = Thermal conductivity through the TIM sheet (W/m-K);

eAir = Diameter of the air gap (mm);

S= Sample area (mm2);

S1 = Upper Aluminum block/TIM contact surface (mm2);

S2 = Lower Aluminum surface/TIM contact surface (mm2).

The global resistivity of the system is represented in figure 3.13,

44

RAl

R int Al/Air

RAir Rint Al/Sample

Rint Air/Sample

RSample

Rint Sample/Air

R Rint Sample/Al Air

Rint Air/Al

RAl

Figure 3.13 Global resistivity of the system

45

With the addition of a TIM at the interface, the resistance at each sub- interface is added to the global resistivity of the system (Figure 3.13). The final schematic of an equivalent thermal circuit with a TIM inserted at the interface

(Figure 3.14)

Upper Aluminum block T4

RAl1= e1 / λAl T3

Tint1 Interface

Rint1 = 1/hAl-global

TIM

e = thickness Rthrough.TIM = e / λTIM

Interface

Rint2 = 1/ hglobal

Tint2

Lower Aluminum block T2

RAl2 = e2 / λAl. T1

Figure 3.14 Thermal resistances at each component

The global heat flow is measured by the temperature difference between two thermocouples of the lower aluminum block (physically known test).

46

And S = 1 In² (contact surface), hence;

– (Eq.3)

To determine the thermal resistivity of the TIM, the temperatures must be calculated between the interfaces Tint 1 and Tint2, determined by the four thermocouples:

Hence, the thermal resistivity of the TIM is defined by:

The global thermal impedance at the interface can be calculated as a function of thermal parameters:

47

This law of the thermal resistivity function of the thickness (e) is determined by a numerical linear fit as shown in Figure 3.15.

0.7

0.6 Slope coefficient

0.5

0.4

0.3 Y-Intercept 0.2

0.1

Thermal Impedance (K/W) ImpedanceThermal 0 0 0.02 0.04 0.06 0.08 0.1 Thickness (mm)

Figure 3.15 Thermal linear fit as a function of thickness

After individually calibrating each ‗T‘ type thermocouple using chilled water bath, each thermocouple reading is stored by using the LabView software. In order to determine the steady state condition, temperature curves are plotted as a function of time (Figure 3.16)

48

Figure 3.16 Calibration of the system

TestWorks 4.1, MTS software, was programmed to control the

Electromagnetic Load Cell, was calibrated by using the ‗Proportional Integral

Derivate‘ (PID) parameters. It was very important to maintain a good proportion between two aluminum blocks in order for crosshead to balance the load for a given time before returning to the initial position (Figure 3.17). Table 1 shows the values of Kp, Ki and Kd, at which the crosshead was stabilized. The maximum load to be applied was set to 180 kg (400 lbf) so as not to damage the load cell and work pieces.

Table 3.2 Configuration parameters of MTS Displacement Unit

PID Parameters

Kp 0.000100

Kl 0.000010 Max Integral (in/min) 2.362200 Derivative Interval 2.000000 Maximum Speed (in/min) 23.62200000

49

Figure 3.17 Accumulation time of the setup

After every test, the aluminum blocks surfaces are cleaned up by Acetone to remove any residual left. It is very essential to maintain same surface roughness during all tests to accurately compare the results. The effect of the heat transfer coefficient h(Al-Sample) is directly related to the roughness and the residual material.

The thermal resistivity calculation (temperature, thermocouple position and area):

50

Error ( ) =

Where:

yT1-yT2 = distance between two thermocouples (Error – 0.02 inch),

T1, T2, T3, T4 = Thermocouple temperatures (Error – 0.05K),

S = Contact area (mm2),

Rexp = Experimental thermal resistivity.

Table 3.3 Final error calculation of thermocouples

y(T1)-y(T2) 0.500 Δy 0.02 T1 27.29 ΔT1 0.05

T2 31.95 ΔT2 0.05

T3 43.64 ΔT3 0.05

T4 47.87 ΔT4 0.05

AL 6.02 S 1.00 ERROR 0.029

51

3.2 Validation of the Test Setup

First step in the process was to validate the test setup with the published results. In order to validate the designed setup, various dry tests were carried out to observe the functionality of the setup. As the next step, Arctic Silver thermal interface material was analyzed at different loads and thicknesses to compare the obtained results with the published ones (Figure 3.18).

Table 3.4 Validation of the setup

Arctic Silver Thermal Interface Material

Load Thickness Thermal Impedance (kg) (mm) (K/W)

Published 135 0.025 0.0136

Experimental 135 0.025 0.0137

Thermal Impedance - Arctic Silver 0.024 Published 0.022 Experimental 0.02

0.018

0.016

0.014

0.012 Thermal Impedance (K/W) ImpedanceThermal 0.01 0 50 100 150 200 250 Load (kg)

Figure 3.18 Thermal impedance of Arctic Silver as a function of load

52

The behaviour of thermal interface materials cannot be analyzed without taking into account the change in thickness of the sample. In order to determine the change in TIM thickness at different loads, the change in cross position is monitored. During the process, the crosshead position changes due to two reasons, the density of the TIM (the conformability) and due to thermal expansion of aluminium surfaces. At high temperatures, the thermal expansion of aluminium surfaces was monitored by carrying out dry tests (without any TIM). The initial and final positions of the upper aluminium block were determined at various loads. In order to analyze buckypaper thermal interface materials under various loads, the sample thickness e can be determined by;

(Eq)

Where;

Ccp = Current crosshead position

Icp = Initial crosshead position

53

2.6 45 kg 90 kg 2.55 135 kg y = -0.00053740x + 2.61086576 180 kg R² = 0.99961953 2.5

y = -0.0004888229x + 2.5403334025 R² = 0.9997716259 2.45

y = -0.000546381x + 2.477263591 R² = 0.999980846 2.4

Sample thickness Current position 2.35 y = -0.00056826x + 2.39990861 Crosshead position (mm) position Crosshead R² = 0.99957014 Initial position 2.3

2.25 0 50 100 150 200 250 Sum Temperatures (ºC)

Figure 3.19 Thickness calculations of samples by using crosshead position

In figure 3.19, four initial crosshead position lines (45, 90, 134, and 180) are shown during the direct contact calibration. Initially the residual values are approximately 1. When the TIM is inserted into the interface, the current crosshead position is 2.374 mm at 45 kg. The sum of four temperature readings is 150 ºC. Therefore the initial crosshead position for this temperature is 2.319 mm. The thickness of the sample is 0.055 inches (2.374 -2.319). The error for thickness calculation is less than 5%.

54

3.3 Experimental Characterization

3.3.1 Thermal conductivity calculation

An experimental characterization of heat conduction between two contacting surfaces has been developed which enables the measurement of thermal impedance as a function of load and thickness. Four buckypaper samples with different thicknesses were tested in order to find the influence of those on global thermal resistance. It was seen that as the load increases, both slope coefficient

A and Y-intercept B (Eq.7) decrease (Figure 3.20). Thus, thermal conductivity through the TIM (λTIM) and at the interface (hglobal,Al) increase with the load.

CNF Buckypaper TIM 0.9 45 kg 0.8 90 kg 135 kg 0.7 180 kg R² = 0.9991 0.6

0.5 R² = 0.9994

0.4 R² = 0.9997

0.3 R² = 0.9999

0.2 Thermal Impedance (K/W) ImpedanceThermal 0.1

0 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 Thickness (mm)

Figure 3.20 Thermal resistance of buckypaper TIMs as a function of thickness at various loads

55

Table 3.5 Thermal conductivity values at various loads

Thermal Conductivity Applied Load (kg)

45 90 135 180

Through the TIM (W/m-K) λTIM 0.71 0.83 0.91 1.02

4 4 4 4 At the Interface (W/m²-K) hglobal,Al 3.99x10 7.74x10 8.41x10 9.09x10

Thermal conductivity λTIM increases with the load due to significant decrease in TIM density. The hTIM-Al also increases due to better conformability of the TIM at higher loads. Figure 3.21 shows the graphical representation of the

TIM in terms of λTIM and hglobal,Al.

1.05 2.00E+05

K) K)

1.80E+05 - - 1

0.95 1.60E+05 1.40E+05 0.9 1.20E+05 0.85 1.00E+05 0.8 8.00E+04 0.75 6.00E+04 0.7 λTIM 4.00E+04 0.65

hglobal,Al 2.00E+04

Thermal Conductivity (W/m Conductivity Thermal Thermal Conductivity (W/m² Conductivity Thermal 0.6 0.00E+00 0 50 100 150 200 Load (kg)

Figure 3.21 Plot showing increase in thermal conductivity of the sample as load increases

56

The initial porosity of the buckypaper was approximately 0.4 but at higher loads, because of compression, carbon fibers tend to get closer allowing better heat transport. It is well known that the thermal conductivity of air is lower than that of carbon nanofibers and contact points increase with the pressure, the thermal conductance through the thermal interface material (λTIM) increases.

Figure 3.22 shows the decrease in porosity with increasing load.

Air gaps LOAD

Compressed Sample Initial Sample

The porosity decreases, hence the thermal conductivity λTIM increases

Figure 3.22 Decrease in porosity of the sample as the load increases

3.3.2 Thermal Impedance comparison

The performance of buckypaper TIM with different heat treated CNFs was investigated in comparison with the direct contact to analyze the advantage of using this compliant material. For this purpose, the TIM of approximately 0.004 inch thickness was selected. It was seen that the thermal impedance values with high heat treated CNF (HHT-CNF) buckypaper were significantly better than low

57 heat treated (LHT) and PS (pyrolytically stripped). For the buckypaper with HHT-

CNF as compared to the direct contact between two surfaces at maximum load, the thermal impedance was reduced by 54% from 0.151 K/W to 0.069 K/W

(Figure 3.23). Since the thickness of the TIM decreases with increasing load, the final bond line thickness attained was 0.075 mm at applied load of 180 kg (Table

3.7).

Table 3.6 Thickness comparison of three types of CNF buckypapers

Load Thickness (mm) (kg) CNF-HHT CNF-LHT CNF-PS 45 0.119 0.117 0.118 90 0.089 0.097 0.092 135 0.084 0.089 0.087 180 0.082 0.082 0.082

Table 3.7 Thermal impedance comparison of direct contact and CNF TIMs

Load Thermal Impedance (K/W) (kg) Direct Contact CNF-HHT CNF-LHT CNF-PS 45 0.391 0.156 1.4662 1.7242 90 0.264 0.105 1.0371 1.2501 135 0.212 0.083 0.8176 0.9756 180 0.151 0.069 0.7235 0.8332

58

2 BP (CNF-HHT) 1.8 BP (CNF-LHT) 1.6 BP (CNF-PS) 1.4 Dry Contact 1.2 1 0.8 0.6 0.4

0.2 Thermal Impedance Impedance Thermal (K/W) 0 0 50 100 150 200 Load (kg)

Figure 3.23 Comparison of thermal impedance of various heat treated CNF buckypapers with the dry contact

Table 3.8 – Thermal Performance comparison of various TIMs

Thermal Impedance Product (K/W) Enhancement Dry Contact 0.151 Reference Aremco 640 0.071 52.98% Shin-Etsu MicroSi® G751 0.098 35.10% Dow Corning® TC-5121 0.096 36.42% Omegatherm 201 0.081 46.36% PowerFilm 51 0.474 -213.91% Thermal Grease TIC-1000A 0.072 52.32% HighFlow 565U 0.221 -46.36% CNF-PS 0.833 -451.66% CNF-LHT 0.723 -378.81% CNF-HTT 0.069 54.30%

59

Table 3.8 shows thermal performance of various commercially used TIMs in comparison with that of CNF based buckypapers. Most conventional TIMs (e.g. thermal greases) are liquid based and messy in installation. They also require flawless application to ensure uniform distribution at the interface. The applied pressure can cause the TIM to leak out of the interface. In contrast, the buckypaper TIM is dry in nature and provides competitive thermal performance.

3.3.3 PVA/CNF Ratio

Four samples of CNF-HTT buckypaper were prepared to analyze the ratio of binder material PVA to the CNF-HHT. By varying the PVA weight and keeping the CNF weight constant (0.5g), different thermal impedance behaviors were seen. For each composition, two samples with different thicknesses were analyzed to find the thermal impedance as a function of thickness. All compositions were finally compared at a sample thickness of 0.1 millimeter.

60

BP1 (1g PVA,0.5g CNF-HTT) 2.5 45 kg 90 kg 2 135 kg 180 kg y = 3.902113x + 0.205368

1.5

y = 3.323697x + 0.126485 1 y = 2.991193x + 0.104678 y = 2.636377x + 0.104728

0.5 Thermal Impedance (K/W) ImpedanceThermal 0 0 0.1 0.2 0.3 0.4 0.5 Thickness (mm)

Figure 3.24 Thermal resistance of BP1 as a function of thickness

BP2 (2g PVA,0.5g CNF-HTT) 1.4 45 kg 1.2 90 kg 135 kg 1 180 kg y = 3.629070x + 0.114869

0.8 y = 2.936704x + 0.091285 0.6 y = 2.727238x + 0.066560 0.4 y = 2.555381x + 0.059062

0.2

Thermal Impedance (K/W) ImpedanceThermal 0 0 0.1 0.2 0.3 0.4 Thickness(mm)

Figure 3.25 Thermal resistance of BP2 as a function of thickness 61

BP3 (3g PVA, 0.5g CNF-HTT) 1.4 45 kg 1.2 90 kg y = 2.846760x + 0.266185 1

0.8 y = 2.413033x + 0.185199

0.6 y = 2.342858x + 0.123765 y = 2.128743x + 0.103730 0.4

0.2 Thermal Impedance (K/W) ImpedanceThermal 0 0 0.1 0.2 0.3 0.4 Thickness (mm)

Figure 3.26 Thermal resistance of BP3 as a function of thickness

BP4 (4g PVA, 0.5g CNF-HTT) 1.4 45 kg 1.2 90 kg 135 kg 1 180 kg y = 2.494422x + 0.317886 0.8 y = 2.217727x + 0.189347 0.6 y = 1.981286x + 0.163762 0.4 y = 1.791196x + 0.138824

0.2 Thermal Impedance (K/W) ImpedanceThermal 0 0 0.1 0.2 0.3 0.4 Thickness (mm)

Figure 3.27 Thermal resistance of BP4 as a function of thickness

62

From above figures (3.24-3.27), it was noticed that two samples with same

PVA to CNF ratio but different global weights like (BP1 and BP2) show identical thermal behavior. It is hence important to define the samples by its PVA to CNF ratio and not its global weight. However, the increase in the weight percent of

PVA with constant CNF weight has minimum effect on the thermal resistance at the interface.

Hence the thermal impedance graph in terms of PVA to CNF ratio was analyzed (Figure 3.28) and it was noticed that at low pressure, as the ratio between two materials increased, the thermal impedance increased. However, the thermal impedance remains constant at high temperatures, regardless of the ratio between PVA and CNF.

Table 3.9 Thermal impedance ratio of PVA to CNF Buckypaper

Thermal Impedance (K/W) PVA/CNF CNF 45 kg 90 kg 135 kg 180 kg Ratio Buckypaper 2 BP1 0.443 0.427 0.375 0.343 4 BP2 0.411 0.356 0.313 0.29 6 BP3 0.523 0.403 0.335 0.296 8 BP4 0.543 0.389 0.342 0.3

63

PVA/HTT Comparison 0.6 45 kg 0.55 90 kg 135 kg 0.5 180 kg

0.45

0.4

0.35

0.3

Thermal Impedance (K/W) Impedance Thermal 0.25

0.2 0 1 2 3 4 5 6 7 8 9 PVA/HTT Ratio

Figure 3.28 Plot showing the thermal impedance of PVA/CNF ratio

Similarly, the ratio of PVA to CNF is plotted in terms of thermal conductivity through the TIM thickness (Figure 3.29) and overall heat transfer coefficient

(Figure 3.30).

64

0.6 45 kg 90 kg 0.55

K) 135 kg - 180 kg 0.5

0.45

0.4

0.35

0.3

Thermal Conductivity (W/m ConductivityThermal 0.25

0.2 0 2 4 6 8 10 PVA/HTT Ratio

Figure 3.29 Thermal conductivity comparison of PVA to CNF ratio

40 45 kg 90 kg 35 135 kg 180 kg 30

25

20

15

10 Heat Transfer Coefficient Transfer Heat

5

0 0 2 4 6 8 10 PVA/HTT Ratio

Figure 3.30 Heat transfer coefficient as a function of ratio of PVA to CNF

65

3.3.3.1 PVA Effects

In order to understand the effects of water soluble binder material

Polyvinyl Alcohol (PVA), seven buckypaper samples were fabricated by varying the weight of PVA by 1 gram while keeping the CNF weight constant to 0.5 grams (Table 3.10). The thermal impedance was calculated for each sample in terms of a) PVA weight (Figure 3.31) and b) applied load (Figure 3.32).

Table 3.10 Weight ratios of various buckypaper TIMs

Material Composition CNF (g) PVA (g) BP1 0.5 1.0 BP2 0.5 2.0 BP3 0.5 3.0 BP4 0.5 4.0 BP5 0.5 5.0 BP6 0.5 6.0 BP7 0.5 7.0

66

PVA Effects 1.4 45 kg 1.2 90 kg 135 kg 180 kg 1

0.8

0.6

0.4

0.2 Thermal Impedance (K/W) ImpedanceThermal 0 0 2 4 6 8 PVA Weight (g)

Figure 3.31 Thermal impedance values as a function of PVA weight

PVA Load Effects 1.4 BP1 1.2 BP2 BP3 BP4 1 BP5 BP6 0.8

0.6

0.4

0.2 Thermal Impedance (K/W) ImpedanceThermal 0 0 50 100 150 200 Load (kg)

Figure 3.32 Thermal impedance of various TIMs as a function of load

67

3.3.3.2 CNF Effects

To study the effects of carbon nanofibers (CNF) on the conformability of the TIM, six buckypaper samples were created by varying the CNF weight by 0.5 grams while keeping PVA weight constant to 3 grams (Table 3.11). The thermal impedance was calculated for each sample in terms of a) CNF weight (Figure

3.33) and b) applied load (Figure 3.34).

Table 3.11 Weight ratios of various buckypaper TIMs

Material Composition CNF (g) PVA (g) BP1 0.5 3.0 BP2 1.0 3.0 BP3 1.5 3.0 BP4 2.0 3.0 BP5 2.5 3.0 BP6 3.0 3.0

68

CNF-HTT Effects 0.7 45 kg 0.6 90 kg 135 kg 180 kg 0.5

0.4

0.3

0.2

0.1 Thermal Impedance (K/W) ImpedanceThermal 0 0 0.5 1 1.5 2 2.5 HTT weight (g)

Figure 3.33 Thermal impedance values as a function of CNF weight

CNF-HTT Load Effects 0.9 BP1 0.8 BP2 BP3 0.7 BP4 0.6 BP5 BP6 0.5 0.4 0.3 0.2

0.1 Thermal Impedance (K/W) ImpedanceThermal 0 0 50 100 150 200 Load (kg)

Figure 3.34 Thermal impedance of various buckypaper TIMs as a function of load

69

The main purpose of the test setup was to characterize the TIMs in terms of effect of the ratio of CNF and PVA (binder material) on the global thermal resistance, and thermal impedance as a function of load and thickness. TIMs with high heat treated fibers minimized the global resistance of the system by

54% compared to the direct contact.

3.4 Modeling

The Thermal Impedance test setup was modeled on Fluent CFD to simulate the applied test conditions. Numerical analysis was performed to compare the results with those from the experimental characterization. A numerical approach offers the advantage that a lot of parameters, including thickness, dimensions and thermal parameters can be varied. The interface resistance can also be tuned for a wide range of materials and applications. This chapter shows a detail description of geometrical design, modeling, simulation and analysis.

In order to conduct the numerical study for heat dissipation between a heat source and heat sink, two aluminum blocks were observed in contact in three ways;

1) Dry Contact – To understand the ideal contact condition with zero

resistance.

2) Using Arctic Silver as a TIM – To validate the test model with

experimental and published results.

3) Using CNF buckypaper TIM – To understand the effect of CNF on the

global thermal resistance of the system.

70

In the model shown in Figure 3.35, the path of heat removal involves conduction across the interface of the upper aluminum block, through a TIM, into the lower aluminum block and then convection to the environment. For the system temperatures close to the ambient, natural convection and radiation terms were neglected.

Heater

Average thickness “e”

T4 Upper aluminum block

T3

Int1 TIM

Int2 T2 Lower Aluminum block

T1

Air gaps Cooler

Figure 3.35 Schematic of the test model

The numerical analysis is based on the exact simulation of the test stand. The simulation is carried out based on the governing equations for the overall system.

The heat flow is calculated as;

(1)

71

Where;

K is the thermal conductivity and

A is the surface area

As shown in Fig 3b, the interface resistance is composed of three resistances in series. The total interface resistance is expressed as;

(2)

Where;

R1 is the contact resistance between TIM and the lower block,

Rcond is the conduction resistance across the thickness of the TIM,

R2 is the contact resistance between TIM and the upper block, and

(3)

Where; t is the thickness of the TIM,

The total R can be expressed as R = R1+R2. One thing to notice is that resistances R1 and R2 are most likely to be same since each of the two contacting surfaces are similar (aluminum). Hence the equation reads,

72

(4)

In order to carry out the model simulation in Fluent for the solid surfaces, the energy equation is given by;

(5)

Where;

is the density, h is the enthalpy,

K is thermal conductivity,

T is temperature and

H is the volumetric heat source.

For the heating plate, the input heat flux is given by;

(6)

And the temperature boundary condition for the cooling plate is given by;

(7)

73

Where;

= Thermal conductivity of the solid,

= Local solid temperature, and

= distance between wall surface and the solid cell center.

3.5 Validation of the model

As a primary step, the numerical model was validated by comparing the results of the present numerical method using thermal paste TIM (Arctic Silver 5) with experimental results available in the literature. The values of total thermal resistance using this model showed practical agreement with the experimental and published results with an error of less than 6%. Table 1 compares the contact thermal resistance of Arctic Silver 5 at a given load and thickness.

74

Figure 3.36 Validation of the test model

Figure 3.36 shows the heat transport from the heating plate (50 W/m2) through the TIM (AS-5) and to the cooling plate (288 K). Four temperature reading at locations (T1, T2, T3 and T4) were obtained to calculate the total thermal impedance as a function of thickness of the TIM.

75

Table 3.12 Validation of the numerical model with experimental and published data for Arctic Silver 5

Arctic Silver Thermal Interface Material

Load Thickness Thermal Impedance (kg) (mm) (K/W)

Published 135 0.0254 0.0136

Experimental 135 0.0254 0.0137

Numerical - 0.0254 0.0133

Thermal Impedance - Arctic Silver 0.024 Published 0.022 Experimental Numerical 0.02

0.018

0.016

0.014

Thermal Impedance (K/W) ImpedanceThermal 0.012

0.01 0 50 100 150 200 250 Load (kg)

Figure 3.37 Graph showing the thermal impedance comparison of Arctic Silver

using various analysis techniques

76

Figure 3.37 shows the Arctic Silver 5 performance comparison of the values calculated from the numerical model with those obtained from the experimental setup and published results.

3.6 Numerical Simulation

Figure 3.38 Contours of static temperature on the numerical model

Figure 3.38 shows the contours of static temperature change from the heat source to the heat sink. The efficient heat dissipation is simulated with the

CNF buckypaper TIM. The heating plate was set to 60 W/m2 and the cooling plate at 10 ºC. The thickness of the TIM was 0.0254 millimeter. Two temperature

77 readings were obtained at contacting block to measure the thermal performance of the TIM.

Figure 3.39 Static temperature distributions on the numerical model

Figure 3.39 shows the overall temperature distribution in terms of the position (geometry) of the test model. The close transition between the upper and lower block (green and red lines respectively), shows the reduction of air gaps resulting in the efficient heat transport path.

3.6.1 Parametric Study

For the designed model, the parametric study is carried out in terms of two parameters of the prepared thermal interface material.

78 a. Thickness of TIM

In an experimental setup, the thickness of the TIM decreases as the contact pressure increases. In order to duplicate this behavior in the numerical study, thickness of the TIM is taken into account as one of the influencing parameter in optimizing the total thermal performance. Five different thicknesses are simulated for each TIM. It is observed that as the thickness of the TIM decreases, more air gaps are filled and more heat is dissipated through the interface minimizing the total interfacial resistance. Figure 3.40 shows thermal impedance results obtained from the simulation on TIMs of different thicknesses.

CNF Buckypaper TIM 0.45 0.431

0.4 0.383

0.35 0.335

0.3 0.285

0.25

0.237 Thermal Impedance (K/W) ImpedanceThermal 0.2 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14

Thickness (mm)

Figure 3.40 Thermal impedance of various TIMs as a function of thickness

79 b. Thermal Performance of TIM

2 BP (CNF-HHT) 1.8 BP (CNF-LHT) 1.6 BP (CNF-PS) 1.4

1.2

1

0.8

0.6

0.4

0.2 ThermalImpedance (K/W) 0 0 50 100 150 200

Load (kg)

Figure 3.41 Thermal impedance of various TIMs as a function of load

The simulations were performed to carry out a parametric study of CNF-

TIMs in terms of thickness and thermal conductivity of the samples. Unlike experimental analysis, the applied load was duplicated as the thickness of the

TIM. Hence, as the load increases, the thickness decreases and hence thermal impedance decreases (Figure 3.41). Overall results from the numerical analysis show a strong agreement with those obtained from the experimental analysis.

80

Chapter IV

Power Spike Analysis

Faster speeds for space vehicles, high power requirements for automotive engines, and brighter beams for optical devices, has resulted in excessive heating of the respective systems. In case of power electronics this quest for miniaturization has lead to the generation of very high heat fluxes (around

107 W/ m2 ). The modern electronic and electrical systems are characterized by increase in the density of circuits within the chips. The denser the circuitry, the higher the power used to empower the chip and, consequently, the higher the heat generated. Intel-R microprocessors advanced from having 1.5 Am lines on

INTEL386k [135] to 0.13 Am in the current PentiumR-4 [136] within a space of 7 years. This is more than ten-fold reduction in the line width. The increase in heat flow from IC chips is apparent in the trend of power dissipation for high- performance chips reported by the International Technology Roadmap for the

Semiconductor (ITRS) [137–139]. Increasing amount of power dissipated has to be offset by efficient cooling because the maximum junction temperature stays at

85–90 C . A major challenge is the ability to manage the heat in the IC chips without compromising on the performance of the device. Thus, there is a growing

81 need to develop more effective TM techniques to control this high order of heat fluxes in different miniaturized devices. Conversely some of the electronic systems are subjected to transient heat loads. These pulsed heat inputs or sudden surge in heat flux called as ‗power spikes‘, can cause significant temperature fluctuations. Thus the management of this thermal energy is very crucial because heat has many detrimental effects on the device. It has been reported [140] that failure rates have near exponential dependence on device temperature. Also the dominant failure in most power electronics systems is thermal related [141]. When a certain upper critical temperature is reached, important parts of the device may cease to function. Furthermore, temperature cycles that result from switching on and off the device can cause problems, even if the operating temperature does not hit the upper critical temperature. This leads to another problem, the reliability of the device, because temperature cycling of more than 20 C can increase the failure rate of a device by eight-fold

[140]. There is therefore an additional challenge to the thermal manager: a need not only to dissipate heat from the device, but also to do it gradually. Thus, control of these power spikes or pulsed heat loads is one of the critical issues in the design of TM systems.

To meet the requirement on the safe operation area of the device, the thermal resistance of the device should be reduced by some kind of cooling. A number of researchers have used conventional interface materials to efficiently dissipate heat upon a sudden surge of power [142-147]. Baker et al. [142] performed a conduction-only analysis of the thermal performance of a

82 heterogeneous package. Two nonmetallic PCMs (octadecane and pentacosane) were used for heat storage with metal fins to enhance the effective thermal conductivity. Vesligaj and Amon [144] investigated passive thermal control using phase change materials during time-varying workloads on portable electronics.

An epoxy polymer was used as a phase change material. They indicated that the operational performance of the portable electronics improved when such a passive thermal storage device was used. Binet and Lacroix [145] conducted a numerical study of natural convection-dominated melting inside a rectangular enclosure from three discrete heat sources. Evans et al. [146] analyzed power electronic packages and provided design guidelines relating the materials, geometry, power input and junction temperature for steady-state conditions and transient pulses. Krishnan and Garimella [147] carried out a fully transient analysis of the phase change process inside a rectangular enclosure for electronics cooling applications with pulsed power dissipation. The experimental investigation included the melting of a pure PCM, n-eicosane, inside a rectangular aluminum container with multiple discrete heat sources mounted on one side. A copper heat spreader was used in conjunction with the heat sources.

The influence of changes in the frequency of pulses, heat source location and aspect ratio of the containment volume on the thermal performance of the PCM unit was studied. It was reported by the authors that the heat source location and container aspect ratio played a very significant role in the performance of the

PCM system with pulsed heat input.

83

Carbon nanotubes (CNTs), because of their size and excellent thermal properties, have proven to be one of the most efficient materials for heat transport [148-152]. Xu and Fisher [148] tested thermal resistance of one sided

CNT array interface (Si-CNT-Cu) under a pressure of 0.414 MPa with a pressure bar method and reported a thermal resistance of 20 mm2 K/W. Thermal resistances between the free CNT array tips and an experimental contact were

17 and 15 mm2 K/W at pressures of 0.040 and 0.100 MPa, respectively. Xu and

Fisher [149] also combined thin layers of phase change material (PCM) with CNT arrays, and the composite produced a resistance of 5 mm2 K/W under moderate pressures for Si–Cu interfaces. Wang et al. [150] used a photothermal technique to measure the thermal resistance between a CNT array and its growth substrate

(Si-CNT). The resistance was relatively large, 16 mm2 K/W, as the CNT array was of poor structural quality and no pressure was applied to the interface. Using a transient, Tong et al. [151] used a transient thermoreflectance technique to measure thermal resistance of one-sided CNT interface (Si-CNT-glass) and measured 18 mm2 K/W. The authors concluded that the interface between the free CNT array tips and their opposing glass substrate (CNT-glass) dominated the total thermal interface resistance and suggested that this resistance could be further decreased by the application of pressure to the interface. Xu and Fisher

[152] fabricated and experimentally studied two-sided CNT interfaces with CNT arrays directly synthesized on Si wafers and Cu blocks. With vertically aligned

CNT arrays, an interface resistance less than 5 mm2 K/W for two-sided CNT interfaces was measured. However, due to CNT growth and fabrication

84 constraints e.g., high-temperature substrate requirement, difficulties in fabricating samples with identical arrays on sides of a test chip [152], and high costs, there is an increased fabrication and measurement uncertainty, which sometimes can be larger than the magnitude of thermal performance.

Carbon nanofibers (CNF), though not that perfect in structure and relatively less conductive in nature, are fairly economical and do not require high temperature substrate for growth. These materials have been increasingly used for aircraft structures and need to be engineered with efficient lightning strike protection to achieve tolerance comparable to metallic components. Gao et al.

[153] developed a paper made of carbon nanofibers and nickel nanostrands as a surface layer on the composite panels to explore their efficiency in transient power spikes. The porous, flexible, non-woven papers of nanofibers and nanostrands were prepared and incorporated onto the surface of carbon fiber reinforced polymer composites through resin transfer molding. Authors reported that the power strike tolerance correlated to the surface conductivities of composite materials and the surfaces of the materials were damaged due to the existence of organic binder. Ngo et al. [154] used electrodeposited Cu as gap filler to enhance the stability and thermal conductance of CNF arrays with a one- dimensional steady state reference bar method. The authors reported a thermal resistance of 25 mm2 K/W under a pressure of 0.414 MPa for Si-Cu interfaces.

These results show that CNFs have competitive thermal properties to those of

CNTs with superiority in terms of cost and fabrication constraints. This section

85 explores the potential of self standing CNF buckypapers to sustain transient power spikes by efficiently dissipating maximum heat to the sink.

4.1 Experimental Setup

In order to test the thermal properties of various carbon-based TIMs, it was necessary to design a test setup. Many previous experimental set ups have been created to measure the thermal properties of TIMs [106-108]. However many of these set ups are expensive and time consuming. A test set up that was accurate, cost effective, and capable of rapidly generating usable data was created for this experiment. A similar setup was reconstructed and reprogrammed to meet the specifications required for the power spike analysis.

Figure 4.1 shows the visual representation of the original test rig. The original design consisted of two main subassemblies; top plate with actuator and heater subassembly, coolant plate with sample holder and guide subassembly. The two subassemblies were connected to each other through two threaded aluminum rods. The upper heating plate consisted primarily of a strip heater with dimensions of 75 mm wide x 100 mm length attached to a Plexiglas plate and an aluminum plate with dimensions of 100 mm x 40 mm x 25 mm. This was connected to a 12 kg force actuator and a 120 watt heater. The lower assembly consisted of a cooling plate designed to maintain a constant ambient temperature. A ―U‖ shaped copper pipe was inserted under the cooling plate which was attached to plastic tubing at each end for cold water supply.

86

Modifications were made in place of the sample holder

Figure 4.1 Original design of the experimental setup [75]

The sample holder (figure 4.2) was designed with two aluminum blocks

(heating and cooling plates) as contacting surfaces. Three holes were drilled down the center of each block. These thermocouple holes were drilled 10 mm deep. It was necessary to mount the blocks to the test rig. Two holes were made in each block, on the faces directly to the left and right of the thermocouples.

These holes were drilled 8 mm deep in order to insure no interference with the thermal couples. The upper block was attached directly to the corresponding hot plate. For the bottom block, two more brackets were created and screwed into the aluminum block. These brackets were then attached to a support block located at the back of the test rig. The two blocks were aligned to insure an accurate contact area possible. To further insure optimal thermal conduction through the blocks, a high thermal conductivity paste (Omegatherm 201) was

87 spread at the interface of the surfaces of the blocks and their corresponding plates.

Figure 4.2 Schematic of the sample holder

88

The test setup was designed for thermal analysis of carbon nanofiber based thermal interface materials. The compliant buckypaper was placed on the lower block in the sample holder. The electromechanical actuator was used to apply and remove the strip heater on the top of the upper block as required.

Thus, a thermal analysis of different samples could be performed by applying heat load at the top of block and dissipating the heat absorbed within the sample by a coolant flow at its bottom through the coolant plate. The above procedure was used for different techniques analyses such as measuring the temperature distribution within samples to analyze their heat transfer ability, TIM‘s ability to handle maximum temperature, and the global thermal resistance. Figure 4.3 shows the experimental setup with the sample holder.

89

Figure 4.3 a) Experimental setup b) Enlarged view of the sample holder

Before stepping into the actual analysis, the temperatures sensors were tested under dry test at ambient temperature. Six thermocouples were connected to the DAQ system and then inserted in two aluminum blocks (three in each) to verify the temperature. The thermocouple readings were noted for 10 minutes exposed to ambient conditions. Figure 4.4 shows the calibration analysis of six thermocouples as a function of time.

90

Thermocouples Calibration 40 38 36 34 32 30 28 T1 T2 26 T3 24 T4 Temperature (ºC) Temperature T5 22 T6 20 1 51 101 151 201 251 301 351 401 451 Time (sec)

Figure 4.4 Calibration of the six thermocouples

The following equations are used to calculate thermal conductivity of a TIM assuming a one-dimensional heat flow through the thickness of the sample. The governing equation is obtained by using Fourier‘s Law:

(Eq)

Where:

k = thermal conductivity

q = heat flow (or power)

t = interface thickness

91

A = area

= temperature difference across the interface

However, it is important to know that the thermal conductivity of a TIM is solely a material property. In other words, the thermal conductivity reflects the type of material used; it does not account for any thermal effects caused by the interfacial resistance or the geometry of the sample. Hence, it is important to take into account the thermal effects on the interface of TIM and contacting surfaces. The global resistance can be determined by the below equation

Where:

R= thermal resistance

= component temperature

= heat sink temperature

Q= heat flow (or power)

By using the equation of the thermal resistance, we can calculate each TIM‘s ability to resist sudden change in temperature. Therefore the significance of heat treatment on carbon structures can be highlighted by analyzing the change in temperature from heat source to the heat sink.

92

4.2 Results

Two tests were set up to examine two different types of heat transfer through each buckypaper TIM; a) the change in temperature between upper and lower surface which determines the ability of a TIM to effectively transport heat and b) maximum temperature reached at the upper surface which determines the ability of a TIM to handle high temperatures. To maximize the accuracy of the results, two thermocouple (closest to the interface) readings were examined.

These tests were designed to imitate typical power demands put upon modern electronics.

4.2.1 Uniform Power Pulse - 24 watts of power was applied at the upper block

for 60 minutes period. To obtain this, maximum power (120 watts) of

power was applied for 5 seconds followed by a 20 second period with zero

watts being applied. This gives 20% of the total power that is 24 watts.

The heater was then turned off, applying zero power for an additional 25

seconds. This test simulates the power usage of a device that uses a

constant power for a short period of time and then switch to stand by

condition (does not utilize any power). This is common for many devices

that are frequently turned on and off.

4.2.2 High Power Spikes - 12 watts of power was applied at the upper block for

60 minutes period. To obtain this, 120 watts of power was applied for 5

seconds followed by a period of 45 seconds with zero watts being applied.

This gives 10% of the 120 watt heater that is 12 watts. This procedure was

followed by a high power heat spike of 96 watts. To obtain that, 120 watts

93

of power was applied for 40 seconds followed by a period of 10 seconds

with zero watts being applied. This gives 80% of the 120 watt heater that

is 96 watts. These high power heat spikes simulate various demands put

upon electronics. With most devices today, power fluctuates constantly

and can experience instances of sudden surge of power.

In order to understand the effects of carbon nanofiber TIM compared the dry contact, and heat treatment effects of these TIM, both tests were run with;

A) Direct contact (no TIM)

Due to the rough nature of contacting surfaces (aluminum blocks in this

case), the dry contact has an immediate disadvantage in terms of interfacial

thermal resistance. Thermal resistance in dry contact therefore varies with

respect to the evenness of contact surfaces which can never be perfect.

Hence, in any case it can be very misleading and inaccurate to assume that

the dry contact performs better than a TIM in real world applications.

Figures 4.5 and 4.6 show the direct contact plots for uniform pulse and

power spikes respectively.

94

Direct Contact - 24 Watt Uniform Pulse 58 56 Dry Contact - TC 1 54 Dry Contact - TC 2 52

50 C) 48 Tmax = 48.85 ºC 46 44 42 40 ∆T = 8.60 ºC 38 36

Temperature ( Temperature 34 32 30 28 26 24 1 201 401 601 801 1001 1201 1401 1601 1801 2001 2201 2401

Time (sec)

Figure 4.5 Direct contact analyses for uniform heat pulse

Direct Contact - 96 Watt Power Spike

155 Dry Contact - TC 1 145 Dry Contact - TC 2 135 º

125 Tmax = 132.82 C C)

° 115 105 95 ∆T = 43.44 ºC 85 75

Temperature ( Temperature 65 55 45 35 25 1 201 401 601 801 1001 1201 1401 1601 1801 2001 2201 2401 2601 2801

Time (sec)

Figure 4.6 Direct contact analyses for transient power spikes

95

In order to compare the results of the direct contact analysis, compliant

TIMs of various heat treatment temperatures are used. Although the insertion of a compliant TIM adds a resistance to the system, but high conformability of the

TIM leads to the enhancement of the global resistance of the system.

B) Buckypaper TIM (CNF-HTT)

Buckypapers (CNF-HHT) - 24 Watt Uniform Pulse

58 BP2 - TC1 56 BP2 - TC2 ∆T = 12.77 ºC 54 BP3 - TC1 Tmax = 53.62 ºC 52 50 BP3 - TC2 48 BP1 - TC1 º ∆T = 9.71 ºC C) Tmax = 50.61 C ° 46 BP1 - TC2 ∆T = 4.45 ºC 44 Tmax = 46.35 ºC 42 40 38 36 Temperature ( Temperature 34 32 30 28 26 24 1 201 401 601 801 1001 1201 1401 1601 1801 2001 2201 2401 2601

Time (sec)

Figure 4.7 CNF-HHT buckypaper analyses for uniform heat pulse

96

Buckypapers (CNF-HHT) - 96 Watt Power Spike

155 BP2 - TC1 BP2 - TC2 145 BP3 - TC1 Tmax = 141.29 ºC 135 BP3 - TC2 BP1 - TC1 125 ∆T = 57.93 ºC

C) BP2 - TC2 ° 115 Tmax = 128.64 ºC 105 º ∆T = 42.60 ºC 95 Tmax = 109.12 C ∆T = 19.95 ºC 85

Temperature( 75 65 55 45 35 25 1 201 401 601 801 1001 1201 1401 1601 1801 2001 2201 2401 2601 2801

Time (sec)

Figure 4.8 CNF-HHT buckypaper analyses for transient power spikes

Three buckypaper TIMs composed of high heat treated (HHT) carbon nanofibers were tested at uniform heat pulse (figure 4.7) and transient power spikes (figure 4.8). The heat dissipation rate was observed in terms of temperature difference from the source to the sink. The maximum temperature values obtained for each TIM show the sustainability of a TIM at uniform and transient temperatures.

97

C) Buckypaper (CNF-LHT)

Buckypapers (CNF-LHT) - 24 Watt Uniform Pulse

58 BP1 - TC1 Tmax = 57.64 ºC 56 BP1 - TC2 54 BP2 - TC1 52 Tmax = 56.73 ºC 50 BP2 - TC2

48 C)

° 46 44 ∆T = 19.27 ºC 42 40 ∆T = 20.38 ºC 38 36 Temperature( 34 32 30 28 26 24 1 201 401 601 801 1001 1201 1401 1601 1801 2001 2201 2401 2601 Time (sec)

Figure 4.9 CNF-LHT buckypaper analyses for uniform heat pulse

Buckypapers (CNF-LHT) - 96 Watt Power Spike

155 BP1 - TC1 145 BP1 - TC2 135 BP2 - TC1

125 BP2 - TC2 C)

° 115 T = 155.62 ºC Tmax = 151.21 ºC 105 max 95 85 ∆T = 86.11 ºC ∆T = 80.0 8ºC 75 Temperature ( Temperature 65 55 45 35 25 1 201 401 601 801 1001 1201 1401 1601 1801 2001 2201 2401 2601 2801

Time (sec)

Figure 4.10 CNF-LHT buckypaper analyses for transient power spikes

98

Two buckypaper TIMs composed of low heat treated (LHT) carbon nanofibers were tested at uniform heat pulse (figure 4.9) and transient power spikes (figure 4.10). The heat dissipation rate was observed in terms of temperature difference from the source to the sink. The maximum temperature values obtained for each TIM show the sustainability of a TIM at uniform and transient temperatures.

D) Buckypaper (CNF-PS)

Buckypapers (CNF-PS) - 24 Watt Uniform Pulse

T = 63.64 ºC 62 BP1 - TC1 max 60 58 BP1 - TC2 T = 59.37 ºC 56 max BP2 - TC1 54

52 BP2 - TC2 ∆T = 20.01 ºC C) ° 50 48 46 44 42 ∆T = 21.72 ºC 40

Temperature ( Temperature 38 36 34 32 30 28 26 24 1 201 401 601 801 1001 1201 1401 1601 1801 2001 2201 2401 2601

Time (sec)

Figure 4.11 CNF-PS buckypaper analyses for uniform heat pulse

99

Buckypapers (CNF-PS) - 96 Watt Power Spike

155 BP1 - TC1 145 BP1 - TC2 135 BP2 - TC1 125 BP2 - TC2

C) 115 Tmax = 162.11 ºC Tmax = 153.47 ºC ° 105 95 85 ∆T = 91.82 ºC ∆T = 86.11 ºC 75 Temperature ( Temperature 65 55 45 35 25 1 201 401 601 801 1001 1201 1401 1601 1801 2001 2201 2401 2601 2801

Time (sec)

Figure 4.12 CNF-PS buckypaper analyses for transient power spikes

Two buckypaper TIMs composed of pyrolitically stripped (PS) carbon nanofibers were tested at uniform heat pulse (figure 4.11) and transient power spikes (figure 4.12). The heat dissipation rate was observed in terms of temperature difference from the source to the sink. The maximum temperature values obtained for each TIM show the sustainability of a TIM at uniform and transient temperatures.

100

Comparison - 24 Watts Uniform Pulse 70 CNF-HTT (T1) CNF-HTT (T2) 65 CNF-LHT (T1) CNF-LHT (T2) 60 CNF-PS (T1) CNF-PS (T2) 55 Dry Contact (T1) Dry Contact (T2) ºC) 50 45 40

Temperature ( Temperature 35 30 25 20 1 6 11 16 21 26 31 36 41

Time (min)

Figure 4.13 Temperature comparison of various heat treated CNF buckypapers

and direct contact subjected to uniform heat pulse

Comparison - 96 Watts Power Spike

200 CNF-HHT (T1) CNF-HHT (T2) 180 CNF-LHT (T1) CNF-LHT (T2) CNF-PS (T1) 160 CNF-PS (T2) Dry Contact (T1) ºC) 140 Dry Contact (T2) 120 100

80 Temperature ( Temperature 60 40 20 1 6 11 16 21 26 31 36 41 46 51 Time (min)

Figure 4.14 Temperature comparison of various heat treated CNF buckypapers

and direct contact subjected to transient power spikes 101

Figure 4.13 and 4.14 show the overall temperature comparison of direct contact and CNF based TIMs for uniform heat and power spikes respectively.

Two thermocouple readings were taken for each test to calculate the temperature difference.

Table 4.1 Maximum temperature results for uniform pulse analyses

24 Watt Uniform Pulse over 60 min. period

TIM Maximum Temperature (°C) ∆T (°C)

Dry Contact 48.85 Baseline

CNF-HHT 46.35 2.50

CNF-LHT 56.73 -7.88

CNF-PS 59.37 -10.52

Table 4.2 Temperature difference results for uniform pulse analyses

24 Watt Uniform Pulse over 60 min. period

TIM Temperature Difference (°C) ∆T (°C)

Dry Contact 8.60 Baseline

CNF-HHT 4.45 4.15

CNF-LHT 19.27 -10.67

CNF-PS 21.72 -13.12

102

Table 4.3 Maximum temperature results for power spikes analyses

96 Watt Power Spike over 60 min. period

TIM Maximum Temperature (°C) ∆T (°C)

Dry Contact 132.82 Baseline

CNF-HHT 109.12 23.70

CNF-LHT 151.21 -18.39

CNF-PS 153.47 -20.65

Table 4.4 Temperature difference results for power spikes analyses

96 Watt Power Spike over 60 min. period

TIM Temperature Difference(°C) ∆T (°C)

Dry Contact 43.44 Baseline

CNF-HHT 19.95 23.49

CNF-LHT 80.08 -36.64

CNF-PS 86.11 -42.67

The most striking and notable inference made from the tables 4.1-4.4 is that buckypaper with high heat treated (HHT) nanofibers reduced the thermal resistance and enhanced heat transport in comparison to the dry contact. For the uniform heat pulse test, the maximum temperature of CNF-HHT-TIM is 46

Celsius and temperature difference between two surfaces was seen to be less than 5 Celsius. This shows that high heat treatment of carbon nanofibers allows the fibers to align in a columnar way and minimizes spacing between fibers. This allows an efficient path for heat to transfer from the source to the sink. For the 103 power spike test, the maximum temperature of CNF-HHT TIM is far better than that of the dry contact. Similarly, significant reduction in temperature difference states the ability to HHT-CNFs to resist sudden surge of power and efficient dissipation to the heat sink.

Thermal interface materials with low heat treated (LHT) and pyrolitically stripped (PS) fibers, however, show increase in thermal resistance and hence do not provide efficient heat path. Significant increase in thermal resistance is seen in these samples due to the randomness of the grown carbon nanofibers.

4.3 Summary

An experimental setup was designed and constructed to perform thermal analysis of carbon nanofiber based thermal interface materials. The setup consisted of four basic components; Main assembly, sample holder, data acquisition device and a control unit. The main assembly includes a heater plate with actuator, a coolant plate with sample holder. The sample holder consisted of two aluminum blocks as contacting surfaces attached to each side of the main assembly. The sample (TIM) was sandwiched between the two contacting blocks under test. A high performance analog/digital data acquisition system, Keithley

2700, was used to measure the power output from strip heater along with the required output data. Efficient data-logging software, ExceLINX, was used as an interface program for transferring the measured data from the DAQ. The pre- installed software was used to directly stream data from the DAQ board into the

Excel. Software was programmed to analyze the data as it was received in Excel

104 with Excel's graphics, charting, and mathematical capabilities. An electronic control unit with the capability to control the speed of actuator along with its dwell time at the initial and retracted positions was used. The control circuit was programmed on MATLAB. Two detailed numerical codes were written for uniform heat pulse condition and spike power condition. The complete setup was calibrated and validated before actual tests.

Various heat treated carbon nanofiber based materials discussed in chapter 3 were used for detailed thermal analysis. All TIM samples were subjected to two types of heat loads; uniform pulse and transient power spikes.

Uniform heat load of 24 Watts was applied over a period of 60 minutes. Transient power spikes of 96 Watts were applied over a period of 60 minutes. Results from all samples were compared to those of direct contact (without a TIM on the interface).

The above analysis emphasized the importance of using high heat treated carbon nanofiber based thermal interface materials. CNF-HTT sample outperformed all others with major enhancement in thermal transport from heat source to heat sink. This was due to the graphitization of fibers and its aligned nature at extremely high temperatures.

105

Chapter V

Conclusions

The present research is focused on the thermal management issues in electronics. These issues are typically caused by thermal interface resistance and pulsed heat loads. A comprehensive literature review provided in chapter 2, highlights these issues with existing solutions. Those solutions however have number of shortcomings that limit the overall performance of a system. By keeping those in view, a step wise process of carbon nanofiber based compliant thermal interface material (TIM) fabrication, experimental and numerical characterization and, analysis of real world applications is carried out. These processed materials are expected to overcome existing limitations by escalating the efficiency of the heat transport.

A detailed analysis was performed in chapter 3 to find the best material that could conform to the mating surface and minimize the thermal interfacial resistance. Various kinds of heat treated Carbon nanofibers (CNFs) were used namely 1) Pyrolytically stripped (PS), Light heat treated (LHT) and High heat treated (HHT). These nanofibers were mixed with a water soluble binder

Polyvinyl Alcohol (PVA). The performance of CNFs was related to their degree of graphitizability. The higher the degree of graphitization and the better the heat transport across it. Various types of thermal interface materials were fabricated to

106 understand an efficient process to achieve pure carbon nanofiber materials with high thermal properties. The thermal properties of the final sample depend on the ratio of CNF to PVA, ultrasonic process time and number of agglomerates present in the solution.

The next step in the research was to design and construct an experimental setup that provides the characterization of the processed materials.

The complete setup was validated using commercially available TIM (Arctic

Silver 5) subjected to various thicknesses and loads. Actual experiments were performed without a TIM first (dry contact), and with fabricated TIM samples. The extension in the displacement cell measures the thickness of the TIM under test which is used to calculate the effective thermal conductivity of the TIM.

The main purpose of the test setup was to characterize the TIMs in terms of effects of the ratio of CNF and PVA on the global thermal resistance, and thermal impedance as a function of load and thickness. Four buckypaper samples with different thicknesses were tested in order to find the influence of those on global thermal resistance. It was found that as the load increases, TIM thermal conductivity increases with the load. This is due to significant decrease in

TIM density. The heat transfer coefficient increases due to better conformability of the TIM at higher loads. The initial porosity of the buckypaper was about 40% but at higher loads, because of compression, carbon fibers tend to get closer allowing better heat transport. It is well known that the thermal conductivity of air is lower than that of carbon nanofibers and contact points increase with the

107 pressure, the thermal conductance through the thermal interface material increases.

For the thermal impedance comparison, the TIMs of approximately 0.1 mm thickness were selected. It was seen that the thermal impedance values with high heat treated CNF (HHT-CNF) buckypaper were significantly better than low heat treated (LHT) and PS (pyrolytically stripped). For the buckypaper with HHT-

CNF as compared to the direct contact between two surfaces at maximum load, the thermal impedance was reduced by 54% from 0.151 K/W to 0.069 K/W.

In order to analyze the CNF to PVA weight ratio, four samples of CNF-

HTT buckypaper were prepared. By varying the PVA weight and keeping the

CNF weight constant (0.5g), different thermal impedance behaviors were seen.

For each composition, two samples with different thicknesses were analyzed to find the thermal impedance as a function of thickness. It was noticed that two samples with same PVA to CNF ratio but different global weights show identical thermal behavior. It was hence important to define the samples by its PVA to

CNF ratio and not its global weight. Hence the thermal impedance graph in terms of PVA to CNF ratio was analyzed which showed that at low pressures, as the ratio between two materials increased, the thermal impedance increased.

However, the thermal impedance remains constant at high temperatures, regardless of the ratio between two materials.

The test setup was modeled using Fluent software package. The simulations were performed to carry out a parametric study of CNF-TIMs in terms of thickness and thermal conductivity of the samples. Overall results from the

108 numerical analysis show a strong agreement with those obtained from the experimental analysis.

The fabricated samples were then tested for real world power application using designed experimental setup. From the analysis, it was found that buckypaper with high heat treated (HHT) nanofibers reduced the thermal resistance and significantly enhanced heat transport in comparison to the dry contact. For the uniform heat pulse test, the maximum temperature of CNF-HHT-

TIM was 47 Celsius and temperature difference between two surfaces was seen to be less than 5 Celsius. This shows that high heat treatment of carbon nanofibers allows the fibers to align in a columnar way and minimizes spacing between fibers. This allows an efficient path for heat to transfer from the source to the sink. For the power spike test, the maximum temperature of CNF-HHT TIM was far better than that of dry contact, taking high contact points into account.

Similarly, major reduction in temperature difference states the ability to HHT-

CNFs to resist sudden surge of power and efficient dissipation to the heat sink.

Thermal interface materials with low heat treated (LHT) and pyrolytically stripped

(PS) fibers, however, show increase in thermal resistance and hence do not provide efficient heat path. Significant increase in thermal resistance is seen in these samples due to the randomness of the grown carbon nanofibers.

The above analysis emphasized the importance of using high heat treated carbon nanofiber based thermal interface materials. CNF-HTT sample outperformed all others with more than 50% enhancement in thermal transport

109 from heat source to heat sink. Because of their graphitic nature, these are highly conductive and hence can be used as a potential TIM for aerospace applications.

The overall research showed that in comparison to existing thermal interface materials like thermal pastes and greases, carbon nanofiber based materials are dry in nature and compliant enough to conform to the mating surfaces. Because of high thermal properties of carbon nanofibers, these materials provide better thermal transport from heat source to heat sink. The added advantages of these materials are these are installation friendly and can perform at extremely high temperatures without breaking or leaking out of the interface.

110

Chapter VI

Recommendations

Carbon Nanotubes (CNTs) have received significant attention among the researchers [87-98] because of its small diameter and ultra high thermal properties. The properties of these nanostructures can further be enhanced by altering its size and orientation of the tubes. As a next step to this research,

CNFs can be replaced by aligned CNTs to maximize the heat dissipation across the interface. CNTs can be grown on a substrate than can sustain temperatures of 700 ºC or more, and peeled off to be used as a self standing CNT buckypaper.

The aligned CNTs generate a fairly good adhesive force which can act as an efficient heat transport from source to sink. Figure 6.1 shows the forces required for a gecko to climb up the vertical wall. The toes of the gecko have a special adaptation that allows them to adhere to most surfaces without the use of liquids or surface tension. The spatula tipped setae on gecko footpads demonstrate that the attractive forces that hold geckos to surfaces are Van der Waals interactions between the finely divided setae and the surfaces themselves [155]. These Van der Waals interactions are strong enough to support up to 100 times the animals‘ bodyweight. Liehui [156] fabricated macroscopic patches of compliant vertically aligned CNT array. The author explained that the fabricated patches of uniform

111 array had adhesive strength similar to that of geckos (10 N/cm2) on a variety of substrates and could be easily removed by peeling. When the array was patterned to mimic the hierarchical structures of gecko foot-hairs, strength was increased up to four times. The author reported that those patches were self- cleaning, non-visco-elastic and gave good strength in vacuum. Since these properties are highly desired in thermal management, microelectronics and space operations, CNT array can be a good candidate for efficient heat dissipation.

Figure 6.1 Adhesive forces of a single gecko foot [157]

112

Bibliography

[1] Shaikh S, Lafdi K, Silverman E. The effect of CNT interface on thermal resistance of contacting surfaces. Carbon 2007; (45):695–703.

[2] Li Y, Kumaran AR. Effect of carbon overcoat on the temperature rise of head/disc interface. Adv. Info. Storage 1991; (3):65–82.

[3] Kristiansen H. Thermal Management in Electronics (Goteborg, Sweden:

Calmers University of Technology, 2001).

[4] Bar-Cohen A, Kraus AD, Thermal Analysis and Control of Electronic

Equipment, Hemisphere Publishing, USA, 1983.

[5] Xu Y, Luo X, Chung DDL. Sodium Silicate Based Thermal Interface Material for High Thermal Contact Conductance. J. Electron. Packaging 2000;

122(2):128–31

[6] Leong CK, Chung DDL. Carbon black dispersions and carbon–silver combinations as thermal pastes that surpass commercial silver and ceramic pastes in providing high thermal contact conductance. Carbon 2004; (42):2323-7.

[7] Leong CK, Chung DDL. Carbon black dispersions as thermal pastes that surpass solder in providing high thermal contact conductance, Carbon 2003;

41(13):2459-69

[8] Leong CK, Aoyagi Y, Chung DDL. J. Electron. Mater. 2005; (34):1336.

113

[9] Howe TA, Leong CK, Chung DDL. Comparative Evaluation of Thermal

Interface Materials for Improving the Thermal Contact between an Operating

Computer Microprocessor and Its Heat Sink. Journal of Electronic Materials

2006; 35(8).

[10] Luo X, Chugh R, Biller BC, Hoi YM, Chung DDL. J. Electron. Mater. 2002;

31:535.

[11] Smalc M. Int. Electronic Packaging Technical Conference Exhibition New

York: ASME, 2003; (2):253.

[12] Savija I. Int. Electronic Packaging Technical Conference Exhibition New

York: ASME, 2003; (2):567.

[13] Marotta E, IEEE Trans. Compon. Packaging Tech. 2005; 28:102.

[14] Liu Z, Chung DDL. Calorimetric evaluation of phase change materials for use as thermal interface materials. Thermochimica Acta 2001; 366(2):135-47.

[15] Nurmawati MH. Int. J. Polymer Anal. 2004; 9:213.

[16] Webb R. Int. Electronic Packaging Technical Conference Exhibition New

York: ASME, 2003; (2):537.

[17] Luo X, Chung DDL. Int. J. Microcircuits Electron Packaging 2001; 24:141.

[18] Ngo Q. Surf. Eng. Mater. Sci. III 2005; 75–82.

[19] Zhang K. Electron. Compo. Trans. Conference 2005; (1):60.

[20] Lee T. Electron. Compo. Trans. Conference 2005; (1):55.

[21] Wu Y. Appl. Phys. Lett. 2005; 87:213108.

[22] Chung DDL. J. Mater. Sci. 1987; 22:4190.

114

[23] Walsh E, Grimes R. Low profile fan and heat sink thermal management solution for portable applications. International Journal of Thermal Sciences

2007; 46:1182–90.

[24] Kandasamy R, Wang X, Mujumdar AS. Application of phase change materials in thermal management of electronics, Applied Thermal Engineering

2007; 27:2822–32.

[25] Bardon JP. Introduction a‘ l'e‘ tude des resistances thermiques de contact.

Rev. Ge‘n. Therm. 1972; 125:429-46.

[26] Zhao JF, Wang AL, Yang CX. Prediction of thermal contact conductance based on the statistics of the roughness profile characteristics. International

Journal of Heat and Mass Transfer 2005; 48:974–85

[27] Defer D, Antezak E, Duthoit B. Measurement of low-thermal effusivity of building materials using thermal impedance method. Measurement Science and

Technology 2001; 12:549-56.

[28] Degiovanni A, Remy B, Andre S. Thermal resistance of a multi-constrictions contact: A simple model. International Journal of Heat and Mass Transfer 2003;

46:3727–35

[29] Tio KK, Toh KC. Thermal resistance of two solids in contact through a cylindrical joint. International Journal of Heat Mass transfer 1998; 41(13):2013-

24.

[30] Zhu D, Wang S. Experimental investigation of contact resistance in adsorber of solar adsorption refrigeration. Solar Energy 2002; 73(3):177–85.

115

[31] Parihar SK, Wright NT. Thermal contact resistance at elastomer to metal interfaces, Int. Comm. Heat Mass Transfer 1997; 24(8):1083-92.

[32] Wahid SMS, Madhusudana CV, Leonardi E. Solid spot conductance at low contact pressure. Experimental Thermal and Fluid Science 2004; 28:489–94.

[33] Bolger JC. Prediction and Measurement of Thermal Conductivity of Diamond

Filled Adhesives. Proceedings of the 42nd Electronic Components and

Technology Conference San Diego, CA 1992; 219–29.

[34] Mirmira SR, Marotta EE, Fletcher LS. Thermal Contact Conductance of

Adhesives for Microelectronic Systems. AIAA Journal of Thermophysics and

Heat Transfer 1997; 11(2):141-5.

[35] Kilik R, Davies R, Darwish SMH. Thermal Conductivity of Adhesive Filled with Metal Powders, International Journal of Adhesion and Adhesives 1989;

9(4):219–23.

[36] Lelea D. Effects of inlet geometry on heat transfer and fluid flow of tangential micro-heat sink. International Journal of Heat and Mass Transfer 2010; 53:3562–

9.

[37] Elshafei EAM. Natural convection heat transfer from a heat sink with hollow/perforated circular pin fins. Energy 2010; 35:2870-7

[38] Barba A, Musi B, Spiga M. Performance of a polymeric heat sink with circular microchannels, Applied Thermal Engineering 2006; 26:787–94.

[39] Chiang KT, Chou CC, Liu NM. Application of response surface methodology in describing the thermal performances of a pin-fin heat sink, International

Journal of Thermal Sciences 2009; 48:1196–205.

116

[40] Khor YK, Hung YM, Lim BK. On the role of radiation view factor in thermal performance of straight-fin heat sinks, International Communications in Heat and

Mass Transfer 2010; 37:1087–95

[41] Kim TH, Do KH, Kim DK. Closed form correlations for thermal optimization of plate-fin heat sinks under natural convection. International Journal of Heat and

Mass Transfer 2011; 54:1210–6.

[42] Huang RT, Sheu WJ, Wang CC. Orientation effect on natural convective performance of square pin fin heat sinks, International Journal of Heat and Mass

Transfer 2008 51:2368–76.

[43] Jouhara H, Axcell BP. Modelling and simulation techniques for forced convection heat transfer in heat sinks with rectangular fins, Simulation Modelling

Practice and Theory 2009; 17:871–82.

[44] Lorenzini M, Fabbri G, Salvigni S. Performance evaluation of a wavy-fin heat sink for power electronics. Applied Thermal Engineering 2007; 27:969–75.

[45] Fabbri G, Lorenzini M, Salvigni S. An innovative design of heat sinks for microelectronics. Proc. 4th European Thermal Sciences Conference,

Birmingham, 2004.

[46] Ismail MA, Abdullah MZ, Mujeebu MA. A CFD-based experimental analysis on the effect of free stream cooling on the performance of micro processor heat sinks, International Communications in Heat and Mass Transfer 2008; 35:771–8.

[47] Wang JC. Superposition method to investigate the thermal performance of heat sink with embedded heat pipe, International Communications in Heat and

Mass Transfer 2009; 36:686–92.

117

[48] Laraqi N, Bairi A. Theory of thermal resistance between solids with randomly sized and located contacts. International Journal of Heat and Mass Transfer

2002; 45:4175–80

[49] Wolfl EG, Schneider DA. Prediction of thermal contact resistance between polished surfaces, International Journal of Heat and Mass Transfer 1998;

41:3469-82.

[50] Zhang XZ, Zhang LW, Xing L, Study of thermal interfacial resistance between TC11/glass lubrication/K403 joint. Experimental Thermal and Fluid

Science 2010; 34:48–52.

[51] Incropera FP, DeWitt DP. Fundamentals of Heat and Mass Transfer, Fourth ed., Wiley, New York, 1996.

[52] Gwinn JP, Webb RL. Performance and testing of thermal interface materials,

Microelectronics Journal 2003; 34:215–22.

[53] Viswanath R, Wakharkar V, Watwe A, Lebonheur V. Thermal performance challenges from silicon to systems. Intel Technol. J. 2000; Q3.

[54] Tzeng JJW, Webber TW, Draqssowski DW. Technical review on thermal conductivity measurement techniques for thin thermal interfaces, Proc. 16th IEEE

Semi-Therm. 2000; 174–81.

[55] Salerno LJ, Kittel P, Spivak AL. Thermal conductance of pressed metallic contacts augmented with indium foil or Apiezon grease at liquid helium temperatures. Cryogenics 1994; 34(8):649–54

118

[56] He Y. Rapid thermal conductivity measurement with a hot disk sensor: Part

2. Characterization of thermal greases. Thermochimica Acta 2005; 436(1-2):130-

4

[57] Sartre V, Lallemand M. Enhancement of thermal contact conductance for electronic systems, Applied Thermal Engineering 2001; 21(2):221-35

[58] Jamison W, Sears G, Larsen G, Hunadi R. Thermally conductive, water cleanable greases In: 1991; 190–203

[59] Chung DDL. Materials for thermal conduction. Applied Thermal Engineering

2001; 21(16):1593-605

[60] Ouelette T, Sorgo M. Proc. Power Electron, Des. Conf. Power Sources

Users Conference, Cerritos, CA 1985; 134.

[61] Galli E. Carbon blacks. Plastics Compounding 1982; 5(2):22–32.

[62] Hu K, Chung DDL. Flexible graphite modified by carbon black paste for use as a thermal interface material. Carbon 2011; 49(4):1075-86.

[63] Stockwell Elastomerics and Inc, Thermal gap pads, thermally conductive pads 2006.

[64] Bergquist, Complete sil-pad selection guide 2008.

[65] Chomerics, thermal products selector guide 2008

[66] Lane GA. Solar Heat Storage: Latent Heat Material 1983; (I).

[67] Chen JCH, Eichelberger JL. Encapsulated Phase Change thermal Energy- storage Materials. US Patent 4,504,402 (1985).

[68] Lane GA, Rossow HE. Reversible Liquid/Solid Phase Change Composition for Storing Energy, US Patent 4,637,888 (1987).

119

[69] Kurklu A. Renewable Energy 1998; 13:89.

[70] Abhat A. Solar Energy 1983; 30:313.

[71] Himran S, Suwon A, Manosoori GA. Energy Sources 1994; 16:117.

[72] Krishnan S, Garimella SV. Thermal management of transient power spikes in electronics, phase change energy storage or copper heat sinks, ASME J Elect

Packag 2004; 126(3):308–16.

[73] Pal D, Joshi Y. Melting in a side heated tall enclosure by a uniformly dissipating heat source, Int J Heat Mass Transfer 2001; 44:375–87.

[74] Alawadhi EM, Amon CH. Performance analysis of an enhanced PCM thermal control unit. Thermal and thermomechanical phenomena in electronic systems. In: Proc ITherm 2000; 1:283–9.

[75] Shaikh S, Lafdi K. C/C composite, carbon nanotube and paraffin wax hybrid systems for the thermal control of pulsed power in electronics Carbon 2010;

48(3):813-24.

[76] Nurmawati MH, Siow KS, Rasiah IJ. Thermal performance of a phase change material on a nickel-plated surface. Thin Solid Films 2004; 462:481-6.

[77] Mey GD, Pilarski J, Wójcik M, Lasota M, Banaszczyk J, Vermeersch B,

Napieralski A et al. Influence of interface materials on the thermal impedance of electronic packages, International Communications in Heat and Mass Transfer

2009; 36:210–2.

[78] Jong KPD, Geus JW. Catalysis Reviews. Science and Engineering 2000;

(42):481–510.

120

[79] Mordkovich VZ. Carbon nanofibers: a new ultrahigh-strength material for chemical technology. Theoretical Foundations of Chemical Engineering 2003;

37(5):429–38.

[80] Tibbetts GG, Lake ML, Strong KL, Rice BP. A review of the fabrication and properties of vapor-grown carbon nanofiber/polymer composites. Composite

Science and Technology 2007; 67(7-8):1709–18.

[81] Choi YK, Sugimoto KI, Song SM, Endo M. Mechanical and thermal properties of vapor-grown carbon nanofiber and polycarbonate composite sheets, Material Letters 2005; 59(27):3514–20.

[82] Cortes P, Lozano K, Barrera EV, Bonilla-Rios J. Effects of nanofiber treatments on the properties of vapor-grown carbon fiber reinforced polymer composites. Journal of Applied Polymer Sciences 2003; 89(9):2527–34.

[83] Higgins BA, Brittain WJ. Polycarbonate carbon nanofiber composites.

European Polymer Journal 2005; 41(5):889–93

[84] Yang YL, Gupta, Dudley KL, Lawrence RW. A comparative study of EMI shielding properties of carbon nanofiber and multi-walled carbon nanotube filled polymer composites. Journal of Nanoscience and Nanotechnology 2005;

5(6):927–31.

[85] Rutkofsky M, Banash M, Rajagopal R, Chen J. Using a carbon nanotube additive to make electrically conductive commercial polymer composites. SAMPE

J 2005; 41(2):54–5.

[86] Trottier CM, Glatkowski P, Wallis P, Luo J. Properties and characterization of carbon-nanotube-based transparent conductive coating. J Soc Inform Display

121

2005; 13(9):759-63.

[87] Yu X, Rajamani R, Stelson KA, Cui T. Carbon nanotube based transparent conductive thin films. J Nanosci Nanotechnol 2006; 6(7):1939–44.

[88] Smith JG, Delozier DM, Connell JW, Watson KA.Carbon nanotube- conductive additive-space durable polymer nanocomposite films for electrostatic charge dissipation. Polymer 2004; 45(18):6133–42.

[89] Mu M, Walker AM, Torkelson JM, Winey KI. Cellular structures of carbon nanotubes in a polymer matrix improve properties relative to composites with dispersed nanotubes. Polymer 2008; 49(5):1332–7.

[90] Du F, Scogna RC, Zhou W, Brand S, Fischer JE, Winey KI. Nanotube networks in polymer nanocomposites: rheology and electrical conductivity.

Macromolecules 2004; 37(24):9048–55.

[91] Zeng JJ, Saltysiak B, Johnson WS, Schiraldi DA, Kumar S. Processing and properties of poly(methyl methacrylate)/carbon nano fiber composites. Compos

Part B-Eng 2004; 35(2):173–8.

[92] Wang Z, Liang Z, Wang B, Zhang C, Kramer L. Processing and property investigation of single-walled carbon nanotube (SWNT) buckypaper/epoxy resin matrix nanocomposites. Composites Part A: Applied Science and Manufacturing

2004 35(10):1225-32

[93] Zaeri MM, Ziaei-Rad S, Vahedi A, Karimzadeh F. Mechanical modelling of carbon nanomaterials from nanotubes to buckypaper. Carbon 2010;

48(13):3916-30

[94] Hong WT, Tai NH. Investigations on the thermal conductivity of composites

122 reinforced with carbon nanotubes. Diamond and Related Materials 2008; 17(7-

10):1577-81

[95] Xu J, Fisher TS. Enhancement of thermal interface materials with carbon nanotube arrays. International Journal of Heat and Mass Transfer 2006;

49:1658–66

[96] Cross R, Cola BA, Fisher T, Xu X, Gall K, Graham S, et al. A metallization and bonding approach for high performance carbon nanotube thermal interface materials. Nanotechnology 2010; 21:445705 (8pp)

[97] Cola BA, Xu J, Fisher TS. Contact mechanics and thermal conductance of carbon nanotube array interface. International Journal of Heat and Mass Transfer

2009; 52(15-16):3490-503

[98] Wang H, Feng JY, Hu XJ, Ng KM. Reducing thermal contact resistance using a bilayer aligned CNT thermal interface material. Chemical Engineering

Science 2010; 65(3):1101-8

[99] Singh IV, Tanaka M, Endo M. Effect of interface on the thermal conductivity of carbon nanotube composites, International Journal of Thermal Sciences 2007;

46:842-7.

[100] Schacht R, May D, Wunderle B, Wittler O, Gollhardt A, Michel B, Reichl H, et al. Characterization Of Thermal Interface Materials to support Thermal

Simulation, Nice, Côte d‘Azur, France, 27-29 September 2006.

[101] Early M, Lee S, Pellilo M. Thermal performance of interface material in microelectronics packaging applications.

123

[102] Konstantin I. Matveev, Gregory W. Swift, Scott Backhaus, Temperatures near the interface between an ideal heat exchanger and thermal buffer tube or pulse tube, Int. Journal of Heat and Mass Transfer 206; 49:868-78.

[103] Rodriguez NM. A review of catalytically grown carbon nanofibers. J Mater

Res 1993; 8(12):3233-50

[104] Dresselhaus MS, Dresselhaus G, Sugihara K, Spain IL, Goldberg HA.

Graphite fibers and filaments. New York: Springer-Verlag; 1998.

[105] Trimm DL. The formation and removal of coke from nickel catalyst. Catal

Rev Sci Eng 1977; 16:155-89

[106] Bartholowmew CH. Carbon deposition in steam reforming and methanation. Catal Rev Sci Eng 1982; 24(4):67-112

[107] Lake ML, Ting JM. In: T.D. Burchell, Editor, Carbon Materials for Advanced

Technologies, Pergamon Press, Oxford, UK (1999), 139–67.

[108] Tibbetts GG, Lake ML, Strong KL, Rice BP. A review of the fabrication and properties of vapor-grown carbon nanofiber/polymer composites, Compos Sci

Technol 2007; 67(7–8):1709–18.

[109] Lee BO, Woo WJ, Kim MS. EMI shielding effectiveness of carbon nanofiber filled poly(vinyl alcohol) coating materials, Macromol Mater Eng 2001;

286(2):114–8.

[110] Pervin F, Zhou Y, Rangari VK, Jeelani S. Testing and evaluation on the thermal and mechanical properties of carbon nano fiber reinforced SC-15 epoxy,

Materials Science and Engineering A 2005; 405:246–253.

124

[111] Mahanta NK, Abramson AR, Lake ML, Burton DJ, Chang JC, Mayer HK,

Ravine JL, et al. Thermal conductivity of carbon nanofiber mats. Carbon 2010;

(48):4457-65

[112] Patton RD, Pittman CU, Wang L, Hill JR. Vapor grown carbon fiber composites with epoxy and poly(phenylene sulfide) matrices. Composites Part A:

Applied Science and Manufacturing 1999; 30:1081–91

[113] Lafdi K, Matzek M. Carbon nanofibers as a nano reinforcement for polymeric nanocomposites. In: 48th International SAMPE symposium proceedings, Long Beach, USA; 2003.

[114] Zhou Y , Pervin F, Rangari VK, Jeelani S. Fabrication and evaluation of carbon nano fiber filled carbon/epoxy composite. Materials Science and

Engineering 2006; 426(1-2):221-8

[115] Chinthaginjala JK, Thakur DB, Seshan K, Lefferts L. How Carbon-Nano-

Fibers attach to Ni foam Carbon 2008; 46(13):1638-47

[116] Tibbetts GG, Kwag C, Howe JY, Lake ML. Heat treating carbon nanofibers for optimal composite performance. In: Proceedings of the forty-eighth

International SAMPE symposium, 2003, Long Beach.

[117] Howe JY, Tibbetts GG, Kwag C, Lake ML. Heat treating carbon nanofibers for optimal composite performance. Journal of Material Research 2006;

21(10):2646–2652.

[118] Zhu H, Li X, Ci L, Xu C, Wu D, Mao Z, et al. Hydrogen storage in heat- treated carbon nanofibers prepared by the vertical floating catalyst method.

Materials Chemistry and Physics 2003; 78(3):670-5

125

[119] Endo M, Kim YA, Hayashi T, Yanagisawa T, Muramatsu H, Ezaka M, et al.

Microstructural changes induced in ―stacked cup‖ carbon nanofibers by heat treatment, Carbon 2003; 41(10):1941–7.

[120] Endo M, Kim YA, Fukai Y, Hayashi T, Terrones M, Terrones H, et al.

Comparison study of semi-crystalline and highly crystalline multiwalled carbon nanotubes, Appl Phys Letts 2001; 79 (10):1531–3.

[121] Parrott EPJ, Zeitler JA, McGregor J, Oei SP, Unalan HE, Tan SC, Milne WI,

Tessonnier JP, Schlogl R, Gladden LF, et al. Understanding the Dielectric

Properties of Heat-Treated Carbon Nanofibers at Terahertz Frequencies: a New

Perspective on the Catalytic Activity of Structured Carbonaceous Materials, J.

Phys. Chem. C 2009; 113:10554–9

[122] Standard Test Method for Thermal Transmission Properties of Thermally

Conductive Electrical Insulation Materials, ASTM Standard D-5470-06, Copyright

ASTM International, Conshohocken, PA, 2007.

[123] Nishino K, Yamashita S, Torii K. Thermal Contact Conductance under Low

Applied Load in a Vacuum Environment. Experimental Thermal and Fluid

Science 1995; 10:258-71

[124] Chung DDL. Materials for thermal conduction, Applied Thermal Engineering

2001; 21:1593-605.

[125] Voller GP, Tirovic M. Conductive heat transfer across a bolted automotive jointand the influence of interface conditioning, International Journal of Heat and

Mass Transfer 2007; 50:4833–44.

126

[126] Cheng WW, Madhusudana C. Effect of electroplating on the thermal conductance of fin–tube interface, Applied Thermal Engineering 2006; 26:2119–

31.

[127] Day RJ, Lovell PA, Wazzan AA. Compos. Sci. Technol. 2001; 61:41–56.

[128] Bagheri R, Pearson RA. Polymer 2001; 41:269–76.

[129] Kawaguchi T, Pearson RA. Polymer 2003; 44:4239–47.

[130] Mahfuz H, Adnan A, Rangari VK, Jeelani S, Jang BZ. Compos. A: Appl.

Sci. Manufact. 2004; 35:519–27.

[131] Vigolo B, Penicaud A, Coulon C, Sauder C, Pailler R, Journet C et al.,

Macroscopic fibers and ribbons of oriented carbon nanotubes, Science 2000;

290(5495):1331–4.

[132] Dalton AB, Collins S, Munoz E, Razal JM, Ebron VH, Ferraris JP et al.,

Super-tough carbon-nanotube fibers, Nature 2003; 423(6941):703.

[133] Miaudet P, Badaire S, Maugey M, Derre A, Pichot V, Launois P et al., Hot- drawing of single and multiwall carbon nanotube fibers for high toughness and alignment, Nano Letters 5 2005; 11:2212–5.

[134] Oberlin A. Rouzaud JN. and Goma J. J. Chime physique 81pp.701-710

(1984).

[135] J. Ameen, R.M. Miller, P.V. Yokimens, 1998. Circuits and Devices, IEEE,

USA, 38– 40.

[136] Intel Inc.,

[137] The International Technology Roadmap for Semiconductor 2000.

[138] The International Technology Roadmap for Semiconductor 2001,

127

[139] The International Technology Roadmap for Semiconductor 2002,

[140] Bar-Cohen A, Kraus AD, Thermal Analysis and Control of Electronic

Equipment, Hemisphere Publishing, USA, Chap. 1, 4, 9, 1983.

[141] Kristiansen H. Thermal Management in Electronics, (Goteborg, Sweden:

Chalmers University of Technology, 2001).

[142] Baker KW, Jang JH, Yu JS. Thermal control of phase change package with periodic pulse heating—a case study,‖ in Proc. ASME/JSME Joint Thermal

Engineering Conference 1985; 4:463–9.

[143] Pal D, and Joshi Y. Application of phase change materials to thermal control of electronic modules: a computational study, J. Electron. Packag. 1997;

119:40–50.

[144] Vesligaj MJ, Amon CH. Transient thermal management of temperature fluctuations during time varying workloads on portable electronics, IEEE Trans.

Comp. Packag. Technol. 1999; 22:541–50.

[145] Binet B. and Lacroix M., 2000, Melting from heat sources flush mounted on a conducting vertical wall, Int. J. Numer. Methods Heat Fluid Flow, 10, 286–303.

[146] Evans AG, He MY, Hutchinson JW, Shaw M. Temperature distribution in advanced power electronics systems and the effect of phase change materials on temperature suppression during power pulses, J. Electron. Packag. 2001;

123:211–7.

[147] Krishnan S, Garimella SV. Analysis of a phase change energy storage system for pulsed power dissipation. IEEE Transactions on components and packaging technologies 2004; 27(1):191–9.

128

[148] Xu J, Fisher TS. IEEE Trans. Composite Packaging Technology 2006,

(29):261.

[149] Xu J, Fisher TS. Int. J. Heat Mass Transfer 2006; 49:1658.

[150] Wang X, Zhong Z, Xu J. J. Appl. Phys. 2005; 97:064302.

[151] Tong T, Zhao Y, Delzeit Y, Kashani A, Majumdar A, Meyyappan M,

Proceeding of the ASME International Mechanical Engineering Congress and

Exposition, Orlando, Florida, 2005

[152] Xu J, Fisher TS. Proceedings of the 18th National and 7th ISHMTASME

Heat and Mass Transfer Conference, IIT Guwahati, India, 2006.

[153] Gou J, Tang Y, Liang F, Zhao Z, Firsich D, Fielding J.

[154] Ngo Q, Gurden BA, Cassell AM, Walker MD, Ye Q, Koehne JE,

Meyyappan M, Li J, Yang CY et al. Nano Lett. 2003; 4:2403.

[155] Huber G. Evidence for capillarity contributions to gecko adhesion from single spatula nanomechanical measurements. Proceedings of the National

Academy of Sciences 2005; 102(45):16293

[156] Liehui G. Synthetic Gecko Adhesives and Adhesion in Geckos. PhD dissertation 2011; University of Akron.

[157] Autumn K, Liang YA, Hsieh ST, Zesch W, Chan WP, Kenny TW, Fearing R,

Full RJ et al. Nature 2000; 405:681-5.

129

APPENDIX

Figure A-1 Cooling Plate

Figure A-2 Aluminum Base Plate

130

Figure A-3 Heater Plate

Figure A-4 Plastic Insulator

131

Figure A-5 1x1 inch Aluminum Block

Figure A-6 2x2 inch Aluminum Block

132

Figure A-7 Test Setup - 1x1 inch blocks installed

Figure A-8 Test Setup - 2x2 inch blocks installed

133