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).

Paleoceanographic and paleoclimatic evolution of the Aegean Sea since the last interglacial (~130 ka BP) with special emphasis on sapropel formation

Five long piston cores collected from different subbasins of the Aegean Sea constitute the primary source of data for this PhD thesis. This study is the first to document a continuous paleoceanographic and paleoclimatic record of the Aegean Sea since the last interglacial. The chronostratigraphic reconstructions of the cored sediments based on organic carbon contents, stratigraphic position of known ash layers and oxygen isotopic curve matching collectively demonstrate the presence of sapropel S1 and MISS sapropels S3, S4 and S5 in the Aegean Sea subbasins. Generally, the organic carbon (TOC wt%) contents in sapropels range between 0.8% and 2% with highest concentrations of 9-13% in sapropels S4 and S5. Average sedimentation rates range between 4.7 and 11.8 cmlka with highest rates being observed in Euboea and North Ikaria basins (9.8 and 11.8 cm lka, respectively). The timing of the onset of sapropels S4 and S5 mostly predate those in the eastern Mediterranean with ages ranging from 106.4-105.6 and 128.6-128.4 ka BP, respectively. On the other hand, the initiation of the onset of sapropel S3 (i.e., 83.2-80.4 ka BP) seems to agree with its Mediterranean counterparts, which highlights the heterogeneity of the Aegean Sea subbasins in terms of rapid vs. lagged response to changing climatic conditions. The sapropel initiations appear to be synchronous across the Aegean Sea; whereas, the terminations display a wider temporal variability implying that the cessation of sapropels is controlled both by the amplitude of paleoclimatic changes and the physiography/location ofthe subbasins. Quantitative variations in the planktonic faunal assemblages exhibit a sequence of bioevents during the last -130,000 years which allow identification of four major biozones. The distributional patterns of the most significant taxa demonstrate similar trends among all core localities suggesting that the major changes in the planktonic foraminifera assemblages have taken place rather synchronously in the Aegean Sea. Sapropels S3, S4 and S5 were deposited under similar hydrographic conditions during which a distinct deep chlorophyll maximum (DCM) layer was established. This situation points to a stratified water column and increased export productivity during times of sapropel formation. On the other hand, the faunal contrast between Sl and older sapropels indicates that the former was developed in the absence of a DCM layer, lacking a deep phytoplankton assemblage. Under such conditions, oxygen advection via intermediate water flow must have been significantly reduced which implies significant stagnation. Sapropels are interpreted to have been deposited under normal marine conditions with temporary establishment of semi-euxinic bottom water conditions. Both marine and terrestrial organic matter contributed equally to MISS sapropels. In addition, organic carbon isotopic values across sapropels are more depleted than those in the eastern Mediterranean which, in tum, suggests enhanced riverine input during their deposition. Primary productivity calculations show that, particularly for sapropels with very high TOC values, both preservation and increased productivity are imperative in order to deposit sapropels with very high organic carbon contents (i.e., up to 13%).