Thermal properties of complex suspension flows are investigated using numerical computations. The objective is to develop an efficient and accurate computational method to investigate heat transport in suspension flows. The method presented here is based on solving the lattice Boltzmann equation for the fluid phase, as it is coupled to the Newtonian dynamics equations to model the movement of particles and the energy equation to find the thermal properties. This is a direct numerical simulation that models the free movement of the solid particles suspended in the flow and its effect on the temperature distribution. Parallel implementations are done using MPI (message passing interface) method. Convective heat transfer in internal suspension flow (low solid volume fraction, φ<10%), heat transfer in hot pressing of fiber suspensions and thermal performance of particle filled thermal interface materials (high solid volume fraction, φ>40%) are investigated. The effects of flow disturbance due to movement of suspended particles, thermo-physical properties of suspensions and the particle micro structures are discussed.

Energy-storing electrochemical batteries are the most critical components of high energy density storage systems for stationary and mobile applications. Lithium-ion batteries have received considerable interest for hybrid electric vehicles (HEV) because of their high specific energy, but face inherent thermal management challenges that have not been adequately addressed. In the present investigation, a fully coupled electrochemical and thermal model for lithium-ion batteries is developed to investigate the impact of different thermal management strategies on battery performance. This work represents the first ever study of these coupled electrochemical-thermal phenomena in batteries from the electrochemical heat generation all the way to the dynamic heat removal in actual HEV drive cycles. In contrast to previous modeling efforts focused either exclusively on particle electrochemistry on the one hand or overall vehicle simulations on the other, the present work predicts local electrochemical reaction rates using temperature-dependent data on commercially available batteries designed for high rates (C/LiFePO4) in a computationally efficient manner. Simulation results show that conventional external cooling systems for these batteries, which have a low composite thermal conductivity (~1 W/m-K), cause either large temperature rises or internal temperature gradients. Thus, a novel, passive internal cooling system that uses heat removal through liquid-vapor phase change is developed. Although there have been prior investigations of phase change at the microscales, fluid flow at the conditions expected here is not well understood. A first-principles based cooling system performance model is developed and validated experimentally, and is integrated into the coupled electrochemical-thermal model for assessment of performance improvement relative to conventional thermal management strategies. The proposed cooling system passively removes heat almost isothermally with negligible thermal resistances between the heat source and cooling fluid. Thus, the minimization of peak temperatures and gradients within batteries allow increased power and energy densities unencumbered by thermal limitations.

During spring, it has been observed that several thousand bees leave their hive, and settle on some object such as a tree branch. Some of the scout bees search for a suitable place where a new hive can be set up, while the rest collect together to form a swarm. Heinrich (J. of Exp. Biology 91 (1981) 25; Science 212 (1981) 565; Scientific American 244:6 (1981) 147) has done some experiments with free and captive swarms. His observations are as follows. (1)The core (centre) temperature is around 35°C irrespective of the ambient temperature. (2)The mantle (outer surface) temperature exceeds the ambient temperature by 2- 3°C, provided the ambient temperature is greater than 20°C. Otherwise the mantle temperature is maintained around 17°C. (3) The temperature gradient vanishes just before take-off of the swarm. The present work attempts to predict temperature profiles in swarms and compare them with the data of Heinrich. A continuum model involving unsteady heat conduction and heat generation within the swarm is used. Heat loss from the outer surface of the swarm by free convection and radiation is accounted for approximately. To simplify the analysis, internal convection within the swarm is neglected. The energy balance equation is solved using the finite element method. The effective thermal conductivity (k) is determined by comparing model predictions with data for a swarm of dead bees. The estimated value of k is 0.20 W/m-K. Both spherical and a non-spherical axisymmetric shapes are considered. Considering axisymmetric swarms of live bees, temperature profiles are obtained using various heat generation functions which are available in literature. The effective thermal conductivity is assumed to be the same as that for the swarm of dead bees. Results based on a modified version of Southwick's heat generation function (The Behavior and Physiology of Bees, pp. 28-47, 1991) are qualitatively in accord with the data. The predicted maximum temperature within the swarm and the temperature at the lower surface of the swarm at the ambient temperature of 5°C are 34°C and 17-20°C, respectively. These are comparable to the measured values of 36°C and 19°C. The predicted maximum temperature within the swarm and the temperature at the lower surface of the swarm at the ambient temperature of 9°C are 36.5°C and 17-22°C, respectively. These are comparable to the measured values of 35°C and 19°C. The predicted oxygen consumption rates are 2.55 ml/g/hr for a swarm of 5284 bees at an ambient temperature Ta = 5°C and 1.15 ml/g/hr for 16,600 bees at Ta = 9°C. These are of the same order as the measured values (2 ml/g/hr for 5284 bees at Ta = 4.4DC and 0.45-0.55 ml/g/hr for 5284 bees at Ta = 10°C). Omholt and Lanvik (J. of Theoretical Biology, 120 (1986) 447) assumed a non-uniform steady state profile and used it to estimate the heat generation function. Using this function in the transient energy balance, it is found that their steady state profile is unstable.

