Abstract

This paper pertains to the study of finite-time control of one dimensional crowd evacuation system. Benefiting from the research of fluid dynamics and vehicle traffic, a one dimensional crowd evacuation system is constructed, whose density-velocity relationship is represented by a diffusion model. In order to deal with the nondirectionality of crowd movement, the free flow speed is chosen as a control variable. Since the control variable is included in a partial derivative, it increases the difficulty of designing the controller. In this paper, finite-time controller is designed, which not only guarantees the effective evacuation, but also obtains the estimation of evacuation time. Then, finite-time tracking problem is solved, which makes the density converge to a given density. Finally, numerical examples illustrate the effectiveness of the controllers.

1. Introduction

In everyday life, crowds would gather in many places, for example, subway stations, stadiums, and cinemas. Effective measures should be taken to ensure the safety and comfort of pedestrians, so how to evacuate people when emergencies occur has been a challenging job. In the early 1990s, the International Conference on Engineering for Crowd Safety [1] has shown the importance of this topic. Physicists, sociologists, psychologists, computer scientists, and traffic scientists have been conducting depth-going researches on this topic from their research fields.

Modeling the crowd dynamics is the primary task of crowd evacuation, but due to the complexity and uncertainty of crowd dynamics, it is very difficult to build a model that suits all situations. Therefore, various models have been developed, such as the continuum model, the network-based model [2], agent-based models [3, 4], game-theoretic models [5], cellular automata models [6, 7], and the fractional model [8].

In this paper, the continuum model, a macroscopic simulation model, is recommended. At medium and high densities, the motion of pedestrian crowds shows some striking analogies with the motion of fluids [9], so the theory describing fluid dynamics is introduced to describe pedestrian dynamics. Pedestrians are treated as a collection rather than individuals. Average density and velocity at a given location are proposed to describe the crowd dynamics. Based on three hypotheses, a first-order continuum model was developed to describe the pedestrian dynamics in [10]. Huang et al. [11] provided an efficient method to solve Hughes’ model. Appert-Rolland et al. [12] extended a macroscopic vehicle traffic model to pedestrian traffic. The cell transmission approach and the continuum approach were combined in [13] to predict densities and travel times. Some continuum models were compared by making use of numerical method in [14]. All above focus on modeling the pedestrian dynamics in different situations, but few literatures present strategies for controlling the crowd dynamics. Wadoo [15] designed advective, diffusive, and advective-diffusive controllers for crowd dynamics in one dimension. Sliding mode control method was applied to crowd dynamic models for the synthesis of robust controllers in [16]. Dong et al. [17] designed feedback control law for two dimensional crowd model. Robin, Neumann, and Dirichlet boundary control laws are designed for a disturbed crowd evacuation system in [18].

In this paper, the problem of finite-time evacuation is studied. In some special circumstances, such as fires and terrorist attacks, pedestrians need to be evacuated as soon as possible. Evacuation strategy cannot be used until its effectiveness is tested, because the cost involves not only property damage but also pedestrian injury or death. Therefore, it is particularly important to estimate the evacuation time. The estimated evacuation time can be used to evaluate the effectiveness of the evacuation strategy, and the corresponding evacuation strategy can be adjusted to minimizing the loss. Benefiting from Orlov’s research [19], finite-time controller is designed to evacuate pedestrians in finite time, and then the finite-time control method is extended to deal with the tracking problem which makes the crowd density follow a given reference density. Some of the latest research on advanced control, such as adaptive control [2022], guaranteed cost control [23], control [24], and tracking control [25, 26], gave us inspiration in the design of the controller. The control and stability problem are formulated directly in the framework of a distributed model of partial differential equations (PDEs), which can avoid errors introduced by spatial discretization. The development of hardware technology has made distributed sensors and actuators realistic, but we mainly explore the crowd management strategies on theoretical interest in this paper, and the implementation needs further study in the future.

The rest of this paper is organized as follows. One dimensional crowd evacuation dynamic is modeled in Section 2. Section 3 presents a feedback controller to evacuate pedestrians in finite time. Section 4 designs a finite-time controller to make the crowd density profile track a reference density. Simulation examples are given in Section 5 to illustrate the effectiveness of the controllers. Conclusions and future work are discussed in Section 6.

