Heat dissipation from carbon nano-electronics

The incorporation of graphitic compounds such as carbon nanotubes (CNTs) and graphene into nano-electronic device packaging holds much promise for waste heat management given their high thermal conductivities. However, as these graphitic materials must be used in together with other semiconductor/insulator materials, it is not known how thermal transport is affected by the interaction. Using different simulation techniques, in this thesis, we evaluate the thermal transport properties - thermal boundary conductance (TBC) and thermal conductivity - of CNTs and single-layer graphene in contact with an amorphous SiO2 (a-SiO2) substrate. First, the theoretical methodologies and concepts used in our simulations are presented. In particular, two concepts are described in detail as they are necessary for the understanding of the subsequent chapters. The first is the linear response Green-Kubo (GK) theory of thermal boundary conductance (TBC), which we develop in this thesis, and the second is the spectral energy density method, which we use to directly compute the phonon lifetimes and thermal transport coefficients. After we set the conceptual foundations, the TBC of the CNT-SiO2 interface is computed using non- equilibrium molecular dynamics (MD) simulations and the new Green-Kubo method that we have developed. Its dependence on temperature, the strength of the interaction with the substrate, and tube diameter are evaluated. To gain further insight into the phonon dynamics in supported CNTs, the scattering rates are computed using the spectral energy density (SED) method. With this method, we are able to distinguish the different scattering mechanisms (boundary and CNT-substrate phonon-phonon) and rates. The phonon lifetimes in supported CNTs are found to be reduced by contact with the substrate and we use that lifetime reduction to determine the change in CNT thermal conductivity. Next, we examine thermal transport in graphene supported on SiO2. The phonon contribution to the TBC of the graphene-SiO2 interface is computed from MD simulations and found to agree well with experimentally measured values. We derive the theory of remote phonon scattering of graphene electrons and compute the heat transfer coefficient dependence on doping level and temperature. The thermal boundary conductance from remote phonon scattering is found to be an order of magnitude smaller than that of the phonon contribution. The in-plane thermal conductivity of supported graphene is calculated from MD simulations. The experimentally measured order of magnitude reduction in thermal conductivity is reproduced in our simulations. We show that this reduction is due to the damping of the flexural (ZA) modes. By varying the interaction between graphene and the substrate, the ZA modes hybridize with the substrate Rayleigh modes and the dispersion of the hybridized modes is found to linearize in the strong coupling limit, leading to an increased thermal conductance in the composite structure.

Bomb disposal in the tropics : a cocktail of metabolic and environmental heat

Bomb technicians perform their work while encapsulated in explosive ordnance disposal (EOD) suits. Designed primarily for safety, these suits have an unintended consequence of impairing the body’s natural mechanisms for heat dissipation. Purpose: To quantify the heat strain encountered during an EOD operational scenario in the tropical north of Australia. Methods: All active police male bomb technicians, located in a tropical region of Australia (n=4, experience 7 ± 2.1 yrs, age 34 ± 2 yrs, height 182.3 ± 5.4 cm, body mass 95 ± 4 kg, VO2max 46 ± 5.7 ml.kg-1.min-1) undertook an operational scenario wearing the Med-Eng EOD 9 suit and helmet (~32 kg). The climatic conditions ranged between 27.1–31.8°C ambient temperature, 66-88% relative humidity, and 30.7-34.3°C wet bulb globe temperature. The scenario involved searching a two story non air-conditioned building for a target; carrying and positioning equipment for taking an X-ray; carrying and positioning equipment to disrupt the target; and finally clearing the site. Core temperature and heart rate were continuously monitored, and were used to calculate a physiological strain index (PSI). Urine specific gravity (USG) assessed hydration status and heat associated symptomology were also recorded. Results: The scenario was completed in 121 ± 22 mins (23.4 ± 0.4% work, 76.5 ± 0.4% rest/recovery). Maximum core temperature (38.4 ± 0.2°C), heart rate (173 ± 5.4 bpm, 94 ± 3.3% max), PSI (7.1 ± 0.4) and USG (1.031 ± 0.002) were all elevated after the simulated operation. Heat associated symptomology highlighted that moderate-severe levels of fatigue and thirst were universally experienced, muscle weakness and heat sensations experienced by 75%, and one bomb technician reported confusion and light-headedness. Conclusion: All bomb technicians demonstrated moderate-high levels of heat strain, evidenced by elevated heart rate, core body temperature and PSI. Severe levels of dehydration and noteworthy heat-related symptoms further highlight the risks to health and safety faced by bomb technicians operating in tropical locations.

