Risk based lifetime costs assessment of a ground source heat pump (GSHP) system design: Methodology and case study

Space heating accounts for a large portion of the world's carbon dioxide emissions. Ground Source Heat Pumps (GSHPs) are a technology which can reduce carbon emissions from heating and cooling. GSHP system performance is however highly sensitive to deviation from design values of the actual annual energy extraction/rejection rates from/to the ground. In order to prevent failure and/or performance deterioration of GSHP systems it is possible to incorporate a safety factor in the design of the GSHP by over-sizing the ground heat exchanger (GHE). A methodology to evaluate the financial risk involved in over-sizing the GHE is proposed is this paper. A probability based approach is used to evaluate the economic feasibility of a hypothetical full-size GSHP system as compared to four alternative Heating Ventilation and Air Conditioning (HVAC) system configurations. The model of the GSHP system is developed in the TRNSYS energy simulation platform and calibrated with data from an actual hybrid GSHP system installed in the Department of Earth Science, University of Oxford, UK. Results of the analysis show that potential savings from a full-size GSHP system largely depend on projected HVAC system efficiencies and gas and electricity prices. Results of the risk analysis also suggest that a full-size GSHP with auxiliary back up is potentially the most economical system configuration. © 2012 Elsevier Ltd.

The influence of stacks on flow patterns and stratification associated with natural ventilation

We investigate the steady state natural ventilation of an enclosed space in which vent A, located at height hA above the floor, is connected to a vertical stack with a termination at height H, while the second vent, B, at height hB above the floor, connects directly to the exterior. We first examine the flow regimes which develop with a distributed source of heating at the base of the space. If hBhB>hA, then two different flow regimes may develop. Either (i) there is inflow through vent B and outflow through vent A, or (ii) the flow reverses, with inflow down the stack into vent A and outflow through vent B. With inflow through vent A, the internal temperature and ventilation rate depend on the relative height of the two vents, A and B, while with inflow through vent B, they depend on the height of vent B relative to the height of the termination of the stack H. With a point source of heating, a similar transition occurs, with a unique flow regime when vent B is lower than vent A, and two possible regimes with vent B higher than vent A. In general, with a point source of buoyancy, each steady state is characterised by a two-layer density stratification. Depending on the relative heights of the two vents, in the case of outflow through vent A connected to the stack, the interface between these layers may lie above, at the same level as or below vent A, leading to discharge of either pure upper layer, a mixture of upper and lower layer, or pure lower layer fluid. In the case of inflow through vent A connected to the stack, the interface always lies below the outflow vent B. Also, in this case, if the inflow vent A lies above the interface, then the lower layer becomes of intermediate density between the upper layer and the external fluid, whereas if the interface lies above the inflow vent A, then the lower layer is composed purely of external fluid. We develop expressions to predict the transitions between these flow regimes, in terms of the heights and areas of the two vents and the stack, and we successfully test these with new laboratory experiments. We conclude with a discussion of the implications of our results for real buildings.

A decoupled whole-building simulation engine for rapid exhaustive search of low-carbon and low-energy building refurbishment options

This paper presents the development of a new building physics and energy supply systems simulation platform. It has been adapted from both existing commercial models and empirical works, but designed to provide expedient exhaustive simulation of all salient types of energy- and carbon-reducing retrofit options. These options may include any combination of behavioural measures, building fabric and equipment upgrades, improved HVAC control strategies, or novel low-carbon energy supply technologies. We provide a methodological description of the proposed model, followed by two illustrative case studies of the tool when used to investigate retrofit options of a mixed-use office building and primary school in the UK. It is not the intention of this paper, nor would it be feasible, to provide a complete engineering decomposition of the proposed model, describing all calculation processes in detail. Instead, this paper concentrates on presenting the particular engineering aspects of the model which steer away from conventional practise. © 2011 Elsevier Ltd.

Transient flows in enclosures: A generalised approach for modelling the effects of geometry, heat gains and wind

