Comparative Study on Occupant Evacuation with Building EXODUS and a Cellular Automaton Model

Building EXODUS software is used to calculate the evacuation times and simulate the evacuation behavior. The results and laws are compared with those from a 2D Cellular Automaton (CA) random evacuation model developed by our group. EXODUS simulation is more reasonable than the CA simulation in the case of evacuation from a simple room, but CA model is more reasonable in the case of evacuation in a long corridor after bottlenecks. As far as the evacuation from a simple room with a single exit is concerned, there is a critical value of exit width. The value of exit width should be bigger than the critical value in order to ensure a dilute pedestrian flow, but the value doesn’t need to be too big. The bigger the original occupant density, the longer the evacuation time is. They can be fitted as a linear relationship. The principle of taking the shortest route is not always useful. If the distribution of occupant density is not uniform at each building part, balancing the use efficiency of each exit should be the main principle in order to improve evacuation efficiency. All the above laws can be obtained both from EXODUS and the CA model.


Introduction
Available Safe Egress Time (ASET) and Required Safe Egress Time (RSET) are two important elements in evacuation evaluation.ASET is determined by the evolvement of disasters, which is not discussed within this paper.RSET is composed of detection and alarm time, pre-movement time [1] [2] [3] and occupant movement time.Pre-movement time is related to the type of building and occu-pant, the mode of fire, the state of occupants, and etc. Fahy

et al. deem that
Weibull distribution is very suitable to describe pre-movement time [4], which ranges from several seconds to a few 10 minutes.This study only focuses on occupant movement time i.e. detection and alarm time, pre-movement time are not considered.
At present, many researchers are committed to develop empirical formulas and models on evacuation movement time.Familiar empirical formulas are Togawa formulas applied on engineering calculation of densely occupied public site [5], Melinek and Booth formulas applied on the calculation of the shortest evacuation time of high-rise buildings [6], Pauls formulas applied on the calculation of the evacuation time of staircases [7], and the method put forward by Nelson et al. (NFPA) based on the research work of Pauls, Proulx, Predtechenskii and Milinskii [8].Familiar models involve continuous models and discrete models.
Jamming is easily happened in evacuation crowd with high density, especially at narrow and long passages, doorways, stairs, etc.The pedestrian flow can change from a lamina flow to a stop-and-go wave, and then to a turbulent flow with the increase of density.Panicking individuals in evacuation tend to show maladaptive behavior like jamming and life-threatening overcrowding.The physical interactions in the jammed crowd add up and cause dangerous pressures up to 4450 N/m which can bend steel barriers or push down brick walls.Escape is slowed by fallen or injured people acting as obstacles.In addition, backtracking and counterflows, overtaking behavior, sub-group phenomenon etc. are also studied by many researchers.
Simulation of occupant behavior during evacuation is very complex and not all the behaviors can be quantified.Apparently models including more rules have become a trend.In this paper, the software of building EXODUS developed by the University of Greenwich is applied to analyze the influence of occupant density, exit width and exits distribution on evacuation.The results are compared with those of a 2D cellular automaton model developed by ourselves.Only simple rooms without obstacles are considered in this study.

Model Description
The common ground on EXODUS and our CA model can be summarized as below: 1) They are discrete models taking into consideration people-people, people-fire and people-structure interactions; 2) They have been written in C++ using Object Orientated techniques and rule-base concepts to control the simulation; 3) Building structure is represented by a two-dimensional grid.Each node represents 0.5 m × 0.5 m, which can hold only one person.

Building EXODUS Model
The rules of building EXODUS (V6.1) [24] can be categorized into five interacting sub-models, the OCCUPANT, MOVEMENT, BEHAVIOUR, TOXICITY and HAZARD sub-models.The last two sub-models are not used in this study.
We only introduce the rules and attributes concerning route selection.The database system of building EXODUS can provide reliable parameters about occupant attribute, building geometry, evacuation environment and so on.For more details, please refer to the user guide and technical manual of EXODUS.
Geometries are represented as two-dimensional grids.Each location on a grid is called a node, and each node may be linked to its nearest neighbors by four or eight arcs as shown in Figure 1.Occupants can move along those arcs from node to node.Walking speed and direction can be adjusted according to the change of environment, such as the crowd density and the distribution of obstacles.Occupants always move from the node of high danger grade to the node of low danger grade and eventually toward the exits.
The OCCUPANT sub-model defines each individual as a collection of attributes which broadly fall into four categories, physical, psychological, experiential and hazard effects.The Target Door attribute allows an occupant to be directed to a specific external door, regardless of the potential map.The values this attribute may take are the names of the exits in the geometry.The default value is "Nearest Door", i.e. follow the potential map, which is applied in this paper.
The MOVEMENT sub-model controls the physical movement of individual occupants from their current position to the most suitable neighboring location, or supervises the waiting period if one does not exist.The movement may involve such behavior as overtaking, sidestepping, or other evasive actions controlled by the BEHAVIOR sub-model.

