Research Article | Open Access

Ying Liu, Cheng Sun, Yiming Bie, "Modeling Unidirectional Pedestrian Movement: An Investigation of Diffusion Behavior in the Built Environment", *Mathematical Problems in Engineering*, vol. 2015, Article ID 308261, 6 pages, 2015. https://doi.org/10.1155/2015/308261

# Modeling Unidirectional Pedestrian Movement: An Investigation of Diffusion Behavior in the Built Environment

**Academic Editor:**Valery Sbitnev

#### Abstract

Unidirectional pedestrian movement is a special phenomenon in the evacuation process of large public buildings and urban environments at pedestrian scale. Several macroscopic models for collective behaviors have been built to predict pedestrian flow. However, current models do not explain the diffusion behavior in pedestrian crowd movement, which can be important in representing spatial-temporal crowd density differentiation in the movement process. This study builds a macroscopic model for describing crowd diffusion behavior and evaluating unidirectional pedestrian flow. The proposed model employs discretization of time and walking speed in geometric distribution to calculate downstream pedestrian crowd flow and analyze movement process based on upstream number of pedestrians and average walking speed. The simulated results are calibrated with video observation data in a baseball stadium to verify the model precision. Statistical results have verified that the proposed pedestrian diffusion model could accurately describe pedestrian macromovement behavior within the margin of error.

#### 1. Introduction

*Pedestrian Movement Theories*. Understanding pedestrian movement is essential for the safety planning and design of facilities. Pedestrian flow theories and models can support infrastructure designers and transport planners to optimize their plans and efficiently control pedestrian traffic and help architects facilitate circulation design in and around assembly occupancies, such as sports venues [1], transportation facilities and interchanges [2], and holy sites and political demonstrations [3]. To meet the increasing safety demands in organizing large scale public events, pedestrian flow theories are often introduced to improve the techniques and effectiveness of crowd management and control. In the design practice of pedestrian facilities, unidirectional movement such as the promotion of one-way traffic in urban road networks and separation of ingress and egress at airports, stations, subways, stadiums, exhibition centers, and supermarkets are considered to be beneficial for efficiently increasing pedestrian flow while maintaining the circulation safety. Unimpeded unidirectional flow keeps pedestrian crowds moving in steady state, which is also the most desirable circumstance in the evacuation process of large crowds.

Pedestrian movement theories are mostly planning related, with large urban and neighborhood scale supporting experimentations. Early studies on pedestrian behavior and movement were mostly empirical, and regression analysis has been widely used for proposing theories, understanding and evaluating design elements of pedestrian facilities, and developing planning guidelines [4]. Fruin proposed the level-of-service (LoS) concept and standard based on a range of pedestrian area occupancies [5]. Pauls studied pedestrian movement in buildings and design solutions of means of egress in cases involving both traditional total evacuation and selective sequential evacuation of large office buildings during fire drills for understanding the performance of the exit geometry, communication, operation, and training [6, 7]. Over the past twenty years, methods for describing and quantifying pedestrian movement attracted much attention from scholars around the world. A number of simulation models have been proposed and widely used for the prediction of pedestrian flow in the planning and design process, such as Simulex [8], Building Exodus [9], Legion [10], and Space Syntax [11]. Up until now, because of the inherent complexity of interpreting collective behavior, research efforts have continuously been evolving to improve the efficiency and accuracy of prediction models for pedestrian movement.

*State-of-the-Art Methods*. A number of computational methods focusing on the pedestrian movement have been introduced and simulation models have been developed in both academia and industry. Existing methods can be divided into two major categories: microscopic and macroscopic. Microscopic models simulate pedestrians as a set of individuals or agents and describe individualized behaviors and interactions in greater detail, which requires greater computational power. Macroscopic models, on the other hand, focus on overall pedestrian movement situations and view pedestrians as flows. These models can be used for minimum travel time estimation and significantly save computational time. In addition, there are some studies using the mesoscopic approach for pedestrian simulation in public buildings [12]. The idea is mainly based on grouping individuals. Physical groups can serve as input data, while virtual or logical groups are used to show similar behaviors or intentions [13].

