Optimization of three-dimensional branching networks of microchannels for thermal management of microelectronics

The aim of this work is to present a methodology to develop cost-effective thermal management solutions for microelectronic devices, capable of removing maximum amount of heat and delivering maximally uniform temperature distributions. The topological and geometrical characteristics of multiple-story three-dimensional branching networks of microchannels were developed using multi-objective optimization. A conjugate heat transfer analysis software package and an automatic 3D microchannel network generator were developed and coupled with a modified version of a particle-swarm optimization algorithm with a goal of creating a design tool for 3D networks of optimized coolant flow passages. Numerical algorithms in the conjugate heat transfer solution package include a quasi-ID thermo-fluid solver and a steady heat diffusion solver, which were validated against results from high-fidelity Navier-Stokes equations solver and analytical solutions for basic fluid dynamics test cases. Pareto-optimal solutions demonstrate that thermal loads of up to 500 W/cm2 can be managed with 3D microchannel networks, with pumping power requirements up to 50% lower with respect to currently used high-performance cooling technologies.

Database of the GLAMAP project

A uniform chronology for foraminifera-based sea surface temperature records has been established in more than 120 sediment cores obtained from the equatorial and eastern Atlantic up to the Arctic Ocean. The chronostratigraphy of the last 30,000 years is mainly based on published d18O records and 14C ages from accelerator mass spectrometry, converted into calendar-year ages. The high-precision age control provides the database necessary for the uniform reconstruction of the climate interval of the Last Glacial Maximum within the GLAMAP-2000 project.

Volume averaging of multiphase flows with hydrate formation in subsea pipelines

In oil and gas pipeline operations, the gas, oil, and water phases simultaneously move through pipe systems. The mixture cools as it flows through subsea pipelines, and forms a hydrate formation region, where the hydrate crystals start to grow and may eventually block the pipeline. The potential of pipe blockage due to hydrate formation is one of the most significant flow-assurance problems in deep-water subsea operations. Due to the catastrophic safety and economic implications of hydrate blockage, it is important to accurately predict the simultaneous flow of gas, water, and hydrate particles in flowlines. Currently, there are few or no studies that account for the simultaneous effects of hydrate growth and heat transfer on flow characteristics within pipelines. This thesis presents new and more accurate predictive models of multiphase flows in undersea pipelines to describe the simultaneous flow of gas, water, and hydrate particles through a pipeline. A growth rate model for the hydrate phase is presented and then used in the development of a new three-phase model. The conservation equations of mass, momentum, and energy are formulated to describe the physical phenomena of momentum and heat transfer between the fluid and the wall. The governing equations are solved based on an analytical-numerical approach using a Newton-Raphson method for the nonlinear equations. An algorithm was developed in Matlab software to solve the equations from the inlet to the outlet of the pipeline. The developed models are validated against a single-phase model with mixture properties, and the results of comparative studies show close agreement. The new model predicts the volume fraction and velocity of each phase, as well as the mixture pressure and temperature profiles along the length of the pipeline. The results from the hydrate growth model reveal the growth rate and location where the initial hydrates start to form. Finally, to assess the impact of certain parameters on the flow characteristics, parametric studies have been conducted. The results show the effect of a variation in the pipe diameter, mass flow rate, inlet pressure, and inlet temperature on the flow characteristics and hydrate growth rates.

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

Coupled transport in rotor models

Acknowledgement One of us (AP) wishes to acknowledge S. Flach for enlightening discussions about the relationship between the DNLS equation and the rotor model.

Investigação experimental das propriedades termofísicas e da convecção forçada de nanofluido de grafeno

This work presents an experimental investigation of thermal hydraulic performance of the nanofluid composed by graphene nanoparticles dispersed in a mixture of water and ethylene glycol at a ratio of 70:30% by volume. The tests were carried out under forced convection inside a circular tube with uniform heat flux on the wall for the laminar-turbulent transition regime. The mass flow rate ranged from 40 to 70 g/s corresponding to Reynolds numbers between 3000 and 7500. The heat flux was maintained constant at values of 11, 16 and 21 kW/m², as well as the inlet temperature of 15, 20 and 25°C. Three samples were produced with the nanofluid volumetric concentration of 0.05%, 0.10% and 0.15%. Thermophysical properties were experimentaly measured for all samples that were critically compared and discussed with theoretical models most commonly used in the literature. Initially, experiments with distilled water confirmed the validity of the experimental equipment for the thermo-hydraulic tests. Therefore, nanofluid samples that showed the highest thermal conductivity, corresponding to the volumetric concentrations of 0.15% and 0.10%, were subjected to the tests. The thermal-hydraulic performance for both samples was unsatisfactory. The heat transfer coefficients for convection of nanofluids reduced 21% in average, for the sample with = 0.15% and 26% and for =0.10%. The pressure drop of the samples was higher than the base fluid. Finally, the pressure drop and heat transfer coefficient by convection of both samples were also compared to theoretical models. The models used for pressure drop showed an excellent agreement with experimental results, which is remarkable considering the transitional flow.