Numerical Study Of Heat Transfer From Pin Fin Heat Sink Using Steady And Pulsated Impinging Jets

The work reported in this thesis is an attempt to enhance heat transfer in electronic devices with the use of impinging air jets on pin-finned heat sinks. The cooling per-formance of electronic devices has attracted increased attention owing to the demand of compact size, higher power densities and demands on system performance and re-liability. Although the technology of cooling has greatly advanced, the main cause of malfunction of the electronic devices remains overheating. The problem arises due to restriction of space and also due to high heat dissipation rates, which have increased from a fraction of a W/cm2to 100s of W /cm2. Although several researchers have at-tempted to address this at the design stage, unfortunately the speed of invention of cooling mechanism has not kept pace with the ever-increasing requirement of heat re- moval from electronic chips. As a result, efficient cooling of electronic chip remains a challenge in thermal engineering. Heat transfer can be enhanced by several ways like air cooling, liquid cooling, phase change cooling etc. However, in certain applications due to limitations on cost and weight, eg. air borne application, air cooling is imperative. The heat transfer can be increased by two ways. First, increasing the heat transfer coefficient (forced convec- tion), and second, increasing the surface area of heat transfer (finned heat sinks). From previous literature it was established that for a given volumetric air flow rate, jet im-pingement is the best option for enhancing heat transfer coefficient and for a given volume of heat sink material pin-finned heat sinks are the best option because of their high surface area to volume ratio. There are certain applications where very high jet velocities cannot be used because of limitations of noise and presence of delicate components. This process can further be improved by pulsating the jet. A steady jet often stabilizes the boundary layer on the surface to be cooled. Enhancement in the convective heat transfer can be achieved if the boundary layer is broken. Disruptions in the boundary layer can be caused by pulsating the impinging jet, i.e., making the jet unsteady. Besides, the pulsations lead to chaotic mixing, i.e., the fluid particles no more follow well defined streamlines but move unpredictably through the stagnation region. Thus the flow mimics turbulence at low Reynolds number. The pulsation should be done in such a way that the boundary layer can be disturbed periodically and yet adequate coolant is made available. So, that there is not much variation in temperature during one pulse cycle. From previous literature it was found that square waveform is most effective in enhancing heat transfer. In the present study the combined effect of pin-finned heat sink and impinging slot jet, both steady and unsteady, has been investigated for both laminar and turbulent flows. The effect of fin height and height of impingement has been studied. The jets have been pulsated in square waveform to study the effect of frequency and duty cycle. This thesis attempts to increase our understanding of the slot jet impingement on pin-finned heat sinks through numerical investigations. A systematic study is carried out using the finite-volume code FLUENT (Version 6.2) to solve the thermal and flow fields. The standard k-ε model for turbulence equations and two layer zonal model in wall function are used in the problem Pressure-velocity coupling is handled using the SIMPLE algorithm with a staggered grid. The parameters that affect the heat transfer coefficient are: height of the fins, total height of impingement, jet exit Reynolds number, frequency of the jet and duty cycle (percentage time the jet is flowing during one complete cycle of the pulse). From the studies carried out it was found that: a) beyond a certain height of the fin the rate of enhancement of heat transfer becomes very low with further increase in height, b) the heat transfer enhancement is much more sensitive to any changes at low Reynolds number than compared to high Reynolds number, c) for a given total height of impingement the use of fins and pulsated jet, increases the effective heat transfer coefficient by almost 200% for the same average Reynolds number, d) for all the cases it was observed that the optimum frequency of impingement is around 50 − 100 Hz and optimum duty cycle around 25-33.33%, e) in the case of turbulent jets the enhancement in heat transfer due to pulsations is very less compared to the enhancement in case of laminar jets.

A proposed qualitative framework for heterogeneous burning of metallic materials : the 'melting rate triangle'

This paper presents a proposed qualitative framework to discuss the heterogeneous burning of metallic materials, through parameters and factors that influence the melting rate of the solid metallic fuel (either in a standard test or in service). During burning, the melting rate is related to the burning rate and is therefore an important parameter for describing and understanding the burning process, especially since the melting rate is commonly recorded during standard flammability testing for metallic materials and is incorporated into many relative flammability ranking schemes. However, whilst the factors that influence melting rate (such as oxygen pressure or specimen diameter) have been well characterized, there is a need for an improved understanding of how these parameters interact as part of the overall melting and burning of the system. Proposed here is the ‘Melting Rate Triangle’, which aims to provide this focus through a conceptual framework for understanding how the melting rate (of solid fuel) is determined and regulated during heterogeneous burning. In the paper, the proposed conceptual model is shown to be both (a) consistent with known trends and previously observed results, and (b)capable of being expanded to incorporate new data. Also shown are examples of how the Melting Rate Triangle can improve the interpretation of flammability test results. Slusser and Miller previously published an ‘Extended Fire Triangle’ as a useful conceptual model of ignition and the factors affecting ignition, providing industry with a framework for discussion. In this paper it is shown that a ‘Melting Rate Triangle’ provides a similar qualitative framework for burning, leading to an improved understanding of the factors affecting fire propagation and extinguishment.

