Fire Safety Journal: Simo Hostikka, Rahul Kallada Janardhan, Umar Riaz, Topi Sikanen
Fire Safety Journal: Simo Hostikka, Rahul Kallada Janardhan, Umar Riaz, Topi Sikanen
Fire Safety Journal: Simo Hostikka, Rahul Kallada Janardhan, Umar Riaz, Topi Sikanen
A R T I C L E I N F O A BS T RAC T
Keywords: Fire-induced pressures have not been considered dangerous in building fires, but the situation may be changing
Pressure as building envelopes become increasingly air-tight. In this study, we investigate whether this can change the
Modelling fire development and pose new risks for structural and evacuation safety. We used experiments to validate the
CFD numerical models, and models for simulating the fire development in buildings with different air-tightness
Smoke
levels. The simulations considered air permeability values typical for traditional, modern and Near-Zero
Airtightness
buildings. Three different smoke damper configurations were studied, and the fire growth rates were varied
Near-zero buildings
High-rise buildings from medium to ultra-fast. The results showed that transitioning from traditional and modern buildings to
Near-Zero buildings can sufficiently increase the peak overpressures from fast-growing fires to cause structural
damage. Conditions were identified for avoiding excessively high overpressures, while preventing smoke from
spreading through the ventilation system.
⁎
Corresponding author.
http://dx.doi.org/10.1016/j.firesaf.2017.04.006
Received 15 February 2017; Received in revised form 27 March 2017; Accepted 3 April 2017
0379-7112/ © 2017 Elsevier Ltd. All rights reserved.
Please cite this article as: Hostikka, S., Fire Safety Journal (2017), http://dx.doi.org/10.1016/j.firesaf.2017.04.006
S. Hostikka et al. Fire Safety Journal xxx (xxxx) xxx–xxx
Nomenclature
S compartment surface area (m2)
Roman Symbols T temperature (K)
u duct flow speed (ms−1)
A crossectional area of the duct (m2) Vleak leakage flow (m3s−1)
AL leakage area (m2) V50̇ volumetric leakage flow at ΔP =50 Pa (m3s−1)
cp spesific heat capacity (kJ kg−1K−1)
Cd discharge coefficient (-) Greek
Δp pressure difference (Pa)
h enthalpy of fluid (W m −2K−1) ρ density (kgm−3)
K loss coefficient (-)
P pressure (Pa) subscripts
Q heat release rate (HRR) (kW)
tg growth time (s) i, k nodes
q50 leakage flow at ΔP =50 Pa (m3h−1m −2 ) j duct segments
The aim of the current research was to examine the influence of 2. Methods
building envelope air tightness on fire pressure and its consequences
in residential buildings. In particular, we intend to determine 2.1 Numerical method
whether new risks can be expected for structural and evacuation
safety and explore the possibility of using the building ventilation FDS is a Large Eddy Simulation (LES) based Computational Fluid
network for reducing pressure without increasing the risk of smoke Dynamics software which solves the low-Mach number combustion
spreading between apartments. First, we use existing [10,11] and equations on a rectilinear grid over time. A dedicated module for
new [12] experiments to validate the modelling capability of fire- modelling Heating, Ventilation and Air-conditioning (HVAC) systems
induced pressures in residential buildings. Next, the validated model is coupled with the gas phase solver. The ventilation network is
is used to perform numerical experiments in hypothetical apartment described as a series of ducts and nodes. The nodes are placed at
buildings with different levels of air-tightness. The severity of the points where ducts intersect each other or the CFD computational
predicted conditions is evaluated in light of two model outcomes: domain. The ducts represent uninterrupted domains of fluid flow which
pressure inside the apartment and smoke spread to neighbouring encompass elbows, expansion/contraction fittings and other fittings.
apartments through the ventilation network. The losses due to friction and duct fittings are assigned as dimension-
less loss numbers to the ducts. The node losses are also attached to the
Fig. 1. Plan drawings and measurements of the FOA (a) and Aalto [12] (b) validation experiments.
2
S. Hostikka et al. Fire Safety Journal xxx (xxxx) xxx–xxx
ducts as loss terms only appear in the duct Eq. (3). However, since the 31–33 a supply network. Schematic diagrams of the ventilation
module does not presently store any mass, mass flux into a duct is networks, simplified for modelling, are shown in Fig. 2. Test 23 had
equal to the mass flux out of the duct. The nodal conservation a slightly different exhaust ventilation network, as explained in [11].
equations for mass, energy and momentum equations are as follows: The ventilation networks consisted of a main duct connecting a fan to
the fire room, and additional branches connecting the main duct to
∑ ρj uj Aj =0
ambient through a damper. The purpose of the branches was to
j (1)
simulate additional rooms of actual ventilation systems. The dampers
∑ ρj uj Aj hj = 0 were used to limit the flow rate to these virtual rooms.
j (2) A detailed description of the Aalto experiments is given in [12], and
only a brief summary is presented here. The experiments are summar-
du ized in Table 2. They were conducted in a single apartment of a 1970's
ρj L j = (Pi − Pk ) + (ρgΔz )j + ΔPj + 0.5Kj ρj |uj |uj
dt (3) apartment building. The apartment (see Fig. 1b) included a living
room, bedroom, kitchen, bathroom, aisle and a closet, having a total
2.2 Leakage modelling floor area of 58.6 m2 and a ceiling height of 2.57 m. It was located on
the first floor of the building. The apartment had two exhaust
Leakage refers to the amount of air flowing through the building ventilation ducts leading to the roof. Three different ventilation
envelope due to the pressure difference between the inside and outside. configurations were created by keeping the original dampers
A standardized method for measuring the leakage (e.g. SFS EN: 13829) (Normal), removing them completely (Open) or closing the ducts
makes use of a variable-speed fan mounted to a door in the envelope tightly (Closed). All the other ventilation paths were closed. Two
being tested [13]. All the ventilation paths are sealed during the test. different fires were used: heptane pools with peak HRR of 800 kW
The building interior is then pressurized or depressurized using the fan were used in experiments 1–3 and flexible polyurethane foam with
and the leakage flows are determined by monitoring the flow through peak HRR of about 1 MW in test 4. Experiments 1–3 were repeated
the fan. Energy efficiency studies commonly measure leakages at an several times.
underpressure of 50 Pa. Overpressures would be more suitable for fire
analyses, though both directions can appear. The leakage can be 2.3 Simulation models
reported as a volumetric flow V50 ̇ , air permeability q50 or air exchange
rate n50. The simulation model of the FOA experiments consists of a single
From the viewpoint of fire CFD, since leakage is a sub grid -scale room with a leakage area connected to an additional duct. The main
phenomenon, a leak path boundary condition cannot be directly components of the HVAC models, shown in Fig. 2, are nodes, ducts,
specified. In FDS simulations, the leaks are modelled as a finite area and a fan. The other components such as dampers and expansion/
using vents on the domain boundaries. FDS utilizes the HVAC module contraction fittings are accounted for in the loss terms. Multiple ducts
to solve the leakage flow have been combined into a single duct with appropriate loss coeffi-
cients for further simplification. Dampers do not completely block the
⎛ ⎞1/2
̇ = AL sign(Δp ) ⎜ 2Δp ⎟
Vleak smoke flow, they only limit the volume flow to 25 l/s.
⎝ ρ ⎠ (4) The exhaust network model consists of 8 nodes, including a room
node connecting the main duct to the CFD domain, four ambient
nodes, and three internal nodes for connecting the branches to the
2.3 Validation studies
main duct. The total length of the main duct is 7.5 m. The exhaust fan
that drives the flow is placed in the duct connecting the main duct to
The FDS validation is performed by simulating two series of
the ambient. The fan curve was defined by specifying the flow rates (0,
compartment fire experiments that were conducted by the Swedish
60 and 120 l/s) at three static pressures (310, 190 and 18 Pa),
FOA Defence Research Establishment and Aalto University. The FOA
respectively. The layout of the supply network is similar to the exhaust
experiments were carried out in a room with dimensions of
network. The fan is obviously operating in different direction, and it is
4.0 m×5.5 m×2.6 m (height), shown in Fig. 1a. The fire room was
connected directly to the ambient without any duct segment.
divided in two parts with a wall, and the wall had a 1.9 m wide opening
In theory, it would be possible to determine the loss coefficients for
from floor to ceiling. The fire source was a heptane pan of 0.73 m×1 m.
A t2-type HRR curve was achieved using a lid that was moved over the
Table 1
pan at a given rate thereby increasing the heptane burning area. FOA experiment configurations, following the original test numbering in [10,11]. Test 15
Assuming a value of 1600 kW/m2 for the HRR per unit area, the had constant pool area of 0.5 m2.
movement speed of the lid was adjusted to yield three different growth
rates: Test Fire type Opening D (m) Exhaust network Supply network
FOA series 1
Type 1: Q̇ = 0.035t 2 kW (between fast and medium). FOA 1 1 0.2 m – –
Type 2: Q̇ = 0.075t 2 kW (between ultra-fast and fast). FOA 2 2 0.2 m – –
FOA 3 3 0.2 m – –
Type 2: Q̇ = 0.085t 2 kW (between ultra-fast and fast).
FOA series 2
The FOA experimental parameters are summarized in Table 1. Expe- FOA 11 1 0.2 – –
riments 1–3 [10] did not include an actual ventilation network, but the FOA 12 1 0.15 – –
fire room had a circular opening (D=0.2 m) at height 0.6 m from the FOA 13 2 0.2 – –
FOA 14 2 0.15 – –
floor, connected to a 2.2 m long and 0.2 m diameter tube. Temperature FOA 15 – 0.2 – –
and flow speed were measured at the end of the tube. FOA 21 2 0.2 Exhaust network 1 –
FOA experiments 11–33 [11] included three different ventilation FOA 22 2 0.15 Exhaust network 1 –
configurations. The leakage openings of different diameters were FOA 23 2 0.15 Exhaust network 2 –
FOA 31 1 0.10 Exhaust network 1 Supply network 1
connected to a 0.32 m diameter and 3.2 m long tube connecting to
FOA 32 2 0.15 Exhaust network 1 Supply network 1
the ambient. The openings were placed at height 0.6 m from the floor. FOA 33 2 0.2 Exhaust network 1 Supply network 1
Experiments 21–33 included an exhaust network, and experiments
3
S. Hostikka et al. Fire Safety Journal xxx (xxxx) xxx–xxx
all the duct sections from the individual losses of the components and
the friction. As all the details cannot be determined at this stage,
effective loss coefficients were estimated first by defining a constant
roughness of 1 mm for the ducts, and then adjusting the loss of the first
duct to match the known volume flow of 25 l/s. The losses from the
damper and the 90 degree bend were thus combined into a single loss
coefficient. Identical losses are used for forward and reverse flow.
For leakages, the bulk leakage method was used over all wall
surfaces. For other boundary conditions, the models of both FOA and
Aalto experiments included only the interiors of the enclosures,
assuming that the exterior sides of the walls are exposed to the
ambient. For the concrete boundaries, thermal conductivity, specific
heat and density were set to 0.7 W/(m. K), 0.75 kJ/(kg. K) and
2200 kg/m3, respectively. Fire boundary conditions were specified as
time dependent heat release rates per unit area (HRRPUA) of the fuel
pan surfaces. Technically, they were converted to time dependent mass
Fig. 2. The ventilation networks of FOA experiments as modelled in FDS. Top: Exhaust flux and species concentration boundary condition of the CFD solver.
network 1. Bottom: Supply network 1. The simulations were carried out using the FDS version 6.3.2 with a
single mesh and 6 OpenMP threads on a personal computer with a 3.2-
GHz Intel Xeon processor and 32-GB ram. The CPU time for a single
Table 2 simulation of 300 s was approximately 8 h.
Parameters of the Aalto experiments. Different numbering from [12] is used here for
simplicity.
3. Validation results
Test Fuel Ducts Roof fan
The experimental and simulated fire room pressures in FOA
Aalto 1 Heptane 3.0 L Open OFF experiments 1–3 are shown in Fig. 3. The positive pressures during
Aalto 2 Heptane 3.0 L Normal ON
the fire growth stage are reproduced by the simulation model with good
Aalto 3 Heptane 3.0 L Closed ON
Aalto 4 Polyurethane foam 3.82 kg Normal ON accuracy. The negative pressures after the fire suppression, in turn, are
not captured as well. The same behaviour was observed in all the
validations tests of this work.
A summary of the gas temperature and peak overpressure predic-
4
S. Hostikka et al. Fire Safety Journal xxx (xxxx) xxx–xxx
5
S. Hostikka et al. Fire Safety Journal xxx (xxxx) xxx–xxx
The influence of the fire/smoke dampers installed in the ventilation The predicted peak pressures are summarized in Fig. 7. These
ducts was studied by analyzing three different compartmentation results were not corrected for the estimated model bias of - 7%. The
damper configurations: trends in the results are clear and consistent. All three parameters - fire
growth rate, damper configuration and air-tightness - are found to be
Damper=Off Both inlet and outlet remain open during the fire. important for the expected peak pressure. Interestingly, the sensitivity
Damper=Inlet The inlet duct of the fire apartment is closed by a da- of the pressure to the parameter values seems to increase when moving
mper 10 s after the ignition. towards a scenario with higher pressure. For instance, the damper
Damper=Both Both inlet and outlet are closed by dampers 10 s from configuration is not very important in traditional or normal buildings,
the ignition. but can become crucially important in very air-tight buildings. The
results in Fig. 7 were found to be independent of the fan operation (on
Additionally the effect of the dampers located at the inlet and outlet or off) and position of the fan unit damper. Therefore, the leakages
fans was investigated. The ventilation network parameters are shown through the other apartments can compensate for the fan pressure
in Table 3. differences and complete closing of the ventilation system. Of course,
Three different levels of the building envelope air-tightness were this leads to smoke spreading to the other apartments through the
examined. These levels were defined using the air permeability values network.
q50, listed in Table 4. The class “Traditional” represents an average of The pressure results can be compared against the simulations of
the required and reference air-tightness value (for heat loss calcula- Fourneau [4] who assumed a fire with a growth rate between medium
tions) described in the current Finnish building code (Part D3: Energy and fast, more airtight envelope (n50 ≈ 0.6 h−1), and two open paths
efficiency, 2012, Ministry of Environment). The “Modern”, corresponds through the mechanical ventilation system. The peak pressure in their
to the measured air-tightness in the concrete element multistorey simulation was about 550 Pa, i.e. between our results obtained for the
buildings [13] and the “Near-zero” represents the current, technically Damper=Off configuration.
achievable target level. In order to evaluate of the occupants’ possibility of escaping from
The volumetric leakage flow rate V50̇ at 50 Pa can be calculated from the fire apartment, we compare the overpressure against an appro-
the air permeability values using Eq. (5) priate threshold. Here, we assume that opening the door would not be
possible if the overpressure were above 100 Pa. For design purposes
q50
V50̇ = S the difference in pressure is usually limited to 50 Pa, but we choose a
3600 (5)
less conservative value for the risk analysis. Based on the predicted
and the air exchange rate as n50 = V50 / V , where V is the building peak pressures, it should be possible to open the door in traditional
volume. The leakage areas AL , through which these air-tightness levels buildings at a medium fire growth rate. For more air-tight buildings
are specified in the FDS models, are then calculated from the and faster fires, opening the door would be challenging. However, it is
volumetric flow rates as important to note that these correspond to momentary peak pressures,
and reveal little about the time duration of the peak pressure.
V50̇
AL = Challenges for the structural integrity can be expected in fast or
⎛ 2Δp ⎞1/2
Cd ⎜ ρ ⎟ ultra-fast fires. If the criterion of failure is chosen based on our own
⎝ ∞⎠ (6) experimental observation (1500 Pa), the fast fires could pose a risk
when the dampers are closed and the envelope is very airtight. This
In the models, these leakage areas are distributed equally at the
could occur in near-zero or high-rise buildings. In ultra-fast fires,
locations of doors and windows.
problems could be expected in all modern buildings using dampers. In
Three different fire scenarios were used in the simulations. Two of
near-zero buildings, the capacity of the ventilation network is insuffi-
them were defined as typical t2 fires with prescribed maximum HRR,
cient to relieve the pressure of an ultra-fast fire even without dampers.
and the third one was based on the experimental heptane pool burning
In such a scenario, alternate means for relieving pressure should be
rate in the Aalto experiments.
investigated.
6
S. Hostikka et al. Fire Safety Journal xxx (xxxx) xxx–xxx
6. Conclusions
7
S. Hostikka et al. Fire Safety Journal xxx (xxxx) xxx–xxx
tightness (reduced leakage), increased use of fire or smoke dampers, exist at fast and ultra-fast fire growth rates.
and the increased rate of fire growth. The 100 Pa limit was reached in Simulations of different damper configurations and fan operation
all the simulations except those made for the traditional buildings with modes showed that the spreading of smoke to neighbouring apart-
medium fire growth rate. In a Near-Zero building, the limit was ments can be avoided if only the inlet ventilation branch is closed with
reached in less than 20 s for ultra-fast fire and in about 100 s for a a damper and exhaust fan is operating. This mode of operation was
medium fire. The duration varied with respect to the fire type, being in found to be the only combination for preventing smoke spread, while
between 100 and 200 s. The risk of structural damage was found to simultaneously maintaining the pressure at an acceptable level.
8
S. Hostikka et al. Fire Safety Journal xxx (xxxx) xxx–xxx