Constructal Design of Energy Systems

This dissertation shows the use of Constructal law to find the relation between the morphing of the system configuration and the improvements in the global performance of the complex flow system. It shows that the better features of both flow and heat transfer architecture can be found and predicted by using the constructal law in energy systems. Chapter 2 shows the effect of flow configuration on the heat transfer performance of a spiral shaped pipe embedded in a cylindrical conducting volume. Several configurations were considered. The optimal spacings between the spiral turns and spire planes exist, such that the volumetric heat transfer rate is maximal. The optimized features of the heat transfer architecture are robust. Chapter 3 shows the heat transfer performance of a helically shaped pipe embedded in a cylindrical conducting volume. It shows that the optimized features of the heat transfer architecture are robust with respect to changes in several physical parameters. Chapter 4 reports analytically the formulas for effective permeability in several configurations of fissured systems, using the closed-form description of tree networks designed to provide flow access. The permeability formulas do not vary much from one tree design to the next, suggesting that similar formulas may apply to naturally fissured porous media with unknown precise details, which occur in natural reservoirs. Chapter 5 illustrates a counterflow heat exchanger consists of two plenums with a core. The results show that the overall flow and thermal resistance are lowest when the core is absent. Overall, the constructal design governs the evolution of flow configuration in nature and energy systems.

Design of Latent Thermal Energy Storage Systems

This dissertation documents the results of a theoretical and numerical study of time dependent storage of energy by melting a phase change material. The heating is provided along invading lines, which change from single-line invasion to tree-shaped invasion. Chapter 2 identifies the special design feature of distributing energy storage in time-dependent fashion on a territory, when the energy flows by fluid flow from a concentrated source to points (users) distributed equidistantly on the area. The challenge in this chapter is to determine the architecture of distributed energy storage. The chief conclusion is that the finite amount of storage material should be distributed proportionally with the distribution of the flow rate of heating agent arriving on the area. The total time needed by the source stream to ‘invade’ the area is cumulative (the sum of the storage times required at each storage site), and depends on the energy distribution paths and the sequence in which the users are served by the source stream. Chapter 3 shows theoretically that the melting process consists of two phases: “invasion” thermal diffusion along the invading line, which is followed by “consolidation” as heat diffuses perpendicularly to the invading line. This chapter also reports the duration of both phases and the evolution of the melt layer around the invading line during the two-dimensional and three-dimensional invasion. It also shows that the amount of melted material increases in time according to a curve shaped as an S. These theoretical predictions are validated by means of numerical simulations in chapter 4. This chapter also shows that the heat transfer rate density increases (i.e., the S curve becomes steeper) as the complexity and number of degrees of freedom of the structure are increased, in accord with the constructal law. The optimal geometric features of the tree structure are detailed in this chapter. Chapter 5 documents a numerical study of time-dependent melting where the heat transfer is convection dominated, unlike in chapter 3 and 4 where the melting is ruled by pure conduction. In accord with constructal design, the search is for effective heat-flow architectures. The volume-constrained improvement of the designs for heat flow begins with assuming the simplest structure, where a single line serves as heat source. Next, the heat source is endowed with freedom to change its shape as it grows. The objective of the numerical simulations is to discover the geometric features that lead to the fastest melting process. The results show that the heat transfer rate density increases as the complexity and number of degrees of freedom of the structure are increased. Furthermore, the angles between heat invasion lines have a minor effect on the global performance compared to other degrees of freedom: number of branching levels, stem length, and branch lengths. The effect of natural convection in the melt zone is documented.

Optimization of Cooling Protocols for Hearts Destined for Transplantation

