# A Social Force Model for Pedestrians’ Movements Affected by Smoke Spreading

Authors: J. Makmul
Journal: Modelling and Simulation in Engineering
 (2020)

---

## Abstract

A social force (SF) model is proposed to simulate the egress of pedestrians while smoke is spreading. The advection-diffusion with source term is used to describe the propagation of smoke. It is incorporated into the SF model. The navigation field in our model is determined by the negative gradient of the solution of the Eikonal equation. It depends on the pedestrian and smoke density. Numerical experiments are performed in a room with multiple exits, and their results are shown.


---

## Body

## 1. Introduction

Many crowd disasters occurred all around the world. This leads to a large number of losses in human life and wealth. Due to these losses, many researchers study pedestrians’ movements and their behaviors in a crowd, especially in evacuation situations. The development of pedestrian dynamic models can help to plan public buildings with respect to safety and evacuation. Nowadays, there are numerous pedestrian models, such as the particle flows [1], the social force model [2], the optimal-velocity model [3], the magnetic force model [4], the cellular automata models [5], and the discrete choice model [6].There is quite a large number of researches focusing on modeling pedestrians’ movements during smoke or fire spreading. Zheng et al. [7] proposed an extended floor-field model to study the influence of fire and smoke on the movement of pedestrians. The results of their simulations demonstrate that the fire location, the spreading rates of the fire, and the smoke highly influence the pedestrian evacuation. Weifeng and Hai simulated the crowd behavior in a smoke-filled compartment [8]. The pedestrians’ visibility is affected by the smoke concentration in their model. Makmul [5] presented a cellular automaton (CA) model to simulate a pedestrian evacuation while gaseous hazardous material is spreading. It is combined with the advection-diffusion equation and the Eikonal equation. The advection-diffusion is used for hazard gas spreading. The Eikonal equation is applied for the arrival times of CA cells.In this work, we are interested in the evacuation of pedestrians during the spread of smoke. A social force model is applied to simulate pedestrian movements, which is one of the most popular choices to model the dynamic behavior of pedestrians. This model produces realistic movements in simulations by taking the individual characteristics into account, such as the desired speed, desired destination, and physical interaction factors. In our model, the social force model is combined with the advection-diffusion equation, which is applied to the smoke density. This equation is used in many applications for fluid motion, heat transfer, flow of gas, or pollutant [9]. The operator splitting method is applied to solve this equation, which is an efficient approach to solve problems in multidimensions. The movement directions of an individual are defined by the negative gradient of the solution of the Eikonal equation in our model. The Eikonal equation is also incorporated into the SF model and solved numerically by the fast marching method [10].Our main objective of this paper is to comprise the advection-diffusion equation [9] and the Eikonal equation [11] into a social force model. The expanded model is used to describe the pedestrian flow while smoke is propagating. The movement direction of a pedestrian is determined by the negative gradient of the solution of the Eikonal equation. The solution of the Eikonal equation depends on the smoke density as well as on the pedestrian density in the individual’s visibility distance. The range of an individual’s visibility varies all the time relying on the smoke density. We perform numerical experiments in a room with multiple exits. The influence of the pedestrian density and the smoke source location on the evacuation time is examined and discussed.The framework of this paper is organized as follows. In Section2, we present the social force model and explain a way to couple it with the advection-diffusion and the Eikonal equation. The numerical methods, which are used to approximate the solutions of the social force model, the advection-diffusion equation, and the Eikonal equation are demonstrated in Section 3. Numerical experiments and results are shown in Section 4. Finally, the discussion and conclusions are made in Section 5.


## 2. Model

