Abstract

In this paper, a lattice hydrodynamic model of two-dimensional bidirectional pedestrian traffic is proposed with consideration of altering path under the cumulative effect of delay time. The stability condition is acquired by linear analysis, and the mKdV equation to describe congestion evolution is derived by reductive perturbation technique. According to the result from the stability analysis, the stability region of pedestrian flow can be divided into stable region, unstable region, and metastable region. On the basis of stable condition, the unstable region is narrowed with the increase of delay time td and the path changing rate γ. It indicates that changing path can effectively improve the stability of the pedestrian flow under the cumulative effect of delay time td. For numerical simulation and analysis of density wave, it is found that the increase of path changing rate γ and the cumulative effect of delay time td are conducive to alleviate pedestrian congestion.

1. Introduction

With the increase of urbanization process, the study of city traffic system has become the hot spot which the scientific and social pays attention. As an important part of the city traffic system, pedestrian transportation system brings great convenience for people to travel. But in recent years, the pedestrian stampedes events happened frequently, which has brought serious disasters to the society and brought pain to the victims’ families. On December 31, 2014, 36 people were killed and 49 injured in a stampede at Chen Yi Square on the Bund in Shanghai [1]. Saudi authorities issue the Hajj stampede in Mecca on September 26, 2015, have killed 769 people, and wounded 934 others who are from different countries [2]. The stampede caused 22 deaths and 27 injured in Mumbai, India, on September 29, 2017 [3]. How to keep evacuate people safely and effectively in limited space has attracted major concern.

In order to better solve the problem of pedestrian traffic congestion, researchers have made a lot of efforts. From about 1980 onwards, research on pedestrian traffic was often divided into two branches: pedestrian simulation and pedestrian model. Pedestrian simulation focuses on reproducing the natural behavior of pedestrian flow. And with the development of computer graphics, virtual pedestrians can be realized.

There are many similar behaviors and different characteristics between the vehicular traffic and the pedestrian traffic. First of all, on the macro level, pedestrian traffic presents the same flow characteristics as vehicle traffic. Similarly, the average velocity and flux can be introduced to describe flow characteristics. However, pedestrian traffic composed of the different individuals shows stochastic properties and multidirection movement due to fewer constraints. Thus, pedestrian traffic is more complex and flexible than vehicular traffic. Henderson firstly put forward the macroscopic model of pedestrian flow in 1971 by applying gas dynamics model and fluid dynamics model to the crowd experience data. He argued that the crowd’s behavior is similar to gas under the freely moving phase with low density and similar to fluids under congested phase with high density. However, the energy conservation in the interaction of pedestrians cannot be satisfied in case of pedestrian conflict [4]. In 1992, Helbing proposed a better fluid dynamics model on the basis of Henderson’s model by taking into account the intention, the desired velocity, and interaction between pedestrians [5]. Hughes’ model of crowd gives the following three assumptions: (1) pedestrian velocity is only determined by the pedestrian density; (2) pedestrian has about the destination of the common sense; and (3) pedestrians seek the shortest travel time [6]. Existing pedestrian models mainly include cellular automata model, lattice gas model, agent-based model, social force model, and hydrodynamic model [4, 720]. Helbing and Molnár presented the social force model to describe pedestrian dynamics in 1995 and assumed that the interaction between pedestrians and pedestrians and between pedestrians and environment is expressed by social forces [21]. But it will cause unexpected effects when the superposition principle of force is satisfied under high density. Cellular automata model (CA model) of pedestrian traffic flow is a micromodel which discretizes time and space. Introducing the cellular automata model of traffic flow into pedestrian flow, it is necessary to refine the discrete space. Each pedestrian to occupy the multiple cells selects the moving direction and speed according to the simple rules [22]. At present, these models have been used for pedestrian facility design, calculation of channel capacity, and emergency evacuation. And many characteristics of pedestrian traffic can be extensively applied to many occasions such as bidirection pedestrian, fire evacuation, aircraft evacuation, T-shaped intersection, bottleneck road section, and slope evacuation [2327].

