Analytical solutions for heat transfer during cyclic melting and freezing of a phase change material used in electronic or electrical packaging

In this paper we develop an analytical heat transfer model, which is capable of analyzing cyclic melting and solidification processes of a phase change material used in the context of electronics cooling systems. The model is essentially based on conduction heat transfer, with treatments for convection and radiation embedded inside. The whole solution domain is first divided into two main sub-domains, namely, the melting sub-domain and the solidification sub-domain. Each sub-domain is then analyzed for a number of temporal regimes. Accordingly, analytical solutions for temperature distribution within each subdomain are formulated either using a semi-infinity consideration, or employing a method of quasi-steady state, depending on the applicability. The solution modules are subsequently united, leading to a closed-form solution for the entire problem. The analytical solutions are then compared with experimental and numerical solutions for a benchmark problem quoted in the literature, and excellent agreements can be observed.

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.

Down-to-earth temperatures: The mechanics of the thermal environment

This lecture describes some recent attempts at unravelling the mechanics of the temperature distribution near ground, especially during calm, clear nights. In particular, a resolution is offered of the so-called Ramdas paradox, connected with observations of a temperature minimum some decimetres above bare soil on calm clear nights, in apparent defiance of the Rayleigh criterion for instability due to thermal convection. The dynamics of the associated temperature distribution is governed by radiative and convective transport and by thermal conduction, and is characterised by two time constants, involving respectively quick radiative adjustments and slow diffusive relaxation. The theory underlying the work described here suggests that surface parameters like ground emissivity and soil thermal conductivity can exert appreciable influence on the development of nocturnal inversions.

Theoretical estimation of electromigration in metallic carbon nanotubes considering self-heating effect

In this paper, we estimate the solution of the electromigration diffusion equation (EMDE) in isotopically pure and impure metallic single-walled carbon nanotubes (CNTs) (SWCNTs) by considering self-heating. The EMDE for SWCNT has been solved not only by invoking the dependence of the electromigration flux on the usual applied static electric field across its two ends but also by considering a temperature-dependent thermal conductivity (κ) which results in a variable temperature distribution (T) along its length due to self-heating. By changing its length and isotopic impurity, we demonstrate that there occurs a significant deviation in the SWCNT electromigration performance. However, if κ is assumed to be temperature independent, the solution may lead to serious errors in performance estimation. We further exhibit a tradeoff between length and impurity effect on the performance toward electromigration. It is suggested that, to reduce the vacancy concentration in longer interconnects of few micrometers, one should opt for an isotopically impure SWCNT at the cost of lower κ, whereas for comparatively short interconnects, pure SWCNT should be used. This tradeoff presented here can be treated as a way for obtaining a fairly well estimation of the vacancy concentration and mean time to failure in the bundles of CNT-based interconnects. © 2012 IEEE.

Natural convection in a series of thermally interacting cavities with wall mounted heaters

Laminar natural convection in a series of thermally interacting cavities is numerically studied. Each cavity consists of a conducting bottom wall with a surface mounted heater. The side walls of the cavities are isothermally cooled. Each cavity thermally interacts with its adjacent cavities through the conducting walls. Flow and heat transfer characteristics are studied in detail for various Rayleigh numbers. The convection characteristics in multiple cavities are compared with those in single independent cavity. The thermal interaction between the cavities results in lower temperatures compared with those in independent cavities. While heat is rejected into the adjacent upper cavity through some portion of the conducting wall, heat is received from the adjacent cavity through the remaining portion of the wall. The influence of substrate conductivity on heat exchange between adjacent cavities are examined. Substrate conductivity shows strong effect on temperature distribution. When cooling at both vertical sides is changed to one side cooling, the heat transfer characteristics are changed drastically and many interesting flow features are observed. Effects of cavity aspect ratio is studied and higher heat transfer rates are observed at higher aspect ratios. Correlations for dimensionless temperature maximum and average Nusselt number are presented in terms of Rayleigh number.

Eulerian two-phase flow simulation and experimental validation of semisolid slurry generation process using cooling slope

Experimental and numerical studies of slurry generation using a cooling slope are presented in the paper. The slope having stainless steel body has been designed and constructed to produce semisolid A356 Al alloy slurry. The pouring temperature of molten metal, slope angle of the cooling slope and slope wall temperature were varied during the experiment. A multiphase numerical model, considering liquid metal and air, has been developed to simulate the liquid metal flow along the cooling channel using an Eulerian two-phase flow approach. Solid fraction evolution of the solidifying melt is tracked at different locations of the cooling channel following Schiel's equation. The continuity, momentum and energy equations are solved considering thin wall boundary condition approach. During solidification of the melt, based on the liquid fraction and latent heat of the alloy, temperature of the alloy is modified continuously by introducing a modified temperature recovery method. Numerical simulations has been carried out for semisolid slurry formation by varying the process parameters such as angle of the cooling slope, cooling slope wall temperature and melt superheat temperature, to understand the effect of process variables on cooling slope semisolid slurry generation process such as temperature distribution, velocity distribution and solid fraction of the solidifying melt. Experimental validation performed for some chosen cases reveals good agreement with the numerical simulations.

