GoZone: A Numerical Model for Travelling Fires Based on Cellular Automata Concept

: Fires in large compartments tend to burn locally and to move across the ﬂoor over a period of time; this particular behaviour has been discovered to challenge the assumption of uniform gas temperature in the ﬁre compartment. Recent studies on ﬁres in large compartments have led to the now widely known concept of “travelling ﬁres”. Several models have been proposed to describe the evolution in time of travelling ﬁres. Although these models represented an innovative step in the ﬁeld of travelling ﬁres, the major drawbacks of these models can be found in the simpliﬁcation of ﬁre dynamics (constant spread rate, 1D imposed ﬁre path) and limited ﬁeld of application (rectangular based geometries). The purpose of this paper is to present a numerical model of travelling ﬁre. The model was based on an improved zone model combined with a cellular automata model. The software GoZone, in which the model was implemented, is intended to be a practical solution to analyse ﬁres in large compartments of potentially any shape. GoZone is aimed to describe the complex dynamics of the ﬁre from ignition to a phase of growing localised ﬁre that may eventually travel in the compartment, possibly followed by a ﬂashover. The main sub models comprising GoZone are presented. A comparison is given with the results of under ventilated ﬁre test 2 of the BST/FSR 1993 test series and with respect to the Vesel ì travelling ﬁre test is shown. GoZone shows a promising capacity to represent ﬁres in a large compartment in both air and fuel controlled ﬁre conditions.


Introduction
Knowledge about fire development in buildings, particularly the concept of compartment fire, was developed during the decades in between 1960 and 1990. An important contribution was made by Kawagoe [1] and thereafter by Thomas [2]. The first author gave a definition of compartment fire, and through observations was able to define the link between gas phase temperature, burning rate, and ventilation. The second author, by extending and formalizing the work of Kawagoe, provided a series of formulations to estimate the fire duration and temperature, the latter was supposed to be uniform in the compartment. From the 70s onward, research was focused on refining the knowledge acquired from Thomas and Kawagoe as well as trying to translate this into design methodologies with an emphasis on structural performance [3,4]. Between the early 90s and in the first decade of the new millennium, uniform burning assumptions have been challenged by experimental fire tests, which have shown that in compartments with larger floor areas, the temperature distribution is not uniform [5]. Combustion does not occur simultaneously through the whole floor of the compartment; it tends to be limited to some portions of the floor and the position of these burning zones changes over time. These fires have been labelled as "travelling fires" [6].
The discipline of structural fire engineering has seen dramatic progress since the early 2000s. Nowadays, it is possible to model the behaviour of complex structures subjected to fire by nonlinear finite element software [7,8]. However, when it comes to the representation • A physical environment, which is made of a discrete lattice of cells of approximately equal area, called grid, in two dimensional problems. The lattice defines the domain in which CA is applied and computed. The physical environment here is the fire compartment, the floor of which is described with triangular and/or quadrilateral cells in a structured or unstructured manner. In GoZone, the discretisation can be made through a graphical user interface based on the open source mesh generator GMSH [28]. The same discretisation is mirrored up and assigned to the lower plane of the smoke layer Figure 1. The versatile character of unstructured triangular meshes makes it possible to model any shape for the plan view of the compartment, while structures of quadrilateral meshes may be used in more regular geometries such as, for example, rectangular compartments.

•
The cell's states, where each cell composing the lattice has a proper state that can change in time. Change in state can be referred to a single cell or to all cells composing the lattice. The cell's state is here: burning or not. Moreover, cells are characterised with properties such as fuel load density, fire growth rate, rate of heat release density, and height of the fuel must be defined in GoZone. While the cell state can be viewed as an "on" or "off" state, a cell property is a given input of the model. • A local transition rule, which acts upon a cell and can generate a change of state from one discrete time step to the next. The transition rule enabling a change in state of a cell from not burning to burning is described by an ignition strategy (Section 3), while a change in state of a cell to air controlled is defined here by the combustion strategy that is presented in Section 6. • The cells' neighbourhoods. For each cell in the lattice, there is a neighbourhood that locally influences the evolution of the cell via the transition rule. The cells' neighbourhoods here are all the burning cells on the floor and the cells of the smoke layer showed that radiation to the cell is not obstructed by walls (eventual internal walls or core and, also, external walls in the compartments of non-convex shapes). Appl. Sci. 2021, 11, x FOR PEER REVIEW 4 of 31 The evaluation of energy and mass balance in the compartment is performed by the well-known 2 zone model OZone [9]. Therefore, quantities such as the time dependent thickness of the smoke layer Zs and its temperature Ts are obtained.