We study pedestrian evacuations in a multiple exits domain with a source of smoke. All pedestrians in our model are assumed to be familiar with the environment, and the smoke density has no effect on the pedestrians’ health. A microscopic social force model without boundary force [2] is introduced to describe the pedestrians’ movements. The moving direction of an individual in our model is guided by the solution of the Eikonal equation, which relies on the smoke density as well as on the pedestrian density in the visibility distance. The linear advection-diffusion equation is applied for the propagation of smoke. The microscopic social force equations of the dynamic interacting particle systems together with the Eikonal equation and the advection-diffusion equation in two dimensions read as
(1)dxidt=vit,(2)midvidt=fidt+∑j≠ifijsoct+fijpht+ξit,with locationxi∈ℝ2 and velocity vi∈ℝ2, i=1,2,⋯,N. N is the number of all pedestrians. mi is the mass of pedestrian i, and ξit is an individual fluctuation term indicating chaotic behavioral changes. For simplification, we set mi=1 and ξi=0 for all i=1,⋯,N.fidt is the desire force that indicates the intention of a pedestrian to reach his destination with a certain desired speed vd in a given desired direction ed. Its mathematical expression for pedestrian i is described by
(3)fidt=vidteidt−vitτi,wherevidt is the intended speed of pedestrian i at time t. In our model, the intended speed of a pedestrian in a smoke region depends on the pedestrian density in the visibility distance. It is set to
(4)vidt=Umax1−ρxi,tρmax,(5)ρxi,t=1πSt2∑j,xi−xj<St1,Umax and ρmax are the maximum speed and the density of pedestrians, respectively. St is the visibility distance of a pedestrian in smoke areas at time t. ρx,t is the pedestrian density in a circle with s radius St at time t. Experiments on human behavior in fire smoke show that the actual visibility distance captured by light reflecting objects can be described through the following equation [7, 12]
(6)St=3Kt,Kt=Kpmst,Kt is the extinction coefficient of smoke at time t, which is measured in m−1. Kp is a constant, and it is set to 7.6 m2/g for light-reflecting objects such as human bodies. mst is the smoke mass at time t.eidt is the unit vector pointing to the desired direction. vit is the actual velocity, and τi is the relaxation time in which the pedestrian adapts his actual velocity to the intended velocity [vidteidt]. In our model, the pedestrians intend to minimize their estimated travel time towards their destination but temper this behavior to avoid areas with high pedestrian and smoke densities. The desired movement direction of pedestrian i in our model is given by
(7)eidt=−∇Txi∣∇Txi∣,whereTxi is the travel cost of the pedestrians to reach their destination at point xi. It is the solution of the following Eikonal equation [11]
(8)∣∇Tx∣=1Fx,x∈Ω,whereΩ is a simulation domain and Tx is the arrival time of the front crossing the point x. Tx is set to 0 for the destination areas. Fx>0 is the moving speed of the front and depends on the position of x. It is set as the follows
(9)Fx=0.01,x∈ΩbUρx,,x∈Ω/Ωb,where Ωb are areas that are obstructed by obstacles [11] or areas with high smoke density. We set the threshold smoke concentration as Cx≥0.05. It is assumed to be high enough to influence the movements of the pedestrians. This value reduces the visibility range of a pedestrian to about 6 m-7 m as described in the references [13, 14]. U is the speed-density function. It describes the relationship between the speed and the density of pedestrians. Presently, there are many speed-density functions available. In our simulation, we choose the following [15]
(10)Uρx,t=Umax1−ρx,tρmax,ρx,t=1πSt2∑j,x−xj<St1.The termfijsoct refers to the repulsive social force. It expresses the psychological tendency of pedestrians that prefer to have a personal space and keep a certain safety distance to other pedestrians. The repulsive social force is a monotonic force depending on the interpedestrian distance. It is modeled as an exponential decaying function
(11)fijsoct=Aiexprij−dijBinijλi+1−λi1+cosφij2.Here,Ai is the interaction strength, Bi is the range of the repulsive interactions, and dij=∣xi−xj∣ is the distance between the centres of mass of the pedestrians i and j. rij=ri+rj is the sum of the pedestrians’ radii ri and rj, and nijt=nij1,nij2=xit−xjt/dijt is the normalized vector pointing in the direction from pedestrian j to pedestrian i. λi is a value in the range [0, 1], and λi<1 reflects the behaviour of a pedestrian to react stronger to an event in front of him than to an event behind him. cosφij=−nijt·eit, where eit=vit/vit is the direction of motion of pedestrian i and φijt denotes the angle between the direction of motion of pedestrian i (eit) and the direction to pedestrian j, which exerts a repulsive force on pedestrian i, i.e., −nijt.fijpht reflects the physical interaction force. It focuses on the aspect that pedestrians can get in physical contact with each other (rij≥dij) in crowded environments. This force separates the pedestrians when collisions occur. The physical interaction force is expressed by
(12)fijpht=knHrij−dijnij+ktHrij−dijΔvjittij,whereH is the Heaviside function, tij=−nij2,nij1 is the unit tangential vector and orthogonal to nij, Δvjit=vj−vi·tij is the tangential velocity difference, and kn and kt are the normal and tangential constants, respectively. The first term on the right-hand side represents a body force for the body compression effect while the second term stands for a sliding friction force that occurs in relative tangential motion.In the proposed model, the gas density is integrated into the SF model and the Eikonal equation. The development of smoke is expressed through the following linear advection-diffusion equation [9]
(13)∂C∂t+w·∇C=κd∇2C+Scs,t∈Ω×ℝ+,with the Dirichlet boundary conditionsC=0 on ∂Ω, the diffusion constant κd>0, the velocity field w=w1,w2∈ℝ2, the source term Scs,t, and the simulation domain Ω⊂ℝ2. We assume that the smoke source emits gas at a constant rate Qcg/s from a single source point cs=xs,ys. Hence, we write the source term as
(14)Scs=Qcδx−xsδy−ys,whereδ is the Dirac delta function given by
(15)δx=1,x=00,x≠0.


