Elsevier

Energy and Buildings

Volume 121, 1 June 2016, Pages 92-103
Energy and Buildings

Validation of a lumped RC model for thermal simulation of a double skin natural and mechanical ventilated test cell

https://doi.org/10.1016/j.enbuild.2016.03.054Get rights and content

Abstract

Most current building thermal codes impose upper limits to the predicted annual building energy demand for heating, ventilation and air conditioning. In the building design phase these predictions are obtained using thermal simulations with variable complexity. The simplest approach uses a single lumped thermal capacitance to model the high thermal mass building elements, combined with five thermal resistances (known as the 5R1C model proposed in EN ISO 13790 standard). This model is used by many European countries as the reference simplified methodology to assess overheating risk and calculate yearly building energy demand. This paper presents a successful extension of this model that allows for its application to the prediction of the internal air temperature of free-running buildings with double skin façades. The extension consists in an increased number of thermal resistances used to model the double skin façade zone. The extended model is validated using a set of detailed thermal measurements obtained in a free-running double skin test cell. For the case analysed the simplifications used in the RC model do not reduce the overall accuracy: the mean absolute error for room air temperature is approximately 1 °C, the same order of magnitude of more detailed EnergyPlus simulations (1.2 °C).

Introduction

The last two decades have seen an increased public awareness of the environmental and operational costs of building energy consumption. As a result, most current building thermal codes impose upper limits to the predicted annual energy demand for heating, ventilation and air conditioning systems (HVAC). In the building design phase these predictions are obtained using thermal simulation models with variable levels of detail and approximations. The most complex buildings require models with several thermal zones and, in some cases, tri-dimensional computational fluid dynamics simulations (CFD). For small buildings with simple indoor climate control systems, such as single-family homes and apartments, there is an increase in use of single thermal zone simulation models. In response to this increased use, the research community is continuously working on improved models with all levels of complexity.

The main challenges of building thermal simulation are in modelling of room airflow and heat transfer in building elements with high thermal mass (floor, walls, etc.). In the majority of simple thermal models room air is considered perfectly mixed and represented as a single thermal node that is connected to all room internal surfaces via thermal resistances [1]. Heat transfer in building elements with high thermal mass can be modelled using variable levels of detail, ranging from detailed finite difference methods to conduction transfer functions (CTF) [2] or, in the simplest approach thermal resistances connected to a thermal capacitance. The thermal resistance and capacitance used in the models have physical meaning and can be calculated approximately from the thermal properties of the building elements. Fig. 1 shows the combination of a perfectly mixed room air approach with a resistance capacitance model for the high thermal mass elements, this approach is known as an RC model. These models are typically labeled according to the number of resistances and thermal capacities used. EN ISO 13790 standard [3], adopted by many European countries as the reference methodology to assess overheating in buildings during summer and calculate building energy demand, uses a single capacitance that represents all high thermal mass elements and five thermal resistances (making it a 5R1C model).

The widespread use of the 5R1C model makes it a preferred target for further development. This paper presents an extension of this model for application to double skin façade (DSF) buildings. The proposed model is validated using a set of detailed thermal measurements obtained in a free-running double skin test cell [4]. The next section presents a review of existing RC models and model validation studies. The following section presents the model methodology followed by the presentation of free-running test cell measurements. The final section presents the results of the model validation.

Section snippets

Literature review

One category of RC models – the lumped parameter construction element models – approximates each building element with high thermal mass (e.g. walls, roofs, floors) by one or more capacitances (Fig. 2). A first-order representation refers to a single capacitance by each building element, a second-order refers to two capacitances. Therefore, the whole room model can have several capacitances and their number increases for higher order representations. It is noteworthy that thermal capacitances

The single thermal zone model

The RC model for a single thermal zone described in EN ISO 13790 standard has five resistances and one capacitance (5R1C). The resistances (or the inverse, the conductances) connect the following temperature nodes: supply air, external air, internal air, ‘star’ and mass. The ‘star’ temperature node results from the delta-star simplification, taking values in the range of internal air, light and heavy elements surface temperatures, although it does not have a specific physical meaning. The

Room and sunspaces thermal envelope

