[go: up one dir, main page]

Academia.eduAcademia.edu
On the effectiveness of multi-feature evacuation systems: an agent-based exploratory simulation study Kashif Zia1 , Umar Farooq2 , Muhammad Shafi1 and Alois Ferscha3 1 Faculty of Computing and Information Technology, Sohar University, Sohar, Oman Department of Computer Science, University of Science and Technology Bannu, Bannu, Khyber Pakhtunkhwa, Pakistan 3 Institute of Pervasive Computing, Johannes Kepler University, Linz, Austria 2 ABSTRACT Evacuation modeling and simulation are usually used to explore different possibilities for evacuation, however, it is a real challenge to integrate different categories of characteristics in unified modeling space. In this paper, we propose an agent-based model of an evacuating crowd so that a comparative analysis of a different sets of parameters categorized as individual, social and technological aspects, is made possible. In particular, we focus on the question of rationality vs. emotionalism of individuals in a localized social context. In addition to that, we propose and model the concept of extended social influence, thereby embedding technological influence within the social influence, and analyze its impact on the efficiency of evacuation. NetLogo is used for simulating different variations in environments, evacuation strategies, and agents demographics. Simulation results revealed that there is no substantial advantage of informational overload on people, as this might work only in those situations, where there are fewer chances of herding. In more serious situations, people should be left alone to decide. They, however, could be trained in drills, to avoid panicking in such situations and concentrate on making their decisions solely based on the dynamics of their surroundings. It was also learned that distant connectivity has no apparent advantage and can be ruled out while designing an evacuation strategy based on these recommendations. Submitted 19 February 2021 Accepted 15 April 2021 Published 19 May 2021 Subjects Adaptive and Self-Organizing Systems, Agents and Multi-Agent Systems, Algorithms Corresponding author Umar Farooq, umar@ustb.edu.pk Keywords Crowd evacuation, Agent based modeling, Rationality vs. emotionalism, Social and Analysis of Algorithms, Scientific Computing and Simulation Academic editor José Manuel Galán influence, NetLogo Additional Information and Declarations can be found on page 21 INTRODUCTION DOI 10.7717/peerj-cs.531 Copyright 2021 Zia et al. Distributed under Creative Commons CC-BY 4.0 OPEN ACCESS Emergency evacuation is often associated with crowd panic, which sometimes turns into a catastrophe. Usually, it happens due to irrational behavior based on factors such as fear, pain, suffocation, or hassle of just a few emotionally intense individuals in the crowd. Disturbance initiated by a small group of individuals transforms into herding, which might cause injuries and deaths (Zia et al., 2016). It is often hard to find empirical evidence for such incidents in terms of factors resulting in them. An experimental examination is generally considered dangerous and unethical. Therefore, it is a normal practice to utilize How to cite this article Zia K, Farooq U, Shafi M, Ferscha A. 2021. On the effectiveness of multi-feature evacuation systems: an agentbased exploratory simulation study. PeerJ Comput. Sci. 7:e531 http://doi.org/10.7717/peerj-cs.531 simulation for modeling and investigation of such systems. Many modeling methodologies have been adopted for conducting evacuation studies in the literature (Pelechano & Malkawi, 2008). One way of modeling such systems focuses on considering global flows (Cova & Johnson, 2003) while the other concentrates only on local interactions (Yuan & Tan, 2007). The majority of practitioners advocate the use of latter (individual-based models) over the former (flow-based models), as it is believed that a global outcome is what emerges out of local interactions, which could be unexpected most of the time and must be dealt with accordingly. However, the understanding of the emergence of global patterns in relation to rules of local decision-making of individuals poses more questions than answers. Furthermore, other factors related to population, environment, and information must also be considered due to their obvious relevance. Ideally, an evacuation strategy should utilize all exits, in such as a way that the distribution of population across the exits is well balanced in terms of time, which makes it a perfect flow control application (Helbing, Farkas & Vicsek, 2000; Helbing, Johansson & Al-Abideen, 2007). Evacuation models have been proposed focusing only on globally scaled collectives (Von Schantz & Ehtamo, 2015) such as clogging (Seyfried et al., 2009), bottlenecks (Kretz, Grünebohm & Schreckenberg, 2006) and faster-is-slower (Garcimartín et al., 2014) effects. However, modeling of an evacuation situation solely based on physics has been greatly criticised (Reilly, Dillon & Guikema, 2018; Batty, 2001; Radianti et al., 2013). The reason, as stated above, is that individual and localized (social) decisions happening at a local scale manifest into global patterns (Kania & Kramer, 2013). Therefore, at an individual level, an agent must be considered as a social entity (Schadschneider et al., 2011), instead of a particle, having its own preferences, biases and emotions, and capability to adapt according to the changing dynamics of the environment. Therefore, more and more researchers have started using agent-based models for evacuation simulation (Wagner & Agrawal, 2014; Almeida, Rosseti & Coelho, 2013). Since people are the central entities in these situations, it is, therefore, important to represent them reasonably well. However, it is hard to get the evidence of human behavior, particularly, in life-critical situations like evacuation. Therefore, in most of the cases, the questions about the efficiency of egress from a troubled environment are answered using behavioral theories and patterns observed in other domains such as crowd mobility, navigation, and finding ways. In addition to that, the dimensions across which such an inquiry can be extended are numerous including physical modalities (such as smoke and visibility, building design, signage, and exits’ designs), crowd dynamics (such as crowd density, lanes and flows, and population diversity), information (such as vision and voice, ICT enabled infrastructure, and the availability of information through social network or local connectivity), individual decision-making (such as emotions, rationalism, and panic behavior), and collective effects (such as leader-following, herding effect, and peer effect). This huge set of aspects could be classified under individual, social and informational domains. Since, most research in evacuation simulation is done at a limited scope, a combinatorial model meaningfully combining these aspects is a challenge and focus of this paper. Zia et al. (2021), PeerJ Comput. Sci., DOI 10.7717/peerj-cs.531 2/25 Zia & Ferscha (2020) addressed this challenge by developing a common agent based model combining individual, social and technological aspects. The model was tested on a grid based hypothetical environment with an aim of getting an insight into the effectiveness of a multi-feature evacuation strategy. Their model, however, has an apparent inconsistency as it transmits the status information of exits to agents using spatial features of the environment instead of social spread. The model proposed in this work helps coping this issue by replacing spatial spread of information based on a mechanism of pure social spread, thus, exploiting a social network of people. The proposed model also explores the impact of long-distance relationships, thereby introducing a new set of exit strategies, tested in an environment of two exits and having different geometric variations (representing complexity). The selection of modeling methods certainly depends on modeler choice and modeling requirements. Agent Based Modeling ABM, (Macal, 2016; Wilensky & Rand, 2015) is one of the most promising techniques for modeling and simulation, which is capable of handling all the multi-facet and multi-layer complexity of crowd panic and evacuation systems. Salient features of (ABM) such as autonomy, re-activity, pro-activity and social interaction make this method a natural choice for scenarios requiring autonomous and adaptive participating agents (Liu, Jiang & Shi, 2016) under the influence of individual, social and informational aspects. With ABM, we started developing models of evacuation simulations based on purely physical acpects (Ferscha & Zia, 2009), which were, later on, extended to incorporate cognitive aspects of human thinking (Zia et al., 2011). Social aspects and neighborhood influence were, then, integrated into our models (Sharpanskykh & Zia, 2014). Furthermore, we developed decision-making models based on game-theoretic approach (Zia et al., 2016). In our most recent work, we investigated the generation and spread of information from exits (Zia & Ferscha, 2020). This paper contributes by proposing a model extension that: 1. removes the inconsistencies of the model proposed in Zia & Ferscha (2020), where the spatial nature of information spread is replaced by a pure social spread. 2. introduces a new behavior in agents’ decision-making repertoire, where the agents are more likely influenced by potentially long-distance ties, when their mobility is hindered. The incorporation of a diversified but related set of features, in a combinatorial way, is itself an achievement, however, the novelty of this work comes from the inclusion of network structure, which enables us to exploit informational flow for the benefit of evacuees. In particular, it would explore the conditions which lead to an efficient evacuation with the help of social structure similar to real life. In other words, this work integrates a new model of social influence with other persisting aspects of interest. The proposed social model is anchored onto the small world phenomenon and the fact that novel information is expected to originate from long distant (or weak) ties, instead of, more frequently encountered strong ties. Moreover, Our work focuses on the question of rationality vs. emotionalism of individuals in a social context while focusing on modeling local interactions. This phenomenon is well researched in social, psychological, and Zia et al. (2021), PeerJ Comput. Sci., DOI 10.7717/peerj-cs.531 3/25 economics domains but relatively less explored in the domain of modeling and simulating evacuation scenarios. This work helps fill this research gap. The rest of the paper is organized as follows. ‘Background and Motivation’ presents the literature review and the contribution of this paper. ‘Models’ provides a detailed description of models including the underlying space and time model. It also presents the decision-making model, which considers individual, social and informational aspects, in detail. The simulation structure and flow using the Overview, Design concepts and Details (ODD) protocol is presented in ‘ODD Protocol’ while ‘Simulation and Results’ is devoted to simulation setup and results. A detailed discussion on the results is presented in ‘Discussion’ while ‘Conclusion’ concludes the paper. BACKGROUND AND MOTIVATION Modelling human behaviour and devising efficient evacuation schemes have been a focus of research, now, for quite some time. Various literature surveys have been carried out summarizing and critically analyzing crowd evacuation models from different perspectives (He et al., 2013; Vermuyten et al., 2016; Bellomo et al., 2016; Shiwakoti, Shi & Ye, 2019; Koponen & Nousiainen, 2016; Zollman, 2012; Murase et al., 2014; Fu, Shi & Wang, 2015; Zheng & Cheng, 2011). ABM has been used extensively in the literature to implement and simulate evacuation scenarios. Researchers have studied the evacuation scenario of crowds in buildings, cities, and classrooms. Some studies presented in Koponen & Nousiainen (2016); Zollman (2012); Murase et al. (2014) have investigated evacuation schemes from the social network aspect while others have studied them using nature-inspired theories and algorithms (Fu, Shi & Wang, 2015; Zheng & Cheng, 2011). Researchers in Wang et al. (2015); Tanimoto, Hagishima & Tanaka (2010); Lo et al. (2006); Altshuler et al. (2005); Wang et al. (2016) have described crowd behaviour and evacuation schemes in detail. Out of this plethora of efforts, we detail a few which are more relevant to our model. Li et al. (2016) studied crowed behaviour in metro using ABM. They studied passengers behaviour in stairs, ticket gates as well as pathways. The goal was to find a safer and faster way to evacuate passengers from metro in case of an emergency. They simulated various evacuation scenarios against different parameters and found that it took an optimal time of three minutes and eight seconds to evacuate all passengers without any injury. However, this goal was achieved with an increased number of ticket gates and the width of the stairs. Game theory was utilized in developing the evacuation model, as passengers think before they move and distribute to different locations. Liu, Jiang & Shi (2016) carried out a detailed study on effectively planning evacuation routes in highly populated buildings such as movie theaters and schools/colleges. They simulated evacuation schemes with four different layouts of classrooms. Changes in the layout of a classroom had an impact on evacuation routes. The students were categorised into two groups, in which one group was provided evacuation instructions while the other was left without any evacuation instructions. The former group of students performed better during the evacuation process compared with the latter group, as they followed quicker and safer evacuation routes. It was Zia et al. (2021), PeerJ Comput. Sci., DOI 10.7717/peerj-cs.531 4/25 also learnt that classrooms with two exits were better for evacuation in an emergency. The overall framework of our model is inspired by the concepts taken from these two papers. Moreover, many research activities come under the domain of ABM such as structural dynamics, agents’ characteristics and decision-making mechanisms, social behavior, and informational flow based on social ties. However, all these efforts are isolated in nature with respect to different aspects of interest. The model, presented in this work, integrates these aspects, which is important to study such scenarios in detail. It is inspired by the work of Wang et al. (2015) and its extension presented in our previous work (Zia et al., 2016; Zia & Ferscha, 2020). Part of our proposed research is also inspired by the work presented in Koponen & Nousiainen (2016), which has demonstrated the impact of appreciation on social interaction and strengthening relationships. The model presented in Wang et al. (2015) differentiated and analyzed the emotional and rational behaviors of agents among a given population. It exploited the generally observed behavior during violent evacuation scenes, represented as emotional behavior that mimics others and the completely opposite rational behavior. The model, we proposed in Zia et al. (2016), questioned Wang et al. (2015) notion of rational behavior and introduced the idea of panic with a hypothesis that a rational agent only, in panic, would do the opposite of mimicking others. It incorporated a formal spatial model with the homophily based social network to capture the panic situation. This model, however, considers only individual characteristics of the evacuees and, thus, the decision about a chosen exit is based on the dynamics of a limited neighborhood. In Zia & Ferscha (2020), we extended this model and incorporated environmental dynamics in it. This extended model realized environmental dynamics using the potential field concept, a collective of the mechanisms used to capture exit events, dealing with content, content spread, and its influence on agents. The extended model, presented in Zia & Ferscha (2020), actually, adopted the concept of guided assistance. In our earlier work, we used a device named, LifeBelt, to provide vibrotactile assistance to the users wearing it, for guiding them towards the safe exit(s) (Ferscha & Zia, 2009). This device is particularly helpful in environments with poor visibility due to blackout or smoke, a potential situation, which might occur during a natural hazard or a technological failure (Ferscha & Zia, 2010). The device was able to actuate appropriate vibra elements based on what was sensed about flow at different exits. It was, thus, producing the locomotive behaviour (short-range mobility) in response to combined effect of directional guidance and social pressure in the neighbourhood. Extensive experimental trials helped us develop models of mobility in the presence of social and cognitive load during evacuations in ambient assisted environments (Zia et al., 2011). The developed model has enough fertility and it was used for city-scale evacuations (Zia & Mitleton-Kelly, 2013). It was also experimented in an industrial setup (Murauer et al., 2019). Reiterating again, in the model presented in Zia & Ferscha (2020), the status information of exits is transmitted to the agents using spatial features of the environment, namely: the potential field. Though, a physical environment, in the current technological advanced era, is certainly capable of transmitting such information. However, this cannot be considered as the social spread of information. To cope with this inconsistency, the model proposed Zia et al. (2021), PeerJ Comput. Sci., DOI 10.7717/peerj-cs.531 5/25 in this paper uses only social spread through people. Also, this paper introduces a new behavioral capability, which tries to capture the strength of long-distance, weak ties supported by the concept of a small-world network, hoping that new information would most probably originate from unfamiliar/unknown corridors. MODELS This section describes various models and the decision-making strategies used in this work. Space and time model This work uses a two-dimensional hypothetical square-shaped simulation world that comprises equal size cells. It is a 51 × 51 cells world, in a grid form, having a total of 2,601 cells, where each cell is represented by an xy-coordinate point. Mobile agents reside on top of cells and they are capable of using the information available at their underlying cells and that of neighboring cells and/or agents, in making their decisions. A neighboring agent can be an agent in the proximity as well as at a distance. Simulation operates in a discrete manner and every agent gets its turn at each timestamp. Agents use incremented state of the world and, thus, make their decisions sequentially. The built-in engine of NetLogo (Wilensky & Rand, 2015) (the software used for simulation in this work), however, provides a fair chance to the agents based on random sequencing at each timestamp. Each simulated space is initially setup with a floor field. As a result of the spread of the floor field, the state variables Hops—having the shortest distance (in terms of hop count) for each exit—and Doms—having the direction towards each of the exits corresponding to the shortest distance—of cells (Netlogo’s patches) are set. For a cell’s neighborhood, we have used Moore’s neighborhood model (Kretz & Schreckenberg, 2007), a cluster of nine neighboring cells around a cell. Figure 1A shows the neighborhood and directions for the cell in the center, exactly the same as used in NetLogo. The spread of floor field starts from exits and disperses outwards following a neighborhood based interaction. When a cell receives information from one of its neighbors, it checks its hop counter of the exit from where it originated. If the new value of hops is smaller than the current value held by the cell, the old value in Hops (for that exit) is replaced by the new value. The direction in Doms (for that exit) is also set to the direction towards the source cell. For example, in the situation shown in Fig. 1A, the information about exit1 to the cell in the center comes from cell at the right and behind and cell at the right. Since the hops to exit1 reported by the cell at the right and behind is less than what is reported by the cell at the right, direction to reach to exit1 is 135 (towards the source cell) and hop count is 1 more than the source cell. The square-shaped simulation world (as shown in Fig. 2) has two exits. The figure also shows the geometric variations for the same two exits. These distributions show an increasing unevenness between exits in terms of the nearness of patches to the exits. Description of different types of patches is needed to completely understand the environmental structure. The patches could be: 1. Normal: A normal patch is walkable and its type is floor, as illustrated as a green patch in Fig. 2. Zia et al. (2021), PeerJ Comput. Sci., DOI 10.7717/peerj-cs.531 6/25 Figure 1 (A) Spread of floor field. (B), (C) Next cell selection mechanism. The textual descriptions of directions assume an agent at the central cell whose current heading (intended direction) is towards the cells at angle labelled as ahead. This diagram also shows Moore’s Neighborhood of size 1 of the cell at the center. Full-size DOI: 10.7717/peerjcs.531/fig-1 Figure 2 Four geometric variations of two exits shown in red with increasing unevenness between exits in terms of the nearness of patches to the exits from (A) to (B) to (C) and to (D) due to obstacles (in gray). The yellow patches are close to the exit on the left and the white patches are close to the exit on right. Full-size DOI: 10.7717/peerjcs.531/fig-2 2. Obstacles: The patches that act as obstacles. The floor field spread does not pass through them and Doms and Hops of these patches, therefore, would have the default values (each −1). Obstacle patches are not walkable and they have type obstacle. Obstacles are shown as gray patches in Figs. 2C and 2D. 3. Exits: The patches that are used as exits. The floor field spread starts from these patches and, therefore, Doms and Hops of these patches have default values (each −1). Exit patches are walkable and they have type exit. These are shown red in Fig. 2. This work uses the microscopic next cell selection strategy based on the concept of Moore’s neighborhood (Kretz & Schreckenberg, 2007). This strategy is about avoiding agents occupying the same cell. An agent chooses one of its eight neighbouring cells depending on the direction, in which it is moving. However, if the desired cell is already occupied, the agent adopts a strategy to choose an unoccupied cell. Figures 1B and 1C shows a dial towards the intended direction labeled with ‘1’. If the cell in this direction is occupied, the agent will choose the direction labeled with ‘2’, and so on. This strategy is based on empirical evidence from the work of Ferscha & Zia (2009). Zia et al. (2021), PeerJ Comput. Sci., DOI 10.7717/peerj-cs.531 7/25 Decision making strategies The decision-making process allows the agents to choose one of the two exits based on the following four strategies: 1. Strategy 1 (S1): which selects the nearest exit. It is based on individual aspect only. 2. Strategy 2 (S2): which selects the exit under the influence of local neighbors. It is based on localized social influence. 3. Strategy 3 (S3): which selects the exit under the influence of exit potential. It is based on technological influence. 4. Strategy 4 (S4): which selects an exit under the influence of both local and distant neighbors. It is based on small-world social influence. These strategies are built upon the space and time model presented before and all of them except S1 support both emotional and rational types of agents. S1 acts as a base case, in which all agents proceed towards the nearest exit irrespective of type. In S2, only rational agents are influenced by a localized social circle. S3 uses exits’ potential along with the localized social circle. A comparison of the strategies S1, S2, and S3 is presented in Zia & Ferscha (2020). However, S3 described in Zia & Ferscha (2020) is updated in this work, and the updated version of this strategy transmits exits’ status information through agent-to-agent communication. The potential field of exits’ still remains a spatial feature, however, the spread in revised S3 is performed by agents instead of patches. S4 is an extended version of strategy S2. S4 has introduced a new behavioral capability, which attempts to capture the strength of long-distance weak ties supported by the concept of a small-world network. Thus, a distant agent, who is successfully moving would help a rational agent in panic. The following subsections describe these strategies in detail. Strategy 1 (S1) All agents in S1 proceed to the nearest exit, which is one of the two exits having a minimum value of the parameter Hops. This information is available at each patch, which is calculated using the floor field spread mechanism when the simulation setup is executed. This is a static strategy and used to compare it with the remaining strategies that allow dynamic changes between exits. Strategy 2 (S2) The strategy S2 is based on game theory. Player 1 is the agent, who is making a decision while Player 2, a random neighbour from Moore’s neighbourhood, is the influencer. The logic behind S2 is implemented in Algorithm 1, where Ae represents the current exit of agent A, Ne represents the current exit of neighbor N, and Ap shows the panic index of agent A. If the behaviour of an emotional agent is the same as the neighbouring agent, then, the payoff is maximum and, thus, the agent retains its behaviour and continues moving towards the current exit. The payoff remains maximum, even, when the behavior of random neighbor is the opposite of the agent behavior. This is the reason that these agents are ‘‘emotional’’. However, there are 10% chances that an agent may change its behavior, if it is in panic, which ensures a random deviation. This means that the game is not really played by emotional but rational agents. Zia et al. (2021), PeerJ Comput. Sci., DOI 10.7717/peerj-cs.531 8/25 Algorithm 1: The Logic behind Strategy 2 (S2) Agent A is making a decision Agent N is a random local agent Ae ← current exit of A Ap ← current panic index of A Ne ← current exit of N if A is emotional then if Ae ! = Ne then if Ap > 0 then Ae = Ne [with probability 0.1] end end end if A is rational then if Ap <= 0 then if Ae == Ne then Ae = Aalternate [with probability 0.1] else Ae = Aalternate [with probability 0.9] end else if Ae == Ne then Ae = Aalternate end end end The strategy for a rational agent is dependent on whether it is in a state of panic or not. If an agent is not in a panic, it will be a source of strengthening its belief on current behavior, when the randomly selected neighbor has the same behavior. This case is represented as a maximum payoff. However, there are 10% chances that the agent may change its behavior, thus ensuring a random deviation. On the other hand, it will be a source of weakening an agent’s belief on the current behavior, when the randomly selected neighbor has the opposite behavior. This case is represented by a minimum payoff. However, there are 90% chances that the Agent may change its behavior leaving 10% chances for random deviation. If the rational agent is in panic, we argue about entirely opposite behavior, in which the same behavior represents the minimum payoff. This directs the agent to change its behavior but strictly this time. Similarly, if the behavior of the randomly selected neighbor is different from the agent’s own behavior, the agent will strictly retain its recent behavior as no change of behavior is required in this case. Strategy 3 (S3) The decision taken by an agent using strategy S3 also depends on the agent being in panic or not. This model is only applicable to the rational agents, and the emotional agents will keep following the nearest exit. S3, however, is based on two additional conjunctures, regarding the spatial change, presented in Tan, Hu & Lin (2015). The spatial change could have a negative or positive effect. The effect is negative when detouring occurs due to the lack of knowledge. On the other hand, when exit utilization is balanced, the spatial change has a positive effect. S3 (like S2) is based on systematic rational thinking, in which the factors of panic, social influence, and exits’ potential plays a major role in deciding an exit. Panic is given more weight (0.5) than the social influence and exits’ potential, which are assigned 0.25 weight each. Social influence and exits’ potential based on positive or negative role both Zia et al. (2021), PeerJ Comput. Sci., DOI 10.7717/peerj-cs.531 9/25 Algorithm 2: The Logic behind Strategy 3 (S3) Agent A is making a decision Agent N is a random local agent Ae ← current exit of A Ap ← current panic index of A Ne ← current exit of N PVL ← Potential of left exit at the current patch PVR ← Potential of right exit at the current patch Ep ← the exit with more potential among the above two if A is rational and Ap <= 0 then if Ae == Ne then if Ae ! = Ep then Ae = Aalternate [with probability 0.25] end else if Ae == Ep then Ae = Aalternate [with probability 0.25] else Ae = Aalternate [with probability 0.50] end end end if A is rational and Ap > 0 then if Ae == Ne then if Ae == Ep then Ae = Aalternate [with probability 0.50] else Ae = Aalternate [with probability 0.75] end else if Ae ! = Ep then Ae = Aalternate [with probability 0.75] else Ae = Aalternate end end end decide these weights. For example, if an agent is not in panic (represented with a weight ’0′) but having a negative input of social influence (0.25) and having a negative input of exits’ potential (0.25), would result in a total weight of 0.50—showing the probability of switching of exit for the agent. The logic behind S3 is given in Algorithm 2 and it is explained as follows. The current behavior of a rational agent A, in terms of left or right exit, is represented as Ae . Agent A chooses a neighbor N, whose current behavior is represented as Ne . The relative potential of the patch underneath the agent A is represented as Ep . Since there are two exits, the value of Ep would be either left or right. The exit selection is based on the greater value among the two for left and right, which represent a more influential exit. When A is not in panic, and Ae is equal to Ne , then, the agent is considered quite satisfied. The agent has nothing to worry about when Ae also equals Ep and no change is needed in this case. However, if Ae is not the same as Ep , then, the agent has something to worry about but two factors (no panic and neighbor’s influence) are still in its favor. Hence, it would only make a change in the exit plan with a 25% probability. When A is not in panic but Ae is not equal to Ne , then, the agent A has something to worry about. The agent is satisfied, if Ae is the same as Ep , as its own behavior and exits’ potential are both pointing towards the same exit. However, the negative neighbor’s influence would force the agent to switch to the other exit but only with 25% probability. If Ae is not equal to Zia et al. (2021), PeerJ Comput. Sci., DOI 10.7717/peerj-cs.531 10/25 Ep , then, the agent has a worry, because both the neighbor’s influence and exits’ potential are against the current behavior. Since the major factor of panic does not exist, the agent would make a switch exit change with a 50% probability. Before the behavior of the agent A in panic is described, it is, worthwhile, stating that A is already in a big worry due to panic. When A is in panic and Ae is equal to Ne , it may exaggerate its worry due to neighbors also moving to the current exit, which seems counter-productive. Now, if Ae is equal to Ep , this may give the agent a little relief and, hence, it would switch an exit with a 50% probability. When Ep is not in favor (means it’s not the same as Ae ), the agent would switch an exit with a 75% probability. If A is in panic and Ae equals Ne , then, there is something to worry about. Now, if Ae is also the same as Ep , then, the agent may decide to change the exit but with a little motivation. In this case, the agent would switch an exit with a 50% probability. Otherwise, the agent would be more motivated about changing the current exit with 75% probability. When A is in panic and Ae is not the same as Ne , then, the agent has to worry a lot. If Ae is also not equal to Ep , the agent would change the current exit with 75% probability, otherwise, it will definitely switch the exit. Algorithm 3: The Logic behind Strategy 4 (S4) Agent A is making a decision Agent N is a random local agent Agent M is agent with maximum mobility index (including both local and distant agents) Ae ← current exit of A Ap ← current panic index of A Ne ← current exit of N Me ← current exit of M if A is emotional then if Ae ! = Ne then if Ap > 0 then Ae = Ne [with probability 0.1] end end end if A is rational then if Ap <= 0 then if Ae == Ne then Ae = Aalternate [with probability 0.1] else Ae = Aalternate [with probability 0.9] end else Ae = Me end end Strategy 4 (S4) Strategy S4 is the same as S2 with one exception that the alternate exit in S4 is the exit of one of the neighbouring agents, which has the highest Mobility Index (MI). The neighbouring agents include both local and distant neighbours. A long distant random link is established for each rational agent in this strategy. The logic of strategy S4 is given in Algorithm 3. Zia et al. (2021), PeerJ Comput. Sci., DOI 10.7717/peerj-cs.531 11/25 Figure 3 Screen-shots of the random simulation outcome at iteration 25, in case of (A) all agents only going to the nearest exit, (B) agents under the influence of locality driven rationality, (C) rational agents under the influence of social as well as potential field, and (D) rational agents under the influence of locality as well as distantly driven rationality. The accompanying graphs indicate the parameters for which the strategies are compared. Full-size DOI: 10.7717/peerjcs.531/fig-3 Differentiating the strategies in graphical form Figure 3 illustrates the differences between the four strategies along with the visualisation of the simulation space. It shows agents’ positions at iteration 25 for all the four situations. Rational agents are represented by circles while emotional by triangles. Different colors are used to represents various types of states, which are described later. In Fig. 3A, all the agents are going only to the nearest exit, as directed by S1. Figure 3B shows the exit preferences of agents towards exits under the influence of locality driven rationality, as specified by S2. Exit preferences of rational agents under the influence of social as well as potential field in a rational way are provided in Fig. 3C, as proposed in S3. The emotional agents in this strategy follow the nearest exit. Figure 3D illustrates the exit distribution of rational agents under the influence of social rationality, as suggested by S4. In this case, a social circle include distant agents and emotional agents use the nearest exit. The comparison of Fig. 3B and 3D and with Figs. 3A and 3C shows that quite a few rational agents moved from the nearest left exit to the right exit under the social influence. Detailed analysis of this will be presented in the following sections. ODD PROTOCOL Overview, Design concepts and Details (ODD) is a scientific protocol (Grimm, Polhill & Touza, 2017), used, for the documentation of ABM simulations. This section is dedicated to present, the work presented in this paper, in terms of ODD protocol. Zia et al. (2021), PeerJ Comput. Sci., DOI 10.7717/peerj-cs.531 12/25 An overview The purpose and patterns The purpose behind the development of the proposed model is to evaluate the performance of an evacuating crowd from a hypothetical square-shaped space with two exits. Entities, state variables and scales Evacuees are the only agents. The state variables used by agents are: current exit, state, and MI. The current exit provides the exit an agent is currently heading to. The current cognitive state of an agent is represented by the ‘‘state’’ variable. Similarly, MI shows the current quality of mobility acquired by an agent. Patches use the state variables named: Structure_type, Doms, Hops, PVL, and PVR. The structure type shows the type of a patch being either normal (walkable), obstacle, or exit. The variable ‘Doms’ provides the directions of nearest paths for both left and right exits while ‘Hops’ maintains the distance (hop count) of nearest paths for both left and right exits. The PVL and PVR represent the potential field values for the left and right exits. Process overview and scheduling Mobility and decision-making happens at each timestamp of the simulation, concurrently, for all agents. When an agent exits through one of the two exits, it is no more part of the simulation. The number of agents in the simulation is used for terminating the simulation. The procedure, therefore, counts the total number of agents in the simulated space at each timestamp before it starts the rest of the activities. The simulation is stopped when this agent count becomes zero. Otherwise, the decision-making process allows an agent to determine the current exit to follow by taking strategy x as an input. Strategy x refers to one of the four strategies, described earlier. Once an exit is selected, the agent uses the procedure heading adjustment to adjust its direction using the floor field information available at the cell level. The agent, then, moves, whose dynamics are explained shortly. If the simulation is run for strategy S3, the potential field information is also updated after the move. Figure 4 provides an abstract view of the flow between various activities of the simulation. The initial state of agents is set to WALK. The agents perform microscopic move (described in ‘Space and time model’), concurrently, at each iteration. In response to a move, three indices (namely the moving, waiting and panic) get updated based on the current state of an agent. If the state of an agent is WALK, the moving index is incremented by 1 while the waiting and panic indices are reset to 0. If the state of an agent is WAIT, the waiting index is incremented by 1 while the rest of the two indices are reset to 0. Similarly, when the state of an agent is PANIC, the panic index is incremented by 1 while the other two indices are reset to 0. When the agent is able to move, it moves and sets its state to WALK for the next iteration. However, when it fails to move, it may enter either into WAIT or PANIC state depending on the waiting index. When the waiting index is less than a threshold, its state is set to WAIT, otherwise, it is set to PANIC. Figure 5 provides the detailed process of move operation and its related transitions. When an agent reaches an exit, it is deleted from the simulation. Global variables (which are detailed in ‘Initializations’ and used for reporting purposes) are updated at each iteration. Zia et al. (2021), PeerJ Comput. Sci., DOI 10.7717/peerj-cs.531 13/25 Figure 4 An abstract view of the simulation process and scheduling executed at each iteration by each agent. Full-size DOI: 10.7717/peerjcs.531/fig-4 Figure 5 The detailed view of move process and state transition system executed at each iteration by each agent. Full-size DOI: 10.7717/peerjcs.531/fig-5 Design concepts Basic principles The proposed model addresses the basic principle of using discrete time and space simulation for evacuation purposes. Different models of agents’ interaction are tested to perform a comparative analysis of various individual, informational and social behaviors in Zia et al. (2021), PeerJ Comput. Sci., DOI 10.7717/peerj-cs.531 14/25 terms of evacuation efficiency. The basic objective behind the adaptive behavior of agents is to avoid panic and keep them moving. Emergence S1 is base strategy for evacuation, in which all the agents are moving towards the nearest exit. Strategies S2 to S4 use social and informational behavior for the emergence of a more distributed and less congested population of agents, with the help of actions of most of the rational agents, while evacuating a space with two exits. We believe that the situations/cases where the people are distributed as much uniform as possible among the two exits and when they are not panicking (actually a by-product of lesser congestion) will have fewer chances of a disaster. The purpose of the three strategies (S2 to S4) is to find conditions, which would let a favourable flow of evacuees to emerge. Adaptation The adaptation of behavior is a central part of the model. Rational agents adapt their behavior, in response to social pressure in the proximity and distant neighbors in S2 and S4. In S3, however, they adapt their behavior based on social proximic pressure as well as potential map diffused by exiting agents. Game-theoretic approach is adopted by all adaptation mechanisms, in which different situations are evaluated based on relative payoffs. Sensing Sensing plays a vital role in making exit decisions. To make rational decisions, agents need to sense the current state of the surrounding neighbors and their exit choices. In case of strategy S3, the agents require to get potential field information. Similarly, strategy S4, further, suggests sensing the MI and exit choice of a distant neighbor. In microscopic move, the agents sense their locality for possible collisions and avoid them by using the corresponding collision avoidance mechanism. Interaction Different types of interactions between the agents and patches such as agent to agent, agent to patch, patch to agent, and patch to patch enable various activities. They not only setup floor field for the nearest direction and distance to both exits but also help spreading and decaying informational content (the potential field) by agents and patches. They materialize social interaction among agents in the neighborhood and enable distant interaction among the agents. Furthermore, they support interactions between agents and environmental structure. Observations Simulation results are evaluated based on three quantitative measures described below: 1. Exit Time: collects the number of iterations (simulation ticks) needed by all agents to exit through the designated exit points. 2. Exit Usage: shows the number of agents, who left the system through two different exits. 3. Agents in Panic: collects the maximum number of agents in panic during the simulation. Zia et al. (2021), PeerJ Comput. Sci., DOI 10.7717/peerj-cs.531 15/25 Details Initialization This work needs to set up patches, agents, and global reporting variables. Therefore, the initialisation is divided into three groups. The default structural type of all the patches’ is normal (walkable). To setup exits, each exit patch is assigned a name: left and right in this work. To setup non-walkable obstacles, the environment type parameter obstacle is used. The last step is to create the floor field value of the patches. Once the patches are initialized, we initialize agents as follows. All agents are randomly assigned to unique walkable patches without overlapping. The moving, waiting and panic indices are, then, all set to 0. The current exit for all agents is set to the nearest exit, which are collected from the ‘Hops’ variable of their underlying patches. The state of all agents is set to WALK and their type is set to emotional while the type of required number of agents is set to rational. This work uses a number of global reporting variables and they are all initialized to 0. They are the tick counter indicating the total simulation time, the number of agents who exit from left and right exits, and the number of rational and emotional agents who enter into the PANIC state at some time during the simulation. Sub models Three types of sub-models are used: update Update Mobility Index,, Spread and Decay, in this work. update MI: Mobility Index (MI) represents the quality of mobility, obtained by an agent, with the passage of time. Initialized with a value of 0, MI is updated at each iteration based on the current state of the agent, as follows. It is incremented with a 1 when the current state is MOVE and decremented with a 1 when the state is WAIT. MI = MI − (MI × sensitivity) when the current state of the agent is PANIC. Sensitivity has a static 0.5 value. MI value is used in S4, in which, it is required to have an agent M with the highest MI value in the neighborhood of a rational agent. spread: This sub-model is responsible for spreading the potential field of an exit (PVX ). When an agent reaches an exit, the patch variable PVX of all neighboring patches is set to 1, where X represents both left and right exits. For other agents, if a patch (say P) on which an agent is residing, has PVX greater than 0, the list of all neighboring patches would be potentially updated. For every neighboring patch N, if its PVX is less than PVX of P, PVX of N will be decremented by a fraction (one-tenth of PVX of P). Hence, starting from the potential field value of 1, the exit announcements would be spread out with lesser and lesser intensity, provided that there is an agent encountered. decay: This sub-model just lets the potential value of patches decay with the passage of time. Zia et al. (2021), PeerJ Comput. Sci., DOI 10.7717/peerj-cs.531 16/25 SIMULATION AND RESULTS Simulation setup This work is simulated with a 2D bounded grid space that comprises 51 × 51 cells. The simulation starts with a number of agents in the space and stops when the space becomes empty. However, the density, in this work, is kept static to a value of 1,000. The simulated space has two exits named left and right. This study investigates different scenarios based on three parameters: Environment Type, Strategy, and % of rational population. Thus, each case is made up of the tuple [Environment type, strategy, and % of rational agent in the population]. It uses four different types of environments namely E1, E2, E3, and E4, which correspond to Figs. 2A to 2D, in the given order. Four different evacuation strategies S1, S2, S3 and S4 are used, which are explained in detail in ‘Models’. This work uses four different percentages (5%, 10%, 15%, and 20%) of rational agents among the total population. Every simulation case is repeated 50 times and their average values are then used for comparison purposes. Simulation results Exit time Figure 6 presents simulation results for different scenarios. It shows that the exit time, generally, increases with an increase in the complexity of the environment. However, the comparative analysis of strategies suggests grouping S1 and S3 as one while S2 and S4 as the second group. The former group is much more efficient (in terms of exit time of the population) than the latter. However, this difference vanishes in more challenging environments, where the number of rational agents is only 5% of the total population. As the percentage of rational agents increases, the former group starts to become more and more efficient. In spite of the fact that S1 is the best (in terms of exit time) as all the agents are directed towards the nearest exit, irrespective of their panic and emotional state, S3 is comparable to it. In fact, it is slightly better in complex environments. Redirection of agents to alternate exit, in the response to social influence, definitely increases the exit time of the population. However, it must be emphasized that in many situations, for a herding crowd, the factor of panic in people can be a more important factor than exit time. When S2 is compared with S4, minor improvements are observed in some settings, however, they can be ignored. Though the effect of population density on exit time is not relevant here, however, it was observed in our previous work (Zia & Ferscha, 2020) that exit time increases with population density. Exit usage Quantitatively, the reason of the exit time, being good or bad, is due to the change of agents’ preferences regarding the exit. Based on the geometric consideration of simulation space, the exit time is good, when many agents change their exit from left to right. The columns in Fig. 7 show the number of agents actually ended up at the left exit for various scenarios, at the end of the simulation. We seek to have a lesser and lesser number here, which is equivalent to equal distribution of agents across both exits. However, the distribution Zia et al. (2021), PeerJ Comput. Sci., DOI 10.7717/peerj-cs.531 17/25 Figure 6 Simulation time (in terms of ticks) for all strategies against four different environments. Illustrating simulation responses for: (A) 5%, (B) 10%, (C) 15%, and (D) 20% rational agents. Full-size DOI: 10.7717/peerjcs.531/fig-6 of exit usage among left and right exits, generally deteriorates, with an increase in the complexity of the environment. The comparative analysis of strategies suggests grouping S1 and S3 into one while S2 and S4 into the second group. The former group is less efficient (in terms of agents’ distribution across the exits) than the latter. The difference between the two groups is more significant in more challenging environments. It can also be seen that S2 performs marginally better than S4. More equal distribution of agents across both exits is expected to decrease the panic in the population, which is analyzed next. Agents in panic Reduction of panic among the crowd was one of the main purposes of this work. Figure 8 illustrates simulation trends for various cases based on different strategies, environments, and rational agents. It is clear that S3 has always less number of panicking agents as compared to S1. However, the number of panicking agents increased with an increase in the complexity of an environment. It can be observed that S1 performs the worst followed by S3 in all cases irrespective of the percentage of rational agents among the population. Strategies S2 and S4 have a similar and improved performance (against both S1 and S3) as they have a reduced number of agents in panic during all the cases. The comparative analysis of these strategies suggest placing S1 and S3 in one group while S2 and S4 in the second group. The former group is less efficient, in terms of panicking agents, than the latter one. Though, the difference between the two groups is consistent in all cases, Zia et al. (2021), PeerJ Comput. Sci., DOI 10.7717/peerj-cs.531 18/25 Figure 7 Number of agents, who exit the space through the left exit for all strategies against four different environments. Illustrating simulation responses for (A) 5%, (B) 10%, (C) 15%, and (D) 20% rational agents. Full-size DOI: 10.7717/peerjcs.531/fig-7 however, it is more obvious in complex challenging environments. In general, S2 performs very similar to S4 in all cases. DISCUSSION The main objective of the development of different strategies was seeking a possible emergence of favorable situation for an evacuating crowd at a global level, where people are either emotional or rational, and the environment has different levels of complexity. The results are evaluated based on three measures: (i) time to exit, (ii) distribution of exits’ usage, and (iii) population in panic. S1 is the base case but it was supposed to perform the best in terms of exit time. Ignoring the case, where the %age of rational agents is really low (5%), S1 was always better than other strategies but specifically in the case of simple environments. S3 was also not that far from S1 in all situations. Moreover, S1 was supposed to perform the worst in terms of exits’ usage and panic. It turned true, especially, for the complex environments. S3 was again very similar to S1 in all these cases. However, S3 performed, consistently, much better than S1 in terms of panic in all cases. Hence, it can be inferred that strategy S3, which relies on rich informational content to be transmitted through potentially peered ICT devices, are not of great use, in general. However, it could be beneficial in those evacuation scenarios, in which the situation is not that dangerous and it is performed in a structured way, where the people know each other. This typically applies to residential buildings, factories, and Zia et al. (2021), PeerJ Comput. Sci., DOI 10.7717/peerj-cs.531 19/25 Figure 8 Maximum %age of rational agents in panic for all strategies against four different environments. Illustrating simulation responses for (A) 5%, (B) 10%, (C) 15%, and (D) 20% rational agents. Full-size DOI: 10.7717/peerjcs.531/fig-8 offices with complex geometries and where there is no disastrous situation such as some chemical leakage, which is not that harmful. In fact, the success of S3 depends on this as people without an apparent herding situation would most probably be able to pay attention to the information available to them. In other words, strategy S3 can be applied in those situations, where we are confirmed that people would behave rationally and would not lose their temper (as a whole). It can be facilitated by providing information such as ’exit to take’, in form of announcements and exit signs, to the evacuees with the help of some technology-based on exit dynamics and knowledge about their impact. Panic is related to the complexity of the environment. When the geometry of the environment is more complex, there would be more panic. It is reiterated here that a person gets panic when he/she is unable to proceed towards one of the exits for a longer period of time. The complexity of the environment, therefore, makes the panic issue severe. The people in the strategies S2 and S4, which are social in nature, get less panicked. S2 is just based on human behavior, whereas, S4 also considers distant connectivity. It has been observed that people do ask for help from their distant friends when they get stuck in an emergency situation. Results have suggested that the response of a simple mechanism of doing against what others are doing, in a panic situation, emerges in a much better shape. However, this is against the psyche of a herd. Still, people may be trained not to be in a panic and act rationally to give them more chances of survival. It was learnt that S4 has no substantial advantage over S2. This might be due to the nature of environmental structures used in this work. It, however, could turn out to be advantageous over S2, when Zia et al. (2021), PeerJ Comput. Sci., DOI 10.7717/peerj-cs.531 20/25 tested against more complex and dynamically changing environments. But, this has yet to be seen and left as a future work. CONCLUSION This paper proposed an agent-based model, providing a framework, to explore individual, social and technological aspects of crowd evacuation. It combined many features, which are related to an evacuating crowd and they are categorized as individual, social and technological features. Models of these (called strategies) were pivoted onto a common agent-based modeling framework and a grid-based hypothetical environment. By simulating these models, an insight into the effectiveness of several interesting evacuation scenarios is provided. Simulation results revealed that there is no substantial advantage of informational overload on a crowd in an evacuation. This may only marginally work in favor of people, and probably only in those situations, in which there are fewer chances of herding. In more serious situations, people should be left alone to decide. But, they can be taught not to panic during the drills so that they are better able to make their decisions, solely based on the dynamics of their surroundings. Also, distant connectivity has no apparent advantage as well. But, this can well be due to the nature of the environmental structures that were used. With more complex and dynamically changing environments, information coming from distant ties may turn out to be vital. Limitations: Although the behavioral assumptions about people in panic during a herding situation have an underpinning on behavioral theories, empirical evidences and models progression but building blocks of human behavior and environmental dynamics used to develop the model, in this paper, are not explicitly verified. ADDITIONAL INFORMATION AND DECLARATIONS Funding This work was supported by FFG funded Pro2Future project under contract No.6112792. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Grant Disclosures The following grant information was disclosed by the authors: FFG: 6112792. Competing Interests The authors declare there are no competing interests. Author Contributions • Kashif Zia conceived and designed the experiments, performed the experiments, performed the computation work, prepared figures and/or tables, authored or reviewed drafts of the paper, and approved the final draft. Zia et al. (2021), PeerJ Comput. Sci., DOI 10.7717/peerj-cs.531 21/25 • Umar Farooq performed the experiments, performed the computation work, authored or reviewed drafts of the paper, and approved the final draft. • Muhammad Shafi analyzed the data, prepared figures and/or tables, and approved the final draft. • Alois Ferscha conceived and designed the experiments, authored or reviewed drafts of the paper, concept, and approved the final draft. Data Availability The following information was supplied regarding data availability: The code of an agent-based model developed using NetLogo is available in the Supplementary File. Supplemental Information Supplemental information for this article can be found online at http://dx.doi.org/10.7717/ peerj-cs.531#supplemental-information. REFERENCES Almeida JE, Rosseti RJ, Coelho AL. 2013. Crowd simulation modeling applied to emergency and evacuation simulations using multi-agent systems. ArXiv preprint. arXiv:1303.4692. Altshuler E, Ramos O, Núñez Y, Fernández J, Batista-Leyva A, Noda C. 2005. Symmetry breaking in escaping ants. The American Naturalist 166(6):643–649 DOI 10.1086/498139. Batty M. 2001. Agent-based pedestrian modeling. Environment and Planning B: Planning and Design 28(3):321–326 DOI 10.1068/b2803ed. Bellomo N, Clarke D, Gibelli L, Townsend P, Vreugdenhil B. 2016. Human behaviours in evacuation crowd dynamics: from modelling to ‘‘big data’’ toward crisis management. Physics of Life Reviews 18:1–21 DOI 10.1016/j.plrev.2016.05.014. Cova TJ, Johnson JP. 2003. A network flow model for lane-based evacuation routing. Transportation Research Part A: Policy and Practice 37(7):579–604 DOI 10.1016/S0965-8564(03)00007-7. Ferscha A, Zia K. 2009. Lifebelt: silent directional guidance for crowd evacuation. In: 2009 international symposium on wearable computers. Piscataway: IEEE, 19–26. Ferscha A, Zia K. 2010. Lifebelt: crowd evacuation based on vibro-tactile guidance. In: IEEE pervasive computing, volume 9, issue 4. Piscataway: IEEE, 33–42 DOI 10.1109/MPRV.2010.83. Fu H, Shi Y, Wang Y. 2015. Emergence of moving pattern in a collective game. In: Journal of physics: conference series. Bristol: IOP Publishing, 012020. Garcimartín A, Zuriguel I, Pastor J, Martín-Gómez C, Parisi D. 2014. Experimental evidence of the ‘‘faster is slower’’ effect. Transportation Research Procedia 2:760–767 DOI 10.1016/j.trpro.2014.09.085. Grimm V, Polhill G, Touza J. 2017. Documenting social simulation models: the ODD protocol as a standard. In: Simulating social complexity. Cham: Springer, 349–365. Zia et al. (2021), PeerJ Comput. Sci., DOI 10.7717/peerj-cs.531 22/25 He G-Q, Yang Y, Chen Z-H, Gu C-H, Pan Z-G. 2013. A review of behavior mechanisms and crowd evacuation animation in emergency exercises. Journal of Zhejiang University Science C 14(7):477–485 DOI 10.1631/jzus.CIDE1301. Helbing D, Farkas I, Vicsek T. 2000. Simulating dynamical features of escape panic. Nature 407(6803):487–490 DOI 10.1038/35035023. Helbing D, Johansson A, Al-Abideen HZ. 2007. Dynamics of crowd disasters: an empirical study. Physical Review E 75(4):046109 DOI 10.1103/PhysRevE.75.046109. Kania J, Kramer M. 2013. Embracing emergence: how collective impact addresses complexity. California: Stanford Social Innovation Review, Leland Stanford Jr. University. Koponen IT, Nousiainen M. 2016. Formation of reciprocal appreciation patterns in small groups: an agent-based model. Complex Adaptive Systems Modeling 4(1):24 DOI 10.1186/s40294-016-0035-6. Kretz T, Grünebohm A, Schreckenberg M. 2006. Experimental study of pedestrian flow through a bottleneck. Journal of Statistical Mechanics: Theory and Experiment 2006(10):P10014 DOI 10.1088/1742-5468/2006/10/P10014. Kretz T, Schreckenberg M. 2007. Moore and more and symmetry. In: Pedestrian and evacuation dynamics 2005. Berlin: Springer, 297–308. Li Z-Y, Tang M, Liang D, Zhao Z. 2016. Numerical simulation of evacuation in a subway station. Procedia Engineering 135:616–621 DOI 10.1016/j.proeng.2016.01.126. Liu R, Jiang D, Shi L. 2016. Agent-based simulation of alternative classroom evacuation scenarios. Frontiers of Architectural Research 5(1):111–125 DOI 10.1016/j.foar.2015.12.002. Lo SM, Huang H-C, Wang P, Yuen K. 2006. A game theory based exit selection model for evacuation. Fire Safety Journal 41(5):364–369 DOI 10.1016/j.firesaf.2006.02.003. Macal CM. 2016. Everything you need to know about agent-based modelling and simulation. Journal of Simulation 10(2):144–156 DOI 10.1057/jos.2016.7. Murase Y, Török J, Jo H-H, Kaski K, Kertész J. 2014. Multilayer weighted social network model. Physical Review E 90(5):052810 DOI 10.1103/PhysRevE.90.052810. Murauer M, Jungwirth F, Anzengruber B, Abbas A, Ahmad A, Azadi B, Cho J, Ennsbrunner H, Ferscha A, Gerhard D, Gollan B, Haslgrübler M, Selymes J, Sopidis G, Stütz M, Weißenbach P. 2019. A task-independent design and development process for cognitive products in industrial applications. In: Proceedings of the 12th ACM international conference on pervasive technologies related to assistive environments, New York: ACM, 358–367 DOI 10.1145/3316782.3322748. Pelechano N, Malkawi A. 2008. Evacuation simulation models: challenges in modeling high rise building evacuation with cellular automata approaches. Automation in Construction 17(4):377–385 DOI 10.1016/j.autcon.2007.06.005. Radianti J, Granmo O-C, Bouhmala N, Sarshar P, Yazidi A, Gonzalez J. 2013. Crowd models for emergency evacuation: a review targeting human-centered sensing. In: 2013 46th Hawaii international conference on system sciences. Piscataway: IEEE, 156–165. Zia et al. (2021), PeerJ Comput. Sci., DOI 10.7717/peerj-cs.531 23/25 Reilly AC, Dillon RL, Guikema SD. 2018. Agent-based models as an integrating boundary object for interdisciplinary research. In: Risk analysis - an official publication of the society for risk analysis. New Jersey: Wiley Online Library DOI 10.1111/risa.13134. Schadschneider A, Klingsch W, Klüpfel H, Kretz T, Rogsch C, Seyfried A.. 2011. Evacuation dynamics: empirical results, modeling and applications. In: Extreme environmental events. New York: Springer, 517–550. Von Schantz A, Ehtamo H. 2015. Spatial game in cellular automaton evacuation model. Physical Review E 92(5):052805 DOI 10.1103/PhysRevE.92.052805. Seyfried A, Passon O, Steffen B, Boltes M, Rupprecht T, Klingsch W. 2009. New insights into pedestrian flow through bottlenecks. Transportation Science 43(3):395–406 DOI 10.1287/trsc.1090.0263. Sharpanskykh A, Zia K. 2014. Understanding the role of emotions in group dynamics in emergency situations. In: Transactions on computational collective intelligence XV. Berlin: Springer, 28–48. Shiwakoti N, Shi X, Ye Z. 2019. A review on the performance of an obstacle near an exit on pedestrian crowd evacuation. Safety Science 113:54–67 DOI 10.1016/j.ssci.2018.11.016. Tan L, Hu M, Lin H. 2015. Agent-based simulation of building evacuation: combining human behavior with predictable spatial accessibility in a fire emergency. Information Sciences 295:53–66 DOI 10.1016/j.ins.2014.09.029. Tanimoto J, Hagishima A, Tanaka Y. 2010. Study of bottleneck effect at an emergency evacuation exit using cellular automata model, mean field approximation analysis, and game theory. Physica A: Statistical Mechanics and its Applications 389(24):5611–5618 DOI 10.1016/j.physa.2010.08.032. Vermuyten H, Beliën J, De Boeck L, Reniers G, Wauters T. 2016. A review of optimisation models for pedestrian evacuation and design problems. Safety Science 87:167–178 DOI 10.1016/j.ssci.2016.04.001. Wagner N, Agrawal V. 2014. An agent-based simulation system for concert venue crowd evacuation modeling in the presence of a fire disaster. Expert Systems with Applications 41(6):2807–2815 DOI 10.1016/j.eswa.2013.10.013. Wang S, Cao S, Wang Q, Lian L, Song W. 2016. Effect of exit locations on ants escaping a two-exit room stressed with repellent. Physica A: Statistical Mechanics and its Applications 457:239–254 DOI 10.1016/j.physa.2016.03.083. Wang T, Huang K, Cheng Y, Zheng X. 2015. Understanding herding based on a coevolutionary model for strategy and game structure. Chaos, Solitons & Fractals 75:84–90 DOI 10.1016/j.chaos.2015.02.008. Wilensky U, Rand W. 2015. An introduction to agent-based modeling: modeling natural, social, and engineered complex systems with NetLogo. Cambridge: MIT Press. Yuan W, Tan KH. 2007. An evacuation model using cellular automata. Physica A: Statistical Mechanics and its Applications 384(2):549–566 DOI 10.1016/j.physa.2007.05.055. Zheng X, Cheng Y. 2011. Modeling cooperative and competitive behaviors in emergency evacuation: a game-theoretical approach. Computers & Mathematics with Applications 62(12):4627–4634 DOI 10.1016/j.camwa.2011.10.048. Zia et al. (2021), PeerJ Comput. Sci., DOI 10.7717/peerj-cs.531 24/25 Zia K, Ferscha A. 2020. An agent-based model of crowd evacuation: combining individual, social and technological aspects. In: Proceedings of the 2020 ACM SIGSIM conference on principles of advanced discrete simulation. New York: ACM, 129–140. Zia K, Ferscha A, Din A, Farooq U. 2016. Quantifying the interplay of emotions and rationality in herding: a game-theoretic simulation study. Adaptive Behavior 24(2):119–129 DOI 10.1177/1059712316635719. Zia K, Ferscha A. 2013. City scale evacuation: a high-performance multi-agent simulation framework. In: Mitleton-Kelly E, ed. Co-evolution of intelligent sociotechnical systems. Understanding complex systems, Berlin: Springer, 239–293 DOI 10.1007/978-3-642-36614-7_10. Zia K, Riener A, Ferscha A, Sharpanskykh A. 2011. Evacuation simulation based on cognitive decision making model in a socio-technical system. In: 2011 IEEE/ACM 15th international symposium on distributed simulation and real time applications. New York: ACM, 98–107. Zollman KJ. 2012. Social network structure and the achievement of consensus. Politics, Philosophy & Economics 11(1):26–44 DOI 10.1177/1470594X11416766. Zia et al. (2021), PeerJ Comput. Sci., DOI 10.7717/peerj-cs.531 25/25