Notation. The notation is used throughout the paper. denotes the infinite-dimensional Hilbert space on interval , represents the interval length, with normFor notational convenience, we denote

2. Mathematical Modeling

The Lighthill-Whitham-Richards (LWR) model [27, 28] based on the conservation law of mass is recommended in this paper to represent the crowd evacuation dynamics in one dimension, implying that the number of pedestrians coming in and going out of a corridor section account for the change of crowd density on that section. The LWR model is given bywhere , and are space and time variables, respectively. is the average crowd density and is the partial derivative of the density with respect to time at position . denotes the flux of the crowd and is the partial derivative of the flux with respect to position at time . The flux is a function of and the average crowd speed , as shown below,

Various models have been developed to mimic the velocity-density relationship, such as Greenshield model, Drew model, Greenberg model, Pipes Munjal model, and Underwood model [29]. Here, the diffusion model is recommended, which is given aswhere denotes the free flow speed, that is, the maximum moving speed when the density is zero. is the maximum crowd density, and denotes diffusion coefficient which is a positive constant, given by , where is a random velocity and is a relaxation parameter.

Remark 1. The diffusion model (5) is an extension of the Greenshield’s model where the speed depends not only on the traffic density but also on the density gradient. The diffusion term demonstrates the fact that pedestrians can adjust their movement speed in real time based on the density ahead. The adjustment makes the change of their speed gradual rather than abrupt in response to the shock wave, which creates the dependence of speed on density gradient [29].

Combining the LWR model (3) with the diffusion model (5) yieldsBy choosing the free flow speed as the distributed control variable denoted by u(x,t), one can derivewhere is the controller and is the maximum velocity. The model is subject to the following initial condition and boundary condition

Remark 2. The paper [30] has suggested that pedestrian traffic can be handled in the similar way as the vehicle traffic. But there is a main difference between pedestrian traffic and vehicle traffic. In the vehicle traffic, the car on a lane is unidirectional, so its speed can be fixed by the traffic density, using the diffusion model. While in the pedestrian traffic people can move in both directions, its density cannot fix the speed with any velocity-density relationship model, so the free flow speed is chosen as the control variable, and, with the actuation system, people can be told to change their speed.

3. Finite-Time Control

In this section, by virtue of the finite-time control theory, a distributed controller is designed to make the state converge to zero in finite time. The stability of the crowd evacuation system under the finite-time controller is analyzed using the Lyapunov method.

The important finite-time control theory ([19], Lemma 4.3) is stated by the following lemma.

Lemma 3 (see [19]). Let an everywhere nonnegative function meet the differential inequalityfor all and for some constants and . Then, for all

The following lemma is an important inequality used in our proof, which can be considered as a special case of Hölder integral inequality.

Lemma 4 (see [31]). Consider an arbitrary real coefficient , and let , where is p-th integrable Banach space defined on interval . Then the following inequality holds,

In order to stabilize the crowd evacuation system (7) in finite time, the following distributed controller is designed,where is a constant control coefficient, is a constant.

Remark 5. The distributed control mentioned in [32, 33] means that the controller of a system can use information of the connected systems to construct the control strategy. However, the distributed controller mentioned in this paper is the controller of distributed parameter systems (as opposed to a lumped parameter system) whose state space is infinite-dimensional.

Theorem 6. The crowd evacuation system (7), subject to the initial and boundary conditions (8), with the distributed controller (12), achieves the attainment of in norm when .

Proof. Consider the Lyapunov functionalComputing the time derivative of for yieldsSubstituting the controller (12) into (14) for and using the Leibniz integral ruleone can deriveIntegrating by parts the first term of (16) and considering the boundary condition (8),As to the second term of (16), by using the same manipulations, one can derive that is,Substituting (17), (19) into (16), the time derivative of for becomesAs to the second term of (20), by virtue of Lemma 4, equation (20) can be rewritten in the formAccording to Lemma 3, when , converges to zero, that is, the distributed controller (12) makes the crowd evacuation system (7)-(8) achieve the attainment of in norm.

4. Finite-Time Tracking Control

In this section, a finite-time tracking controller is designed to make the crowd density follow a given reference density .

