LES of heat transfer in electronics

Numerical methods based on the Reynolds Averaged Navier-Stokes (RANS) and Large Eddy Simulation (LES) equations are applied to the thermal prediction of flows representative of those found in and around electronics systems and components. Low Reynolds number flows through a heated ribbed channel, around a heated cube and within a complex electronics system case are investigated using linear and nonlinear LES models, hybrid RANS-LES and RANS-Numerical-LES (RANS-NLES) methods. Flow and heat transfer predictions using these techniques are in good agreement with each other and experimental data for a range of grid resolutions. Using second order central differences, the RANS-NLES method performs well for all simulations. © 2011 Elsevier Inc.

Finite element analysis of compressor disc heat transfer using the transient conduction method with sparse boundary collocation points

Surface temperature measurements from two discs of a gas turbine compressor rig are used as boundary conditions for the transient conduction solution (inverse heat transfer analysis). The disc geometry is complex, and so the finite element method is used. There are often large radial temperature gradients on the discs, and the equations are therefore solved taking into account the dependence of thermal conductivity on temperature. The solution technique also makes use of a multigrid algorithm to reduce the solution time. This is particularly important since a large amount of data must be analyzed to obtain correlations of the heat transfer. The finite element grid is also used for a network analysis to calculate the radiant heat transfer in the cavity formed between the two compressor discs. The work discussed here proved particularly challenging as the disc temperatures were only measured at four different radial locations. Four methods of surface temperature interpolation are examined, together with their effect on the local heat fluxes. It is found that the choice of interpolation method depends on the available number of data points. Bessel interpolation gives the best results for four data points, whereas cubic splines are preferred when there are considerably more data points. The results from the analysis of the compressor rig data show that the heat transfer near the disc inner radius appears to be influenced by the central throughflow. However, for larger radii, the heat transfer from the discs and peripheral shroud is found to be consistent with that of a buoyancy-induced flow.

A novel theoretical approach to passive control of thermo-acoustic oscillations: Application to ducted heat sources

In this paper, we develop a linear technique that predicts how the stability of a thermo-acoustic system changes due to the action of a generic passive feedback device or a generic change in the base state. From this, one can calculate the passive device or base state change that most stabilizes the system. This theoretical framework, based on adjoint equations, is applied to two types of Rijke tube. The first contains an electrically-heated hot wire and the second contains a diffusion flame. Both heat sources are assumed to be compact so that the acoustic and heat release models can be decoupled. We find that the most effective passive control device is an adiabatic mesh placed at the downstream end of the Rijke tube. We also investigate the effects of a second hot wire and a local variation of the cross-sectional area but find that both affect the frequency more than the growth rate. This application of adjoint sensitivity analysis opens up new possibilities for the passive control of thermo-acoustic oscillations. For example, the influence of base state changes can be combined with other constraints, such as that the total heat release rate remains constant, in order to show how an unstable thermo-acoustic system should be changed in order to make it stable. Copyright © 2013 by ASME.

Transient Ventilation Dynamics Following a Change in Strength of a Point Source of Heat

We investigate the transient ventilation flow within a confined ventilated space, with high- and low-level openings, when the strength of a low-level point source of heat is changed instantaneously. The steady-flow regime in the space involves a turbulent buoyant plume, which rises from the point source to a well-mixed warm upper layer. The steady-state height of the interface between this layer and the lower layer of exterior fluid is independent of the heat flux, but the upper layer becomes progressively warmer with heat flux. New analogue laboratory experiments of the transient adjustment between steady states identify that if the heat flux is increased, the continuing plume propagates to the top of the room forming a new, warmer layer. This layer gradually deepens, and as the turbulent plume entrains fluid from the original warm layer, the original layer is gradually depleted and disappears, and a new steady state is established. In contrast, if the source buoyancy flux is decreased, the continuing plume is cooler than the original plume, so that on reaching the interface it is of intermediate density between the original warm layer and the external fluid. The plume supplies a new intermediate layer, which gradually deepens with the continuing flow. In turn, the original upper layer becomes depleted, both as a result of being vented through the upper opening of the space, but also due to some penetrative entrainment of this layer by the plume, as the plume overshoots the interface before falling back to supply the new intermediate layer. We develop quantitative models which are in good accord with our experimental data, by combining classical plume theory with models of the penetrative entrainment for the case of a decrease in heating. Typically, we find that the effect of penetrative entrainment on the density of the intruding layer is relatively weak, provided the change in source strength is sufficiently large. However, penetrative entrainment measurably increases the rate at which the depth of the draining layer decreases. We conclude with a discussion of the importance of these results for the control of naturally ventilated spaces.

Transient Natural Ventilation of a Room with a Distributed Heat Source

We report on an experimental and theoretical study of the transient flows which develop as a naturally ventilated room adjusts from one temperature to another. We focus on a room heated from below by a uniform heat source, with both high- and low-level ventilation openings. Depending on the initial temperature of the room relative to (i) the final equilibrium temperature and (ii) the exterior temperature, three different modes of ventilation may develop. First, if the room temperature lies between the exterior and the equilibrium temperature, the interior remains well-mixed and gradually heats up to the equilibrium temperature. Secondly, if the room is initially warmer than the equilibrium temperature, then a thermal stratification develops in which the upper layer of originally hot air is displaced upwards by a lower layer of relatively cool inflowing air. At the interface, some mixing occurs owing to the effects of penetrative convection. Thirdly, if the room is initially cooler than the exterior, then on opening the vents, the original air is displaced downwards and a layer of ambient air deepens from above. As this lower layer drains, it is eventually heated to the ambient temperature, and is then able to mix into the overlying layer of external air, and the room becomes well-mixed. For each case, we present new laboratory experiments and compare these with some new quantitative models of the transient flows. We conclude by considering the implications of our work for natural ventilation of large auditoria.

