Advanced thermal analysis of microelectronics using spreading resistance models

Thermal analysis of electronic devices is one of the most important steps for designing of modern devices. Precise thermal analysis is essential for designing an effective thermal management system of modern electronic devices such as batteries, LEDs, microelectronics, ICs, circuit boards, semiconductors and heat spreaders. For having a precise thermal analysis, the temperature profile and thermal spreading resistance of the device should be calculated by considering the geometry, property and boundary conditions. Thermal spreading resistance occurs when heat enters through a portion of a surface and flows by conduction. It is the primary source of thermal resistance when heat flows from a tiny heat source to a thin and wide heat spreader. In this thesis, analytical models for modeling the temperature behavior and thermal resistance in some common geometries of microelectronic devices such as heat channels and heat tubes are investigated. Different boundary conditions for the system are considered. Along the source plane, a combination of discretely specified heat flux, specified temperatures and adiabatic condition are studied. Along the walls of the system, adiabatic or convective cooling boundary conditions are assumed. Along the sink plane, convective cooling with constant or variable heat transfer coefficient are considered. Also, the effect of orthotropic properties is discussed. This thesis contains nine chapters. Chapter one is the introduction and shows the concepts of thermal spreading resistance besides the originality and importance of the work. Chapter two reviews the literatures on the thermal spreading resistance in the past fifty years with a focus on the recent advances. In chapters three and four, thermal resistance of a twodimensional flux channel with non-uniform convection coefficient in the heat sink plane is studied. The non-uniform convection is modeled by using two functions than can simulate a wide variety of different heat sink configurations. In chapter five, a non-symmetrical flux channel with different heat transfer coefficient along the right and left edges and sink plane is analytically modeled. Due to the edge cooling and non-symmetry, the eigenvalues of the system are defined using the heat transfer coefficient on both edges and for satisfying the orthogonality condition, a normalized function is calculated. In chapter six, thermal behavior of two-dimensional rectangular flux channel with arbitrary boundary conditions on the source plane is presented. The boundary condition along the source plane can be a combination of the first kind boundary condition (Dirichlet or prescribed temperature) and the second kind boundary condition (Neumann or prescribed heat flux). The proposed solution can be used for modeling the flux channels with numerous different source plane boundary conditions without any limitations in the number and position of heat sources. In chapter seven, temperature profile of a circular flux tube with discretely specified boundary conditions along the source plane is presented. Also, the effect of orthotropic properties are discussed. In chapter 8, a three-dimensional rectangular flux channel with a non-uniform heat convection along the heat sink plane is analytically modeled. In chapter nine, a summary of the achievements is presented and some systems are proposed for the future studies. It is worth mentioning that all the models and case studies in the thesis are compared with the Finite Element Method (FEM).

Experimental investigation of water, snow and granular ice effects on ice failure processes and impact loads

A large series of laboratory ice crushing experiments was performed to investigate the effects of external boundary condition and indenter contact geometry on ice load magnitude under crushing conditions. Four boundary conditions were considered: dry cases, submerged cases, and cases with the presence of snow and granular ice material on the indenter surface. Indenter geometries were a flat plate, wedge shaped indenter, (reverse) conical indenter, and spherical indenter. These were impacted with artificially produced ice specimens of conical shape with 20° and 30° cone angles. All indenter – ice combinations were tested in dry and submerged environments at 1 mm/s and 100 mm/s indentation rates. Additional tests with the flat indentation plate were conducted at 10 mm/s impact velocity and a subset of scenarios with snow and granular ice material was evaluated. The tests were performed using a material testing system (MTS) machine located inside a cold room at an ambient temperature of - 7°C. Data acquisition comprised time, vertical force, and displacement. In several tests with the flat plate and wedge shaped indenter, supplementary information on local pressure patterns and contact area were obtained using tactile pressure sensors. All tests were recorded with a high speed video camera and still photos were taken before and after each test. Thin sections were taken of some specimens as well. Ice loads were found to strongly depend on contact condition, interrelated with pre-existing confinement and indentation rate. Submergence yielded higher forces, especially at the high indentation rate. This was very evident for the flat indentation plate and spherical indenter, and with restrictions for the wedge shaped indenter. No indication was found for the conical indenter. For the conical indenter it was concluded that the structural restriction due to the indenter geometry was dominating. The working surface for the water to act was not sufficient to influence the failure processes and associated ice loads. The presence of snow and granular ice significantly increased the forces at the low indentation rate (with the flat indentation plate) that were higher compared to submerged cases and far above the dry contact condition. Contact area measurements revealed a correlation of higher forces with a concurrent increase in actual contact area that depended on the respective boundary condition. In submergence, ice debris constitution was changed; ice extrusion, as well as crack development and propagation were impeded. Snow and granular ice seemed to provide additional material sources for establishing larger contact areas. The dry contact condition generally had the smallest real contact area, as well as the lowest forces. The comparison of nominal and measured contact areas revealed distinct deviations. The incorporation of those differences in contact process pressures-area relationships indicated that the overall process pressure was not substantially affected by the increased loads.