## 3. Numerical Methods

This section is devoted to numerical methods that are used to discretise the equations of the SF model (1) and (2), the Eikonal equation (8), and the advection-diffusion equation (13).


### 3.1. Solving the SF Model

To solve the SF model, first, we write the equations (1) and (2) as
(16)u′t=ft,ut,ut0=u0,with initial conditionut0=u0. The equidistant grid Ωt with respect to time t is generated by
(17)Ωt=tk,tk=kΔt,k=0,1,2,⋯,MandΔt=1M.Then, we apply the two-stage second-order Runge-Kutta method as follows(18)k1=ftk,uk,(19)k2=ftk+23Δt,yk+23Δtk1,(20)uk+1=uk+Δt14k1+34k2,whereuk=utk and uk+1=utk+1. The solution of u in the next time step is obtained from equation (20).


### 3.2. Solving the Eikonal Equation

There are a quite number of numerical methods existing to approximate the solution of the Eikonal equation: the fast marching method [10], the fast marching level set method [10], the fast sweeping method [16], and the fast iterative method [17] for example. In our experiments, the fast marching method is employed in all simulations. The details of this method can be reviewed in reference [10].


### 3.3. Solving the Advection-Diffusion Equation

We discretise the equation of the advection-diffusion (13) in order to obtain the smoke density on the simulation domain at each time step. Equidistant grids with respect to time and space are generated as follows
(21)tk=kΔt,k=0,1,2,⋯,M,Δt=1M,(22)xa=aΔx,a=0,1,2,⋯,Nx,Δx=1Nx,(23)yb=bΔy,b=0,1,2,⋯,Ny,Δy=1Ny.Then, we perform the operator splitting on the two-dimensional advection-diffusion equation. Thex-direction and the y-direction are solved separately over two time steps. The splitting method on equation (13) reads as
(24)Ca,b∗=Ca,bk−Δtw1∗Dx−C∗−κdDxxC∗,(25)Ca,bk+1=Ca,b∗−Δtw2t+1Dy−Ct+1−κdDyyCt+1+ΔtSa,bk,where
(26)Dx−C∗=Ca,b∗−Ca−1,b∗Δx,DxxC∗=Ca+1,b∗−2Ca,b∗+Ca−1,b∗Δx2,Dy−Ck+1=Ca,bk+1−Ca,b−1k+1Δy,DyyCk+1=Ca,b+1k+1−2Ca,bk+1+Ca,b−1k+1Δy2.We assume that the smoke density on the boundaries is zero. Thus, we have(27)C0,y,t=0,CNx,y,t=0,Cx,0,t=0,Cx,Ny,t=0,with the initial condition given by
(28)Cx,y,0=Qcδx−xsδy−ys,wherexs,ys is a hazard source point and Qc is the gas concentration at the initial time. First, we solve implicitly equation (24) in the x direction and obtain Ca,b∗ for all grid points a,b. The value Ca,b∗ is used for equation (25) in y direction. Equation (3.6) is solved implicitly. Finally, we receive the smoke density for the next time step Ca,bk+1 for all grid points a,b. For further details and the convergence of this method, we refer to reference [18].