The lattice-gas models, fluid-dynamic models, and network flow models commonly use macroscopic modeling approaches. Lattice gas models consider each pedestrian as an active particle on the grid for studying the features of pedestrian crowds by means of probability and statistics. Helbing et al. proposed a lattice gas model of pedestrian flows and compared the results with experimental data [14]. Fluid-dynamic models describe pedestrian crowds with fluid-like properties that behave similarly to gases or fluids. Helbing presented a model for pedestrian collective movement based on Boltzmann-like gas-kinetic approach [15]. Hughes presented a model based on three hypothesis from well-defined observations [16]. These fluid-dynamic models often describe density and velocity change over time using partial differential equations [17]. The network flow model can be also used to estimate pedestrian travel time and explore optimal routes. Hamacher et al. proposed a model which considered the quickest cluster flow problem based on network flow theories [18]. In addition, other common macroscopic models include the regression model, route choice model, queuing model, and gas-kinetics model.

Microscopic models can be divided into several categories: physical-based models, cellular automata models (CA), rule-based models, and multiagent models. Physical-based models recognize crowds as individuals reacting to surrounding events and have been used to study the emergency and panic situations, where optimal acceleration is controlled by varying physical force. The magnetic force model, NOMAD (Normative Pedestrian Behavior Theory), and social force model are examples of physical-based models [12]. CA models, also known as matrix-based models, divide a building floor into grids or cells. Each individual evacuee occupies one or multiple cells and moves forward when there are empty cells available. CA models can be further divided into either coarse-grid models or fine-grid models based on the size of the grids. Rule-based models have been widely used in computer graphics for modeling groups of animals and virtual crowds. An example of this is Reynolds’ “birds” model which followed three local rules [19]. These rules separation, alignment, and cohesion are used to simulate flocking behavior in the model. Multiagent models allow the environment to be populated with interactive and autonomous agents, which is a preferred approach for large scale and complex evacuation simulation.

*Current Barriers and Proposed Model*. Despite the numerous models that have been developed for simulating pedestrian movement, limitations still exist in interpreting collective behavior. Particularly in macroscopic methods, only a few studies integrate behavioral impact into the model. Hughes developed a continuum model to describe crowds as “thinking fluids” [16]. Colombo and Rosini proposed a model for pedestrian flow to describe typical features of flow under the effect of panic [20]. Helbing et al. investigated a model to simulate the dynamic feature of the disastrous form of collective human behavior caused by panic and jam [21]. The numerically efficient cellular automata (CA) models treat pedestrians as interacting particles. Some CA models allow partial behavior, which is a compromise between the macroscopic and microscopic approaches [22]. Diffusion phenomena have also been observed and specified in models for vehicle transportation [23, 24]. However, the existence of diffusion behavior in unidirectional pedestrian crowd movement has neither been noticed nor addressed yet in current models, which can be important in representing spatial-temporal crowd density differentiation in the movement process.

This paper proposes a new approach that incorporates diffusion behavior into modeling of pedestrian movement. The “spread-out” phenomenon that pedestrian crowds display during the unidirectional movement is illustrated in algorithms for pedestrian flow estimation. The proposed methodology features discrete-time, geometrical speed distribution and degree of diffusion in relation to speed divergence, which can be used to calculate pedestrian arrival state at any forward position and analyze unidirectional pedestrian walking processes from the macroscopic perspective. The model is then validated using video tracking data for the evacuation process of a long egress ramp in a stadium. The minimum difference value between the actual video sampling data and estimated results in various time intervals is used as an evaluation index to achieve the optimized degree of diffusion. The verification results show that the proposed model could accurately describe the pedestrian movement on the passage and makes the macroscopic pedestrian movement modeling approach more realistic by taking diffusion behavior into consideration.

#### 2. Methodology

##### 2.1. Diffusion Model

The diffusion modeling concept is illustrated in Figure 1. As shown in the figure, as a pedestrian crowd moves from left to right along a unidirectional long passage, both crowd density and flow rate are relatively higher at time segment . However, during the movement process towards the right direction, because of walking speed differentiations among individual pedestrians, some faster individuals in the group will overtake the slower ones. The interpersonal distances among each pedestrian will increase gradually, while pedestrian density per unit will decrease. Therefore, pedestrian flow rate at time () segment will be smaller than segment . This observed “spread-out” phenomena of the walkers can be described as crowd diffusion. The crowd diffusion phenomenon commonly exists in buildings and urban streets. Studying crowd diffusion behavior can be beneficial for analyzing pedestrian flow rate and density in buildings or on streets and provides a method for evaluating the rationality of level of service (LoS) and pedestrian facilities.