Transient flows in a confined ventilated space induced by a buoyancy source of time-varying strength and an external wind are examined. The space considered has varying cross-sectional area with height. A generalised theoretical model is proposed to investigate the flow dynamics following the activation of an external wind and an internal source of buoyancy. To investigate the effect of geometry, we vary the angle of the wall inclination of a particular geometry in which a point source of constant buoyancy is activated in the absence of wind. Counter-intuitively the ventilation is worse and lower airflow rates are established for geometries of increasing cross-sectional areas with height. We investigate the effect of the source buoyancy strength by comparing two cases: (1) when the buoyancy input is constant and (2) when the buoyancy input gradually increases over time so that after a finite time the total buoyancy inputs for (1) and (2) are identical. The rate at which the source heat gains are introduced has a significant role on the flow behaviour as we find that, in case (2), a warmer layer and a more pronounced overshoot are obtained than in case (1). The effect of assisting and opposing wind on the transient ventilation of an enclosure of constant cross-sectional area with height and constant heat gains is examined. A Froude number Fr is used to define the relative strengths of the buoyancy-induced and wind-induced velocities and five different transient states and their associated critical Fr are identified. © 2010 Elsevier Ltd.

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.

Thermal stratification produced by plumes and jets in enclosed spaces

The airflow and thermal stratification produced by a localised heat source located at floor level in a closed room is of considerable practical interest and is commonly referred to as a 'filling box'. In rooms with low aspect ratios H/R ≲ 1 (room height H to characteristic horizontal dimension R) the thermal plume spreads laterally on reaching the ceiling and a descending horizontal 'front' forms separating a stably stratified, warm upper region from cooler air below. The stratification is well predicted for H/R ≲ 1 by the original filling box model of Baines and Turner (J. Fluid. Mech. 37 (1968) 51). This model represents a somewhat idealised situation of a plume rising from a point source of buoyancy alone-in particular the momentum flux at the source is zero. In practical situations, real sources of heating and cooling in a ventilation system often include initial fluxes of both buoyancy and momentum, e.g. where a heating system vents warm air into a space. This paper describes laboratory experiments to determine the dependence of the 'front' formation and stratification on the source momentum and buoyancy fluxes of a single source, and on the location and relative strengths of two sources from which momentum and buoyancy fluxes were supplied separately. For a single source with a non-zero input of momentum, the rate of descent of the front is more rapid than for the case of zero source momentum flux and increases with increasing momentum input. Increasing the source momentum flux effectively increases the height of the enclosure, and leads to enhanced overturning motions and finally to complete mixing for highly momentum-driven flows. Stratified flows may be maintained by reducing the aspect ratio of the enclosure. At these low aspect ratios different long-time behaviour is observed depending on the nature of the heat input. A constant heat flux always produces a stratified interior at large times. On the other hand, a constant temperature supply ultimately produces a well-mixed space at the supply temperature. For separate sources of momentum and buoyancy, the developing stratification is shown to be strongly dependent on the separation of the sources and their relative strengths. Even at small separation distances the stratification initially exhibits horizontal inhomogeneity with localised regions of warm fluid (from the buoyancy source) and cool fluid. This inhomogeneity is less pronounced as the strength of one source is increased relative to the other. Regardless of the strengths of the sources, a constant buoyancy flux source dominates after sufficiently large times, although the strength of the momentum source determines whether the enclosure is initially well mixed (strong momentum source) or stably stratified (weak momentum source). © 2001 Elsevier Science Ltd. All rights reserved.

Palm frond biochar production and characterisation

Palm oil has been the world's main source of oil and fats since 2004, producing over 45 million tonnes in 2009. Malaysia alone has over 4·5 million hectares planted with oil palm and, based on common practice, ~300 palm fronds are pruned per hectare per year. This agricultural waste is currently either being used as roughage feed or, more frequently, being left between rows of palm trees to prevent soil erosion, or for nutrient recycling purposes. This paper proposes an alternative use for palm frond as a source of biochar. A traditional method commonly use by gardeners in Malaysia to improve soil fertility was used to produce the biochar. A shallow earth pit was dug in the ground for the carbonisation process. The process is described and the impact of carbonisation on the earth wall is analysed and presented. The process was later re-assessed by using TGA-FTIR. Most of the hemicelluloses had fully disintegrated, but the depolymerisation of the cellulose was still incomplete at the carbonisation temperature. Most of the lignin aromatic structure was still present in the biochar. The carbonisation process was repeated in the laboratory and biochar was characterised by using BET, SEM and FTIR. An adsorption isotherm study was conducted and the experimental data were fitted to the Langmuir model. The model predicted Pb2+ adsorption rates of 83·3 mg/g, Cu2+ 41·4 mg/g, Ni2+ 13·0 mg/g and Zn2+ 19·7 mg/g. Copyright © The Royal Society of Edinburgh 2012.