Natural convection and heat transfer in attics subject to periodic thermal forcing

The heat transfer through the attics of buildings under realistic thermal forcing has been considered in this study. A periodic temperature boundary condition is applied on the sloping walls of the attic to show the basic flow features in the attic space over diurnal cycles. The numerical results reveal that, during the daytime heating stage, the flow in the attic space is stratified; whereas at the night-time cooling stage, the flow becomes unstable. A symmetrical solution is seen for relatively low Rayleigh numbers. However, as the Ra gradually increases, a transition occurs at a critical value of Ra. Above this critical value, an asymmetrical solution exhibiting a pitchfork bifurcation arises at the night-time. It is also found that the calculated heat transfer rate at the night-time cooling stage is much higher than that during the daytime heating stage.

Effect of MHD and heat generation on natural convection flow in an open square cavity under microgravity condition

Purpose - Thermo-magnetic convection and heat transfer of paramagnetic fluid placed in a micro-gravity condition (g = 0) and under a uniform vertical gradient magnetic field in an open square cavity with three cold sidewalls have been studied numerically. Design/methodology/approach - This magnetic force is proportional to the magnetic susceptibility and the gradient of the square of the magnetic induction. The magnetic susceptibility is inversely proportional to the absolute temperature based on Curie’s law. Thermal convection of a paramagnetic fluid can therefore take place even in zero-gravity environment as a direct consequence of temperature differences occurring within the fluid due to a constant internal heat generation placed within a magnetic field gradient. Findings - Effects of magnetic Rayleigh number, Ra, Prandtl number, Pr, and paramagnetic fluid parameter, m, on the flow pattern and isotherms as well as on the heat absorption are presented graphically. It is found that the heat transfer rate is suppressed in increased of the magnetic Rayleigh number and the paramagnetic fluid parameter for the present investigation. Originality/value - It is possible to control the buoyancy force by using the super conducting magnet. To the best knowledge of the author no literature related to magnetic convection for this configuration is available.

MHD natural convection flow from an isothermal horizontal circular cylinder under consideration of temperature dependent viscosity

Magnetohydrodynamic (MHD) natural convection laminar flow from an iso-thermal horizontal circular cylinder immersed in a fluid with viscosity proportional to a linear function of temperature will be discussed with numerical simulations. The governing boundary layer equations are transformed into a non-dimensional form and the resulting nonlinear system of partial differential equa-tions are reduced to convenient form, which are solved numerically by two very efficient methods, namely, (i) Implicit finite difference method together with Keller box scheme and (ii) Direct numerical scheme. Numerical results are presented by velocity and temperature distributions of the fluid as well as heat transfer characteristics, namely the shearing stress and the local heat transfer rate in terms of the local skin-friction coefficient and the local Nusselt number for a wide range of magnetohydrodynamic parameter, viscosity-variation parameter and viscous dissipation parameter. MHD flow in this geometry with temperature dependent viscosity is absent in the literature. The results obtained from the numerical simulations have been veri-fied by two methodologies.

Physiological tolerance times while wearing explosive ordnance disposal protective clothing in simulated environmental extremes