The heavy elements of DSF test cell, with a total area of 58.04 m2, are the external walls, the external roof and the ground floor (Table 2, Table 3). The light elements are the two double glazed windows and vents, that in this experiment are kept closed (Tables 2 and 3). The internal wall partition is also assumed as a light thermal mass element (negligible mass) due to its small area compared to the total heavy elements area. The ratio between total surface area and floor area (At/Af) is

Error parameters

Assuming yi as the measured air temperature at time-step i and yˆi, the corresponding model calculated value, for a total of n time-steps, the following error parameters expressed in temperature units or percentage, by normalisation, considering as normalised factor the measured absolute temperature difference defined by ymax  ymin, with:

  • Mean absolute error (MAE) [°C] and normalised mean absolute error (NMAE) [%]MAE=i=1n|yˆiyi|nNMAE=i=1n|yˆiyi|n(ymaxymin)

  • Mean bias (MB) [°C]MB=i=1n(yˆiyi)n

Conclusions

This paper presented a successful extension of the 5R1C model used in EN ISO 13790 standard that allows for its application to the prediction of the internal air temperature of free-running buildings with double skin facades. The extension consists in an increased number of thermal resistances used to model the double skin façade zone. The predicted room and façade sunspace air temperatures were compared with Energy Plus simulations and detailed measurements performed in a test cell located in

Acknowledgments

The authors acknowledge the support from Fundação Calouste Gulbenkian, through Ph.D. Grant no. 126724.

References (42)

  • H. Buchberg

    Electric analogue prediction of the thermal behavior of an inhabitable enclosure

    Trans. Am. Soc. Air-Cond. Eng.

    (1955)
  • A. Foucquier et al.

    State of the art in building modelling and energy performance prediction: a review

    Renew. Sustain. Energy Rev.

    (2013)
  • European Committee for Standardisation and International Organization for Standardization, EN ISO 13790

    Energy Performance of Buildings, calculation of energy use for space heating and cooling

    (2008)
  • N. Mateus et al.

    Validation of EnergyPlus thermal simulation of a double skin naturally and mechanically ventilated test cell

    Energy Build.

    (2014)
  • G.G.J. Achterbosch et al.

    The development of a convenient thermal dynamic building model

    Energy Build.

    (1985)
  • G. Hudson et al.

    A simple modeling procedure for matlab/simulink

  • C. Underwood

    An improved lumped parameter method for building thermal modelling

    Energy Build.

    (2014)
  • G. Lefebvre

    Modal-based simulation of the thermal behavior of a building: the m2m software

    Energy Build.

    (1997)
  • T.A. Mara et al.

    Empirical validation of the thermal model of a passive solar cell test

    Energy Build.

    (2001)
  • G. Fraisse et al.

    Development of a simpli(r)ed and accurate building model based on electrical analogy

    Energy Build.

    (2002)
  • T. Dewson et al.

    Least squares parameter estimation of a reducer order thermal model of an experimental building

    Build. Environ.

    (1993)
  • M. Gouda et al.

    Building thermal model reduction using nonlinear constrained optimization

    Build. Environ.

    (2002)
  • S. Wang et al.

    Parameter estimation of internal thermal mass of building dynamic models using genetic algorithm

    Energy Convers. Manag.

    (2006)
  • X. Xu et al.

    Optimal simplified thermal models of building envelope based on frequency domain regression using genetic algorithm

    Energy Build.

    (2007)
  • J. Crabb et al.

    A simplified thermal response model

    Build. Serv. Eng. Res. Technol.

    (1987)
  • A. Tindale

    Third-order lumped-parameter simulation method

    Build. Serv. Eng. Res. Technol.

    (1993)
  • J. Schultz et al.

    Winsim: a simple simulation program for evaluating the influence of windows on heating demand and risk of overheating

    Sol. Energy

    (1998)
  • T. Nielsen

    Simple tool to evaluate energy demand and indoor environment in the early stages of building design

    Sol. Energy

    (2005)
  • J.H. Kampf et al.

    A simplified thermal model to support analysis of urban resource flows

    Energy Build.

    (2007)
  • J.-R. Millet

    The simple hourly method of prEN 13790: a dynamic method for the future

  • G. Kokogiannakis et al.

    Comparison of the simplified methods of the ISO 13790 standard and detailed modelling programs in a regulatory context

    J. Build. Perform. Simul.

    (2008)
  • Cited by (0)

    View full text