Comparing with micro-CA model, the macroscopic model is suitable to real-time online simulate pedestrian flow and carry out prediction due to the neglect of pedestrian movement details resulting in rapid calculation. The first microlattice hydrodynamic model was put forward by Nagatani in 1998. He simplified Kerner and Konhäuser's high-order continuous model and referenced to the idea of the car-following model obtained by flow dynamics equation and conservation equation [28]. And for simulating more actual status of traffic flow, Nagatani proposed a two-lane traffic model to consider the problem of vehicle lane change in 1999 [29]. When the motion of a crowd is seen as fluid flow at a macrolevel, the lattice dynamic equations describing the flow can be applied to pedestrian traffic. Based on Nagatani’s lattice hydrodynamic model, Tian et al. extended the two-lane lattice hydrodynamic model of traffic flow to the two-dimensional lattice hydrodynamic model considering path change in bidirectional pedestrian traffic flow and studied the impact of path changing rate on the stability of pedestrian traffic flow through discussing the value of critical point ac [30, 31]. Kuang et al. presented a new lattice hydrodynamic model for two-way pedestrian flow considering pedestrian visual field effect. The phase transitions between the freely moving phase, the coexistence phase, and the jamming phase below which the critical point ac has been found [32]. Zhou et al. studied the effect of memory effect on pedestrian flow and assumed pedestrians may walk by his memory when it is difficult for pedestrians to accurately capture the position [33]. Li et al. extended the original hydrodynamic model of traffic flow to single row pedestrian movement with medium and high density and considered asymmetric interaction [34]. Yu et al. solve the general solution of the mKdV (modified Korteweg–de Vries) equation [35]. In 2015, Redhu and Gupta present a new lattice hydrodynamic model based on delayed feedback control considering the flow difference effect [36]. In 2016, Redhu and Gupta propose a multiphase lattice hydrodynamic model by considering the forward looking effect and discuss the impact of forward leading lattices on multiphase traffic dynamics [37]. Peng et al. consider the lateral effects and extended the lattice hydrodynamic model to the non-lane-based lattice hydrodynamic model [38]. Wang et al. study the impact of the multiple density difference effect on the stability of traffic [39].

Until now, the lattice hydrodynamic model of pedestrian flow with consideration of path change has been little developed. Also, the cumulative effect of delay time is important to adjust current pedestrian flow reaching the optimal status. The movement of pedestrians will be affected by the transverse effect and the longitudinal effect between them. Assumed that if there is more space to allocate pedestrian in the transverse direction based on the optimal effect and comfort level, one can select to continue along the way forward or also select to change path with a faster rate. In the longitudinal direction, Tian et al. have studied that the current grid point flow is adjusted by the optimized flow of the grid point ahead. The pedestrian flow reaches the optimal flow within the delay time td, while the downstream optimal flow can continuously adjust the current flow under the cumulative action of the delay time. Considering both the transverse and longitudinal effects of pedestrian traffic, the analysis of pedestrian traffic evolution characteristics will be more real.

In this paper, a lattice hydrodynamic model of change path pedestrian traffic under the cumulative action of delay time is presented in Section 2. Section 3 carries out the linear stability analysis to obtain the stability condition. In Section 4, the mKdV equation of congestion evolution is derived by applying perturbation theory to perform nonlinear analysis. Section 5 executes the numerical simulation to analyze the evolution of the density wave to discuss the impact of the cumulative effect of delay time. Finally, the main conclusion is simply summarized.

2. Models

In 1999, Nagatani extended the one-lane lattice model to two-lane traffic lattice model and assumed the changing lane happens when the density at site on the second lane is greater than the density at site j on the first lane, and the changing lane rate is , where and , respectively, denote the local density at site j on the second lane and the local density at site j on the first lane. The changing path frequency is denoted as γ, which reflects the strength of the changing path. The total average density is and the corresponding first derivative is . The velocity is called the optimal velocity function of the density , which is defined as follows [28]:where denotes the maximal velocity of traffic flow and represents the critical density which is the inflection point of the function of optimization velocity . Similarly, if the density at site j on the first lane is greater than that at site on the second lane, the changing lane rate from the first lane to the second lane is .

The conservation equation of one-dimensional lattice hydrodynamic model is expressed as follows:where .

We will extend equation (2) to two-dimensional pedestrian flow with consideration of pedestrian changing path. Pedestrians usually change travelling path along the transverse direction in Public Square and so on. As shown in Figure 1, four-way pedestrian traffic is considered: eastbound, westbound, northbound, and southbound pedestrian. Assumed the fraction of pedestrians moving along each direction is different. The ratio of eastbound and westbound pedestrians to total pedestrians is C, in which the ratio of westbound pedestrians to eastbound is C1 and the westbound pedestrians correspond to . The ratio of northbound and southbound pedestrians to total pedestrians is . The ratio of southbound pedestrians to southbound is C2 and the northbound pedestrian is . The eastbound, westbound, northbound, and southbound pedestrian local densities at site (j, m) at time t are represented by , respectively. The average density can be obtained as follows:

The conservation equation of the two-dimensional bidirectional pedestrian traffic flow can be written as follows:where . The eastbound, westbound, northbound, and southbound pedestrian flows are denoted as , respectively.

To simulate more realistic pedestrian traffic conditions with consideration of the impact of downstream pedestrian flow on current pedestrians, the evolution equation of one-dimensional pedestrian traffic flow is obtained as follows:where is the pedestrian flow at site j at time t and a is the sensitivity coefficient of pedestrian flow (the inverse of the response time τ). In effect, the response time is a relaxation time when the pedestrian flow reaches a steady state from a nonequilibrium state, which corresponds to the time interval in response to a stimulus. λ denotes the strength coefficient, and the integral term represents the cumulative effect of downstream optimization flow and current pedestrian flow differences within delay time td. The delay time td is the time pedestrian needs to make decisions about what to travel next and its value is larger than the response time. We employ the following trapezoidal integral formula to approximately estimate integral in equation (5) and extend to four-way pedestrian flow:

The evolution equation of four-way pedestrians flow can be written as follows:

By eliminating velocity in equations (4) and (7), the density evolution equation of four-way pedestrian flow is derived as follows:

3. Linear Stability Analysis

To obtain the stability condition of equation (8), the linear stability analysis of the lattice hydrodynamic model is performed as follows. At first, we consider the uniform pedestrian flow which has constant density and optimization velocity . The solution of the uniform steady state is obtained as . Under disturbance , the density of pedestrian flow is expressed as . Substituting it into equation (8), the linearized equation is derived as follows:

Letting the disturbance be expressed as a Fourier series and substituting it into equation (9), we can gain the following equation:

By performing Taylor expansion for z and keeping the first- and second-order terms of (ik), we can abstract the following equation from equation (10):where . The uniform steady-state flow becomes unstable pedestrian flow for long-wavelength modes in the condition if z2 is a negative value. When z2 is a positive value, the state of the steady flow maintains. The neutral stability condition can be obtained as follows:

When γ = 0.0, the following stability condition is obtained:

The stability condition is identical with the solution of the lattice hydrodynamic model of bidirectional pedestrian traffic presented by Cen et al. [40]. By comparing the stability condition, we discovered that the critical value of the sensitivity coefficient without changing path frequency is larger than the critical value of the sensitivity coefficient considering changing path frequency. It shows the importance of considering changing lane to alleviate traffic congestion. Moreover, if the cumulative effect of delay time td is not taken into account, the stability condition equation (12a) become as follows:

It is the same as the research results in [31]. Thus, it indicates the cumulative effect of delay time td is conducive to suppress traffic jamming.

4. Nonlinear Analysis

Using the reductive perturbation technique to equation (8), we carry out nonlinear analysis by considering long-wave patterns in pedestrian flow on coarse-grained scales. In order to describe the long-wave modes, the slowly varying behavior of long wavelengths near the critical point () is studied. Introducing a small parameter ε (), we defined the slow space variables X and slow time variable T:where b is an undetermined constant. The instantaneous density is obtained as follows:using equations (13) and (14) and expanding each term of equation (8) to the fifth-order of ε, we gainwhere the first derivative and third derivative of the optimal velocity at the critical point () are . By taking and introducing , the second-order term and the third-order term of ε are eliminated. We derived the following equation:where

For obtaining the high-order regularized equation, the transformations are introduced:

Consequently, the mKdV equation with term is obtained as follows:where . The solution of the mKdV equation (19) is given bywhere and .

Thus, the solution of kink-antikink wave of the mKdV equation is yielded as follows:where , and the amplitude A of kink-antikink wave is expressed as follows:

The solution of kink-antikink wave indicates the coexisting phase which is composed of the free flow phase (low density) and the congested phase (high density). The density of the free flow phase is, and the density of the congested phase is . Thereafter, the coexisting curve can be shown as the dotted curve in Figure 2.