Explosive ordnance disposal (EOD) technicians are required to wear protective clothing to protect themselves from the threat of overpressure, fragmentation, impact and heat. The engineering requirements to minimise these threats results in an extremely heavy and cumbersome clothing ensemble that increases the internal heat generation of the wearer, while the clothing’s thermal properties reduce heat dissipation. This study aimed to evaluate the heat strain encountered wearing EOD protective clothing in simulated environmental extremes across a range of differing work intensities. Eight healthy males [age 25±6 years (mean ± sd), height 180±7 cm, body mass 79±9 kg, V˙O2max 57±6 ml.kg−1.min−1] undertook nine trials while wearing an EOD9 suit (weighing 33.4 kg). The trials involved walking on a treadmill at 2.5, 4 and 5.5 km⋅h−1 at each of the following environmental conditions, 21, 30 and 37°C wet bulb globe temperature (WBGT) in a randomised controlled crossover design. The trials were ceased if the participants’ core temperature reached 39°C, if heart rate exceeded 90% of maximum, if walking time reached 60 minutes or due to fatigue/nausea. Tolerance times ranged from 10–60 minutes and were significantly reduced in the higher walking speeds and environmental conditions. In a total of 15 trials (21%) participants completed 60 minutes of walking; however, this was predominantly at the slower walking speeds in the 21°C WBGT environment. Of the remaining 57 trials, 50 were ceased, due to attainment of 90% maximal heart rate. These near maximal heart rates resulted in moderate-high levels of physiological strain in all trials, despite core temperature only reaching 39°C in one of the 72 trials.

Heat transfer analysis of viscous incompressible fluid by combined natural convection and radiation in an open cavity

The effect of radiation on natural convection of Newtonian fluid contained in an open cavity is investigated in this study. The governing partial differential equations are solved numerically using the Alternate Direct Implicit method together with the Successive Over Relaxation method. The study is focused on studying the flow pattern and the convective and radiative heat transfer rates are studied for different values of radiation parameters namely, the optical thickness of the fluid, scattering albedo, and the Planck number. It was found that in the optically thin limit, an increase in the optical thickness of the fluid raises the temperature and radiation heat transfer of the fluid. However, a further increase in the optical thickness decreases the radiative heat transfer rate due to increase in the energy level of the fluid, which ultimately reduces the total heat transfer rate within the fluid.

Inside the ‘Hurt Locker’: The combined effects of explosive ordnance disposal and chemical protective clothing on physiological tolerance time in extreme environments

Background Explosive ordnance disposal (EOD) technicians are often required to wear specialised clothing combinations that not only protect against the risk of explosion but also potential chemical contamination. This heavy (>35kg) and encapsulating ensemble is likely to increase physiological strain by increasing metabolic heat production and impairing heat dissipation. This study investigated the physiological tolerance times of two different chemical protective undergarments, commonly worn with EOD personal protective clothing, in a range of simulated environmental extremes and work intensities Methods Seven males performed eighteen trials wearing two ensembles. The trials involved walking on a treadmill at 2.5, 4 and 5.5 km.h-1 at each of the following environmental conditions, 21, 30 and 37°C wet bulb globe temperature (WBGT). The trials were ceased if the participants’ core temperature reached 39°C, if heart rate exceeded 90% of maximum, if walking time reached 60 minutes or due to volitional fatigue. Results Physiological tolerance times ranged from 8 to 60 min and the duration (mean difference: 2.78 min, P>0.05) were similar in both ensembles. A significant effect for environment (21>30>37°C WBGT, P<0.05) and work intensity (2.5>4>5.5 km.h-1, P< 0.05) was observed in tolerance time. The majority of trials across both ensembles (101/126; 80.1%) were terminated due to participants achieving a heart rate equivalent to greater than 90% of their maximum. Conclusions Physiological tolerance times wearing these two chemical protective undergarments, worn underneath EOD personal protective clothing, were similar and predominantly limited by cardiovascular strain.

Can perceptual indices estimate physiological strain across a range of environments and metabolic workloads when wearing explosive ordnance disposal and chemical protective clothing?

Objective Explosive ordnance disposal (EOD) often requires technicians to wear multiple protective garments in challenging environmental conditions. The accumulative effect of increased metabolic cost coupled with decreased heat dissipation associated with these garments predisposes technicians to high levels of physiological strain. It has been proposed that a perceptual strain index (PeSI) using subjective ratings of thermal sensation and perceived exertion as surrogate measures of core body temperature and heart rate, may provide an accurate estimation of physiological strain. Therefore, this study aimed to determine if the PeSI could estimate the physiological strain index (PSI) across a range of metabolic workloads and environments while wearing heavy EOD and chemical protective clothing. Methods Eleven healthy males wore an EOD and chemical protective ensemble while walking on a treadmill at 2.5, 4 and 5.5 km·h− 1 at 1% grade in environmental conditions equivalent to wet bulb globe temperature (WBGT) 21, 30 and 37 °C. WBGT conditions were randomly presented and a maximum of three randomised treadmill walking trials were completed in a single testing day. Trials were ceased at a maximum of 60-min or until the attainment of termination criteria. A Pearson's correlation coefficient, mixed linear model, absolute agreement and receiver operating characteristic (ROC) curves were used to determine the relationship between the PeSI and PSI. Results A significant moderate relationship between the PeSI and the PSI was observed [r = 0.77; p < 0.001; mean difference = 0.8 ± 1.1 a.u. (modified 95% limits of agreement − 1.3 to 3.0)]. The ROC curves indicated that the PeSI had a good predictive power when used with two, single-threshold cut-offs to differentiate between low and high levels of physiological strain (area under curve: PSI three cut-off = 0.936 and seven cut-off = 0.841). Conclusions These findings support the use of the PeSI for monitoring physiological strain while wearing EOD and chemical protective clothing. However, future research is needed to confirm the validity of the PeSI for active EOD technicians operating in the field.

