# American Institute of Mathematical Sciences

doi: 10.3934/krm.2020044

## An anisotropic interaction model with collision avoidance

 University of Mannheim, School of Business Informatics and Mathematics, 68159 Mannheim, Germany

* Corresponding author: Claudia Totzeck

Received  December 2019 Revised  July 2020 Published  September 2020

In this article an anisotropic interaction model avoiding collisions is proposed. Starting point is a general isotropic interacting particle system, as used for swarming or follower-leader dynamics. An anisotropy is induced by rotation of the force vector resulting from the interaction of two agents. In this way the anisotropy is leading to a smooth evasion behaviour. In fact, the proposed model generalizes the standard models, and compensates their drawback of not being able to avoid collisions. Moreover, the model allows for formal passage to the limit 'number of particles to infinity', leading to a mesoscopic description in the mean-field sense. Possible applications are autonomous traffic, swarming or pedestrian motion. Here, we focus on the latter, as the model is validated numerically using two scenarios in pedestrian dynamics. The first one investigates the pattern formation in a channel, where two groups of pedestrians are walking in opposite directions. The second experiment considers a crossing with one group walking from left to right and the other one from bottom to top. The well-known pattern of lanes in the channel and travelling waves at the crossing can be reproduced with the help of this anisotropic model at both, the microscopic and the mesoscopic level. In addition, the 'right-before-left' and 'left-before-right' rule appear intrinsically for different anisotropy parameters.

Citation: Claudia Totzeck. An anisotropic interaction model with collision avoidance. Kinetic & Related Models, doi: 10.3934/krm.2020044
##### References:

show all references

