Air Tightness
Air Tightness
Air Tightness
Air-tightness
Page....................................7-117
1. Introduction
Uncontrolled ventilation from air leakage has an impact on the indoor environment and the building
energy efficiency. In that context, the air tightness is generally an important issue calling for
evaluation. However, that evaluation is not an easy task (Sherman 2004, Sikander 2008, Gränne 2001,
Bankvall 2013). Evaluation of building airtightness, based on field measurements for new buildings,
buildings under construction, as well as for inhabited buildings is often conducted by trained experts.
For their use there are several available methods. Based on the conditions in question and the demands
on the evaluation, each method demonstrates its strengths and drawbacks, described in the literature
(Allard 2012, Nirvan 2010, Lee 2011, Sherman 1992, Hassan 2013). The reported evaluations were
often based on experimental conditions. It can be assumed that the experts’ experience from field
measurements can differ from what is documented under controlled laboratory conditions. Field
experiences can supplement the experimental findings in general but also be useful for future
development of field measurements.
In this paper we have investigated how a number of Finnish and Swedish experts value their expertise,
the performance and usefulness, and expected future use of a selected number of identified methods
for measuring airtightness. The scope of the paper is to problematize usefulness of the methods, how
and when they are used and their potential for future development. The intension was for the findings
of this pre-study to be used for a future similar, but more qualitative, investigation.
3. Questioner
The questionnaire used for the survey is developed by a Finnish marketing company (ZEF 2013) and
utilizes a graphical 2D-interface.The results are reported in Fig 1-3, below, where the smaller circle
with a number (1..n) is referring to the evaluated questions. The position of the point is based on the
average of all who have answered the question. Around each circle there is a darker area/ellipse,
calculated with a standard deviation formula. It shows, in a graphical way, how much answers from
the respondents (experts) vary from each other. The smaller the ellipse is, the more persons agree on
that topic. The model of the ellipse also tells you in which direction/criteria (horizontal/vertical)
agreement is smaller or bigger.
The questions were categorized into three levels; Background, Recommendations and Forecast, all
further introduced in this section.
3.1 Background: Expert awareness questions
The background question is supposed to give information about expertise and practical experience of
the methods among the participants. This shows which methods are most commonly used within the
expert group. The background question was stated as: “I assess my expertise and practical experience
to be good”. Therefore, both information on expertise (which may be only theoretical) and practical
experience will be gathered. This information may be used for qualitative assessments of the answers
in combination with the other questions.
3.2 Recommendations: Method recommendation questions
The recommendations are supposed to indicate which methods that is appropriate to use under certain
circumstances. The circumstances, or scenarios, depend on building type and life stage. The
recommendation question was stated as: “I would recommend this method by the terms of ease of use
and performance for this building type and life stage”. The question was modified to apply to the
different building types and life stages.
3.3 Forecast: Future method use questions
The questions were asked about the potential and future use of the investigated methods. The forecast
question was stated as: “I believe there is likelihood of use and potential of performance for this
method in the future”. The term ‘future’ is used since the exact date is not important.
FIG 1a. Expertise and practical experience of FIG 1b. Performance and ease of use for
fan pressurization with a blower door, fan qualitative methods: Recommendations for
pressurization using ventilation and the tracer single-zone, multi-zone, high rise, leaky as well
gas methods as air tight and new buildings
FIG 3a. Expertise and experience FIG 3b. Performance and ease of FIG 3c. Potential and likelihood of
of experimental methods use for experimental methods use for experimental methods
4.4 General results
The identification of experts for survey in this pre-study was based on mapping procedure supported
by the national research institutes. It was found difficult to identify Swedish and Finnish experts of the
investigated field measurement methods. For the study, 20 experts were identified and the
1
Building Physics Section, Department of Civil Engineering, KU Leuven, Leuven, Belgium
SUMMARY
A full scale test room was constructed for the analysis of the energy performance of a room with night
ventilation with air supply through a hopper window, supported by an extraction fan. Contrary to
typical climate rooms, the room was exposed to realistic climatic conditions. This paper illustrates by
means of measurements the air flow phenomena during night ventilation. The room includes a
concrete floor and ceiling and detailed monitoring is done of local air and mass temperatures, air
flow rate and local air velocities and surface heat flux. Based on a selection of summer nights, an
assessment is made of the flow pattern occurring during night ventilation. Furthermore, to allow
taking into account such flow changes in BES models, the timing of the flow changes are connected to
the Richardson number. Three air flow typologies were identified: forced flow with short-circuit
between supply opening and outlet (Coandă at ceiling), transition flow with limited Coandă effect, and
buoyancy dominated flow with supply air falling down upon entry. The implications for assumptions
used in BES models are discussed shortly, as the flow pattern will strongly influence the surface heat
transfer.
1. Introduction
Night ventilation is a promising technique, able to reduce the cooling load through use of cold outdoor
air at high air change rates to cool down the internal thermal mass of a building in moderate climates
with sufficiently cold nights. A previous study focussed on the determination of the transient surface
convection during night ventilation through numerical simulations using the commercial CFD package
Fluent (Leenknegt et al. 2013). A simplified 2D geometry was simulated with unsteady RANS, using
constant boundary conditions. Various supply air temperatures and velocities were studied, for a flow
time of up to eight hours. Results from this numerical study showed the prevalence of unstable flow
behaviour during constant high ventilation rates. Three flow regimes were observed, illustrated in FIG
1. Buoyant flow was defined as air falling down upon entry, with the main flow path over the floor
surface. At a critical point, the air detaches from the wall under the inlet, moves upwards and attaches
to the ceiling. The transition flow regime is then defined as a flow adhering to the ceiling with variable
attachment length. During the simulations, this flow regime was only an intermediate regime, typically
lasting only up to half an hour. Finally, the room air flow stabilizes into forced flow, with maximum
ceiling attachment length and supply air short-circuiting via the ceiling towards the outlet. These flow
regimes strongly influence the surface convection on the different surfaces during the intensive
ventilation period. This influence was particularly noticeable at the ceiling where the evolution from
buoyant to forced flow increased the surface averaged convective heat transfer coefficient (CHTC) by
a factor 10. This illustrates the importance of taking into account this unsteady flow behaviour in
rooms with night ventilation. However, this is currently not taken into account in Building Energy
Simulation (BES) models.
2.1 Location
The main orientations of the building are north-east and south-west, with the room located on the
north-western side, façade facing south-west. The most prevalent wind direction in Belgium is south-
west, so the building was oriented to have one long façade exposed to normal wind forces. However,
as the room is located at the corner of the building, the wind climate at the window will be more
unstable. The wind climate during the measurements was summarized in FIG 2b, showing the wind
velocity as a function of the wind direction during the measurements. These were measured using an
ultrasonic anemometer, placed on a mast (M) at a height of 10 m at approximately 20 m south-west of
the building. Remaining outdoor climatic conditions were measured in the weather station (W),
(a) Orientation of Vliet building and room (b) Wind velocity as a function of the wind direction
(black) (M = wind mast, W = weather
station)
FIG 2. Overview of Vliet test building orientation and environmental parameters
FIG 3. Vertical section over the length of the test room with sensor locations (round black:
thermocouples, black triangle: tracer gas injection (Kin) and sampling (K), black diamonds:
anemometers (V)
3. Results
3.1 Classification according to flow type
The same three flow typologies that were defined in the numerical study, namely buoyant flow (BF),
transition flow (TF) and forced flow (FF), can be used as well to classify the flow pattern observed
during the measurement campaign. However, the following differences with regard to timing and
occurrence were noted. Firstly, the forced flow pattern was only observed at near isothermal cases or
when the supply temperature is warmer than the room air temperature, placing the room in heating
modus. Secondly, the transition flow, defined by a ceiling attachment length xcl larger than zero but
smaller than room depth D, will be shown to be a near stationary condition, rather than a short
(a) Air temperatures for different flow types (b) Schematic visualisation of flow types
FIG 4. Illustration of different flow types
TABLE 1. Criteria to sort measurement data per night according to flow development
Flow type Criteria Number of nights
forced (FF) Te > TAFb4 0
transition (TF) Te < TAFb4 & Te ≈ TAFb4 18
buoyant (BF) Te < TAFb4 & Te << TAFb4 20
FF → TF Te > TAFb4 → Te ≈ TAFb4 5
FF (→ TF) → BF Te > TAFb4 →Te << TAFb4 3
TF → BF Te ≈ TAFb4 → Te << TAFb4 14
, −
= (1)
⁄ ℎ
Where g gravitational acceleration (m/s²)
β thermal expansion coefficient of air (K-1)
H room height (m)
Tr,c central room air temperature (°C)
Te outdoor air temperature (°C)
volumetric air flow rate (m³/h)
b, h opening width and height of supply window in façade (m)
FIG 7. Range of Richardson values for which flow pattern changes occur (dataset: 39 nights)
5. Acknowledgements
This research was funded by the Research Foundation Flanders (FWO Vlaanderen). Their financial
contribution is greatly appreciated. Furthermore, Patricia Elsen was a great help in conducting the
measurements.
References
Artmann, N. et al., 2010. Experimental investigation of heat transfer during night-time ventilation.
Energy and Buildings, 42(3), pp.366–374.
Beausoleil-Morrison, I., 2000. The adaptive coupling of heat and air flow modelling within dynamic
whole-building simulation. University of Strathclyde, Glasgow, UK.
Goethals, K. et al., 2012. Experimental investigation of the impact of room/system design on mixed
convection heat transfer. Energy and Buildings, 49, pp.542–551.
Leenknegt, S. et al., 2013. Numerical study of convection during night cooling and the implications
for convection modeling in Building Energy Simulation models. Energy and Buildings, 64,
pp.41–52.
Leenknegt, S. & Saelens, D., 2014. Determination of ventilation rate in a full scale night ventilation
experiment through tracer gas measurements. In submitted to Nordic 2014.
Leenknegt, S. & Saelens, D., 2013. Use of dimensionless numbers to classify room air flow and
convection regimes: part 1 - literature review. Building and Environment, submitted,.
1
Building Physics Section, Department of Civil Engineering, KU Leuven, Leuven, Belgium
SUMMARY
A full scale test room was constructed for the analysis of the energy performance of a room with night
ventilation with air supply through a hopper window, supported by an extraction fan.. Contrary to
typical climate rooms, the room was exposed to realistic climatic conditions. The determination of the
air change rate is paramount in a night ventilation system. This paper compares two methods for the
determination of the air change rate: (1) the pressure difference over the orifice in the exhaust
channel and (2) tracer gas measurements using the averaging method on the tracer gas mass balance.
The latter method allowed accounting for the back flow through the hopper window, which was not
captured via the orifice method. A large underestimation was seen in case of the pressure based ACH,
as this method does not take into account two-way flow via the large supply opening. Therefore, a
correlation was developed to correct the pressure based flow rate, to be able to determine the flow
rate in future measurement campaigns, which cannot be combined with tracer gas experiments.
1. Introduction
Night ventilation is a promising technique, able to reduce the cooling load through use of cold outdoor
air at high air change rates to cool down the internal thermal mass of a building in moderate climates
with sufficiently cold nights. To investigate the energy performance of night ventilation, a full scale
test room was constructed inside the Vliet test building of the Building Physics Section of KU Leuven,
Leuven, Belgium. The room differs from typical full scale test rooms in similar research (Artmann et
al. 2010; Goethals et al. 2012), as it is exposed to real climatic conditions through a south-west
oriented hopper window. High air change rates are established through an extractor fan. As ventilation
is done with outdoor air, the supply air temperature is not controlled and the thermal response of the
room in case of realistic boundary conditions can be investigated. This limitation in ventilation control
is compensated by more than 200 sensors, allowing a detailed monitoring of the air, surface and mass
temperatures, as well as local air velocities, surface heat flux and pressure difference. One of the key
parameters in the analysis of the energy performance is the ventilation rate. Therefore, it is crucial to
have a reasonable estimate of the air change rate during the night ventilation experiments. This is
complicated by the influence of wind on the large supply opening. Therefore, the objective of this
paper is to determine the air change rate in this mechanically ventilated case.
One of the commonly used methods to asses ventilation rate is the use of tracer gas, which is also the
main method discussed here. Alternatives found in literature are hot wire anemometry or pressure
difference measurements either by wind on the façade or over an orifice in a extraction/supply duct
(Van Buggenhout et al. 2009; Shen et al. 2012). Furthermore, an orifice was placed in the exhaust
channel to allow determination of the flow rate through the orifice using the measured pressure
difference. This paper will discuss the determination of the ventilation rate through tracer gas
2.1 Location
The Vliet test building is a 25.2 m long and 7.2 m wide building, designed for the study of wall,
window and roof components under real field conditions. The main orientations of the building are
north-east and south-west, with the room located on the north-western side, façade facing south-west.
The most prevalent wind direction in Belgium is south-west, so the building was oriented to have one
long façade exposed to normal wind forces. However, as the room is located at the corner of the
building, the wind climate at the window will be more unstable. The wind climate during the
measurements was summarized in FIG 1b and c, showing a box plot of respectively the wind direction
and wind velocity. These were measured using an ultrasonic anemometer, placed on a mast (M) at a
height of 10 m at approximately 20 m south-west of the building. The box plot is based on the 10th
and 90th percentile, with conditions occurring only 2 % of the time not considered in the calculation
of the spread. Remaining outdoor climatic conditions were measured in the weather station (W),
located on the roof of the building. FIG 2 shows a vertical section of the test room. The room internal
dimensions are 1.80 m wide, 3.45 m long and 2.40 m high. A window with dimensions 1.25 m wide
and 1.60 m high is placed centrally in the façade, starting at a height of 0.68 m above room floor level.
The top part of the window is a motorized hopper window, with a maximum opening angle of 30°.
Opposite the window, an outlet opening is placed, with dimensions of 1.00 m wide and 0.10 m high
and located at a height of 2.10 m and at 0.20 m distance from the ceiling. In front of the outlet, a grid
is placed with horizontal lamellae with 2 cm spacing. The outlet opens into a metal plenum box, which
is 40 cm deep. The plenum is connected with an extractor fan via a 1.50 m long flexible channel with
125 mm diameter. A maximum air change rate of up to 18 h-1 could be realized. An exhaust opening is
made in the NW-façade of the building, provided with wind shield. A flexible duct leads 3.20 m
upwards, connecting fan to exhaust.
(a) Orientation of Vliet building (b) Spread on wind direction (c) Spread on wind velocity
and room (black) (M = wind mast, during analysed nights at 10 m high during
W = weather station) analysed nights
FIG 1. Overview of Vliet test building orientation and environmental parameters
V = − (1)
Where V zone volume (m³)
C tracer gas concentration in the zone (10-6 mol/mol)
t time (h)
QT injected tracer gas volumetric flow rate (m³/h)
Q volumetric flow rate (m³/h)
Sometimes however, the standard techniques are not usable, for example if the air change rate is
unstable and the controller is too slow to follow the changes. When no assumptions are made
regarding constant injections and concentrations apart from the continuity equation itself, Sherman
describes three different approaches for determination of the ventilation directly from the
concentration data: (1) regression, (2) integral and (3) averaging techniques, of which the latter
technique is used here (Sherman 1990). The continuity equation cannot be used to determine
instantaneous ventilation, but it can be averaged over a measurement period T to get the average
ventilation, as given in equation (2).
= − ln (2)
3. Results
The following results will start with a description of the ventilation rate based on the averaging
method applied to the measured tracer gas concentrations. The influence of wind climate is shortly
discussed, resulting in unstable flow rates. These results can then be compared to the air change rate
based on the pressure difference over the orifice in the exhaust channel. The problems concerning this
latter method are identified and addressed.
FIG 3. ACH based on averaging technique from equation (2) over 30 minute periods as a function of
orifice opening and fan setting (wind box present: black bars, no wind box: grey bars)
For the measurements with fan power 1, the effect of the wind box is clearly visible, with much more
stable flow rates. At fan power 2, all measurements are conducted with the wind box in place. Most of
the calculated ACH display only a very limited range, indicated by the box plot. Nevertheless, the
minima and maxima, indicated by the vertical lines, are still quite far apart. The measurements with
the widest ACH-range correspond again to the moments with the most unstable wind direction.
At fan power 3, measurements were conducted with and without wind box. At this high fan power, the
correlation between wind direction and ACH range is less strong. Some measurements show a very
stable ACH, despite a large variation in wind direction. Nevertheless, the wind box again results in
much more stable air flow rates. One case at fan power 3 and orifice opening 1 is an exception, with
the air change rate ranging from 17 to nearly 23 h-1 in one night, illustrated in more detail in FIG 4a,
displaying the ACH, along with current wind direction and velocity. The tracer gas based air change
rate nTr is calculated for a 30 minute averaging period (nTr30m), which is indicated on the graphs by
horizontal bars. The wind velocities (WindVel) are rather low, i.e. equal to or lower than 2 m/s.
However, the wind direction (WindDir) is continuously nearly perpendicular to the window. This
could explain why the moderate wind still has such a high influence on the air flow rate. However, the
increase in ACH over the presented period seems not correlated with the wind climate and could not
be explained.
Another interesting case is shown in FIG 3 at fan power 1, orifice opening 5 and with wind box, but
nevertheless with a large spread. The main dataset is very stable at 8.5 to 9 h-1, though a peak up to 10
h-1 is seen. This day is shown in more detail in FIG 4b. Cross-referencing with wind climate explains
this variation. During the stable phase, the wind velocity is rather low, i.e. equal to or lower than 2
m/s. However, the local peak in ACH coincides with a sudden increase in wind velocity up to 5.5 m/s.
(a) Tracer gas based ACH on day 62 (F3V1) (b) Tracer gas based ACH on day 48 (F1V5)
FIG 4. Detail of tracer gas based ACH, correlation to climatic data
FIG 5. Comparison of pressure based ACH (nΔP) and tracer gas based ACH (nTr30m), see equation (3)
for Δnrel (filled markers = with wind box, not filled markers = without wind box)
During the analysed periods of intensive ventilation, the pressure difference is always 10 Pa or higher,
resulting in an uncertainty on the calculated ACH of approximately ± 7 %, as the uncertainty on the
measured pressure difference is very small. FIG 4 shows for both days also the pressure based ACH,
indicated by nΔP. The uncertainty is indicated by the vertical error bars, based on the uncertainty on the
pressure difference measurement, combined with the pressure-flow characteristics. Also the outdoor
temperature is indicated by Text. Two conclusions can be drawn: firstly, the pressure based air change
rates nΔP are significantly lower than the tracer gas based air change rates nTr, and secondly, they are
independent of the wind climate. They may fluctuate somewhat by the indoor-outdoor temperature
difference, although no clear correlation was found. Furthermore, the ventilation rate is driven by the
extraction fan in all measurements. The pressure based ACH for the full measurement period are
shown in box plots in FIG 5. The tracer gas based median values nTr30m of each measurement are
shown as well. To illustrate the difference between both methods, the relative difference Δnrel
according to equation (3) has been plotted on the right axis, with nTr30m based on the 30 minutes
(a) Relative difference between tracer gas based (b) Correlation between pressure-based on tracer
and pressure based ACH gas based ACH
FIG 6. Correlation to predict actual ACH based on pressure-based ACH
Future research is aimed at predicting the air flow pattern in the room during intensive night
ventilation, as this strongly influences the local surface convection. This was shown in previous
research through numerical simulations (Leenknegt et al. 2013) . The presence of mixing fans would
disrupt the flow pattern, which is a crucial part of the research. Therefore, it is impossible to perform
tracer gas measurements during the night ventilation experiments and the ACH must be determined
based on the pressure difference over the orifice. It was shown before that this value only takes into
account the flow rate leaving through the extraction duct and ignores the additional air leaving through
the window opening. As such, the pressure based ACH is an underestimation of the actual value and a
correlation was sought and found between nTr and nΔP to obtain a corrected value nΔP,corr. This
correlations is illustrated in FIG 6b and given by equation (4). Despite the good correlation, there is
still some scatter, and the corrected value will have an uncertainty which is estimated at ± 15 %, based
on the typical difference between the measured values and the correlation.
.'(
∆ ,!" = 2.35 ∆ (4)
5. Acknowledgements
This research was funded by the Research Foundation Flanders (FWO Vlaanderen). Their financial
contribution is greatly appreciated. Furthermore, Patricia Elsen was a great help in conducting the
tracer gas measurements.
References
Artmann, N. et al., 2010. Experimental investigation of heat transfer during night-time ventilation.
Energy and Buildings, 42(3), pp.366–374. Available at:
http://linkinghub.elsevier.com/retrieve/pii/S0378778809002412.
Van Buggenhout, S. et al., 2009. Influence of sampling positions on accuracy of tracer gas
measurements in ventilated spaces. Biosystems Engineering, 104(2), pp.216–223.
Goethals, K. et al., 2012. Experimental investigation of the impact of room/system design on mixed
convection heat transfer. Energy and Buildings, 49, pp.542–551.
Leenknegt, S. et al., 2013. Numerical study of convection during night cooling and the implications
for convection modeling in Building Energy Simulation models. Energy and Buildings.
Shen, X., Zhang, G. & Bjerg, B., 2012. Comparison of different methods for estimating ventilation
rates through wind driven ventilated buildings. Energy and Buildings, 54, pp.297–306.
Sherman, M.H., 1990. Tracer-gas techniques for measuring ventilation in a single zone. Building and
Environment, 25(4), pp.365–374.
SUMMARY
Leakages in the building envelope may lead to increase of primary energy demand and carbon foot-
print of buildings. Additionally, particle and bacteria infiltration or exfiltration are possible. Fur-
thermore, the building shell may be damaged by hygric enveloping surface infiltration. For investiga-
tion of the impact of leakages in the building envelope, different measurements and simulations are
used. The disadvantage of the both is that the knowledge on discharge coefficient of leakages in real
buildings is limited. Therefore, the discharge coefficient is estimated or based on sharp edged, circu-
lar standard orifice. But real leakages are different, e.g. the shape is not comparable to a standard
orifice. Due to these, the flow through different leakages in perforated air tightness layers of light-
weight timber frame constructions was analysed. The applied pressure differences were real meas-
ured values of long term investigations on building shells. At these specific pressure differences, the
flow profile through leakages is different to studies of standard orifice. The volumetric flow through
leakages is significantly influenced by this difference. The results show the impact of the different
leakage parameters onto the discharge coefficient. Thus, simulations should not be based on standard
orifice investigations.
Introduction
In walls of new or refurnished buildings, it is common to integrate different insulation. Thus, a
reduction of energy losses through the building envelope is possible. However this climate separation
of interior and environment represents also a potential danger. If the water vapour diffusion resistance
of the room side materials is too small, condensation may arise in the wall construction. For
prevention, room-side diffusion inhibiting or blocking materials e.g. vapour barriers are arranged.
According to (Schulze H. 2011), awareness regarding diffusive moisture transport processes has
increased. Therefore, moisture problems caused by diffusion in constructions of new buildings are
infrequent at present time. A far bigger problem is given by the deficient air tightness of the building
envelope. According to (Biskop R. 2008), these leakages are caused by mistakes in planning and
execution. Especially during the heating season, bacteria, particle and moisture may be transported by
convection from the interior to the environment. This transport of matter may lead to deposits in the
wall structure. This might result in growth of mould and mildew or the formation of heat bridges e.g.
in case of condensation in non-capillary active components. Additionally, the convective heat
exchange leads to increased energy losses. Light-weight timber constructions are particularly at risk
concerning leakages. Damages in the air tight and protection layer are easily possible and the used
insulation materials are frequently air permeable. In addition, the convective leakage flow is promoted
by the non-obligatory air tightness of the weather protection layer. These findings led to the
development of simulation software with the aim of analysing the disadvantages of leakage flows on
buildings. The best-known algorithms are implemented in the software of Delphin and WUFI. These
algorithms are based on leakage flows caused by density differences (Häupl P. et al. 1997, Künzel
H.M. 2010). Only a few calculation models are developed for leakage flow analyses under real
Basics
The pressure difference between interior and the surrounding area is the potential of flows through
leakages in the building envelope. In naturally ventilated buildings, these are influenced by wind and
buoyancy. The change of pressure differences occurs in intervals of seconds (Schmidt J. et al. 2010b).
According to the Bernoulli equation, the leakage flow is from the higher to the lower pressure level.
In order to analyse and simulate the exfiltration flow through a room-side layer with a single leak and
adjacent permeable insulation in a light-weight timber structure, knowledge on the pressure loss
caused by the leak is required. Commonly, the cross-section of a leak in the building shell is much
smaller than the cross-section area of the room, which forces a constriction of the exfiltration flow.
This results in contraction and velocity losses. The contraction factor Ψ is affected by the ratio of
stream and leakage cross-section. In addition, the shape of leakage edge is crucial (fig. 1, 2). Ψ of e.g.
sharp edged standard orifice is 0.61 … 0.64 (Bohl W. et al. 2005). The velocity coefficient Φ of the
leakage flow is influenced by the macroscopic and microscopic geometry of the opening, the viscosity
and the surface tension (if not air) of the fluid. Φ is in case of sharp edged standard orifice 0.97 (Bohl
W. et al. 2005). The product of Ψ and Φ gives the discharge coefficient ζ. “Caused by the difficulty to
investigate Ψ and Φ experimentally separated from each other, values of these coefficients are only
scarcely documented in the literature.” (Bohl W. et al. 2005) Theoretical considerations were already
carried out by Torricelli and Borda in 17th and 18th century (Borda J. C. 1769).
2 ⋅ ∆p
AL ⋅ ζ ⋅ = Vreal (1)
r ⋅ (1 − N ² )
Vreal
⋅ (1 − N ² ) = ζ (2)
V
theo
f ⋅L ρ ⋅u ²
+ ΣC ⋅ = ∆p (3)
dh 2
In eq. (3) ΣC is representing the sum of local losses. f is the Darcy-Weisbach friction factor. In case of
fully developed laminar flow it is equal 96/Re (Langmans J. 2013). The other parameters in eq. (3) are
L the length, dh the hydraulic diameter and the mean square velocity u ² in the channel. Re denotes the
Reynolds number, which defines the stream profile (laminar, Re ≤ 2300). DIN EN ISO 5167 provides
additional opportunities for estimation of ζ under conditions of 0.05 ≤ N ≤ 0.64 and 5 · 104 < Re <
107. A tabulation of some ζ-values of standard orifice is implemented in this standard. Furthermore
the Reader-Harris/Gallagher equation may be used for estimation under knowledge of Re at N ≥ 0.1
(Böswirth L. 2007). It is only applicable to pipes with circular standard orifice. Caused by the
approach of N 0 and the laminarity of flows through common leakages in buildings, (Langmans J.
2013) values that these models are unsuitable for forecasts of pressure losses. Hence, the discharge
coefficient was analysed using eq. (2).
Measurements
The investigations were carried out on vapour barrier foils of about 1 x 1 m in the “Measurement
Setup for research on convective moisture transfer (MCMT)” (Schmidt J. et al. 2010a). Every
investigated leakage was sharp edged by use of special tools. The climate conditions between the over
and under pressure chamber were equal at all measurements. The perforated vapour barrier was the
separation of these two connected chambers. For the early investigations, the foil was stiffened
without influencing the flow. The volumetric flows through the leaks were measured parallel by use
of Tracer-gas-system “TGS” and calorimetric flow meter “MSD”. According to (Böhle M. et al. 2002,
Leick Ph. 2008), the flow-through was under open outflow wherein the air flow after overcoming the
obstacle could open conically. First, the pressure difference between the inflow and outflow-side were
calculated by weather data of a German test reference year (Christoffer J. et al. 2004). Later, these
values were arranged based on long term measurements at real building facades. According to the
design of the MCMT (Schmidt J. et al. 2012), the maximum total leak size was 10 cm ², at a
maximum pressure difference of 350 Pa. Under real measured pressure differences (Schmidt J. 2013)
the leakage flow is usually laminar. Due to poor execution of real buildings, cavities between the
insulation and airtightness layer may occur. In order to analyse their influence on the leakage flow of
perforated vapour barrier foils, double chamber measurements were carried out under free arching of
Results
The first part of the study related to the analysis of the flow of differently shaped leaks. Therefore,
leakages of the same cross-section and shape of edge were investigated in stiffened PE-foil.
FIG 3. Volumetric flow through different shaped FIG 4. Comparison of volumetric flow through
leaks, AL = 0.2 cm² different shaped leaks to the calculated
leakage flow, AL = 0.2 cm²
Figure (3) shows the leakage flow through 0.2 cm² leakages of different shape at different pressure
differences. In this context, the flow may be approximated by a leak according to (DIN EN 13829,
ASTM E779-10, ASTM E1186-03) as a function of the applied pressure difference by an exponential
function. As leakages in figure (3, 4) were used: an edged rectangular leak (0.1 x 2 cm), a squared
leak (0.45 x 0.45 cm) and a circular leak (d = 0.5 cm). The measurement error between the carried out
parallel measurements was about 3%. The association for air tightness in the building trade “Flib”
(2008) suggests using a ζ-value of 0.61 to analyse the equivalent leakage area of a building.
Calculated flows, using this ζ-value, were compared with the measurement results, see figure (3, 4).
According to figure (4), the calculated flows had an average difference of about 20 % compared to the
measured flows. When using the hydraulic diameter to calculate the flow rate, the error would
increase further. For e.g.the rectangular leak, it would yield dh = 0.19 cm. Thus, a much smaller
calculation cross-section would be given.
The discharge coefficient is calculable by use of eq. (2) and neglecting N. According to figure (5), a
linear relationship exists between theoretical and measured volumetric flow. The inclination of the
linear fit defines the correlation between the two flow rates. Thus, the inclination resembles a
resistance coefficient respectively the discharge coefficient. The properties of the ζ-value were
investigated on sharp edged, circular leaks of different cross-section. As shown in figure (6), the
discharge coefficient ζ is hyperbolic, approaching a limit value. This relationship can be
approximated by an exponential function. The result for circular holes was similar (figure 7). Both,
figure 6 and 7 show that the discharge coefficient approaches the theoretical, calculated contraction
coefficient according to Torricelli and Borda (1769) with increasing leakage size. This suggests that
the influence of the velocity on Φ decreases when the leakage area increases.
FIG 7. Change of the discharge coefficient in comparison of the change of the leakage area of
annular, circular, sharp edged leaks
In the next step of this investigation, the deformation behaviour of a 1 m² large, 0.2 mm thick,
perforated PE-foil was examined. In figure (8) is shown that the PE-foil is arching like a nozzle due to
the one-sided pressure load. The maximum of arching was 7.4 ± 0.1 cm at a leak size of 0.2 cm² and
350 Pa of pressure difference (figure 9).
FIG 8. Arching of perforated PE-foil at ∆p = 2Pa, FIG 9. Maximum of arching of perforated PE-
AL = 0.2 cm², circular, sharp edged leak foil at different ∆p, AL = 0.2 cm²,
circular, sharp edged leak
FIG 10. Impact of foil arching onto the leakage FIG 11. Comparison of leakage flows through a
flow compared to calculated flow with single leak with the total cross-section
ζ =0.61, AL = 0.2 cm², circular, sharp of 5 small leaks of the same properties
edged hole in deformable PE-foil
In figure (10) is shown that the smallest variation between the flows through arched and stiffened PE-
foil occurs for circular leaks. The difference of leakage flow between deformed and undeformed foil
increases with increasing pressure difference (figure 10) for leakages of non-circular shape. The ratio
between calculated and measured volumetric flows was never 1 for all measurements (figure 10).
These results show that only one value for ζ is not generally valid for every single leakage.
Further investigations should clarify if the ζ-value of the volumetric flow of one leak is equal to the ζ-
value of several leakages of the same total cross-section and properties. For the experiment, a leakage
area of 1 cm² was chosen. The reference was the flow through a single, sharp edged, centred leakage
in deformable PE-foil. The total flow at different pressure differences for the single leak was
compared to 5 leakages of AL = 0.2 cm² (total 1 cm²) with different positions.
In figure (11), circle means a uniform distribution of the small holes with a 30 cm radius around one
central leak. The definition star describes a hole distribution with one centre leak and 4 leaks in a
distance of 5cm on the ordinate and abscissa. Coincidence “coinc.” denotes a freely chosen, non-
definable distribution of the small leakages. Figure (11) shows, that the total leakage flows are similar
for all 5-leak-tests. The relationship of the flow through the single leak and the 5 holes of the same
total cross-section is linear. However, the inclination of the flow fit between the 5 small to the one
large leak shows a deviation of 11 to 13 %. In order to be able to define only one discharge coefficient
for existing leakages in the building envelope, e.g. in Zone models or to determine the equivalent
leakage rate, accurate knowledge about the existing leaks is essential. Even the size of the leakage can
be critical.
Conclusions
The investigations of leakage flows through PE-foil show that the volumetric flow through a leak in
lightweight timber structures and the adjacent layers is significantly influenced by the discharge
coefficient ζ of the room-side leak (exfiltration). It is shown, that the assumption of a discharge
coefficient of standard orifices to describe a leakage flow under real pressure differences (wind +
buoyancy) may contain errors. The ζ-values of building leakages deviate to ζ-values of standard
References
ASTM E1186-03, 2009, Standard Practices for Air Leakage Site Detection in Building Envelopes and
Air Barrier Systems.
ASTM E779-10, 2010, Standard Test Method for Determining Air Leakage Rate by Fan
Pressurization.
ASHRAE Standards Committee, American National Standards Institute, 2009, ASHRAE Handbook -
Fundamentals. Atlanta: ASHRAE Customer Service.
Biskop R., 2008, Bauphysikalische Probleme - Fallbeispiele aus der Praxis des mo-dernen Bauens, 1.
Workshop InnoProfile / nuBau, CD, Bauhaus-Universität Weimar (editor), Germany.
Bohl W. & Elmendorf W., 2005, Technische Strömungslehre, Stoffeigenschaften von Flüssigkeiten
und Gasen, Hydrostatik, Aerostatik, inkompressible Strömungen, kompressible Strömungen, Strö-
mungsmesstechnik, 13.th ed., Vogel, Würzburg, Germany.
Böhle M., Etling D. et al., 2002, Prandtl - Führer durch die Strömungslehre, Grundlagen und Phäno-
mene, 11th ed., Vieweg, Braunschweig, Germany.
Borda J. C., 1769, Mémoire sur l'Ecoulement des. Fluides par les Orifices des Vases, Mémoires de
l'Académie, Royale des Sciences, Paris, France.
Böswirth L., 2007, Technische Strömungslehre: Lehr- und Übungsbuch, 7th ed., Vieweg, Braun-
schweig, Germany.
Buckingham E., 1931, Beitrag zur Berechnung der Kontraktionszahl, Mitteilung des Bureau of
Standards, Washington, Rundschau, Forschungsband, vol. 2., 5th ed., Fritz W. (editor), Berlin,
Germany, pp. 185–192.
Christoffer J., Deutschländer Th. et al., 2004, Testreferenzjahre von Deutschland für mittlere und
extreme Witterungsverhältnisse TRY, Selbstverl. des Deutschen Wetter-dienstes, Offenbach a.
Main, Germany.
KEYWORDS: Air flow, tracer gas measurement, comfort, façade systems, cooling ceiling, cooling
capacity
SUMMARY:
The aim of the work is to clarify how the façade system and the cooling ceiling impacts the air
movement within a room. Therefore tracer gas measurements with different inject locations and
different room set-ups were done.
One part of the results show the impact of different façade systems – a single skin façade and a single
storey double skin façade – on the air movement. The measurements show an increase of
concentration within the ceiling zone which is attributed to a higher air flow rate over the ceiling.
Furthermore different set ups of the cooling ceilings were measured – one with natural ventilation
over the suspended cooling ceiling and one without natural ventilation. The measurement results
show that there is a decrease of concentration due to missing natural ventilation.
The measurement results show a clear difference in the rate of concentration on different room points
due to solar radiation, the location of impact and the room set up.
1. Introduction
This paper presents the impact of different room boundary conditions (façade system and ventilation
rate over cooling ceiling) on the air movement inside an office room. The motivation for this work is
that the calculation of comfort in the common areas of buildings plays a key role for the acceptance of
new concepts. In the completeness the calculation of complex flow processes in closed and open
façades, the dynamic of building elements' performance in combination with the control system is at
the moment, if possible, very time consuming. The aim is to develop a simplified room model for the
evaluation of building concepts in terms of energy efficiency and comfort.
Measurements for the model validation are one part of the development. A series of measurements
according to the operative temperature and the cooling capacity have been carried out and the results
show that there is an impact on the cooling capacity with respect to different façade systems. There is
an increase of the cooling capacity of about 20 W/m²cooling area in dependency of the temperature
difference between the operative temperature in the middle of the room and the mean medium
temperature due to a second internal glazing (Eder, 2012).
For verification of the assumption, that the capacity increase is due to a higher air flow rate over the
ceiling and therefore to a higher convection load, a mathematical model is built. The validation of this
model is done by a series of tracer gas measurements and the results of these measurements are
presented in this paper. Chen (2010) describes an overview of ventilation assessment methods. The
methods are divided into analytical, empirical, small-scale experimental, full-scale experimental,
FIG 1. Floor plan (left side) / section (right side); location of tracer gas injection is highlighted; floor
plan: direct injection to the room (1); concrete ceiling (2) cooling ceiling (3); supply air area (4);
façade area (5)
The U-Values of the surrounding areas of the room are summarized in the following list:
glazing: U = 1.10 W/m²K
frame (percentage 8%): U = 2.95 W/m²K
ceiling / floor: U = 1.47 W/m²K
Interior wall (gypsum cardboard): U = 0.65 W/m²K
Interior door: U = 1.80 W/m²K
There is the possibility to change the façade system from a single skin façade to a single storey double
skin façade by adding a second internal glazing element. These two façade systems were tested, both
with an internal shading element. The difference between these two façade variations are shown in
figure 3, there the difference is highlighted between picture 1 and 3, and between 2 and 4 respectively.
2.1 Cooling
The cooling of the rooms is done by a suspended cooling ceiling with capillary pipe mats, which are
placed in the ceiling panels. The cooling area is splitted into a façade area and a room area, these are
overlapping so that there is a higher capacity in the façade area. For the room area of the measurement
room the active area of the cooling ceiling is 5.56 m² and for the façade area it is 3.19 m². Overall
2.2 Ventilation
There is a mechanical ventilation system. The design of the ventilation system is as full air
conditioning – heating, cooling / dehumidification and humidification.
The supply air for the room is realized via the raised floor, the whole floor is air leading. The air
outlet is over “slot diffusors” next to the façade. There is a flow rate of about 50 m³/h for the room,
which means an air change rate of about 2.0 1/h with a supply air temperature of about 21 °C.
The exhaust air is realized by exhaust air lights (light fixture combined with extract air terminals).
There are flexible ducts between the lights and the suspended ceiling in the corridor. The idea of
doing the return air in this way is to remove the heat were it is generated. The suspended ceiling in the
corridor is air leading as well.
3. Measurement Set up
The goal of the measurement is to paint a picture of the air flow movement in the room and to see if it
is changing with different façade situations. Therefore a series of tracer gas measurements was done.
The ASHREA Handbook of fundamentals (2009) describes three different tracer gas methods to
determine the air change rate:
concentration decay or growth method
constant concentration method
constant injection method
The concentration decay method, measuring the decay of tracer gas is the simplest and most often
used method. The following three conditions should strictly be satisfied using single point
measurements. For single point measurements the tracer gas should be uniformly mixed in the
enclosure, there should be no unknown tracer gas sources and the ventilation flow should be the
dominant means of the tracer gas removal from space to get reliable results (according ASHREA
Handbook of fundamentals (2009)).
According to Charlesworth (1988) it is very difficult to achieve a uniformly mixture of the tracer gas
in the entire enclosure of in-situ measurements. To solve this problem tracer gas concentrations
should be measured at several points of the considered areas and it is assumed that the mean value is
the average concentration for the entire enclosure.
As tracer gas carbon dioxid (CO2) was used, as described in Laussmann (2011) it fulfils a number of
the specification of a good tracer gas.
FIG 2. Measurement set up – floor plan (left side) and section (right side) including the sensor
numbers, where the first 4 numbers are the maximum concentration rate and the last 2 are the
numbers referring to the evaluation diagrams
FIG 3. Schematic overview of the four different room types; from left to right: single skin
façade+closed suspended cooling ceiling; single skin façade+open suspended cooling ceiling, single
2. Results
2.1 Decay curves – CO2 concentration
As a result of the measurement there are curves of concentration for each variant and each sensor. The
measurement period for one variant is 1 hour, but the main air movement is within the first 10
minutes, then the concentration is equal within the room (see figure 4 and 5).
Figure 4 represents the measurement results for room type 4 (single storey double skin façade+open
ceiling) for the different locations in the room. The results are for a clear sunny day and the CO2
impulse was injected directly into the room. In the diagram there are results for different room
heights, where the black line represents the measurement points at 1m height, the black dotted line the
measurement results at 2m room height and the grey line the measurement results for the ceiling. The
location of the measurement is marked in the diagram at the right corner. The figure shows the change
of concentration for different locations within the room and the measurement results indicate that
there is no big difference of concentration with respect to the room location, the existing difference
can be explained by the supply air location next to the façade.
FIG 4. Measurement result of CO2 concentration a different location in the room with the impact in
the room for room type 4, the numbers shown in the legend are referring to the numbers of the sensor
(last two numbers in fig 2). The measurement results of the air temperatures (t) for the different room
zones are: troom= 24.5 °C, tfacade = 31.2 °C, tceiling = 23.5 °C and tfloor= 21.8 °C.
FIG 5. Measurement result of CO2 concentration with respect to different impact location for room
type 4 (dskin+openc); with solar radiation; measurement period 1 hour for each variant; different
concentrations in the specified areas within the first 10 minutes, after the first 10 minutes the
concentration is equal in each area independent of the impact location
FIG 6. Measurement results of CO2 concentration for the different room types with solar radiation,
the tracer gas injection directly into the room (floor area). The diagram shows an increase of
concentration in the ceiling area due to the single stores double skin façade and a decrease due to the
lack of natural ventilation for the cases with the closed suspended ceiling.
FIG 7. Measurement results of CO2 concentration for the different room types without solar
radiation, the tracer gas injection directly into the room (floor area). The diagram shows that there is
almost no difference of concentration between the room types
4. Acknowledgements
The doctoral program “Energy Systems 2030 (ENSYS 2030)” is an internal research project funded
by the Vienna University of Technology, which has the main goal to provide doctoral students within
the research area "Energy and Environment," a remarkably structured and interdisciplinary education.
10 PhD positions are advertised and supported as part of the College of the Vienna University of
Technology.
Energy and Environment is part of the development plan of the Vienna University of Technology
2010+ as one of the five research priorities. The doctoral program ENSYS 2030 has the following key
themes: solar energy, energy storage, electro mobility and energy active building.
References
ASHRAE. Handbook of fundamentals.2009. Atlanta, USA: American Society of Heating,
Refrigerating and Air-Conditioning Engineers, Inc.
Charlesworth P.S.1988. Air exchange rate and airtightness measurement techniques - an applications
guide. Coventry, UK: Air Infiltration and Ventilation Centre
Chen, Q., Lee, K., Mazumdar, S., Poussou, S., Wang, L., Wang, M., and Zhang, Z. 2010. “Ventilation
performance prediction for buildings: Model Assessment,” Building and Environment, 45(2), 295-
303.
Eder K., Steininger C., Bednar T. 2012 “Enhancing the performance of a cooling ceiling by an
innovative façade system”; Proceedings of the 5th International Building Physics Conference, IBPC,
2012, Kyoto
Laussmann D., Helm D. 2011. Air Change Measurements Using Tracer Gases: Methods and Results.
Significance of air change for indoor air quality. Chemistry, Emission Control, Radioactive Pollution
and Indoor Air Quality, Dr. Nicolas Mazzeo (Ed.), ISBN: 978-953-307-316-3, pp. 365-399
KEYWORDS: façade systems, natural ventilated double skin façade, computational fluid dynamics,
operative temperature, cooling capacity, comparison measurement and CFD calculation
SUMMARY:
One of the key issues for the calculation of a room’s condition and the resulting thermal comfort is
the knowledge of the air flow distribution in the room. A simplified but still accurate enough CFD
model should be developed to determine the effect of different façade types on the cooling capacity of
a cooling ceiling. A parametric study using the program COMSOL Multiphysics 4.3a was done. The
output of the calculations was compared with the results of a series of in-situ measurements and is
presented in this paper.
The measurement and the calculation results show an impact on the cooling capacity by changing air
flow conditions within the room. The design of the façade system can increase the capacity of a
cooling ceiling and therefore decrease the operative temperature in the room. Neglecting the flow
characteristics in the room can lead to incorrect results, especially for complex façade systems in
combination with component activation (e.g. cooling ceiling).
Therefore an enhanced building model for the design of future buildings, taking into account the
interaction of building envelope and building service systems should be developed using the results of
the CFD calculation. Recommendations for energy active buildings related to room conditions, room
comfort and the energy demand can be formulated by simulations over yearly periods.
1 Introduction
This study presents results of the cooling capacity of a cooling ceiling in dependence of two different
façade systems – a single skin façade and a single storey double skin façade. The influence on the
cooling capacity due to the façade systems was measured by an intensive set of measurement and the
results are presented in Eder et al. (2012) and are used for the validation of the mathematical model.
Due to multiple coupled physical phenomena (air flow, convection, conduction and radiation) the
behaviour of double skin façades is very complex. Both, experimental and numerical models for the
optimization of the performance of double skin façades have been developed e.g. network models
(Tanimoto et al., 1997), zonal models (Jiru et al., 2008), energy simulation with coupled air flow
models (Stec et al. 2005) and detailed computational fluid dynamics studies (Manz, 2004).
Most of the developed models are for external double skin facades. This paper studies the effect of a
(internal) single storey double skin façade on the HVAC system. The diving force is natural buoyancy,
therefore the air flow is not easy to control and the dependence on weather conditions (mainly solar
radiation) makes it discontinuous. The validated results of the air flow calculation should be
2 Theoretical background
The software program COMSOL Multiphysics version 4.3a was used, with the aim to calculate the
impact of air movement with different façade systems. As physical model non isothermal flow was
used to describe the problem. The main equations are the Navier Stokes (momentum balance and mass
balance), heat balance and energy transport (convection and conduction).
One of the main topics is the modeling of the façade. Because the simulation program COMSOL
Multiphysics version 4.3a consider all radiation as the same, it is not possible to make a distinction
between large infrared radiation and the solar spectrum (UV, light & short IR). Within this study the
interaction between the façade system and the cooling ceiling is the main topic, therefore the surface
to surface radiation has to be considered. For this work the incident solar radiation was realized by
heat impacts to the individual layers of the façade according to its absorption and solar transmission
coefficient.
3 2-D Model
3.1 Geometry and mesh properties
The geometry of the simulation model is according to the architectural plans, a section through the
room. A simplification of the model, from three dimensions to two was done. The symmetry is along
the depth of the room (see figure 1).
The main geometric data of the room is summarized in table 1.
FIG 1. Simplification of the model geometry from three dimensions to two dimensions
shading shading
2nd glazing
FIG 2. Façade type 1 (left side): single skin façade with internal shading element / façade type 2
(right side): single storey double skin façade with a shading element between the glazing elements.
4 Results
4.1 Operative Temperature
In figure 3 the results of the operative temperature in the middle of the room are shown for the two
façade systems. The results for the single skin façade show an increase of the operative temperature
from 24.7 to 28.5 °C with an increasing solar impact. The mean capacity of the cooling ceiling variate
between 40 and 80 W/m²cooling area. The effect of the single storey double skin façade is a smaller
increase of the operative temperature in the middle of the room due to a higher efficiency of the
cooling ceiling. The diagram shows, that the operative temperature stays around 26 °C and the cooling
capacity increases with respect to the solar impact.
FIG 3. Comparison of operative temperature and the cooling capacity for façade type 1 (left side) and
façade type 2 (right side); the temperature is below 27°C with façade type 2, the cooling capacity is
increasing with increasing solar radiation; with façade type 1 the operative temperature is getting
above 27°C and the cooling capacity is lower.
FIG 5. Cooling capacity in dependence of the depth of the room with an incident solar radiation of
800W/m² and a temperature difference of 9 K between the mean medium temperature and the
operative temperature in the middle of the room.
Figure 6 shows the resulting mean cooling capacity with different impacts due to the incident solar
radiation for the two façade systems. The temperature difference between the mean medium
temperature and the operative temperature in the middle of the room is 9 K.
FIG 7. Measurement results on operative temperature (left side) and cooling capacity (right side;
(M)..measurement; (C)..CFD calculation) – there is a good accordance between the measurement
results and the CFD calculation; the difference of operative temperature is about 0.2 K;
6 Acknowledgements
The doctoral program "Energy Systems 2030 (ENSYS 2030)" is an internal research project funded
by the Vienna University of Technology, which has the main goal to provide doctoral students within
the research area "Energy and Environment," a remarkably structured and interdisciplinary education.
10 PhD positions are advertised and supported as part of the College of the Vienna University of
Technology
Energy and Environment is part of the development plan of the Vienna University of Technology
2010+ as one of the five research priorities. The doctoral program ENSYS2030 has the following key
themes: solar energy, energy storage, electro mobility and energy active building.
References
Eder K., Steininger C., Bednar T. 2012 “Enhancing the performance of a cooling ceiling by an
innovative façade system”; Proceedings of the 5th International Building Physics Conference, IBPC,
2012, Kyoto
EN 14240. 2004. Ventilation for buildings - Chilled ceilings - Testing and rating
Jiru T.E., Haghighat F., 2008 “Modeling ventilated double skin façades – a zonal approach”, Energy
and Buildings 40, 1567–1576.
Manz H. 2004 “Total solar energy transmittance of glass double façades with free convection”,
Energy and Buildings 36, 127–136.
Safer, N., Woloszyn, M., Roux, J.J., 2005 "Three-dimensional simulation with a CFD tool of the
airflow phenomena in single floor double-skin façade equipped with a venetian blind", Solar Energy
79, 193–203.
Stec W.J., Van Paassen A.H.C. 2005, “Symbiosis of the double-skin façade with the HVAC system”,
Energy and Buildings 37, 461–469.
Tanimoto J., Kimura K., 1997 “Simulation study on an airflow window system with an integrated roll
screen”, Energy and Buildings 26, 317–325.
KEYWORDS: air leakage, internal insulation, log house, test-house, corner notch
SUMMARY:
In this study the air leakage of four different log corners and hygrothermal performance of three
different interior insulation materials are studied in a log test house. Field measurements of the air
leakages were analysed in a test house built with different types of log junctions: corner post, dovetail
notch, Scandinavian saddle notch, and double notch with wind lock. One wall made of logs (average
thickness of 270 mm) in a test house was internally insulated with three different materials: cellulose
fibre, mineral wool and reed mat. The air leakage rate of the overall house was also measured using
the standardized building pressurization technique. Hygrothermal performance of walls was measured
with t&RH sensors and heat flow plates.
The mean air leakage rate of the overall house at the pressure difference of 50 Pa was q50=2.8
m3/(h·m2) and mean air change rate was n50=5.0 h-1. The air leakage rate lowered slowly over time.
The lowest air leakage was through the Scandinavian saddle notch - 6.7 m3/(h⋅m).
Simulation models on the hygrothermal performance of the studied test walls were validated based on
the measurement results. Temperature, relative humidity and heat flux showed good agreement
between the measured and the calculated results.
1. Introduction
Log houses have a long history and they represent a variety of building techniques employed in
Estonia and in other Nordic countries. Results of measurements of old log houses have shown that
this building type is typical of the highest air leakage (Alev & Kalamees 2013). In a new log house,
the quality of the envelope has to be higher than a century ago, because the requirements for comfort,
function, and energy-efficiency of today’s residents are different.
Uncontrolled air movement through a building envelope leads to problems related to the hygrothermal
performance, health, energy consumption, performance of the ventilation systems, thermal comfort,
noise and fire resistance. Air leakage through the building envelope depends on the results of the air-
pressure differences across the envelope, the distribution of air leakage places and the airtightness of
the building envelope.
Many studies have analyzed the possibilities to use the internal thermal insulation for improving the
thermal resistance of external walls, including those focused on stone walls (Stopp et al. 2001; Häupl
et al. 2004; Toman et al. 2009) and only few have studied the internal insulation of log walls (Ojanen
2007; Alev et al. 2012; Arumägi & Kalamees 2012; Arumägi et al. 2011). Ojanen (2007) studied the
2. Methods
2.1 Tested house
The field measurements were carried out in a small test house (one room with a net area of 18 m2)
specially designed and built for current study (FIG 1). The house was made of square logs with an
average thickness of 200 mm (except the back wall with half round logs and an average thickness of
270 mm, insulated from the internal side). Every corner had a different type of log junction (FIG 2).
A B C
FIG 1. View of the test house from the back (left) and internally insulated wall (right).
1 2 3 4
FIG 2. Different corner notches of the log wall: double notch with wind lock (1), Scandinavian saddle
notch (also known as Norwegian notch) (2), dovetail notch (3) and corner post (4).
2.3 Simulations
The measurement results were compared with a complex hygrothermal simulation model, WUFI 5.1
Pro. The comparison was made to validate the simulation model for future simulations with different
initial and climatic conditions as well as with different dimensions of the building envelope layers.
FIG 3. Air leakage rate of the total house (left) and the corners (right).
Two methods may be used to compare the thermal performance of the test walls: to compare
temperatures on the inner surface of the log, or to compare the measurement results of the heat flux
plate. The simulation model and this experiment were validated in both ways. The figures below show
the difference between the measured and the calculated temperatures (FIG 4), RH (FIG 5) and water
vapour pressure (FIG 6) between the log wall and the insulation layer for three materials. Thin lines on
both figures represent the measured values and thick lines the calculated values. There is a good
correlation between the measured and the calculated results in addition to the temperature and the RH,
also with the heat flux. Differences in the temperatures reflect the different thermal resistances of the
insulation materials. Reed mat has higher thermal conductivity and therefore the temperatures between
the insulation layer and the log wall were higher during the winter period.
FIG 4. Measured and calculated temperatures on the inner surface of the log wall.
The log wall with an average thickness of 27 cm had a thermal transmittance of U=0.49 W/m2K. After
adding internal insulation, the thermal transmittance decreased to U=0.31 W/m2K in the case of
cellulose fibre, U=0.28 W/m2K in the case of mineral wool, and U=0.36 W/m2K in the case of reed
mat.
The RH level of different materials was different due to different construction methods. Insulation
with cellulose fibre was most critical, because it was installed by a wet spray method. While the
FIG 5. Measured and calculated RH on the inner surface of the log wall.
FIG 6. Water vapour pressure on the inner surface of the log wall.
The required moisture excess was guaranteed over the first half year without the use of air humidifier,
caused by the drying out of the logs and the absence of ventilation. The vertical temperature gradient
was small due to air movement in the room caused by the air-air heat pump.
4. Discussion
Comparison of airtightness measurements of this specially built house (2.8 m3/(h·m2)) and previous
measurements in the log houses in Estonia reveals substantial improvements using the new sealing
method and the quality of work: the average airtightness of 12 log dwellings was 9.2 m3/(h·m2)
(Kalamees 2008); the average of 35 measurements in wooden apartment buildings (made of logs) 10
5. Conclusions
The average air leakage rate of all measurements of the entire envelope at the pressure difference of 50
Pa was q50=2.8 m3/(h⋅m2). The average air leakage rate of all measurements in every corner was 10
m3/(h⋅m). The lowest air leakage was through the Scandinavian saddle notch - 6.7 m3/(h⋅m), almost
the same was through the dovetail notch - 6.8 m3/(h⋅m)), the corner post had much higher leakage rate
of 10.4 m3/(h⋅m) and the double notch with the wind lock was the leakiest - with 17.3 m3/(h⋅m).
The internal insulation parts were constructed differently from the hygrothermal point of view. The
mineral wool part was constructed as a dry wall, cellulose fibre was installed using a wet method and
the reed mat was covered with clay plaster that added also moisture to the reed mat. Cellulose fibre
and mineral wool were covered with water vapour barrier, which prevented these wall parts to dry to
the room side. The reed mat with clay plaster dried out within less than a month. The RH in the
cellulose fibre part was over 80% for over 6 months and the RH level in the mineral wool part had
increased to the same level as in the cellulose fibre part after 6 months and both started to dry out to
the external side during the winter period. During the warm period the RH was about 10% lower than
in other parts, during the cold period the RH level increased at the same rate as in other parts.
The simulation models of three internally insulated log walls were validated using long-term field
measurements. The WUFI was selected for the hygrothermal performance simulations. A good
correlation between the calculated results and the measured values was achieved after the modification
of the material properties and adding a factor as the air change rate in the material layers inside the
wall. The thermal transmittance decreased by 27..43% in different internally insulated wall parts as
compared to an uninsulated wall. Drying and wetting are determined more accurately if the convective
air flow is included in the hygrothermal simulation model. The drying out moisture from fresh, but a
kiln-dried log with an average RH level of 21.6% caused significant moisture excess to the room and
increased the RH in the internal insulation layer. The validated model will be used in our further
studies to focus on the performance of the internally insulated log wall in cold climates.
6. Acknowledgements
The research has been conducted as part of the projects IUT1-15 “Nearly-zero energy solutions and
their implementation on deep renovation of buildings” financed by Estonian Research Council.
Alev, Ü. et al., 2012. Comparison of thermal performance of mineral wool and reflective insulation on
internally insulated log wall. In Proceedings of Healthy Buildings 2012. Brisbane: Queensland
University of Technology.
Alev, Ü. & Kalamees, T., 2013. Field Study of Airtightness of Traditional Rural Houses in Estonia. In
M. L. Karel Kabele, Miroslav Urban, Karel Suchý, ed. Proceedings of CLIMA 2013. p. 6882.
Arumägi, E. et al., 2011. Field study of hygrothermal performance of log wall with internal thermal
insulation. In International Conference on Durability of Building Materials and Components.
Proceedings. pp. 811 – 819.
Arumägi, E. & Kalamees, T., 2012. Validation of a Simulation Model for Hygrothermal Performance
of Log Wall with Internal Thermal Insulation in Cold Climate. In Proceedings of 5th IBPC.
Kyoto, pp. 345–352.
Häupl, P., Fechner, H. & Petzold, H., 2004. Interior retrofit of masonry wall to reduce energy and
eliminate moisture damage: Comparison of modelling and field performance. In Thermal
Performance of the Exterior Envelopes of Buildings IX. Florida.
Hukka, A. & Viitanen, H.A., 1999. A mathematical model of mould growth on wooden material.
Wood Science and Technology, 33, pp.475–485.
Kalamees, T., Arumägi, E., et al., 2011. Eesti eluasemefondi puitkorterelamute ehitustehniline
seisukord ja prognoositav eluiga, Tallinn.
Kalamees, T., Alev, Ü., et al., 2011. Maaelamute sisekliima, ehitusfüüsika ja energiasääst I, Tallinn:
Tallinna Tehnikaülikool.
Ojanen, T., 2007. Low Energy Log Walls Under Cold Climate Conditions. In Proc. Thermal
Performance of the Exterior Envelopes of Whole Buildings X International Conference.
Clearwater Beach, Florida, p. 9.
Stopp, H. et al., 2001. The hygrothermal performance of external walls with inside insulation. In
Thermal Performance of the Exterior Envelopes of Buildings VIII. Clearwater Beach, Florida.
Toman, J., Vimmrová, A. & Černý, R., 2009. Long-term on-site assessment of hygrothermal
performance of interior thermal insulation system without water vapour barrier. Energy and
Buildings, 41(1), pp.51–55.
1 Introduction
It is important that the building envelope is airtight, to prevent outside air to leak into the house and
cause draught, excessive energy use and moisture convection which affects people's health and
satisfaction. Increased energy use is caused by cold outside air that blows into the building which must
be heated to the indoor temperature. Moisture convection can cause damage to the exterior walls and
roof when warm moist air condenses inside the wall and ceiling of the exterior parts that are colder.
Air leakage as the result of poor air tightness can be measured by a pressurization test where a fan is
controlled to keep a 50 Pa pressure difference between indoors and outdoors. The airflow through the
fan is measured and presented with a functional unit that can be for example envelope area or floor
area. The EN-13829 standard states that the resulting value should be weighted from one test with
under pressure and one with over pressure. Simultaneously with the pressurization test, it is common
practice to perform thermal imaging of the interior building envelope to find locations of thermal
bridging and air leakage.
2 Method
For the pressurization tests, according to the standard EN 138 229, all air inlets and devices and other
intentional openings through the building envelope must be sealed. Before the measurements were
carried out, all supply and exhaust devices were taped and, sewage traps were filled with water and all
doors and windows were closed. A Minneapolis Blowerdoor with fan model 4 was used. The fan was
mounted in an outgoing balcony door. The fan was mounted in a Blowerdoor custom canvas that was
stretched on adjustable aluminum rails with a rubber strip on the side facing the door frame to get it
completely tight in the opening.
A pressure meter of model DG-700 measured the pressure difference between indoor and outdoor. A
pc application, Tectite Express 3.1, controled the fan to obtain a set point pressure difference and
another pressure meter gives the airflow through the fan. The standard specifies a number of set point
pressure differences that will be tested up to 50 Pa.
3 Results
During the testing day the weather was sunny, the wind light and the temperature -6 ° C. The
pressurization test showed that apartment number seven which is an apartment on the second floor of a
two-floor house, has an airflow at 60 Pa of 499 l / s, compared to the nine other apartments that have a
maximum airflow of between 82 and 159 l / s. The high value of apartment seven was due to a
construction defect close to the shaft penetration to the attic space. Because of this large leakage in
apartment seven, its result was excluded from further analysis. However, it is included in Table 2 to
show the difference in airflow. Table 2 shows the leakage flows normalized to quantities depending on
the different surfaces and volumes in the apartments. This is to compare the apartments with each
other but also with the results of other studies. Table 1 gives the quantity description used in Table 2.
Table 1 Quantities of Table 2
The study did not measure airflow to adjacent apartments which makes it impossible to compare the
survey results with the requirement from the former building regulation, which specified airflow per
envelope area for the entire building. Leakage airflow in apartments can also enter adjacent
apartments. The flow at 50 Pa per exterior envelope area (Q50, ext) varies between 0.59 (APT 8) to 2.29
B50 / (l/(s·m))
APT 1
7
6 APT 3
5 APT 4 APT 8
APT 10
4
APT 2 APT 5 APT 9
3 APT 6
2
1
0
0 20 40 60 80 100
Figure 1. Leakage airflow B50 per concrete floor slab edge length.
Figure 1 shows apartments’ leakage airflow per concrete floor slab edge length to the outdoor sorted in
order. The apartments belonging to two-floor buildings is marked blue, italic, and the label is
100,00 APT10
APT4
APT8
80,00 APT3
60,00
40,00
20,00
0,00
0 5 10 15 20 25 30 35
Concrete floor slab edge/ (m)
n50 / (m3/(h·m3) )
2,50
APT5 APT6
2,00 APT9 APT2
APT3 APT1
1,50
APT4 APT10
APT8 Two storey houses
1,00
Three storey houses
0,50
0,00
0 20 40 60 80 100
Ratio of apartments %
14.7°C
Figure 6. Right side of the bedroom in APT 1 towards outdoors.
5 Acknowledgements
This study is funded by the Research Council Formas.
6 References
Stein J. (2008). Air tightness in multifamily - measurements and analysis. THID-08/5098 (In Swedish)
Kalamees T. (2007). Airtightness and air leakages of new lightweight single-family detached houses
in Estonia. Building and Environment, 42 (6), pp2369 – 2377
Nevander, L.E. and Elmarsson, B. (2006). Moisture Manual - practice and theory. ISBN 9173331562
Stockholm: Swedish Building Centre Ltd
Kiruna a 2013, Climate. [www]. Retrieved from <http://www.kiruna.se/kommun/Kommun-
politik/Kommunfakta/Geografi/Klimat/>
Kiruna b 2013, Event Handling new city hall. [www] . Retrieved from
<http://www.kiruna.se/PageFiles/6614/134368%20t%C3%A4vlingsprogram_uppslag_low.pdf?epslan
guage=sv>
The Swedish National Board of Housing, Building and Planning 1 2013,
<http://www.boverket.se/Bygga--forvalta/Regler-om-byggande/Boverkets-byggregler-BBR/>
SUMMARY:
The conditions in one zone are often treated as ‘fully mixed’ in Building Energy Simulation programs
(BES) e.g. TRNSYS. By contrast, Computational Fluid Dynamics (CFD) are too complex and time-
consuming to predict temperature and relative humidity in a building for a longer time period, such as
one year. Nevertheless, several building applications require a prediction of the vertical distribution
of the indoor relative humidity and temperature under transient boundary conditions. This paper
presents the development and possibilities of a coupled TRNSYS-zonal model, which allows
accounting the stratification in a room. The zonal model is based on the so-called block model
proposed by Togari[1] and describes a one-dimensional heat and mass transfer between horizontal
layers. The goal of this preliminary study was to examine the effect of parameters such as the number
of layers, the value of the heat convection coefficient and the value for the time step. To validate the
model, the case of natural convection in a test room (3m x 3m x 2.5m ) described by Arai and Togari
has been used .
1. Introduction
‘Energy use’ is a concept which we could not ignore nowadays: every day we are flooded with
information and publicity about energy saving measures and cutting down on power consumption.
However, when dealing with historical buildings, architectural and deontological criteria e.g., the rules
of the Flemish monumental Guard complicate the idea of sustainability. Due to increased thermal
comfort expectations from churchgoers, in recent decades heating systems have been installed in the
historical church buildings. Even today the demand for new or additional heating is still large. So
despite the idea of energy saving and despite the fact that a church building represents a large volume
to be heated, one can ask whether heating a church is ecologically and financially wise.
There also is the additional problem that a lot of churches contain valuable and historic pieces of art.
Hence, the heating systems, which often only operate during service, significantly alter the
microclimate in the church and these fast-changing indoor conditions may result in a faster
deterioration of the artworks. To avoid the damage of the artworks, guidelines can be found in
literature for the acceptable indoor temperature and relative humidity variations (Anon 2011). To
reconcile these preservation needs for the artworks with the heating demand, computer simulations
provide an alternative for experiments and measurements. Because heating often leads to stratification
in this type of buildings, the simulation study must correctly estimate the stratification. Therefore, the
calculation of the airflow in the space is also necessary, next to the calculation of the energy exchange.
To predict this temperature and humidity distribution in building, different modelling approaches have
been used; namely the CFD-method and a zonal airflow model. The CFD-method is a widespread
approach to simulate the airflow in a building. The models based on the CFD-method predict
The mass transfer along the wall is modelled using a wall current model based on the boundary layer
theory for natural convection for a vertical wall (Schlichting and Gersten 1979). The model assumes
that heat convection drives mass flow mout(I,K) with an average temperature TD(i,K) from layer i to its
related boundary layer (
FIG 1). To calculate mout(I,K) with temperature TD(i,K), the following assumptions were made:
TD(i,K) = 0.75Ti + 0.25Tw(i,K) (2)
αC(i,K) .Aw(I,K)
mout(i,K) = 4 Ci
(3)
Where mm(i,K) the mass of the wall current of wall K adjacent to layer i, kg/s
mmd(i+1,K) the mass of the wall current of wall K adjacent to layer i+1 or i-1 to layer i, kg/s
Tm(i,K) the temperature of the wall current of wall K adjacent to layer i, °C
Some of the air returns to the air layer i (𝑚𝑖𝑛(𝑖,𝐾) ) and some continues to the cell down/up(𝑚𝑚𝑑(𝑖,𝐾) ).
The splitting of the mass 𝑚𝑚(𝑖,𝐾) into 𝑚𝑖𝑛(𝑖,𝐾) and 𝑚𝑚𝑑(𝑖,𝐾) is calculated by the ratio P(i,K).
ZONAL Model
- Initial values: T, RH -additional
- Twall wall currents convective heat gain
-Vapour flux of the Mass balance -additional moisture
wall Twall Temperature balance Heat gain gain
Tair node Moisture
RHair node gain
Qconv,air node
TRNSYS
-Geometric model Air node heat and Air node temp
-Boundary conditions moisture balance Air node RH
-radiation model
FIG 2. Schematic overview of the coupling between the multizone building model in TRNSYS and the
developed thermal zonal model.
3. Initial results
3.1 Test case
The case studied in this paper was the case that can be found in the report of Togari et al. (Togari et al.
1993) and that of Arai et al.(Arai et al. 1994) which was used for evaluating the implemented model.
The geometrically simple test room had a ground plane of 3m x 3m and measures 2,5m in height. The
room consisted of insulation boards (three vertical walls, a ceiling and a floor) and one glass wall. In
the wall opposite to the glass wall, two openings were foreseen in the symmetry plane: a supply inlet
at 0,625 m above the floor and a return outlet at 0,250m above the floor. Several configurations were
measured in this test room e.g. air heating or cooling and natural convection.
To validate the calculated results of the newly implemented model, the natural convection case (N10)
has been used in this paper. In this case, the temperature outside the room was kept at about 12°C and
was then raised to 42°C, while 24 hours later it was again lowered to 12°C. The response of the test
room was monitored, while no air was supplied into the room. The moisture transfer between the air
wall surface and the air was not considered in this case study.
The first calculation with the zonal model was a steady-state calculation of the layers temperature
using the measured interior surface temperatures at two hours of testing as boundary conditions.
TABLE 1.Measured temperature for the wall and glass surface and the zone temperature at different
heights (0,25m – 0,75m – 1,25m 1,75m 2,25m) in the test case presented by Togari and the calculated
temperature in the zonal model.
The next step was to perform a dynamic calculation taking transient terms in eq.6 and eq.7 into
account. For that, 48 hours were simulated with coupled TRNSYS - zonal model. The goal of this
study was to examine the effect of physical and numerical parameters such as the thermal resistance of
the walls, the value of the heat convection coefficient, the number of layers and the value for the time
step. In this case, the full test room was modelled.
Because material parameters were unknown, several compositions for the walls were tested. Good
agreement was found for the thermal resistance R=1.6 m²K/W for the insulated walls and with the
glass wall as a massless layer. However, the surface temperature of glass was slightly under predicted.
On FIG 3, results were depicted for the first seven hours in which outdoor temperature was raised
from 12 to 42°C. The measured outdoor temperature served as boundary conditions, while the
temperatures for the insulated walls and the glass surface were calculated by TRNSYS. This was
opposite to the calculated results of Togari in which all interior surface temperatures at each time step
(one hour) were given by the experiment (steady-state calculation for every hour). The calculated
results were compared to the measurements and the calculations of Togari (FIG 3). The zone
temperature was lower than was seen in the calculations of Togari, which was probably due to the
lower surface temperatures of glass and the dynamic character of the calculation (Togari adapted the
boundary conditions every hour based on his measured results).
On FIG 4, the effect of the number of layers on the zone temperature and relative humidity was
visualised. The temperature and the relative humidity in the zone represented the values for whole the
volume. The more layers, the more the stratification was detailed. Further, it could be noticed that the
number of layers has an impact on the zone temperature and the surface temperature of the floor and
the ceiling. The more layers, the higher the temperature and the lower the relative humidity became. In
future, further research needs to be done on the allowed maximal number of layers. Next, the effect of
the time step has also been investigated, which had no effect on the calculated results, and the value of
the heat convection coefficients (FIG 5). In the first calculation, heat transfer coefficients were
considered as simple constant values as described earlier. In the second calculation, heat transfer
coefficients were adapted for every time step based on the temperature in the layer. In first instance
2,5
Height above floor level [m]
measured
2
1,5
zonal model
1 + TRNSYS
0,5 Model
Togari
0
10 15 20 25 30 35
Temperature [°C]
FIG 3:Comparisons between calculated values, measured values and the calculations of Togari in
case where outside was heated for every hour, during seven hours(N10).
2,5 2,5
L3
Height above floor level [m]
2,0 2,0
L5
1,5 1,5
L7
1,0 1,0
L10
0,5 0,5
0,0 0,0
30,0
32,0
34,0
36,0
38,0
40,0
16,0
16,5
17,0
17,5
18,0
18,5
19,0
19,5
20,0
20,5
21,0
T[°C] RH[%]
FIG 4: Comparison between the calculated temperature and relative humidity after two hours of
heating for a different number of layers.
45
40 Outside temperature
Temperature [°C]
35
Ti 0,01 H fixed
30
25 Ti 0,15 H fixed
20
Ti 0,01 H TRNSYS
15
10
4
0
12
16
20
24
28
32
36
40
44
48
Time [h]
FIG 5. Comparison between the calculated values with a time step of 0,01h and 0,15h and between a
constant heat transfer coefficients and the heat transfer coefficients calculated by the algorithms from
TRNSYS
References
Anon (2011) 'ASHRAE handbook: Heating, ventilating, and air-conditioning applications, SI edition'
in, American Society of Heating, Refrigerating and Air-Conditioning Engineers, 23.1-23.22.
Arai, Y., Togari, S. and Miura, K. (1994) 'Unsteady-state thermal analysis of a large space with vertical
temperature distribution', ASHRAE Transactions, 100(part 2), 396-411.
Heiselberg, P., Murakami, S. and Rulet, C.-A. (1998) Annex 26: Ventilation of large spaces in
buildings. Part 3. Analysis and Prediction Techniques, Denmark, Aalborg, Aalborg
University
Megri, A. C. and Haghighat, F. (2007) 'Zonal Modeling for Simulating Indoor Environment of
Buildings: Review, Recent Developments, and Applications', HVAC&R Research, 13(6), 887-
905.
Schlichting, H. and Gersten, K. (1979) Boundary-Layer Theory, Seventh Edition ed., New York and
London: MacGraw-Hill.
SEL, U. o. W.-M. (2013) 'A TRaNsient SYstems Simulation Program TRNSYS 17', [online],
available: http://sel.me.wisc.edu/
Solar Energy Laboratory, U. o. W.-M. (2010) TRNSYS 17: A Transient System Simulation Program,
Solar Energy Laboratory, U. o. W.-M. (2012) Mult izone Building modeling with Type56 and
TRNBuild, Volume 5, Solar Energy Laboratory, University of Wisconsin-Madison.
Steeman, M. (2010) Hygrothermal modelling for building energy simulation applications, unpublished
thesis Ghent University.
Takemasa, Y., Togari, S. and Arai, Y. (1996) 'Application of an unsteady-state model for predicting
vertical temperature distributions to an existing atrium', ASHRAE Transaction, 102(part 1),
239-247.
Togari, S., Arai, Y. and Milura, K. (1993) 'A Simplified Model for Predicting Vertical Temperature
Distribution in a Large Space', ASHRAE Transaction, 99(part 1), 84-90.
Wang, X., Huang, C. and Cao, W. (2009) 'Mathematical modeling and experimental study on vertical
temperature distribution of hybrid ventilation in an atrium building', Energy and Buildings,
41(9), 907-914.
KEYWORDS: air leakage, building envelope, numerical model, permeability, heat loss, defect
SUMMARY:
Air leakage through the building envelope can result in an increase in energy consumption and in
potential moisture damages. In this paper, the development of a simplified coupled heat-airflow model
is presented. The model is applied to a ceiling section insulated with blown-in cellulose, separating an
attic space from a heated indoor space. This ceiling section, part of a full scale test-hut built in an
environmental chamber, is tested experimentally with and without air leakage, where temperature
sensors in the cellulose provide a map of the temperature profile at various locations about the air
leakage point. A comparison of the experimental data and the simulation results is performed.
1. Introduction
Ineffective design and poor workmanship can both lead to a dramatic decrease in building’s air
tightness. Occupants can also inadvertently create additional openings in the building envelope after
the building is delivered. Excessive air leakage can result in significant heat losses, as well as
condensation formation and mould growth in the building envelope. Statistical surveys carried out in
France over the past decade identified the areas of the building envelope where air leakage is more
often observed (Litvak, 2005). Excluding electrical conduits and thin gaps between window frames
and the walls, airflow due to air leakage through building defects generally flows through the
insulation material. The model developed here aims to assess the effect of the airflow on the
temperature profile in porous insulation materials. Following a brief overview of the existing research
work on heat-airflow coupling, the governing equations, geometry, and boundary conditions of the
model are presented. A comparison between the experimental data and the simulation results
demonstrates the ability of the model to predict the temperature conditions in the insulation.
2. Numerical model
2.1 Current state
The influence of air leakage on the performance of porous insulation has been investigated since the
late eighties. An experiment conducted at the Fraunhofer Institute of Holzkirchen showed that a 1 mm
crack in an air membrane placed on fibrous insulating material could reduce its effective R-value by a
factor of 4.8 (Wagner, 1989). Langlais et al. (1990) presented a numerical model to assess natural
The momentum conservation for laminar airflow is described by the Navier-Stokes equation, which is
highly nonlinear and often solvable only with numerical methods.
∂v
ρ air + ρ air (v ⋅ ∇) v = −∇P + µ Δv − ρ air g (4)
∂t
The pore Reynolds number provides information on the flow regime in the porous medium and is
expressed as
Lc ρ air u
Re p = (7)
µ air
The characteristic length (Lc) for the flow can be approximated by the square root of the material
permeability (Nield and Bejan, 2006).
Lc = kmat (8)
In building physics, the infiltration/exfiltration air velocity driven by pressure differences typically do
not exceed 0.5 m/s (Chan et al., 2003). For our experiment, given that the permeability of cellulose
insulation is 7.67 × 10-10 m2 (Trechsel, 2011) and the inlet airflow rate is 2 or 5 LPM (see section 3):
Re p = 4.8 → 1 < Re p < 10 (9)
For Rep>10, the flow regime is still laminar but the friction drag becomes quadratic as described in the
Forchheimer equation:
ρ air $ ∂u u' µ
& + (u ⋅ ∇) ) = −∇P −
−1/2
u − cF kmat ρ air u u + µ Δu (10)
ε % ∂t ε( kmat
3. Experimental Setup
A test hut built inside an environmental chamber at Concordia University provides the experimental
setup for this study. The test hut contains an unvented attic space insulated at the ceiling level with
380 mm of blown-in cellulose above 38 mm of rigid polyisocyanurate (PIR). The PIR also acts as the
air and vapor barrier of the ceiling. To simulate air leakage, sampler pumps installed in the indoor
space are used to continuously deliver indoor air into the attic space at a controlled rate. The air is
supplied at the bottom of the cellulose insulation through an orifice in the PIR using a tube with a 6.4
mm inner diameter. In the case where there is no airflow into the attic, the orifice is sealed with tape at
bottom of the PIR.
FIG 1. Location of temperature sensors FIG 2. Temperature sensors held in place with a
across the ceiling low-profile wooden support system
Figure 1 shows the location of the temperature sensors in the indoor, ceiling, and attic space. Tatt and
Tint are the attic and indoor space temperature sensors, respectively. A way to indirectly map air
leakage through building components is to measure temperature in the vicinity of the hole (Desmarais
et al., 2000), so the remaining sensors are located within the cellulose insulation at three different
heights and at various horizontal distances away from the air leakage orifice in the PIR. Tatt, Tint, T3A
and T6 are resistive temperature detectors embedded in relative humidity sensors with an accuracy of
±0.3ºC. The other temperatures are measured by thermocouples (30AWG, NSB special limits of error
4. Numerical simulation
The model in this study uses COMSOL's built-in physics of "heat transfer in porous media" and
"Brinkman equations" to implement the conservation equations in the different domains as described
in section 2. To reduce computational time, the transfer phenomena are considered to be rotationally
symmetric about the longitudinal axis that starts at the air orifice and points upward, neglecting
therefore the thermal bridge induced by the truss and wooden structures supporting the sensors (Figure
2). The simulation is therefore performed on a 2D axisymmetric plane, creating the 3D rotational polar
coordinate system.
TABLE 1. Input Parameters
Parameter Symbol Units Value Source/Note
Thermal conductivity of cellulose λcel W/(m·K) 0.038 Manufacturer
Thermal conductivity of PIR λpoly W/(m·K) 0.022 ASHRAE
Thermal conductivity of air λair W/(m·K) 0.026 At 20°C
Dynamic viscosity air µair Pa·s 1.8 × 10-5 At 20°C
Surface film coefficient int. hint W/(m2·K) 9.26 ASHRAE
Surface film coefficient attic hattic W/(m2·K) 4.32 ASHRAE
Permeability of cellulose kcel m2 7.67 × 10-10 (Trechsel, 2011)
Porosity of cellulose εcel - 0.95 CSTB (Optimob)
Air density ρair kg/m3 1.2 At 20°C
Thermal capacity air cpair J/(kg·K) 1006 Assumed constant
Boundary conditions for air and heat must be chosen on each boundary. The measured temperatures
Tint and Tatt are used as boundary conditions for the model. To simplify the geometry, polyisocyanurate
insulation has been integrated as thermal boundary condition with the interior surface film coefficient
in an equivalent global heat transfer coefficient (Figure 3):
1
heq = (13)
1
Rpolyiso +
hint
The interior temperature is imposed at the entrance of the hole, and both open boundary and natural
convection are assigned to the top of the cellulose. The heat flux brought by air can both enter and
leave the domain on interfaces with the thermal open boundary condition (COMSOL, 2011):
%T = Tatt if n ⋅ u < 0
& (14)
'−∇T ⋅ n = 0 if n ⋅ u ≥ 0
FIG 3. Domain physics and boundary FIG 4. 3D plot of the temperature field for a flow
conditions of simulation rate of 5 LPM
TABLE 2. Temperature data from the experiment and simulation, values are in ºC
Experiment Simulation
Data Point No airflow 2 LPM 5 LPM No airflow 2 LPM 5 LPM
Tint 21.9 21.8 21.8 21.9 21.8 21.8
T1 16.6 20.3 20.8 17.7 21.5 21.8
T2A 12.1 16.1 19.4 13.1 15.9 18.9
T2B 12.1 15.7 19.0 13.1 15.5 18.3
T3A 6.4 6.8 7.9 7.2 7.5 8.0
T3B 6.4 6.6 7.5 7.2 7.5 8.0
T3C 6.4 6.6 7.2 7.2 7.4 7.8
T3D 6.3 6.4 6.8 7.2 7.4 7.6
T4 16.3 16.9 17.3 17.7 17.9 18.2
T5 12.0 12.6 13.3 13.1 13.4 13.9
T6 6.4 6.5 6.6 7.2 7.3 7.4
Tatt 6.7 6.8 6.8 6.7 6.8 6.8
6. Conclusion
Air leakage is a major contributor to moisture problems and energy losses in buildings. Airflow
through the building envelope often occurs through porous insulation materials, and understanding its
effect on the enclosure performance is of importance. This study presents a simplified, coupled heat-
airflow model, applied to an experimental setup. The experiment focuses on the temperature profile of
cellulose insulation in an attic of test hut, subjected to a 15ºC temperature difference and various air
leakage rates. The results from the experiment and simulation show good agreement. Ongoing
research aims to expand this model to couple moisture transport, which can further increase its
applicability, especially for greater temperature differences and transient processes. Parametric and
sensitivity analyses with multiple orifice areas, thicknesses of insulation, and flow rates could also be
performed.
7. Acknowledgements
This work is financially supported by ADEME (Agence De l'Environement et de la Maîtrise de
l'Énergie), CSTB (Centre Scientifique et Technique du Bâtiment), and the Région Rhône-Alpes. The
experimental data is taken from a project under the NSERC Smart Net-Zero Energy Buildings
Strategic Research Network (SNEBRN), sponsored by the Natural Sciences and Engineering Research
Council of Canada (NSERC) and 14 industrial partners including KOTT Group. The project is also
supported through a NSERC discovery grant. The materials and installation of the test hut are supplied
by KOTT Group.
SUMMARY:
The airtightness of office and school buildings influences energy use and thermal comfort. A leaky
building is likely to have a high use of energy and thermal discomfort. The knowledge of real
airtightness levels of entire buildings and their impact on the energy use is very low, except for a study
carried out in the USA. Therefore two different methods of airtightness testing were applied to six
entire Swedish office and school buildings built since 2000. The first method involves using the
ventilation system of the building and the second one to use a number of blower doors. Information on
30 other airtightness tests of similar buildings was collected. During the airtightness testing the air
leakage paths were detected using infrared scanning and smoke sticks. Most of the tested buildings
showed a very good airtightness level. All previously tested office buildings in the USA, Canada and
the UK are much leakier. The tested buildings had some leakage paths, which could easily have been
taken care of during construction, but are rather difficult to stop now. The infiltration rate for these
buildings was estimated by a simplified calculation procedure to be low, implying a low energy use
caused by infiltration. For most of the buildings stringent airtightness requirements had been applied.
1. Introduction
It is well-known that the building sector plays an important role in the work towards sustainable
development. The sector represents extensive economic, social and cultural values, at the same time as
it causes extensive environmental impact due to its high use of energy and materials. An important
part of the energy use within the building sector is related to office and school buildings. The total
energy use of an average Swedish office building is 220 kWh/(m²year) (heated usable floor area) of
which electricity stands for 108 kWh/(m²year). Of this 108 kWh/(m²year), 57 kWh/(m²year) is due to
office equipment, of which 23 kWh/(m²year) is lighting. This was shown in a study of 123 office and
administrative buildings of different ages (Persson 2007). Of the floor area in all office buildings, 69
% is heated by district heating and the average use of district heating energy is 110 kWh/(m²year)
(SCB 2006). Both new and old office buildings have a substantial potential for energy savings and
improvement of indoor climate. While many new office buildings may have a low energy use for
heating compared with older office buildings, they may have a higher electricity use. This is due to a
high use of electricity for ventilation, cooling, lighting and office equipment. The situation for school
buildings is similar, with the exception of cooling which is rare in school buildings. An important
parameter affecting the energy use for space heating and cooling, and thus the indoor climate, is the
airtightness of the building envelope. In a leaky building the energy use increases due to uncontrolled
infiltration/exfiltration. The air leaking in and out through the building envelope increases the energy
use as it, for example, does not pass through a heat recovery unit. The uncontrolled air leakage can
contribute to discomfort such as draught, which can result in the indoor temperature being raised to
improve the comfort, causing a further increase in energy use from the user’s behaviour.
2. Method
The hypothesis is that, in many cases, the airtightness can be measured using the ventilation system of
the building instead of using a number of blower door units. Two different methods were used:
- Airtightness testing using the ventilation system of the building. Canadian standards were
applied (CGSB 1986, CGSB 1996).
- Airtightness testing using a number of blower doors (portable fans),
www.energyconservatory.com. European standard 13829, Method B was applied (CEN
2000).
The measurements involve pressurizing and depressurizing the entire building and measuring the
corresponding air flow to maintain the different pressure differences between inside and outside.
Ventilation openings and lead-throughs are sealed before the measurements. Thus the airtightness of
the building envelope is determined. The location of leakage paths are determined using thermography
and smoke.
When using the ventilation system of the building the following has to be investigated before
proceeding:
- Assess the building automation system to ensure that the ventilation air flows can be
controlled and that it has the capacity to provide the necessary air flows. It is usually easier if
the building has a demand controlled ventilation system.
- Ensuring that the air flows can be measured and that it can be done with adequate accuracy.
Within this project three schools and three office buildings were tested. Additional tests had been
carried out before by other Swedish organizations.
To determine the air infiltration/exfiltration rate from the results of pressurization tests there are
different ventilation models. The ventilation models can be divided into: “air change” methods,
reduction of pressurization test data, regression techniques, theoretical network methods, simplified
theoretical methods (Liddament 1986). The first three models are empirical techniques, which tend to
be loosely based on the physical principles of air flow. The other models are theoretical models, which
are based on a much more fundamental approach involving the solution of the equations of flow for air
movement through openings in the building envelope. Empirical methods are usually straightforward
to use, but tend to be unreliable and have a limited field of application. On the other hand, theoretical
models have a potentially unrestricted applicability but are often demanding in terms of data and
computer execution time. Theoretical calculation techniques can be divided into: single zone network
models, multi zone network models and simplified theoretical techniques. These models require a lot
of information e.g. wind pressure coefficients, air leakage distribution for the building envelope, local
wind speed, geometry of the building. Due to the limited amount of information on the tested
3. Tested building
The aim was to test school and office buildings built after the year 2000 with a floor area preferably
larger than 1000 m². They should be a mix of buildings both with and without specific airtightness
requirements.
4. Results
All the buildings tested outside the project were very airtight (see table 1). The average airtightness
was 0.3 l/sm² @ 50 Pa which is equivalent to the voluntary Swedish requirement for passive houses
(Nollhus 2012). The best building had a value of 0.1. For most of the buildings airtightness
requirements ranged from 0.2 to 0.8 l/sm² @ 50 Pa, which can be compared with the requirement of
the previous Swedish building code (before year 2006), of 1.6 l/sm² @ 50 Pa. The current building
code does not have any specific requirement. Only two buildings did not meet their requirement. All
TABLE 2. Measured air leakage and calculated energy use for heating infiltrating air.
Type of Year Measured Measured Infiltration/e Energy use for heating
building airtightness airtightness xfiltration, infiltration,
l/sm² @ 50 ach @ 50 ach kWh/m²year
Pa Pa
Shop 2011 0.18 0.20 0.01 2
Office 2008 0.34 0.60 0.03 3
Office 2010 0.43 0.43 0.02 3
Food store 2011 0.62 1.11 0.06 10
School 2010 0.09 0.13 0.01 1
School 2011 0.17 0.34 0.02 2
School 2011 0.48 1.18 0.06 7
School 2010 0.4 1.04 0.05 7
School 2011 0.16 0.24 0.01 2
School 2010 0.88 1.28 0.06 9
Home for the 2012 0.20 0.20 0.01 1
elderly
Home for the 2011 0.14 0.18 0.01 1
elderly
Average 0.34 0.58 0.03 4
Also the recently tested six buildings were fairly airtight, but not as airtight as the previously tested
buildings (see table 3). One contributing factor might be that there were only two buildings which had
a specified airtightness requirement.
5. Conclusion
Two different methods of measuring the airtightness of entire buildings have been used, using the
building’s ventilation systems or using a number of blower doors. Both methods can be used
separately or combined. Choice of method depends on the prerequisites of the individual building. For
large buildings using the ventilation system can be preferable. This presupposes that the air flow can
be controlled within the desired range and measured accurately enough. Usually this is the case for
demand controlled ventilation. Tests during construction, which are recommended to ensure good
airtightness, can often only be carried out using blower doors. The two methods can be applied to
office buildings, apartment buildings, industrial buildings and other premises. For apartment buildings
the blower door technique is often the only method as the ventilation system often has insufficient
capacity, unless the building is very airtight. Complete testing includes determination of the location
of leakage paths.
This study clearly shows that it is possible to build very airtight school and office buildings i.e. which
fulfils the Swedish criterion for passive houses at 0.3 l/sm² @ 50 Pa. This seems to be possible with a
building envelope which is prefabricated or framework construction or curtain wall. Most likely, the
energy use for infiltration in these buildings is almost negligible i.e. in the order of magnitude of a
couple of kWh/m²year. This number can be compared with the total energy use for space heating for a
typical average Swedish office building of 110 kWh/m²year, where infiltration might account for 10-
20 kWh/m²year if only the airtightness requirement of the previous building code is fulfilled, which is
likely.
An important prerequisite for arriving at a very good airtightness is obviously that requirements
regarding airtightness are set at an early stage. These requirements have to be specified as a maximum
air leakage value, which should be tested during the building process.
6. Acknowledgements
The project was funded by SBUF (Development Fund of the Swedish Construction Industry), NCC
and WSP.
7. References
ASHRAE 2009. ASHRAE Handbook Fundamentals. ASHRAE, Atlanta, USA.
Blomsterberg, Å., 1990. Ventilation and airtightness in low-rise residential buildings - Analyses and
full-scale measurements. Swedish *Council for Building Research, D10:1990, ph.d. thesis.
KEYWORDS: Air tightness, air leakage, blower door, pressure test, building envelope
SUMMARY:
It is important to approach issues concerning building air tightness systematically and implement this
throughout the building process to ensure that air leakage requirements are met for large buildings.
Air leakage testing needs to be applied early on during construction to better understand how to
increase the quality of the building and efficiency of the construction process. This paper explains how
to work with air tightness management in the different construction phases and gives a practical
example of a large building that met high ambitions for air tightness.
1. Introduction
Several benefits can be derived from testing air tightness in the early stages of constructing large
buildings. The key to making this possible is to use a systematic approach to managing issues
concerning air leakage throughout the construction process.
Difficulties that occur concerning air leakage testing during the production phase of large buildings
are foremost of the practical types. Examples include limitations of the equipment capacity, the scale
of the provisional sealing methods and time consumption of the capsuling techniques and other
operation steps. Different parties may participate in different ways to contribute towards an air-tight
end result, through defining requirements, technical collaboration meetings or workshops,
coordination for difficulties in practice, storage and use of sustainable materials and well thought-out
modifications.
This paper aims to answer how to efficiently manage air tightness issues from experience of working
in projects for the construction of large buildings. First we present the authors’ general method of
practice for implementing air-tightness management in the different stages of the building process.
Furthermore a project with a high ambition level is presented where the actions taken concerning air
tightness are mapped out and how the issued air leakage requirements were met.
2.1 Air tightness management during the design and development phases
A project developer can set an ambition level for air tightness that can be regulated as a part of the
company’s standard for construction or it can be project specific. Air tightness in this phase is often
related to energy efficiency, which has a strong connection to the attractiveness of the building as a
product, i.e. if it can be promoted as a “green” building by itself or with an environmental
classification system.
The areas of interest for testing were marked out on overview plans and the contractor was given
instructions on the method for constructing the provisional encapsulations, which were also explained
in detail beforehand on site together with the carpenters.
FIG 3. Blower-door equipment attached to a provisional wall with a smoke machine placed at the fan
inlet.
3.5 Skylights
The skylights are situated over the atrium in the core of the building with a vast open space
underneath, see figure 5. The provisional encapsulation that was constructed around one of the
skylights was made up of a customised polyethylene tent that was sealed to the floor of the
scaffolding. The air tightness engineer and representatives from the contractor and subcontractors
enclosed themselves inside the testing area during the testing. The fan was mounted to the
polyethylene foil and the test area was pressurised and inspected for air leakage with the use of
thermography, anemometer and a smoke machine. It turned out that the attachment of the capsulation
to the scaffolding could not ensure an air tight test result during pressurisation but air leakages could
still be found and repair work was prescribed.
4. Conclusions
There are several benefits for air tightness management throughout the building phases and to scale
down planning for manageable testing methods to be applied during the course of the construction.
One of the foremost reasons for these benefits is that the air tightness engineer or consultant may
implement systematic routines for spreading information to the different organisational parts of the
building process, so that the end installer or carpenter is up to date with accurate knowledge on how to
effectively work with these issues and why they are important. There are vast possibilities to engage
these issues early on to affect the building process and have an impact on the end result.
Work that is needed to ensure a comprehensive grip on air tightness is implemented during the design
phase, when construction designs are made that are critical to air tightness issues and need to be
identified and analysed in an Air Tightness Description, for issues that concern both the construction
and operational phase of the building. Before commencing the production, a Plan for Air Tightness
needs to be established, which describes the execution of critical steps that are practically explained
through technical meetings. The choice of suitable building sections to be tested is of great
importance, especially for buildings where it is not applicable to pressurise the total area of the
building envelope at once, for results from these sections will represent the whole building. It is
important to work with air leakage surveillance and apply pressurised testing to strategically chosen
parts of the building during the course of construction, to assure the quality of execution so that the
end result amounts to less than the design value. Through a systematic approach to air tightness
management throughout the construction process and documenting the path taken, lessons may be
learnt and a more widespread knowledge of the importance of air tightness may infiltrate the building
industry.
5. Acknowledgements
The authors wish to thank Fabege for involving us early in the project Uarda 5 and Peab Sverige AB
for the partnership surrounding air tightness testing.
References
CAN/CGSB-149.15-96. Determination of the Overall Envelope Airtightness of Buildings by the Fan
Pressurization Method Using the Building’s Air Handling System, Canadian General Standards
Board, Ottawa, Canada, 1996.
SS-EN 13187. Thermal performance of buildings - Qualitative detection of thermal irregularities in
building envelopes - Infrared method. Swedish Standards Institute SIS, 1999.
SS-EN 13829. Thermal performance of buildings - Determination of air permeability of buildings -
Fan pressurization method (ISO 9972:1996, modified). Swedish Standards Institute SIS, 2000.
SUMMARY:
Air-tightness of constructions is essential when evaluating a building’s energy efficiency and
durability. Even in connection with refurbishment air-tightness is becoming increasingly important.
Especially when using facing formworks with vapour barriers as interior insulation to thermally
renovate an existing building, air-tightness is an essential parameter. It has a significant influence on
the proper functioning of the construction because convective moisture entry can lead to considerable
damage.
In this paper an in-situ measuring carried out on a real object is being presented, which allows to
assess and to proof the air-tightness of facing formworks. This method was developed on a testing
wall in the laboratory of the Research Centre for Building Physics and Sound Protection of Vienna
University of Technology.
The measurement starts after the completion of the construction by using a tracer-gas (CO2) and by
adjusting excess pressure in the testing room. At the same time the increase of CO2 in the facing
formwork is recorded. Measurements are carried out with at least two different pressure ranges. From
the results gained, the air flow volume coefficient C and the leakage exponent n are calculated in
accordance to EN 12114.
1. Introduction
A construction’s air-tightness substantially influences its functionality as convective moisture entries
contain a high risk of causing damage. The quality of the airtight layer often does not meet the planned
requirements. This has been pointed out in many publications. For example, Sandberg (2005)
described typical problems caused by poor air-tightness of constructions and specified experience
from site visits. Wahlgren (2010) presented a method detecting air leakages at an early construction
stage by using negative pressure. However, the given examples of joints and connections solely refer
to the entire building envelope as an airtight structure.
A building’s air-tightness can be determined by negative- or excess-pressure-measurement using
standardised measuring methods (Blower-Door-measurement according to EN 13829). Thereby, the
total air-tightness at a difference in pressure of 50 Pa is given as a n50-value in air exchange rate per
hour or as a q50-value in m³/(m²h) according to the formulas (1). There is, however, no differentiation
whether the leakiness occurs in the area of the windows, doors or joints and connections. Because of
this, size and form of leakages cannot be quantified with this method. Furthermore, a judgement
concerning a construction’s durability cannot be drawn from these measurement results.
V50 V50
n50 q50 (1)
V AE
If, on the other hand, the air-tightness of a single building component or a single construction layer is
to be determined, a different measurement concept must be used. In this case the focus is not on
FIG. 1 left: reveal with facing formwork as interior insulation; right: ceiling-to-wall connection with
facing formwork as interior insulation
It is noticeable that the air-tightness of the facing formwork depends on at least two layers. In the one-
dimensional case these are the vapour barrier within the construction and the exterior load-bearing
wall, e.g. the masonry. With multi-dimensional component connections flow-paths to neighbouring
rooms or storeys may exist. As a prerequisite for a convective moisture entry into the facing formwork
there has to be an air-path from the facing formwork to the outside or a neighbouring room
additionally to the crack in the vapour retarder. Thus, to determine the airflow through the facing
formwork, all leakages in the construction’s interior and exterior cladding must be scrutinized. The
volume flow rate into the construction in essence depends on the leakage size and the difference in
FIG. 2 facing formwork in laboratory test rig; left without foil, right with foil and indication of fields
The facing formwork is lined with a foil which is stuck to the existing wall around the test
construction in a manner as to provide an airtight space. Furthermore, the foil is attached to the facing
formwork’s vertical uprights with wooden strips to simulate planking with plasterboard. Deliberately,
no full planking was provided and see-through foil was used in the laboratory test rig for us to be able
FIG. 3 left: hose leading from facing formwork to the outside; right: schematic depiction of the test
set-up (layout)
During the test excess pressure is set up in the measuring room using a Blower-Door device. At the
same time, CO2 is released. The concentration of CO2 in the room is set to approximately 2.000 ppm
and kept on a constant level. Because of the difference in pressure between the test room and the
facing formwork resp. the exterior the CO2-air-mixture runs into the facing formwork through leaks in
the foil, thus increasing CO2-concentration in all nine measuring fields. In this test, the connection
hose attached to the facing formwork simulates the leakage towards the exterior and thus favours the
mass through cracks on the inside.
If one looks at a period of time with approximately constant CO2-concentration in the test room, the
mass entry ṁ into the facing formwork can be calculated in kg/d. This is effected by determining a
linear regression from the CO2-concentration measured in the facing formwork and by averaging over
the period of time with approximately constant CO2-concentration in the test room.
c
m V C p
b
cR cc (3)
Where ṁ massflow rate of air into cavity in kg/s
c rise of CO2-concentration in test field in facing formwork in ppm/d
cR CO2-concentration in test room in ppm
cc CO2-concentration of test field in facing formwork construction in ppm
V volume of the air space in test field in facing formwork in m³
air density in kg/m³
p Pressure difference in Pa
C air flow coefficient non-dimensional
b air flow exponent non-dimensional
FIG. 6. Leakage in field 5; left with insulation in field 5; right without insulation in field 5
Based on the CO2-progressions it becomes clear that the measurement of hollow facing formworks
provides a solid mean regarding the actual airtightness of a wall-area. When measuring facing
formworks with insulation, however, heavily varying results are gained. This is presumably due to
flowpaths and cavities within a construction filled with insulation material causing an uneven spread
of the tracer gas.
FIG. 7. Schematic view of wall area with three test fields and possible flowpaths both taken and not
taken into account. Each test field contains multiple sensors to detect flowpaths.
In FIG. 8 some results of the in-situ measurements are presented. In this case three test fields – one in
the reveal and two in the wall left and right to the window – were analysed. The results show that the
increase of CO2 in the window reveal (field 3) and in the neighbouring test field 1 are approximately
the same. It stands to reason that a connected plenum exists between these two test fields, whereas test
field 2 shows no connecting flowpath to the other test fields. The measurements were done at three
stages of pressure differences. Thus, the coefficients C and n could be determined.
testfield 2
5 E -7 1,587
testfield 3
2 E -6 1,786
FIG. 8. CO2-concentration in three fields of facing formwork (left); Coefficients C and n (right)
4. Conclusions
Measurements in the laboratory test rig were conducted with and without insulation in the facing
formwork. The results clearly show that the insulation material must be removed from the construction
to obtain a defined volume of air which then serves as basis for evaluation. Furthermore, it was shown
the air-volume in the test field examined in the facing formwork has to be thorough mixed with small
fans to achieve an average airtightness. The measurements in the laboratory test rig have also shown
that CO2 entering the facing formwork does not only depend on leakage size and difference in
pressure, but also on the dimension of the connection hose between facing formwork and outside air
(i.e. the outer leakage). Because of this connection the mass-flow into the construction can be depicted
in dependence of the relation between inner and outer leakage. This thesis will be the topic of future
laboratory tests.
The in-situ measurement on a real object has shown that all test fields of a cohesive facing formwork
provide similar results. It may thus be assumed that the tracer gas does spread immediately after
entering the construction and forms flowpaths to neighbouring rooms. This means that more test fields
must be arranged for more air-tight constructions to gather a statistic figure for the construction’s air-
tightness. The sizing of the artificial outer leakage in the form of a hose makes for further need for
research.
References
Bednar T. & Deseyve C. & Jung M. & Nusser B. & Teibinger M. 2010. Impact of Airflow on the Risk Assessment of Flat
Roofs. in: 7th International Conference on Indoor Air Quality, Ventilation & Energy Conservation in Buildings -
IAQVEC 2010, Syracuse University, (2010)
Hall M. & Hauser G. 2003. In situ Quantifizierung von Leckagen bei Gebäuden in Holzbauart.
Sandberg P.I. & Sikander E. 2005. Airtightness issues in the building process. Proceedings of the 7th Nordic Symposium on
Building Physics NSB 2005
Wahlgren P. & Sikander E. 2010. Methods and Materials for an Airtight Building. Buildings XI, 2010
Wegerer P. & Neusser M. & Bednar T. 2012. Auswirkungen der Luft(un)dichtheit auf die Feuchtebelastung von
Konstruktionen mit Innendämmsystemen. 7th International Symposium on Building and Ductwork Air Tightness in
Practice May 11 – 12, 2012, Stuttgart, Germany
EN 12114:2000 Thermal performance of buildings – Air permeability of building components and building elements –
Laboratory test method
EN 13829:2001 Thermal performance of buildings – Determination of air permeability of buildings – Fan pressurization
method
SUMMARY:
An airtight thermal envelope is important to achieve buildings with high energy efficiency and
moisture safety. The airtightness is commonly measured using the fan pressurization method, and
international studies have shown large variations in the measured airtightness at different seasons. In
this paper the airtightness variation at different seasons is studied using measurements on two
wooden frame, one-family residential buildings. Numerical simulations are presented that show the
variation in airtightness due to varying air densities of the air.
The airtightness measurements in the two buildings have been performed during seven months, from
summer to winter. The trend in the measurements is that the airtightness is lower (more leaky
envelope) when the indoor air is drier (low relative humidity). Consequently, the air leakage is
largest during the winter measurements. The decrease in airtightness during the measurement period
(from summer to winter) is in the order of 8-10%.
1. Introduction
The airtightness of a building has an impact on the energy use and on the moisture safety of a
building. It also affects the thermal comfort, the air quality in a building, sound insulation and the
spread of fire gases (Sandberg et al. 2007a and 2007b). Measuring the airtightness of a building has
become more common lately, much due to the increased energy use in leakier buildings. The
airtightness of a building can be measured in order to attain a certification or on demand from a
developer. The consequences if failing the target can sometimes be severe. Therefore, it is of great
importance to obtain a correct and representative measure of the airtightness.
The airtightness in a building is created by a continuous and airtight thermal envelope. The airtight
layer in a thermal envelope can be either a thin layer, such as a polyethylene foil, a board, such as
plywood, a homogeneous construction (e.g. a concrete component) or an outer coating, such as
rendering. In all examples it is of great importance that the joints are properly sealed.
Airtightness measurements are usually performed in accordance with EN 13829:2000 (Fan
pressurization method). In this standard there are limitations with respect to the climatic conditions
during measurements. There is for example a limit on the maximum allowed wind speed and the
maximum allowed temperature difference over the thermal envelope. The purpose of the limitations is
to assure a correct measured airtightness. Nevertheless, measurements have shown that there is a
variation in the measured airtightness with respect to the time of year for the measurement. Yoshino
(2012) described variations of ±20% over the year. Boorsboom et al. (2012) analyzed airtightness
measurements, from the 80ies, made on 21 window frames mounted in masonry or concrete walls.
The average difference in air tightness between summer and winter was about 30% (higher leakage
during winter) and the maximum seasonal difference was 120%. To be noted, some window frames
showed a lower leakage rate during winter. Boorsboom et al. suggest measurements during three
subsequent seasons in order to obtain correct values. Also Kim and Shaw (1986) showed increased
leakages during winter time. The highest leakages occurred in winter and early spring, and the lowest
2. Methods
The air tightness’ variation over season and climate will be studied by measuring the air tightness of
three buildings, during one year, and by performing numerical simulations on the climate and the
effect on airtightness. The first measurements, on two one-family wooden buildings, started in June
2013. Measurements on a multi-family concrete building will start March 2014. The measurements
have been performed by SP, Technical Research Institute of Sweden. Initial numerical calculations
have been made on the influence of wind, and on air properties.
The first airtightness measurements were performed on two residential one family houses, both
located in the south west part of Sweden, one house in Landvetter and one in Sevred, located outside
Borås. The houses are light weight wooden houses in two floors (plain wood in Landvetter and light
weight wooden beams/joists in Sevred) and they both have slab on ground and cold attics. Both
buildings have mineral wool insulation and polyethylene foil on the inside (between insulation and
board) as air barrier and moisture barrier. The house in Landvetter is built in 2004 and the house in
Sevred is built in 1993. Both houses have mechanical exhaust ventilation systems.
The airtightness quantity used is air permeability, q50 (l/sm²). It is the amount of air that passes
through the thermal envelope at a pressure difference of 50 Pa, per area of thermal envelope (l/sm²).
The air flow is measured both when the building is pressurized and depressurized and the mean value
is used. The airtightness measurements are made according to EN 13829:2000, using a Minneapolis
BlowerDoor. Temperature and relative humidity, indoor and outdoor, is measured at each
measurement occasion, as well as the outdoor wind speed. Air leakage search was performed during
the first airtightness measurements (June 2013) using air velocity meter and thermal camera, and will
be performed again during winter conditions. The airtightness measurements are performed
approximately every 1.5 month in order to study the different seasons and climate conditions.
3. Numerical simulations
Simulations have been performed in order to study how the measured air flow is affected by different
densities of the air. The density of the air is different due to variations in temperature and relative
12 µ ⋅ L
Sg = (2)
b2 ⋅ A
Here S g (Pa∙s/m³) is the air flow resistance, L is the length of the air gap (m), µ (Ns/m²) is the
dynamic viscosity, b (m) is the height (or width) of the gap and A (m²) is the entrance (or exit) area. In
order to include the entrance and exit pressure losses, the air flow resistance at entrance and exit S e
needs to be known. The air flow resistance S e depends on the air flow which requires the use of a
parameter S e' (Pa/(m³/s)²).
1.8 ⋅ δ a
S e' = (3)
2 ⋅ A2
The density of the air δ a (kg/m³) is included. Total air flow through a crack is determined by
1
Ra = ( S g2 + 4 ⋅ ∆P ⋅ S e' − S g ) (4)
2 ⋅ Se
'
As previously mentioned, the leakages are evenly distributed over the building. By investigating
extreme temperatures and relative humidities, the factors that possibly influence the airtightness
measurements are determined.
The investigation shows that the temperature of the air can affect the measurements by affecting the
air flow through the air gaps, see Table 1. Using -20°C, the measured airtightness is 0.76 l/m²s and at
a temperature of 30°C the measured airtightness is 0.81 l/m²s. Consequently, the difference between
the two extreme measurement situations is 5.7%.
v2
∆P = (C p − C pi ) ⋅ (5)
2
The building is first simulated with equally distributed leakages, and then with a windward side that is
twice as leaky as the other sides.
The results from the simulations show that there in many cases is a small difference in the measured
airtightness values for pressurization and depressurization when wind is present. However, the
average value is not affected until the wind speed increases. For example, at a wind speed of 9 m/s,
the building is estimated 2% more airtight with wind than without wind. At a wind speed of 9 m/s, the
zero flow pressure difference is most likely exceeding the value accepted in EN13823:2000. In the
simulations, a higher wind speed resulted in lower calculated air permeability (more airtight building).
For the case of a non-uniform air leakage distribution, simulations were made for a wind speeds up to
9 m/s. The results are similar to those of the equal leakage distribution, i.e. unless the wind speeds are
high there is little error due to wind.
Both buildings have the lowest airtightness at the winter measurement (January-February). This
coincides with the lowest indoor air relative humidity. The measured airtightness (expressed as air
permeability) as a function of time is shown in Figure 1, and as a function of relative humidity indoor
in Figure 2.
There is an increase in air leakage over time. The increase from the first measurement is 10% for the
building in Sevred and 7.7% for the building in Landvetter. For both buildings, the airtightness is less
when the indoor relative humidity is the lowest. There is no clear correlation with indoor or outdoor
temperature, or outdoor relative humidity. The variations in airtightness can be noticeable for
stakeholders aiming for a certification.
The measurements will continue and will also include a concrete building (from March 2014). In the
two measured wooden buildings shown above, the reason for the increase in permeability could be
that the wooden construction dries when the relative humidity indoor decreases. The main leakages in
the leakiest building (Landvetter) are found around the windows in the bottom floor, around the attic
hatch (see Figure 3) and at the connection between the top and bottom floor (see Figure 4). The
reason is probably a poor connection between window and polyethylene for the windows, a poor
connection between hatch and polyethylene foil plus a leaky hatch (poor seal) for the attic hatch and a
discontinuous polyethylene foil in the exterior wall.
FIGURE 4. Thermographic image of the connection between the upper and lower floor, the brighter
diagonal part in front is the hand rail in the stairs.
5. Conclusions
Airtightness has been measured in two residential wooden buildings. The trend in the airtightness
measurements is that the airtightness is lower when indoor air is drier. The winter measurements (last
reported measurements) have the lowest airtightness of all measurements during the seven months
that the measurements have been ongoing. The decrease in airtightness during the measurement
period is in the order of 8-10% (from July to February).
The simulations show a small change in air flow through the leakages due to high or low air
temperatures (affecting the air density of the air that flow in the leakages), but no change in air flow
due to different relative humidities. The effect of the air density on the measurement equipment (not
leakages) is compensated for when measuring according to EN 13823:2000.
The measurements will continue to cover a whole year, and in March 2014 a concrete building will be
included.
References
Borsboom W. and de Gids W. 2012. Seasonal variation of facade airtightness: field observations and
potential impact in NZEB. Proceedings of the AIVC-TightVent International Workshop, Bryssel,
Belgium, March 2012
Bracke W. Laverge J. Van Den Bossche N. Janssens A. 2013, Durability and measurement uncertainty
of airtightness in extremely airtight dwellings. Proceedings of the AIVC-TightVent International
Conference, Athens, Greece, September 2013
Dickinson J.B, Feustel H.E. 1986. Seasonal variations in effective leakage area, Thermal performance
of the exterior envelops of buildings III, Atlanta, ASHRAE, 144-160.
Kim A.K. Shaw C.Y. 1986. Seasonal Variation in Airtightness of Two Detached Houses, Measured Air
Leakage of Buildings, ASTM STP 904, 1986
Sandberg P-I. Sikander E. Wahlgren P. Larsson B. 2007a. Lufttäthetsfrågorna i byggprocessen- Etapp
B. Tekniska konsekvenser och lönsamhetskalkyler, SP Report 2007:23 (in Swedish only)
Sandberg P-I. Bankvall C. Sikander E. Wahlgren P. Larsson B. 2007b. The effects and cost impact of
poor airtightness- Information for developers and clients. Proceedings of the Thermal Performance
of the Exterior Envelopes of Whole Buildings X, Florida, USA, December 2007
EN 13829:2000 (Thermal performance of buildings - Determination of air permeability of buildings -
Fan pressurization method (ISO 9972:1996, modified))
Yoshino. H. 2012. System for ensuring reliable airtightness level in Japan. Proceedings of the AIVC-
TightVent International Workshop, Bryssel, Belgium, March 2012
Page....................................119-320
1 Introduction
In England, approximately 6.5 million homes are built of solid wall – 31% of the total housing stock,
of which around 60% have been built before 1920. Solid wall dwellings are considered “hard-to-treat-
homes”, since they cannot be upgraded with easy or cost-effective fabric energy efficiency measures
(BRE, 2008). Improving the energy efficiency of these dwellings becomes even harder in conservation
areas, listed buildings, or building with decorative façades where the only feasible solution is internal
wall insulation (IWI); planning permission for external wall insulation (EWI) is often denied.
However, the installation of IWI may affect the interstitial temperature and vapour permeability of the
building envelope leading to moisture accumulation and the reduction of the building durability; high
interstitial relative humidity is ideal for mould growth and timber decay.
This paper describes the experimental test carried out for the validation of a numerical tool for heat,
air and moisture transport. The experiment was designed to help understand the hygrothermal
behaviour of internal wall insulation exposed to transient boundary conditions of relative humidity and
temperature, and to validate a simulation tool commonly used to estimate moisture movement within
building elements and the likelihood of mould growth and fabric decay in buildings. Two walk-in
environmental chambers are utilised for the experiment; wall samples were exposed to climate
conditions set independently in each chamber.
a) b)
FIG 1. Test wall built between the environmental chambers a) view of eight insulation systems exposed
to internal boundary conditions b) view of clay block wall exposed to external boundary conditions
1. 2.
Gypsum plasterboard 12.5 mm Gypsum plasterboard 12.5 mm
Light woodfibre 80 mm Vapour Control Layer (VCL 1) 1 mm
Wood-clay board 20 mm Light woodfibre 100 mm
Lime plaster 5 mm Gypsum plaster 10 mm
Gypsum plaster 10 mm Aerated clay blocks 175 mm
Aerated clay blocks 175 mm
3. 4.
5. 6.
Gypsum plasterboard 12.5 mm Gypsum plasterboard 12.5 mm
Air layer 25 mm Vapour Control Layer (VCL 2) 1 mm
Aluminium foil 0.1 mm Mineral wool 90 mm
Polyisocyanurate (PIR) 60 mm Gypsum plaster 10 mm
Aluminium foil 0.1 mm Aerated clay blocks 175 mm
Gypsum plaster 10 mm
Aerated clay blocks 175 mm
7. 8.
2.1.3 Monitoring
The relative humidity at the interface between the insulation and the substrate wall was monitored
using six temperature and relative humidity sensors applied to the gypsum plaster and subsequently
covered by the insulation system (total of 96 sensors). Thermocouples and capacitive sensors were
used for temperature and relative humidity respectively. Also, room temperature and relative humidity
were monitored in each environmental chamber. Data were collected every 5 minutes, for 43 days, and
averaged every hour. The collected data were then used as the boundary conditions input in the
simulations for the tool validation.
FIG 5. Comparison between experimental data (solid line) and simulation results (dashed line).
Acknowledgments
The authors would like to thank Natural Building Technologies and Technology Strategy Board who
funded the research project.
References
BRE. 2008. A Study of Hard-to-treat Homes Using the English House Condition Survey. Part 1 —
Dwelling and Household Characteristics of Hard-to-treat Homes. Watford, UK, Building Research
Establishment.
Carmeliet J. & Derome D. 2012. Temperature driven inward vapor diffusion under constant and cyclic
loading in small-scale wall assemblies: Part 1" Building and Environment. 48-56.
Fossa M. & Petagna P. 2004. Humidity measurements inside Atlas and CMS: notes on sensor
calibration. [http://proj-jcov.web.cern.ch/proj-
JCOV/EBmeeting_22/040422_Rh_Sensor_calibration.pdf accessed on 13/12/2013]
Vereecken E. & Roels S. 2011. Hygric performance of different interior insulation systems: an
experimental comparison. 12th International conference on Durability of Building Materials and
Components, Porto, Portugal.
SUMMARY:
In recent years there has been an increased use of wood-frame construction in Belgium, which
likewise requires an increased focus on both hygrothermal analysis as well as best practice
guidelines to ensure the longevity of this construction type in a temperate maritime climate. This
experimental research investigates both the airtightness and watertightness of window-wall
interfaces, and the interrelation between both. Three different approaches to secure the
watertightness on the exterior side were evaluated, more specifically the use of a watertight layer,
spray-in-place polyurethane foam, and self-expanding tapes. The sensitivity of these installation
methods was analyzed in respect to the airtightness of the window-wall interface. Based on 3
airtightness performance levels, it was found that the airtightness, and hence the pressure
equalization, has a significant effect on the risk for water ingress. Furthermore, the amount of water
ingress was also quantified under a range of pressure differences under both static and dynamic
boundary conditions. This quantitative approach is also of importance to evaluate the moisture
sensitivity and moisture tolerance of wood-frame constructions. Finally, the paper highlights typical
installation errors, and provides practical guidelines for airtight window installation in wood-frame
construction.
1. Introduction
The market share of wood-frame construction in Belgium is increasing. Studies show that the market
share of wood-frame construction for newly built residential dwellings was 5.6% in 2004, it rose to
11% in 2009, and it is expected to increase up to a market share of 20% by 2020 (BBRI, 2020). In
Belgium, residential buildings are typically constructed with brick cavity walls with about 8 to 20cm
of insulation in the cavity. Consequently, the majority of technical guidelines and sector documents
refer to that building practice, and little to no information is available on wood frame construction.
Furthermore, Viitanen et al. (2010) report the modelled mass loss of pine sapwood for all of Europe,
which shows that the climatic conditions in Belgium are among the worst in respect to wood decay.
In Belgium, the BBRI (2007) reported that 24% of all interventions for buildings relate to moisture,
and water infiltration is the most important cause of building damage (condensation, absorption due to
capillary action and initial moisture content are excluded). Brick cavity walls are the predominant
construction type in Belgium, and this type of wall has the specific capacity of absorbing significant
amounts of water without leading to premature deterioration or failure. Commonly, there is an exterior
brick leaf with solid bricks, a ventilated air cavity, insulation, an extruded and perforated interior
brick wall, and a stucco finish. A limited amount of inadvertent water ingress past the insulation layer
may be stored by the interior brick wall, and due to the high vapour permeability of the construction
3. Experimental setup
A non-operable wooden window (1.55m high by 1.23m wide) was installed in a typical wood frame
wall of 2.28m high and 1.94m wide. This wall comprised an interior airtight sheeting of taped
oriented strand boards, an air-open exterior impregnated fibre board, but for reasons of visual
inspection by means of an endoscope, no insulation was installed between the studs. On the exterior
side, the cladding consists of slender horizontal planks, which was assumed the be rather permeable
for wind driven rain. At the sill and top of the window, cross cavity flashing was installed by means of
an aluminium profile, with end dams at the sill. The wood frame wall was installed perfectly airtight
in a steel test box.
At the bottom side of the rough window opening a water collection tray and drainage system was
installed, to measure the amount of water ingress into the window wall interface. Infiltration water
from the drainage was directed to a hermetically sealed gravimetric trough and weighing system. The
drainage system was designed in such a way that the water that was drained did not occlude the
drainage tube, and in order to avoid artificial pressure gradients caused by the weighing system itself,
the collection through was connected to the window wall interface with a pressure tube. That allowed
for perfect pressure equalization, omitting any inadvertent pressure gradients due to perturbations in
the collection system.
3.1 Test protocol
There is no specific test standard to evaluate the watertightness of joints or interfaces. For windows, a
static test sequence is provided in EN 1027. The European standard EN 12865:2001 describes a test
protocol to determine the resistance of external wall systems to driving rain under pulsating air
pressure. According to this standard, water is sprayed at two locations: 72l/h-m at the top as run-off,
and 90L/h-m² on the whole area to simulate direct Wind driven rain. The pressure difference is
applied in a step-wise approach to determine at which level water entry occurs. Furthermore, the
standard offers two procedures: one for qualitative testing (A - steps of 10 minutes after 20 minutes of
initial wetting), and the other for quantitative testing (B - steps of 60 minutes). Although there is no
standard for window-wall interfaces, the interface is perhaps expected to perform well according to
both test sequences (windows and walls). Consider that the EN 12865 only provides a test sequence,
without information on the required performance levels as a function of exposure to climate
conditions. Consequently, the samples will be tested according to the static test protocol up to 750Pa
(EN 1027) and the cyclic test protocol up to 600Pa (EN 12865).
To assess the effect of the airtightness of the window-wall interface on the risk for water ingress, the
airtightness of every setup was also determined based on EN 12114. Please refer to (Van Den Bossche
et al., 2012) for more information on the airtightness measurements and error analysis.
3.2 Test parameters
Four different parameters were varied during the test sequence: mechanical installation of the window
to the wall, installation method to ensure the watertightness, position of the window in respect to the
wall, and the level of airtightness of the interface.
3.2.1 Mechanical installation
In Belgian brick construction practice, windows are typically fixed to the wall using steel brackets.
Direct fixation with screws through the window frame itself – typical for German window installation
in ETICS – is not advocated in Belgium, and flanged windows are not available. The brackets allow
for a practical and flexible installation, but conversely hamper the easy installation of airtightness
membranes at the window perimeter. The alternative is a plywood casing around the window frame,
which in turn replaces the brackets. The plywood provides a solid ground for the finish of the window
reveal, and allows to secure the airtightness easily.
3.2.2 Position of the window in respect to the wall
The position of the window frame might depend on the alignment with the thermal insulation,
esthetical or mechanical boundary conditions, but by consequence will also affect the exposure to
wind driven rain. Four positions are considered (see figure 2): behind the impregnated fibre board (a),
onto the back of the fibre board (b), aligned with the fibre board (c), or projected outwards (d). Mind
that the critical pathways for water ingress also shift from the bottom to the top as the position of the
window shifts to the exterior. The highest risk for infiltration is located at the corners, where –
depending on the position of the window – complex folding may be required to ensure the continuity
of the foil. When the window is aligned with the fibre board, the foil can easily be installed with a low
risk for installation errors.
Figure 3. Installation with self-expanding tape (left) and foil (middle and right)
4. Results
4.1 Mechanical installation
The tests indicate that the type of mechanical installation did not affect the watertightness
performance when self-expanding sealant tape or watertightness foil was applied. This was most
evident for installations with a good airtightness, whereas installations with significant air leakages
tend to show small discrepancies. Nevertheless, differences in water ingress could be attributed to
executing performance. In general, installations using a plywood frame around the window frame
proved to be less susceptible to errors in the airtightness: the use of metal brackets typically
introduces more difficulties to ensure the airtightness, but for equal airtightness similar performances
were found. When spray in place polyurethane foam was applied, a difference in failure behaviour
was noticed. The foam itself did not show any leakage, but local deficiencies at obstacles may lead to
water ingress. When installing metal brackets, it is important to make sure that the bracket itself does
to protrude to the exterior, where it is exposed to runoff. Even though the foam was installed
carefully, it does not provide a continuous barrier. For the installation with plywood frames, small
wood blocks are used as spacer to position the window frame in respect to the rough window opening.
At these blocks, the continuity of the foam is interrupted, which also gave rise to water ingress.
4.2 Position of the window in respect to the wall
Only for the case where the watertightness was ensured by a water resistive barrier as flashing, the
effect of the position of the window was evaluated. As the window is located more inwards, the
flashing at the bottom is more exposed, and inadequate drainage may lead to local accumulation,
hydrostatic pressure and infiltration. Likewise, as the window is located more towards the exterior,
the flashing at the top is similarly more exposed. However, the results show that more problems can
be expected at the sill, as the window itself causes a high runoff rate at the sill, whereas at the top
most of the water is diverted by the cross-cavity flashing.
4.3 Installation method and airtightness
The installation with sealant tape did not perform well. At the sill, the installation of the sealant tape
apparently led to accumulation of water, which caused a hydrostatic pressure. Next to that, additional
repetition tests showed that the specific sealant tape used in this study displayed shrinkage and
swelling behaviour under wetting and drying cycles. With each cycle the sealant tape shifted slightly
downwards due to gravity. Consequently, horizontal joints should be avoided with this type of tape,
and otherwise these must be shielded from rain and runoff as much as possible. The compression of
the tape should be adequate to prevent any sinkage of the tape. For a good airtightness the amount of
water infiltrating to the interior was rather limited, whereas the water ingress increased as the
airtightness decreased. The results reported in table 1 are those for window location b. The
installation with a self-adhering flashing membrane to the window frame proved to be very
susceptible to execution errors. The adhesion was found inadequate, and the joint between window
and flashing was henceforth sealed by means of caulking (reported in the table below). Without
caulking, water ingress rates of 14,7 up to 34,6 gram/min were measured under static conditions, and
also under dynamic conditions the infiltration rates were higher. Table 1 again indicates that a good
5. Conclusions
The watertightness of window-wall interfaces in wood-frame constructions was evaluated in lab
conditions, under static and cyclic wind pressure loadings. Two types of mechanical fixing to the
interior wall were evaluated: steel brackets and plywood casing around the window. The steel
brackets require more attention when installing the airtightness foil, but no significant effect was
found on the level of watertightness. When polyurethane foam is used, one should ensure a
continuous layer of foam, avoiding any crossing brackets or other obstructions. Tests on the position
of the window were not reported in detail, but for all positions the sill was most prone to water
ingress probably due to higher local runoff rates. The use of sealant tape led to high infiltration rates
References
BBRI, 2010. Research on new building concepts, and the impact on natural resources. (in dutch).
Belgian Building Research Institute, VITO, KULeuven.
Viitanen H., Toratti T., Makkonen L., Peuhkuri R., Ojanen T., Ruokolainen L., Räjsänen J., 2010.
Towards modelling of decay of wooden materials. European Journal of Wood and Wood Products
68, 3: 303-313.
BBRI, 2007. Annual report. www.wtcb.be
Ruest K., Wall moisture problems in Alberta dwellings. CMHC Research Highlights (2000) CMHC,
Ottowa, Canada.
RDH, 2000a Water penetration resistance of windows – study of manufacturing, building design,
installation and maintenance factors, RDH Building Engineering Limited, Vancouver, 2002.
RDH, 2000b Water penetration resistance of windows –study of codes, standards, testing and
certification. RDH Building Engineering Limited, Vancouver, 2002.
Rousseau M., An Overview of the Survey of Building Envelope Failures in the Coastal Climate of
British Columbia, performed by Morrisson-Hershfield Limited for CMHC (1996). Journal of
Building Physics (1999) 22: 364-367.
Lstiburek J.W., Rainwater management performance of newly constructed residential building
enclosures during august and September 2004. (2005) Building Science Corporation, MA.
Nelson C., Norris R.E., Mock-Up Water Test Results of Sample Flashing Systems for Storefront
Windows in Stucco Walls. Interface (2010) RCI inc., Raleigh, NC.
Knowles S., WANZ Window Installation System. Window Association of New Zealand. Wellington,
NZ, 2002.
House of Representatives, Weathertightness of Buildings in New Zealand. Forty-Seventh Parliamant,
March 2003.
Lacasse M.A., Rousseau M.Z., Cornick S.M., Manning M.M., Nicholls M., Nunes S., Performance
evaluation of Wall-Window Interface Details Phase 1 – Watertightness, Air Leakage and Rainwater
Management of CMHC Specified Assemblies. Report B-1229.1. National Research Council, 2008.
Lacasse M.A., Rousseau M.Z., Cornick S.M., Manning M.M., Ganapathy G., Nicholls M., Williams
M.F., Assessing the effectiveness of window-wall interface details to manage rainwater – selected
results from US window installation practice. Journal of ASTM International (2009) 6, 8: 1-35.
Salzano C.T., Masters F.J., Katsaros J.D., Water penetration resistance of residential window
installation options for hurricane-prone areas. Building and Environment (2010) 45: 1373-1388.
Birkeland Ø. The design of multi-layer walls. Report No 24, Norges Byggforskningsinstitutt. Oslo,
1957.
Edgar J., Performance of Source Drainage Exterior Insulation Finish System at the window/wall
junction. Journal of Building Physics (1999) 23: 57-77.
Finn D., Rain penetration control of the window-wall joint. Institute for Research in Construction,
NRC. NRCC-38747, 1991.
Van Den Bossche N., Huyghe W., Moens J., Janssens A., Depaepe M., 2012. Airtigthness of the
window-wall interface in brick cavity walls. Energy and Buildings 45: 32-42.
SUMMARY:
Evaluating the thermal performance of novel advanced building enclosures is challenging, and there is
a need for better and more cost effective methods for doing so. Laboratory testing of components
provides accurate results under special conditions, however, these may fail to represent complex
dynamic interactions between system components. Field testing can offer more realistic performance
evaluation, however, studies involving occupied buildings suffer from many confounding factors, such
as occupant interactions and unknown building material properties that make it difficult to generalize
results. Whole building energy modeling is the ultimate goal for evaluating and designing new
building systems, however, without proper validation, the models cannot be trusted. In this paper we
present a flexible approach for cost effective evaluation field testing of building advanced enclosures.
1. Introduction
New building components or systems must be evaluated and tested before going to market. From an
energy and durability standpoint, it is desirable to measure performance on several dimensions,
including thermal performance (heat gain or loss through components), occupant comfort (temperature
and humidity in occupied zones), moisture performance (moisture accumulation and material
durability), air tightness (and its impact on energy consumption), and energy performance (HVAC
energy reduction, peak load shifting). Various testing approaches have been developed and tried in the
last century. This paper is not meant to be a comprehensive review of such methods; however, it does
illustrate some different philosophies and why the proposed approach may be appropriate for specific
testing goals.
We present a flexible, nearly-calorimetric testing platform that combines the advantages of outdoor
field testing with the controllability of a laboratory experiment. The platform is comprised of multiple
small testing structures, or huts. One hut is built to standard construction practices (e.g., typical
residential construction) and acts as a control or baseline case. The other hut(s) are constructed
identically, with the exception of the novel building materials or components. All huts are conditioned
with identical HVAC systems and control strategies, and HVAC energy consumption is metered
during the experiment. Air tightness of each hut is measured prior to each experiment. Detailed
temperature and heat flux measurements are made across assemblies of each hut. A weather station
records outdoor conditions. Resulting data may be used to support the validation of energy models for
advanced systems or to rapidly evaluate prototype building systems or components.
This paper presents the results of testing measurements made during the fall in three huts used to
evaluate an advanced insulation material. Here we focus on description of the platform, construction
and experimental technique, and the measurements obtained during testing. Through a six month
2.1.3 Examples
Two example of a single-structure test is the test room of (Athienitis et al. 1997), which compared
phase change material board temperature with that of ordinary wall board and developed mathematical
models of material performance, and the test room of (Shilei et al. 2007), which tested phase change
wallboard energy performance by measuring electric heating energy and temperature.
A more sophisticated, but slightly different type of single-structure test is a calorimetric testing
chamber that exposes one side of a test element to the exterior environment and controls the interior
side precisely. All remaining sides of the test chamber are heavily insulated and instrumented with
heat flux sensors to obtain a careful heat balance of the space. This kind of apparatus has been used
and tested extensively throughout research laboratories in Europe under various names DYNASTEE,
PASLINK, and PASYS (e.g., Baker 2008). This approach is useful for obtaining assembly-specific
parameters that can be used in simulation, but it may not capture the full system-level interactions
between multi-zone spaces.
Examples of multi-structure testing include (Szymocha et al. 2005), which tested an advanced thermal
solar system with heat storage and compared surface temperature, surface irradiance, and space
heating with three test structures; (Medina et al 2007) tested a phase change material structural
insulated panel using arrays of thermistors and heat flux sensors with two test structures; (Khudhair
and Farid 2007) tested phase change material load shifting in two test structures
In some (intentionally unnamed) cases the researchers have mistakenly relied on insufficient data to
produce conclusions about energy savings. For instance, measuring only interior temperature and
cumulative heating or cooling energy consumption is normally not sufficient to convincingly show
that (1) savings do exist, and (2) savings originate from the element in question. Accurate and
sufficient measurement is the cornerstone of any evaluation. Since, in many cases, the effect of
improving a building component can be small (say 5% improvement), it is essential to use sensitive
equipment and sample frequently enough to measure these effects reliably.
2.2.1 Experiment
In this experiment, we set out to measure the thermal and energy performance of two proprietary
insulation systems against an uninsulated baseline structure. Each hut was conditioned to an identical
setpoint temperature, and allowed to run for several months during the summer and fall of 2013 in the
desert climate of Albuquerque, NM. We monitored temperature and heat flux through the assemblies,
energy consumption of the cooling equipment, and climate conditions. We used the results to analyze
energy reduction potential, peak load shifts and reductions, and zone temperature behavior.
The insulation systems were identical other than the presence of a reflective surface foil. Ordinarily,
the baseline condition would include some standard insulation materials; however, in this case, the
uninsulated case was more representative of the particular end-use application.
Figure 1: Three identical test huts, pictured from the North-east corner of the lot.
Figure 2: Left to right: uninsulated hut interior without ceiling, non-reflective insulation installed on
the roof interior of Huts 2 and 3, reflective insulation installed on the walls of Hut 3 with ceiling
installed.
The test insulation was then nailed to the wall and roof deck studs of Huts 2 and 3, while Hut 1
remained uninsulated. Afterwards, ceilings were installed in each hut to separate the interior
conditioned zone from the unconditioned attic space. Hut 2 received the insulation packaged in a non-
reflective facing surface, while Hut 3 received insulation packaged with a reflective foil surface on
both sides of the batts. The purpose of the reflective facing on the Hut 3 insulation was to reduce
radiative heat transfer across the stud cavities created by the insulation. After the insulation was
installed in Huts 2 and 3, OSB ceilings were installed on all three huts to create a separate,
unconditioned attic space.
3.2 Temperature
We measured temperature across all key wall assembly components, as well as the interior zone
temperatures of all huts. By measuring temperature across lightweight construction materials of
known properties, it becomes possible to estimate the heat flux through those materials. This provides
a backup method for calculating heat flux through surfaces where sensors fail or where cost prohibits
their use.
The temperature measurements were obtained using NIST Calibrated type T thermocouples with
specified accuracy of ±0.5°C or 0.4%. Further laboratory tests were conducted to determine the
accuracy by measuring thermocouples in conjunction with a calibrated thermistor in a metrology well
from a range of 0°C to 60°C. To more accurately measure temperatures with our data logger, an
external reference junction box was used. The isolated isothermal external reference junction box
(zone box) is used instead of the data logger’s internal reference junction, minimizing fluctuations that
could be caused by rapid changes in the outdoor temperature.
Figure 3: Installing temperature and heat flux sensors. Left to right: temperature sensor on exterior
wall, heat flux sensor on interior side of OSB and thermocouple array through assembly, and
installing heat flux sensor on reflective insulation with reflective tape.
3.4 Electricity
We monitored electricity consumption of each hut in one-minute resolution to approximate the cooling
energy required by each hut. In this experiment, the instrumentation power was very low (<5W), and
there were no additional internal building loads.
4. Sample Results
The following results from a sample week illustrate the sample data results and relative differences in
performance among the three test huts.
Figure 6: Wall heat flux from different surfaces: EW=east wall, WW=west wall, etc.
Figure 7: Zone Load (calculated from surface heat flux measurements) and AC Power draw.
5. Conclusion
We have described and demonstrated a flexible platform for accurately testing the thermal
performance of advanced building systems through the use of small outdoor test structures. By
carefully measuring key building parameters, we have found that is possible to study and characterize
the performance of building systems in scaled-down versions of realistic building constructions.
In the materials we tested we found that insulated structures consumed about 20% of the energy
needed by the uninsulated structure. Future work will use the detailed thermal measurements to
develop calibrated energy models for the test huts, and energy savings will be simulated with buildings
that have similar characteristics.
6. Acknowledgements
This material is based upon work supported by The Small Business Innovative Research (SBIR)
program and the Engineer Research and Development Center – Construction Engineering Research
Laboratory (ERDC-CERL) under Contract No. W9132T-12-C-0012. Special thanks to Bob Andrews
and Alliston Watts for their many important contributions to this project.
References
ASTM. 2007. C1046-95(2007). “Standard Practice for In-Situ Measurement of Heat Flux and
Temperature on Building Envelope Components.” American Society for Testing and Materials.
ASTM. 2012. C1130-07(2012). “Standard Practice for Calibrating Thin Heat Flux Transducers.”
American Society for Testing and Materials.
Athienitis, A.K., C. Liu, D. Hawes, D. Banu and D. Feldman. 1997. “Investigation of the Thermal
Performance of a Passive Solar Test-room with Wall Latent Heat Storage.” Building and
Environment, Vol. 32, No.5, pp. 405-410.
Baker, P.H. and H.A.L. van Dijk. 2008. “PASLINK and dynamic outdoor testing of building
components.” Building and Environment. 43(2008)143-51.
Khudhair, A.M. and M.M. Farid. 2007. “Use of Phase Change Materials for Thermal Comfort and
Electrical Energy Peak Load Shifting: Experimental Investigations. In: GOSWAMI, D. Y. &
ZHAO, Y., eds. Solar World Congress 2007, September 18-21. Beijing, China. Solar Energy and
Human Settlement.
Medina, M.A., J.B. King, and M. Zhang. 2008. “On the Heat Transfer Rate Reduction of Structural
Insulated Panels Outfitted with Phase-change Materials” Energy. Vol. 33, Issue 4:667-678.
Shilei Lv., F. Guohui, Z. Neng, and D. Li. 2007. “Experimental Study and Evaluation of Latent Heat
Storage in Phase Change Materials Wallboards.” Energy and Buildings. Vol. 39 (2007), pp 1088-
91.
Szymocha, K. 2005. “Advanced Thermal Solar System with Heat Storage for Residential House Space
Heating.” Proceedings of SESCI 2005 Conference British Columbia Institute of Technology,
Burnaby, British Columbia, Canada, August 20-24.
Summary
The modeling of hygrothermal performance using computer programs has increasingly become a
common task in building envelope design. Typically, the deterministic approach is used in
hygrothermal modeling to assess the Heat, Air and Moisture responses and moisture risks of various
designs. However, in reality factors influencing the hygrothermal responses are stochastic in nature
such as the variability of material properties, boundary conditions, construction quality, and
environmental loads, especially the wind-driven rain and runoff loads. To improve the reliability of
hygrothermal modeling, stochastic modeling approach can be used to establish the correlation between
input uncertainties with the output uncertainty and identify the most influential parameters. However,
the importance of these parameters is not constant and varies with wall configurations, boundary
conditions, and moisture loads. To investigate the importance of various design parameters, stochastic
modeling approach is applied to a brick veneer wood-frame wall assembly. Four scenarios with
different cavity ventilation rates and rain leakage deposition positions are tested. For each scenario,
Monte Carlo simulation and multiple regression analysis are performed to rank the influential
parameters. The results indicate that the importance of material properties is largely influenced by the
cavity ventilation and the deposition location of rain leakage.
1. Introduction
Hygrothermal performance modeling has increasingly become a common tool used in building
envelope design in recent years, especially for innovative and energy efficient buildings. Typically, the
deterministic approach is used in hygrothermal modeling to assess the Heat, Air and Moisture
responses and moisture risks of various designs. However, in reality factors influencing the
hygrothermal responses are stochastic in nature such as the variability of material properties, boundary
conditions, construction quality, and environmental loads, especially the wind-driven rain and runoff
loads. The stochastic approach has been applied in several studies in hygrothermal performance
simulations to investigate the influence of material properties (Salonvaara et al. 2001; Holm and
Kuenzel, 2002; Defraeye et al. 2012). Zhao et. al. (2011) also looked into the influence of boundary
conditions and ranked the importance of each parameter using multiple regression analysis. The
influence of environmental conditions including rain loads and internal moisture gains on the
hygrothermal performance of wall assemblies was also investigated by parametric studies and scenario
analysis (Cornick, et al. 2009 and Pallin, 2013).
Previous studies mainly focused on the influence of the stochastic nature of material properties. As
indicated in the study by Zhao et al. (2011), the importance of each influential parameter was not
constant and varied with the weather conditions. To provide a reliable hygrothermal performance
assessment of building envelopes, especially for those highly insulated or with innovative materials, a
stochastic approach taking into account the combined influence of material properties, boundary
conditions, and moisture loads is necessary. The objective of this paper is to develop a methodology to
assess the moisture damage risks of wood-frame construction using a stochastic approach. A brick
Air Cavity
Plywood
Glass Fiber
Polyethylene
Gypsum Board
According to Zhao et al. (2011), the thermal properties seem having insignificant influence when MC
is taken as the performance indicator. Therefore, only the hygric properties are considered as
stochastic variables. To examine the effect of cavity ventilation and rain leakage, the material
properties of red brick, SBP sheathing membrane and plywood are selected for the stochastic analysis.
Both the deterministic and stochastic parameters are listed in Table 1. The material properties are
taken from the database in WUFI Pro 5.0. The properties of stochastic variables, mean value with a
standard deviation in the bracket, are determined based on the measurement results reported by
Kumaran (2002). All the stochastic variables are assumed to follow a normal distribution.
Figure 2 Stochastic Results of Scenario 1 (no rain leakage with cavity ventilation)
Figure 3 Stochastic Results of Scenario 2 (1% rain leakage on SBP membrane with cavity ventilation)
Figure 4 Stochastic Result of Scenario 3 (1% rain leakage on Plywood with cavity ventilation)
Figure 5 Stochastic Results of Scenario 4 (1% rain leakage on SBP membrane without cavity
ventilation)
PCC in Feb. 5th
1
0.5
PCC
0
-0.5
-1
PCC NO Rain Leakage with Cavity Ventilation PCC 1% Rain Leakage on SBP with Cavity Ventilation
PCC 1% Rain Leakage on Plywood with Cavity Ventilation PCC 1% Rain Leakage on SBP without Cavity Ventilation
u_brick-vapour diffusion resistance factor of brick
mr_brick-reference moisture content of brick
mf_brick-free water saturation of brick
A_brick-water absorption coefficient of brick
u_SBP-vapour diffusion resistance factor of Spun Bound Polyolefin
u_plywood-vapour diffusion resistance factor of plywood
mr_plywood-reference moisture content of plywood
mf_plywood-free water saturation of plywood
A_plywood-water absorption coefficient of plywood
FD-rain deposition factor
VR-air change rate
0
-0.5
-1
PCC No Rain Leakage with Cavity Ventilation PCC 1% Rain Leakage on SBP with Cavity Ventilation
PCC 1% Rain Leakage on Plywood with Cavity Ventilation PCC 1% Rain Leakage on SBP without Cavity Ventilation
u_brick-vapour diffusion resistance factor of brick
mr_brick-reference moisture content of brick
mf_brick-free water saturation of brick
A_brick-water absorption coefficient of brick
u_SBP-vapour diffusion resistance factor of Spun Bound Polyolefin
u_plywood-vapour diffusion resistance factor of plywood
mr_plywood-reference moisture content of plywood
mf_plywood-free water saturation of plywood
A_plywood-water absorption coefficient of plywood
FD-rain deposition factor
VR-air change rate
6. Conclusion
Stochastic modeling can be used in the sensitivity analysis of hygrothermal performance of building
envelopes. The sensitivity of the performance indicator in relation to input parameters often varies
with wall configurations, environmental conditions, and moisture loads. To examine the combined
effect of material properties and environmental conditions, both stochastic analysis and scenario
analysis are performed using hygrothermal simulations. A typical wood framed building envelope
with cladding ventilation is used for the stochastic modeling. The influence of ventilation rates and
deposition rain factor is investigated by performing stochastic simulations under four different
scenarios. The results indicate that the importance of the material properties is largely influenced by
the cavity ventilation and the deposition location of rain leakage. Further research should be focused
on developing wall configurations that are less sensitive to the variation of local environmental
conditions and rain loads.
References:
ASHREAE. 2009. ASHRAE Standard 160-2009 Criteria for Moisture-Control Design Analysis in Building.
American Society of Heating, Refrigerating and Air-Conditioning Engineers, Inc.
Cornick, S.M., Dalgliesh, W.A. & Maref, W. 2009. Sensitivity of Hygrothermal Analysis to Uncertainty in Rain
Data. Journal of ASTM International, 6 (4): 1-17.
Defraeye, T., Blocken, B. & Carmeliet, J. 2013. Influence of uncertainty in heat-moisture transport properties on
convective drying of porous materials by numerical modelling. Chemical Engineering Research and Design ,
91(1): 36-42.
Holm, A. & Kunzel, H.M. 2002. Practical application of an uncertainty approach for hygrothermal building
simulation-drying of an AAC flat roof. Building and Environment. 37(8-9): 883-889.
Kumaran, K., Lackey, J., Normandin, N., van Reenen, D. & Tariku, F. 2002. Summary Report from Task3 of
MEWS Project. Institute for Research in Construction, National Research Council Canada. Ottawa, Canada,
(NRCC-45369), pp.1-68.
Lomas. K.L., & Eppel, H. 1992. Sensitivity Analysis Techniques for Building Thermal Simulation Programs.
Energy and Building, 19(1): 21-44.
Pallin, S. 2013. Risk Assessment of Hygrothermal Performance-Building Envelope Retrofit. Ph.D Thesis.
Chalmers University of Technology. Goteborg, Sweden.
Salonvaara, M., Karagiozis, A., & Holm, A. 2002. Stochastic Building Envelope Modelling-The Influence of
Material Properties. Proceedings for Performance of Exterior Envelopes of Whole Building Ⅷ: Integration of
Building Envelope, Clearwater Beach, Florida.
Simpson, Y. 2010. Field evaluation of ventilation wetting and drying of rainscreen walls in coastal British
Columbia. MASc. Thesis. Concordia University. Montreal, Canada.
Zhao, J., Plagge, R., Nicolai, A. Grunewald, J. & Zhang, J. 2011. Stochastic Study of Hygrothermal Performance
of a Wall Assembly- The influence of Material Properties and Boundary Coefficients. HVAC&R Research,
17(4): 591-601.
KEYWORDS:
CDM, Cities, Thermal Energy Simulation, CO2-Emissions, Energy Demand
SUMMARY:
Increasing urbanisation and climate change are one of the greatest challenges in the 21st Century.
Many regions already face different negative impacts on cities such as growing local pollution, solid
and liquid waste, traffic congestion and noise. A growing consumption of fossil fuels in cities leads to
increasing CO2-emissions accelerating climate change. All over the world and also in China, cities
are trying to reduce their carbon footprint. In that context, the concept of so called “low carbon
cities” is promoted. Currently, 79 cities in China are considering low carbon concepts. To build a low
carbon city is not an easy task as there are various barriers to be removed, such as financial,
traditional and institutional barriers. Building belong to the most traditional and slow-changing
innovative areas in industry. Further solutions are different for hot and warm climates than for cold
climates. Due to a very fragmented and divided structure in smaller enterprises and in different areas, a
holistic approach is hardly to get.
1. Introduction1
The Chinese government decided to reduce the domestic emission reduction goals in the 12th Five
Year Plan (2011-2015) according to the goals of international emission reduction. The reduction of
carbon intensity per unit of GDP (Gross Domestic Product) in China is supposed to be reduced by 40-
45% in 2020 against the intensity of 2005 levels.
The implementation of CDM2 in the United Nations is in the beginning; therefore China has the
chance to implement these systems as the first and the most important Nation worldwide in the
building sector for LOW Carbon Cities. Germany, as one of the leaders in Low Carbon projects, has
the chance, in collaboration with China, to increase the leadership in this working field and to become
the master in the implementation of CDM projects in buildings worldwide. An excellent possibility to
reach this goal is therefore the implementation on financial incentives by pCDM (Programmatic
1
Corresponding author email: [email protected]
1
This report is part of the final report of the study for the German Federal Ministry for the Environment, Nature
Conservation and Nuclear Safety, Chapter 1602 Title 896 05, 2009 as granted by the notification dated
December 3rd, 2009
2
The Clean Development Mechanism (CDM), defined in Article 12 of the Protocol, allows a country with an
emission-reduction or emission-limitation commitment under the Kyoto Protocol (Annex B Party) to implement
an emission-reduction project in developing countries. Such projects can earn saleable certified emission
reduction (CER) credits, each equivalent to one tonne of CO2, which can be counted towards meeting Kyoto
targets. [http://unfccc.int/2860.php]
4.1 TRNSYS
TRNSYS is a complete and extensible simulation environment for the transient simulation of thermal
systems including multi-zone buildings.
4.2 LEC
LEC (Low Energy Certificate) is the result of two years of research to develop a planning tool to
evaluate the energy performance of buildings. Thanks to the evaluation programme it is possible to
evaluate nearly all building types and parts of a building with regard to their energetic quality
separately after the heating period and after the cooling period. The examination of the buildings with
regard to the regional climate conditions is based on a pure physical basis. In this context it is
important to mention that the calculations are exclusively based on results of analysis of the thermal
equations.
The basis for the evaluation of the heating period is a comparison with similar buildings (Reference–
building) which is in accordance with a building according the standard of the 80ies. As far as the
cooling periods are concerned, the evaluation is based on comparisons with an “optimal” envelope of a
building. The evaluation of the energy for cooling and heating are re-evaluated, classified and shown
with regard to certain criteria. The result of the energetic verification is presented with a simple star
system. An increase of stars clearly shows the energetic quality of the building, which means that the
user can immediately recognize the energetic quality of the building thanks to a simple illustration.
4.3 BEED
Building Energy Efficiency Design, Calculation and Economy Analysis Software was developed by
the Centre of Science & Technology of Construction Minister of Construction; P.R. China cooperates
with Beijing E-house building science and technology development Co. Ltd. to finish the task of the
World Bank "economy analysis about building energy efficiency measures and stimulation calculation
module in BEED". The result has passed the MOC evaluation and the level is in the top in China.
BEED was designed according as "Thermal design code for civil building" and "Energy conservation
design standard for new heating residential building.
Table 2. – CO2 emissions Jimei, different building standards according to the building standard LEC
(t CO2/a)
Building type Floor area (m²) Baseline 3* 4* 5*
Single family houses 125969 4624 3081 2977 2664
Multi storage houses 1142439 47464 29545 27139 25180
Office buildings 946974 68184 61800 54332 46194
Multi storage houses /
Office buildings (mix) 1247028 64288 43795 39312 36264
Shopping Centers 647986 52158 43418 40041 37955
Hotels 142016 16831 8503 7854 7083
Congress / Fair 11192 1078 834 772 733
Schools 68966 3493 1476 1249 1187
Other buildings 535384 40473 35351 31715 28292
TOTAL 4867954 298594 227805 205390 191837
Total (in % of baseline) 100% 76% 68% 64%
CO2-emissions for heating, cooling and dehumidification MIN values.
Standards according to: Baseline = applicable Chinese standard
3* acc. Chinese Green Building Standard; 4* German ENEV 2001; 5* German ENEV 2009
Even considering regional characteristics, it is shown that there is a possibility to determine the CO2
potential for new „Low Carbon Cities” of China. The results will deliver the base to make decisions by
implementing new CDM-Baselines to limit the energy demand in a special „Low Carbon Standard“, a
standard quite above the actual energy standard of China.
References
PROPOSED NEW BASELINE AND MONITORING METHODOLOGIES - (CDM-NM)
CDM project activity categories – III.AE. Energy efficiency and renewable energy measures in new
residential buildings
ASHRAE5 Guideline 14-2002, Measurement of Energy and Demand Savings, Whole Building
Calibrated Simulation Performance Path6
GB 50189-2005 Chinese Guideline - public buildings
GB/T 50378-2006 Chinese Guideline - Green Building Standard
GB 50176-93 Chinese Guideline - residential buildings
Energieeinsparverordnung–EnEV2009 and DIN V 18599–Energetische Bewertung von Gebäuden
(07/2005)
TRNSYS - Transient Energy System Simulation Tool; University of Wisconsin, Madison
LEC- Low Energy Certificate - http://www.lowenergycertificate.com
KEYWORDS: Energetic retrofitting, façade integration, central ventilation system, air ducting
SUMMARY:
For successful transition of the energy system being sustainable, the energetic retrofitting of the
building stock is a crucial factor. About 40% of the current final energy consumption in Germany is
caused by the building sector for heating purposes. Energetic retrofitting can reduce the heat demand
of buildings considerably but extensive refurbishments often fail due to high costs, leading to only
individual measures on the façade being carried out. New retrofitting concepts are required which
are both energy efficient and cost-effective to increase the rate of retrofitting significantly.
A promising solution, developed by the Fraunhofer Institute for Building Physics, provides the
implementation of the air duct system in the building envelope. Hereby, the air ducts for the
installation of a central residential ventilation system are integrated directly in the mounted ETICS
plate. Such a system has already been successfully tested in one and two family houses. However, the
challenge of implementing the system in a multi-storey building is much larger due to structural
dependencies on the façade, fire protection regulations and other limiting factors. These constraints
have been identified to develop appropriate solutions for the planning of the air ducting on the
façade.
1. Introduction
1.1 Importance of the building sector in our future energy system
In Germany, the building sector causes 40 % of the current final energy consumption. The main part
of that energy is used in buildings for space heating and domestic hot water obtained from fossil
energy sources. The energy concept of the German Federal Government (BMWi/BMU, 2010) assigns
the energetic retrofitting of the existing building stock a significant role for the upcoming energy
transition. Until the year 2050, the entire energy supply in Germany which is today mainly supplied
by fossil and nuclear energy sources should be replaced by renewable energies. Furthermore these
enormous efforts to rehabilitate the existing building stock are crucial to achieve the climate change
objectives of the European Union.
Currently, the legislature in European countries has already reached a very high level of energy
efficiency for new building, like the directive on the energy performance of buildings (EPBD) in 2003
or national ordinances like the Energy Conservation Act (EnEV) in Germany. New energy standards
for buildings have been developed like the passive house concept, a building without a conventional
heating system due to a high level of thermal insulation and optimized passive solar energy gains or
the plus energy house, which produces throughout the year due to a large photovoltaic system more
energy than it consumes. These concepts are proven, technically mature and well-established in
building practice (Erhorn, 2012). Despite these technical processing, throughout all European
countries the existing building stock contributes about 80 % of the overall energy consumption. The
FIG 1: ETICS panel with integrated air duct (FAW) FIG 2: Air ducting on the façade
For every building class, a typical building was selected and examined regarding possible solution for
air ducting on the façade with the FAW system. Class 1 building is a detached residential house with
two dwellings and class 2 building is a terraced house with one dwelling. Furthermore, the example
building for class 3 is a three-storey residential building which you can find mostly in larger cities.
Also class 4 building is a typical multi-storey building in urban area with four floors.
FIG 5: Connection between the inner distribution channels and the air ducting integrated in the
thermal insulation in the attic through the jamb wall (left) and in the basement through the cellar wall
(right)
FIG 6: Front view of class 3 building with lateral inner and outer walls in dark lines (left) and the
associated floor plan (right)
Buildings from the reconstruction years in Germany rarely have projections and recesses on the
façade but on the rear side often recessed balconies or loggias can be found. The rooms located
behind the balconies cannot be connected to a vertical air duct without using deflections on the
façade. These constraints are shown in figure 6 (right side) to illustrate that challenge. Therefore, the
air ducting in the area of the balcony can be realized for example by breaking through the ceiling
selective. However, the sealing of the balcony should be maintained and the process of retrofitting
should be preferably with little effort. Another possibility is to provide the lateral surfaces for air
ducting, when they are not needed otherwise.
FIG 7: Solution for the air ducting on the façade with limited space between the windows: connection
of two rooms with one duct (left) and with air ducting through the attic and the basement (right)
For large buildings (building class 4) additional fire protection measures are required e.g. fire barriers.
The concept of separation of floors, described above, can be applied (figure 8). The fire barrier then
splits two floors and the rooms are either supplied via air ducts from the attic or the basement. If only
one ventilation device is desired, the fire barrier will be crossed by air ducts. This is technically
possible with nonflammable ducts made of metal and laminated with an intumescent coating. For this
purpose, further investigations are necessary to evaluate the performance in case of fire and smoke.
FIG 8: Fire barrier and separated supply air ducting through the attic and the basement, exhaust air
ducts installed subsequently in existing chimneys
6. Acknowledgements
The author wants to thank the German Federal Environmental Foundation (DBU) for the PhD
scholarship and the associated promoting of the research project. Furthermore, the authors are grateful
to ABG Nova for the providing of the demonstration building that is to be renovated within the EU
Project (EU RetroKit, 2012).
References
Erhorn H. 2012. Das Effizienzhaus Plus. Fachtagung „Die neue Plus-Energie-Welt“. Berlin
http://www.forschungsinitiative.de/PDF/Plusenergie_IBP.pdf
European Union. 2002. Directive on the energy performance of buildings (EPBD). http://eur-
lex.europa.eu/LexUriServ/LexUriServ.do?uri=OJ:L:2003:001:0065:0065:EN:PDF
European Union. 2012. RetroKit project. http://www.retrokitproject.eu/web/guest
Federal Government of Germany. Energiekonzept für eine umweltschonende, zuverlässige und
bezahlbare Energieversorgung. Press information Office of the Federal Government. Berlin.
Germany
Fraunhofer Institute for Building Physics. 2010. Innovativ saniert! – mit integrierten Lüftungskanälen
in der Außendämmung, press release. Kassel. Germany
Hauser G. &Kaiser J. 2013. Dämmstoffintegrierte Kanäle für zentrale Lüftungsanlagen mit
Wärmerückgewinnung. Article Bauphysik 6/2013. Ernst & Sohn. Berlin.
IS-Argebau. 2002. Musterbauordnung (MBO). Model Building Code. Berlin. Germany
Stiegel H. & Krause M. 2012. Minimalinvasives Sanierungssystem mit vorgefertigtem,
multifunktionalem WDVS-Fassadenmodul. GI Gesundheits-Ingenieur. 133(2012) Heft 6.pp. 290-
302
WDVSysteme. Fachverband Wärmedämmverbundsysteme e.V. 2009. Technische Systeminfo –
„Brandschutz“. specialized information. Baden-Baden. Gemany
Ziegler M. 2012. Thermische und strömungstechnische Simulation einer Fassadendämmung mit
integrierter Luftführung für die Bestandssanierung. conference BauSim 2012. Berlin. Germany
Ziegler M. & Krause M. 2013. Innovative façade refurbishment with integrated air ducting for the
existing building stock. conference SB13. Munich. Germany
SUMMARY:
The discussion on sustainability is now over a decade long. How this concept relates to roof systems is
unique in that the variety of potential solutions is almost endless. If one accepts the conclusion of the
CIB/RILEM Joint Committee on Sustainable Low Slope Roofing (2001), that established the 21 Tenets
of Low Slope Roofing (which have become the international standard bearer for sustainable low slope
roofing), that the essence of sustainability is ‘Long Term Service Life’, then designers must move
beyond current fad and use empirical experience and innovative methods to achieve sustainable roof
systems. This paper will present two roof system design case studies in which the author called upon
his 25 years of roof system design experience, empirical in-field service performance, and
architectural background to design roof systems whose long-term service life expectations are 50
years or more. Design concepts, details and photographs from these projects will be presented.
1. Introduction
Prior to any new concept implementation, such a Sustainable Roof Systems, the ultimate end goal
should be clearly and definitively stated and understood through the definition of terms. Make no
mistake about it: Sustainability is about long-term performance. Additionally, in combination with the
concept that roofs are systems that are only as good as their weakest component, the idea of
sustainability must be holistic in approach, involving the entire roof system. The inclusion of green
parameters, which are only one part of sustainability, and without the other key parameters will not
result in the realization of long-term performance. Thus, the keys to sustainability are not independent,
but must act synergistically so that they result in a system that results in long-term performance,
achieving the owner’s and designer’s desired sustainability goals.
Keys to sustainability that have been developed over years and have been proven by the test of time
are as follows (Hutchinson 2007):
1. Achieving the owner’s buy-in.
2. Be aware of budget appropriateness. Quality cannot be achieved with under-developed budgets.
3. A roofing system designed for long-term performance.
4. Use of quality materials that support long-term performance.
5. Inclusion of thermal insulation performance appropriate for the region in which the roof is
installed.
6. Environmental considerations.
7. Use of qualified contractors.
8. Quality assurance.
9. Management of rooftop equipment.
10. Maintenance.
The following is a discussion of each of the above-mentioned 10 keys to sustainability.
1.3 Materials
It continues to amaze this author that new materials (note: new material are not only those coming on
the market, but also those whose formulation are revised, changed and modified in an attempt to
correct in field performance issues). In the U.S.A. the primary example is TPO (thermoplastic
polyolefin membrane), which are so readily and easily accepted by the building design community,
when there is no proven history of performance for these materials in the roofing industry, nor in the
region or locality in which the materials are being installed. Therefore, a key of all materials
implemented on a new roof system looking to achieve sustainability is that each material has an
established history of long-term in situ performance for the system in which it will be installed. The
roof installation methodology should also be appropriate for the roof system involved.
Additionally, if possible, roof system materials should also be restorable, recyclable or reusable after
their service life. What is the availability of roof system materials and their recycling programs in the
locale in which the system it is being installed? Can the roof be restored? One fine example is the
restoration of EPDM roofs after they are either out of warranty or their service life, resulting in an
extended service life.
2. Case Studies
In the past 5 years this authors has been involved with educational, healthcare and fortune 50
companies whose building stock is substantial, and for whom the concept of long-term service life is
especially appealing. The decided upon long-term service life goal was to exceed 30 years. Design,
detailing, material selection, contactor qualifications, and on site observation raised the authors
prediction of long-term service life to 50 years.
A couple of the main parameters (concerns) of designing roof systems for 50 years are: 1.) How will
the material perform under UV radiation, and its potential deleterious effects on the roofing material,
and 2.) What type of in situ performance does the material have. In regard to the case study projects,
the material selected was EPDM. This author’s experience with the material is excellent, with
numerous projects installed over 25 years ago still performing. In fact, studies by SKZ (Süddeutsche
Kunststoff Zentrum) in Würzburg, Deutschland (Zahn 2003), declared that EPDM should have a
service life of 50 years. Furthermore, this author and the EPDM Roofing Association (ERA) undertook
a study looking at EPDM roofs over 30 years old, tested the material for the physical characteristics,
and found that 90% held characteristics at levels above the minimum level for new materials. A second
study heat aged the 30 year in situ material for an additional 30 year equivalency, and then performed
physical characteristics tests. Not surprisingly the results were positive (a paper on the results and
conclusion of these findings entitled Characteristics of In-Service, 30-Year Old EPDM Roof
Membrane will be presented at the 2014 Durability of Building Materials and Components Conference
in Sao Paulo, Brasil, in October 2014).
2.1 Case Study 1: Hawthorn School District 73: Middle School North
Client: Hawthorn School District 73
Building: Middle School North
Location: Vernon Hills, IL (20 miles north of Chicago)
Building Use: Middle School for grades 6-8
Size: 145,000 ft2 (13,471 m2)
Height: 1 and 2 Story: Max height: 28 feet (8.5 m)
Project Goals:
• 30 Year Service Life: Increased to 50 by author after design
• Increase thermal value by 200%
• Budget: $1,650,000 (approx. 1,218,200 € / 10,816,500 kr)
• Roofing to be completed during summer break: June to mid-August
A B
FIG 2A -B. Roofing crews set the 90-mil EPDM into the full coverage two component spray
polyurethane adhesive, which will be rolled into place with a 90 # (41 kg) water roller. Note the two
layers of insulation with offset joints and the roof edge sheet metal installation on the roof edges
beyond. In the right photo the 45-mil vapor retarder has been installed, and the raising of the roof
curbs and edges, replacing roof drains, and modifying mechanical equipment can take place.
2.1.3 Construction:
The project was bid in early spring 2012 with the main construction commencing in early June. The
roof installation was observed by a Registered Roof Consultant (RRC) and took approximately 4
months. Roof edge sheet metal took another 2 months.
Contentedly, the owner reported an approximate drop in the total building energy usage in the first
year.
2.2.3 Construction:
The construction commenced in mid-March 2013 and was completed in late June. The author observed
the construction and several anomalies were observed. The existing roof edge wood blocking was
found to be in good condition and was allowed to remain. Even with a high caliber contractor, their
project organization and crew management was less than desirable and lengthened the construction
time. Undulating concrete roof decks challenged the tapered insulation, and the required detailing on
what many would consider the simplest of roof systems took considerable time to install.
3. Conclusion
In over twenty-five years of roof system design, this author has had the opportunity to design
numerous long-term performance roof systems, even before the concept of sustainability was popular.
To this author’s knowledge, 90% of all the roof systems that he has designed are still in place and
functioning properly. Through empirical observations and discussions with roofing contractors,
building scientists, and others in the building industry such as manufacturers and owners, the keys to
sustainability have actually been known intuitively by all of us for years. Now they are stated and need
to be put into universal practice. Building owners, architects, the roofing industry, ‘Green Activists’
need to move forward with enthusiasm, to embrace the concept of long-term service life and thus
sustainability (long-term performance), and to accept the challenges to achieving roof systems that
will perform for half a century: for our clients, for ourselves, and for our children.
References
CIB W083 / RILEM 166-MRS Joint Committee on Roofing Materials and Systems. 2001. Towards
sustainable roofing. CIB Publication No. 271. http://cibworld.xs4all.nl/dl/publications/Publi271.pdf
Desjarlais A.O. et al. 2008. Evaluating the energy performance of ballasted roof systems. ORNL
Report Number UF-04-396
Hutchinson T. W. 2007. Keys to sustainable roofing. Roofing Canada, p. 16-30.
Hutchinson T.W. 2014. Characteristics of in-service, 30 year old EPDM roof membrane. Proceedings
of 2014 Durability of Building Materials and Components Conference in Sao Paulo, Brasil.
U.S. Department of Energy. 2011. Building Energy Data Book.
Zahn A. et al. 2003. Final Results of Study to Estimate the Service Life of Hertlan EPDM roof
sheeting. Life expectancy of EPDM roofing membrane over 50 years. SKZ (Süddeutsche
Kunststoff Zentrum) Final report no. 37236/99-IV. SKZ - TeConA GmbH, Würzburg, Germany.
KEYWORDS: External insulation, Cellular plastic material, Fire properties, Thermal insulation,
Moisture, Long term stability, ETICS
SUMMARY:
External thermal insulation composite systems (ETICS) can be used as extra insulation of existing
buildings. The system can be made of cellular plastic materials or mineral wool. There is a European
Technical guideline, ETAG 004, that describe the tests that shall be conducted on such systems. This
paper gives a comparison of systems with mineral wool and cellular plastic, based on experience
from practice and literature. It is important to look at the details in the system and at long time
stability of the properties such as thermal insulation, moisture and fire. Investigation of fire
properties must be done before utilisation of the system, including the risk of fire spread from one
storey to the next for practical solutions. An elaboration of fire spread risks require thermo physic
knowledge about ignition temperatures, critical radiation, upward flame spread velocities etc. of the
actual insulation.
1. Introduction
External insulation is a good solution for many buildings as it reduces the effect of cold bridges and in
most cases does not need to disturb the inhabitants in the building. It can be used for both single
family houses and blocks of flats. In most cases it will change the architecture of the building as we
get a new outer layer. This can be a problem for old houses, where we will keep the exterior look. The
external insulation can be mineral wool or cellular plastic material covered of an external layer of
plate material or plaster. There is a European Technical guideline (ETAG 004 from 2011) for external
insulation systems. The systems are called ETICS, an abbreviation for External Thermal Insulation
Composite Systems. The guideline describes the tests for each component and the system that has to
be done to get a CE-mark. The acceptance of the system is based on the whole system, so it is not
allowed just to change for instance the insulation material without a new approval. The European
Association for External thermal insulation (EAE) has made a report (EAE 2011) that describes a
quality system for external extra insulation. In Denmark it has not been allowed to use plastic
insulation as external insulation and an overview of the experience from other countries is interesting.
We will specially focus on the use of cellular plastic material and the problems related to fire. But we
will also discuss the use of mineral wool and compare systems with plastics and mineral wool, as they
have different properties on some important parameters.
2. Components
The systems consist of:
Adhesive (binds the insulation to the existing wall surface)
For anchors a European approval (ETAG 014) defines the properties needed for plastic anchors. For
each of the insulation materials are also found European standards that define the tests that have to be
done.
An external insulation can be performed on existing walls. The outer surface material of the wall will
influence the solution.
Heavy constructions. This includes walls of lightweight concrete, concrete and masonry of
different types. In these cases it will normally be easy to fix the insulation systems with anchors
and/or glue.
Light constructions. This can be wood, plywood and gypsum. The mechanical strength of the
material can be limited and there is a higher risk for movements of the surface. So the fixing is
more complicated. In this case it is important to keep water away from the light construction.
In ETAG 004 are described test methods that must be used to evaluate the system for different outer
surfaces. In the EAE guide is a check list to avoid mistakes in design and execution.
The thermal insulation material must have a CE-mark and the most important factor is the thickness
and the thermal conductivity. The thermal conductivity shall be declared by the producer. For external
insulation we can normally use the declared value as the insulation should be kept dry. The declared
values for mineral wool are 0.034 -0.040 W/mK. For cellular plastic materials there is more variation
depending on the type – polystyrene, polyurethane, polyisocyanurate or phenol foam. For polystyrene
the value will also depend on the production method – expanded polystyrene (EPS) 0.034 – 0.041
W/mK and extruded polystyrene (XPS) 0.025 – 0.038 W/mK. For phenol, polyurethane and
polyisocyanurate foam, the values can go down to 0.022 W/mK. This is caused by gases used in the
production process. The declared value takes into account the change in gas concentration during the
lifetime.
The insulation must be fixed to the wall without air gab between wall and insulation material, as this
will reduce the heat resistance of the system. For EPS or mineral wool it is normal to glue all over the
surface or along the borders of each plate and some points in the centre. It is very important to
consider details as corners and connections, so thermal bridges are reduced. The insulation layer must
be fixed with anchors if the glue is not sufficient. In the description of the ETICS it must be described
which type of anchors (and how many) that are needed. Insulation of plastics does not tolerate
ultraviolet-radiation as it can degrade the material, so is it important that the surface plaster is put on
3. Thermal properties
The systems must keep the thermal insulation properties over time, so we still get the energy savings
we calculated for the system, also after 10-20 years. We will normally expect that if we calculate the
transmission heat loss after the European standards then this is a good estimate for the energy savings.
It is possible to measure U-values before and after the installation, but it is very difficult in real
buildings to be sure that the indoor (and outdoor) climate is the same and the inhabitants behave in the
same way. So a theoretical calculation is sufficient, but we have to look at possible problems.
For most materials there will be no change in thermal conductivity over time. Only cellular plastics
with a foam gas that can diffuse out will increase the thermal conductivity over time, but this should
be included in the declared conductivity.
With air gaps between the insulation material and the wall or between the insulation panels, then cold
outer air can flow into the insulation layer, and we do not achieve the expected U-value. This is
important to avoid in the construction phase. Furthermore we must ensure that the fixing of the panel
not loosen during the lifetime of the system.
Moisture in the insulation will increase the U-value. The moisture can come from the outside in the
form of driving rain or from the inside, if the moisture barrier is poor. Most critical is driving rain. If
the extra moisture only comes in short periods it will probably dry out.
4. Fire properties
Cellular plastics are made of plastic which is foamed by injecting a gas. The final foam material
consists of approximately 95-98% gas. Cellular plastic insulations are flammable, unlike mineral
wools which are non-combustible. Fire characteristics for cellular plastics are, of this reason,
important to investigate. The risk of fire spread, especially vertically from one storey up to the next, is
important for practical solutions. An elaboration of fire spread risks requires thermo physic
knowledge about ignition temperatures, critical radiation, flame spread velocities etc. Additives as fire
retardants could be relevant to slow down the fire speed or reducing the flammability. This is
discussed further below.
Cellular plastic insulation materials have very different fire characteristics, and there exist a number
of opportunities to improve these through the addition of varying amounts of fire retardants. The
We do not have one common European test method for facades. Some countries apply their own
national tests, other utilize test methods, which in principle are developed for other purposes. This
applies to the use of the EN 13823 single burning item (SBI) test for construction products` reaction
to fire. Sweden and Norway use a Swedish full-scale test (SP Fire 105).
Cellular plastic has a large surface area relative to body mass, i.e. a large specific surface, and the
material has great access to oxygen. These factors result in a fast and clean combustion of the
material. The cell structure results in a low density of the material and consequently a relatively low
energy production by combustion, measured by volume unit (Troitzsch 2004). (DBI 2012) provides an
overview of the typical cellular plastic type density, thermoplastic properties and energy, and is
summarized below in Table 1.
Table 1. Overview of selected physical properties of typical plastic-based insulation materials (DBI,
2012).
The ignition temperature is dependent on whether there is a pilot flame present or not. If the foam
only is affected by heat radiation, the possibility of ignition depends on the exposure time, the
radiation intensity and the size of the exposed area.
Table 2 from (DBI 2012) shows examples of ignition temperatures and critical thermal radiation of
various cellular plastic insulation materials. The critical thermal radiation is the lowest thermal
radiation which may cause ignition. Due to the materials very product specific properties, there may
be a discrepancy between the table and some product values.
PF 490 °C 450 °C - -
Table 2. Examples of ignition temperatures and critical thermal radiation of the typical cellular
plastic insulation materials (DBI 2012).
By the addition of a fire retardant (additive) during production, the fire performance of the foam is
changed. An additive can be used to alter cellular plastics ability to shrink or melt, char, slow down
the fire speed or to reduce the flammability (DBI 2012). For example, the bromide flame retardants
are widely used.
It may be noted that, for example ISOBYG`s EPS burn different than normal EPS, due to adding of
flame retardant substances, which make the material does not burn in a self-sustaining reaction. That
is, it turns off automatically when there is no direct flame exposure (ISOBYG 2012).
There are a large variety of additives. This coupled with a myriad of different combinations and
proportions of different types of cellular plastics, makes it very difficult to accurately describe the
additives effect on cellular plastics fire performance. Therefore, a fire-related classification of a foam
material must be determined on the basis of a test of the exact composition of the foam and
additive(s) (DBI 2012).
Fire and flame retardant additives have the greatest impact by initial heating, from small heat sources,
i.e. corresponding to fire effects in the early course of the fire. In a fully developed fire, flame
retardants have no effect on the fire performance (NFPA 2003).
The fire risk is not yet fully documented for cellular plastic insulation on facades. The risk of external
flame spread seems to be the most critical parameter in this context. It is recommended that flame
spread from one storey to the one above, should be further investigated in order to prove a fire safety
level according to the Building Regulations. At the moment, there exists only the Nordic SP Fire 105
test, which in fact is developed for other types of facades. By the way the test is not suitable for
handling molten material dripping from the burning facade.
The lifetime for plaster is important as it will give the lifetime for the ETICS systems. The
information from literature is that the lifetime is similar to a masonry wall or a wooden panel. In all
cases it is assumed that ordinary maintenance is carried out.
In Sweden it has been popular to build high insulated wooden wall with insulation and plaster on the
outside as seen in Figure 1. An investigation of 800 buildings (Samuelson and Jansson 2009) shows
that moisture will come into the construction. There is damage in 55 % of the buildings, mostly in
houses with EPS as insulation. Using mineral wool as insulation reduces the damage to 32%. This is
still a very high damage percentage. The problems are related to details as around windows, doors,
shading roofs, balconies and cracks in the plaster at fixation of light or sunblinds. Testing the systems
before it is used is important and also to reduce the number of special details, where there is a risk of
damage. Later reports (Jansson 2011) confirm the results with EPS giving a high risk of moisture
problems for light walls.
FIG 1. Sketch of system on a light wall from Samuelson 2009, 1- plaster layer 2- external insulation 3
gypsum plate 4 internal heat insulation 5 vapour barrier 6 gypsum plate
The Swedish problems are related to water in the insulation. Drying out of moisture in the external
insulation layer is important. The most efficient way is to use a ventilated space behind the plaster or
outer plate in the system. Swedish calculations from Falk and Sandin 2013 has shown large variations
6. Conclusion
External insulation of facades has, in Denmark, typically been done with systems of mineral wool as
insulation material. This has worked very well, but it is interesting to look at the use of systems with
cellular plastic material. Danish experience with plastic material in external insulation system is
rather limited, primarily due to limiting building regulations, especially in relation to fire. Experience
is taken from Sweden and Germany where it is widespread.
The fire risk is not yet fully documented for cellular plastic insulation on facades. The risk of external
flame spread is found to be the most critical parameter in this context. It is recommended that flame
spread from one storey to the one above, should be further investigated in order to prove a fire safety
level that fulfils the building regulations.
There is no extensive documentation showing that houses with external placed insulation should have
a higher general fire risk than buildings without external insulation.
References
Babrauskas, V. 2003. Ignition Handbook. Interscience Communications Ltd. 2003.
Davies, J.M. 2001. Lightweight Sandwich Construction. Blackwell Science Ltd. 2001.
EAE 2011; European guidelines for the application of ETICS (external thermal insulation composite
system) - Quality with system, EAE – European association for External thermal insulation
composite systems, March 2011.
EN 13823:2010 Reaction to fire tests for building products. Building products excluding floorings
exposed to the thermal attack by a single burning item. 2010
ETAG 004 2011. Guideline for European technical approval of external thermal insulation composite
systems with rendering, European Organisation for Technical Approvals, Brussels, 2011.
ETAG 014 2011. Guideline for European technical approval of plastic anchors for fixing of external
thermal insulation composite systems with rendering. European Organisation for Technical
Approvals, Brussels, 2011.
Falk, J. and Sandin, K.; 2013. Ventilated rain screen cladding: Measurements of cavity air velocities,
estimation of air change rates and evaluation of driving forces, Building and Environment vol. 59
164-176. 2013
Künzel, H., Künzel, H.M. and Sedlbauer, K. 2006; Long-term performance of external thermal
insulation systems (ETICS), Architectura 5 (1) 2006, 11-24.
NFPA 2003. Fire Protection Handbook, 19. edition, Volume II. National Fire Protection Association.
2003.
Rakic, J. 2003. Fire rated insulated (sandwich) panels. Fire Australia. May 2003.
SP Fire 105. Swedish standard for large scale testing of facades system. SP Technical Research
Institute of Sweden, Borås, Sweden.
Sørensen, L.S., Nielsen, A. and Wittchen, K.B. 2013. Udvendig efterisolering med plastbaserede
materialer (external insulation with plastic based materials). Statens Byggeforskningsinstitut,
Aalborg Universitet, København.
Troitzsch, J. 2004. Plastics Flammability Handbook. Published by Carl Hanser Publishers, Munich.
2004
SUMMARY:
Single leaf masonry walls with the application of thermal insulating mortar can be an interesting
solution for façades, providing the performance of its components can lead to a competitive thermal
resistance. The determination of the walls U-value will therefore be of most importance.
Laboratory tests were conducted in a wall composed of a thermal insulating mortar applied on a
single leaf masonry of complex shape units. The objective was the evaluation of the U-value and the
validation of calculation methodologies. The tests included material properties evaluation and U-
value determination of the wall using a climate chamber and heat flow meters.
Numerical simulations were performed with the focus on the calculation of the thermal resistance of
masonry using two approaches. One was a simplified 2D model based on EN ISO 6946 calculation
methodology, and the other one a 3D model based on the finite element method (FEM).
The numerical results obtained from the simulations were compared with laboratory measurements
leading to a difference of about 5%. The differences observed and the results obtained in the material
characterization phase supported recommended strategies that can be used to improve the feasibility
of U-value calculations for the studied wall, including manufacturer declared moisture dependence of
key material properties.
1. Introduction
A growing interest in envelope solutions based on thermal insulating mortars applied on masonry
walls of optimized geometry can be observed nowadays. The effective U-value of these solutions
should however be carefully evaluated. Practitioners often use reference values available in literature
or declared by manufacturers to determine U-values that will be computed in building thermal
calculations. The case of masonry walls is typically sensitive in the sense that manufacturers
frequently base the declared values on calculations using (EN-ISO-6946, 2007) methodology. For the
case of masonry using complex shape units and thermal rendering the definition of the expected U-
value raises additional problems that have to be tackled so that manufacturers declared values can be
trusted.
2. Studied System
A specimen representing the wall system under analysis was installed in a climatic chamber. It
includes a thermal insulation mortar applied as external rendering. The specimen is 3.61 m2 (1.9 m x
1.9 m) and its constitution is described in TABLE 1 and FIG. 1
3. Material Properties
The thermal conductivity (λ) of the masonry units base materials and of the thermal rendering was
determined. Conductivity can be affected by different factors. The λ is moisture dependent (Freitas,
1994), especially for these materials. Also, according to (Stahl et al., 2012) the pressurization used on
the plastering machine can affect the actual performance of the insulating mortar. Hence, the actual
insulating mortar applied on the wall was studied. Two 30 x 30 cm2 panels were extracted from the
specimen (FIG. 4). The λ of the masonry units material was determined from shell samples removed
from the units. The determination of the thermal conductivity was been carried out using a thermal
conductivimeter by Hot Wire “CT-mètre”. The dry value was calibrated with measurements using the
guarded hot plate according to (EN-12664, 2001).
The λ values of the samples of thermal insulating mortar (A) (bulk density ~400 kg/m3) and masonry
units (B) (bulk density ~440 kg/m3) was obtained for different moisture content values, with results
displayed in TABLE 2 and FIG. 3. In FIG. 3, the thermal conductivity as a function of moisture
content of two types of lightweight concrete (Kumaran, 1996) is also included for comparison. The
lightweight concrete and aerated concrete values from literature demonstrate an upgraded
performance of the material applied in the masonry units of the studied system. The thermal insulating
mortar performance is similar to the polystyrene concrete from literature, although its moisture
dependence is higher.
TABLE 2 – Thermal conductivity values of the samples in various conditions.
W λ
Sample Sample conditions
[%-kg/kg] [W/m.ºC]
A Dried at 70ºC 0% 0.095
A Laboratory conditions (RH ≈ 70%) 2% 0.107
A Intermediate moisture content 28% 0.188
A Saturated 47% 0.273
B Dried at 70ºC 0% 0.178
B Laboratory conditions (RH ≈ 70%) 1% 0.221
B Intermediate moisture content 12% 0.320
B Saturated 19% 0.360
0.4
0.3
0.2
0.1
0.0
0% 10% 20% 30% 40% 50%
W [kg/kg - %)
FIG. 3 –Thermal conductivity values obtained for the thermal insulating rendering sample.
FIG. 4 – Extraction of samples of insulating mortar from the specimen and shells from the masonry
units
4. Numerical Simulations
Numerical simulations were performed to determine the thermal resistance of the masonry and the
resulting U-value of the wall system, with two different procedures:
• 2D model based on a simplified geometry of the unit and vertical joints, and on the
calculation methodology of (EN-ISO-6946, 2007) for the determination of the thermal
resistance of the masonry unit with air spaces/voids and (EN-1745, 2012) used on the
calculation of the overall resistance of the wall, i.e. with the influence of the mortar bed joints
and rendering layers;
• 3D model based on the finite element method (Simulia), in which micro modelling of the
masonry system was performed with an equivalent thermal conductivity for the air spaces
(voids) calculated trough the method indicated in (EN-ISO-6946, 2007).
5. Laboratory Tests
The specimen was conditioned by the climatic chamber Fitoclima EDTU 1000 (FIG. 5). The
temperature provided by the device ranges from -50 to 180° C with an accuracy of ± 0.5°C while the
relative humidity range is set from 10 to 98% with an accuracy of ± 2%. The test specimens were
constructed in a detachable part of the climatic chamber (FIG. 5). The wall surface with thermal
rendering is conditioned by the climatic chamber and the opposite surface is exposed to laboratory
conditions. For the determination of the heat flow through the specimen, the equipment presented on
the FIG. 6 was used. This equipment has an accuracy of 5%. The temperatures on the surfaces were
determined using thermocouples of the type T as shown on FIG. 6. The temperature and the relative
humidity near the specimen inside and outside of the climatic chamber was measured with Onset
HOBO data loggers with a temperature range from -20 to 70 ºC with an accuracy of 0.35ºC while the
relative humidity range is from 10 to 90% with an accuracy of 3% (FIG. 6).
6. Conclusions
The main conclusions of this study are the following:
• The 2D and 3D numerical simulations demonstrated good accuracy with the experimental
results, although 3D gave a better approach;
• A difference of about 5% was found between the numerical simulations (3D) and the
experimental determination of the U value;
• Considering the penetration of laying mortar into the voids of the masonry units, reduces the
difference between measured and calculated (3D) U-Value to 4%.
• A numerical procedure based on a 3D FEM model, in which the thermal conductivity values
used in the model are obtained from laboratory tests made with samples removed from the
masonry components, can be useful to determine the U values for several moisture conditions
without using too many laboratory resources;
• The U-value of this wall system is very sensitive to moisture content even in the hygroscopic
range and the manufacturers should clearly state it;
• This study is part of a work in progress and additional analysis will be performed on the
hygrothermal behaviour of these systems.
7. Acknowledgements
The authors would like to thank Secil and Secil Argamassas for their technical support.
References
EN-1745 2012. Masonry and masonry products. Methods for determining thermal properties.
EN-1934 1998. Thermal performance of buildings. Determination of thermal resistance by hot box
method using heat flow meter. Masonry. 1998-09-30: CEN/TC 89.
EN-12664 2001. Thermal performance of building materials and products - Determination of thermal
resistance by means of guarded hot plate and heat flow meter methods - Dry and moist
products of medium and low thermal resistance.
EN-ISO-6946 2007. Building components and building elements -- Thermal resistance and thermal
transmittance - Calculation method.
FREITAS, V. P. 1994. A influência da humidade na resistência térmica de elementos de construção.
GHAZI WAKILI, K. & TANNER, C. 2003. U-value of a dried wall made of perforated porous clay
bricks: Hot box measurement versus numerical analysis. Energy and Buildings, 35, 675-680.
ISO-9869 1994. Thermal insulation -- Building elements -- In-situ measurement of thermal resistance
and thermal transmittance. ISO/TC 163/SC 1.
KUMARAN, M. K. 1996. IEA ANNEX 24 - Heat, Air and Moisture Tranfer Through New and
Retrofitted Insulated Envelope Parts (Hamtie) - Task 3: Material Properties, Laboratorium
Bouwfysica - Department Burgerlijke Bouwkunde - K.U. Leuven Belgium.
SUMMARY
An experimental facility consisting of a wall exposed to the external environment and installed on a
cell with indoor controlled conditions is studied. The wall is timber framed with wood fibre panels for
insulation. The first experimental results are presented. The results show the predominant effect in the
short term of the temperature and insolation on the distribution of moisture in the walls, when for the
long term it is rather the water content of inside and outside atmospheres that directs the walls’
moisture distribution.
NOMENCLATURE
Symbol Definition Unit
λ Thermal conductivity mW/(m.K)
θm Moisture content kg/kg
HR Relative Humidity %
µ Vapour diffusion equivalent air layer thickness -
Cp Specific heat capacity J/(kg.K)
T Temperature °C
C Water adsorption coefficient kg/(m².min0.5)
1. Introduction
To study Heat, Air and Moisture (HAM) coupled transfers in building envelopes, many numerical
studies have been done and many experimental facilities under controlled conditions or outdoor
climate have been established [Kêdowidé et al., 2012]. Most of these studies concern realistic
constructions, with high air and water tightness together with a vapour barrier, thereby limiting the
impact of mass transfers on energy performance of the envelope. Yet moisture movements influence
the temperature fields inside the walls [Labat et al., 2012] and indoor comfort [Teodosiu et al., 2003].
There is also a persisting need for precise assessment and understanding of the impact of mass transfer
on the performance of highly hygroscopic wall assemblies, as well as for detailed validation of
dedicated simulation tools. In order to better assess the effects of coupled transfers within wall
assemblies, a new experimental work with highly vapour-permeable and very hygroscopic materials
was used in the work presented below.
The walls are then coated with a 3mm roughcast to protect the insulation from liquid loads such as
rain (C<0.20 kg/m².min0.5) [EN 998-1] without being a barrier to vapour diffusion (µ<25).
During the first period of the experiment, reported in the present paper, the 2 walls are identical, and
exposed to the same indoor climate.
Interior
FIG 1: Position of the sensors: General view of the tested façade from the outside (left); Sensors in the
thickness of the wall (right), dimensions in cm
3. Implemented tests
Tests on presented experimental facility began in March 2013, which corresponds to the end of winter
- early spring. The figure 2 shows an overview of the experimental indoor (int) and outdoor (Ext)
conditions over 5 months: temperature, relative humidity and vapour pressure. According to the
weather data, recorded results can be divided into two main sessions.
4. Results
4.1 Spring session
Here will be presented the first relevant test done with the experimental facility. Comparing pairs of
sensors placed at the same position in the two tested wall showed that the responses were very similar;
therefore, only one wall’s results will be presented.
FIG 4: Temperatures, Vapour pressures and radiation on the walls during the steps experiment
FIG 5: Temperatures, Vapour pressures in the wall and radiation on it when turned to the south
FIG 7: Temperatures, Vapour pressures and RH on the walls when vapour pressure dropped
5. Conclusions
A test made on a wall exposed to external environment was studied. The wall is a timber frame with
wood fibre insulation, and a vapour permeable outer coating. The results show that for a very
hygroscopic and vapour permeable wall, on the one hand the movement of vapour are directed by
temperature fluctuations in the short term and the vapour pressure gradient between the exterior and
interior in the long term, and secondly that the vapour pressure variation reciprocally influence the
temperature of the hygroscopic material in the wall. Numerical studies should help better understand
and assess these effects.
6. Acknowledgments
This study is done for the HYGROBAT (ANR-10-HABI-0005) project, financed by the National
Agency of Research of France (ANR).
References
Desta, T. Z., Langmans, J., & Roels, S. , 2011. Experimental data set for validation of heat, air and
moisture transport models of building envelopes. Building and Environment, 46(5), 1038-1046.
Kêdowidé Y., Woloszyn M., Le Pierrès N., 2012, Synthèse sur les dispositifs expérimentaux d'étude
des transferts couplés de chaleur et d'humidité dans les parois des bâtiments ; Colloque International
Francophone d’Energétique et Mécanique (CIFEM), Ouagadougou.
Labat M., Noel M., Woloszyn M., Piot A., Garnier G., Roux J.-J.. , 2012, Comparison of moisture and
temperature measurements for six envelope types with different moisture transfer properties. 5th
International Building Physics Conference. Kyoto, Japan, 28-31 May, 2012.
NF EN 998-1 (Décembre 2010) : Définitions et spécifications des mortiers pour maçonnerie – Partie 1
Teodosiu C., Hohota C., Rusaouën G., Woloszyn M., 2003, Numerical prediction of indoor air
humidity and its effect on indoor environment, Building and Environment, Volume 38, Issue 5, Pages
655–664.
1
KU Leuven, Building Physics Section, Department of Civil Engineering, Belgium
SUMMARY:
Several studies show that the actual ‘as-built’ thermal performance of the building envelope can differ
significantly from the theoretical, calculated value. Characterisation of building envelope
performances based on in situ dynamic measurements can help to bridge this gap between ‘designed’
and ‘as-built’ performances. A common method to evaluate the thermal performance of a building in
situ is the co-heating test, which is a quasi-stationary method based on linear regression analysis of
dynamic measurement data.
After a short state-of-the-art on the co-heating test methodology, the limitations and opportunities
associated with the use of the co-heating test method to characterise the thermal performance of
buildings are investigated on the basis of a real full-scale experiment: a co-heating test performed on
a terraced house in Herstal, Belgium. Renovation induced drops in the overall heat loss coefficient of
the dwellings are characterised by the test method and compared with calculated values.
1. Introduction
In order to reduce the energy use of buildings, several countries have put forward more stringent
demands on energy performance of new and renovated buildings. Without exception, these supervised
buildings are characterised or awarded a label in the design phase. A theoretical energy use calculated
on the basis of building plans and specifications determines the performance category. An important
distinction needs to be made, however, between this theoretical energy performance and the actual as-
built performance. Several studies have shown that these can differ rather significantly (Bell et al.
2010, Lowe et al. 2007).
The energy performance of a building is essentially determined by the (1) thermal characteristics of
the building envelope, (2) installed services and (3) building usage. As the latter is not easily
predicted nor controlled, the first two are decisive in achieving the envisaged building energy
performance, both for new buildings and renovations. Hence, the thermal performance
characterisation of a building envelope represents a crucial first step to bridge the gap between its
designed and as-built energy performance. A common method to evaluate the thermal performance of
a building envelope in situ is the co-heating test.
To the knowledge of the authors, (Sonderegger and Modera 1979) is first to mention the use of
thermostatically controlled electric heating experiments to assess a building’s energy performance.
Real full-scale dwellings were alternately heated using the building’s own services and electric heaters
with known efficiency. Hence the name co-heating. Ever since it’s conception, the co-heating test
method has been used to (Sonderegger and Modera 1979):
where Asw,* now represents the solar aperture coefficient corresponding to the sole global solar
radiation component selected as input. For instance, in a next section, measurements of global solar
radiation on a horizontal surface, qsw,hor, is selected.
Where c describes the discrepancy between the regression model fit and the actual aggregated
measurement data.
Essentially, three options can be discerned:
As illustrated in Figure 1 (b), an estimate of HLC is then given by the intercept. Note that this
mathematical transformation implicitly forces the above described multiple linear regression
through zero. In both of the earlier mentioned methods, a non-zero intercept is possible due to
discrepancies between the measurement data and the assumed stationary model to which it is
fitted. In the third method, these errors are included in the HLC-estimate.
(a) Basic multiple linear regression: ΔT as (b) Multiple linear regression when only
descriptive variables, thermal lag not considering winter months, ΔT and
taken into account, intercept c. qsw,* as descriptive variables, taking
into account thermal lags and forcing
intercept through zero (c = 0).
FIG 2. U-value estimation through multiple linear regressions applied on data acquired during a
simulated co-heating test on an insulated cavity wall component. Different subsets of the
measurement data are considered: with different starting dates (vertically aligned points) and
various durations (along the x-axis). The data points are coloured according to the average Qh
during the considered measurement data subset: red points indicate a high average Qh, blue
points indicate a low average Qh.
FIG 3. Investigated terraced house in Herstal, FIG 4. Co-heating test equipment: heat sources
Belgium. (a) Front façade; (b) Garden façade. controlled by thermostats, ventilators,
temperature sensors and pulse meters spread
throughout the investigated dwelling. Sensor
data (including weather station data) is
transmitted to a central logger located on site.
FIG 5. HLC estimates through multiple linear regressions with Qh as a function of ΔT and qsw,hor,
taking into account thermal lags and forcing regression surface through the origin, applied on
measurement data acquired during co-heating test on a terraced house in Herstal, Belgium. Different
subsets of the measurement data are considered: with (1) different starting dates and (2) different
durations. The data points are coloured according to the average Qh over the course of the considered
measurement data subset: red points indicate a high average Qh, blue points indicate a low average
Qh.
TABLE 2. Results of applying simple and multiple linear regression on daily averaged co-heating test
data, before and after renovation step 1 and 2: coefficient of determination R2, estimates for HLC and
Asw,*, p-value Pr(>|t|). The p-value indicates the probability the variable is not significant.
R2 HLC Std.Error Pr(> |t|) Asw,* Std.Error Pr(> |t|)
- (W/K) (W/K) - (m2) (m2) - ΔHLC
Before 0.9983 296.55 13.51 3.98e-09 3.15 3.85 0.433
After step 1 0.9935 185.63 6.27 2.14e-15 4.76 0.69 3.70e-06 110.92
After step 2 0.9899 157.84 9.62 4.54e-10 2.48 0.68 3.06e-03 27.79
Both ΔT and qsw,* prove to be significant variables in all cases but one. Before renovation step 1, the
average solar radiation was rather limited, which explains the fact that the solar radiation as a
descriptive variable proves to be less significant in that case (p-value > 0.05).
6. Acknowledgements
The authors wish to thank Knauf Insulation, for setting up an expert network concerned with
investigating the gap between the ‘designed’ and ‘as-built’ thermal performances of building
envelopes
References
Sonderegger, R. C., Modera, M.P. (1979). Electric Co-heating: A method for Evaluating Seasonal
Heating Efficiencies and Heat Loss Rates in Dwellings. In Second CIB Symposium on Energy
Conservation in the Built Environment. Lawrence Berkeley National Laboratory, March 1979.
Everett, R. (1985). Rapid thermal calibration of houses, ERG, Open University Energy Research
Group, Milton Keynes UK.
Roels, S. (2011). Reliable Building Energy Performance Characterisation Based on Full Scale
Dynamic Measurements. Annex text, November 2011.
(http://www.kuleuven.be/bwf/projects/annex58/index.htm)
Bell, M., Wingfield, M. and Miles-Shenton D. (2010). Low carbon Housing. Lessons from Elm
Treed Mews. Technical report, Joseph Rowntree Foundation.
Lowe, R.J., Wingfield, J., Bell, M. and Bell J.M. Evidence for significant heat losses through party
wall cavities in load-bearing masonry construction. Building Services Engineering Research and
Technology, 28 2:161-181, 2007.
Bauwens, G., Standaert, P., Delcuve, F., and Roels, S. (2012). Reliability of co-heating measurements.
First Building Simulation and Optimization Conference. IBPSA, September 2012.
Max Sherman, Estimation of Infiltration for Leakage and Climate Indicators, in Energy and Buildings,
10, 1987, p.81.
SUMMARY:
Unventilated wooden roofs typically need a vapor barrier at the warm side to avoid interstitial
condensation due to vapor diffusion and air leakages from the interior. So-called smart vapour
barriers (SVB) or retarders, could allow condensed moisture, built-in-moisture or moisture from
minor leakages to dry to the interior. The concept of SVBs or retarders have been known for some
time, consisting of a material that changes its vapour resistance according to the level of relative
humidity in the surrounding air. SVB’s should be relatively vapour open during summer conditions to
allow drying to the interior. During winter conditions however, the vapour resistance should
preferably be as high as possible to avoid condensation problems. Various commercial SVB’s have
different properties in regard to the level of their vapour resistance and variation with relative
humidity.
In this laboratory study the level of moisture accumulation during winter conditions have been
investigated for unventilated wooden roofs applying four various types of SVB’s and compared with a
standard roof with polyethylene foil as vapour barrier. Four levels of indoor air humidity were
applied during the test period of three months. The measurement results were compared with
hygrothermal simulations. The results showed that the various SVB’s had very different robustness
against condensation problems during winter conditions, and thus the type of SVB should be chosen in
regard to the expected level of indoor air humidity.
1. Introduction
In Nordic climates a vapour barrier is typically used on the warm side of the building envelope to
avoid interstitial condensation during the heating season due to vapour diffusion and air leakages from
the interior. Often a polyethylene foil is used, which has a very high vapour resistance. Such a vapour
tight product does however not allow moisture to dry to the interior. A more vapour open product
could allow condensed moisture, built-in moisture or moisture from minor leakages to dry to the
interior. This could be especially useful for constructions with reduced or no possibility to dry to the
exterior, such as unventilated flat roofs. It might allow cheaper solutions, making it possible to use
organic (wooden) materials in such constructions or skip the use of wooden preservatives.
The term “vapour retarder”, as opposed to “vapour barrier”, is often used on products that have a
lower vapour resistance than recommended for vapour barriers. Such products may allow for some
drying to the interior. While a vapour retarder has a given constant vapour resistance, some vapour
barriers or retarders are sold (Note: “barrier” is used in the rest of this paper for simplification) on the
European and North American market with adaptable vapour resistance in regard to what is actually
needed. Popular terms for these products are “smart”, “intelligent”, “moisture adaptive” or “humidity
dependent” vapour barriers. The physical behaviour of these products varies, but the main principle is
that the vapour barrier should function as an ordinary vapour tight vapour barrier most of the time,
preventing vapour diffusion into the construction from the indoor air – especially during winter
conditions. If, on the other hand, the construction is wet, for example due to built-in moisture or
2. Laboratory tests
2.1 General
The purpose of these tests was to investigate the risk for undesirable vapour diffusion from interior air
during winter conditions in unventilated wooden roofs with SVB’s, possibly leading to condensation,
moisture accumulation or other moisture problems. SVB’s should be relatively vapour open during
summer conditions to allow drying to the interior. During winter conditions however, the vapour
resistance should preferably be as high as possible to avoid condensation problems.
In this laboratory study the level of moisture accumulation during winter conditions have been
investigated for unventilated wooden roofs applying four various types of SVB’s and compared with a
standard roof with polyethylene foil as vapour barrier. Four levels of indoor air humidity were applied
during the test period of three months. The effect of having some part of the insulation at the interior
of the SVB was also investigated.
The following type of compact unventilated wood frame roof was tested (from exterior side):
Roofing membrane
200 mm glass wool
SVB or PE-foil
13 mm gypsum board
FIG 1. Test box with RH/T-sensor before mounting of SVB (left) and after mounting of SVB and
gypsum board.
To simulate a cold external climate on one side and a warm indoor climate on the other side the boxes
were installed in a wall separating two different climate chambers with controlled temperature and
relative humidity, as shown in Figure 2. A horizontal section through two of the boxes as installed in
the wall is shown in Figure 3. This means the specimens were installed vertical instead of horisontal,
that may potentially give some internal convection in the insulation layer. The insulation was however
installed with no air gaps, so it is assumed that the effect of internal convection is insignificant. Since
the main focus of this test was on the total moisture accumulation, any internal convection is anyhow
considered of little importance. On the cold side a continuous wood fibre board layer was installed to
FIG 2. The experimental set-up showing the 10 test boxes installed in the wall, seen from the warm
side.
FIG 3. Horizontal section through two of the test boxes installed in the wall.
2.3 Materials
The four types of SVBs are;
Delta-Novaflexx
Intello
AirGuard Smart
Hygrodiode 200A
Delta-Novaflexx is a polyamide with an Sd-value ranging from 0,2–5 m. This product is manufactured
under many different names, including MemBrain in the US, and Difunorm Vario in Germany. It was
originally developed in Germany in the late 90’s (Künzel, 1996), and was the first commercial
available relative humidity dependent vapour barrier product for building applications. Most scientific
studies on use of SVBs in building applications are limited to this product only. Intello is a
polyethylene copolymer, where the copolymer is an acrylic with hygroscopic properties. The Sd-value
ranges from 0,25-26 m (Pro Clima, 2006). AirGuard Smart is a polyvinyl alcohol film (with spun bond
polypropylene as reinforcement and protecting layer) with an Sd-value ranging from 0,05 - 102 m
(DuPont, 2012). Hygrodiode 200A is composed of a fabric of synthetic fibres with good capillary
properties laminated with stripes of perforated PE-foil and polyamide-foil. It should be noted that this
product is a further development of the original Hygrodiode mentioned earlier, i.e. the polyamide film
giving the product drying possibilities also for conditions below 100% RH. The Sd-value ranges from
60
Novaflexx
40
20
0
0 20 40 60 80 100
Relative humidity (%)
2.5 Measurements
The rig and test boxes were constructed so that the boxes could be dismounted and weighed regularly
to follow the total drying of the configurations. In addition the temperature and RH at the interface
between the SVB and the 200 mm insulation were measured (logged each sixth minute) in all test
boxes. The RH and temperature of the indoor climate were also continuously logged. The RH/T-
sensors were controlled before and after the test, and the results were adjusted according to this
calibration.
3. Hygrothermal simulations
The laboratory test setup were simulated for all the described configurations, except for the
configuration with Hygrodiode where detailed material properties were missing, using WUFI 1D Pro
5.2 (WUFI 2013). WUFI 1D Pro 5.2 is a simulation program for coupled non-stationary, one-
dimensional heat and moisture transport. The program includes the moisture transport by vapour
diffusion and capillary transport, and takes into account the moisture capacity of the materials. The
35 300
Total accumulated water (g/m2)
30 PE (1) 250
Total accumulated water (g)
Intello (1)
25
Hygrodiode (1) 200
PE
20 Novaflexx (1)
150 Intello
15 AirGuard Smart (1)
Hygrodiode
PE (2)
100 Novaflexx
10 Intello (2)
AirGuard Smart
5 Hygrodiode (2) 50
Novaflexx (2)
0 AirGuard Smart (2) 0
0 21 42 63 84 0 21 42 63 84
Days Days
FIG 5. Total accumulated water in the five roof configurations, all 10 specimens to the left and the
average of two specimens to the right.
In Figure 6 is shown the comparison of measured moisture accumulation with the hygrothermal
simulations with WUFI. For the Novaflexx configuration the simulation compare very well with the
measurements through all cycles. Also for the PE-foil and AirGuard Smart configurations the
simulations and measurements compare very well, with the exception for some deviation for the last
cycle. For AirGuard Smart the simulations give higher moisture accumulations than the measurements
during the last cycle, while the opposite applies for the PE-foil. One explanation for the deviation in
regard to AirGuard Smart could be that the Sd-value in reality is higher for the given RH levels than
we see in Figure 4. New cup-measurements made for AirGuard Smart seem to confirm that
explanation. For Intello the deviation between the simulation and measurements is however bigger, the
simulation showing much lower moisture accumulation than the measurements. Due to the little
deviation between the two test specimens of the Intello configurations we observe in Figure 5, the
most obvious explanation may be that the real Sd-value for Intello is lower than given in Figure 4. It
250
Total accumulated water (g/m2)
200 PE
Intello
Novaflexx
150
AirGuard Smart
PE‐sim
100 Novaflexx‐sim
Intello‐sim
50 AirGuard Smart‐sim
0
0 21 42 63 84
Days
FIG 6. Comparison of measured and simulated (“sim”) total accumulated water in four of the roof
configurations.
22
20
PE (1)
Intello (1)
Relative humidity (%)
18
Hygrodiode (1)
Novaflexx (1)
16
AirGuard Smart (1)
PE (2)
14 Intello (2)
Hygrodiode (2)
12 Novaflexx (2)
AirGuard Smart (2)
10
0 21 42 63 84
Days
FIG 7. Measured RH at the interface between the glass wool and SVB/PE-foil for all test specimens.
The RH at the interface between the glass wool and SVB/PE-foils were measured and is shown in
Figure 7. We observe that even if the various vapour barriers have very different vapour resistances,
this do not influence the RH at the external side of the vapour barrier too much. For the PE-foil
configuration the RH varies between approximately 12 – 14 %, while for the other SVB’s the RH
varies between approximately 14 – 21% - increasing somewhat as the indoor air humidity and water
accumulation in the test specimens increases during the test. According to Figure 7 the average RH at
the interface between glass wool and the SVS are approximately 16%, 17%, 19% and 20% for Cycle
5. Conclusions
In this laboratory study the level of moisture accumulation during winter conditions have been
investigated for unventilated wooden roofs applying four various types of SVB’s and compared with a
standard roof with polyethylene foil as vapour barrier. Four levels of indoor air humidity were applied
during the test period of three months. The measurement results were compared with hygrothermal
simulations.
The results showed that the various SVB’s had very different robustness against condensation
problems during winter conditions, and thus the type of SVB should be chosen in regard to the
expected level of indoor air humidity. However, whether a certain level of moisture accumulation
during winter conditions is a problem, also depends on the drying during summer conditions. It is for
example interesting to observe that while built-in moisture drying experiments for summer conditions
(Geving, Stellander and Uvsløkk 2013; Geving, Thorsrud and Uvsløkk 2013) show faster drying for
Novaflexx and AirGuard Smart compared to Intello and Hygrodiode, Novaflexx show the highest
moisture accumulation and AirGuard Smart the lowest for this experiment for winter conditions.
The simulations of the test conditions compared well with the measurements, except for a few cases –
where one explanation for the deviations could be that the given Sd-values for those products may not
be sufficiently correct or detailed.
6. References
DuPont. 2012. Tyvek ©. Product Brochure. http://construction.tyvek.co.uk.
Fraunhofer IBP. 2010. WUFI 1D (Version 5.1) [Computer Program]. Fraunhofer IBP, Holzkirchen,
Germany.
Geving, S., Stellander, M. and Uvsløkk, S. 2013. Smart vapour barriers in compact wood frame roofs.
Proceedings of Buildings XII Conference – Thermal Performance of Ecterior Envelopes of Whole
Buildings, December 1-5, 2013, Clearwater Beach, Florida.
Geving, S., Thorsrud, E. & Uvsløkk, S. 2013. Smart vapour barriers in unventilated wooden roofs in a
Nordic climate – laboratory study of drying effect under shaded conditions. Proceedings of 2nd
Central European Symposium on Building Physics, September 9-11, 327-334.
Icopal, 2005. Hygrodiode ©. Product Brochure. www.icopal.dk.
Korsgaard, V. 1985. Hygro diode membrane: A new vapor barrier. ASHRAE/DOE/BTECC
Conference: Thermal Performance of the Exterior Envelopes, Clearwater Beach, Florida,
December 2-5.
Künzel, H.M. 1996. Humidity controlled vapour retarders reduce risk of moisture damages.
Proceedings of the 4th Symposium on Building Physics in the Nordic Countries, Espoo, Finland,
September 9-10.,447-454.
WUFI. 2013. http://www.wufi-pro.com
KEYWORDS: Exterior air barrier, HAM-simulation, DELPHIN 5, air flow, natural convection
SUMMARY:
This paper studies the hygrothermal performance of timber frame wall configurations with exterior air
barrier systems in various European climates. First, the impact of different insulation materials will
be discussed for Belgium climate conditions. The simulation results demonstrate that the application
of mineral wool insulated timber frame walls in combination with exterior air barriers result in
increased moisture loads. Moreover, unavoidable (small) air channels between the mineral wool and
the adjacent exterior air barrier significantly increase natural convection and add up to harmful
moisture levels. Yet the use of cellulose insulation can avoid these issues. The second part of the paper
studies the performance of these wall elements for 12 European climates, focussing on Nordic
countries in specific. The simulation results indicate that the technique of exterior air barrier is more
suitable for continental climates rather than for cold and moderate sea climates in Europe.
1. Introduction
The execution of exterior air barrier systems in timber frame construction has gained increasing
interest in the aim to meet the severe airtightness requirements of Passive houses (n50<0.6 ACH). In
practice, fewer joints and perforations of exterior air barriers are an advantage compared to the
traditional interior air barrier systems (Langmans, 2010). Also in Norway the practical advantages of
exterior air barriers are mentioned (Myhre & Aurlien, 2005 and Holøs & Relander, 2010). Yet for cold
and moderate climates, such as the Northern European countries, the hygrothermal impact of exterior
air barriers remains unstudied. Consequently, the execution of wall systems with an exterior air barrier
only is currently not recommended for these climate types.
The present paper performs numerical heat, air and moisture (HAM)-simulations on light weight walls
with exterior air barrier systems to explore the risks involved. Yearly simulations under realistic
climate conditions were conducted with an adjusted version of DELPHIN 5, capable to model forced
and natural convection in interaction with detailed heat and moisture transport ( Nicolai, 2007 and
Langmans et al. 2012). The model applies a two-domain approach; (1) porous building materials with
(2) adjacent air channels assuming fully-developed laminar flow. The applied model is extensively
documented and evaluated with experimental results in previous work (Langmans, 2013).
The current article performs numerical simulations to investigate the hygric response of exterior air
barrier systems in timber frame construction in a European climates. In a first step, interest lies in the
impact of different insulation materials. Thereafter, simulation are conducted for various European
climate conditions. Herein, focus is on the performance in Nordic countries. The transient simulations
were conducted for yearly climate conditions and the mould index (M) (Hukka & Viitanen, 1999) and
accumulated interstitial condensation (AC) were applied to assess the performance of these timber
frame elements (Mmax < 3 and ACmax < 0.1 kg/m²).
The following section will discuss the used wall configurations and climate conditions.
(1)
in which V (m³) corresponds to the volume of the room and HIR* (kg/m3%RV) to its hygric inertia, n
(1/h) is the ventilation rate and Gvp(kg/s) is the vapour production in the room. The used parameters
for the single zone model are chosen rather conservative: a small volume of 50 m3 with a high
moisture load of two active persons (120 gram/h) between 8h-22h. The nominal ventilation rate in
living spaces is 3.6m3/m2/h according to the Belgium standard, corresponding to a ventilation rate of
1.51/h for this room. Several studies, however, indicated that the actual ventilation rate is often much
lower in reality. Yet to include the effect of a realistic ventilation system in operation, the nominal
ventilation rate was reduced to 0.5 1/h. Finally, a value of 1.5gram/m3%RV was chosen for the hygric
inertia.
3. Simulation results
3.1 Insulation material
First, the impact of the insulation material is investigated. Herein, two commonly applied groups of
insulation materials are selected for the numerical investigation: (1) mineral wool and (2) cellulose
insulation. Most often mineral wool insulation with a density of 20 kg/m3 (MW20) is used in Belgium.
Densities of 30 kg/m3 (MW30) have only a limited market as they are more expensive. Nevertheless,
both densities (MW20 and MW30) are included in the investigation here. For inflated loose-fill
cellulose insulation, however, densities of around 60 kg/m3, commonly applied in practice, are
adopted here. In addition, cellulose insulations with lower densities (40 kg/m3 and 50 kg/m3) are
examined as well. Table 2 presents the maximum mould index and accumulated condensate on the
exterior air barrier (BIFB) for the reference configuration (Fig.1) with the discussed insulation
materials. The predicted maximum mould index never exceeds the proposed limit of 3 and no
condensation occurred for this type of insulation for the wall elements insulated with cellulose. For the
mineral wool insulation, however, mould growth problems are predicted for the standard density of 20
kg/m3. A mould index of 4.3 was found for this density, corresponding to a visual mould covering
percentage between 10-50%. Mineral wool with increased density (30 kg/m3), in contrast, meets the
proposed moisture limit state. Finally, this table shows that no condensation on the BIFB exterior air
barrier occurred for both the mineral wool and cellulose.
The results of Table 2 correspond to an ’ideal’ installation of the insulation layer, in that, perfect
contact between the insulation and the adjacent layers is assumed. On the contrary, several
experimental studies proved the importance of small discontinuities along the interface of mineral
wool insulation layers (Brown et al., 1993). Such small air cavities are responsible for increased levels
of natural convection and consequently higher moisture loads. As a result the presence of small air
cavities around mineral wool insulation should be considered. Therefore, the simulations with mineral
wool insulated elements were repeated, introducing imperfections along the interface for the insulation
and the adjacent layers. Herein, five configurations are considered, varying the position of these air
cavities: an air channel along the (a) cold side of the insulation, (b) on top of the insulation layer (c)
both cold and warm sides, (d) both the cold and warm side and at the top of the insulation, and (e) at
all four sides of the insulation. For these configurations the impact of channels of 1 mm, 3 mm and 5
mm have been studied. All combinations are simulated for both mineral wool with a density of 20
kg/m3 and 30 kg/m3 for which the corresponding simulation results are summarised in Table 3. This
table confirms that the risk for moisture problems is highly influenced by the presence of small air
channels along the mineral wool layer. Even a small air channel of 1mm at both sides of the insulation
increases the mould index from 2.3 to 3.7 for MW30. As a consequence, the above-mentioned
conclusion that higher mineral wool densities decrease the risk for moisture problems only holds for
perfect contact conditions. From the moment imperfections around the mineral wool insulation layer
(which are unavoidable in practice) are considered, the density of the mineral wool becomes of minor
importance. In addition, Table 3 shows that most critical moisture limit state again corresponds to
mould growth.
TABLE 3:Iimpact of air cavities (1mm, 3mm and 5mm) around insulation layer on the maximum
mould index M (-) and maximum accumulated condensation MC (kg/m²) for mineral wool.
MC (kg/m²)
M (-)
Density (kg/m³) 20 30 20 30 20 30 20 30 20 30
1 mm 4.5/- 2.8/- 4.5/- 2.9/- 4.6/- 3.7/- 4.8/- 4.1/- 4.8/- 4.1/-
3 mm 5.5/- 4.6/- 4.7/- 3.6/- 5.8/- 5.1/- 6/1.1 6/0.6 5.9/1.4 6/1.2
5 mm 5.9/- 4.9/- 4.7/- 3.7/- 6/0.1 5.6/- 6/>5 6/4.5 6/>5 6/>5
Fig. 2 depicts the hourly condensation along the height of the exterior air barrier for MW20 with air
channels of 5mm at both sides. This graph demonstrates that condensation clearly occurs during
winter periods and is concentrated at the upper height of this layer, indicating the importance of
natural convection within the component.
FIG. 3: Overview of the European locations for which simulations are performed.
First, the averaged outdoor winter temperatures (θwinter), provides a direct indication of the potential for
buoyancy driven air flow within the component, and thus, the increased moisture loads. Second, the
Πy -factor expresses the severity of the exterior climate. This factor (kg/m³), proposed by Hagentoft &
Harderup (1993) for the selection of reference years, calculates the yearly averaged difference between
the vapour concentration in the outdoor air and the maximum moisture concentration on the building
envelope surface on hourly data:
𝑣 𝜃 𝑣 | year (2)
in which vout,sat (θs) represents the saturation vapour concentration at the buildings surface for a North
oriented wall and vout corresponds to the moisture concentration of the outdoor air. As this value is
calculated on a yearly basis, it represents the drying potential of the climate, and is thus, a relevant
FIG. 4: Top: maximum mould index on BIFB and accumulated condensation on foil as a function of
θwinter (left) and Πy (right) for MW30 insulated elements, bottom: maximum mould index on BIFB and
maximum moisture content of cellulose insulation as a function of θwinter (left) and Πy (right) for CL60
insulated elements.
4. Conclusions
This paper investigated the hygric response of a timber frame walls with an exterior air barrier in a
Belgian climate. The simulation results identified a large impact of the applied insulation material on
the hygrothermal performance of this component. No moisture problems were found for the walls
insulated with loose-fill cellulose insulation. For mineral wool (glass wool) insulated walls, however,
excessive mould growth levels were obtained on the upper parts of the exterior air barrier as a result of
natural convection within the insulation layer. Furthermore, the simulations emphasised the
importance of small air channels along the mineral wool and the adjacent layers. Such limited
deficiencies increased the risk for mould growth and harmful condensation amounts against the
exterior air barrier. The second part of the paper verified the performance of highly insulated walls
with an exterior air barrier in 12 European climates. From the limited amount of locations tested, it
followed that Northern and Western Europe regions were most vulnerable to mould problems. No
mould problems were found for the (dryer) continental climates, such as Berlin and Vienna. For
interstitial condensation on spunbonded foils, however, the decisive climate parameter appeared to be
the winter temperatures. Excessive amounts of interstitial condensation were found for the coldest
climates, such as Helsinki, Oslo, Trondheim and Stockholm.
5. Acknowledgements
Research funded by a Ph.D. grant (grant number 81153) of the Institute for the Promotion of
Innovation through Science and Technology in Flanders (IWT-Vlaanderen).
References
Brown, W. C., Bomberg, M. T., Ullett, J. M., & Rasmussen, J. (1993). Measured thermal resistance of
frame walls with defects in the installation of mineral fibre insulation. Journal of Building Physics,
16(4), 318–339.
Hukka, A & Viitanen, H (1999), ‘A mathematical model of mould growth on wooden material’, Wood
Science and Technology 33(6),475–485.
Holøs, S. B., & Relander, T.-O. (2010). Airtightness Measurements of Wood Frame Low Energy Row
Houses. In BEST conference (pp. 1–11). Portland.
Langmans, J., Klein, R., De Paepe, M. & Roels, S. (2010), ‘Potential of wind barriers to assure
airtightness of wood-frame low energy constructions’, Energy and Buildings 42(12), 2376–2385.
Langmans, J., Andreas Nicolai, Ralf Klein, and Staf Roels. (2012). “A quasi-steady state
implementation of air convection in a transient heat and moisture building component model.”
Building and Environment 58 (12) 208-218.
KEYWORDS: building envelope, hollow brick, computational simulation, heat and moisture
transport, hygrothermal performance
SUMMARY:
Computational analysis of hygrothermal performance of several types of hollow-brick based building
envelopes is presented, aimed at the optimal hygric properties of hollow brick and exterior plaster. In
the simulations, Künzel’s mathematical model of coupled heat and moisture transport is used.
Dynamic climatic data in a form of Test Reference Year for Prague are utilized as boundary
conditions. The computational results show that for a suitable hygrothermal performance of the
analyzed envelopes the key factor is the hygric properties of the exterior plaster, whose value of
moisture diffusivity has to be reduced to at least 1 · 10-10 m2/s, e.g. using hydrophobic agents. In that
case, moisture content can be held under the hygroscopic threshold for almost the whole reference
year, which ensures the elimination of the presence of liquid moisture inside the envelope.
1. Introduction
Hollow brick blocks have found a widespread use in building industry during the last decades. The
increasing requirements to the thermal insulation properties of building envelopes (Uygunoglu and
Kecebas 2011) given by national standards in Europe led the brick producers to reduce drastically the
production of common solid bricks. The brick blocks with complex systems of internal cavities (da
Silva Almeida et al. 2013, Alhazmy 2010, Arendt et al. 2011) replaced the traditional bricks and
became dominant on the building ceramics market. However, the development of the new types of
bricks was very fast and their properties are often not known with a sufficient accuracy. Consequently,
in the current durability assessment of hollow brick-based building envelope systems a complex view
is often missing. Precise and serious analyses of hygrothermal performance based on sound scientific
knowledge are not performed very frequently. However, hollow brick-based building envelopes should
be designed as a system consisting of hollow bricks, internal and external finishes and possibly also
thermal insulation layer, and in this design process, the details have to be solved using suitable
methods. In this paper, influence of the effective hygric parameters of hollow brick and exterior plaster
on the hygrothermal performance of building envelope is analyzed, in order to identify possible
problems related to the coupled heat and moisture transport. Main objective of this paper is to find
such combination of hygric parameters of exterior plaster and hollow brick, which will reduce the
moisture content in the envelope as much as possible.
Thermal Lime-metakaolin
insulating plaster plaster
Hollow brick
30 500 10
5.0E-10
Moisture diffusivity (m2/s)
4.5E-10
4.0E-10
3.5E-10
3.0E-10
2.5E-10
2.0E-10
1.5E-10
1.0E-10
5.0E-11
0.0E+00
0.00 0.05 0.10 0.15 0.20 0.25
Moisture content (m3/m3)
dH ∂T
dT ∂t
[ ]
= div(λgradT ) + Lv div δ p grad (ϕp s ) , (4)
1 · 10-8 A1 A2 A3
Moisture
diffusivity of
1 · 10-9 B1 B2 B3
thermal insulating
plaster (m2/s)
1 · 10-10 C1 C2 C3
The results of hygric and thermal behaviour of studied building envelopes are presented in the form of
temperature and relative humidity profiles for November 22 of the reference year. This day appeared
statistically as one of the most critical from the point of view of hygric performance. According to the
reference climatic data for Prague, a 3-day period with highly increased relative humidity of the air (>
97 %) due to the effect of rainfalls ends in this day. Such high value of the relative humidity
corresponds to the overhygroscopic (liquid) moisture content, thus the presence of liquid moisture in
the pore space of exposed building materials is to be expected. Although in the subsequent days the
relative humidity of the air is still very high, it remains in the hygroscopic range and the liquid
100
A1
A3
95
90
Relative humidity [%]
85
80
75
70
65
60
55
0 0.1 0.2 0.3 0.4 0.5
Distance [m]
100
A3
C3
95
90
Relative humidity [%]
85
80
75
70
65
60
55
0 0.1 0.2 0.3 0.4 0.5
Distance [m]
4. Conclusions
The effect of hygric properties of hollow brick and external plaster on the hygrothermal performance
of hollow-brick based envelopes was analyzed in the paper. According to the results, hygric
parameters of hollow brick did not play as significant role as the hygric parameters of the external
plaster. This confirmed that the properties of surface layers of building envelopes are very important
from the point of view of hygrothermal performance of whole building envelopes. In order to reduce
moisture intake, it is necessary to decrease the value of moisture diffusivity of external plaster. Among
all investigated variations, the apparent moisture diffusivity of external plaster equal to 1 · 10-10 m2/s
provided the best results. This value can be easily reached using external plasters with hydrophobic
admixtures.
5. Acknowledgements
This research has been supported by the Ministry of Industry and Trade of the Czech Republic, under
project No FR-TI3/085.
References
Alhazmy M.M. 2010. Numerical investigation on using inclined partitions to reduce natural
convection inside the cavities of hollow bricks, International Journal of Thermal Sciences 49 2201-
2210.
SUMMARY:
Transient hygrothermal simulations of two-dimensional detailed wall constructions set high demand
on computational capacity. When using implicit time integration combined with spatial Finite-Volume
discretization the resulting system of equations contains a large number of unknowns and the
Jacobian matrices are sparse. This motivates the application of iterative solvers, particularly Krylov
Subspace methods, to the resulting linear equation system. We compare two representative
algorithms for direct and iterative linear equation system solvers: the Gaussian elimination with LU-
decomposition versus the ILU-preconditioned GMRES. The numerical study covers three classes of
problem types that may occur in the hygrothermal simulation: diffusion-dominated transport
problems with highly nonlinear coupling of heat and moisture, diffusion-dominated transport inside
constructions with extremely varying material properties, mixed convection-diffusion hygrothermal
transport problems. The aim of the study is to evaluate the applicability of Krylov Subspace methods
to different problem types and to provide advice for the choice of the linear equation system solver.
1. Introduction
Building component simulation considers the dynamic thermal and hygric transport inside wall
constructions. We use the transport equations based on porous media theory by (Grunewald 1997) that
include both liquid water convection and vapor diffusion. The transport equations are highly
nonlinear and stiff. Therefore, only implicit time integration methods combined with established
nonlinear solution methods may solve these equations sufficiently stable and efficient.
The numerical solution requires a time and spatial discretization. When using implicit time integration
methods, many linear equation systems have to be solved. With increasing complexity of the
construction the solution of the linear equation system dominates the simulation performance
(Vogelsang 2013, Nicolai 2014). Therefore, the demand for an acceptable simulation time limits the
problem dimension. Consequently, efficient numerical solution methods are of high interest.
The discretization of partial differential equations (PDE) results in matrices with characteristic
pattern. This fact offers a high potential for optimizing the linear equation solver. The classical
approach focuses on banded matrices. Alternatively the matrix can be considered as sparse with most
elements of the matrix equal to zero. These matrix properties motivate the application of sparse matrix
linear solvers.
3. Simulation examples
Typical for practical applications are two-dimensional simulation problems with a high geometrical
complexity. Sufficiently refined spatial discretization is necessary and leads to critical demand on
computational capacity. Further, for hygrothermal transport problems the intensity of nonlinear
interactions varies highly for each simulation case. Consequently, we take into consideration different
levels of physical model complexity for the current study.
We start with a thermal simulation, continue with a moderately nonlinear hygrothermal simulation,
and finish with a discussion of a highly nonlinear hygrothermal simulation of a steel detail in the
presence of real climate conditions.
4. Conclusions
The study considers three simulation cases that appear in practical applications. Requirements on
sufficient accuracy lead to high grid detail and critical simulation times. Choosing suitable numerical
methods may drastically reduce the computational effort and therefore give access to simulation
problems with an increased complexity.
In the case of a linear problem (case 1) the Band Solver combined with a Modified Newton method
proves to be a high-performance method. However, hygrothermal simulations are of nonlinear type.
We have shown that Krylov Subspace methods prove to be a very efficient alternative to direct banded
linear equation system solvers when applied to nonlinear transport problems. As the equation system
setup does not demand significant simulation time the Modified Newton strategy is no longer
necessary. We suggest the combination of the GMRES method with a Standard Newton procedure for
large nonlinear problems.
However, the GMRES method is only applicable when problem-specific preconditioned. In our tests
the standard ILU preconditioner provided a base-level for convergence acceleration. The development
of suitable preconditioning strategies will be subject of the future work.
References
Thermal bridges in building construction - Heat flows and surface temperatures - Detailed
calculations (ISO 10211:2007). German version EN ISO 10211:2007
Grunewald J. 1997. Diffusiver und konvektiver Stoff- und Energietransport. Ph.D thesis. Dresden
University of Technology.
Hindmarsh A. C. et al. 2005. SUNDIALS: Suite of Nonlinear and Differential/Algebraic Equation
Solvers. ACM Transactions on Mathematical Software. Vol. 31(3). pp. 363-396.
Meister A. 2008. Numerik linearer Gleichungssysteme. Vieweg & Sohn Verlag, Wiesbaden.
Nicolai A. 2008. Modeling and Numerical Simulation of Salt Transport and Phase Transitions in
Unsaturated Porous Building Materials, Ph.D thesis, Dresden University of Technology.
Nicolai A. 2014. Performance improvement of HAM simulations through an optimized grid
numbering technique, 10th Nordic Symposium on Building Physics. Lund.
Paepcke A. 2014. Building constructions for the performance test of dynamic hygrothermal transport
solvers arising from practical application. In prep.
Saad Y. 2003. Iterative Methods for Sparse Linear Systems, SIAM Society for Industrial and Applied
Mathematics. Philadelphia.
Vogelsang S. 2013. Parallel Hardware Architectures. EnTool 2013 – Symposiuum Dresden.
SUMMARY:
The certification of building energy performance according to the European Building Performance
Directive is in many EU countries based on monthly quasi-steady-state calculation of the expected
energy use. But does this calculated energy use correspond to the actual energy use? The aim of this
paper is to contribute to the minimisation of the performance gap between calculated and actual
energy use by improving the precision of the calculated energy use. It is common practice to use the
declared U-value of a glazing in calculations of the expected energy use despite the fact that this U-
value is varying with outdoor temperature, wind speed and direction. Data presented in this paper
suggests that the dynamic changes in U-value due to weather conditions should be taken into account
to obtain a more accurate calculation of the expected energy use of a building in the Nordic
countries. Consequently, three alternative approaches to minimise the performance gap is suggested:
1) abandon the monthly quasi-steady-state method for certification in favour of a dynamic method, 2)
use the approach suggested in this paper to reduce the performance gap and then visualise the
remaining performance gap as an uncertainty of the overall calculation result, or 3) visualise the full
uncertainty of the overall calculation result due to the performance gap.
1. Introduction
The European Building Performance Directive (EPBD 2010) has since 2006 been the framework for
certification of building energy performance in European Union (EU) member states. The certification
process varies from country to country but is very often based on a calculation of the expected energy
use (Lausten et al. 2010). Various stakeholders in the building industry use this calculated energy use
for estimating operating budgets, investment security, policy making etc. Consequently, an increasing
number of stakeholders are raising the question: does the calculated energy use correspond to the
actual energy use?
EPBD prescribes the use of a standardised energy calculation method for calculating the expected
energy use. ISO 13790 (2008) describes two basic types of methods for this purpose where the
monthly quasi-steady-state method is the most widespread method in the EU member states. This
method has its benefits but it is also in many aspects a simplified representation of reality. In relation
to the above raised question, it is therefore relevant to investigate the significance of these simplified
representations. This paper focuses on the simplified representation of heat loss through glazings.
The monthly quasi-steady-state energy calculation relies on the so-called declared U-values of
glazings in its calculation of the annual heat loss. The purpose of the declared value is product
comparison and it is obtained by calculating the U-value according to EN 673 (1997) using a set of
fixed standardised boundary conditions. But EN 673 (1997) also makes it clear that the U-value of a
glazing is not a static value as it varies with the environmental temperatures on each side of the
window, outdoor wind speed and direction, and indoor air flows at the window. This paper therefore
investigate the theoretical effect on building energy calculations when taking into account the
2. Method
This section describes a procedure for calculating the dynamic U-value of a glazing due to changes in
hourly weather conditions. The procedure is used in two types of energy calculations:
1. Simplified dynamic calculation.
The procedure is integrated in a simplified hourly-based tool called iDbuild (Petersen &
Svendsen 2010). The purpose is to investigate the theoretical effect of a dynamic U-value on
energy performance calculations compared to the use of the static declared U-value.
2. Monthly quasi-steady-state calculation.
The procedure is used to calculate the mean monthly U-values. These values are used for
calculating the energy use of an office space in the program Be10 (SBi 2011). Be10 is the
Danish implementation of the monthly quasi-steady-state energy calculations method
described in ISO 13790 (2008). The purpose is to investigate whether the use of mean
monthly U-values in monthly quasi-steady-state energy calculations is a reasonable approach
to take into account any effect of weather data on glazing U-value identified in the simplified
dynamic calculations.
The results from the calculations and comparisons can be found in section 3.
Where k factor depending on number of panes, k=0.5 for two-layer glazings, k=0.25 for the
cavity closest to outside in a three-layer glazing and k=0.75 for the cavity closest to
inside in a three-layer glazing
These guesses are used to calculate the initial U-value of the glazing. The initial U-value is used to
calculate the temperatures of the glazing pane surfaces facing the gas space (t s,x). ts,x of a certain
surface x can be expressed as:
x
h z
1
t s , x ti ti to z i 1 (6)
U
Where ti indoor air temperature (°C)
to outdoor air temperature (°C)
z index referring to the heat transfer coefficients in equation 1 starting with hi
A proportion of the total solar radiation is absorbed in the panes and thus increases the pane
temperature. For example, an incident solar radiation of 500 W/m2 increases the temperature of an
FIG 1. Hourly U-value of two-layer west-facing glazing in Kiruna, Sweden. The vertical dotted line is
the declared U-value using the boundary conditions in EN 673 (1997).
hi
t s , x 4 t rm 4 (4)
t s , x t rm
3. Results
The following case was used to investigate the impact of dynamic glazing U-value on energy need for
heating and cooling using the formulas specified in section 2 compared to the use of the static
declared U-value. Assume a single-sided, single-zone two-person office with a west-facing window.
FIG 3. Primary energy use for heating and cooling, respectively, for the three-layer glazing at the
four locations. Dynamic energy calculations.
U T
is
i i
U month n
(6)
T
is
ref Tout,i
Where s the annual hour number in the start of a certain month (-)
n the annual hour number in the end of a certain month (-)
ref a reference indoor temperature (°C), e.g. 20 °C
Results are depicted in FIG 4 and FIG 5. The heating demand for both glazing solutions shows the
same tendency in deviations as the dynamic calculation in FIG 1 but not the same relative or absolute
magnitude. The tendency for the cooling demand in the dynamic calculation is not repeated in the
quasi-steady-state calculation.
FIG 5. Primary energy use for heating and cooling, respectively, for the three-layer glazing at the
four locations. Monthly mean quasi-steady-state calculations.
4. Conclusions
The aim of this paper is to minimise the performance gap between calculated and actual energy use by
improving the precision of the calculated energy use. Focus is on the gap caused by the use of the
declared U-value of a glazing as a static all-year input in energy calculations instead of the U-value as
a function of hourly changes in weather conditions.
Data from dynamic hourly building simulations suggest that the dynamic changes in U-value due to
weather conditions should be taken into account to obtain a more accurate calculation of the expected
energy use of a building in the Nordic countries. The use of declared U-value in performance
simulation of two-layer glazing solutions leads to an underestimation of the heating demand for all
investigated locations and an overestimation of the cooling demand. For a three-layer glazing solution
the use of dynamic U-value results in less heating demand whereas the cooling demand is practically
the same. The use of declared U-values instead of dynamic U-values in energy performance
calculations may therefore be a reason for performance gaps between calculated and actual energy
use. Based on this data, the paper investigates whether the use of monthly mean temperature-weighted
U-values is a reasonable approach to represent this effect when calculating the expected energy use
according to the monthly quasi-steady-state method in ISO 13790. The conclusion is, assuming that
the dynamic hourly calculation represents actual conditions, that the suggested approach is not
sufficient to represent the identified effect even though the approach does compensate somewhat for
the performance gap.
Three alternative approaches to minimise the performance gap is hereby suggested: 1) abandon the
monthly quasi-steady-state method for certification in favour of a dynamic method, 2) use the
approach suggested in this paper to reduce the performance gap and then visualise the remaining
performance gap as an uncertainty of the overall calculation result, or 3) visualise the full uncertainty
of the overall calculation result due to the performance gap.
References
Cooper P.I. 1981. The effect of inclination on the heat loss from flat-plate solar collectors. Solar
Energy 27 (5) 413-420.
EN 673. 1997. Glass in building – Determination of thermal transmittance (U value) – Calculation
method. European Committee for Standarization. Brussels Belgium
EPBD. 2010. Directive 2010/31/EU of the European Parliament and of the Council of 19 December
2010 on the energy consumption of buildings (recast).
ISO 13790. 2008. Energy performance of buildings – calculation of energy use for space heating and
cooling. International Organization for Standardization. Geneva Switzerland
ISO 15099. 2003. Thermal performance of windows, doors and shading devices – detailed
calculations. International Organization for Standardization. Geneva Switzerland
Mirsadeghi M., Cóstola D., Blocken B. & Hensen J.L.M. 2013. Review of external convective heat
transfer coefficient models in building energy simulation programs: Implementation and
uncertainty. Applied Thermal Engineering 56 (1–2) 134-151.
Ito N., Kimura K., & Oka J. 1972. A field experimental study on the convective heat transfer
coefficient on exterior surface of a building. ASHRAE Transactions 78 (1) 184-191
Lausten J., Brophy V. & Ryan L. 2010. Energy performance certification of buildings. IEA Energy
Conservation in Buildings and Community Systems Annex 53 – Total energy use in buildings:
analysis and evaluation methods. Paris France
Palyvos J.A. 2008. A survey of wind convection coefficient correlations for building envelope energy
systems' modelling. Applied Thermal Engineering 28 (8-9) 801-808.
Petersen S. & Svendsen S. 2010. Method and simulation program informed decisions in the early
stages of building design. Energy and Buildings 42 (7), 1113–1119.
Rowley F.B., Algren A.B. & Blackshaw J.L. 1930. Surface conductance as affected by air velocity,
temperature and character of surface, ASHRAE Transactions 36 p. 429.
SBi 213. 2011. Bygningers energibehov. 5th ed. Copenhagen, Statens Byggeforskningsinstitut. 120 p.
U.S. Department of Energy. 2013. EnergyPlus Energy Simulations Software website,
www.eere.energy.gov/buildings/energyplus/
KEYWORDS: vertical greening, green walls, living walls, moisture conditions, evaporative cooling,
energy balance, mould growth.
SUMMARY:
Vertical greening on external walls has lately been the subject of increasing interest, where much
hope is spent upon the benefits of enhanced energy performance of buildings. However, there is a lack
of knowledge and experiences of the technical performance of such walls, especially their hygro-
thermal and energy related behaviour under different climatic conditions.
This paper presents a long-term, full-scale field experiment study on a masonry wall in the southern
part of Sweden. Different kinds of wall greening have been included in the experiment, mainly
solutions based on commercially available cassettes, simple home-built pocket-like plant pots and
different climbers.
Monitored parameters in the technical part of the program include outdoor and indoor climate, air
gap and surface temperature and relative humidity, driving rain and thermal flux.
The wall greening systems studied quite strongly influence the hygro-thermal behaviour of the original
walls; greening may even contribute to dryer conditions by protecting from driving rain. Caution must
be paid to protect from water leakage, which may cause mould growth. The energy balance of poorly
insulated walls is slightly improved in winter by raised air gap temperatures, and in summer by
shading. No evaporative cooling effect has been found in this study.
1. Introduction
Vertical greening on external walls has during the last years been the subject of increased interest from
city planners, architects and real estate developers (Köhler 2008). Much hope is spent upon the
benefits of green walls, not only from an esthetical point of view, but also technically: e.g. enhanced
energy performance with vertical greening, cooler street canyon climate during hot days and traffic
noise reduction. Vertical greening concepts have been developed in a number of countries during the
last decades, especially in Germany, Austria, the Netherlands and the US (Perini et al. 2013, Enzi
2011, Mir 2011, Tilley et al. 2012). In the northern parts of Europe and North America there exists
however much less practice and thus experience of exterior wall greening and so little is known today
about different greening systems and adequacy of different plants as well as the technical behaviour of
wall greening in relation to the building under influence of e.g. more severe winter conditions.
In a collaborative research project with the Swedish University of Agricultural Sciences (SLU),
Malmö University and the building developer and contractor PEAB as partners, research on the hygro-
thermal and energy related performance forms one part which is performed by Malmö University and
CEC Design AB. The results from this part of the project form the basis for this paper. The most
extensive part of the project, however, is related to horticultural issues, such as planting systems, plant
choice and adequacy, irrigation solutions etc. This research is mainly performed by SLU.
3. Experimental setup
The study is using two experimental sites, both located in an industrial area in the city of Malmö in
southern Sweden (N 55.6108, E 12.9896), see FIG 1.
FIG 2. Site 1 with driving rain meters. FIG 3. The wall-greening systems (Site 2).
Left: Humulus lupulus (Hops) Left: Felt pocket module
Right: Fallopia baldschuanica (Bukhara fleeceflower) Right: Modular panel
14 4
12 3
10 2
8 1
6 0
4 -‐1
2 -‐2
0 -‐3
-‐2 -‐4
00:00 06:00 12:00 18:00 00:00 06:00 12:00 18:00
K Climber Pockets RefArea Module
(NoIns) Module
(Insul) WindSpeed m/s
16 5
14 4
12 3
10 2
8 1
6 0
4 -‐1
2 -‐2
0 -‐3
-‐2 -‐4
00:00 06:00 12:00 18:00 00:00 06:00 12:00 18:00
FIG 5. Air temperature differences in the air spaces or close to the wall compared to outdoor air
during; 17–18 January (top); and 17–18 July (bottom).
It is often argued that vegetated walls and roofs generate a considerable evaporative cooling effect that
reduces the need for active cooling of the building. The backside temperatures of the systems were not
measured, but in FIG 6 we see the differences between wall surface and air gap temperatures of the
test areas. A negative balance could thus indicate a radiative cooling effect from the green systems.
However, this effect seems to be very marginal. Probably most of the evaporative heat exchange will
take place between the plants and the surrounding air outside the plants, thus cooling the air in e.g. a
street canyon. Later in the project this will be studied by means of simulations.
K Climber Pocket RefArea Module Module-‐Insul
6
4
Difference:
Tsurf
— Tair
3
0 Missing
data
-‐1
okt-‐12 nov-‐12 dec-‐12 jan-‐13 feb-‐13 mar-‐13 apr-‐13 maj-‐13 jun-‐13 jul-‐13 aug-‐13 sep-‐13 okt-‐13
FIG 6. Monthly average temperature differences between wall surfaces and air gaps.
1
00:00 06:00 12:00 18:00 00:00 06:00 12:00 18:00
g/m³ Climber Pockets RefArea Module
(NoIns) Module
(Insul)
15
12,1 13,4 11,8 12,2 12,0
Above: Monthly
average
(g/m³)
14
13
12
11
00:00 06:00 12:00 18:00 00:00 06:00 12:00 18:00
FIG 7. Vapour concentrations in the air spaces or close to the wall during;
17–18 January (top); and 17–18 July (bottom). Monthly averages also included.
There are two reasons for the higher moisture level behind the pockets; 1) the above mentioned steel
sheet cover on top of the space reduces ventilation; 2) at the same time water transport through the felt
structure into the space was discovered. Conditions suitable for mould growth occur from time to time
and a preliminary investigation performed with an endoscope seems to confirm the existence of cob
web and signs of mould growth. High moisture levels found in the wall surface could also indicate
water leakage from the irrigation system. These issues will be further investigated.
0
okt-‐12 nov-‐12 dec-‐12 jan-‐13 feb-‐13 mar-‐13 apr-‐13 maj-‐13 jun-‐13 jul-‐13 aug-‐13 sep-‐13 okt-‐13
FIG 9. Driving rain. Two climbers compared to bare reference area on a wall facing prevailing wind.
6. Acknowledgements
The Swedish Research Council FORMAS and the building developer and contractor PEAB have
financially supported the research study. The support is gratefully acknowledged.
References
Enzi V. et al. 2011. Leitfaden Fassadenbegrünung. (Guide book for vertical greening.) ÖkoKauf Wien,
Arbeitsgruppe 25 Grün- und Freiräume. www.oekokauf.wien.at
Köhler M. 2008. Green facades – a view back and some visions. Urban Ecosystems, 2008, Vol.
11:423-436. Springer.
Mir M.A. 2011. Green facades and building structures. Master thesis CIE5060, TU Delft.
Nevander L.E. & Elmarsson B. 1994. Fukthandbok, Praktik och teori. (Moisture handbook, Practise
and theory) 2nd ed. Stockholm, AB Svensk Byggtjänst and the authors. 538 p.
Perini K. et al. 2012. Vertical greening systems, a process tree for green facades. Urban Ecosystems,
2012, Vol. 16:265-277. Springer.
Sandin K. 1987. Fukttillstånd i autoklaverade lättbetongväggar. (The moisture condition in aerated
lightweight concrete walls.) Report TVBM-3026. Div. of Building Materials. Lund Institute of
Technology, Lund, Sweden
Tilley D. et al. 2012. Vegetated walls: Thermal and growth properties of structured green facades.
Final report to Green Roofs for Healthy Cities – Green Walls Group. UM-09040836. Toronto.
www.greenroofs.org
KEYWORDS: actual energy performance, in situ characterisation, dynamic data analysis, grey box
models
SUMMARY:
Several studies have shown that actual thermal performance of buildings after construction may
deviate significantly from that anticipated at design stage. As a result, there is growing interest in full
scale testing of components and whole buildings. The IEA EBC Annex 58-project ‘Reliable Building
Energy Performance Characterisation Based on Full Scale Dynamic Measurements’ is developing the
necessary knowledge, tools and networks to achieve reliable in-situ dynamic testing and data analysis
methods that can be used to characterise the actual thermal performance and energy efficiency of
building components and whole buildings. The research within this project is driven by case studies.
As a first simple case, an experiment on testing and data analysis is performed on a round robin test
box. This test box can be seen as a scale model of a building, built by one of the participants, with
fabric properties unknown to all other participants. Full scale measurements have been performed on
the test box in different countries under real climatic conditions. The obtained dynamic data are
distributed to all participants who have to try to characterise the thermal performance of the test
box’s fabric based on the provided data.
This paper describes the rationale, aim and objectives of IEA EBC Annex 58 and presents more in
depth the first results obtained on the round robin experiment. It is shown how different techniques
can be used to characterise the thermal performance of the test box, ranging from a simple stationary
analysis to advanced dynamic data analysis methods.
1. Introduction
The rise of living standards, the scarcity of natural resources and the awareness of climate change
resulted in an international pressure to significantly reduce the energy consumption of buildings and
communities. In several countries more stringent requirements are imposed by energy performance
legislation and also an increased awareness for environmental issues in building codes can be noticed.
Mostly, requirements and labelling of the energy performances of buildings is done in the design
phase by calculating the theoretical energy use. Several studies showed however that the actual
performance after realisation of the building may deviate significantly from this theoretically designed
performance. Part of the deviations can be explained by user behaviour, but the other part has to be
attributed to the physical features of the building and its systems. For the latter, building performance
characterisation based on full scale testing – testing of building components or whole buildings under
FIG 1. Left: infrared pictures of the outer leaf of two full scale highly insulated cavity walls in the
VLIET-test building at KU Leuven. The impact of workmanship on the thermal performances of the
walls is clearly visible. Right: Measured versus predicted whole house heat losses (W/K) for 18 new
build dwellings in the UK. None of the houses is able to reach the designed values.
Examples as those mentioned above, explain why at present several in situ testing activities are going
on. A recent international workshop showed the interest for full scale testing from all over the world
(Janssens et al. 2011). A growing activity is observed in both full scale testing on building components
(as e.g. in Paslink-cells or in situ on components of real buildings) and on whole buildings (to
characterise thermal performance and energy efficiency of either test buildings or real buildings). So it
is clear that, contrary to what was expected, the numerical building component and building energy
simulation models did not make full scale testing of building (components) redundant. On the
contrary, together with an increased application of numerical simulations, a renewed interest in full
scale testing can be observed. This is not so strange, because dynamic full scale testing showed not
only to be of interest to study building (component) performances under different real conditions – and
as illustrated, quite often a huge difference is observed between predicted and realised performances –,
it is also a valuable and necessary tool to integrate advanced components and systems into simulation
models. Jiménez et al. (2008a) among others used full scale tests on BIPV-façade (BIPV: Building
Integrated Photovoltaic Cells) to develop a so-called grey box-model. A grey box model is based on a
combination of prior physical knowledge and statistics by identifying the unknown parameters of the
system with dynamic data analysis. Once identified, the grey box model is able to predict the thermal
dynamic response of ventilated photovoltaic double skin facades under different climatic conditions.
This way it can be ensured that the behaviour of new advanced building components is integrated in a
correct way in building energy simulation (BES) models. A similar approach of parameter
identification based on dynamic measurements can be used to identify suitable models to describe the
thermal dynamics of whole buildings including building systems (Bacher and Madsen 2011).
Characterising the dynamic behaviour of buildings is an essential and very valuable input e.g. when
optimising energy grids for building communities.
ST1.
State
of
the
art
on
full
scale
testing
and
dynamic
data
Use
of
results
Data
analysis
ST
1.
State
of
the
art
on
full
scale
testing
and
dynamic
data
analysis
analysis
ST
2.
Optimising
full
scale
dynamic
testing
Subtask
2 Subtask
3
ST4.
Application
of
the
developed
framework
Network
of
Excellence ST5.
Setting
up
a
network
of
excellence
Subtask
5
Subtask 1 is a short introductory subtask. Based on a literature review an overview and evaluation is
made of previous and ongoing in situ test activities. An inventory is made both on full scale test
facilities available all over the world as well as on common methods to analyse dynamic data.
Subtask 2 establishes the procedure how to realise a good test environment and test set-up. A decision
Winter 2012-2013 the test box has been tested at the premises of the Belgian Building Research
Institute in Limelette, Belgium (50°41’ N, 4°31’ E). Afterwards the box has been shipped to Spain,
where it was measured under summer conditions in Almeria (37.1° N, 2.4° W). In general, the weather
conditions in Belgium are temperate, with a mild, but rainy, humid and cloudy winter. The weather at
Almeria on the other hand is dry and extremely hot in summer, with large temperature amplitudes
FIG 4.Test box during winter at the measuring site at BBRI. Belgium (left) and during summer at the
Plataforma Solar de Almeria, Spain (right).
At both sites, different experiments have been performed, ranging from co-heating tests with constant
indoor temperature, over free floating temperature runs, to imposed dynamic heating sequences
(ROLBS-signals). During the experiments, heat fluxes on all internal surfaces, together with internal
and external surface temperatures, indoor temperature and delivered heating energy within the box
have been measured. In addition, both test sites are equipped with an outdoor weather station,
measuring all relevant boundary conditions (temperature, relative humidity, wind direction and speed,
diffuse and direct solar radiation, long wave radiation,…). The measured data has been provided to all
participants in the Annex 58-project. They are requested to characterise the thermal performance of
the round robin test box as good as possible based on the provided dynamic data. Both stationary
properties, e.g. the overall heat loss coefficient, and dynamic properties of the test box are aimed for.
TABLE 1. Determined overall heat loss coefficient (W/K) of the round robin test box by different
modelling teams and making use of different data analysis methods
Team Winter data Belgium Summer data Spain
1 Averaging method 3.77-3.92
State space model (RC using LORD) 3.07-3.42
2 Averaging method 2.86-4.15
Linear regression (5’-data) 2.84-4.11
Linear regression (daily averaged data) 3.68-4.12 4.32-4.48
AR(MA)X-models 4.06 4.07-4.32
State space models (RC using LORD) 3.98-4.04 4.23
3 Multiple linear regression (hourly data) 4.77-5.24
Multiple linear regression (daily data) 3.73-4.39
4 State space models 4.27-4.56
5 Linear regression (daily averaged data) 3.99-4.08
State space models (RC using CTSM-R) 3.99
6 State space models (RC using Matlab) 3.97
7 ARX-models 3.95 4.05-4.10
State space models (RC using CTSM-R) 3.84 3.96
8 Averaging method 3.72-3.99
Linear regression (5’-data) 2.98-3.94
AR(MA)X-models 4.01-4.08
State space models (RC using CTSM-R) 4.48
4. Conclusions
The Annex 58-project of the IEA EBC-programme shows that there is currently a larger international
interest in full scale testing and dynamic data analysis. This can be explained by the fact that full scale
testing remains necessary for several reasons. It is for instance the only way to verify our numerical
BES-models. Furthermore, full scale testing allows evaluation and characterisation of the thermal
5. Acknowledgements
The construction of the IEA EBC Annex 58 round robin test box was financially supported by Knauf
Insulation. This as well as the input from the different Annex 58 members who participated in the
common exercises is greatly acknowledged.
References
Bacher P and Madsen H, 2011. Identifying suitable models for the heat dynamics of buildings. Energy
and Buildings, 2011, 43(7): 1511-1522.
Bacher P. and Delff P. 2013. IEA EBC Annex 58 Common Exercise 4: ARX, ARMAX and grey-box
models for thermal performance characterisation of the test box. Internal report.
Bauwens G., Standaert P., Delcuve F. & Roels S. 2012. Reliability of co-heating measurements.
Building Simulation and Optimization Conference 2012, Loughborough, UK. 348-355 p.
Hens H., Janssens A, Depraetere W, Carmeliet J. and Lecompte L. 2007. Brick cavity walls: a
performance analysis based on measurements and simulations. Journal of Building Physics, 31(2):
95-124.
ISO 9869. 1994. Thermal insulation. Building elements. In-situ measurements of thermal resistance
and thermal transmittance (ISO 9869:1994(E)). International Organization for Standardization ISO.
Janssens A., Roels S., Vandaele L. 2011. Full scale test facilities for evaluation of energy and
hygrothermal performances. Proceedings of the International Workshop, Brussels, Belgium, 30-31
March 2011.
M.J. Jimenez and M.R. Heras. 2005. Application of multi-output ARX models for estimation of the u
and g values of building components in outdoor testing. Solar Energy, 79(3):302–310.
M.J. Jiménez, H. Madsen, H. Bloem, B. Dammann. 2008a. Estimation of Non-linear Continuous Time
Models for the Heat Exchange Dynamics of Building Integrated Photovoltaic modules. Energy and
Buildings. 40, pp. 157-167.
M.J. Jimenez, H. Madsen, K. Andersen. 2008b. Identification of the main thermal characteristics of
building components using matlab. Building and Environment, 43(2): 170–180.
Strachan P.A., Baker P.H. 2008. Editorial. Outdoor testing, analysis and modelling of building
components. Building and Environment 43: 127-128.
Wingfield J., Miles-Shenton D., Bell M. 2011. Comparison of Measured versus Predicted Heat Loss
for New Build UK Dwellings – Unpublished Data, Leeds Metropolitan University, Leeds, UK
SUMMARY:
The applicability of the co-heating test for the estimation of a building’s thermal dynamic
characteristics, i.e. the heat loss coefficient, time constants, and effective thermal capacity, is
investigated. Different grey-box models are identified based on a simulated co-heating test while the
heating power supplied to the building is controlled using: (1) a constant indoor air temperature set-
point, (2) multi-sinusoidal signals, and (3) a pseudo random binary sequence signals.
The study shows that measurement data obtained from a co-heating test with the indoor air
temperature maintained constant is not suitable for estimation of the building’s thermal dynamics,
while multi-sine and PRBS based data provide a better basis for model identification. Multi-sine or
PRBS signals covering a larger range of the frequency spectrum provided the best data set for
identification heat dynamics, since the influence of the outdoor environmental conditions is reduced.
While the advantage of the PRBS signal compared to a multi-sine signal is that its frequency spectrum
is continuous, a PRBS is more stringent with respect to the selection of the sampling frequency for
data collection. However, a multi-sine signal demands a well-designed frequency spectrum in order to
excite the building properly.
1. Introduction
The consideration of the heat losses through the building envelope is crucial in order to develop
energy efficient buildings and to reduce carbon emissions. During the last decades, the thermal
insulation of building components has increased a lot, resulting in lower U values. However,
experimental investigations and field studies (Hens et al. 2007) (Wingfield et al. 2009) indicated that
the U values, and hence, the heat loss to the environment, in practice are often much higher compared
to the theoretical U values resulting from calculations (EN ISO 6946). Moreover, the U value in
practice depends strongly on the design and the installation of insulation layers.
A valuable test method to measure the heat loss of a building (both the building fabric and background
ventilation) to the environment is a co-heating test. During the test, the building (often a dwelling) is
heated electrically, using electric resistance point heaters, to an elevated mean internal temperature
(typically 25oC) over a specific period of time, typically 3 winter weeks. The electrical energy
consumption which is required for maintaining the building at the elevated temperature is measured,
resulting in the heat loss coefficient (HLC) by plotting the daily heat input against the daily
temperature difference between the inside and outside air. The development of the co-heating method
started in the late 1970s and early 1980s and has been further developed, resulting in the current
experimental guidelines available in the UK (Wingfield et al. 2011). Moreover, improvements to the
co-heating test for characterization of the dynamic thermal properties of a building are currently
studied in the European project IEA ECBCS Annex 58.
2. Analysis
A case study building is selected and modelled using the TRNSYS software for building energy
performance simulation. The software is used to simulate different co-heating experiments.
FIG 2. Indoor and outdoor air temperature (Ti and Te [oC]), solar radiation (Q [kJ/(hr.m2)]), and
heating power (Q [W]) during the co-heating experiment with constant indoor air temperature in the
heavyweight building during two weeks in January and April
3. System identification
The following iterative methodology has been applied for selecting, identification and validation of a
suitable model of the heat dynamics of the building:
1. The first seven days of each data set, comprising respectively data from January 5 to January
12, and data from April 1 to April 8, are selected for model identification.
2. A sampling period containing sufficient information regarding the heat dynamics, is chosen.
Average values of the signals are calculated for the sampling period.
3. The model’s dependent variable(s) are selected, being the outdoor air temperature, total solar
radiation on the horizontal plane and the heating power as input variables, while the indoor air
temperature is selected as the output variable.
4. An ARMAX model is identified and validated with time series plots of the residuals and inputs,
based on the auto-correlation function (ACF) of the residuals and cross-correlation functions
(CCFs) from the residuals to the inputs. The data periods used for validation comprise January
12 to January 20 and April 8 to April 14.
5. If necessary, the mathematical order of the model is increased until the model’s output does not
contain any significant autocorrelation and cross-correlation.
6. The heat loss coefficient, time constants and effective thermal capacity of the building are
determined.
where τ is a time constant and the term αheat,1/ αheat,0 could be interpreted as the sum of time
constants for Q. With respect to the identified ARMAX models and the corresponding coefficients A,
B1, B2, and B3 obtained. The coefficients αi are replaced by the corresponding ARMAX parameters.
It should be noticed that the thermal capacity Ceff is only as exact as the coefficients of the model,
which in practice are determined by means of the data set obtained. The normal diurnal cycling does
not give enough time for a heat pulse to penetrate all parts of a massive building, which means that Ceff
is an effective thermal capacity. It is a measure of the heat storable under normal cycling. For slower
cycling Ceff would increase, and so would the time constants of the data fit. This illustrates to what
extent the parameters of the model are dependent on the frequencies represented in the data.
4. Results
This section presents the simulation results and the physical parameters estimated by the different
ARMAX models. Figure 3 presents a comparison of the measured indoor air temperature and the
predicted indoor air temperature by a second order ARMAX model, while using a multi-sinusoidal
heating power signal (S4). In addition, the figure presents plots of the ACF of the output’s residuals
and the CCFs between the residuals and the inputs, respectively the outdoor air temperature (Te), the
heating power (Q) and the solar radiation (Gv). No significant auto-correlation and cross-correlation is
observed while the residuals lie near the confidence band of 95%. With respect to the indoor air
temperature and the solar radiation, it could be argued that significant cross-correlation is observed.
Additional investigations, which are not presented in this paper, showed that using a separate input for
the direct solar radiation on each vertical facade instead of the total horizontal radiation on the
FIG 3. Comparison of the measured indoor air temperature and the indoor air temperature predicted
by a second order ARMAX model using a multi-sinusoidal (S4) heating power signal and the
corresponding autocorrelation and cross-correlation functions for the lightweight building
Figure 4 presents a comparison of the estimated HLC by the models based on the heating power
signals (Table 2) both for the lightweight and the heavyweight building and for January and April. The
ability of the models to predict the steady-state thermal properties of the building, i.e. HLC and the
solar aperture are analysed. The models are compared using a relatively small sample time of 5
minutes (including both fast and slow system dynamics) and larger sample times of 15, 60, and 120
minutes focusing on the slow dynamics and the determination of the steady-state thermal parameters.
First of all, figure 4 shows that the models identified based on data obtained from the standard co-
heating tests showed to be incapable of always giving a reliable prediction of the building’s HLC.
Since the indoor air temperature is constant during the measurement period, a reliable model
describing the building’s thermal dynamics is difficult to identify. Information regarding heat
dynamics is only obtained from over-heating periods. If overheating does not occur or seldom occurs,
which is the case in the heavyweight building, the capturing the building’s thermal dynamics is
impossible.
Second, the figure shows that the predicted HLCs for the lightweight building lie close to the
theoretical value of 105.6 W/K, while the predicted standard deviation is relatively small, both for
January and April based estimates. Moreover, model identification and prediction of the HLC was
possible for every signal except for signal S3 in April. This signal contains too little information to
excite the system sufficiently, comprising only 5 sines in a relatively large range of frequencies.
References
EN ISO 6946. 2007. Building components and building elements. Thermal resistance and thermal
transmittance. Calculation method. (EN ISO 6946:2007). International Standards Organisation.
EN ISO 13790. 2008. Energy performance of building – Calculation of energy use for space heating
and cooling. (ISO 13790:2008). International Standards Organisation
Hens H, Janssens A, Depraetere W, Carmeliet J, Lecompte J. 2007. Brick Cavity walls: performance
analysis based on measurements and simulations. Journal of building physics 31(2), pp. 95-124.
Jiminez MJ, Madsen H, Andersen KK. 2008. Identification of the main thermal characteristics of
building components using MATLAB. Building and Environment, vol. 43, pp. 170-180.
Madsen H. 2008. Time Series Analysis. Chapman and Hall.
Rabl. 1988. Parameter Estimation in Buildings: Methods for Dynamic Analysis of Measured Energy
Use. Journal of Solar Energy Engineering, vol. 110, pp. 52-66.
Wingfield, J. Bell M, Miles-Shenton D, Seavers J. 2011. Elm Tree Mews Field Trial – Evaluation and
Monitoring of Dwellings Performance Final Technical, Leeds Metropolitan University
SUMMARY:
The EU Directive on the Energy Performance of Buildings requires that energy certification of
buildings should be implemented in Denmark so that houses that are sold or let should have an energy
performance certificate. The result is that only a small part of existing houses has an energy
performance certificate. The Danish Building Research Institute has described a method that can be
applied for estimating the energy demand of dwellings. This is based on the information in the Danish
Building and Dwelling Register and requirements in the Danish Building Regulations from the year of
construction of the house. The result is an estimate of the energy demand of each building with a
variation. This makes it possible to make an automatic classification of all buildings. The paper
discusses the uncertainties and makes a sensitivity analysis to find the important parameters. The
variations are compared with measured energy demand. The method can be applied in other countries
with modifications for local building requirements, climate and building registers.
1. Introduction
In Denmark, energy performance certification has been in use since 1997 with energy performance
rating and for new houses (since 2006), including sold or let buildings. The energy rating classifies the
building on an efficiency scale ranging from A (high energy efficiency) to G (poor efficiency) (EPBD
2011). The rating is important when buying and selling houses. The problem is that many houses still
do not have an energy rating, and that it will take many years to cover the entire building stock. Until
ultimo 2013, only 22% of single-family houses had an energy rating. Energy ratings are given in
connection with a visit by an energy expert and the result is a rating and a description of possible
energy-saving measures. For houses that have not been energy rated, we do not know the energy
demand and potential energy savings. An automatic calculation of the energy rating would thus be
very useful. It is possible that some owners would make energy saving measures, if they were
informed about the energy rating of their property. It would also make it possible for authorities to
give an economic subsidy or tax deduction for the buildings with the poorest energy efficiency. The
result should be more energy saved for the same money.
It is possible to apply a normal certification method based on monthly values and details on heating
system and constructions to make energy demand calculations of existing houses. The problem is that
we have to visit the houses and we will not always be able to obtain information about the
constructions at the site or in drawings. Another important point is that this would be a very time-
consuming and expensive method. The Danish Energy Performance Certification Scheme for existing
buildings includes energy demand calculations based on a monthly calculation method (Aggerholm
and Grau 2005). This can be used at the visit as you can find some information and have additional
information from a guideline with typical constructions. The cost of this is paid by the owner before
selling the building. As it will take many years before all buildings have been officially rated, the
Danish Building Research Institute (SBi) has described a method (Bertelsen et al. 2011) that can be
used without physically visiting the house but only based on available information.
3. Geometric model
BBR contains no information on building height, length or width. So we have to make an estimate of
the geometry. The simplest form is a box-shaped building. We know the roof and floor area as it
equals the building footprint. For buildings without attic and basement, we use a box-shaped building.
We need the building width. Based on typical Danish buildings, we estimate the building width to be 9
m for single-family houses. The building length is the building footprint divided by the width.
Calculations with other length/width proportions show that it does not influence the surface area very
much. The storey height is assumed to be 2.8 m as a typical value. The volume and the external wall
area can then be calculated.
If the building has a useful attic area, we get an extra wall area based on the attic which is also a box.
The box is placed on top of the building box and does not change the floor or roof area. The length of
the attic is estimated to be the length of the building. The width and volume of the attic and wall area
can now be calculated.
For houses with a basement, this is calculated as a box placed under the building. As it is placed
underneath the building, there is no change in floor or roof area. The walls of the box give the
basement wall area. Note: We calculate as if the whole basement area is heated, as most people use it
like that. The length of the basement is estimated to be equal to the length of the building. The width
and volume of the basement and the basement wall area can now be calculated.
4. U-value model
U-values for the different building parts are taken from the Danish Building Regulations in force when
the building was constructed. Energy requirements were changed in 1961, 1967, 1972, 1977, 1982,
1985, 1995, 1998 and 2008. The five U-values used in the calculation are outer wall, basement wall,
floor, roof and window. We do not take into account that we can have different outer wall
constructions in the same house.
TABLE 1. Typical U-values and requirements in different periods
Heavy outer wall Roofs Windows
(> 100 kg/m2)
Before 1930 1.2 W/m2K Before 1950 0.97 W/m2K Before 1930 3.7 W/m2K
1931-1950 1.5 W/m2K 1951-1960 0.6 W/m2K 1931-1994 2.9 W/m2K
1951-1960 1.5 W/m2K 1961-1976 0.4 W/m2K 1995-2007 1.8 W/m2K
1961-1971 1.1 W/m2K 1977-1994 0.2 W/m2K From 2008 1.5 W/m2K
1972-1976 1.0 W/m2K From 1995 0.15 W/m2K Note: We assume that most old
1977-1994 0.4 W/m2K windows are upgraded.
1995-2007 0.3 W/m2K
From 2008 0.2 W/m2K
BBR information on the exterior wall material is used to decide whether it is a light or heavy wall as
the Danish Building Regulations define different U-values for these wall types. Data about older
houses are supplemented by the experience of SBi. Table 1 is an example of the U-value variation
over the years for heavy outer walls, roofs and windows. Data for other construction types and periods
are given in a Danish report (Bertelsen et al. 2011).
If we compare the official energy rating for the five houses that have an official energy certificate with
the result of the automatic rating method based on BBR it is seen that the officially calculated energy
demand by the energy expert lies between the upper and lower limits of the automatic method. As
expected, the value is closer to the lower limit, as some energy-saving measures were made. These are
typically additional insulation, tightening and retrofitting of windows and new heating system.
FIG 2. Measured energy consumption (dots) in Danish houses for different construction year. The red
line is the requirement in the Danish Building Regulations.
Figure 3 shows that old houses have a high energy demand and large variation and new houses has
lower energy demand and less absolute variation (Jensen 2004). The arrows indicate the average
values. This shows that our assumptions for variations are realistic.
References
Aggerholm, S. and Grau, K. (2005). Bygningers energibehov - PC-program og beregningsvejledning.
(Building energy demand – PC program and user guide) SBi-Anvisning 213. Statens
Byggeforskningsinstitut (SBi), Hørsholm, Denmark
Bertelsen, N. H., Nielsen, A., Sørensen, N.L. and Wittchen, K.B. (2011) Automatisk
energirammeberegning for den eksisterende bygningsmasse (Automatic energy performace
calculation for the existing buildings), SBi-rapport 2011-20, Statens Byggeforskningsinstitut (SBi),
Hørsholm, Denmark
Danish Building Regulations (in Danish: BygningsReglementet. Erhvervs- og Byggestyrelsen,
København. Use www.ebst.dk – to find older building regulations
EPBD (2011). Implementation of the Directive on the Energy Performance of Buildings - Country
reports 2010. EPBD Buildings Platform, Brussels, European Commission ISBN: 978-92-9202-090-
3.
Jensen, O. M. (2004). Barrierer for realisering af energibesparelser i bygninger. (Barriers for
realisation of energy saving in buildings), Statens Byggeforskningsinstitut (SBi), Hørsholm.
Petersen, T.D. (1997). Uncertainty analysis of energy consumption in dwellings, NTNU, Dr.avh.
1997-122, Trondheim, Norway
Wittchen, K. B. (2009). Potentielle energibesparelser i det eksisterende byggeri. (Potential energy
savings in existing buildings) Forskningsprojekt SBi 2009:05. Statens Byggeforskningsinstitut,
Hørsholm, Danmark
Wittchen, K.B, and Kragh, J. (2012) Danish building typologies, Danish Building Research Institute,
SBi 2012:01, Hørsholm, Denmark
Keywords: Windows, glazing, solar shading, thermal transmittance, U-value, Hot box, laboratory,
simulation
SUMMARY:
Shading systems are widely used in conjunction with glazed facades, especially in office buildings.
Solar gains in combination with high internal gains will often lead to a cooling need during
operational hours. Thus, shading devices are traditionally used in order to reduce cooling demands
during the time of day when the sun is up.
The possible effects of a solar shading device on the thermal transmittance value (U-value) of the
glazing unit has, however, not been given much attention.
In order to assess energy saving potentials by applying shading devices, detailed information about the
component level characteristics is necessary. Thermal and optical properties of these components are
important input data for simulation tools and the calculation of cooling and heating demands as well as
daylight accessibility in buildings.
Shading systems can be categorized into three main types; external, internal and in-between panes
shading devices. When an in-between pane shading device is deployed (i.e. closed shut), it can act as
an additional layer in the insulating glazing unit (IGU) and thus reduce the U-value of the IGU. In a
cold climate, such a system could be utilized in order to reduce heat losses during periods of time
when a room is not occupied and when no view to the outdoors is necessary.
In this article, the thermal properties of such a system have been explored. Hot-box measurements
have been carried out for a window with 4-panes. The glazing unit is built up as a 3-pane IGU with an
external coupled fourth pane. A venetian blind shading device is placed between the outermost,
coupled pane and the IGU. U-values have been measured for various slat-angles of the blinds. These
results have been compared with calculated values using the THERM and WINDOW simulation
software. This is software developed for calculations of thermal and optical properties of entire
windows.
1 Introduction
Some of the previous work pertaining to louvered blinds has been done for a horizontal blind sealed
between the panes of an Insulated Glazing Unit (IGU). Garnet et al. (1995) carried out an
In spite of a large amount of theoretical work carried out in order to assess thermal and optical
performance of solar shading systems, few measurement campaigns have been carried out. In 2001,
Breitenbach et al. (2001) presented thermal resistance measurements for a two-pane IGU with
integrated venetian blinds. The authors found that the thermal transmittance value (U-value) varied
almost linearly with the slat angle from 2.92 ± 0.14 W/(m²K) with blinds retracted to 1.44 ±
0.21 W/(m²K) with the blinds shut and slats in vertical, closed position. This contradicts the values
found by Tzempelikos (2005) during a measurement campaign performed on a double-pane window
with in-between venetian blinds. Here it was found that a blind-tilt angle of 120° gave the lowest
thermal resistance and that blinds closed with the slats in vertical (i.e. 180°) position gave the highest
thermal resistance.
None of the two aforementioned publications present values for a situation with the blinds in a
retracted position and only the centre-of-glass U-values are presented. Presenting only centre-of-
glazing U-values does not take into account any effects caused by the mounting system and any motor
of the venetians on the total U-value of the glazing unit. Such effects could be substantial as mounting
systems and motors may contain large amounts of continuous metal spanning from the warm to the
cold side.
In 2006 Huang et.al (2006) presented measurements performed using a Guarded Heater Plate
apparatus (GHP). Here, centre-of-glass U-values (Ucog) for double-glazing units with and without
shades, were measured using Heat Flow Meters (HFM). This was carried out by controlling the
temperature on each side of the sample using the GHP. 12 sample configurations were measured. The
measured values were used for benchmarking/validation of CFD simulations. Values are extrapolated
to also include solar heat gain coefficients (SHGC) and visual transmittance (Tvis) values.
A study performed by Shahid H. and Naylor D. (2005) presents a numerical method for a window
with internal shading. The method is validated by use of previous measurement results for tall vertical
cavities. No validation using measurements on real windows was performed.
In 2009, Laouadi (2009) presented an overview of existing studies related to modelling of glazing
units coupled with solar shading systems. A model for calculating properties of the centre-of-glazing
area were proposed and found to correspond to measurements performed by Huang within a 7 %
margin of error. The model is only verified for double pane glazing units with a limited (lower than
25 mm) cavity thickness. The U-value of IGU's with low-e coated glazings and higher cavity
thicknesses (40 mm was modelled) were underestimated with the model. Furthermore, the authors
point out that at that time, the models regarding the thermal performance of shading systems are based
on simple algorithms and that further validation work needs to be carried out.
Figure 1. Graphical illustration of the 4-pane window. Shading device, in dark blue, shown in
retracted position. The shading device is integrated in the cavity behind an exterior coupled glass
pane. The frame is insulated with polyurethane foam (coloured yellow).
Table 1 shows an overview and description of the window configurations that have been measured.
Generic data as described in ISO 15099 (ISO 2003) has been used as material data for the shading
slats.
2.2 The hot box test facility – test procedure and instrumentation
Measurements have been carried out in a guarded hot box apparatus, as pictured in Figure 2, according
to procedures described in ISO 8990:1994 Determination of steady-state thermal transmission
properties Calibrated and guarded hot box (ISO 1994) and ISO 12567-1:2010 Thermal performance
of windows and doors - Determination of thermal transmittance by the hot-box method - Part 1:
Complete windows and doors (ISO 2010).
Figure 2. The guarded hot box test facility used for the measurements
The tests were performed at steady state conditions of +20 °C interior and 0 °C exterior temperatures.
Window U-values are calculated based on the measured heat flows, surrounding temperatures and
window area.
The window was mounted in a template constructed as a sandwich element consisting of a 100 mm
thick Extruded Polystyrene (XPS) core, clad with 12 mm plywood on the faces exposed to the hot and
cold sides of the hot box. The joints between the window and the surround panel are sealed with tape
on both sides to ensure an air tight seal. Minor gaps between window frame and template were filled
with Expanded polystyrene (EPS).
The metering area of the hot box is 2.45 m x 2.45 m. The window is placed in the middle of the
metering area template wall at a distance of 1.0 m from the floor to the lower edge of the frame.
Surface thermal resistances were adjusted close to the standardised ones prior to the tests by adjusting
air flow velocities adjacent to the template surface on the cold side. Natural convection driven air-flow
were maintained on the warm side. However, during the measurements the surface thermal resistances
differ slightly from the standardized values. Corrections have thus been made for these deviations, so
that all U-values are stated with normalised surface thermal resistance values as specified in ISO
8994:1994 (ISO 1994). The standardized conditions are; interior surface thermal resistance,
Rsi = 0.13 W/m²K and external surface thermal resistance Rse = 0.04 W/m²K.
Centre-of-glazing U-values (Ucog) were measured using two PU 43 T heat flow meters, from
Hukseflux (Hukseflux), with a declared accuracy of ± 5 % at 20°C. These were mounted mid-height
approximately 10 cm from the vertical symmetry axis of the IGUs on the warm side of the samples.
2 2 2
Δ PU w ⎡ Δ P Φw ⎤ ⎡ Δ P Α w ⎤ ⎡ Δ P δθie ⎤
= ⎢ ⎥ + ⎢ ⎥ + ⎢ ⎥ (1)
Uw ⎣ Φw ⎦ ⎣ Α w ⎦ ⎣ θie ⎦
The uncertainty in the sample heat flow is based on the heat balance equation for the metering
chamber. The uncertainties in test sample specimen heat flow, ΔPΦw/Φw, is expressed using Eq. 2.
2 2 2 2
Δ P Φw ⎛ Δ P ΦIN ⎞ ⎛ Δ P Φsur ⎞ ⎛ Δ P ΦEXTR ⎞ ⎛ Δ P ΦFL , w ⎞
= ⎜ ⎟ + ⎜ ⎟ + ⎜ ⎟ + ⎜⎜ ⎟⎟ (2)
Φw ⎝ Φw ⎠ ⎝ Φw ⎠ ⎝ Φw ⎠ ⎝ Φw ⎠
A calibration experiment for the thermocouples was carried out prior to the hot box measurements
using a reference temperature bath. The relative scattering in measured temperatures between the
thermocouples were found to be lower than 0.02 °C. Thus it can be concluded that the influence from
the factor ΔPδθie as described in Eq.1 is negligible.
The uncertainties stated in this work, are given with a coverage factor of two standard deviations and
the corresponding 95 % confidence interval.
Simulations of the thermal properties of the window with the integrated shading device have been
carried out using the THERM 7.1 and WINDOW 7.1 software developed at Lawrence Berkeley
National Laboratories in Berkeley (LBNL 2013; LBNL 2013). THERM is a finite element method
program suitable for calculating U-values of window frames, whereas the WINDOW software
calculates thermal and optical properties of glazing units. Thus, THERM has been used for calculation
of frame U-values and linear edge transmission losses caused by the spacer used. Cavities in the
frames have been modelled according to the NFRC 100-2001 definitions in THERM. WINDOW has
been used for calculating the centre-of-glazing U-values including the shading device. The results are
compared to the values measured in the hot box.
3 Results
3.1 Measurement results
From Figure 3 and Table 2, one can see that the effect of operating the blinds has minor or no effect
on the measured mean U-value. Deploying the blinds with horizontal slats, increases the U-value
slightly, whereas closing the slats shut (vertical slats) will give a slight decrease in U-value. The
alternations of the U-value as function of blind positions were found to be negligible in terms of
statistical significance, which means that the individual measured mean U-values lies inside the
measurement error boundaries.
Figure 3. Measured U-values for the window. Upper and lower uncertainty values shown with lines.
The measured U-values for the whole window, Uwindow, and the center-of-glazing U-values, Ucog, are
shown in Table 2.
Table 2. Measured mean U-values, Uwindow, and center of glazing U-values, Ucog, for the window.
retracted deployed blinds, deployed blinds, deployed blinds,
Measured value
blinds horizontal slats 45°deg slats vertical slats
Uwindow (W/m²K) 0.84 ± 0.04 0.85 ± 0.04 0.83 ± 0.04 0.81 ± 0.04
Ucog (W/m²K) 0.61 ± 0.03 0.60 ± 0.03 0.58 ± 0.03 0.57 ± 0.03
The outermost cavity in the glazing unit is an air filled cavity. Tests carried out using THERM to
model this cavity treating the cavity as both a slightly ventilated and un-ventilated cavity according to
the NFRC definition used in THERM (LBNL 2013), gave U-value-results comparable to the values
calculated assuming an air-filled cavity using WINDOW. Assuming a slightly ventilated cavity gave
center-of glazing U-values approximately 4 % higher than an un-ventilated cavity. However, based on
the minor openings constricting air cross-flow in the cavity, it has been treated as an un-ventilated
If one looks at the measured center-of-glazing U-values (Ucog), the same behavior as for the Uwindow-
values can be seen. Thus the two independent measurments support each other in that the effects on
the window U-value as function of shading device slat angle and position is negligible.
The U-values found from simulations using THERM and WINDOW were consistently lower than the
measured values but also here, the effects of blind posistion and angle ion the U-value were found to
be negligible. The deviation between the measured and calculated values can likely be explained by
differences between declared (from the producer) and actual properties of the low-e coatings as well as
the amount of argon in the gas-filled cavities of the IGU.
As opposed to finding any potential gains in U-value by the introduction of such a shading device, it
was found to have a negative effect on the U-value of the window compared to one without the
shading device. The protruding aluminum components of the shading device and venetian blinds, act
as thermal bridges. Numerical simulations showed that the U-value of the window with the shading
device will increase with approximately 5 % compared to a window without an integrated shading
device. Thus, any beneficial effects expected to be achieved by using the venetians as an additional
layer in the IGU, was found to be more than counteracted by the thermal bridging.
5 Conclusions
Measurements have been carried out in order to investigate the U-value of a window with a 4-pane
glazing unit incorporating an in-between panes venetian blind shading device. The aim of the study
has been to assess the effect of operating the blinds on the window U-value. Various slat angles and
blind positions have been studied. It was found that no statistically significant effect could be found on
the U-value as function of various slat angles.
Numerical simulations showed that the U-values were underestimated compared to the measured
values. This can have several reasons, including deviations in gas-filling levels and quality of the low-
emissivity coatings in the IGU. Identifying the reasons for these deviations is the aim of further studies
that are currently being carried out within the Research Centre on Zero Emission Buildings (ZEB).
6 Future work
There are several possibilities which should be explored in order to achieve more effective shading
devices, and the following should be explored
• Reduce slat thermal conductivity in order to reduce thermal bridging effects
7 Acknowledgements
The authors gratefully acknowledge the support from the Research Council of Norway and several
partners through the Research Centre on Zero Emission Buildings (ZEB).
8 References
Breitenbach, J., S. Lart, et al. (2001). "Optical and thermal performance of glazing with integral
venetian blinds." Energy and Buildings 33(5): 433-442.
Garnet, J. M., R. A. Fraser, et al. (1995). Effect of internal Venetian blinds on center-glass U-values.
Window Innovations Toronto, Canada.
Huang, N. Y. T., Wright, J.L., Collins, M.R. (2006). "Thermal resistance of a window with an
enclosed Venetian blind: Guarded Heated Plate measurements." ASHRAE Transactions
112(2): 13-21.
Hukseflux. "PU Series." Retrieved 06.08.2013, from
http://www.hukseflux.com/sites/default/files/product_brochure/PU%20series%20v0717.pdf.
ISO (1994). ISO 8990:1994 - Determination of steady-state thermal transmission properties Calibrated
and guarded hot box.
ISO (2003). NS-ISO 15099 – Thermal performance of windows, doors and shading devices - Detailed
calculations, International Organization for Standardization.
ISO, N.-E. (2010). NS-EN 12567:2010 Thermal performance of windows and doors - Determination
of thermal transmittance by the hot-box method - Part 1: Complete windows and doors.
Laouadi, A. (2009). "Thermal modeling of shading devices of windows NRCC – 51121." NRCC-
51121, ASHRAE Transactions 115(2): 1-20.
LBNL. (2013). "THERM 7.1 beta." Retrieved 06.08, 2013, from
http://windows.lbl.gov/software/therm/7/index.html.
LBNL. (2013). "WINDOW 7.1 beta." Retrieved 06.08, 2013, from
http://windows.lbl.gov/software/window/window.html.
Rheault, S. and E. Bilgen (1989). "Heat transfer analysis in an automated Venetian blind window
system." Solar Energy Engineering 111(1): 89-95.
Shahid, H. and D. Naylor (2005). "Energy performance assessment of a window with a horizontal
Venetian blind." Energy and Buildings 37(8): 836-843.
Tzempelikos, A. (2005). A methodology for integrated daylighting and thermal analysis of buildings.
Faculty of Engineering and Computer Science - Building, Civil and Environmental
Engineering, Concordia University.
SUMMARY:
The development of complex physical models in building simulation creates a great demand for
powerful and flexible transient simulation programs, and thus numerical solvers. Especially when
turning to three dimensional models computational speed and memory use strongly limits the problem
size of an application. The calculation complexity increases further more when taking anisotropic
transport processes into account. This article describes a solution for a transient anisotropic three
dimensional solver for combined heat, air, and moisture transport in porous materials of building
envelop systems, utilizing the Delphin solver framework. After describing the physical model,
boundary conditions, and certain aspects of internal solver optimisation regarding memory use and
calculation speed, three initial academic verification examples will be presented. The influence of
anisotropic material properties onto three dimensional simulation results is researched, targeting an
evaluation of simulation performance and correctness while comparing to isotropic simulation results
of the same experiments. The paper closes with an outlook onto parallelisation of the introduced
solver solution.
1. Introduction
Prediction of damage in constructions (e.g. timber beam heads, interior insulations, etc.) is recently
done by performing transient transport simulations. It is state of the art to utilise continuum physics
for describing conserved quantities coupled through systems of partial differential equations. Our
Delphin solver framework (Nicolai 2011) implements continuum physics models for porous media
(construction materials) by defining a representative elementary volume (REV) following Whitaker
(1969), Bear and Bachmat (1992) and Whitaker (1986). Grunewald (1997) and Nicolai (2007) defined
essential conservation equations for hygrothermal, pollutant as well as salt transport. Moreover a
lower bound (infimum) as a function of the pore radius spectrum was defined as a constraint for an
accurately selected REV size.
The numerical solution of resulting partial differential equations is done using a spatial discretisation
that introduces location-fixed control volumes and leads to a large set of ordinary differential
equations to be integrated. Numerical integrators developed to solve such problems implement some
form of time discretisation and solve the problem in a step-wise manner. Each integration step (time
step) involved hereby a solution of a non-linear equation system, which in turn requires solving large
sparse linear systems of equations. The actual implementation is done in several cascaded
computational loops (see Listing 1).
Increasing dimension of simulation space (e.g. 1D, 2D to 3D) introduces new degrees of freedom, and
thus complexity to all calculation procedures executed during a simulation. It basically increases
To emphasis additional effort of multidimensionality a vector arrow marks all divergence symbols,
usually this vector arrow is implicitly defined in a divergence and therefore omitted. In isotropic
models the transport coefficient is a scalar. In anisotropic models the transport coefficient tensor L is
FIG 1. Example for a material and simulation reference system that are not aligned
The relation between material reference system and simulation reference system may change relative
to the spatial location. In FIG. 1 the detail (right) shows the MRS for a specific location timber beam
end. The direction W marks the tangential direction, V marks radial direction with respect to the
material (MRS). Moving to a different location will change the orientation of the MRS within the
global simulation coordinate system.
The mathematical transformation process responsible for mapping these two reference systems is thus
location dependent. For completeness we briefly review the derivation of rotation matrices involved.
Detailed information can be found in standard literature on the matter.
A starting position ( uP , vP ) may now be substituted, following equations (1 - 2) into the found
equations (3 - 4) , which drops out all unknowns. The rotation equation, written in matrix notation, for
a rotation by an angle around the w-axis reads:
0 0 1 0 0
The rotation axis u and v can be derived equivalently:
1 0 0 cos 0 sin
Ru 0 cos sin (5) Rv 0 1 0 (6)
0 sin cos sin 0 cos
A three dimensional point in a u,v,w coordinate system is transformed by stepwise execution of single
direction rotations until it matches with the target x,y,z coordinate system.
pxyz RT puvw
RT Ru Rv Rw
cos v cos w cos v sin w sin v
RT sin u sin v cos w cos u sin w ( sin u )sin v ( sin w ) cos u cos w sin u cos v
cos u sin v cos w sin u sin w cos u sin v sin w sin u cos w cos u cos v
In analogy a material parameter/transport coefficient tensor should be used in SRS, this coefficient
tensor is first represented in MRS as Lk , and afterwards multiplied with the transformation matrix
k
RT . L RT Lk
With respect to the numerical solution we have to keep in mind that the rotation matrix is
location dependent.
For large 2D and generally all 3D simulation setups efficient memory handling alongside a fast flux
calculation is essential.
FIG 3. Simulation Time Comparison of the Direct (left) and the Iterative GMRES Solver (right) both
utilizing the CVODE integrator.
6. Conclusions
An anisotropic three dimensional simulation model was developed and implemented in C++ within
the Delphin 6 modelling framework. The framework was tested with scientific test cases to illustrate
the functionality of the anisotropic property transformations with different mapping assignments (FIG
4).
FIG 4. Anisotropic Transport for Spruce, Oak and Beech (24h) sorption in a corner experiment
The computational speed for the 3D tests was enhanced by selecting an iterative GMRES solver with
ILU preconditioning. But large hygrothermal 3D cases take still too much time while running as serial
iterative code, future work will be devoted to a parallel implementation of the Delphin 6.1 model. The
preconditioner ILU must thus be replaced since its code is hard to parallelise. Another approach will
be applying domain decomposition techniques to split large problems onto multiple processors.
For an accurate comparison of different solution algorithms a standardised set of three dimensional
test benchmarks is required. Those test shall not only specify results, constructions, materials, and
conditions but also truncation error control parameters to ensure a valid comparison of a simulations
wall clock times or other convergence properties after valid simulation results are meet.
7. Acknowledgements
The authors gratefully acknowledge funding support by the German Federal Ministry of Economics
and Technology under the Contract No. 0329663N.
SUMMARY:
This text studies the required temperature accuracy criterion for mould growth index calculations in
crawlspace. First outdoor air conditions data and the updated VTT-TUT mould growth model are
used to evaluate the accuracy criterion for temperature field. Then heat conduction and air change
simulations on crawlspace are made and the temperature accuracy criterion is used to give guidelines
on modelling a crawlspace structure.
When the ground under the crawlspace is initially at outdoor air average temperature, it is
recommended to have a minimum of three years pre-conditioning time. Based on the literature and
simulations, it is recommended to have the cut-off planes in the ground at a distance of at least three
times the periodic penetration depth or the width of the building. The temperature accuracy limit
±0.13 °C doesn’t include many essential parameters, but gives direction on how sensitive the updated
VTT-TUT mould growth model is for the change in temperature and humidity conditions.
1. Introduction
1.1 General
According to the European directive 2010/31/EU all new buildings must be near zero energy buildings
by the end of year 2020. To design and build truly energy efficient and durable buildings it is essential
to understand the hygrothermal behaviour of building structures. Numerical simulations provide a
powerful way to contribute to this knowledge.
Calculation of temperature field is an essential part of evaluating the hygrothermal behaviour of
building components. For example the relative humidity is a key factor in material moisture content
and many degradation phenomena. Relative humidity depends essentially on temperature and because
of that the temperature conditions should be evaluated as accurately as possible. Different authors
however use different criterion for evaluating the applicability of a simulation model. This text aims to
contribute to that by suggesting a method to determine limits for temperature accuracy.
Crawlspace is a building foundation type where there is an air space between the soil and living areas
of the building. The structure allows making the building on several types of soils and for example on
top of piles. Crawlspace structure is also suitable when making the building from volume elements
and the crawlspace is also recommended in radon areas. However the major drawback of the structure
is its inadequate hygrothermal behaviour in the summertime. The temperature accuracy criterion is
used to give reference values when modelling crawlspace structure.
2. Methods
2.1 Temperature accuracy criterion
The temperature accuracy criterion was determined with the following steps:
1. Hourly T/RH data for 30 years was acquired from the Finnish Meteorological Station from
four different localities (Sodankylä, Jyväskylä, Jokioinen and Vantaa).
2. For every hour of the year, new relative humidity was calculated by assuming a constant
temperature change and a constant amount of water vapour in that time step. When
temperature was decreased the relative humidity was limited to a maximum of 100 % RH.
3. The updated VTT-TUT mould growth model (Ojanen et al. 2010) was used to calculate the
maximum value of mould index for every year and as a function of temperature change.
Mould index calculation was done for the material class: “Very sensitive” (although no actual
material surface was present).
4. Cumulative distribution function of the impact of temperature change was plotted for the set
of 30 values (ΔM(ΔT)) at a time. Two maximum values (ΔM = 0.6 and ΔM = 1) were used
and the temperature change where less than 90 % of the years stay under the limit was chosen.
2.2 Simulations
The crawlspace geometry and other simulation model information are shown in the next drawing.
Figure 2.1 The geometry of studied crawlspace. Floor insulation is 200 mm and the thickness of
foundation wall layers is 100 mm. The outside ground is 500 mm lower than the floor inside surface.
For the 300 mm length of wall, adiabatic boundary condition is used. The amount of ground in the
simulations (parameter z) is given values of 1, 2, 4, 8, 16 and 32 m. Monitoring point Tairmiddle is
also shown.
As initial condition the yearly average of outdoor air temperature was used (5.63 °C). Long-wave
radiation in the y-direction in crawlspace was taken into account where it was possible, but not
between the wall and horizontal surfaces.
In Comsol Multiphysics the “Finer” Physics-controlled mesh was used for all calculations. For Wufi
calculations the default discretization was used without modifications. In Delphin calculations first the
For yearly variations the periodic penetration depth in soil is δsoil = sqrt(1·10-6 m2/s ·8760·3600 s/π) =
3.2 m
Figure 3.1. The impact of temperature change on the calculated mould index (sensitivity class: very
sensitive). One line is associated with a one year of outdoor air T/RH conditions. The outdoor air
moisture content is kept constant while the air temperature is decreased or increased. The maximum
value of mould index is plotted against the temperature change.
Big differences exist between years. For example during 1980-2009 in Jokioinen the maximum value
of the yearly outdoor air mould index has varied from 0.6 to 5.8 (ΔT = 0 °C). In Jyväskylä, Jokioinen
and Vantaa the variability between years becomes smaller the better (warmer) or worse (colder) the
conditions become. Sodankylä is located in Northern Finland and the average temperature there is
lower than in the other three cities. It seems that the outdoor air conditions in Sodankylä are not as
critical for mould growth as they are in other cities and even with some average temperature drop, the
maximum mould index is not very high for all years.
The temperature has a big impact on the resulting mould index. A temperature change of 0.5 °C can
change the mould index value over two units (scale 0...6). This leads to the conclusion that the
determination of temperature field accurately is essential for accurate evaluation of mould growth and
the hygrothermal behaviour of building envelope structures. A second conclusion is that it is possible
Figure 3.2 Results from crawlspace simulations. Left: Impact of simulation time and the choice of
simulation software on the yearly average crawlspace temperature (z = 8 m, adiabatic). Right:
Maximum difference in the hourly values compared to the last simulation year (Crawlspace air,
Delphin, q = 0.03 W/m2).
Temperature changes faster in the beginning of the simulation than after a few years. The conditions
approach asymptotically their periodical equilibrium. In the beginning of each year the conditions
differ more from the last simulation year than the end of the year. In many cases two years pre-
simulation time would have been enough but three years pre-simulation was enough in all calculated
cases.
The next figures show the impact of increasing the amount of soil in the simulation model.
References
Adjali, M.H., Davies, M., Ni Riain, C. & Littler, J.G. (2000) In situ measurements and numerical
simulation of heat transfer beneath a heated ground floor slab, Energy and Buildings 33 (2000) 75-
83.
Cheng, Dong (2013) Dynamic three-dimensional heat transfer calculation for uninsulated slab-on-
ground constructions. Energy and Buildings 60 (2013) 420-428.
Hagentoft, Carl-Eric (2001) Introduction to Building Physics. Studentlitteratur. Lund, Sweden. ISBN
91-44-01896-7
Janssen, Hans; Carmeliet, Jan & Hens, Hugo (2004) The influence of soil moisture transfer on
building heat loss via the ground. Building and Environment 39 (2004) 825-836.
Ojanen, Tuomo; Lähdesmäki, Kimmo; Viitanen, Hannu; Vinha, Juha; Peuhkuri, Ruut & Salminen,
Kati (2010) Mold growth modeling of building structures using sensitivity classes of materials.
Presented at the Thermal Performance of the Exterior Envelopes of Whole Buildings XI
International Conference.
Pallin, Simon & Kehrer, Manfred (2012) Hygrothermal simulations of foundations: Part I: Soil
material properties. Journal of Building Physics 37(2) 130-152.
Rees, S. W., Zhou, Z. & Thomas, H.R. (2007) Ground heat transfer: A numerical simulation of a full-
scale experiment. Building and Environment 42 (2007) 1478-1488.
SFS-EN ISO 10211. 2008. Thermal bridges in building construction – Heat flows and surface
temperatures – Detailed calculations (ISO 10211:2007). Finnish Standards Association SFS.
SFS-EN ISO 13370. 2008. Thermal performance of buildings. Heat transfer via the ground.
Calculation methods (ISO 13370:2007). Finnish Standards Association SFS.
Zoras, Stamatis; Georgakis, Chrissa; Kosmopoulos, Panagiotis & Dimoudi, Argiro (2011) Multi-year
application of the three-dimensional numerical generation of response factors (NGRF) method in
the prediction of conductive temperatures in soil and passive cooling earth-contact components.
Solar Energy 85 (2011) 2275-2282.
SUMMARY:
Overheating, glare and high energy consumption for air conditioning and ventilation are recurrent
problems in office buildings in Santiago of Chile over the whole year. Solar heat gains are the major
cooling load. Large variety of complex fenestration systems (CFS) such as perforated screens,
louvers and motorized rollers are being implemented on modern office buildings. However, it is
unknown their effectiveness to control short-wave solar radiation before it is transmitted through the
glazed building façades. This paper aims to study the performance of CFS currently used in existing
office buildings in Santiago of Chile by measuring the transmission of short-wave solar radiation in
four office buildings at different façade´s orientations. The main results show that perforated and
undulated exterior screens and louvers and motorized rollers reduce peaks of solar heat transmission
up to 2% of incident solar radiation, while solar transmissions can be as low as 8%. Also, this study
reveals that CFS are more effective to reduce solar heat transmission than spectrally selective
windows with low/medium SHGC. However, non-proper architectural arrangements of CFS can
significantly diminish their performance.
1. Introduction
Overheating, glare and high energy consumption are recurrent problems in office buildings in
Santiago of Chile over the whole year. This situation is consequence of the combination of a warm
climate and modern architecture of office buildings. Figures 1a and 1c show the hourly temperature
and beam solar radiation of an average day in January, March, July and October, while Figure 1b
shows the global horizontal solar radiation during the whole-year (ASHRAE 2005). It is observed that
Santiago’s climate (S 33° 22'; W 70° 46') is characterized by high temperatures and solar irradiance
for around 8 months of the year. This climatic condition exposes office buildings to large heat gains
due to solar radiation and heat conduction.
This situation would not be a problem if office buildings would be properly designed for this weather
condition. Otherwise, modern architecture imposes buildings with completely or large glazing façades
as shown in Figure 2. Modern architecture, warm weather and internal gains (people, equipment and
lighting) cause excess of solar heat gains and daylighting transmission to the indoor environment. For
instances, it is observed in Figure 2b that internal rollers in Titanium Tower are down in the north
FIG 1: Weather data of Santiago of Chile. a) Hourly temperature. b) Global horizontal solar
radiation. c) Hourly beam solar radiation.
Solar heat gains through windows are the major component of the building energy performance and
occupant’s comfort. Several authors have reported the large contribution to cooling loads of solar heat
gains through fenestration in warm climates (Reilly and Hawthorne 1998; Li and Lam 2000;
Winkelmann 2001; Kuhn 2006). Exterior shading devices are the most effective way to reduce solar
heat gains through glazed façades because they intercept direct solar radiation before it reaches the
glass. Fully shaded glazing façades can reduce solar heat gains up to 80% (ASHRAE, 2005).
Building C: San Agustín Building is a five-story university building with teaching labs, computer
labs, classrooms and offices. The east façade is composed of single-pane clear glass, while the
west façade consist of double glazed clear glass and exterior perforated and undulated screens on
the 3rd and 4rd floors of the building.
Building D: MideUC Building is a university building mainly composed of offices and meeting
rooms. The envelope consists of double glazed clear glass with different CFS on each façade. The
north façade includes perforated and undulated louvers, while undulated and perforated screens
are placed on the east and west façades.
The CFS implemented in these buildings are motorized exterior rollers in Building A, perforated and
and ondulated louvers in Building D, and ondulated and perforrated screens in Buildings C and D.
Building B does not count with CFS but it has double glazing windows with medium SHGC. Finally,
the east façade of Building C is composed of windows with a single clear glass.
Building A Building B
Building C Building D
FIG 3: Monitored buildings
(a) (b)
FIG 5: a) Mounting diagram of piranometers (dimensions are in millimeters). b) Comparison of
vertical solar radiation measurements of piranometers SP Lite 2 and CMP 11.
SP Lite2 sensor has a spectral range from 400 to 1100 nm. Since this range covers a fraction of the
whole solar spectral range, thus they measure a range much smaller than more precise piranometers, it
was necessary to evaluate if SPLite2 were suitable for this study. This was done by means of
comparing the measurements of solar radiation of three sensors SP Lite2 against the measurements of
sensor CMP11 of Kipp&Zonen, which has a spectral range from 285 to 2800 nm. Measurements were
carried out for four days during winter in a vertical plane oriented to north. Figure 5b shows that
FIG 6: Measured SW solar radiation transmission through different fenestration systems, buildings,
façade orientations, and season.
The following sections presents detailed results of incident solar radiation on building façades, solar radiation
between the CFS and glass, and short-wave solar radiation transmitted through the fenestration system towards
the indoor environment.
3.1 Building A
Figure 7 shows the measurements of solar radiation on the west façade in July 2012 (winter) and
January 2013 (summer). In both figures is clear the effectiveness of exterior rollers to control the
transmission of SW solar radiation. Despite of the high incident solar radiation on the northwest
façade, the peak of transmitted SW solar radiation is below 20 W/m2 in winter and summer. This
(a) (b)
FIG 7: Incident and transmitted solar radiation through fenestration systems on northwest façade of
Building A: a) July 2012. b) January 2013.
3.2 Building B
Figure 8 deploys the measurements of incident solar radiation on the northeast façade in August
(winter) and October (spring). Figure 8a shows peaks of incident solar radiation around 800 W/m2 and
600 W/m2 in sunny days of winter and summer, respectively. In summer, incident solar radiation is
lower than that in winter due to higher incident sun angles. Similarly to the measurements in Building
A, it is shown that the transmitted SW solar radiation is higher in winter than that in spring and
summer. This fact, in addition to that the people occupancy ratio of buildings is usually higher in
Chile than that in North America and Europe, causes high energy consumption for cooling in office
buildings even though in winter.
Figures 8a and 8b show that the double glazed window only transmits around 25% of the incident
solar radiation. This indicates that spectrally selective windows with medium to low SHGC allow
controlling solar heat gains effectively during the whole year. However, they still allow transmitting
20 to 25% of incident solar radiation, which is significant in climates with high solar radiation even in
winter.
3.3 Building C
Measurements were carried on the east and west façades in an unoccupied office and a small library,
respectively, that are located on the third floor. Main results are presented below.
3.4 Building D
Figure 10 shows that the CFS on the east and north façades significantly blocks the incident solar
radiation. On the east façade (Figure 10a), the peak of transmitted SW solar radiation is about 380
W/m2 but it occurs for very short time (< 30 min). The rest of time when sun rays strike the east
façade, the transmitted SW solar radiation is lower than 70 W/m2, which is only 10% of the incident
(a) (b)
FIG 9: Incident and transmitted SW solar radiation in San Agustín Building in November: a) East
façade (single glazed clear glass). b) West façade (double glazed clear glass with perforated and
undulated exterior screen).
4. Conclusions
This paper presented a study of the effectiveness of typical glazing façades and CFS to control solar
heat gains in existing buildings in Santiago of Chile. The main conclusions that could be drawn from
this study are:
CFS are more effective than double glazed windows with low/medium SHGC to reduce
transmitted SW solar radiation. Perforated and undulated exterior screen/louvers and motorized
rollers reduce peaks of SW solar radiation transmission up to 2% of incident solar radiation,
while SW solar transmissions can be as low as 8%.
Architectural design of CFS can extremely influence the effectiveness of CFS to control solar
heat gains. When CFS are not arranged properly, their effectiveness to control solar heat gains
can be significantly diminished. Therefore, it is important that designers take care how CFS are
installed to maximize their effectiveness to control solar heat gains.
Since high internal gains in Chilean offices buildings due to high occupancy ratios, significant
SW solar radiation transmission during wintertime might cause the need for cooling even in
winter. In consequence, designers should take into consideration the implementation of CFS to
control solar heat gains over the whole-year.
References
Appelfeld, D., A. McNeil & S. Svendsen. 2012. An hourly based performance comparison of integrated micro-
structural perforated shading screen with standard shading systems. Energy and Buildings, 50: 166-176.
ASHRAE. 2005. Handbook Fundamentals. Atlanta: American Society of Heating, Refrigerating and Air-
Conditioning Engineers, Inc.
Khun, T.E., S. Herkel, F. Frontini, P. Strachan & G. Kokogionnakis. 2011. Solar control. A general method for
modeling solar gains through complex façades in building simulation programs. Energy and Buildings, 43
(1): 19-27.
Khun, T.E. 2006. Solar control: A general evaluation method for façades with venetian blinds or other solar
controls systems. Energy and Buildings, 38 (6): 648-660.
Li, D.H.W. & J.C. Lam. 2000. Solar heat gain factors and the implications to building designs in subtropical
regions. Energy and Buildings, 32 (1): 47-55.
Reilly, S. & W. Hawthorne. 1998. The impact of windows on residential energy use. ASHRAE Transactions, 104
(2): 791-798.
Winkelmann, F.C. 2001. Modeling windows in EnergyPlus. Proceedings of international IBPSA Conference,
Building Simulation, Rio de Janeiro, august 13-15, pp. 457-464.
SUMMARY
Energy efficient buildings are an essential factor to reduce the energy consumption by 2020. New
buildings have to meet severe requirements, whereas older buildings need renovation to reduce the
heat losses through the building envelope. Infrared thermography (IRT) might be an improvement
over existing methods to assess the thermal performance of an existing wall in a non-destructive way,
or to check upon the as-built quality, specifically in the case of window glazing. The technique
instantly visualises the surface temperature of a whole building part, and in turn might allow to
deduce the thermal transmittance accordingly. However, many parameters can influence the surface
temperature and lead to distorted conclusions. This paper reports on the impact of different indoor
and outdoor boundary conditions for the assessment of the U-value of glazing, using the results from
a numerical simulation model. After an analysis of 6 types of windows, it is concluded that for
specific conditions, IRT might allow to estimate the U-value with an acceptable accuracy, based
upon the instantaneous indoor surface temperature. In future research, experiments will be
performed to validate this conclusion and the assessment methodology will be improved.
1. Introduction
Energy efficiency and sustainability are major concerns of our time. In 2011, the European
Commission decided that the energy consumption in general should decrease with 20% by 2020
(compared to 1990) (EC, 2011). To reach this goal, the building industry has to comply with stringent
regulations in respect to insulating performance, airtightness quality, renewable energy sources, etc.
However, new buildings only constitute a minority of the building stock. In Belgium, 62 % of the
buildings have been constructed before 1970 (WTCB, 2005), and typically do not comprise any
insulation in the building envelope. Consequently, deep energy renovation and refurbishment has
become essential to hit the target by 2020.
To verify whether the thermal performance of a building envelope meets the requirements, or to
determine where renovation is necessary, the as-built state has to be evaluated. To this end, it is
generally accepted to consider the thermal transmittance (U-value). This value can be calculated for
steady-state conditions from technical standards (EN673, 1997, ISO6946, 2007) or be measured on
site by means of heat flux sensors (ISO9869(E), 1994). The latter technique determines the U-value of
a building element by measuring the heat flux on one side and the surface temperature on both sides.
It is important to collect data during a considerable time span (preferably about 2 weeks) as in reality
the boundary conditions are always fluctuating. With a substantial dataset it is typically possible to
derive a precise U-value. This method is time consuming, and in principle different points on the wall
should be analysed to exclude singularities, which renders this approach rather elaborate. Infrared
The heat flux can also be expressed in function of the indoor-outdoor air temperature difference,
Q = U ⋅ (θ i − θ e ) (3)
By definition, Eq. 4 is only valid under static boundary conditions, but it will be applied to transient
conditions to assess to what extent it might provide useful information. Since the value for he is
function of the strongly varying wind velocity (Emmel et al., 2007), this study emphasises on the
determination of the U-value from the interior side. A sensitivity analysis of θsi to variations in θi is
performed by subjecting six types of windows to 3 different types of indoor temperature regimes.
Additionally, the influence of the most critical outdoor parameters is analysed. Note that specific
constraints in respect to thermography on glazing units are not addressed here. It is assumed that the
surface temperature can be determined accurately with an IR camera. Evidently, uncertainties in
emissivity and background temperature will propagate in the uncertainty interval of the estimated U-
value, as well as the accuracy of the camera itself.
2.1 Cases
FIG. 2 illustrates the different glazing assemblies that are studied. These configurations are
representative of those found in the Belgian building stock, and the broad variety of thermal
performance levels ensures the wide applicability of the results. The properties of each layer are
FIG. 2 Overview of the 6 investigated glazing types. Diversity in glass thickness, cavity fill and
coating broadens the scope of the analysis.
FIG. 3 Variable outdoor climate data for the period of December17-26, Uccle, Belgium, and 3
different temperature regimes for the indoor air (0= Dec 17 0h00, 10= Dec26 24h00 )
FIG. 4 U-value assessment based on θsi for 3 θi- regimes and 6 glazing types. Note: in some cases the
estimated U-value became negative, but these values are excluded from the graph.
For each indoor temperature regime, the U-value for the 6 types of glazing is estimated every 10
minutes, based on the internal surface temperature, using Equation 4 with hi equal to 7.7 W/(m².K). In
clear overcast
FIG. 5 Influence of the boundary conditions for the assessment of the U-value, illustrated for window
type 3. During sunshine hours and clear skies, the deviation to the theoretical value increases.
FIG. 6 Distribution of the assessed U-value(based on θsi ) for 3 θi - regimes and 6 glazing types. The
results are either only excluding sunny periods (no sun) or also moments of clear sky(cloudy) .
To increase the reliability of the estimation, events that cause extreme aberrations should be excluded
from the data. From FIG. 5 it is clear that the largest deviations of the estimated value (Utsi) come to
4. Conclusions
This paper discusses the application of quantitative infrared thermography on window glazing. Based
on numerical simulations for 6 window assemblies with different thermal performances, a sensitivity
analysis has pointed out the most critical restrictions to the boundary conditions for a reliable
assessment of the U-value from the interior surface temperature. Assuming that this temperature can
be determined correctly, the limited simulation results show that the thermal transmittance can be
estimated within a minor confidence interval of 0.65 W/(m².K) when direct sunlight and clear skies
are avoided (except for single glazing). The regime of the indoor air temperature turned out to be of
minor importance for the quality of the results. These findings are promising for the application of IR-
thermography in assessing the insulation quality of glazing, but further research is needed to extend
them to other types of windows and other periods of the year. Of course, the confidence interval will
enlarge for on-site measurements because θsi can only be determined within a certain range, e.g. due to
the accuracy of the camera, the variability in emissivity, and the influence of the background. Future
research will extend the error estimation by means of experiments, to verify the presented conclusions
and include additional noise. Furthermore, the use of external surface temperatures for the assess-
ment of the U-value will be investigated. Up to now, this data is ignored in the analysis, but including
them might improve the methodology as information of both sides of the assembly is used. Also,
consecutive measurements in a short time span may add to the accuracy of the U-value calculation.
References
ASDRUBALI, F., BALDINELLI, G. & BIANCHI, F. 2012. A quantitative methodology to evaluate
thermal bridges in buildings. Applied Energy, 97, 365-373.
BALARAS, C. A. & ARGIRIOU, A. A. 2002. Infrared thermography for building diagnostics. Energy
and Buildings, 34, 171-183.
BURN, K. & SCHUYLER, G. 1980. Applications of infrared thermography in locating and
identifying building faults. journal of international institute for Conservation 4, 3-14.
DALL'O, G., SARTO, L. & PANZA, A. 2013. Infrared Screening of Residential Buildings for Energy
Audit Purposes: Results of a Field Test. Energies, 6, 3859-3878.
EC 2011. Europe 2020 targets. EU EUROPEAN COMMISSION.
EMMEL, M. G., ABADIE, M. O. & MENDES, N. 2007. New external convective heat transfer
coefficient correlations for isolated low-rise buildings. Energy and Buildings, 39, 335-342.
EN673 1997. Glass in building. Determination of thermal transmittance (U value). Calculation
Method. European Standard EN.
FOKAIDES, P. A. & KALOGIROU, S. A. 2011. Application of I thermography for the determination
of the overall heat transfer coefficient in building envelopes. Applied Energy, 88, 4358-4365.
ISO6946 2007. Building components and building elements. Thermal resistance and thermal
transmittance - calculation method. ISO International Organization for Standardization
ISO9869(E) 1994. Thermal insulation - Building elements. In-situ measurement of thermal resistance
and thermal transmittance. International Organization for Standardization ISO.
ISO10456 2007. Building materials and products - Hygrothermal properties - Tabulated design values
and procedures for determining declared and design thermal values. ISO International
Organization for Standardization.
KALAMEES, T. 2007. Air tightness and air leakages of new lightweight single-family detached
houses in Estonia. Building and Environment, 42, 2369-2377.
LEHMANN, B. & GHAZI WAKILI, K. 2013. Effects of individual climatic parameters on the
infrared thermography of buildings. Applied Energy, 110, 29-43.
LUCCHI, E. 2011. Non-invasive method for investigating energy and environmental performances in
existing buildings. PLEA Conference on Passive and Low Energy Architecture. Belgium.
PHYSIBEL 2011. Voltra 7.0 Manual. Physibel Software.
STRAUBE, J. & BURNETT, E. F. P. 1999. Rain Control and Design Strategies. Journal of Building
Physics, 23, 41-56.
TAKI, A. H. & LOVEDAY, D. L. 1996. External convection coefficients for framed rectangular
elements on building facades. Energy and Buildings, 24, 147-154.
TAYLOR, T., COUNSELL, J. & GILL, S. 2013. Energy efficiency is more than skin deep: Improving
construction quality control in new-build housing using thermography. Energy and Buildings,
66, 222-231.
THOMAS, F., GLANZMANN, J., KELLER, B., QUEISSER, A. & RAGONESI, M. 1990. Element
29: Wärmeschutz im Hochba, Zürich, Faktor Verlag.
WTCB 2005. Visietekst 'Bouwen en Innoveren'- Bouwen aan de toekomst. Belgisch Bouwplatform.