### 3.4. Main Update Algorithm

A way to obtain an individuals’ trajectory in each time step is described in Algorithm1. The algorithm, that is used to update an individuals’ position and velocity in each time step, is shown below.Algorithm 1
Main update algorithm.input: xi0;vi0;Cx,0;Tx,0;ρmax;N;tendoutput: xit;vit;Cx,t;Tx,tfort=1totenddofors=1# grid points doSolve the advection-diffusion equation (13)Find grid points with high smoke density, i.e.,Cx≥0.05fori=1 to Ndo1. Compute the visibility distance of personi (St);2. Compute the pedestrian density (ρx,t) around person i in a ball radius St;3. SetTx=0, where x is a grid point on an exit;SetFx=Umax1−ρx,t/ρmax for walkable areas;SetFx=0.01 for obstacle cells or cells with high smoke density Cx≥0.05;4. Solve the Eikonal equation (8);5. Seteidt=−∇Txi/∣∇Txi∣;6. Computefidt;forj=1 to NdoComputefijsoct;Computefijphyt);Solve the social force model (1) and (2);Update the positions and the velocities of all persons;Remark 1.
We set the threshold smoke concentration asCx≥0.05 which is assumed to be high enough to influence the movements of the pedestrians. This value reduces the visibility range of a pedestrian to about 6 m-7 m, and his speed is affected significantly, as described in the references [13, 14].


## 3.1. Solving the SF Model

To solve the SF model, first, we write the equations (1) and (2) as
(16)u′t=ft,ut,ut0=u0,with initial conditionut0=u0. The equidistant grid Ωt with respect to time t is generated by
(17)Ωt=tk,tk=kΔt,k=0,1,2,⋯,MandΔt=1M.Then, we apply the two-stage second-order Runge-Kutta method as follows(18)k1=ftk,uk,(19)k2=ftk+23Δt,yk+23Δtk1,(20)uk+1=uk+Δt14k1+34k2,whereuk=utk and uk+1=utk+1. The solution of u in the next time step is obtained from equation (20).


## 3.2. Solving the Eikonal Equation

There are a quite number of numerical methods existing to approximate the solution of the Eikonal equation: the fast marching method [10], the fast marching level set method [10], the fast sweeping method [16], and the fast iterative method [17] for example. In our experiments, the fast marching method is employed in all simulations. The details of this method can be reviewed in reference [10].


## 3.3. Solving the Advection-Diffusion Equation