Performance assessment in a heat exchanger tube fitted with double counter twisted tape inserts

The present study explored the effects of the double counter twisted tapes on heat transfer and fluid friction characteristics in a heat exchanger tube. The double counter twisted tapes were used as counter-swirl flow generators in the test section. The experiments were performed with double counter twisted tapes of four different twist ratios (y = 1.95, 3.85, 5.92 and 7.75) using air as the testing fluid in a circular tube turbulent flow regime where the Reynolds number was varied from 6950 to 50,050. The experimental results demonstrated that the Nusselt number, friction factor and thermal enhancement efficiency were increased with decreasing twist ratio. The results also revealed that the heat transfer rate in the tube fitted with double counter twisted tape was significantly increased with corresponding increase in pressure drop. In the range of the present work, heat transfer rate and friction factor were obtained to be around 60 to 240% and 91 to 286% higher than those of the plain tube values, respectively. The maximum thermal enhancement efficiency of 1.34 was achieved by the use of double counter twisted tapes at constant blower power. In addition, the empirical correlations for the Nusselt number, friction factor and thermal enhancement efficiency were also developed, based on the experimental data.

Compressible MHD boundary layer in the stagnation region of a sphere

The steady laminar compressible boundary layer flow of an electrically conducting fluid in the stagnation region of a sphere with an applied magnetic field has been studied. The effects of the induced magnetic field, mass transfer, and viscous dissipation have been taken into account. Both isothermal and adiabatic wall conditions have been considered. The governing equations have been solved numerically using a shooting method. The skin friction and heat transfer are found to be strongly affected by the magnetic field, mass transfer, wall temperature and Mach number. It is found that the magnetic field reduces the heat transfer. This is a significant result which can be used in controlling the heat transfer rate. The boundary layer solutions break down as the magnetic parameter tends to a certain critical value.

Computational and experimental studies of fluid flow and heat transfer in a Calandria based reactor

CFD investigations are carried out to study the heat flux and temperature distribution in the calandria using a 3–Dimensional RANS code. Internal flow computations and experimental studies are carried out for a calandria embedded with a matrix of tubes working together as a reactor. Numerical investigations are carried on the Calandria reactor vessel with horizontal inlets and outlets located on top and the bottom to study the flow pattern and the associated temperature distribution. The computations have been carried out to simulate fluid flow and convective heat transfer for assigned near–to working conditions with different moderator injection rates and reacting heat fluxes. The results of computations provide an estimate of the tolerance bands for safe working limits for the heat dissipation at different working conditions by virtue of prediction of the hot spots in the calandria. The isothermal CFD results are validated by a set of experiments on a specially designed scaled model conducted over a range of flows and simulation parameters. The comparison of CFD results with experiments show good agreement.

Effective thermal conductivity of a heat generating rod bundle dissipating heat by natural convection and radiation

A numerical study of conjugate natural convection and surface radiation in a horizontal hexagonal sheath housing 19 solid heat generating rods with cladding and argon as the fill gas, is performed. The natural convection in the sheath is driven by the volumetric heat generation in the solid rods. The problem is solved using the FLUENT CFD code. A correlation is obtained to predict the maximum temperature in the rod bundle for different pitch-to-diameter ratios and heat generating rates. The effective thermal conductivity is related to the heat generation rate, maximum temperature and the sheath temperature. Results are presented for the dimensionless maximum temperature, Rayleigh number and the contribution of radiation with changing emissivity, total wattage and the pitch-to-diameter ratio. In the simulation of a larger system that contains a rod bundle, the effective thermal conductivity facilitates simplified modelling of the rod bundle by treating it as a solid of effective thermal conductivity. The parametric studies revealed that the contribution of radiation can be 38-65% of the total heat generation, for the parameter ranges chosen. Data for critical Rayleigh number above which natural convection comes into effect is also presented. (C) 2011 Elsevier B.V. All rights reserved.