Source-linked AI summary

Pedestrian flows in bounded domains with obstacles

Benedetto Piccoli, Andrea Tosin

arXiv:0812.4390v1math-phphysics.soc-ph

TL;DR

The paper addresses how to model pedestrian flows in complex, obstacle-filled domains while representing target-seeking and crowd avoidance. It develops a discrete-time Eulerian model using time-evolving Radon measures and push-forward dynamics. The framework reproduces reported pedestrian-flow phenomena and is presented as easier to handle analytically and numerically than nonlinear hyperbolic conservation-law models, while remaining limited by its discrete-time formulation and difficulties constructing smooth desired velocities around obstacles.

  • Problem

    Pedestrian flows require models that capture target-directed motion, crowd avoidance, and complex two-dimensional domains while overcoming limitations of existing approaches.

  • Method

    The paper models pedestrian occupancy with discrete-time Radon measures evolved by push-forward motion mappings whose velocity combines desired target-directed motion with nonlocal crowd-avoidance interactions.

  • Results

    The model reproduces several experimentally documented pedestrian-flow features and offers unified analytical and numerical treatment across spatial dimensions, including two-dimensional applications.

  • Takeaways & Limitations

    Time-evolving measures provide a practical macroscopic framework for studying pedestrian motion in complex domains with obstacles.

  • Takeaways & Limitations

    The framework lacks a continuous-time limit as ∆t → 0, and constructing smooth desired velocities becomes difficult in domains with many or changing obstacles.

Abstract

from arXiv · show

In this paper we systematically apply the mathematical structures by time-evolving measures developed in a previous work to the macroscopic modeling of pedestrian flows. We propose a discrete-time Eulerian model, in which the space occupancy by pedestrians is described via a sequence of Radon positive measures generated by a push-forward recursive relation. We assume that two fundamental aspects of pedestrian behavior rule the dynamics of the system: On the one hand, the will to reach specific targets, which determines the main direction of motion of the walkers; on the other hand, the tendency to avoid crowding, which introduces interactions among the individuals. The resulting model is able to reproduce several experimental evidences of pedestrian flows pointed out in the specialized literature, being at the same time much easier to handle, from both the analytical and the numerical point of view, than other models relying on nonlinear hyperbolic conservation laws. This makes it suitable to address two-dimensional applications of practical interest, chiefly the motion of pedestrians in complex domains scattered with obstacles.

1. Introduction

Pedestrian flows are complex multi-agent systems whose behavior matters for both scientific understanding and infrastructure design. Experimental observations reveal characteristic walking patterns, but mathematical models are needed to move beyond description toward analysis of unsteady dynamics.

  • 1. Introduction: Crowd-structure coupling can create dangerous lateral oscillations, as illustrated by the London Millennium Footbridge closure.The phenomenon arose from synchronous lateral excitation exerted by pedestrians crossing the bridge.
  • 1. Introduction: Mathematical and numerical tools are needed to formalize pedestrian behavior and support engineering design and optimization.The paper situates this need within broader challenges in applied mathematics and crowd-related applications.
  • 1. Introduction: Pedestrians balance reaching targets with avoiding crowding, using mainly frontal visibility and consequently forming parallel lanes.The passage contrasts this behavior with flocking in birds and fish, whose lateral eye placement provides rear visibility.
  • 1. Introduction: Experimental data provide fundamental behavioral insights but are mainly descriptive and scarcely predictive for far-from-equilibrium crowd dynamics.Unsteady dynamics can strongly affect the overall evolution of pedestrian flows.
  • 1. Introduction: The paper applies time-evolving measures to macroscopic pedestrian-flow modeling and organizes the work around theory, applications, conclusions, and research perspectives.The introduction identifies subsequent sections covering prior models, the measure framework, applications, and conclusions.

2. Overview of mathematical models of human crowds