We discretise the equation of the advection-diffusion (13) in order to obtain the smoke density on the simulation domain at each time step. Equidistant grids with respect to time and space are generated as follows
(21)tk=kΔt,k=0,1,2,⋯,M,Δt=1M,(22)xa=aΔx,a=0,1,2,⋯,Nx,Δx=1Nx,(23)yb=bΔy,b=0,1,2,⋯,Ny,Δy=1Ny.Then, we perform the operator splitting on the two-dimensional advection-diffusion equation. Thex-direction and the y-direction are solved separately over two time steps. The splitting method on equation (13) reads as
(24)Ca,b∗=Ca,bk−Δtw1∗Dx−C∗−κdDxxC∗,(25)Ca,bk+1=Ca,b∗−Δtw2t+1Dy−Ct+1−κdDyyCt+1+ΔtSa,bk,where
(26)Dx−C∗=Ca,b∗−Ca−1,b∗Δx,DxxC∗=Ca+1,b∗−2Ca,b∗+Ca−1,b∗Δx2,Dy−Ck+1=Ca,bk+1−Ca,b−1k+1Δy,DyyCk+1=Ca,b+1k+1−2Ca,bk+1+Ca,b−1k+1Δy2.We assume that the smoke density on the boundaries is zero. Thus, we have(27)C0,y,t=0,CNx,y,t=0,Cx,0,t=0,Cx,Ny,t=0,with the initial condition given by
(28)Cx,y,0=Qcδx−xsδy−ys,wherexs,ys is a hazard source point and Qc is the gas concentration at the initial time. First, we solve implicitly equation (24) in the x direction and obtain Ca,b∗ for all grid points a,b. The value Ca,b∗ is used for equation (25) in y direction. Equation (3.6) is solved implicitly. Finally, we receive the smoke density for the next time step Ca,bk+1 for all grid points a,b. For further details and the convergence of this method, we refer to reference [18].


## 3.4. Main Update Algorithm

A way to obtain an individuals’ trajectory in each time step is described in Algorithm1. The algorithm, that is used to update an individuals’ position and velocity in each time step, is shown below.Algorithm 1
Main update algorithm.input: xi0;vi0;Cx,0;Tx,0;ρmax;N;tendoutput: xit;vit;Cx,t;Tx,tfort=1totenddofors=1# grid points doSolve the advection-diffusion equation (13)Find grid points with high smoke density, i.e.,Cx≥0.05fori=1 to Ndo1. Compute the visibility distance of personi (St);2. Compute the pedestrian density (ρx,t) around person i in a ball radius St;3. SetTx=0, where x is a grid point on an exit;SetFx=Umax1−ρx,t/ρmax for walkable areas;SetFx=0.01 for obstacle cells or cells with high smoke density Cx≥0.05;4. Solve the Eikonal equation (8);5. Seteidt=−∇Txi/∣∇Txi∣;6. Computefidt;forj=1 to NdoComputefijsoct;Computefijphyt);Solve the social force model (1) and (2);Update the positions and the velocities of all persons;Remark 1.
We set the threshold smoke concentration asCx≥0.05 which is assumed to be high enough to influence the movements of the pedestrians. This value reduces the visibility range of a pedestrian to about 6 m-7 m, and his speed is affected significantly, as described in the references [13, 14].


## 4. Numerical Experiments and Results