On the transition from displacement to mixing ventilation with a localized heat source

We investigate the steady state natural ventilation of a room heated at the base and consisting of two vents at different levels. We explore how the air flow rate and internal temperature relative to the exterior vary as a function of the vent areas, position of the vents and heat load in order to establish appropriate ventilation strategies for a room. When the room is heated by a distributed source, the room becomes well mixed and the steady state ventilation rate depends on the heating rate, the area of the vents and the distance between the lower and upper level vents. However, when the room is heated by a localised source the room becomes stratified. If the effective ventilation area is sufficiently large, then the interface separating the two layers lies above the inlet vent and the lower layer is comprised of ambient fluid. In this case the upper layer is warmer than in the well mixed case and the ventilation rate is smaller. However, if the effective area for ventilation is sufficiently small, then the interface separating the two layers lies below the inlet vent and the lower layer is comprised of warm fluid which originates as the cold incoming fluid mixes during descent from the vent through the upper layer. In this case both the ventilation rate and the upper layer temperature are the same as in the case of a distributed heat load. As the vertical separation between lower and upper level vents decreases, then the temperature difference between the layers falls to zero and the room becomes approximately well mixed. These findings suggest how the appropriate ventilation strategy for a room can be varied depending on the exterior temperature, with mixing ventilation more suitable for winter conditions and displacement ventilation for warmer external temperatures.

Heat release rate measurement in turbulent flames

The dominant industrial approach for the reduction of NO x emissions in industrial gas turbines is the lean pre-mixed prevaporized concept. The main advantage of this concept is the lean operation of the combustion process; this decreases the heat release rate from the flame and results in a reduction in operating temperature. The direct measurement of heat release rates via simultaneous laser induced fluorescence of OH and CH 2O radicals using planar laser induced fluorescence. The product of the two images correlated with the forward production rate of the HCO radical, which in turn has correlated well with heat release rates from premixed hydrocarbon flames. The experimental methodology of the measurement of heat release rate and applications in different turbulent premixed flames were presented. This is an abstract of a paper presented at the 7th World Congress of Chemical Engineering (Glasgow, Scotland 7/10-14/2005).

Convective heat dissipation with lattice-frame materials

This paper presents experimental results on heat transfer and pressure drop for a compact heat sink made of fully triangulated, lightweight (porosity∼0.938), aluminum lattice-frame materials (LFMs). Due to the inherent structural anisotropy of the LFMs, two mutually perpendicular orientations were selected for the measurements. Constant heat flux was applied to the heat sink under steady state conditions, and dissipated by forced air convection. The experimental data were compared with those predicted from an analytical model based on fin analogy. The experimental results revealed that pressure drop is strongly dependent upon the orientation of the structure, due mainly to the flow blockage effect. For heat transfer measurements, typical local temperature distributions on the substrate under constant heat flux conditions were captured with infrared camera. The thermal behavior of LFMs was found to follow closely that of cylinder banks, with early transition Reynolds number (based on strut diameter) equal to about 300. The Nusselt number prediction from the fin-analogy correlates well with experimental measurements, except at low Reynolds numbers where a slightly underestimation is observed. Comparisons with empty channels and commonly used heat exchanger media show that the present LFM heat sink can remove heat approximately seven times more efficient than an empty channel and as efficient as a bank of cylinders at the same porosity level. The aluminum LFMs are extremely stiff and strong, making them ideal candidates for multifunctional structures requiring both heat dissipation and mechanical load carrying capabilities. © 2003 Elsevier Ltd. All rights reserved.

Modelling soot formation in a DISI engine

In this work, the formation of soot in a Direct Injection Spark Ignition (DISI) engine is simulated using the Stochastic Reactor Model (SRM) engine code. Volume change, convective heat transfer, turbulent mixing, direct injection and flame propagation are accounted for. In order to simulate flame propagation, the cylinder is divided into an unburned, entrained and burned zone, with the rate of entrainment being governed by empirical equations but combustion modelled with chemical kinetics. The model contains a detailed chemical mechanism as well as a highly detailed soot formation model, however computation times are relatively short. The soot model provides information on the morphology and chemical composition of soot aggregates along with bulk quantities, including soot mass, number density, volume fraction and surface area. The model is first calibrated by simulating experimental data from a Gasoline Direct Injection (GDI) Spark Ignition (SI) engine. The model is then used to simulate experimental data from the literature, where the numbers, sizes and derived mass particulate emissions from a 1.83 L, 4-cylinder, 4 valve production DISI engine were examined. Experimental results from different injection and spark timings are compared with the model and the qualitative trends in aggregate size distribution and emissions match the exhaust gas measurements well. © 2010 The Combustion Institute. Published by Elsevier Inc. All rights reserved.