Pedestrian-flow models range from microscopic agent dynamics to macroscopic conservation laws, each capturing crowd behavior with different analytical and computational trade-offs. The paper motivates a time-evolving-measures approach as a macroscopic Eulerian alternative for intrinsically discrete pedestrian systems.

  • 2. Overview of mathematical models of human crowds: Microscopic models track individual pedestrians and their interactions, but many coupled ODEs hinder analytical treatment and global system-level understanding.The models represent each pedestrian as an agent with a position, velocity, and interaction-related dynamics.
  • 2. Overview of mathematical models of human crowds: Existing pedestrian models encode preferred motion toward destinations together with discomfort or repulsion arising from nearby pedestrians and obstacles.Microscopic social-force and kinematic approaches use desired velocities, repulsive effects, and obstacle-related potentials.
  • 2. Overview of mathematical models of human crowds: More refined microscopic models may add attractive effects from environmental elements and stochastic fluctuations in individual behavior.These extensions are described within social-force formulations of pedestrian motion.
  • 2. Overview of mathematical models of human crowds: Related macroscopic approaches include multidimensional conservation laws and one-dimensional models with density-dependent flux behavior.The overview contrasts these formulations with the measure-based framework developed in the paper.
  • 2. Overview of mathematical models of human crowds: Macroscopic conservation-law models provide density-based descriptions but face intrinsic hyperbolic difficulties, especially in the two-dimensional settings natural for pedestrian flows.The overview presents continuity- and momentum-based formulations alongside density-dependent velocity closures.
  • 2. Overview of mathematical models of human crowds: The proposed time-evolving-measures approach evaluates pedestrian space occupancy with measures, yielding a macroscopic Eulerian description despite the system’s Lagrangian granularity.This framework is presented as a way to address difficulties associated with both microscopic models and multidimensional conservation laws.

3. Mathematical modeling by time-evolving measures

The paper models pedestrian occupancy with time-evolving positive Radon measures transported by discrete-time motion mappings. It establishes density-preservation properties, adapts the framework to outgoing flows and obstacles, and analyzes a grid-based approximation with stability and localization-error bounds.

  • Modeling framework: Positive Radon measures represent macroscopic localization and crowding without assigning physically meaningful pointwise pedestrian densities.Densities can still be introduced through absolute continuity because they are useful for analysis and numerical approximation.
  • Modeling framework: Discrete-time measure dynamics update pedestrian occupancy by pushing each measure forward through a Borel motion mapping.The mapping transports the distribution from time n to n+1, yielding µn+1 = γn#µn and a conservation-law interpretation.
  • Analytical properties: Under suitable motion-map conditions, absolutely continuous initial measures retain unique L1 densities, with mass conservation and boundedness preservation.The framework requires a non-clustering condition on γn; a Lipschitz velocity field with L < ∆t^-1 is sufficient, and the resulting density problem is well posed.
  • Numerical treatment: The numerical scheme uses grid-based piecewise translations and piecewise-constant densities, while allowing obstacle unions of grid cells under stated simplifying assumptions.The construction can be extended by dropping most simplifying assumptions, at the cost of additional implementation technicalities.
  • Numerical treatment: Stability and localization-error results relate approximate measures to prior-step or initial approximation quality and grid size, provided cell displacement is controlled by the time-step/grid-size condition.The analysis includes one-step and multistep stability, boundedness preservation, and error estimates for the approximate measures λn_h.

4. A model for pedestrian motion

The model combines a geometry-dependent desired velocity toward targets with a density-dependent interaction velocity that lets pedestrians avoid crowding. This measure-based construction supports obstacle-aware dynamics while retaining analytical regularity and well-posedness under stated assumptions.

  • 4.1. The desired velocity.: The desired velocity guides pedestrians around obstacles toward boundary targets through a scalar potential solving the Laplace equation.Boundary conditions encode attraction to targets and repulsion from walls or alternative tangential sliding along obstacles.
  • The velocity field combines target-directed motion with a crowding-avoidance correction, yielding the actual pedestrian velocity.The desired component depends on domain geometry, while the interaction component depends on the current pedestrian measure.
  • 4.1. The desired velocity.: The obstacle-aware potential construction avoids spurious aggregation points, but complex obstacle layouts make intermediate-target constructions nonsmooth and laborious.The model also requires parameter control to prevent the resultant velocity from becoming retrograde relative to the desired direction.
  • 4.2. Interaction velocity.: The interaction neighborhood is an anisotropic forward sector, modeling pedestrians as looking ahead within a finite visibility range.Its angular restriction is controlled by the maximum visibility angle, while the interaction velocity depends on the current distribution of pedestrians.
  • 4.2. Interaction velocity.: Option (i) produces a Lipschitz-continuous interaction velocity, enabling inductive existence and uniqueness results for pedestrian densities with L1 and L∞ estimates.The alternative normalized construction lacks Lipschitz continuity because of the measure in its denominator.

