key: cord-0045626-a255svyx
authors: Wąs, Jarosław; Karp, Artur; Łukasik, Szymon; Pałka, Dariusz
title: Modeling of Fire Spread Including Different Heat Transfer Mechanisms Using Cellular Automata
date: 2020-05-26
journal: Computational Science - ICCS 2020
DOI: 10.1007/978-3-030-50371-0_33
sha: 036cbd09cf8eaa90342c63fae61ee080276c6852
doc_id: 45626
cord_uid: a255svyx

The article presents a new method of modeling the spread of fire using Cellular Automata based on 3D mesh. The models based on different heat transport mechanisms were used in the study, including conductivity, convection and thermal radiation. General mechanisms of fire spread are exemplified by a fire inside a building. The tests of the created simulator for fire in rooms were conducted, and the reference point was the known simulator based on CFD, namely FDS Fire Dynamics Simulator.

Fire is a complex phenomenon which includes thermal degradation or pyrolysis of solid fuel into gas (of a volatile character) and combustion of the volatile fuel. It should be noted here that the generated heat influences the solid. Heat transport is carried out through thermal conductivity, convection (natural and forced) and thermal radiation. Modern computational sciences provide various methods and IT tools that allow for modeling this complex process [8] in a simplified way.

The most common models are high-fidelity CFD models, also known as Field models, like e.g. Fire Dynamics Simulator (FDS). based on the physical principles of conservation of: mass, energy, and momentum. In the field model, the motion of heat and smoke is calculated in particular cells in consecutive time intervals. In such models we assume that in each cell the properties like temperature, density, flow, etc. are constant. Another popular approach are Zone models, in which the whole space (for instance, a room) is divided into areas with different properties (like a 'combustion zone' or a 'convection zone'), and for each of them a set of algebraic equations is solved in order to calculate relevant flow properties. The most recent approach assumes looking for fire models that are realistic and computationally efficient.

The paradigm of Cellular Automata (CA) can be considered an effective and efficient solution. In the literature, CA are used to model simplified fire spreading schemes [1] , often in relation to evacuation [13] or specific scenarios like forest fires [6, 14] etc.

We propose a new model of fire spread based on Cellular Automata in which different mechanisms of heat transfer are taken into account, namely: conductivity, convection and thermal radiation. We propose the application of the cuboid structure of space using three-dimensional CA.

We developed a theoretical model, implemented it, and performed preliminary validation based on widely used CFD model FDS, created by NIST (National Institute of Standards and Technology), which is often applied as a reference model [4, 10] 2 Heat Transfer Mechanisms