CA Model
Von Neumann neighborhood is applied to consider four possible movement directions, i.e. up and down, left and right.In each time step, there are two basic problems to be resolved, i.e. route selection and conflict resolution.Each occupant checks all the cells within visual field and selects a proper target cell according to the value of ij p determined by the following formula: ( ) ( ) ( ) ( ) ( )( ) where ij p denotes the probability of moving to the cell ( ) that the cell ( ) , i j is occupied by an occupant at the time of t , otherwise, ij n = 0; the parameter N is introduced: where ij S and S k denote position attraction and its influence coefficient of cell ( ) , i j , respectively; ij R and R k denote repulsive force of occupants around and its influence coefficient of cell ( ) , i j , respectively; ij A and A k denote at- traction of occupants around and its influence coefficient of cell ( ) , i j , respec- tively; ij D and D k denote attraction of movement direction and its influence coefficient of cell ( ) , i j , respectively; ij F and F k denote fire repulsive force and its influence coefficient of cell ( ) going with the crowd or the kin behavior.The adjacent occupants form ij R in order to avoid colliding with each other.ij S is determined by the choice of safety exits.It reflects the geometry of the structure and the locations of the exits.In general, the nearer to the exit, the greater the value of the cell is.The simple formula is: , max min where ( ) ie je denotes the coordinates of each exit.The second item on the right of the equation denotes the minimum of the distance between ( ) , i j and each exit.The first item is the maximum of the above values.
A set of intelligent local rules are introduced into the model: 1) Determine the probability of each adjacent cell based on Equations ( 1)-(3).2) Each occupant chooses one of the adjacent cells as its object at the next time step depending on the probability of each cell in its visual field.3) For each cell where more than one occupant want to enter, randomly assign it to one of them with a certain probability; other occupants still stay where they are.In this paper, we consider that every competitor has the same probability.4) In order to avoid a determinate model, a random slowdown rule is introduced: Give each occupant who has decided to move a probability to stay.An experiential value, 5%, is used here because it's most close to real data.5) After parallel updating each occupant in each time-step, compute each cell's probability if necessary.6) Do the above until all the occupants finish their evacuation.For more details, please refer to our previous work [25] 1.For the accuracy of the data, all the results in this paper are the mean values of 10 times calculation in the same condition.

Intelligence
The different intelligence and pedestrian flow configuration of building EXO-DUS and our CA model are shown by three examples.

Case One
Occupant distribution states of EXODUS and the CA model during evacuation process are very different as shown in Figure 2. In our CA model, all occupants form an assembling state around the exit very quickly, while only the occupants around the exit form an assembling state and the occupants away from the exit almost keep still at the initial stage in EXODUS.EXODUS is able to represent the delay that evacuees might experience before commencing their evacuation.This is attributed to each agent who then delays their response accordingly.By default, this ranges from 0 -30 seconds.During the middle stage, occupants form a semicircle assembling state at the exit in EXODUS, and there are scarcely any interspaces.While occupants form a prolate assembling state and there are more interspaces in the CA model.EXODUS can more reasonably simulate the characteristic of assembling in the case of evacuation from a simply room because there are eight and four possible movement directions for each occupant in EXODUS and the CA model, respectively.

Case Two
Occupant evacuations from a room to a long corridor are simulated.The whole structure is a 20 × 100 area with a room of 20 × 40 and a total number of 371 occupants.Single door (d = 4) and two doors (each d = 2) are considered.As shown in Figure 3, compared with the CA simulation, only the middle one of three parts of the corridor is occupied after the bottleneck (door).In Figure 3(b), the occupants form two separate pedestrian flows in the corridor in EXODUS simulation; while the occupants are almost evenly distributed at the whole width of the corridor and form a single pedestrian flow after the bottlenecks in the CA simulation.Our CA model is more reasonable in this case.