5. Application to some cases study

The model is applied to two-dimensional pedestrian-flow cases involving bottlenecks, obstacles, interaction patterns, and crossing streams. Across these cases, it reproduces lane formation, crowd redistribution between passages, clustering, and alternating walking lanes.

  • 5.1. Pedestrian flow through a narrow passage: The model predicts lane formation and bottleneck self-organization as pedestrians merge, obstruct passage entrances, and eventually flow through narrow corridors.High-density lanes are separated by less crowded regions, with bifurcation appearing near the bottleneck entrance.
  • 5.3. Boundary conditions: Different boundary conditions at obstacle edges produce different pedestrian-density configurations under the same initial condition and desired velocity field.Figure 7 compares Dirichlet and Neumann conditions for the potential at obstacle boundaries at the same evolution time.
  • 5.2. Pedestrian flow through two adjacent narrow passages: Pedestrians redistribute between adjacent passages: sideways crowd pressure sends some individuals to the second bottleneck, producing branching and lower density in the farther passage.The model captures the observed preference for using the passage closest to each pedestrian’s starting point.
  • 5.4. Lane formation versus clustering: Anisotropic interactions produce parallel walking lanes, whereas isotropic interactions produce clusters drifting with the desired velocity field.The characteristic spacing between lanes or clusters is comparable to the interaction-neighborhood size R = 0.1 in the simulation.
  • 5.5. Crossing flows: Oppositely directed pedestrian groups first break left-right symmetry and then form stable alternating uniformly walking lanes.The result reproduces a characteristic self-organization pattern observed in crossing pedestrian flows.
  • 5. Application to some cases study: These behavioral patterns emerge as by-products of general nonlocal-interaction principles rather than being imposed as specific modeling objectives.The applications use the time-evolving-measures model with desired motion and interaction mechanisms to generate the observed configurations.

6. Conclusions and research perspectives

The paper applies a discrete-time, measure-based framework to pedestrian flows, modeling target-seeking and crowd-avoidance behavior while supporting unified analysis and computation in arbitrary spatial dimensions. Applications reproduce experimentally observed flow features, although continuous-time dynamics remain unavailable.

  • Modeling framework: The model represents pedestrian occupancy with Radon positive measures updated by a push-forward recursion and defines motion through velocity mappings.Under non-singularity conditions, the measures admit bounded densities when the initial measure has a density.
  • Numerical approximation: Numerical approximation uses piecewise constant measures on a spatial partition, with localization error controlled under technical assumptions and a CFL-like relation between grid size and time step.The approximation is constructed for the push-forward operation and requires choosing h with respect to ∆t.
  • Analytical and numerical scope: The framework supports analytical and numerical treatment in any dimension d ≥1, avoiding additional technical difficulties for two-dimensional applications.This addresses obstacles and boundary conditions that are difficult to handle with nonlinear hyperbolic conservation-law models.
  • Pedestrian behavior: Pedestrian velocity combines target-directed motion, congestion avoidance, and boundary or obstacle corrections.The desired velocity depends on domain geometry and obstacles, while interactions represent the tendency to avoid congested areas.
  • Applications: Representative obstacle-domain applications test whether the model reproduces complex pedestrian-flow features documented in experimental research.The study includes domains with obstacles and cases involving more than one target.
  • Limitations and perspectives: A central limitation is that the time-evolving-measures framework currently lacks a continuous-time limit as ∆t → 0.The authors identify control, optimization, momentum-balance extensions, and possible turbulence modeling as directions for further research.
Loading 0812.4390v1…