The amount of energy which is transferred in thermal conduction can be computed using the law of heat conduction (Fourier's law). The differential form of this law is as follows:

The heat flux that can be emitted by radiation coming from the black body surface is given by Stefan-Boltzmann law:

where: q -is radiative heat flux (SI unit [ W

, T -is the temperature of the body surface (SI unit [K])

The flux emitted by the real surface is lower than the flux emitted by the black body surface:

where: -a factor called emissivity, which is a property characterising a surface ( < 1)

Convective heat transfer is transfer of heat from one place to another caused by the movement of fluid. The rate equation for this form of heat transfer is described by Newton's law of cooling:

where: q -is convective heat flux (SI unit [ W m 2 ]), h -is convection heat transfer coefficient (SI unit [ W m 2 K ]), (T w − T f ) -is the temperature difference between the surface of an object and fluid (SI unit [K]).

The assumptions of the model used in the implementation of the CAFire programme which simulates fire spread are presented below.

The whole space under consideration can be divided into three-dimensional elements (e.g. cubes) small enough to be treated as homogeneous -i.e. all physical and chemical properties (e.g. temperature, specific heat capacity) in a given element are the same at a given time.

Three basic types of neighborhoods exist in a three-dimensional array of a cubic cell:

-6 cells sharing faces with the reference cell -18 cells sharing edges with the reference cell -26 cells sharing vertices with the reference cell Because the simulation results obtained with the use of all the above types of neighborhoods are comparable and the first variant has the lowest calculation costs, it is used in the model (Fig. 1 ). 

Thermal Conduction. At a given temperature gradient, the heat transfer rate is proportional to the surface through which the heat flows. Generally, the amount of heat dQ that will flow L over dt with the temperature difference dT across the transverse surface S is as follows:

where: Q -heat, λ -thermal conductivity coefficient, S -transverse surface of the cell, L -distance between cells, dt -time, dT -temperature difference.

In the proposed model, the cells have the same surface area. Consider the following one-dimensional situation shown in Fig. 2 .

-warmer cell -cooler cell In this type of analyses, it is assumed that the heat exchange occurs between the cells, and so the λ factor coefficient indicates how the heat passes from one cell to another.

With two cells only one factor is determined. However, if a whole cell chain is considered, it may happen that for a given cell n the conduction coefficient from/to cell (n − 1) is different than from/to (n + 1).

In the discrete one-dimensional case, the formula for temperature T of the cell with index n and of a unit size can be written as: (6) where: λ = λ m·c , T n (t) -Cell n temperature at a given time t, λ -thermal conductivity coefficient, c -specific heat capacity, m -cell mass.

The phenomenon of convection is characteristic of liquids and gases. The model includes natural convection.

where: α = α m·c , T u (t) -the temperature of the upper cell at a given time t, T d (t) -the temperature of the lower cell at a given time t, α -convection heat transfer coefficient, m -cell mass, c -specific heat capacity (Fig. 3 ).

-upper cell -lower cell 

Thermal radiation is a process in which energy is emitted by a heated surface in all directions in the form of electromagnetic waves (electromagnetic radiation).

In the fire safety literature, thermal radiation has been identified as important and even dominant in terms of heat transfer in medium to large fires. Its mechanisms are analysed at various levels of accuracy ranging from simplified analyses to sophisticated methods of calculating radiation [9] .

Stefan-Boltzmann law refers to radiation in all directions, whereas in the model only some of the energy radiated by one surface will reach another. The heat exchange from a warmer (1) to a cooler (2) surface depends on the following factors [5] :

-temperatures T 1 and T 2 -surface areas A 1 and A 2 -spatial configurations between surfaces -surface radiation characteristics -additional surfaces in the environment -medium between surfaces (if it is air, it can be assumed as non-invasive)

To calculate the amount of energy that will reach from one surface to another, we need to calculate the integral, taking into account the density of the radiation at an angle under which fragments of the radiating surface are visible for the other surface. Then, an additional factor view factor appears in Stefan-Boltzmann law.

where: Q net -total radiation energy, F 1−2 -view factor, A 1 -surface area, σ -Stefan -Boltzmann constant, T 1 , T 2 -temperature in Kelvins. The building is divided into radiation areas (for example, walls or tables) with a direction from which they can receive heat (Figs. 4 and 5) . 

then:

where: F 1−2 -view factor, a -length of surface, b -width of surface, cdistance between surfaces (Fig. 6) . -Perpendicular surfaces let:

then:

where: F 1−2 -view factor, h -width of the first surface, w -width of the second surface, l -length of the surface For each wall of a fire cell parallel radiation surfaces are considered where no fire cells or solid body cells are located between them in a straight line.

T surf ace was calculated as the average temperature of the cells contained in this surface. Of course, the following condition must be met: T cell > T surf ace .

As an obstacle, we took into account a cell which was a solid body or a fire cell. Knowing surface dimensions a i b and the distance from the cell c coefficient, the view factor was calculated on the basis of Table 1 .

The last step was to calculate the transferred energy Q net according to Eq. 8. Energy is transferred equally to all cells belonging to the radiation surface.

In case of radiation surface ⇒ surface, the main radiation surfaces were determined at the beginning of the simulation; they can be generally treated as walls (in the model for each room, 6 such surfaces were designated for the floor, ceiling and each wall). For each surface other surfaces in its "field of visions" and its view factors were determined on the basis of Eq. 10 or 12. The same procedure was repeated in case of cell ⇒ surface.

The proposed model takes into account the properties of relevant materials, because each of them has a different autoignition and ignition temperature. Ignition is possible if a cell receives temperature proper to a given kind of material (e.g. wood 250 • C, carpet 100 • C, fabrics 170 • C) and one of the neighboring cells remains on fire long enough (Fig. 7) .

The only condition for autoignition is to reach a sufficiently high temperature (e.g. wood 2000 • C, carpet 500 • C, fabrics 650 • C).

The linear fire spread rate was calculated for each material. By transforming the formula to take into account the speed, we calculated how long the cell must burn (to be able to ignite its neighbor) (Fig. 8 ).

set in the cell "fire" (red) state. 

In the presented model smoke is generated in the cell located directly above the flame, which is shown in Fig. 9 . The amount of smoke emitted depends on the type of material that is burning. For example, less smoke will be released from a wooden cell than from a textile cell.

-cell with smoke -cell with fire The principle of mass conservation and the convection mechanism were retained in smoke modeling.

It was assumed that every cell tends to get rid of all its smoke. To ensure that the cell is not filled above its maximum capacity, the amount of smoke that it can take from one neighbor was limited to 1 6 of the remaining free space [12] . In order to ensure that the cell does not transmit more smoke than it has, the amount of smoke that a given cell can transmit was limited to 1 6 of the amount of smoke.

The amount of smoke in the cell will be reduced by the amount the cell will pass on to its neighbors and increased by the amount it will receive from its neighbors ( Fig. 10) :

where:

D 0→n -the amount of smoke passed to a neighbouring cell D n→0 -the amount of smoke received from a neighbouring cell d n -quantitative factor, according to Table 2 n -neighbor index When the cell above C U is not air, or the amount of smoke in it D U is close to 100% (that is, it cannot accept more smoke), the d n factor will be modified (see Table 3 ).

This modification adequately demonstrates pushing smoke out of the ceiling zones.

While calibrating fire scenarios, we took into account bibliographical data: standard PN-EN ISO 6946 and [3] , as well as some data from reference fires [2] . The simulation was based on materials that are usually found in a typical household - Fig. 11 (Table 4) . of fire and smoke was visually compared, while in case of FDS, the distribution of the temperatures in selected fire areas was additionally validated. Fire Dynamics Simulator [3, 4, 7, 11] is an application that, according to many experts in fire engineering, best reflects the course of fire in buildings. In order to compare our CA-based model with the same model in FDS, we compared two identical scenarios of a building fire. The size of the building: 7.5 m × 2.9 m × 7.5 m, including a living-room, a bathroom, a wardrobe, a study, a kitchen and a hall. Ignition temperature was 150 • C. It should be noted that other data sets were used for the calibration and validation procedures. Both, fire development scenario and basic parameters are similar in our model and the reference FDS model. Sample comparison results are presented in Fig. 12 and 13 . 

The main purpose of this study was to create an effective and efficient method for simulating the spread of fire. The model and simulations were based on nonhomogeneous Cellular Automata to present the distribution of fire and smoke in a residential building, taking into account the temperature fields in the rooms. The current model of fire called CAFire includes three main phenomena accompanying the spread of heat during a fire, i.e. conductivity, convection and thermal radiation. The simulations showed relatively reliable convergence of the results compared to the reputable FDS tool.

In the presented model, only in a very simplified form, the fluid flow phenomena (air flow, smoke flow) are simulated. The interactions between the air streamlines (due to Bernoulli's principle) and the influence of pressure caused by the difference of temperatures on the formation of laminar and turbulent flows etc. are not taken into account.

The idea presented in the paper requires further development and testsand additional emphasis should be placed on the chemical aspects of fire and exhaustive validation of the final results.

However, it should be emphasized that we consider the proposed direction of work very promising. Creating simulations in the FDS program requires time, even up to several dozen minutes, while calculations obtained using the proposed Cellular Automata method take seconds to be completed.

Cellular automata based simulation for smoke and fire spreading in large buildings

Living room fires

SFPE Handbook of Fire Protection Engineering

Risk assessment of main control board fire using fire dynamics simulator

A heat transfer textbook

Spread vector induced cellular automata model for real-time crown fire behavior simulation

Modelling real fire by FDS and 2-zone model for structural post-fire assessment

A comparison of three fire models in the simulation of accidental fires

Radiative heat transfer in fire safety science, radiation IV proc

Zone computer fire models for enclosures

Assessing the influence of fuel geometrical shape on fire dynamics simulator (FDS) predictions for a large-scale heavy goods vehicle tunnel fire experiment. Case Stud. Fire Saf

Cellular automata model for simulation of effect of guiders and visibility range

Evacuation dynamics with fire spreading based on cellular automaton

Forest fire spread simulating model using cellular automaton with extreme learning machine