Design and analysis of conceptually different cooling systems for the human heart preservation are numerically investigated. A heart cooling container with required connections was designed for a normal size human heart. A three-dimensional, high resolution human heart geometric model obtained from CT-angio data was used for simulations. Nine different cooling designs are introduced in this research. The first cooling design (Case 1) used a cooling gelatin only outside of the heart. In the second cooling design (Case 2), the internal parts of the heart were cooled via pumping a cooling liquid inside both the heart’s pulmonary and systemic circulation systems. An unsteady conjugate heat transfer analysis is performed to simulate the temperature field variations within the heart during the cooling process. Case 3 simulated the currently used cooling method in which the coolant is stagnant. Case 4 was a combination of Case 1 and Case 2. A linear thermoelasticity analysis was performed to assess the stresses applied on the heart during the cooling process. In Cases 5 through 9, the coolant solution was used for both internal and external cooling. For external circulation in Case 5 and Case 6, two inlets and two outlets were designed on the walls of the cooling container. Case 5 used laminar flows for coolant circulations inside and outside of the heart. Effects of turbulent flow on cooling of the heart were studied in Case 6. In Case 7, an additional inlet was designed on the cooling container wall to create a jet impinging the hot region of the heart’s wall. Unsteady periodic inlet velocities were applied in Case 8 and Case 9. The average temperature of the heart in Case 5 was +5.0oC after 1500 s of cooling. Multi-objective constrained optimization was performed for Case 5. Inlet velocities for two internal and one external coolant circulations were the three design variables for optimization. Minimizing the average temperature of the heart, wall shear stress and total volumetric flow rates were the three objectives. The only constraint was to keep von Mises stress below the ultimate tensile stress of the heart’s tissue.

Chemical analysis of water samples of station M1_1207, Red Sea (Table 1-3)

Two water samples and two sediment samples taken in 1965 by the R. V. "Meteor" in the area of the hot salt brine of the Atlantis II-Deep were chemically investigated, and in addition the sediment samples were subjected to X-ray and optical analysis. The investigation of the sulfur-isotope-ratios showed the same values for all water samples. This information combined with the Ca-sulfate solubility data leads us to conclude that, for the most part, the sulfate content of the salt brine resulted from mixing along the boundary with the normal seawater. In this boundary area gypsum or anhydrite is formed which sinks down to the deeper layers of the salt brine where it is redisolved when the water becomes undersaturated. In the laboratory, formation of CaS04 precipitate resulted from both the reheating of the water sample from the uppermost zone of the salt brine to the in-situ-temperature as well as by the mixing of the water sample with normal Red Sea water. The iron and manganese delivered by the hot spring is separated within the area of the salt brine by their different redox-potentials. Iron is sedimented to a high amount within the salt brine, while, as evidenced by its small amounts in all sediment samples, the more easily reducible manganese is apparently carried out of the area before sedimentation can take place. The very good layering of the salt brine may be the result of the rough bottom topography with its several progressively higher levels allowing step-like enlargements of the surface areas of each successive layer. Each enlargement results in larger boundary areas along which more effective heat transfer and mixing with the next layer is possible. In the sediment samples up to 37.18% Fe is found, mostly bound as very poorly crystallized iron hydroxide. Pyrite is present in only very small amounts. We assume that the copper is bound mostly as sulfide, while the zinc is most likely present in an other form. The sulfur-isotope-investigations indicate that the sulfur in the sediment, bound as pyrite and sulfides, is not a result of bacterical sulfate-reduction in the iron-rich mud of the Atlantis II-Deep, but must have been brought up with the hot brine.

(Table 1) Age control and temporal resolution of paleoceanographic records of the North Atlantic

