FIRE MODELLING OF ENERGY-EFFICIENT APARTMENT BUILDINGS
Consideration of air-tightness and mechanical ventilation
Simo Hostikka
Aalto University
Rahul Kallada Janardhan
Aalto University
Umar Riaz
Aalto University
Topi Sikanen
VTT Technical Research Centre of Finland Ltd.

Abstract. In this work, we wanted to study how the construction trends aiming at energy efficient and high-rise buildings are changing the fire modelling practices. Through experiments, FDS validation and a simulation case-study, we investigate the reliability and modelling practices of the mechanical ventilation systems and air-tight building envelopes. The simulation results indicate that the new, very air-tight building envelopes can pose a risk for both occupant and structural safety in fires.

1. INTRODUCTION

In the modeling of enclosure fires, it has been common to assume that the buildings that we analyze are not air-tight, with a few exceptions, such as the special industrial facilities like nuclear power plants. Pressure build-up in the simulations has been avoided because (i) it has been deemed sufficient to focus on the hydrostatic pressure driving the flows in the space, (ii) accurate information about the air-tightness has not been generally available, and (iii) the means of calculating the leakage flows have been difficult to integrate with the numerical fire models. For long, it was a common practice to define a 10 by 10 cm leakage hole to the rooms when modeling them with zone or CFD models. With the current trends of energy-efficiency and high-rise buildings, the situation has changed. Both trends pose much more strict requirements for the building envelope's air-tightness than what has been the common practice before. As a result, the development of over- and under-pressures has become an essential feature of building fires and their simulations.

At the same time with building envelopes becoming more air-tight, the ventilation systems are getting more complicated. Meeting the energy-efficiency requirements necessitates the use of mechanical supply and exhaust ventilation with heat recovery. Integration with building services and active fire protection technologies makes the systems complicated and challenging for modeling.

In FDS, the technical capability to simulate HVAC systems and leakages was introduced by the implementation of a dedicated HVAC module [4].

2. FEATURES OF MODERN BUILDINGS

2.1. Ventilation systems

The traditional ventilation systems in residential buildings have been based on mechanical or buoyancy driven exhaust, with the supply air provided through the building envelope either as an uncontrolled leakage or through valves. In modern HVAC systems, there are usually separate networks for supply and exhaust air. Both networks are typically equipped with a fan unit to control the flow rate and to implement the heating/cooling, as well as heat recovery and air filtration for indoor air quality control. Single fan unit may serve an entire building or a single floor of a multi-storey apartment building. The coils and filters introduce drag to the flow. As a result, the fan unit will cause pressure losses even when the fans are turned off. For fire compartmentation and smoke control, fire and/or smoke dampers are typically installed to the ducts entering and leaving the apartments. In addition, the modern ventilation fans are often equipped with dampers that automatically close the ambient connection when the fan is turned off. In fire situation, turning off the fan can therefore lead to complete closing of the ventilation system.

In principle, the ventilation system could be modelled starting from the drawing describing the design or ‘as built’ status. This is possible if we know the exact dimensions and characteristics of the various components of the system. The losses were determined using based on individual pressure loss of each duct fitting and the friction loss of the duct. The volume flow or velocity of flow through the ducts is required to find the loss coefficient $K_{j}$ appearing in the duct momentum equation

(1)
\[\rho_jL_j\frac{du}{dt} = (P_i - P_k) + (\rho g \Delta{z})_j + \Delta{P_j} + 0.5K_j\rho_j|u_j|u_j \]

For example, assuming a duct of length 1m and diameter of 0.1m with an expansion fitting which leads to a duct of 0.125m and a damper regulating the flow to 25 L/s. The loss coefficient for this duct is calculated as follows:

This kind of process requires a lot of information that can be difficult or impossible to obtain for a fire engineer. Additionally, the process is sensitive to errors.

Practical experiences from the two validation studies and two case studies has shown us that the practical HVAC modelling procedure in the near future will go along the following lines:

  1. Gather the basic design information about flow rates, known pressure losses (HVAC valves and fan unit components), and the fan operating characteristics (fan curve).
  2. Model the HVAC network in the simplest possible way, maintaining the physical dimensions (diameters, lengths and heights) of the system.
  3. Assign some realistic value for the duct roughness to introduce main duct pressure losses.
  4. Specify additional losses to the ducts and adjust them to reach the expected flow rates through the system.
  5. Verify the pressure levels and flow rates in a normal state by adding measurements for the flows and node pressures.

So, the actual system information is limited to the design flow rates and the fan characteristics,as the fan pressure levels dominate the system's performance.

2.2. Air-tight envelopes

The building envelope air-tightness is important for many building performance aspects, including the energy-efficiency of buildings and for the control of indoor environment and ventilation in high-rise buildings. In energy-efficiency, the goal is to reduce the energy losses by convection through the cracks and gaps of the envelope. The convection is driven by the pressure difference between the interior and exterior, which in small houses can vary from few tens of Pa underpressure to a few Pa overpressure. For the moisture control, it is, at least in Finland, common to maintain the buildings at small negative pressure. In high-rise buildings, the pressure differences can be much higher.

The envelope air-tightness is commonly measured using a blower-door test [1] where a powerful fan is attached to a door leading outside (see Fig. 1), and the other vents of the building are closed tightly. The fan is used to create a pressure difference between the interior and the ambient. The amount of air flow through the fan is measured at several pressure differences. The normal practice in the energy efficiency studies is to report the leakage rate at 50 Pa underpressure in a form of a volumetric flow rate $\dot{V}_{50}$, air permeability $q_{50}=\dot{V}_{50}/A_{\rm env}$ or air exchange rate $n_{50}=\dot{V}_{50}/V$. Here, $A_{env}$ is the envelope area and $V$ is the building volume.