In this section, we investigate numerical experiments which are performed in a room of size16m×20m. Two exits of width 2 m are located at the bottom and on the right side of the room. They are labeled as Exit 1 and Exit 2, respectively, see Figure 1. Since our model is nondeterministic, different runs will present different results. Thus, to obtain more general results, we simulate ten trial runs for the same example and their average is used. The computations are executed on a HP Intel(R) Core(TM) i7-7700CPU, 3.6GHz. All programs are implemented in MATLAB.Figure 1
Modeling area for the numerical experiments.To examine the effect of smoke spreading on an individual’s movements, we limit our study to only one source and assume that the smoke density at the source point is rather high at the initial time. Furthermore, the source emits a constant smoke density of 0.01 g/s afterward, i.e.,(29)Qc=10g/s,t=00.01g/s,t>0.The velocity fieldw1,w2 in equation (13) is chosen randomly out of the interval [-0.5, 0.5] at each time step. The maximum velocity Umax is set to 3 m/s, and the maximum density ρmax is 10 ped/m2. At initial time, the pedestrians are randomly distributed in the room, and their velocities are set to zero.The model parameters that are used in the evacuation simulations are shown in Table1. We simulate the individuals’ temptation to move away from a high smoke density area and to simultaneously move towards an exit to evacuate. Two scenarios are regarded to examine the influence of the smoke source position on the pedestrian evacuation. In scenario (i), the smoke source’s location is in the middle of the room, whereas it is located in front of Exit 1 in scenario (ii). There are 100, 300, and 500 individuals in experiments.Table 1
Parameters for evacuation simulations.
ParameterSymbolValueEquationRelaxation timeτi0.5(3)Maximum desired speedUmax3(4)Maximum pedestrian densityρmax10(4)Human-human interaction strengthAi2(11)Human-human range of repulsive interactionBi0.21(11)Contact distancerij0.5(11)Anisotropic parameterλi0.61(11)Body force coefficientkn2(12)Friction force coefficientkt2(12)Velocity field inx-directionw1[-0.5, 0.5](13)Velocity field iny-directionw2[-0.5, 0.5](13)Diffusion constantκd0.05(13)Time step sizeΔt0.02(20), (21)Space grid size inxΔx0.4(22)Space grid size inyΔy0.4(23)The average evacuation time and the percentage of pedestrians who move out through Exit 1 in scenario (i) and in scenario (ii) are presented in Table2. We see that the average evacuation time in scenario (ii) is longer than in scenario (i) in all cases. Therefore, the location of the smoke source affects the pedestrian’s evacuation time. It takes more time to evacuate when the source is placed in front of Exit 1 compared to when the source is in the middle of the room. A high number of evacuees move out through Exit 2 in scenario (ii) leading to a large jamming and crowd around this exit. In scenario (ii), the pedestrians reject to move out through Exit 1 because the source is in front of this exit and it is too dangerous to move out.Table 2
Average evacuation time and percentage of pedestrians who move out through Exit 1. The source is in the middle of the room in scenario (i) and in front of Exit 1 in scenario (ii).
Number of pedestriansEvacuation time (s)Exit 1 (%)Scenario (i)Scenario (ii)Scenario (i)Scenario (ii)1007.4360 s13.3660 s56.702.6030013.0800 s22.3920 s170.004.9050020.0180 s32.5400 s285.209.30The plot of the average number of outside pedestrians of each exit versus the time in scenario (i) and scenario (ii) with 300 individuals are demonstrated in Figure2. The usage of Exit 1 is slightly higher than the usage of Exit 2 in scenario (i), in which the smoke source is placed in the center of the room. The pedestrians evacuate less out of Exit 1 in scenario (ii).Figure 2
Number of outside pedestrians in scenario (i) and scenario (ii). The total number of pedestrians is 300.
(a)
Source is in the middle(b)
Source is near Exit 1The comparisons of the average number of outside pedestrians in scenario (i) and scenario (ii) with 100, 300, and 500 pedestrians in the experiments are displayed in Figure3. Both scenarios show that as the number of pedestrians in experiments rises, the average number of outside pedestrians also increases. The time that is needed to evacuate all individuals out of the room increases when the number of pedestrians in simulations get higher. Scenario (ii) takes more time to evacuate all people compared to scenario (i).Figure 3
Average number of outside pedestrians in scenario (i) and scenario (ii). The number of pedestrians is 100, 300, and 500.
(a)
Source is in the middle(b)
Source is near Exit 1Table3 shows the computation time of ten trial runs in scenario (i) and in scenario (ii) with 100, 300, and 500 pedestrians. In all cases, the simulations in scenario (ii) take more computation time compared to scenario (i). This is because the evacuation times in scenario (ii) are larger than in scenario (i), see Table 2.Table 3
The computation time of ten trial runs in scenario (i) and in scenario (ii) with 100, 300, and 500 evacuees.
Number of pedestriansCPU time (10 rounds)Scenario (i)Scenario (ii)1002.9956 h3.3485 h30010.6187 h14.8713 h50022.5628 h34.3668 hFigure4 demonstrates the movement directions of 300 evacuees during smoke spreading, where the smoke source is in the center of the room. Red circles refer to the positions of pedestrians, and green arrows show the moving directions. At time t=0, all pedestrians are placed randomly inside the room and their velocities are zero. There is no green arrow points out from any pedestrian. As time goes on, the individuals start to move. The moving direction of each pedestrian points in the direction towards the exits and points out from the smoke region. Clogging and arching effects are observed around the two exits as we can see in Figures 4(d)–4(f). The number of individuals who evacuate through Exit 1 is slightly higher than the number of individuals who move towards Exit 2.Figure 4
Movements of 300 pedestrians during smoke spreading, where the source is in the middle of the room. The velocity fieldw1,w2 of the advection-diffusion equation (13) is chosen randomly out of the interval [-0.5,0.5].
(a)(b)(c)(d)(e)(f)The movements of 300 pedestrians during the spreading of smoke in scenario (ii), where the source is at Exit 1, are illustrated in Figure5. At the beginning of the simulation, the smoke density is less and has no effect on the individuals’ movements. The moving directions of some pedestrians point to Exit 1 where the smoke source is situated. Pedestrians, who are near Exit 1, go out through this exit, see Figure 5(b). While the time increases, the smoke density around Exit 1 increases. The moving directions of the individuals change. They point out of the smoke source region and point mostly towards Exit 2, see Figure 5(c). Most of the pedestrians reject to evacuate through Exit 1 since it is dangerous. Large crowds and jammings occur around Exit 2, see Figure 5(f).Figure 5
The movements of the pedestrians during smoke spreading, where the source is at Exit 1. The velocity fieldw1,w2 of the advection-diffusion equation (13) is chosen randomly out of the interval [-0.5,0.5].
(a)(b)(c)(d)(e)(f)The trajectories of the 300 evacuees during propagation of smoke are presented in Figure6, in which the smoke source is in the middle of the room and the velocity field of the advection-diffusion is w1,w2=0.5,0 in equation (13).Figure 6
The movements of 300 pedestrians during smoke spreading, as the velocity field of the advection-diffusion isw1,w2=0.5,0 in equation (13) and the source is in the middle of the room.
(a)(b)(c)(d)(e)(f)Figure7 displays the contour plot of the traveling time (arrival time) to reach a destination of the 300 individuals in the simulations considering scenario (i) and scenario (ii) at t=5s. The traveling cost is obtained by solving the Eikonal equation (8) with the help of the fast marching method. To plot the traveling time, we set Tx=50 for grid space that receives traveling cost higher than 50 (Tx>50), such as high smoke density grids or walls grids. In Figure 7, we see that the arrival times are relatively high in areas where the smoke source or the wall grids are placed. This is on account of combining the smoke density and obstacles to the speed of the front Fx in equation (2.6).Figure 7
Contour plots of the dynamic potential (estimated remaining travel time) of 300 pedestrians in scenario (i) and scenario (ii) att=5s.
(a)
Scenario (i)(b)
Scenario (ii)Negative gradients of the arrival times (−∇Tx/∣∇Tx∣) of 300 individuals for the example in scenario (i) and scenario (ii) at t=5s are shown in Figure 8. They demonstrate the movement directions at each point on the simulation domain which guides the way to move for the evacuees. In Figure 8, we see that negative gradients of the arrival times point out of the smoke source and wall areas. They point towards the exits.Figure 8
Negative gradient of the arrival times of 300 pedestrians in scenario (i) and scenario (ii) att=5s.
(a)
Scenario (i)(b)
Scenario (ii)


## 5. Discussion and Conclusions

In this paper, we present a model that is used to describe pedestrians’ movements during the propagation of smoke. The social force model is applied and coupled with the advection-diffusion equation and the Eikonal equation. The advection-diffusion equation is employed for smoke spreading. The Eikonal equation is used for the moving directions of pedestrians. It depends on the smoke as well as on the pedestrian density in the visibility range. In our model, we observed the arching and the clogging effect which are pedestrian collective phenomena occurring in pedestrian crowds around exits [2, 19]. We performed different experiments and showed their results.



---

*Source: 8819076-2020-11-16.xml*