**(a) Pedestrian movement status at time**

**(b) Pedestrian movement status at time**

##### 2.2. Modeling Procedure

The proposed model partitions time into multiple intervals. The length of each time interval is seconds. Time 0 is defined by the time when the first pedestrian in the group walks through segment , hence represents the moment when the last pedestrian walk through segment , and the total number of time intervals is :where Ceiling is defined as a top integral function. For example, both* Ceiling (4.2)* and* Ceiling (4.7)* are equal to 5.

Then, the number of pedestrians that walk through segment is counted. Taking time interval () as an example, the number of pedestrians walking through segment is . (m/s) represents the average pedestrian walking speed. However, there are differences between the walking speeds of various pedestrians, assuming all pedestrian walking speeds follow geometric distribution.

For the pedestrians in one interval, their speeds can be calculated according towhere is the smoothing coefficient and is the maximum speed. The determination method for will be illustrated in the following paragraphs.

When the pedestrians arrive at the downstream segment , the arrival flow rate should follow geometric distribution as well. The predicted number of pedestrians that will arrive at segment at the th time interval equalswhere represents the walking time of the pedestrian with fastest speed from segment to segment , expressed by time intervals:where is the diffusion coefficient; is the travel time coefficient; and is the average walking time for all pedestrians from segment to segment , which equals , is the distance between segment and segment in meters.

is an important parameter in this model, as it can be used to depict the diffusion level of the pedestrians. It should be calibrated in practice by considering various contributing factors, such as the road width, road grade, pedestrian flow volume, and pedestrian types (e.g., male, female, teenagers, children, and the elderly). In this study, is affected by and . Different values of and will result in a different value. In Section 3.2, a calibration method is brought forward to determine the optimal values for and , which will also produce the optimal value for .

In (4) and (5), represents the degree of pedestrian diffusion, which is a positive number less than 1. Equation (4) shows that, as increases, the value decreases, meaning that the larger the diffusion coefficient, the lower the degree of pedestrian diffusion.

is a positive number smaller than 1. This parameter is introduced because there are differences in pedestrian walking speeds. Those who walk faster will arrive at segment earlier than other pedestrians. Parameter is the walking time of the fastest pedestrian from segment to segment . In other words, the pedestrian with the largest speed at segment during time interval will arrive at segment in the () time interval. In addition, the fastest walking speed is () times the average walking speed.

Thus, the maximum walking speed (m/s) equals

In (6), is the number of seconds in one unit interval and equals the minimum travel time between segments and .

To further develop the formula, (7) can be as follows:

Equation (7) estimates the number of pedestrians at each time interval in the downstream segment based on the number of pedestrians and the average pedestrian walking speed in each time interval from the upstream segment . Segment can be set up at any location downstream from segment . Therefore, (6) can be used to calculate pedestrian arrival state at any location in front of segment and to analyze the walking process of the pedestrian stream on the passage from macroperspectives.

##### 2.3. Model Features

The proposed crowd movement diffusion model has the following features.

###### 2.3.1. Discretization of Time Variable

In the modeling process, time is discretized into multiple intervals, the unit length of which is defined in seconds. Because of the individual differences in gait parameters and walking speeds, diffusion phenomena exists in each time interval. Diffusion behavior of the whole pedestrian stream consists of crowd diffusion behavior in each time interval. When compared to the continuous model, the discrete model is more easily translated into programming and the computation resource requirement is relatively low. This allows the model to satisfy the requirements of large scale software development.

###### 2.3.2. Walking Speed in Geometric Distribution

Pedestrian stream diffusion occurs because each pedestrian in the crowd has his/her own walking velocity. If all pedestrians walk at the same speed, then there will be no diffusion behavior. Therefore, the greater the individual pedestrian velocity differences are, the more serious the pedestrian stream diffusion can be, which means that as the value of increases, the value of decreases accordingly. This study describes the pedestrian crowd walking speed distribution by a geometrical distribution model. It cannot only describe the pedestrian walking speed differences, but also adjust the modeling structure by parameter to make the geometrical model more practical and realistic.