blower_door Closed_vent


Figure 1. Blower-door test (left) and sealed exhaust ventilation duct (right). (Photographs by Mikko Yli-Piipari, Vertia Oy).

FDS utilizes the HVAC module to solve the leakage flow

(2)
\[ \dot{V}_{leak} = A_{L} \mathrm{sign}(\Delta p) \left(\frac{2\Delta p}{\rho} \right)^{1/2} \label{eq_leakage_FDS} \]

The main input for the leakage modelling is therefore the effective leakage area, which can now be determined from the air-tightness measurements as

(3)
\[A_{\rm L} = \frac{\dot{V_{50}}}{C_d \left(\frac{2\times 50}{\rho_{\infty}} \right)^{1/2}} \; {\rm [m^2]} \]

Of course, the area can be calculated using the leakage flow at some other pressure than 50 Pa, if available.

3. VALIDATION

3.1. FOA Test series

The Swedish FOA Defence Research Establishment conducted two sets of experiments in the late 1990's to investigate the fire pressure and duct flows [2, 3]. The first series [2] consists of three tests with $t^2$ fires of different growth rates between medium and ultra-fast. The fire room was 4.0 m $\times$ 5.5 m $\times$ 2.6 m (high) and the fire source was a heptane pan of 0.73 m $\times$ 1 m (c.f. Fig. 2). These tests did not include actual ventilation network, but the fire room had a circular opening ($D$ = 0.2 m) connected to a 2.2 m long tube. Temperature and flow speed were measured at the end of the tube. The opening was located at 0.6 m from the floor. The fire room was divided in two parts with a wall, and the wall had a 1.9 m wide opening from floor to ceiling. The $t^2$ behaviour of the HRR was achieved using a lid that was moved over the pan at a given rate thereby increasing the heptane burning area. The quantitative value of HRR was obtained by assuming 1600 kW/m$^2$ for the HRR per unit area of pool surface.


FOA_Geometry


Figure 2. FOA test geometry.

The second series [3] consists of three groups of experiments with the same room and same fire types as the first set, but with different ventilation configurations. The leakage openings of different diameters were connected to a $0.32$ m diameter and $3.2$ m long tube connecting to the ambient. The opening was located at $0.6$m height from the floor on one wall of the room. In the first group, there was not ventilation system, only leakage opening. In the second group, the room was equipped with an exhaust system, shown in the top part Fig.3. A supply network (Fig. 3 lower part) was added in the third group.

The HVAC model components are nodes, ducts, and a fan. The other components such as dampers and expansion/contraction fittings are accounted for in the loss terms. Multiple ducts have been combined into a single duct with appropriate loss coefficients for further simplification. Dampers only limit the volume flow through the ducts to 25 l/s.

The exhaust network was modelled using 8 nodes, starting from Node 1 which connects the computational domain to the main duct which is also connected to three fictive compartments maintained at ambient pressure. The total length of the duct from Node 1 to 5 was about 7.5 m. Nodes 5,6,7 and 8 are maintained at ambient pressure. The exhaust fan that drives the flow is placed in the last duct connecting nodes 4 and 5. The fan curve was defined by specifying the flow rates (0, 60 and 120 Ls) at three static pressures (310, 190 and 18 Pa), respectively.


Exhaust1 Supply


Figure 3. FOA exhaust and supply ventilation systems.

3.2. Aalto University test series

The third series of validation experiments consists of experiments carried out in 2015 in a 1970's apartment building in western Finland. A detailed explanation of the experiments is given in [5]. The apartment (see Fig. 4) had two exhaust ventilation ducts leading to the roof. All the other ventilation paths were closed during the experiments. The originally natural ventilation had been enforced by post-installed fan in the bathroom exhaust. Two different fires (Fig. 5) and three exhaust damper configurations (Fig. 6) were used. Open configuration means that the dampers were removed completely, normal configuration that the original dampers were in place, and closed configuration that the exhaust ducts were tightly sealed.

Kurikka_geometry


Figure 4. Geometry of the Aalto experiments.

ignition foam


Figure 5. Heptane pool (top) and Polyurethane foam fires.

The ventilation ducts were modelled as combinations of two or three duct segments. The first duct was installed for the measurements, and the third segment of the bathroom duct was used for including the fan.

Ventilation


Figure 6. Ventilation system configurations. Open, normal, and closed.

3.3. Validation results

The experimental and simulated fire room pressures in FOA Series 1 are shown in Fig. 7. The positive pressures during the fire growth stage are reproduced by the simulation model with good accuracy. The negative pressures after the fire suppression, in turn, are not captured as well. The same behaviour was observed in all three validation series. Figure 8 shows the corresponding comparisons for the heptane pool tests of the Aalto University's series. The simulations were performed using both bulk and localized leakage method. The with the bulk method, the peak overpressures are overpredicted by 21 % but the underpressure peaks are captured at least qualitatively. With localized leakages, the peak pressure predictions varied from overprediction (Test 3) to underprediction for Test 8. The difference between the bulk and localized methods is at least partially explained by the different pressure value used for driving the flows: In bulk method, the background pressure drives the flow and in localized method the driving pressure is the local pressure at the vent, including the hydrostatic and perturbation components.

FOA_results


Figure 7. Measured and simulated gas pressures in FOA series 1.