Assumption 7. The reference density is smooth enough and its spatial derivatives up to the second order are square integrable in norm. Also, it satisfies the initial and boundary conditions

Define the tracking error as . By using (7), the error dynamic is given asand the initial and boundary conditions are

In order to stabilize the error dynamic systems (23) in finite time, the distributed controller is designed aswhere is a constant control coefficient and is a constant. Then, the following result is gotten.

Theorem 8. Consider the crowd evacuation system (7) with the initial and boundary conditions (8) and the reference density satisfying Assumption 7. Then, the error dynamic systems (23) subject to the boundary conditions (24) can be stabilized to zero in norm with the distributed controller (25) when .

Proof. Consider the Lyapunov functional Computing the time derivative of for yieldsSubstituting the distributed controller (25) into (27) for , and using the Leibniz integral rule (15), one can deriveAs to the first term of (28), integrating it by parts and considering the boundary condition (24), the equation (28) can be written as By virtue of Lemma 4, the time derivative of for is According to Lemma 3, the error dynamic system (23) with boundary conditions (24) can be stabilized to zero when .

Remark 9. In the application of crowd evacuation, the crowd density can be stabilized to different values to achieve different control objectives, such as maximizing the evacuation flow and maximizing the pedestrian movement speed. Therefore, the research of tracking control is of great practical significance. Meanwhile, estimating the time when the crowd density stabilizes to the reference density can evaluate the effectiveness of the control strategy.

5. Simulation Results

In this section, numerical results are given to illustrate the effectiveness of the finite-time controller (12) and the finite-time tracking controller (25), respectively. The numerical method is the finite volume method. For simulation, the initial density is given by , with being the center of the Gaussian distribution and being the highest magnitude of the distribution. The constant reference density is and a general reference density is given by . The main parameters are .

When the controller , the mathematical model of the crowd evacuation system (7) shown in Figure 1 is a diffusion model .

Figure 2 illustrates the density response of the crowd evacuation system with the finite-time controller (12). Because of the effect of the advection term, the density profile moves towards the exit () and converges to zero in finite time. For a clearer demonstration, the density evolutions at and are shown in Figure 3. The density at becomes zero at about 6.1 seconds, the density at becomes zero at about 7.1 seconds, and the crowd evacuation process ends. The evacuation time calculated by Theorem 6 is It can be easily seen that the calculated evacuation time and the simulated evacuation time are almost equal, and the error is within a reasonable range.

Figure 4 demonstrates the density response of the crowd evacuation system with the finite-time tracking controller (25), where is chosen as the reference density. The density evolutions of and are shown in Figure 5. It can be seen that the density at converges to 2.5 at about 1.9 seconds and the density at converges to 2.5 at about 2.2 seconds. The evacuation time calculated by Theorem 8 isThere is a small error between the calculated evacuation time and the simulated evacuation time, which may be caused by the discretization of simulation, but it is within a reasonable error range.

Next, a more general reference density is selected to show the effectiveness of the finite-time tracking controller, as shown in Figure 6.

The density response of the crowd dynamic system under the finite-time tracking controller is demonstrated in Figure 7. Figure 8 illustrates the density evolution at and with reference . The density at reaches the reference density at about 2.0 seconds and the density at reaches the reference density at about 2.9 seconds. Comparing the evacuation time calculated by Theorem 8, and they are almost equal.

To sum up, the effectiveness of the designed controllers has been shown by the comparisons, and the calculated evacuation time is almost equal to the simulated evacuation time, so the estimated evacuation time mentioned in the theorem is feasible.

6. Conclusion

In this paper, the crowd dynamic model was constructed by combining the LWR model and the diffusion model. Then, finite-time controllers were designed for the crowd evacuation system, which solved the problem of nondirectionality of crowd movement and got the estimation of evacuation time. This theoretical research can promote the improvement of practical application, but the effect of time delay and disturbance in implementation needs to be further studied.

Data Availability

All the data used to support the findings of this study are available from the corresponding author upon request. The email address of the corresponding author is [email protected].

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by National Natural Science Foundation of China (grant number 61807016) and Postgraduate Research & Practice Innovation Program of Jiangsu Province (grant number KYCX).