Figure 2 reveals the neutral stability lines with solid lines and the coexisting curves with dotted curves for a fixed value of the ratio C1, C2, and C, the delay time td, the cumulative coefficient λ, and the changing path rate γ. The phase diagram can be divided into three regions: the stable region above the neutral line, the metastable region between the neutral stability line and the coexisting curve, and the unstable region below the neutral line. In the metastable region, if the disturbance strength exceeds its critical amplitude, the disturbance will become larger with the time. The pedestrian flow is in the unsteady state, which causes pedestrian congestion. When the disturbance is less than the critical amplitude, pedestrian flow is still in the steady state. In the stable region, pedestrian flow remains the steady state due to the attenuation of disturbance. On the contrary, the disturbance in the unstable region gradually increases with time and eventually develops into kink-antikink wave.

Figure 2(a) plots that the neutral stability lines obviously reduce with increasing value of the ratio C. It is shown that the distribution of pedestrian density impacts the stability of the pedestrian traffic system, and reasonable distribution of pedestrian density can improve the stability of the traffic system. By comparison of Figures 2(a) and 2(b), the increase of delay time td leads to the contraction of the unstable region. It indicates that the delay time td can significantly enhance the pedestrian traffic stability. Moreover, the stable region is accordingly expanded with the increase of the changing path frequency γ. It means the increase of the changing path frequency γ has a favorable influence on suppressing pedestrian traffic congestion.

5. Simulation and Analysis

In order to verify the results in the stability analysis and nonlinear analysis, we carry out simulation and comparison with the consequence of the theoretic analysis. The periodic boundary condition is adopted in the process of the numerical simulation. The initial density and the critical density are chosen as and the sensitivity of driver a = 1.8. The four-way traffic on square is divided into L × L lattice with L = 100. The small perturbations are added to a two-dimensional four-way pedestrian traffic: and . The unit time step is equal to the inverse of the sensitivity a.

The evolution pattern of density wave is depicted in Figure 3 at 2000 time steps for C1 = C2 = 0.1 without delay time td. As we can see that the fluctuation of density wave in Figure 3(a) is great without changing path frequency. It is easy to cause congestion and is not conducive to evacuation. In Figures 3(b) and 3(c), the fluctuation of density wave is gradually weakened under the increase of path change frequency γ, it indicates that changing path can improve the stability of the traffic system on certain level, and it is consistent with the stability condition equation (12c). When the ratio C is less than or equal to 0.5, Figures 3(b) and 3(d) exhibit that there is no congestion of pedestrian traffic with the increase of the ratio C, which is consistent with the theoretical analysis of system stability.

Figure 4 shows the spatial distribution of density at 2000 time steps for C1 = C2 = C = 0.1. It is clearly observed that the cumulative effect of delay time td can significantly improve the pedestrian traffic stability and eliminate jams of pedestrian traffic from the fluctuation of density wave in Figures 3 and 4. If pedestrians do not think about changing path (γ = 0), Figure 4(a) shows the great changes in crowd density. When the cumulative effect of delay time td enhances, the density change in Figure 4(b) is very small. It indicates the more impact of the greater strength coefficient λ on suppression of traffic congestion. The cumulative term of delay time in equation (5) is suitable to use as a feedback control one.

6. Conclusions

In this paper, a two-dimensional lattice hydrodynamic model of pedestrian traffic flow is proposed considering path change under the cumulative effect of delay time. Using the linear stability analysis, the neutral stability condition is obtained. By performing the nonlinear analysis, the mKdV equation describing the density wave of pedestrian flow is derived and the coexistence curve is gained. The phase diagram is divided into three regions: stable region, unstable region, and metastable region. The theoretical analysis shows that changing path γ and the cumulative effect of delay time td can reduce the value of critical point ac and significantly improve the stability of pedestrian flow. The unstable region can be shrunken by the increase of path change frequency, which may be conducive to crowd evacuation. Numerical simulations further illustrate the improvement of changing path rate and delay time td on pedestrian traffic flow. Moreover, increasing ratio C (C < 0.5) of eastbound pedestrians and westbound pedestrians to total pedestrians can also improve the stability of pedestrian traffic flow. It indicates pedestrian density must be properly distributed in order to improve the stability of the traffic system.

Data Availability

The data used to support the findings of this study are included within the article’s references.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This project was supported by the National Natural Science Foundation of China (Grant nos. 11962002 and 11672176), the Natural Science Foundation of Guangxi, China (Grant no. 2018GXNSFAA138205), and the Doctoral Research Project of Guangxi University of Finance and Economics (Grant no. BS2019027).