##### References:
Illustration of isotropic and anisotropic pairwise interaction induced by different values of $\lambda$. The initial positions of the pedestrians are highlighted with a point. Top: Two pedestrians are approaching each other. Positive $\lambda$ leads to right-stepping, negative $\lambda$ results in left-stepping. The isotropic case corresponds to $\lambda = 0.$ Bottom: Two pedestrians meet at a crossroad. In the isotropic case, both depart from their desired trajectory. For $\lambda \ne 0$ the two return to their desired velocity after avoiding the collision. In case of the crossing with $\lambda = -0.25$ the blue pedestrians slows down, while the red pedestrian accelerates to avoid the collision. Then, the pedestrians relax their velocities towards the desired velocity. For $\lambda = 0.25$ the two pedestrians change their roles
Visualization of the influence of $\lambda.$ For $\lambda = 0$ we are in the case of the grey vector. The interaction force reduces the force resulting for the desired velocity. The particles may stop in front of each other. For $\lambda > 0$ we are in the yellow region. The force vector resulting from the interaction is turned. Adding this rotated vector and the vector resulting from the desired velocity, yields the evasive behaviour. The red particle moves to the bottom and the blue particles moves to the top. For $\lambda < 0$ the roles of the to particles change. We are then in the blue region and the red particle moves to the top and the blue particle to the bottom
Visualization of the influence of $\lambda$ in a crossing scenario. Top: For $\lambda = 0$ we are in the case of the grey vector. The particles push each other diagonally and leave the path given by the desired velocity, see Figure 3 in the middle. For $\lambda > 0$ we are in the yellow region. The force vector resulting from the interaction is turned. Adding this rotated vector and the vector resulting from the desired velocity, yields the evasive behaviour. The red particle slows down and moves to the top before accelerating and adjusting the velocity to get back onto its path and the blue particle does not slow down as much as the red one and moves to the right. For $\lambda < 0$ the roles of the to particles change. We are then in the blue region and the red particle moves to the top and the blue particle moves to the right before being push back onto its path. Intrinsically, we see here the 'right-before-left' rule for $\lambda > 0$ and 'left-before-right' for $\lambda <0.$
Illustration of the reflective boundary conditions. The black particle with velocity depicted with the dashed vector is about to leave the domain in the next time step. Due to the reflecting boundary conditions, it is projected into the domain and the y-component of its velocity is reflected (black vector)
Illustration of the boundary setting. The green and yellow parts of the boundary refer to periodic boundary conditions. Walls are indicated by black lines. In case of the crossing the particles that leave the domain through a green boundary are flowing in at the other green boundary and analogous for the yellow boundary parts
Left: Illustration of artificial agents with fixed position and artificial velocities modelling a circular obstacle. Right: Influence of a circular obstacle on the trajectories of the pedestrians
Lane formation in a channel - particle simulation with only few particles involved, we see a formation of multiple horizontal lanes as stationary state. The parameters are $N_b = 10 = N_r, t = 150.$ The arrows show the velocities of the particles
Lane formation in a channel - particle simulation. The blue pedestrians go from right to left, the red from left to right, i.e. $u_b = (-0.2,0)^T$ and $u_r = (0.2,0)^T.$ The arrows show the velocities of the pedestrians. As $\lambda = 0.25,$ pedestrians prefer to step to the right to avoid a collision. The snapshots are made at the times $t = 0,\; 50,\; 100,\; 150,\; 200,\; 250.$
Lane formation in a channel - mean-field simulation. Initially the red and the blue group are uniformly distributed in the domain. The difference of the densities vanishes as shown in the plot on the top-left. Then we see the formation of diagonal stripes at time t = 50 and finally a formation of lanes as time proceeds. The plots correspond to times t = 0,250,500,750, 1000, 1250
Density distribution in a channel - mean-field simulation. The densities $\rho_i^2, i \in \{r, b\}$ are integrated w.r.t. $x^1$ in order to extract the information along the $x^2$-axis
Lane formation in a channel - mean-field simulation. Initially the red and the blue group are uniformly distributed in the velocity domain $\Omega_{v_r}$ and $\Omega_{v_b}$, respectively, as shown in the plot on the top-left. Then the relaxation towards the desired velocities $u_r = (0.2,0)^T$ and $u_b = (-0.2,0)^T$ starts, see plot on the top-right. Only small changes in the $y$-coordinates lead the crowd to the lane configuration as can be seen when comparing the top-left and the bottom plot of the velocity density. The time instances are $t = 0,\; 250,\; 1500.$
Travelling waves at a crossing - particle simulation. The blue go from bottom to top, the red from left to right, i.e. $u_b = (0,0.2)$ and $u_r = (0.2,0).$ The arrows show the velocities of the pedestrians. As $\lambda = 0.25,$ pedestrians prefer to step to the right to avoid a collision. From top-left to bottom-right: $t = 0$, $t = 50$, $t = 100$, $t = 150$, $t = 200$ and $t = 250$. The simulation was done with $N_b = 150 = N_r.$ At $t = 250$ the pattern allows most pedestrians to walk with their desired velocities
Travelling waves at a crossing - mean-field simulation. Initially the red and the blue group are uniformly distributed in the spatial domain. The difference of the densities vanishes as shown in the plot on the top-left. At time t = 50 the groups start to separate, see the plot in the top-right. Afterwards the line pattern is forming (t = 150,375 and 750). Finally, we see the stationary configuration of travelling waves at time t = 1500 in the bottom-right plot
Travelling waves at a crossing - mean-field simulation. Initially the red and the blue group are uniformly distributed in the velocity domain $\Omega_v$ as shown in the plot on the top-left. Then the relaxation towards the desired velocities $u_r = (0.2,0)^T$ and $u_b = (0, 0.2)^T$ starts, see plot on the top-right. Only small changes in the $x$- and $y$-coordinate lead the crowd to the travelling waves configuration as can be seen when comparing the top-left and the bottom plot of the velocity density. The time instances are $t = 0,250, 1500.$
 [1] Michael Herty, Mattia Zanella. Performance bounds for the mean-field limit of constrained dynamics. Discrete & Continuous Dynamical Systems - A, 2017, 37 (4) : 2023-2043. doi: 10.3934/dcds.2017086 [2] Martin Burger, Marco Di Francesco, Peter A. Markowich, Marie-Therese Wolfram. Mean field games with nonlinear mobilities in pedestrian dynamics. Discrete & Continuous Dynamical Systems - B, 2014, 19 (5) : 1311-1333. doi: 10.3934/dcdsb.2014.19.1311 [3] Rong Yang, Li Chen. Mean-field limit for a collision-avoiding flocking system and the time-asymptotic flocking dynamics for the kinetic equation. Kinetic & Related Models, 2014, 7 (2) : 381-400. doi: 10.3934/krm.2014.7.381 [4] Franco Flandoli, Enrico Priola, Giovanni Zanco. A mean-field model with discontinuous coefficients for neurons with spatial interaction. Discrete & Continuous Dynamical Systems - A, 2019, 39 (6) : 3037-3067. doi: 10.3934/dcds.2019126 [5] Gerasimenko Viktor. Heisenberg picture of quantum kinetic evolution in mean-field limit. Kinetic & Related Models, 2011, 4 (1) : 385-399. doi: 10.3934/krm.2011.4.385 [6] Seung-Yeal Ha, Jeongho Kim, Jinyeong Park, Xiongtao Zhang. Uniform stability and mean-field limit for the augmented Kuramoto model. Networks & Heterogeneous Media, 2018, 13 (2) : 297-322. doi: 10.3934/nhm.2018013 [7] Charles Bordenave, David R. McDonald, Alexandre Proutière. A particle system in interaction with a rapidly varying environment: Mean field limits and applications. Networks & Heterogeneous Media, 2010, 5 (1) : 31-62. doi: 10.3934/nhm.2010.5.31 [8] Seung-Yeal Ha, Jeongho Kim, Peter Pickl, Xiongtao Zhang. A probabilistic approach for the mean-field limit to the Cucker-Smale model with a singular communication. Kinetic & Related Models, 2019, 12 (5) : 1045-1067. doi: 10.3934/krm.2019039 [9] Young-Pil Choi, Samir Salem. Cucker-Smale flocking particles with multiplicative noises: Stochastic mean-field limit and phase transition. Kinetic & Related Models, 2019, 12 (3) : 573-592. doi: 10.3934/krm.2019023 [10] Seung-Yeal Ha, Jeongho Kim, Xiongtao Zhang. Uniform stability of the Cucker-Smale model and its application to the Mean-Field limit. Kinetic & Related Models, 2018, 11 (5) : 1157-1181. doi: 10.3934/krm.2018045 [11] Joachim Crevat. Mean-field limit of a spatially-extended FitzHugh-Nagumo neural network. Kinetic & Related Models, 2019, 12 (6) : 1329-1358. doi: 10.3934/krm.2019052 [12] Max-Olivier Hongler. Mean-field games and swarms dynamics in Gaussian and non-Gaussian environments. Journal of Dynamics & Games, 2020, 7 (1) : 1-20. doi: 10.3934/jdg.2020001 [13] Franco Flandoli, Matti Leimbach. Mean field limit with proliferation. Discrete & Continuous Dynamical Systems - B, 2016, 21 (9) : 3029-3052. doi: 10.3934/dcdsb.2016086 [14] Patrick Gerard, Christophe Pallard. A mean-field toy model for resonant transport. Kinetic & Related Models, 2010, 3 (2) : 299-309. doi: 10.3934/krm.2010.3.299 [15] Thierry Paul, Mario Pulvirenti. Asymptotic expansion of the mean-field approximation. Discrete & Continuous Dynamical Systems - A, 2019, 39 (4) : 1891-1921. doi: 10.3934/dcds.2019080 [16] Doron Levy, Tiago Requeijo. Modeling group dynamics of phototaxis: From particle systems to PDEs. Discrete & Continuous Dynamical Systems - B, 2008, 9 (1) : 103-128. doi: 10.3934/dcdsb.2008.9.103 [17] Andreas Schadschneider, Armin Seyfried. Empirical results for pedestrian dynamics and their implications for modeling. Networks & Heterogeneous Media, 2011, 6 (3) : 545-560. doi: 10.3934/nhm.2011.6.545 [18] Diogo A. Gomes, Gabriel E. Pires, Héctor Sánchez-Morgado. A-priori estimates for stationary mean-field games. Networks & Heterogeneous Media, 2012, 7 (2) : 303-314. doi: 10.3934/nhm.2012.7.303 [19] Illés Horváth, Kristóf Attila Horváth, Péter Kovács, Miklós Telek. Mean-field analysis of a scaling MAC radio protocol. Journal of Industrial & Management Optimization, 2019  doi: 10.3934/jimo.2019111 [20] Yufeng Shi, Tianxiao Wang, Jiongmin Yong. Mean-field backward stochastic Volterra integral equations. Discrete & Continuous Dynamical Systems - B, 2013, 18 (7) : 1929-1967. doi: 10.3934/dcdsb.2013.18.1929

2019 Impact Factor: 1.311

## Metrics

• HTML views (36)
• Cited by (0)

• on AIMS