###### 2.3.3. Key Parameter Values and Impact

The degree of diffusion and walking time of the fastest person in each time interval are the key parameters. In the equations, represents the diffusion degree of the pedestrian stream, which is influenced by the divergence degree of each pedestrian’s walking speed. equals the walking time of the pedestrian with the fastest speed in each time interval from the upstream segment to downstream segment , represented by the number of time intervals. is a positive number less than or equal to 1.0. The higher the value is, the closer the fastest walking speed is to the average pedestrian walking speed; the lower the value is, the greater the fastest walking speed is compared to the average pedestrian walking speed. If is equal to 0.8, this indicates that the greater fastest walking speed is 1.25 (1.0/0.8) times larger than the average walking speed .

#### 3. Experimentation

##### 3.1. Data Collection

To verify the validity of the proposed model, we collected data of pedestrian movements on a field passage to verify the model precision. Since the values of parameters and have a large fluctuation range, this section will also verify the impact of these parameters on model precision.

Video recordings were conducted on the long egress ramp of Citizens Bank Park in Philadelphia (Pennsylvania, United States) during the spectators’ evacuation process. Citizens Bank Park is a baseball stadium built in 2004, the capacity of which is 43.500. The playing field is 7 m lower than the ground level. Spectators from the upper stands could evacuate to the ground by a ramp, escalators, and elevators. The observed match was very popular. All the seats had been occupied since the game started, and there were some spectators standing on the northwest open platform. The match (Phillies versus Yankees) started at 7:30 PM and finished at 10:15 PM on the 7th of August, 2010. The artificial illumination after the match could fully satisfy the spectators’ evacuation requirements. There were no difficulties in either way-finding or signage recognition. After the match, spectators immediately started to egress, and the vast majority of them chose to walk on the main ramp. The width and length of the ramp could facilitate many pedestrians and the continuous walking surface made walkers feel relaxed. Therefore, this main ramp was chosen as the observation and study site.

The evacuation process was video-recorded on site for 4 minutes, while a majority of spectators from the upper stands passed through the main ramp (Figure 2). We then analyzed video data in the laboratory to get individual flow data of each evacuee. Over 500 valid data plots on pedestrian movement and travel time were obtained. Crowd diffusion behavior on the ramp was also observed. When pedestrians walked towards the same direction, their desired speed was easier to achieve. The relaxed state encouraged greater differences in individual walking speed. In each time interval, pedestrians with higher walking speed showed overtake behavior. Also small group gathering in large crowds was noticed. Those family and/or friend groups tended to keep approximately the same walking speeds and closer personal distances, so that all group members could watch out for one another and stay all together along the way. These small groups of pedestrians often move at the speed of the slowest individual member in the group.

##### 3.2. Model Validation

Under the condition of and taking different values, the estimated results of flow rate in the downstream segment based on the diffusion model will be different as well. This study uses the minimum difference value between actual video sampling data on a number of pedestrians and estimated results in various time intervals through the downstream segment as the evaluation index to achieve the optimized and value as shown inwhere is the number of time intervals with pedestrian walk through downstream segment ; is the actual number of pedestrians at time interval in downstream segment ; is the predicted number of pedestrians at time interval in downstream segment .

Equation (8) shows that the prediction error indicator is the function of and . This study chooses and values which make achieve its minimum value as the best value of these two parameters. The value intervals of these two parameters are , starting from 0.1 and taking 0.1 as a time interval to determine the value of these two parameters. The set of value of each parameter is . This set includes 9 values in total. By using these 2 parameter values to develop verification plan, there are 81 verification plans in total.

##### 3.3. Result Analysis

As shown in Figure 2, upstream segment and downstream segment are chosen as study segments. The distance between these two segments is 100 m, m. The average pedestrian walking speed is equal to 1.4 m/s. This is slightly larger than the average pedestrian walking speed on roads with a plane surface, because this passage was a downward walking ramp. The phenomenon of pedestrians surpassing people in front was commonly recognized during data collection and individual pedestrian walking speed varied greatly. Therefore, the level of diffusion in the pedestrian stream was relatively higher on the long passage ramp.