Ignition Strategy
Unlike what has been done in previous travelling fire models [16,18,20], the fire dynamics in GoZone is not imposed by assigning a predefined fire velocity or fire path; instead, it is regulated by the ignition strategy that is presented here.
Development and validation of this strategy was based on experimental data obtained on wood, first because wood-based products still represent a substantial fraction of the fire load in buildings, but also because of the wide availability of scientific data for this specific material. Conceptually speaking, it should be possible to adapt the model for other materials, assuming the availability of data on these materials. Application of the model in a compartment containing different materials would then require defining a combination rule, taking into account the proportion of each material present in the compartment. It will be shown in Section 5 that the proposed ignition strategy can be applied to replicate the global results in terms of ignition and propagation velocity as defined by a time constant linked to each type of occupation in EN1991-1-2,without any reference to the types of material that are present in the compartment.
Theories and models about the ignition of wood have evolved since the first studies started more than one century ago [29] from a concept of critical oven temperature [30] to the introduction of the time of exposure [31] and surface temperature [32][33][34], then to the idea of a constant incident heat flux [14,15,[35][36][37][38][39][40]. To study materials under constant incident heat flux exposure, a concept of critical heat flux " " was introduced: a constant flux level below which ignition cannot be obtained in practical conditions. However, this concept is based on a constant applied heat flux, not representative of a real fire. An ignition strategy considering a non-constant heat flux would thus probably be more appropriate to represent a real fire. One of the first works on time varying heat flux is from Bilbao et al. [41], who were still using a time decreasing heat flux, which is not representative of the first growing phase in a compartment fire. In 2006, Ji et al. [42] proposed a simple method to derive the time of wood ignition based on linearly increasing heat flux exposition.
The ignition strategy proposed here is based on the notion of critical energy (i.e., the minimum energy needed to be accumulated by a specific material to get ignited) and has been developed from a material characterisation in terms of ignitability, accounting for several kinds of time varying heat flux exposures. The concept is similar to the one used in the cumulative flux method [43] or in the flux time product (FTP) developed mainly for constant heat flux exposure [37]. The evaluation of energy and mass balance in the compartment is performed by the well-known 2 zone model OZone [9]. Therefore, quantities such as the time dependent thickness of the smoke layer Z s and its temperature T s are obtained.

Ignition Strategy
Unlike what has been done in previous travelling fire models [16,18,20], the fire dynamics in GoZone is not imposed by assigning a predefined fire velocity or fire path; instead, it is regulated by the ignition strategy that is presented here.
Development and validation of this strategy was based on experimental data obtained on wood, first because wood-based products still represent a substantial fraction of the fire load in buildings, but also because of the wide availability of scientific data for this specific material. Conceptually speaking, it should be possible to adapt the model for other materials, assuming the availability of data on these materials. Application of the model in a compartment containing different materials would then require defining a combination rule, taking into account the proportion of each material present in the compartment. It will be shown in Section 5 that the proposed ignition strategy can be applied to replicate the global results in terms of ignition and propagation velocity as defined by a time constant linked to each type of occupation in EN1991-1-2,without any reference to the types of material that are present in the compartment.
Theories and models about the ignition of wood have evolved since the first studies started more than one century ago [29] from a concept of critical oven temperature [30] to the introduction of the time of exposure [31] and surface temperature [32][33][34], then to the idea of a constant incident heat flux [14,15,[35][36][37][38][39][40]. To study materials under constant incident heat flux exposure, a concept of critical heat flux "q cr " was introduced: a constant flux level below which ignition cannot be obtained in practical conditions. However, this concept is based on a constant applied heat flux, not representative of a real fire. An ignition strategy considering a non-constant heat flux would thus probably be more appropriate to represent a real fire. One of the first works on time varying heat flux is from Bilbao et al. [41], who were still using a time decreasing heat flux, which is not representative of the first growing phase in a compartment fire. In 2006, Ji et al. [42] proposed a simple method to derive the time of wood ignition based on linearly increasing heat flux exposition.
The ignition strategy proposed here is based on the notion of critical energy (i.e., the minimum energy needed to be accumulated by a specific material to get ignited) and has been developed from a material characterisation in terms of ignitability, accounting for several kinds of time varying heat flux exposures. The concept is similar to the one used in the cumulative flux method [43] or in the flux time product (FTP) developed mainly for constant heat flux exposure [37].

Methodology
A series of 1D numerical analyses were made using SAFIR [7] with a wood semiinfinite media exposed to several laws of radiation exposure until the moment when the surface temperature reaches 350 • C, assumed as the ignition temperature [34,44]. The material was wood from EN 1995-1-2 [45] with a moisture content of 10%, a density of 450 kg/m 3 , a coefficient of convection of 25 W/m 2 K, and an emissivity of 0.8. For each simulation, a critical energy E cr is computed according to Equation (1).
where t ign is the time of ignition (surface temperature = 350 • C); t is the time; q is the applied heat flux; and q cr is the critical heat flux.

Results
In Figure 2, some of the flux exposure laws used are shown up to the computed ignition time. For a constant heat flux, exposure of 12.5 kW/m 2 ignition has not been observed, thus the critical heat flux was set at 12.5 kW/m 2 . This is in line with what has been observed by Spearpoint [46].

Methodology
A series of 1D numerical analyses were made using SAFIR [7] with a wood semiinfinite media exposed to several laws of radiation exposure until the moment when the surface temperature reaches 350 °C, assumed as the ignition temperature [34,44]. The material was wood from EN 1995-1-2 [45] with a moisture content of 10%, a density of 450 kg/m 3 , a coefficient of convection of 25 W/m 2 K, and an emissivity of 0.8. For each simulation, a critical energy Ecr is computed according to Equation (1).
where tign is the time of ignition (surface temperature = 350 °C); t is the time; q′ is the applied heat flux; and qcr is the critical heat flux.

Results
In Figure 2, some of the flux exposure laws used are shown up to the computed ignition time. For a constant heat flux, exposure of 12.5 kW/m 2 ignition has not been observed, thus the critical heat flux was set at 12.5 kW/m 2 . This is in line with what has been observed by Spearpoint [46].   Table 2 shows the same results for a linear variation in the form of q′(t) = m t.   Table 2 shows the same results for a linear variation in the form of q (t) = m t.  Values of the critical energy computed on the base of a critical heat flux of 12.5 kW/m 2 ranged from 350 to 600 kJ/m 2 . The value of 400 kJ/m 2 was implemented in the code GoZone, considering this value closer to critical energies obtained with a "quadratic" heat flux exposure Table 2, assumed to better represent a growing fire. Moreover, it can be seen that the time of ignition does not vary significantly if a value of 600 kJ/m 2 is used ( Figure 3).   Table 3 presents the results for a quadratic heat flux exposure law similar to the proposal of Eurocode 1 [47]: ′( ) = 1000 . Values of the critical energy computed on the base of a critical heat flux of 12.5 kW/m 2 ranged from 350 to 600 kJ/m 2 . The value of 400 kJ/m 2 was implemented in the code GoZone, considering this value closer to critical energies obtained with a "quadratic" heat flux exposure Table 2, assumed to better represent a growing fire. Moreover, it can be seen that the time of ignition does not vary significantly if a value of 600 kJ/m 2 is used ( Figure  3).

Comparison against Fire Test
A series of natural fire tests have been performed at Liege University in the framework of the European TRAFIR project. Incident heat flux history was measured at different locations from the centre of the fire through vertical and horizontal thermos-plate flux meters [44]. The proposed ignition strategy was tested by comparing the time of ignition evaluated using the recorded heat fluxes with the time of ignition observed during the test. Table 4 shows the results in terms of time of ignition observed in the test, time of ignition evaluated by the model using a critical energy of 400 kJ/m 2 , and the error between the model and the tests.

Comparison against Fire Test
A series of natural fire tests have been performed at Liege University in the framework of the European TRAFIR project. Incident heat flux history was measured at different locations from the centre of the fire through vertical and horizontal thermos-plate flux meters [44]. The proposed ignition strategy was tested by comparing the time of ignition evaluated using the recorded heat fluxes with the time of ignition observed during the test. Table 4 shows the results in terms of time of ignition observed in the test, time of ignition evaluated by the model using a critical energy of 400 kJ/m 2 , and the error between the model and the tests.  Figure 3 shows the value computed by the model for a critical energy of 400 kJ/m 2 and 600 kJ/m 2 as a function of the observed time of ignition, showing a low dependency of the result of the model on this parameter. It can be observed that the times of ignition were significantly longer than those mentioned in Tables 1-3. This is due to the fact that in these fire propagation tests, it took a significant amount of time for the fire to approach the points of measurements, which means that the received flux was below the critical heat flux until then.

Radiation toward the Cells
This section describes the models for computing the heat flux received by each cell, essentially by radiation [48] from the hot smoke or from the fire source, in order to apply the ignition strategy described in Section 3.

Radiation from the Smoke Layer
Radiation from the smoke layer was evaluated considering the smoke layer as an isothermal surface with temperature and elevation provided by the zone model (OZone) [49]. Absorption of radiation by air below the smoke layer is considered as negligible [50].
The flux received by a fuel cell j from a smoke cell i is evaluated accounting for the view factor between the two cells [47] and considering a smoke emissivity coefficient ε f equal to 1.
The contributions from all smoke cells that are visible by the fuel cell are summed up to compute the total flux received by the fuel cell.

Radiation from the Flames
The semi-empirical model known as the point source model (PSM) [51] was used here to compute the radiation from each burning cell to the fuel cells. PSM assumes the fire to radiate isotropically from a point "P" at the centre of the flame located at half of the flame length ( Figure 4).
The angle θ between the normal of the receiving surface and the line of sight (target to point source) allows for considering the orientation of the receiving target with respect to the point source. The incident radiant heat flux on the target surface is evaluated according to Equation (3), where RHR(t) is the fire power as a function of time and χ r is the radiative fraction coefficient set as 0.3 [52].  is the radiative fraction coefficient set as 0.3 [52].
Opinions on PSM are controversial in the literature. Casal stated that PSM overestimates the intensity of the thermal radiation for a target close to the fire [53]. Drysdale, in contrast, wrote that PSM underestimates thermal radiation for close targets [52]. Fleury, after having tested several semi-empirical models for the incident heat flux evaluation toward a wide range of target orientations and fire powers, concluded that PSM was the most accurate model among the ones analysed in his work [54]. PSM was tested here with respect to real experimental tests on localised pool fires and to CFD simulations. Figure 5 shows a comparison between the PSM and 18 different pool fire tests conducted within the framework of the RFCS project LOCAFI [13]. The flux measurements during the tests were taken at a fixed distance of approximately 3 m from the centre of the fire and each test had a different steady state fire power created by different diameters.  Opinions on PSM are controversial in the literature. Casal stated that PSM overestimates the intensity of the thermal radiation for a target close to the fire [53]. Drysdale, in contrast, wrote that PSM underestimates thermal radiation for close targets [52]. Fleury, after having tested several semi-empirical models for the incident heat flux evaluation toward a wide range of target orientations and fire powers, concluded that PSM was the most accurate model among the ones analysed in his work [54]. PSM was tested here with respect to real experimental tests on localised pool fires and to CFD simulations. Figure 5 shows a comparison between the PSM and 18 different pool fire tests conducted within the framework of the RFCS project LOCAFI [13]. The flux measurements during the tests were taken at a fixed distance of approximately 3 m from the centre of the fire and each test had a different steady state fire power created by different diameters. The angle θ between the normal of the receiving surface and the line of sight (target to point source) allows for considering the orientation of the receiving target with respect to the point source. The incident radiant heat flux on the target surface is evaluated according to Equation (3), where RHR(t) is the fire power as a function of time and is the radiative fraction coefficient set as 0.3 [52].
Opinions on PSM are controversial in the literature. Casal stated that PSM overestimates the intensity of the thermal radiation for a target close to the fire [53]. Drysdale, in contrast, wrote that PSM underestimates thermal radiation for close targets [52]. Fleury, after having tested several semi-empirical models for the incident heat flux evaluation toward a wide range of target orientations and fire powers, concluded that PSM was the most accurate model among the ones analysed in his work [54]. PSM was tested here with respect to real experimental tests on localised pool fires and to CFD simulations. Figure 5 shows a comparison between the PSM and 18 different pool fire tests conducted within the framework of the RFCS project LOCAFI [13]. The flux measurements during the tests were taken at a fixed distance of approximately 3 m from the centre of the fire and each test had a different steady state fire power created by different diameters.

Lf/2
Lf θ T P Figure 5. PSM compared with 18 pool fires; data from the LOCAFI project.
Gas burners of different constant powers (200, 500, 750, 1000 kW) have been modelled in FDS [15]. The incident radiative heat flux has been computed at 1 m, 2 m, and 3 m from the fire centre by varying the target orientation. Figure 6 shows the results for vertical orientation facing the fire source, represented with black dots, whereas the white filled dots are for the horizontal orientation. Gas burners of different constant powers (200, 500, 750, 1000 kW) have been modelled in FDS [15]. The incident radiative heat flux has been computed at 1 m, 2 m, and 3 m from the fire centre by varying the target orientation. Figure 6 shows the results for vertical orientation facing the fire source, represented with black dots, whereas the white filled dots are for the horizontal orientation. Figures 5 and 6 tend to confirm the conclusion drawn by Fleury [54] that PSM is an appropriate tool for the radiative heat flux estimation.
Fire is started by choosing a fire origin cell as most fires start as localised. In several analytical TFM, the fire is assumed to start as a line perpendicular to the path that the fire will follow, with the line as long as the width of the compartment. In this model, the fire is started by choosing a fire origin cell according to the fact that, except in the case of arson or explosion, most fires start as a local fire. The RHR of a burning cell is defined according to the indication provided by Eurocode 1-2 [47]. Assuming a fuel controlled fire, there are three main fire stages for the burning cells: (1) an initial growing phase, depending on the fire growth rate assigned; (2) a stationary regime in which the cell is burning entirely and the RHR is a function of the area of the fire and the RHRf; and (3) a linear descending phase starting when 70% of the fire load is consumed. In this context, if the fire propagates, the sum of RHR(t) produced by each cell is the total fire power and the total fire area is evaluated the same way by summing the fire area of each burning cell.
The PSM model cannot be used here with a representation of the whole fire engulfed area by a single point. This is due to the complex shape that the fire can take in a plan view depending on the fire origin and considering that the model should be applicable in compartments of any shape. In compartments with a concave plan view, the fact that a single point is visible for a cell or not would result in dichotomist yes/no possibilities, whereas a cell may see part of the fire source but not all of it. The PSM model has to be applied cell by cell. This yields the question of the evaluation of the flame length for each burning cell and of the mesh sensitivity of the discretised model.
Heskestad formulation (Equation (4)) is widely adopted to evaluate the flame length of localised fires [47,55]. Figures 5 and 6 tend to confirm the conclusion drawn by Fleury [54] that PSM is an appropriate tool for the radiative heat flux estimation.
Fire is started by choosing a fire origin cell as most fires start as localised. In several analytical TFM, the fire is assumed to start as a line perpendicular to the path that the fire will follow, with the line as long as the width of the compartment. In this model, the fire is started by choosing a fire origin cell according to the fact that, except in the case of arson or explosion, most fires start as a local fire. The RHR of a burning cell is defined according to the indication provided by Eurocode 1-2 [47]. Assuming a fuel controlled fire, there are three main fire stages for the burning cells: (1) an initial growing phase, depending on the fire growth rate assigned; (2) a stationary regime in which the cell is burning entirely and the RHR is a function of the area of the fire and the RHR f ; and (3) a linear descending phase starting when 70% of the fire load is consumed. In this context, if the fire propagates, the sum of RHR(t) produced by each cell is the total fire power and the total fire area is evaluated the same way by summing the fire area of each burning cell.
The PSM model cannot be used here with a representation of the whole fire engulfed area by a single point. This is due to the complex shape that the fire can take in a plan view depending on the fire origin and considering that the model should be applicable in compartments of any shape. In compartments with a concave plan view, the fact that a single point is visible for a cell or not would result in dichotomist yes/no possibilities, whereas a cell may see part of the fire source but not all of it. The PSM model has to be applied cell by cell. This yields the question of the evaluation of the flame length for each burning cell and of the mesh sensitivity of the discretised model.
Flame length in Equation (4) relates to the fire central axis for a given fire diameter D and power Q. It is not appropriate for evaluating the flame length of individual cells that form an otherwise greater fire because it would amount to considering a large number of fires with a small diameter located side by side. An approach inspired by a non-local variable concept is proposed here. The non-local variable concept is used, for example, as a solution for failure processes characterized by strain localisation into shear bands where related classical continuum theories reach the limits of application [56].
In general, all non-local variable methods rely on a characteristic length that governs the size of the region of influence on the point that is being considered. A characteristic length L CH of 2 m was found appropriate here through calibration of the model presented hereafter with circular shaped fires with a radius varying from 1 to 5 m and rate of heat release densities of 250 and 500 kW/m 2 .
An equivalent fire area was first evaluated for each burning cell i with central coordinates (x i , y i ) using Equation (5) where the surface integration was conducted on the whole burning domain.
An associated equivalent diameter was computed from the equivalent fire area (see Equation (6)).
where Q eq,i is computed by the product of the rate of heat release density by the equivalent area of the cell. Figure 7 shows the results of the non-local variable methodology in terms of the shape of the global fire as reconstructed from the length of the flame in each cell, for two different shapes of the fire area and two very different mesh refinements.
A comparison between the flame length evaluated with this methodology and the one provided by Heskestad can be conducted within the range of applicability of the Heskestad formulation, thus in the central axis of a circular fire area limited to a maximum radius of 5 m. Several circular fires of radius ranging from 2 m to 5 m were discretised. In Figure 8 solution for failure processes characterized by strain localisation into shear bands related classical continuum theories reach the limits of application [56].
In general, all non-local variable methods rely on a characteristic length that g the size of the region of influence on the point that is being considered. A charac length LCH of 2 m was found appropriate here through calibration of the model pre hereafter with circular shaped fires with a radius varying from 1 to 5 m and rate release densities of 250 and 500 kW/m 2 .
An equivalent fire area was first evaluated for each burning cell i with centra dinates (xi, yi) using Equation (5) where the surface integration was conducted whole burning domain.
An associated equivalent diameter was computed from the equivalent fire ar Equation (6)).
The flame length for the cell is then computed using Equation (7). where Qeq,i is computed by the product of the rate of heat release density by the equ area of the cell. Figure 7 shows the results of the non-local variable methodology in terms shape of the global fire as reconstructed from the length of the flame in each cell, f different shapes of the fire area and two very different mesh refinements.
(a)   A comparison between the flame length evaluated with this methodology and t one provided by Heskestad can be conducted within the range of applicability of t Heskestad formulation, thus in the central axis of a circular fire area limited to a maximu radius of 5 m. Several circular fires of radius ranging from 2 m to 5 m were discretised.  In Figure 9a the flame length at the centre of the circular area from the non-local va iable method was plotted as a function of the number of discretised elements; this show that the result was not mesh sensitive. In Figure 9b, a comparison between the analytic solution at the centre of the circular area and the non-local variable method showed error between 6 and 13%, which was also not mesh sensitive. In Figure 9a the flame length at the centre of the circular area from the non-local variable method was plotted as a function of the number of discretised elements; this showed that the result was not mesh sensitive. In Figure 9b, a comparison between the analytical solution at the centre of the circular area and the non-local variable method showed an error between 6 and 13%, which was also not mesh sensitive.
In Figure 10, the flame length at the centre was evaluated for different square fire areas, varying the cell size discretised with quadrilateral structured elements. The behaviour of the function was in line with that observed for a circular shape area in Figure 9a, appearing to be quasi mesh insensitive.

Fuel Control Situations
The rate of heat release (RHR) of each burning cell is pre-defined according to the recommendations of Eurocode 1 part 1-2 [47] for a fuel controlled fire, with three main stages: (1) an initial growing phase, depending on the assigned fire growth rate; (2) a stationary regime during which the RHR is equal to the product between the area of the cell and the value of given rate of heat release density; and (3) a linear descending phase starting when 70% of the fire load of the cell is consumed. When the fire propagates, the sum of RHR(t) produced by all cells makes the total fire power in the compartment used by the zone model. The advantage in defining the fire power evolution according to the recommendations of EN1991-1-2 is in having input parameters that are relatively easy to understand, familiar to practitioners (RHR density, fuel load, and fire growth rate). A 20 × 20 m 2 well ventilated compartment with the ceiling at 2.7 m was taken to check the growing phase evolution in terms of the rate of heat released evaluated by GoZone through a comparison with the analytical solution of EN1991-1-2. Table 5 shows the three different compartment typologies that were tested with regard to the fuel load density, rate of heat release density, and fire growth rate described in Tables E.4 and E.5 of EN1991-1-2. As has been conducted for the evaluation of the flame length, the temperature of the flame along the central axis of a circular fire was computed and is shown in Figure 11, first considering the fire as a whole (full line), and then considering the temperature evolution in the central cell of the discretised fire (dotted line).

Fuel Control Situations
The rate of heat release (RHR) of each burning cell is pre-defined according to the recommendations of Eurocode 1 part 1-2 [47] for a fuel controlled fire, with three main stages: (1) an initial growing phase, depending on the assigned fire growth rate; (2) a stationary regime during which the RHR is equal to the product between the area of the cell and the value of given rate of heat release density; and (3) a linear descending phase starting when 70% of the fire load of the cell is consumed. When the fire propagates, the sum of RHR(t) produced by all cells makes the total fire power in the compartment used by the zone model. The advantage in defining the fire power evolution according to the recommendations of EN1991-1-2 is in having input parameters that are relatively easy to understand, familiar to practitioners (RHR density, fuel load, and fire growth rate). A 20 × 20 m 2 well ventilated compartment with the ceiling at 2.7 m was taken to check the growing phase evolution in terms of the rate of heat released evaluated by GoZone through a comparison with the analytical solution of EN1991-1-2. Table 5 shows the three different compartment typologies that were tested with regard to the fuel load density, rate of heat release density, and fire growth rate described in Tables E.4 and E.5 of EN1991-1-2.

Fuel Control Situations
The rate of heat release (RHR) of each burning cell is pre-defined according to the recommendations of Eurocode 1 part 1-2 [47] for a fuel controlled fire, with three main stages: (1) an initial growing phase, depending on the assigned fire growth rate; (2) a stationary regime during which the RHR is equal to the product between the area of the cell and the value of given rate of heat release density; and (3) a linear descending phase starting when 70% of the fire load of the cell is consumed. When the fire propagates, the sum of RHR(t) produced by all cells makes the total fire power in the compartment used by the zone model. The advantage in defining the fire power evolution according to the recommendations of EN1991-1-2 is in having input parameters that are relatively easy to understand, familiar to practitioners (RHR density, fuel load, and fire growth rate). A 20 × 20 m 2 well ventilated compartment with the ceiling at 2.7 m was taken to check the growing phase evolution in terms of the rate of heat released evaluated by GoZone through a comparison with the analytical solution of EN1991-1-2. Table 5 shows the three different compartment typologies that were tested with regard to the fuel load density, rate of heat release density, and fire growth rate described in Tables E.4 and E.5 of EN1991-1-2. The fire origin was set in the centre of the compartment and the comparison was made up to the moment when the smoke layer was predominantly influencing the fire dynamics or when the fire completely engulfed the compartment. This is because the equation of Eurocode 1 does not stand when any of these two criteria is met. Figures 12 and 13 show the evolution of the power released during the growing phase of a fast and a medium fire. Outcomes provided by GoZone seem to be reasonably in line with the analytical solution.  The fire origin was set in the centre of the compartment and the comparison was made up to the moment when the smoke layer was predominantly influencing the fire dynamics or when the fire completely engulfed the compartment. This is because the equation of Eurocode 1 does not stand when any of these two criteria is met. Figures 12 and 13 show the evolution of the power released during the growing phase of a fast and a medium fire. Outcomes provided by GoZone seem to be reasonably in line with the analytical solution.   Figure 14 shows the fire propagation predicted by GoZone, whic of a medium fire from the time of first ignition to a moment just before fast and medium fire development were similar in terms of the shape o only one case is shown and is considered to be representative of the tw Figure 15 shows the evolution of the power for a slow fire scen general trend of the numerical model is still the one of the analytical s behaviour could be observed compared to Figures 13 and 14, where t did not increase monotonously, especially for cruder meshes.  Figure 14 shows the fire propagation predicted by GoZone, which is representative of a medium fire from the time of first ignition to a moment just before the flashover. Since fast and medium fire development were similar in terms of the shape of the fire evolution, only one case is shown and is considered to be representative of the two. Figure 15 shows the evolution of the power for a slow fire scenario. Although the general trend of the numerical model is still the one of the analytical solution, a different behaviour could be observed compared to Figures 13 and 14, where the power released did not increase monotonously, especially for cruder meshes.  This is due to the fact that, in this case, a slow fire produces some kind of r elling fire, as can be seen in Figure 16. In this figure, red cells are in the increasing state period of combustion; pink cells are in the descending phase; and blac completely burnt. Due to the discretised nature of the floor in cells of finite si the regular geometry, it may occur that a significant amount of cells enter in the ing branch or are completely consumed, whereas for a certain amount of time, n is ignited. This phenomenon is amplified by the larger size of the meshes and b ble symmetry of the problem, which makes groups of several cells, thus chan behaviour at the same time.  Figures 12, 13 and 15 showed that the simple adaptation rameter talpha in each cell led to an ignition strategy that globally reproduced th propagation velocity depending on the occupation of the compartment, bypassi cessity to define the type and proportion of the different materials present in the This is due to the fact that, in this case, a slow fire produces some kind of radial travelling fire, as can be seen in Figure 16. In this figure, red cells are in the increasing or steady state period of combustion; pink cells are in the descending phase; and black cells are completely burnt. Due to the discretised nature of the floor in cells of finite sizes and to the regular geometry, it may occur that a significant amount of cells enter in the descending branch or are completely consumed, whereas for a certain amount of time, no new cell is ignited. This phenomenon is amplified by the larger size of the meshes and by the double symmetry of the problem, which makes groups of several cells, thus changing their behaviour at the same time. This is due to the fact that, in this case, a slow fire produces some kind of radial travelling fire, as can be seen in Figure 16. In this figure, red cells are in the increasing or steady state period of combustion; pink cells are in the descending phase; and black cells are completely burnt. Due to the discretised nature of the floor in cells of finite sizes and to the regular geometry, it may occur that a significant amount of cells enter in the descending branch or are completely consumed, whereas for a certain amount of time, no new cell is ignited. This phenomenon is amplified by the larger size of the meshes and by the double symmetry of the problem, which makes groups of several cells, thus changing their behaviour at the same time. Comparison of Figures 12, 13 and 15 showed that the simple adaptation of the parameter talpha in each cell led to an ignition strategy that globally reproduced the desired propagation velocity depending on the occupation of the compartment, bypassing the necessity to define the type and proportion of the different materials present in the compartment, each with its specific ignition propensity.

Strategy for Air Control
The mass of oxygen in the compartment was considered as 23% of the total mass of air present in the compartment. The oxygen fraction in the air was equal to 21%, however, Comparison of Figures 12, 13 and 15 showed that the simple adaptation of the parameter t alpha in each cell led to an ignition strategy that globally reproduced the desired propagation velocity depending on the occupation of the compartment, bypassing the necessity to define the type and proportion of the different materials present in the compartment, each with its specific ignition propensity.

Strategy for Air Control
The mass of oxygen in the compartment was considered as 23% of the total mass of air present in the compartment. The oxygen fraction in the air was equal to 21%, however, OZone has been extensively validated with 23%, so this value was taken for sake of consistency [57]. The zone model considers the oxygen mass balance equation according to Equation (10), where a stoichiometric combustion factor of wood equal to 1.27 was considered for the combustion process [52].
The oxygen distribution in the compartment is assumed to be uniform according to the hypothesis of a well steered reactor and the fire is assumed to be in a fuel controlled regime as long as the oxygen content obtained by time integration of Equation (10) is positive. When the oxygen content drops to zero, the strategy in air-controlled regime is inspired from an extended fire duration concept, in the sense that the pyrolysis rate and, as a consequence, the power released in the compartment are proportional to the mass of oxygen mass entering from the openings in the compartment (Equation (11)).

RHR air
where H c,e is the effective heat of combustion taken as 14 MJ/kg. The total power available is distributed to the cells that are closest to the opening's location among those that were burning before the switch to the air-controlled regime. Ignition of new cells and propagation of the fire during the air controlled regime is still possible provided that the critical energy criteria is met in some cells and if the limit in terms of power given by Equation (11) is not exceeded. This will happen typically when either breakage of a new opening is detected (according to the different strategies implemented in the zone model [9]) or if some burning cells enter in the decreasing phase. In this case, the cell among the non-burning cells, which is the closest to an active opening, will ignite or re-ignite first.

Comparison with BST/FSR 1993: Test 2
The BST/FSR test series performed in 1993, better known as Gordon Cooke's fire tests, consist of nine different compartment fire tests performed in one compartment with varying fuel load, opening size, and ignition methodology [21]. The test compartment was 22.8 m long, 5.6 m wide, and 2.75 m high, with an opening on one short end of the compartment. The fuel was made of wood crib piles uniformly distributed on the whole floor area, ignited at the opposite side with respect to the opening position. The authors consistently observed in the nine tests performed that the fire travelled from the back of the compartment toward the opening, seeking air, then once all the fuel was consumed near the opening region, the fire travelled back again toward the ignition region.
This test series is emblematic because it may be regarded as an early example of what is now termed as "travelling fires", where the position of the burning zones moved and was severely influenced by the ventilation. Recently, Dai et al. [22] tried to model and represent the "BST/FSR Test 2" with FDS [15], obtaining a good agreement with the fire dynamic observed in the test. However, the rate of heat release assigned to each row of crib in the model was extrapolated based on the crib mass loss rate measured during the test and assigned crib by crib according to their position using several ramp functions. The authors thus showed how FDS is able to replicate travelling fire behaviour in under-ventilated condition but were not able to predict its complex behaviour a priori.
GoZone was compared with the BST/FSR test 2 in order to evaluate its ability to predict this particular fire dynamic. For test number 2, Kirby et al. provided the evolution of the temperature distribution in the compartment for three specific crib lines (Figure 17). These information were used here to extrapolate the fire dynamics and to compare the test results with the model's outputs. The compartment was modelled in GoZone with a mesh size of 1 m and discretised with 205 unstructured quadrilateral cells. Each cell was characterised with the parameters in Table 6. The fuel load density is given in Kirby et al. [21]; however, to characterise the fuel cell, GoZone also needs the rate of heat release density and fire growth rate. These input parameters were chosen according to the suggestion provided by Eurocode 1 part 2 where a fuel load density of 380 MJ/m 2 can be associated with a Hotel Room compartment according to Table E.4. (Annex E of EN1991-1-2) Therefore, the rate of heat release density and fire growth rate were chosen consistently according to Table E.5 (Annex E of EN1991-1-2). A comparison with the tests was thus made in a situation that is as close as possible to a design situation in which the fire action is not known a priori. Critical energy of ignition and critical heat flux were not tuned for this particular test but taken as determined in Section 3. Figure 18 shows the results in terms of fire evolution. The fire was ignited on the left hand side and initially propagated slowly toward the opening (Figure18a), then the increasing temperature in the hot gases and the associated radiation helped the fire to spread faster toward the opening. The compartment at 19 min was completely engulfed in flames. After a few minutes, a transition to an air-controlled situation due to a lack of oxygen took place. Cells on the opposite side of the opening, depicted in yellow ( Figure  18b), cannot burn due to a lack of oxygen. The air entering from the opening feeds the burning cells (red) and the fire travels from the opening to the opposite side (Figure 18c). At 64 min, the fire is completely extinguished.
Comparison between Figures 17 and 18 show that GoZone can predict the complex fire behaviour observed in this specific test in terms of fire position and total duration. The total simulation time on a laptop was 54 s. The compartment was modelled in GoZone with a mesh size of 1 m and discretised with 205 unstructured quadrilateral cells. Each cell was characterised with the parameters in Table 6. The fuel load density is given in Kirby et al. [21]; however, to characterise the fuel cell, GoZone also needs the rate of heat release density and fire growth rate. These input parameters were chosen according to the suggestion provided by Eurocode 1 part 2 where a fuel load density of 380 MJ/m 2 can be associated with a Hotel Room compartment according to Table E.4. (Annex E of EN1991-1-2) Therefore, the rate of heat release density and fire growth rate were chosen consistently according to Table E.5 (Annex E of EN1991-1-2). A comparison with the tests was thus made in a situation that is as close as possible to a design situation in which the fire action is not known a priori. Critical energy of ignition and critical heat flux were not tuned for this particular test but taken as determined in Section 3. Figure 18 shows the results in terms of fire evolution. The fire was ignited on the left hand side and initially propagated slowly toward the opening (Figure 18a), then the increasing temperature in the hot gases and the associated radiation helped the fire to spread faster toward the opening. The compartment at 19 min was completely engulfed in flames. After a few minutes, a transition to an air-controlled situation due to a lack of oxygen took place. Cells on the opposite side of the opening, depicted in yellow (Figure 18b), cannot burn due to a lack of oxygen. The air entering from the opening feeds the burning cells (red) and the fire travels from the opening to the opposite side (Figure 18c). At 64 min, the fire is completely extinguished.
Comparison between Figures 17 and 18 show that GoZone can predict the complex fire behaviour observed in this specific test in terms of fire position and total duration. The total simulation time on a laptop was 54 s.
Comparisons in terms of evolution of gas temperature among the compartment length just above crib lines 2, 6, and 10 in the centre line of the compartment are shown in Figure 19, as measured during the test, as evaluated by GoZone, and extrapolated from Dai et al. [22] as a result of FDS analysis. Figure 19a,b shows that the two separate peaks of temperature when the fire travels back and forth can be reproduced by the CFD replication of the test as well as by the GoZone a priori prediction, with the level of maximum temperature better predicted in the latter. Both models correctly predicted the maximum temperature in crib line 10, which was close to the opening (see Figure 19c). The cooling down phase was better predicted by GoZone. Comparisons in terms of evolution of gas temperature among the compartment length just above crib lines 2, 6, and 10 in the centre line of the compartment are shown in Figure 19Fig, as measured during the test, as evaluated by GoZone, and extrapolated from Dai et al. [22] as a result of FDS analysis.
(a) Comparisons in terms of evolution of gas temperature among the compartment length just above crib lines 2, 6, and 10 in the centre line of the compartment are shown in Figure 19Fig, as measured during the test, as evaluated by GoZone, and extrapolated from Dai et al. [22] as a result of FDS analysis.
(a)  Figure 19a,b shows that the two separate peaks of temperature when the fire travels back and forth can be reproduced by the CFD replication of the test as well as by the GoZone a priori prediction, with the level of maximum temperature better predicted in the latter. Both models correctly predicted the maximum temperature in crib line 10, which was close to the opening (see Figure19c). The cooling down phase was better predicted by GoZone.

Comparison with Veselì Travelling Fire Test
The Veselí fire test was aimed at investigating the fire dynamics and temperature distribution of a fuel controlled fire in a large open compartment in a two-storey building (Czech Republic, 2011) [58]. The compartment size was 10.4 m × 13.4 m on the floor and 4 m in height. The fire load was made of wood cribs with sticks of dimensions of 50 mm × 50 mm × 1000 mm with a moisture content of 12%. A linear channel filled with "petroleum" was used for the ignition of the fire. There was a 5 m wide and 2 m high opening in one of the long walls ( Figure 20). The gas temperature in the compartment was recorded during the test at different locations whereas the fire behaviour was recorded on camera. Information not directly measured during the test was evaluated numerically by Horova et al. by a CFD based simulation performed with FDS [58]. The test compartment as well as the GoZone model are as shown in Figure 20.

Comparison with Veselì Travelling Fire Test
The Veselí fire test was aimed at investigating the fire dynamics and temperature distribution of a fuel controlled fire in a large open compartment in a two-storey building (Czech Republic, 2011) [58]. The compartment size was 10.4 m × 13.4 m on the floor and 4 m in height. The fire load was made of wood cribs with sticks of dimensions of 50 mm × 50 mm × 1000 mm with a moisture content of 12%. A linear channel filled with "petroleum" was used for the ignition of the fire. There was a 5 m wide and 2 m high opening in one of the long walls ( Figure 20). The gas temperature in the compartment was recorded during the test at different locations whereas the fire behaviour was recorded on camera. Information not directly measured during the test was evaluated numerically by Horova et al. by a CFD based simulation performed with FDS [58]. The test compartment as well as the GoZone model are as shown in Figure 20. For the model in GoZone, the fire load area was discretised by a structured mesh with 24 cells of 1 m 2 (blue cells in Figure 20b) according to Section 5. The green cells do not contain any fuel and the red cell is the cell of origin of the fire. The comparison with GoZone is meant to be a priori, thus the values of the fuel load density, the rate of the heat release density, and the fire growth rate assigned to the fuel were taken according to the insights given by the authors of the test [49] (see Table 7). The critical energy and the critical heat flux were taken consistently with the outcomes of Section 3. The power released by the fire as a function of time was not directly measured during the test, however, it was evaluated a posteriori by means of the CFD based model with the software FDS [49]. Figure 21 shows the comparison between the results from FDS and GoZone. The peak value of the power as well as the peak time and the total duration were in the same order of magnitude with differences compared to a FDS of −8%, −10%, and −4%, respectively. A discrepancy can be noticed in the first phase, when FDS accounts for the petroleum burned as means of ignition, which was not the case in GoZone. For the model in GoZone, the fire load area was discretised by a structured mesh with 24 cells of 1 m 2 (blue cells in Figure 20b) according to Section 5. The green cells do not contain any fuel and the red cell is the cell of origin of the fire. The comparison with GoZone is meant to be a priori, thus the values of the fuel load density, the rate of the heat release density, and the fire growth rate assigned to the fuel were taken according to the insights given by the authors of the test [49] (see Table 7). The critical energy and the critical heat flux were taken consistently with the outcomes of Section 3. The power released by the fire as a function of time was not directly measured during the test, however, it was evaluated a posteriori by means of the CFD based model with the software FDS [49]. Figure 21 shows the comparison between the results from FDS and GoZone. The peak value of the power as well as the peak time and the total duration were in the same order of magnitude with differences compared to a FDS of −8%, −10%, and −4%, respectively. A discrepancy can be noticed in the first phase, when FDS accounts for the petroleum burned as means of ignition, which was not the case in GoZone. For the model in GoZone, the fire load area was discretised by a structured mesh with 24 cells of 1 m 2 (blue cells in Figure 20b) according to Section 5. The green cells do not contain any fuel and the red cell is the cell of origin of the fire. The comparison with GoZone is meant to be a priori, thus the values of the fuel load density, the rate of the heat release density, and the fire growth rate assigned to the fuel were taken according to the insights given by the authors of the test [49] (see Table 7). The critical energy and the critical heat flux were taken consistently with the outcomes of Section 3. The power released by the fire as a function of time was not directly measured during the test, however, it was evaluated a posteriori by means of the CFD based model with the software FDS [49]. Figure 21 shows the comparison between the results from FDS and GoZone. The peak value of the power as well as the peak time and the total duration were in the same order of magnitude with differences compared to a FDS of −8%, −10%, and −4%, respectively. A discrepancy can be noticed in the first phase, when FDS accounts for the petroleum burned as means of ignition, which was not the case in GoZone.  The fire evolution has been described by Horova et al. and Dai et al. [49,58,59]. The fire development was characterised by three main phases: phase I, from 0 to 15 min, when the fire leading edge spread along the fuel bed; phase II, from 20 to 25 min, when post-flashover conditions were observed; and phase III, from 25 to 40 min when the fire trailing edge progressed to the end of the fuel bed, characterised by a decay of the fire. In Figure 22, the fire evolution is shown through a comparison between pictures taken during the test (from Horova et al. [9]) and images from the GoZone Graphical User Interface (GUI). Pink cells represent a fire in its decay phase and black cells are painted where the fuel has been consumed completely.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 23 o The fire evolution has been described by Horova et al. and Dai et al. [49,58,59]. fire development was characterised by three main phases: phase I, from 0 to 15 min, w the fire leading edge spread along the fuel bed; phase II, from 20 to 25 min, when p flashover conditions were observed; and phase III, from 25 to 40 min when the fire trai edge progressed to the end of the fuel bed, characterised by a decay of the fire. In Fig  22, the fire evolution is shown through a comparison between pictures taken during test (from Horova et al. [9]) and images from the GoZone Graphical User Interface (G Pink cells represent a fire in its decay phase and black cells are painted where the fuel been consumed completely.  Moreover, it has been observed that the fire spread rate varied during the first ph Dai et al. [49] provided an estimation of the position of the leading edge of the fire u the moment when the fire covers the whole fuel bed surface (end of phase I). Figure 23 shows a comparison between the test and the model for the fire front p tion in the first phase of the fire; the average fire spread rate in the test, computed f the positions at 5, 10 and 15 min, nearly tripled in the third period of 5 min compare the initial period, from 3 to 8.9 mm/s. This figure shows that GoZone is capable of acco ing for a variation of fire spread rate depending on the changing fire conditions, in c trast to the travelling fire frameworks proposed in the literature so far, which are all ba on a constant spread rate.  Moreover, it has been observed that the fire spread rate varied during the first phase. Dai et al. [49] provided an estimation of the position of the leading edge of the fire up to the moment when the fire covers the whole fuel bed surface (end of phase I). Figure 23 shows a comparison between the test and the model for the fire front position in the first phase of the fire; the average fire spread rate in the test, computed from the positions at 5, 10 and 15 min, nearly tripled in the third period of 5 min compared to the initial period, from 3 to 8.9 mm/s. This figure shows that GoZone is capable of accounting for a variation of fire spread rate depending on the changing fire conditions, in contrast to the travelling fire frameworks proposed in the literature so far, which are all based on a constant spread rate.  Moreover, it has been observed that the fire spread rate varied during the first phase. Dai et al. [49] provided an estimation of the position of the leading edge of the fire up to the moment when the fire covers the whole fuel bed surface (end of phase I). Figure 23 shows a comparison between the test and the model for the fire front position in the first phase of the fire; the average fire spread rate in the test, computed from the positions at 5, 10 and 15 min, nearly tripled in the third period of 5 min compared to the initial period, from 3 to 8.9 mm/s. This figure shows that GoZone is capable of accounting for a variation of fire spread rate depending on the changing fire conditions, in contrast to the travelling fire frameworks proposed in the literature so far, which are all based on a constant spread rate.  Gas temperatures at different locations below a steel beam located above the fuel load were also measured. A comparison is given for the two thermocouples indicated in Figure 24. The gas temperature measured in the test, provided by the GoZone simulation and obtained by the CFD based model of Dai et al. [49] are compared in Figure 25.
The gas temperature measured in the test, provided by the GoZone simulation and obtained by the CFD based model of Dai et al. [49] are compared in Figure 25. Gas temperatures at different locations below a steel beam located above the fuel load were also measured. A comparison is given for the two thermocouples indicated in Figure 24. The gas temperature measured in the test, provided by the GoZone simulation and obtained by the CFD based model of Dai et al. [49] are compared in Figure 25. Overall, it is possible to observe a reasonable agreement in terms of the predicted gas temperature by both software (at a very different computational cost). However, both FDS and GoZone seem to underestimate the gas temperature in the compartment during the decay phase. The reason for this trend, which has been observed in different occasions, has yet to be identified.

Conclusions
A numerical model was presented to predict the fire development in compartments. While being based on simplifying assumptions, it proved to have the capability to predict several of the important features of the development of a fire in a compartment from ignition to flashover such as fire propagation and temperature development, in fuel control as well as air control regime. Having the capability to represent several features of the compartment (real shape of the plan view of the floor, real position and surface of any opening) and of the fire load characteristics (density, rate of heat release density, time constant for fire propagation, elevation above ground level, all possibly varying across the compartment, plus any position possible for the fire origin), whereas the run time of current problems on a laptop is in the order of a minute, this model has the ambition to occupy an empty niche between the geometry constrained analytical models presented thus far for travelling fires and the more sophisticated CFD analyses.
The basic idea is to combine an existing two zone model with a discretisation of the floor by a mesh of triangular or quadrilateral cells of finite size, with each cell possibly having its own characteristics with regard to fire load and its own status with regard to combustion. Fire propagation is driven by radiation from burning cells as well as from the hot smoke layer.
Different sub models required particular attention for the development of this tool. The point source model was adopted to represent the fire source and to compute the radiation that it will generate to the zones of the compartment where combustion has not yet started. Discretisation of the floor and the associated fire load into cells of finite size led to the development of a non-local model concept for computing the flame length in the burning cells. It can be applied to a burning zone of any shape, which is an essential requirement for a model that has the ambition to not be geometrically constrained. This model was tested against an experimental campaign and CFD analyses, which showed good agreement compared with the Heskestad equation in terms of flame length at the centre line of the fire (for shapes of the fire where this equation applies) and it appeared to be almost mesh insensitive for a range of practical mesh sizes. Overall, it is possible to observe a reasonable agreement in terms of the predicted gas temperature by both software (at a very different computational cost). However, both FDS and GoZone seem to underestimate the gas temperature in the compartment during the decay phase. The reason for this trend, which has been observed in different occasions, has yet to be identified.

Conclusions
A numerical model was presented to predict the fire development in compartments. While being based on simplifying assumptions, it proved to have the capability to predict several of the important features of the development of a fire in a compartment from ignition to flashover such as fire propagation and temperature development, in fuel control as well as air control regime. Having the capability to represent several features of the compartment (real shape of the plan view of the floor, real position and surface of any opening) and of the fire load characteristics (density, rate of heat release density, time constant for fire propagation, elevation above ground level, all possibly varying across the compartment, plus any position possible for the fire origin), whereas the run time of current problems on a laptop is in the order of a minute, this model has the ambition to occupy an empty niche between the geometry constrained analytical models presented thus far for travelling fires and the more sophisticated CFD analyses.
The basic idea is to combine an existing two zone model with a discretisation of the floor by a mesh of triangular or quadrilateral cells of finite size, with each cell possibly having its own characteristics with regard to fire load and its own status with regard to combustion. Fire propagation is driven by radiation from burning cells as well as from the hot smoke layer.
Different sub models required particular attention for the development of this tool. The point source model was adopted to represent the fire source and to compute the radiation that it will generate to the zones of the compartment where combustion has not yet started. Discretisation of the floor and the associated fire load into cells of finite size led to the development of a non-local model concept for computing the flame length in the burning cells. It can be applied to a burning zone of any shape, which is an essential requirement for a model that has the ambition to not be geometrically constrained. This model was tested against an experimental campaign and CFD analyses, which showed good agreement compared with the Heskestad equation in terms of flame length at the centre line of the fire (for shapes of the fire where this equation applies) and it appeared to be almost mesh insensitive for a range of practical mesh sizes.
Radiation from the hot smoke to the fire load was obtained by the discretisation of the lower plane of the hot smoke that was mirrored by the discretisation of the floor and computation of the view factors from all cells of the hot smoke to all cells on the floor (with due consideration of possible obstruction in the case of the ground floor with a concave shape or in the case of obstacles present in the compartment such as internal walls).
An ignition strategy was developed based on the concepts of critical energy of ignition and critical heat flux. This strategy has shown reasonably good prediction capabilities including for transient heat flux exposures.
Air control regime is triggered by the total oxygen content in the compartment, and the combustion or not of the different cells is then determined by their proximity to individual openings that deliver air in the compartment.
The model was implemented in the existing OZone model and thus benefits from all capabilities of this model such as the possibility of openings that appear at predefined times or according to defined temperature criteria, the presence of vertical as well as horizontal openings eventually with forced ventilation, and heat loss through the boundaries of the compartment were explicitly considered.
It also had the same inherent limitations such as a horizontal floor and ceiling and a ratio between the characteristic horizontal dimension of the compartment and a floor to ceiling distance that was not excessive.
The new code in which these developments were implemented was given the name of GoZone. The next step will be the automatic transmission of information regarding the results of this fire model to a thermo-mechanical software such as SAFIR ® , in order to allow for the simulation of the behaviour of a load bearing structure located in the fire compartment.