Model-free simulations for compressible mixing layer

Confined supersonic mixing layer is explored through model-free simulations. Both two- and three-dimensional spatio-temporal simulations were carried out employing higher order finite difference scheme as well as finite volume scheme based on open source software (OpenFOAM) to understand the effect of three-dimensionality on the development of mixing layer. It is observed that although the instantaneous structures exhibit three-dimensional features, the average pressure and velocities are predominantly two-dimensional. The computed wall pressures match well with experimental results fairly well, although three-dimensional simulation underpredicts the wall pressure in the downstream direction. The self-similarity of the velocity profiles is obtained within the duct length for all the simulations. Although the mixing layer thicknesses differ among different simulations, their growth rate is nearly the same. Significant differences are observed for species and temperature distribution between two- and three-dimensional calculations, and two-dimensional calculations do not match the experimental observation of smooth variations in species mass fraction profiles as reported in literature. Reynolds stress distribution for three-dimensional calculations show profiles with less peak values compared to two-dimensional calculations; while normal stress anisotropy is higher for three-dimensional case.

A model of coupled thermal, mechanical, and electrostatic field effects in III-N thin film heterostructures

A one-dimensional coupled multi-physics based model has been developed to accurately compute the effects of electrostatic, mechanical, and thermal field interactions on the electronic energy band structure in group III-nitrides thin film heterostructures. Earlier models reported in published literature assumes electro-mechanical field with uniform temperature thus neglecting self-heating. Also, the effects of diffused interface on the energy band structure were not studied. We include these effects in a self-consistent manner wherein the transport equation is introduced along with the electro-mechanical models, and the lattice structural variation as observed in experiments are introduced at the interface. Due to these effects, the electrostatic potential distribution in the heterostructure is altered. The electron and hole ground state energies decrease by 5% and 9%, respectively, at a relative temperature of 700 K, when compared with the results obtained from the previously reported electro-mechanical model assuming constant and uniform temperature distribution. A diffused interface decreases the ground state energy of electrons and holes by about 11% and 9%, respectively, at a relative temperature of 700 K when compared with the predictions based on uniform temperature based electro-mechanical model. (C) 2013 AIP Publishing LLC.

Analytical solutions for movement of cold water thermal front in a heterogeneous geothermal reservoir

In the present study an analytical model has been presented to describe the transient temperature distribution and advancement of the thermal front generated due to the reinjection of heat depleted water in a heterogeneous geothermal reservoir. One dimensional heat transport equation in porous media with advection and longitudinal heat conduction has been solved analytically using Laplace transform technique in a semi infinite medium. The heterogeneity of the porous medium is expressed by the spatial variation of the flow velocity and the longitudinal effective thermal conductivity of the medium. A simpler solution is also derived afterwards neglecting the longitudinal conduction depending on the situation where the contribution to the transient heat transport phenomenon in the porous media is negligible. Solution for a homogeneous aquifer with constant values of the rock and fluid parameters is also derived with an aim to compare the results with that of the heterogeneous one. The effect of some of the parameters involved, on the transient heat transport phenomenon is assessed by observing the variation of the results with different magnitudes of those parameters. Results prove the heterogeneity of the medium, the flow velocity and the longitudinal conductivity to have great influence and porosity to have negligible effect on the transient temperature distribution. (C) 2013 Elsevier Inc. All rights reserved.

Studies on rheocasting using cooling slope

In the present work, a cooling channel is employed to produce semi-solid A356 alloy slurry. To understand the transport process involved, a 3D non-isothermal, multiphase volume averaging model has been developed for simulation of the semi-solid slurry generation process in the cooling channel. For simulation purpose, the three phases considered are the parent melt, the nearly spherical grains and air as separated but highly coupled interpenetrating continua. The conservation equations of mass, momentum, energy and species have been solved for each phase and the thermal and mechanical interactions (drag force) among the phases have been considered using appropriate model. The superheated liquid alloy is poured at the top of the cooling slope/channel, where specified velocity inlet boundary condition is used in the model, and allowed to flow along gravity through the channel. The melt loses its superheat and becomes semisolid up to the end of cooling channel due to the evolving -Al grains with decreasing temperature. The air phase forms a definable air/liquid melt interface, i.e. free surface, due its low density. The results obtained from the present model includes volume fractions of three different phases considered, grain evolution, grain growth rate, size and distribution of solid grains. The effect of key process variables such as pouring temperature, slope angle of the cooling channel and cooling channel wall temperature on temperature distribution, velocity distribution, grain formation and volume fraction of different phases are also studied. The results obtained from the simulations are validated by microstructure study using SEM and quantitative image analysis of the semi-solid slurry microstructure obtained from the experimental set-up.

Long way to go: how outflows from large galaxies propagate through the hot halo gas