When predicting the numbers of pedestrians at each time interval in the downstream segment , the value of is set as 5 seconds. The predicted results are shown in Figure 3.

As shown in Figure 3, the corresponding value of scenario 34 is the smallest, only 0.03. This indicates that the corresponding and value of test scenario 34 could achieve the preferable prediction effect. Moreover, the prediction error is only 0.03 at this point. This indicates that the difference between the predicted number of pedestrians at each time interval through downstream segment and the actual number of pedestrians was only 0.03 persons, which is a very small margin of error. Therefore, the proposed pedestrian diffusion model could accurately describe pedestrian macromovement behavior on the passage.

In scenario 3.4, is equal to 0.4, which indicates that the level of diffusion at this moment is not very high. is equal to 0.7, meaning that in each time interval the walking time of the fastest individual pedestrian is 0.7 times that of the average pedestrian walking time. Thus, his/her walking speed is equal to .

#### 4. Conclusion

Diffusion behavior in pedestrian flow commonly exists in unidirectional crowd movement such as during the evacuation process of sports venues, transportation facilities, large public exhibition buildings, and pedestrians walking on roads and pavements in urban and neighborhood scale. Quantifying and describing such diffusion behavior more precisely could allow for analysis of pedestrian movement speed and density effectively in the above facilities, thus having an important impact on evaluating the circulation performance and egress safety of the facility planning and design. Meanwhile, the proposed model can be used in pedestrian simulation software as well for simulating unidirectional pedestrian macroscopic movement in buildings and urban environments and improving the accuracy and applicability of current pedestrian simulation tools.

In the modeling process, this study is based on the hypothesis that all pedestrian walking speeds obey geometrical distribution and discretize time to build the pedestrian flow rate prediction model for the downstream segments. The proposed model has the following advantages: (1) the modeling structure is simple and could satisfy the requirements of real-time and rapidity; (2) the proposed model is a discrete model which is suitable for computer programming and is also easier to be embedded in pedestrian simulation software; and (3) obtaining a database does not require the movement trajectory of each pedestrian. The database only requires calculating the number of pedestrians in each time interval and the average pedestrian walking speed in the upstream segment, which could overcome the problem of complex structure in former macroscopic pedestrian simulation models.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

This work was financially supported by the National Natural Science Foundation of China (NSFC), Grants nos. 51308142 and 61304198, China Postdoctoral Science Foundation Funded Project (Grant no. 2014M550193), supported by “The Fundamental Research Funds for the Central Universities” (Grant no. HIT. NSRIF. 2014072).

#### References