Case Three
Evacuation from a room with two exits and initially 130 occupants is simulated.s respectively.In the CA model, the bigger the radius of visual field, the more reasonable route choice can be made and thus a shorter evacuation time is needed.

Influence of Exit Width
The physical distance between nodes is 0.5 m, which is assumed as one unit.
Each dimensionless size can be obtained by actual size divided by 0.5 m.First of all, the parameters d, q, ρ, t are interpreted as below: d is dimensionless exit width, q is flux per dimensionless unit exit width (people/s), ρ is occupant density (the number of occupants divided by the number of nodes of the structure) and t is evacuation time (s).separated by two rows of unfilled nodes to prevent the premature occurrence of jamming [25].Figure 6 shows the results of EXODUS and the CA model.
We can see that q decreases with the increase of d.By observing the dynamic play of EXODUS, we find that occupants form a stable and regular dense pedestrian flow state when d is smaller than 4, while the pedestrian flow takes on a dilute state when d is bigger than 4.
Evacuation time decreases nonlinearly with the increase of d. viz.when d is small, increasing d is very helpful for decreasing the evacuation time.When d is bigger than 4, the evacuation time decreases very slowly; especially when d is bigger than 8, evacuation time is almost unchanged, and the use ratio of the door (q) is very low.In this case, d = 4 is regarded as the critical value.

Case Five
The room is sized 50 × 50 with 1365 occupants initially randomly distributed.
The results are shown as Figure 7. Although, the size of the room and the density of the occupants are different with those of case four, the laws we obtained are the same.In this case, d = 10 is regarded as the critical point.
As far as the evacuation form a simple room with a single exit is concerned, there is a critical value of exit width.In practice, the value of exit width should be bigger than this critical value in order to ensure a dilute pedestrian flow and a safe evacuation.Occupants need not to queue up at the exits when d is bigger than the critical value.Considering the other factors such as the economy, the appearance, the structure safety and etc., it's no need to make the exits too big.
We can obtain the same rules form EXODUS and the CA model, but the evacuation times of the CA model are slightly greater than those of EXODUS.

Influence of Initial Occupant Density
A room sized 18 × 14 with a single exit (d = 2) is calculated.As shown in Figure 8(a), when occupant density ρ is smaller than 0.4, q increases with the increase of ρ; While q changes a little (the flux of exit takes on a saturated state) when ρ is greater than 0.4.We can get the same law from EXODUS and the CA model, but the values of q in the CA simulation are bigger than those in EXODUS.
As shown in Figure 8(b), the bigger the occupant density, the longer the evacuation time is.The results can be linearly fitted as: t = 205.263ρ+ 3.635 (EXODUS).The results of CA are slightly smaller than those of EXODUS.In practice, occupant density should be controlled to be less than 0.4, i.e.Corresponding to the four cases in Figure 9, the evacuation times simulated by EXODUS are 61.5 s, 75.6 s, 63.8 s, and 79.1 s, respectively.By the CA simulation, they are 55.8 s, 62.9 s, 53.1 s, and 67.0 s, respectively.Compared with EXODUS, the evacuation times by CA simulation are reduced by 9.3%, 16.8%, 16.8%, and 15.3%, respectively.In the CA simulations, the principle of taking the shortest routes is not always been applied, and occupants choose the evacuation routes depending on not only the distance from the exits, but also the distribution of the occupants around in their visual field.In Figure 9(b), some occupants originally located at the right part choose the left exit as their target.In these cases, the principle of taking the shortest route is not always useful when the distribution of occupants is not uniform at each area, while balancing the use ratio of each exit is benefit for improving evacuation efficiency.In practice, the distribution and size of each exit should be adjusted according to the actual distribution of occupant density.More exits and greater exit width should be designed at the areas with higher occupant density.

