Ventilation effectiveness measures based on heat removal: Part 1. Definitions

The effectiveness of ventilation flows is considered from the perspective of buoyancy (or heat) removal from a space. This perspective is distinct from the standard in which the effectiveness is based on the concentrations of a neutrally buoyant contaminant/passive tracer. Three new measures of effectiveness are proposed based on the ability of a flow to flush buoyancy from a ventilated space. These measures provide estimates of instantaneous and time-averaged effectiveness for the entire space, and local effectiveness at any height of interest. From a generalisation of the latter, a vertical profile of effectiveness is defined. These measures enable quantitative comparisons to be made between different flows and they are applicable when there is a difference in density (as is typical due to temperature differences) between the interior environment and the replacement air. Applications, therefore, include natural ventilation, hybrid ventilation and a range of forced ventilation flows. Finally, we demonstrate how the ventilation effectiveness of a room may be assessed from simple traces of temperature versus time. © 2006 Elsevier Ltd. All rights reserved.

Heat source modelling and natural ventilation efficiency

We compare natural ventilation flows established by a range of heat source distributions at floor level. Both evenly distributed and highly localised line and point source distributions are considered. We demonstrate that modelling the ventilation flow driven by a uniformly distributed heat source is equivalent to the flow driven by a large number of localised sources. A model is developed for the transient flow development in a room with a uniform heat distribution and is compared with existing models for localised buoyancy inputs. For large vent areas the flow driven by localised heat sources reaches a steady state more rapidly than the uniformly distributed case. For small vent areas there is little difference in the transient development times. Our transient model is then extended to consider the time taken to flush a neutrally buoyant pollutant from a naturally ventilated room. Again comparisons are drawn between uniform and localised (point and line) heat source geometries. It is demonstrated that for large vent areas a uniform heat distribution provides the fastest flushing. However, for smaller vent areas, localised heat sources produce the fastest flushing. These results are used to suggest a definition for the term 'natural ventilation efficiency', and a model is developed to estimate this efficiency as a function of the room and heat source geometries. © 2006 Elsevier Ltd. All rights reserved.

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.

Shroud heat transfer measurements from a rotating cavity with an axial throughflow of air

The paper discusses measurements of heat transfer obtained from the inside surface of the peripheral shroud. The experiments were carried out on a rotating cavity, comprising two 0.985-m-dia disks, separated by an axial gap of 0.065 m and bounded at the circumference by a carbon fiber shroud. Tests were conducted with a heated shroud and either unheated or heated disks. When heated, the disks had the same temperature level and surface temperature distribution. Two different temperature distributions were tested; the surface temperature either increased, or decreased with radius. The effects of disk, shroud, and air temperature levels were also studied. Tests were carried out for the range of axial throughflow rates and speeds: 0.0025 ≤ m ≤ 0.2 kg/s and 12.5 ≤ Ω ≤ 125 rad/s, respectively. Measurements were also made of the temperature of the air inside the cavity. The shroud Nusselt numbers are found to depend on a Grashof number, which is defined using the centripetal acceleration. Providing the correct reference temperature is used, the measured Nusselt numbers also show similarity to those predicted by an established correlation for a horizontal plate in air. The heat transfer from the shroud is only weakly affected by the disk surface temperature distribution and temperature level. The heat transfer from the shroud appears to be affected by the Rossby number. A significant enhancement to the rotationally induced free convection occurs in the regions 2 ≤ Ro ≤ 4 and Ro ≥ 20. The first of these corresponds to a region where vortex breakdown has been observed. In the second region, the Rossby number may be sufficiently large for the central throughflow to affect the shroud heat transfer directly. Heating the shroud does not appear to affect the heat transfer from the disks significantly.

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.

Turbine blade tip heat transfer in low speed and high speed flows

In this paper, high and low speed tip flows are investigated for a high-pressure turbine blade. Previous experimental data are used to validate a CFD code, which is then used to study the tip heat transfer in high and low speed cascades. The results show that at engine representative Mach numbers the tip flow is predominantly transonic. Thus, compared to the low speed tip flow, the heat transfer is affected by reductions in both the heat transfer coefficient and the recovery temperature. The high Mach numbers in the tip region (M>1.5) lead to large local variations in recovery temperature. Significant changes in the heat transfer coefficient are also observed. These are due to changes in the structure of the tip flow at high speed. At high speeds, the pressure side corner separation bubble reattachment occurs through supersonic acceleration which halves the length of the bubble when the tip gap exit Mach number is increased from 0.1 to 1.0. In addition, shock/boundary-layer interactions within the tip gap lead to large changes in the tip boundary-layer thickness. These effects give rise to significant differences in the heat-transfer coefficient within the tip region compared to the low-speed tip flow. Compared to the low speed tip flow, the high speed tip flow is much less dominated by turbulent dissipation and is thus less sensitive to the choice of turbulence model. These results clearly demonstrate that blade tip heat transfer is a strong function of Mach number, an important implication when considering the use of low speed experimental testing and associated CFD validation in engine blade tip design. Copyright © 2009 by ASME.

Decay heat in fast reactors with transuranic fuels

In this work, we performed an evaluation of decay heat power of advanced, fast spectrum, lead and molten salt-cooled reactors, with flexible conversion ratio. The decay heat power was calculated using the BGCore computer code, which explicitly tracks over 1700 isotopes in the fuel throughout its burnup and subsequent decay. In the first stage, the capability of the BGCore code to accurately predict the decay heat power was verified by performing a benchmark calculation for a typical UO2 fuel in a Pressurized Water Reactor environment against the (ANSI/ANS-5.1-2005, "Decay Heat Power in Light Water Reactors," American National Standard) standard. Very good agreement (within 5%) between the two methods was obtained. Once BGCore calculation capabilities were verified, we calculated decay power for fast reactors with different coolants and conversion ratios, for which no standard procedure is currently available. Notable differences were observed for the decay power of the advanced reactor as compared with the conventional UO2 LWR. The importance of the observed differences was demonstrated by performing a simulation of a Station Blackout transient with the RELAP5 computer code for a lead-cooled fast reactor. The simulation was performed twice: using the code-default ANS-79 decay heat curve and using the curve calculated specifically for the studied core by BGCore code. The differences in the decay heat power resulted in failure to meet maximum cladding temperature limit criteria by ∼100 °C in the latter case, while in the transient simulation with the ANS-79 decay heat curve, all safety limits were satisfied. The results of this study show that the design of new reactor safety systems must be based on decay power curves specific to each individual case in order to assure the desired performance of these systems. © 2009 Elsevier B.V. All rights reserved.