- Y. Liu, D. Liu, N. Badler, and A. Malkawi, “Analysis of evacuation performance of merging points in stadiums based on crowd simulation,” in
*Proceedings of the 12th Conference of the International Building Performance Simulation Association (IBPSA, BS '11)*, pp. 2651–2658, 2011. View at: Google Scholar - S. B. Young, “Evaluation of pedestrian walking speeds in airport terminals,”
*Transportation Research Record*, vol. 1674, no. 1, pp. 20–26, 1999. View at: Publisher Site | Google Scholar - C. McPhail and R. T. Wohlstein, “Individual and collective behaviors within gatherings, demonstrations, and riots,”
*Annual Review of Sociology*, vol. 9, no. 1, pp. 579–600, 1983. View at: Publisher Site | Google Scholar - Transportation Research Board, “Highway capacity manual,”
*Special Report*204, TRB, Washington, DC, USA, 1985. View at: Google Scholar - J. J. Fruin,
*Pedestrian Planning and Design*, no. 206, 1971. - J. L. Pauls and B. K. Jones, “Building evacuation: research methods and case studies,” in
*Fires and Human Behavior*, pp. 251–275, 1980. View at: Google Scholar - J. Pauls, “The movement of people in buildings and design solutions for means of egress,”
*Fire Technology*, vol. 20, no. 1, pp. 27–47, 1984. View at: Publisher Site | Google Scholar - P. A. Thompson and E. W. Marchant, “Testing and application of the computer model ‘SIMULEX’,”
*Fire Safety Journal*, vol. 24, no. 2, pp. 149–166, 1995. View at: Publisher Site | Google Scholar - E. R. Galea and J. M. Perez Galparsoro, “A computer-based simulation model for the prediction of evacuation from mass-transport vehicles,”
*Fire Safety Journal*, vol. 22, no. 4, pp. 341–366, 1994. View at: Publisher Site | Google Scholar - G. K. Still,
*Crowd dynamics [Ph.D. thesis]*, University of Warwick, Coventry, UK, 2000. - A. Penn and A. Turner, “Space syntax based agent simulation,” in
*Proceedings of the 1st International Conference on Pedestrian and Evacuation Dynamics*, Essen, Germany, 2001. View at: Google Scholar - N. Shiwakoti, M. Sarvi, and G. Rose, “Modelling pedestrian behaviour under emergency conditions—state-of-the-art and future directions,” in
*Proceedings of the 31st Australasian Transport Research Forum (ATRF '08)*, pp. 457–473, 2008. View at: Google Scholar - J. Tolujew and F. Alcalá, “A mesoscopic approach to modeling and simulation of pedestrian traffic flows,” in
*Proceedings of the 18th European Simulation Multi-Conference*, pp. 123–128, 2004. View at: Google Scholar - D. Helbing, M. Isobe, T. Nagatani, and K. Takimoto, “Lattice gas simulation of experimentally studied evacuation dynamics,”
*Physical Review E*, vol. 67, no. 6, Article ID 067101, 2003. View at: Google Scholar - D. Helbing, “A fluid-dynamic model for the movement of pedestrians,”
*Complex Systems*, vol. 6, no. 5, pp. 391–415, 1992. View at: Google Scholar | Zentralblatt MATH | MathSciNet - R. L. Hughes, “The flow of human crowds,”
*Annual Review of Fluid Mechanics*, vol. 35, pp. 169–182, 2003. View at: Publisher Site | Google Scholar - X. Zheng, T. Zhong, and M. Liu, “Modeling crowd evacuation of a building based on seven methodological approaches,”
*Building and Environment*, vol. 44, no. 3, pp. 437–445, 2009. View at: Publisher Site | Google Scholar - H. W. Hamacher, K. Leiner, and S. Ruzika, “Quickest cluster flow problems,” in
*Pedestrian and Evacuation Dynamics*, pp. 327–336, Springer, New York, NY, USA, 2011. View at: Publisher Site | Google Scholar - C. W. Reynolds, “Flocks, herds and schools: a distributed behavioral model,” in
*Proceedings of the 14th Annual Conference on Computer Graphics and Interactive Techniques (SIGGRAPH '87)*, pp. 25–34, ACM, 1987. View at: Publisher Site | Google Scholar - R. M. Colombo and M. D. Rosini, “Pedestrian flows and non-classical shocks,”
*Mathematical Methods in the Applied Sciences*, vol. 28, no. 13, pp. 1553–1567, 2005. View at: Publisher Site | Google Scholar | MathSciNet - D. Helbing, I. Farkas, and T. Vicsek, “Simulating dynamical features of escape panic,”
*Nature*, vol. 407, no. 6803, pp. 487–490, 2000. View at: Publisher Site | Google Scholar - G. J. Perez, G. Tapang, M. Lim, and C. Saloma, “Streaming, disruptive interference and power-law behavior in the exit dynamics of confined pedestrians,”
*Physica A: Statistical Mechanics and its Applications*, vol. 312, no. 3-4, pp. 609–618, 2002. View at: Publisher Site | Google Scholar - Y. Bie, Z. Liu, D. Ma, and D. Wang, “Calibration of platoon dispersion parameter considering the impact of the number of lanes,”
*Journal of Transportation Engineering*, vol. 139, no. 2, pp. 200–207, 2013. View at: Publisher Site | Google Scholar - Y. Bie, D. Wang, and X. Qu, “Modelling correlation degree between two adjacent signalised intersections for dynamic subarea partition,”
*IET Intelligent Transport Systems*, vol. 7, no. 1, pp. 28–35, 2013. View at: Publisher Site | Google Scholar

#### Copyright

Copyright © 2015 Ying Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.