Using hydrodynamic simulations, we study the mass-loss due to supernova-driven outflows from Milky Way type disc galaxies, paying particular attention to the effect of the extended hot halo gas. We find that the total mass-loss at inner radii scales roughly linearly with total mass of stars formed, and that the mass loading factor at the virial radius can be several times its value at inner radii because of the swept up hot halo gas. The temperature distribution of the outflowing material in the inner region (similar to 10 kpc) is bimodal in nature, peaking at 10(5) K and 10(6.5) K, responsible for optical and X-ray emission, respectively. The contribution of cold/warm gas with temperature <= 10(5.5) K to the outflow rate within 10 kpc is approximate to 0.3-0.5. The warm mass loading factor, eta(3e5) (T <= 3 x 10(5) K) is related to the mass loading factor at the virial radius (eta(v)) as eta(v) approximate to 25 eta(3e5) (SFR/M-circle dot yr(-1))(-0.15) for a baryon fraction of 0.1 and a starburst period of 50 Myr. We also discuss the effect of multiple bursts that are separated by both short and long periods. The outflow speed at the virial radius is close to the sound speed in the hot halo, less than or similar to 200 km s(-1). We identify two `sequences' of outflowing cold gas at small scales: a fast (approximate to 500 km s(-1)) sequence, driven by the unshocked free-wind; and a slow sequence (approximate to +/- 100 km s(-1)) at the conical interface of the superwind and the hot halo.

COUPLED OPTICAL-THERMAL-FLUID MODELING OF A DIRECTLY HEATED TUBULAR SOLAR RECEIVER FOR SUPERCRITICAL CO2 BRAYTON CYCLE

Recent studies have evaluated closed-loop supercritical carbon dioxide (s-CO2) Brayton cycles to be a higher energy density system in comparison to conventional superheated steam Rankine systems. At turbine inlet conditions of 923K and 25 MPa, high thermal efficiency (similar to 50%) can be achieved. Achieving these high efficiencies will make concentrating solar power (CSP) technologies a competitive alternative to current power generation methods. To incorporate a s-CO2 Brayton power cycle in a solar power tower system, the development of a solar receiver capable of providing an outlet temperature of 923 K (at 25 MPa) is necessary. The s-CO2 will need to increase in temperature by similar to 200 K as it passes through the solar receiver to satisfy the temperature requirements of a s-CO2 Brayton cycle with recuperation and recompression. In this study, an optical-thermal-fluid model was developed to design and evaluate a tubular receiver that will receive a heat input similar to 2 MWth from a heliostat field. The ray-tracing tool SolTrace was used to obtain the heat-flux distribution on the surfaces of the receiver. Computational fluid dynamics (CFD) modeling using the Discrete Ordinates (DO) radiation model was used to predict the temperature distribution and the resulting receiver efficiency. The effect of flow parameters, receiver geometry and radiation absorption by s-CO2 were studied. The receiver surface temperatures were found to be within the safe operational limit while exhibiting a receiver efficiency of similar to 85%.

Study of Gas Combustion in a Vented Cylindrical Vessel

Complicated interaction of a flame front with a turbulent flow induced by venting is studied during combustion of the stoichiometric propane/air mixture in a relatively large vented cylindrical vessel. Flame position, its shape, and combustion pressure were measured as a function of time and vent parameters. The experimental data were used to verify numerical simulation of the combustion process. The proposed numerical model satisfactorily simulates the main features of combustion in a closed and vented vessel such as flame configuration, flow and temperature fields, and pressure variation pattern. Simulated velocity and temperature distribution are very useful pieces of information because they are not available from experiments.

Thermocapillary flow in an annular liquid layer painted on a moving fiber

In the present paper, a theoretical model is studied on the flow in the liquid annular film, which is ejected from a vessel with relatively higher temperature and painted on the moving solid fiber. A temperature gradient, driving a thermocapillary flow, is formed on the free surface because of the heat transfer from the liquid with relatively higher temperature to the environmental gas with relatively lower temperature. The thermocapillary flow may change the radii profile of the liquid film. This process analyzed is based on the approximations of lubrication theory and perturbation theory, and the equation of the liquid layer radii and the process of thermal hydrodynamics in the liquid layer are solved for a temperature distribution on the solid fiber.

Numerical simulation of magnetic field design for damping thermocapillary convection in a floating half-zone

The magnetic fields produced by electrical coils are designed for damping the the thermocapillary convection in a floating half-zone in microgravity. The fields are designed specially to reduce the flow near the free surface and then in the melt zone by adjusting the longitudinal coil positions close to the melt zone. The effects of the designed magnetic fields on reducing the flow velocity and temperature distribution non-uniformity in the melt zone are stronger than those of the case of an uniform longitudinal magnetic field obtained by numerical simulation, particularly at the melt-rod interface. It brings fundamental insights into the heat and mass transfer control at the solidification interface by the magnetic field design for crystal growth by the floating full-zone method.