Summary
In this paper, the software of building EXODUS is used to analyze the influence of some physical factors on evacuation including exit width, layout of exits, occupant density and distribution, and etc.The results and laws by EXODUS are compared with those by our CA simulations.The following conclusions can be obtained: 1) The evacuation pedestrian flow states in EXODUS simulation are very different from those in the CA simulation.In the beginning, some occupants away from the exits almost keep still in EXODUS simulation, while all occupants try to move toward the exits in the CA simulation.As far as the evacuation from a simple room is concerned, EXODUS is more reasonable than the CA model to exhibit the assembling state at the exit.It takes on a semicircular shape in EXODUS and a prolate shape in the CA model, respectively.As far as the evacuation in a long corridor after bottlenecks is concerned, the CA model is more reasonable than EXODUS.The occupants evenly distribute at the whole width of the corridor in the CA model, while the widths of pedestrian flow are limited to the width of the bottlenecks in EXODUS.
2) There is a critical value of exit width.If the value is smaller than the critical value, the pedestrian flow takes on a dense state and the evacuation time is very long.If the value is too big, evacuation time can only be reduced a little, and the use ratio of exit is reduced a lot.In practice, the value of exit width should be bigger than the critical value in order to ensure a dilute pedestrian flow and a safe evacuation.For example, the critical value is 4 for a room sized 18 × 14, and it is 10 for a room sized 50 × 50.The same laws can be obtained from both EXODUS and the CA model, but the evacuation times of the CA model are slightly greater than those of EXODUS.
3) The bigger the original occupant density, the longer the evacuation time is.
They can be fitted as a linear relationship.This law can be obtained from both EXODUS and the CA simulation, but the evacuation times of the CA model are slightly smaller than those of EXODUS.

Figure 1 .
Figure 1.Connection mode of cells in EXODUS.

3. 1
. Randomicity A case of a 30 × 30 (i.e. 15 m × 15 m) room with a single exit located at the mid-dle of the wall are simulated by both EXODUS and the CA model.Exit width d equals 3, i.e. 1.5 m, and 600 occupants randomly located in the beginning.The results are shown in Tbl.1.In the CA model, each time step (ts) represents different real time based on walking velocity and the size of a cell.By experience, 1ts represents 0.66 s (walking velocity is 0.76 m/s, the length of a cell is 0.5 m, 0.5/0.76≈ 0.66 s).The mean values of evacuation time by EXODUS and the CA model are 325.67 s and 325.45 s (i.e.493.1 × 0.66 = 325.45s), respectively.The relative deviations of the CA model are slightly bigger than those of EXODUS as shown in Table

(a) 2 Figure 2 .
Figure 2. Occupant distribution states at different times with initially 300 people in a square room.

Figure 3 .
Figure 3. Occupant distribution states in a long corridor after the bottlenecks.

Figure 4 .Figure 5 .
Figure 4. Occupant distribution states at different times of case three.3.3.1.Case FourA room (18 × 14) that initially contains 200 randomly distributed occupants is calculated.As shown in Figure5, the occupants and the exit door are purposely

Figure 6 .Figure 7 .
Figure 6.Results of case four, (a) Relations between q and d; (b) Relations between t and d.

AFigure 8 .
Figure 8. Results of evacuation with different density, (a) Relations between q and ρ; (b) Relations between t and ρ.

Figure 9 .
Figure 9. Original states with different layout of exits and density distribution, (a) opposite layout，the density of the upper half part is twice than that of the lower half part; (b) opposite layout，the density of the right half part is twice than that of the left half part; (c) adjacent layout, the density of the right half part is twice than that of the left half part; (d) adjacent layout, the density of the left half part is twice than that of the right half part.

4 )
Both EXODUS and our CA model have certain intelligence.EXODUS simulation corresponds to the CA simulation with a visual field radius of 2. The principle of following the shortest route (potential map) is applied both in EXODUS and the CA model.The CA model also considers the distribution of occupants and the visual field.If the distribution of occupant density is not uniform in each building part, balancing the use ratio of each exit should be the main principle.40900), the National Natural Science Foundation of China (Grant No. 51006072) and the develop project of master pilot of management science and Engineering (1021ZK14000100502, 4521ZK130059047).The authors deeply appreciate the supports.
exit (i.e. the total evacuation time) and the left exit are 75.6 s and 39 s, respectively.The former is almost twice than the latter.In the CA simulation, the evacuation times of the right exit and the left exit are 62.9 s and 42.2 s, respectively. right