Short-term changes in sea surface conditions controlling the thermohaline circulation in the northern North Atlantic are expected to be especially efficient in perturbing global climate stability. Here we assess past variability of sea surface temperature (SST) in the northeast Atlantic and Norwegian Sea during Marine Isotope Stage (MIS) 2 and, in particular, during the Last Glacial Maximum (LGM). Five high-resolution SST records were established on a meridional transect (53°N-72°N) to trace centennial-scale oscillations in SST and sea-ice cover. We used three independent computational techniques (SIMMAX modern analogue technique, Artificial Neural Networks (ANN), and Revised Analog Method (RAM)) to reconstruct SST from planktonic foraminifer census counts. SIMMAX and ANN reproduced short-term SST oscillations of similar magnitude and absolute levels, while RAM, owing to a restrictive analog selection, appears less suitable for reconstructing "cold end" SST. The SIMMAX and ANN SST reconstructions support the existence of a weak paleo-Norwegian Current during Dansgaard-Oeschger (DO) interstadials number 4, 3, 2, and 1. During the LGM, two warm incursions of 7°C water to occurred in the northern North Atlantic but ended north of the Iceland Faroe Ridge. A rough numerical estimate shows that the near-surface poleward heat transfer from 53° across the Iceland-Faroe Ridge up to to 72° N dropped to less than 60% of the modern value during DO interstadials and to almost zero during DO stadials. Summer sea ice was generally confined to the area north of 70°N and only rarely expanded southward along the margins of continental ice sheets. Internal LGM variability of North Atlantic (>40°N) SST in the GLAMAP 2000 compilation (Sarnthein et al., 2003, doi:10.1029/2002PA000771; Pflaumann et al., 2003, doi:10.1029/2002PA000774) indicates maximum instability in the glacial subpolar gyre and at the Iberian Margin, while in the Nordic Seas, SST was continuously low.

Planktic foraminiferal distribution and stable isotope ratios of sediment core MSM05/5_712-1 from the Arctic Ocean

The Arctic is responding more rapidly to global warming than most other areas on our planet. Northward flowing Atlantic Water is the major means of heat advection towards the Arctic and strongly affects the sea ice distribution. Records of its natural variability are critical for the understanding of feedback mechanisms and the future of the Arctic climate system, but continuous historical records reach back only ~150 years. Here, we present a multidecadal scale record of ocean temperature variations during the last 2000 years, derived from marine sediments off Western Svalbard (79°N). We find that early-21st-century temperatures of Atlantic Water entering the Arctic Ocean are unprecedented over the past 2000 years and are presumably linked to the Arctic Amplification of global warming.

Preparation of Polymer Nanoparticles by Ring Opening Metathesis Polymerization in Aqueous Dispersed Systems with Water-Soluble Ruthenium-based Catalyst

Ring opening metathesis polymerization (ROMP) is a variant of olefin metathesis used to polymerize strained cyclic olefins. Ruthenium-based Grubbs’ catalysts are widely used in ROMP to produce industrially important products. While highly efficient in organic solvents such as dichloromethane and toluene, these hydrophobic catalysts are not typically applied in aqueous systems. With the advancements in emulsion and miniemulsion polymerization, it is promising to conduct ROMP in an aqueous dispersed phase to generate well-defined latex nanoparticles while improving heat transfer and reducing the use of volatile organic solvents (VOCs). Herein I report the efforts made using a PEGylated ruthenium alkylidene as the catalyst to initiate ROMP in an oil-in-water miniemulsion. 1H NMR revealed that the synthesized PEGylated catalyst was stable and reactive in water. Using 1,5-cyclooctadiene (COD) as monomer, we showed the highly efficient catalyst yielded colloidally stable polymer latexes with ~ 100% conversion at room temperature. Kinetic studies demonstrated first-order kinetics with good livingness as confirmed by the shift of gel permeation chromatography (GPC) traces. Depending on the surfactants used, the particle sizes ranged from 100 to 300 nm with monomodal distributions. The more strained cyclic olefin norbornene (NB) could also be efficiently polymerized with a PEGylated ruthenium alkylidene in miniemulsion to full conversion and with minimal coagulum formation.

Parametric Analysis of Concrete Solar Collectors

Concrete solar collectors offer a type of solar collector with structural, aesthetic and economic advantages over current populartechnologies. This study examines the influential parameters of concrete solar collectors. In addition to the external conditions,the performance of a concrete solar collector is influenced by the thermal properties of the concrete matrix, piping network andfluid. Geometric and fluid flow parameters also influence the performance of the concrete solar collector. A literature review ofconcrete solar collectors is conducted in order to define the benchmark parameters from which individual parameters are thencompared. The numerical model consists of a 1D pipe flow network coupled with the heat transfer in a 3D concrete domain. Thispaper is concerned with the physical parameters that define the concrete solar collector, thus a constant surface temperature isused as the exposed surface boundary condition with all other surfaces being insulated. Results show that, of the parametersinvestigated, the pipe spacing, ps, concrete conductivity, kc, and the pipe embedment depth, demb, are among those parameterswhich have greatest effect on the collector’s performance. The optimum balance between these parameters is presented withrespect to the thermal performance and discussed with reference to practical development issues.