Not available / text

The following research contributes original concepts to the fields of deterministic neutron transport modeling and reactor power excursion simulation. A deterministic neutron transport code was created to assess the value of new methods of determining neutron current, fluence, and flux values through the use of view factor and average path length calculations. The neutron transport code is also capable of modeling the highly anisotropic neutron transport of deuterium-tritium fusion external source neutrons using diffusion theory with the aid of a modified first collision source term. The neutron transport code was benchmarked with MCNP, an industry standard stochastic neutron transport code. Deterministic neutron transport methods allow users to model large quantities of neutrons without simulating their interactions individually. Subsequently, deterministic methods allow users to more easily couple neutron transport simulations with other physics simulations. Heat transfer and thermoelastic mechanics physics simulation modules were each developed and benchmarked using COMSOL, a commercial heat transfer and mechanics simulation software. The physics simulation modules were then coupled and used to simulate reactor pulses in fast burst and externally driven nuclear assemblies. The coupled system of equations represents a new method of simulating reactor pulses that allows users to more fully characterize pulsed assemblies. Unlike older methods of reactor pulse simulation, the method presented in this research does not require data from the operational reactor in order to simulate its behavior. The ability to simulate the coupled neutron transport and thermo-mechanical feedback present in pulsed reactors prior their construction would significantly enhance the quality of pulsed reactor pre-construction safety analysis. Additionally, a graphical user interface is created to allow users to run simulations and visualize the results using the coupled physics simulation modules. / text

Not available / text

No description available.

No description available.

This thesis employs computational and experimental methods to explore hotspot cooling and high heat flux removal from a 3-D stacked chip using thermoelectric and microfluidic devices. Stacked chips are expected to improve microelectronics performance, but present severe thermal management challenges. The thesis provides an assessment of both thermoelectric and microfluidic technologies and provides guidance for their implementation in the 3-D stacked chips. A detailed 3-D thermal model of a stacked electronic package with two dies and four ultrathin integrated TECs is developed to investigate the efficacy of TECs in hotspot cooling for 3-D technology. The numerical analysis suggests that TECs can be used for on demand cooling of hotspots in 3-D stacked chip architecture. A strong vertical coupling is observed between the top and bottom TECs and it is found that the bottom TECs can detrimentally heat the top hotspots. As a result, TECs need to be carefully placed inside the package to avoid such undesired heating. Thermal contact resistances between dies, inside the TEC module, and between the TEC and heat spreader are shown to significantly affect TEC performance. TECs are most effective for cooling localized hotspots, but microchannels are advantageous for cooling large background heat fluxes. In the present work, the results of heat transfer and pressure drop experiments in the microchannels with water as the working fluid are presented and compared to the previous microchannel experiments and CFD simulations. Heat removal rates of greater than 100 W/cm2 are demonstrated with these microchannels, with a pressure drop of 75 kPa or less. A novel empirical correlation modeling method is proposed, which uses finite element modeling to model conduction in the channel walls and substrate, coupled with an empirical correlation to determine the convection coefficient. This empirical correlation modeling method is compared to resistor network and CFD modeling. The proposed modeling method produced more accurate results than resistor network modeling, while solving 60% faster than a conjugate heat transfer model using CFD. The results of this work demonstrate that microchannels have the ability to remove high heat fluxes from microelectronic packages using water as a working fluid. Additionally, TECs can locally cool hotspots, but must be carefully placed to avoid undesired heating. Future work should focus on overcoming practical challenges including fabrication, cost, and reliability which are preventing these technologies from being fully leveraged.

In modern microprocessors, thermal management has become one of the main hurdles in continued performance enhancement. Cooling schemes utilizing single phase microfluidics have been investigated extensively for enhanced heat dissipation from microprocessors. However, two-phase fluidic cooling devices are becoming a promising approach, and are less understood. This study aims to examine two-phase flow and heat transfer within a pin-fin enhanced micro-gap. The pin-fin array covered an area of 1cm x 1cm and had a pin diameter, height and pitch of 150μm, 200μm and 225μm, respectively, (aspect ratio of 1.33). This study covers both uniformly and partially heated scenarios. The working fluid used was R245fa. The average heat transfer coefficient and high speed flow visualization results indicated a rapid transition to the annular flow regime with a strong dependence on heat flux. Also, unique, conically-shaped two-phase wakes were observed, demonstrating the lateral spreading capability of the pin-fin array geometry.