Pedestrian Dynamics – Pedestrian Dynamics

Calibrating a pedestrian model against real data, part 2: where it stops

Mon, 07 Sep 2026 00:00:00 +0000

Abstract. Part 1 calibrated JuPedSim’s Collision Free Speed model with Dakota on the Hermes bottleneck experiment and found two parameter sets with near-equal residual norms. This part takes them to two other open experiments without refitting: a 0.5 m gate with the motivation varied between runs, and an entrance without guiding barriers. At the gate one set stalls in most seeds and the other passes it 20 to 40 % too fast. Recalibrating on the gate makes passage reliable but not every observable accurate, and two optimizer starts again disagree. Separate timing parameters absorb part of the low-motivation condition, in-sample, but not the speed. At the unguided entrance every set drains the crowd two to three times faster than the people did. A joint calibration over both experiments is worse than either specialist. The note closes with what we would do differently.

A calibration is only validated within the domain it was tested in. Part 1 ended with two parameter sets, A and B, that reproduce flow, density and speed at five bottleneck widths of 2.4 to 5.0 m mostly within ten percent, though neither passes the stated tolerance everywhere, and that disagree with each other on what the model gets wrong at the widest opening. The desired speed is fixed at the 1.55 m/s measured on the Hermes participants in every set below; carrying that value to other experiments is an assumption. The sets found along the way:

parameter	default	set A	set B	set C	set D	joint
desired speed [m/s]	1.2	1.55	1.55	1.55	1.55	1.55
radius [m]	0.20	0.13	0.15	0.124	0.101	0.115
time gap [s]	1.0	0.81	0.55	1.038	0.958	0.962
neighbor repulsion strength	8	2.1	9.4	9.01	1.70	8.62
neighbor repulsion range [m]	0.10	0.25	0.10	0.062	0.336	0.189
wall repulsion strength	5	5 (fixed)	5 (fixed)	2.96	1.39	4.69
wall repulsion range [m]	0.02	0.02 (fixed)	0.02 (fixed)	0.072	0.105	0.031
fitted to		Hermes 2.4, 3.6, 5.0 m	Hermes 2.4, 3.6, 5.0 m	CrowdQueue h0	CrowdQueue h0	both

Sets A and B, and sets C and D, are two optimizer starts on the same data. The tools, observables and tolerances are those of part 1: JuPedSim, PedPy on both experiment and simulation, Dakota’s surrogate-based optimizer, and a decision rule of 6 % on flow and 10 % on density and speed, chosen as calibration weights and adopted as tolerances after the fact.

Step 6 — a second experiment, and where the model stops

The CrowdQueue experiment, Wuppertal 2018: a run in the 5.6 m corridor seen from above, and the setup with the 0.5 m gate at the origin. Images: Pedestrian Dynamics Data Archive, Forschungszentrum Jülich.

The second open dataset is the CrowdQueue experiment from Wuppertal 2018 (Adrian et al. 2020): a crowd in front of an entrance again, but through a 0.5 m gate instead of a 2.4 to 5 m opening, in corridors from 1.2 to 5.6 m wide, with 11 to 75 participants and the motivation varied between runs. The archive files carry the full geometry and every trajectory, and each simulation is seeded from the data: agents in view at the first frame start there, and agents that enter the tracked corridor later are injected at the time and place of their first observation. Every tracked participant is simulated. An injection that finds its spot occupied is retried for 10 s and then dropped; this happened to one to four people in a few seeds of the two densest 1.2 m runs, and such a seed cannot count as emptied. Density and speed use the same window in experiment and simulation, the frames between the tenth and ninetieth percent of that run’s own crossings. Density is the time average over every frame; speed is averaged only over frames in which the area is occupied. Whether an observable is conditional on occupancy is part of its definition, and sparse runs expose the difference: averaging the zeros of empty frames had once reported 0.12 m/s instead of 0.56 m/s for the sparsest run. The flow is measured at a gate line and density and speed in a 2.2 m² upstream area, so the three values are not a flux identity and should not be forced to satisfy $q=\rho v$.

Each run is simulated three times and every seed is shown, because averages hide the thing that matters here. A seed counts as emptied if every expected participant was injected and passed the gate within 120 s, twice the longest experiment. Otherwise it is not emptied, and if it also had an interval of 20 s or more without a crossing while agents remained, it counts as stalled. Flow is the slope of N(t) between the tenth and ninetieth percent of crossings for emptied runs; for runs that did not empty it is the number who passed divided by the time from the first crossing to the end of the run, which is lower by construction.

N(t) at the gate for three baseline runs, experiment against the three seeds of CrowdQueue set C. The simulated slopes are steeper than the measurements. In run 110 each seed passes 62 of 63 tracked participants because one late entrant could not be placed; none is classified as stalled.

The transfer test comes first: the two Hermes parameter sets applied to the 21 runs without refitting. Set A, the long time gap with weak long-range repulsion, stalls at the funnel in 46 of 63 seeds and empties 16. Set B, the short gap with strong short-range repulsion, empties 59 of 63 and stalls in none, but where it flows it overshoots the measured gate flow by 22 % on average in the baseline runs and by 37 % in the low-motivation runs. Of the two Hermes optima with near-equal residual norms on the wide bottleneck, one mostly stalls at a half-metre gate and the other passes it too fast. That is the practical meaning of the parameter compensation in part 1: sets with the same fit quality on one experiment behave differently on another.

Recalibrating on the seven baseline runs at 1.2, 3.4 and 5.6 m, with the wall parameters free because the 1.2 m corridors exercise them, was done from two optimizer starts, and again they disagree. Set C has a three-seed norm of 11.3 on its calibration runs and empties 59 of 63 seeds over all conditions. Set D has a norm of 16.2 and stalls in 14 seeds. Set C’s baseline flows are 12 % high on average, with individual runs from 2 % low to 30 % high. Its radius is 0.124 m; set D puts the radius almost at its 0.10 m lower bound. A second optimizer start changed both the inferred parameters and whether the narrow gate remained passable.

CrowdQueue: experiment (bars) against Hermes set B transferred without fitting (blue), the CrowdQueue calibration, set C (orange), and the joint calibration (green), one marker per seed, columns ordered by corridor width with the run number; asterisks mark the runs used in the calibration. Filled circle: the run emptied within 120 s; open circle: not emptied, discharge continuing; open square: stalled. Top: baseline motivation, bottom: low motivation.

The low-motivation runs show what a flow-only validation would miss. Between the baseline and the low-motivation condition the measured flow at the same gate is about 10 % lower; the model has no input for motivation, so with parameters fixed at the baseline calibration it predicts the same or a higher flow, and overshoots the low-motivation runs by 24 % on average. In the 1.2 m corridor with 24 people, set C gives a flow 29 % high, a front area 44 % denser than measured and people moving 21 % slower: throughput alone hides two compensating errors in the upstream state.

What static parameters can absorb

Can a different parameter set stand in for motivation? Freeing all seven parameters is a poor diagnostic because every interaction can compensate for every other one. We therefore freeze the five interaction parameters at set C and profile only desired speed v0 and time gap T, separately for the baseline condition h0 and the low-motivation condition h−. These profiles are fitted and evaluated on the same runs, with one-seed grids and Dakota runs that stopped at the iteration cap with their convergence criteria unmet. They are in-sample diagnostics, not a validated motivation model.

For h0, Dakota selects (v0, T) = (1.17 m/s, 0.79 s) with a three-seed norm of 9.19 over 33 residuals. For h− it selects (1.50, 1.73) with 11.74. Nearby low cells follow a diagonal valley in each grid, so the pair is constrained jointly more clearly than either coordinate is identified alone, the same degeneracy between desired speed and headway that Kretz et al. (2026) describe for the social force model. Neither profile point transfers to the other condition: the h0 point scores 15.65 on h− data and the h− point 13.95 on h0 data. The first profile had stopped at a bound of T = 1.2 s and appeared to show that no pair could represent h−; widening the box to 2.0 s removed that conclusion. A bound hit is a prompt for a sensitivity check, not a result.

The larger time gap reduces the residuals in this conditional slice. At the h− profile point, flow is within 13 % in every run and density within 11 % in seven of nine. It pays for that with speed: seven runs are 24 to 53 % too slow. Part of that residual is the observable rather than the model. The speed is a PedPy individual speed over a 0.4 s window, and people in the front area who move less than 10 cm in 2 s still register 0.09 to 0.13 m/s, a standing-frame floor that plausibly comes from head tracking and body sway. In the 4.5 and 5.6 m corridors the measured speeds of 0.12 to 0.21 m/s are within a factor of two of that floor, so a 10 % tolerance there is tighter than the observable resolves. We have not built an error model for it: the floor is an empirical warning, not a distribution, and we do not know how it varies with density or width or how it covaries with the density measurement. A post-hoc check with the radial speed towards the gate lowers the h− norm from 11.7 to 9.8 and leaves five runs more than two assumed standard deviations too slow. The static v0–T pair absorbs part of motivation, but it does not reproduce flow, density and speed together across corridor widths, and the speed residuals near the floor are poorly resolved.

The densities become tangible as head counts because the front area is 2.2 m². In the two 1.2 m pairs, the experiment averages 6.1 people for h0 against 3.9 for h− with 24 participants, and 12.2 against 5.8 with 63. With one common set C, simulation compresses those contrasts to 5.9 against 5.7 and 9.7 against 7.2. After the separate v0–T fits, the means become 5.6 against 5.0 and 12.3 against 5.6. The simulations replay each run’s observed initial positions and arrivals, so this is conditional evidence: it tests whether the model preserves the supplied positioning difference, not whether it generates that difference from motivation.

People inside the 2.2 m² area immediately upstream of the gate, experiment against three simulation seeds at the condition-specific v0–T profile points. The paired runs have the same corridor width and nominal population. Each simulation is conditioned on that run’s observed initial positions and arrivals; the figure is not a generative test of motivation.

The two 63-person runs in the 1.2 m corridor, tracked experiment beside one simulation seed at the condition-specific v0–T profile point. The dashed rectangle is the front measurement area; the counts below each panel are the people inside it and the people who have passed the gate. Illustration only: the simulation starts from the observed positions, so the sparser low-motivation queue is partly inherited.

We also tested a spacing explanation by fixing v0 and T at the h0 profile point and varying radius and neighbor range. Its best point, radius 0.150 m and range 0.083 m, has a three-seed norm of 14.4, worse than the timing profile, and most of the surface is not a fit landscape but a stall cliff: from a neighbor range of 0.337 m upward all nine runs stall or fail at every radius, because radius and range govern both upstream separation and access through the same 0.5 m gate. This negative slice does not by itself identify the missing state variable.

One-seed residual profiles with all unshown parameters fixed. Left and centre: desired speed v0 against time gap T for h0 and h−; right: radius against neighbor range for h− with v0 and T fixed at the h0 profile point. Open circles mark grid minima and white crosses the set-C coordinates. On the spacing panel, black dots mark grid points with some stalled or failed runs and black crosses points where all nine were stalled or failed. The common colour scale is clipped at norm 30.

The one high-motivation run, 11 people through the 1.2 m corridor at 2.1 persons per second, gets its own sweep over the time gap. A time gap near 0.3 s reaches the measured flow; there the density is about three times the measurement and the speed 30 % low. This condition has one usable run, so the sweep is a diagnostic, not a calibration.

The high-motivation run: flow, density and occupied-frame speed in the corridor against the time gap, all other parameters at set C, three seeds per point shown individually. Filled markers and the line: runs that emptied. Grey band: the measured value with the uncertainty assumed in the calibration weights.

With the stated tolerance, recalibration on CrowdQueue makes passage reliable without bringing every observable within tolerance, and parameter sets that fit the wide bottleneck equally well behave very differently at the narrow gate. Why the remaining discrepancies are there, whether the round body, the interaction rules or something else, these tests do not say.

Step 7 — a third experiment, without barriers

The BaSiGo entrance experiment, Düsseldorf 2013: the crowd in front of the unguided entrance, top right. Image: Pedestrian Dynamics Data Archive, Forschungszentrum Jülich.

The last test uses the BaSiGo entrance experiment from 2013 (Sieben et al. 2017): an entrance with two half-metre lanes and no guiding barriers, 319 people told that their favourite artist is playing and they want to be first in, 273 of them tracked. Nothing was fitted to it, and the comparison is a spatial profile rather than a number, the kind of comparison Liao et al. (2017) asked for, here conditional on the observed arrivals. The simulation injects each of the 273 tracked agents at the time and place where the data first see them. The 46 participants who were never tracked are absent, so the simulated crowd is about 15 % smaller than the real one, which changes the interactions in ways we did not quantify.

Time-averaged density in front of the entrance, 20 to 110 s into the run, on a 0.5 m grid: the experiment and each simulation seed separately, with its entrance flow over the same window and the number of people who passed within the 120 s run. Grey: the entrance barriers. Pale yellow is zero density; white lies outside the mapped grid.

Hermes set A produces a semicircle, but the wrong one. Its density is concentrated at the entrance and falls off regularly, with peaks around 6 per square metre; the measured crowd is dense to 10 per square metre over a region that extends two to three metres upstream and is visibly asymmetric, heavier to the left of the entrance. The model drains the crowd at 1.3 persons per second where the people, pushing to be first, achieved 0.57; the joint parameters give 1.36 to 1.39 persons per second, with the same regular shape. What the model lacks is a hypothesis, not a finding of this figure: we suspect the pressure of a crowd that wants to be first, and the shoulder rotation that gets real people through half a metre. The maps only show that the shape, the extent and the flow are all wrong, in the same direction for both parameter sets.

The joint calibration

The last set in the table comes from one Dakota run over both experiments: the three Hermes calibration widths and the seven CrowdQueue baseline runs, 30 residuals in all, each divided by its assumed uncertainty and entering with equal weight, so CrowdQueue carries about 70 % of the objective by count. One seed per run, one optimizer start, efficient_global with at most 40 iterations; both stopping criteria were met after 17. The result is compared with the specialists on post-hoc residual norms from stored seed means and on how many CrowdQueue seeds empty:

parameter set	Hermes norm (9 residuals)	CrowdQueue norm (21 residuals)	CrowdQueue seeds emptied / stalled / otherwise not emptied
Hermes set A	2.7	30.0	16 / 46 / 1
Hermes set B	2.4	14.9	59 / 0 / 4
CrowdQueue set C		11.3	59 / 0 / 4
CrowdQueue set D		16.2	47 / 14 / 2
joint	6.2	13.6	57 / 3 / 3

The joint point is worse than the relevant specialist in both experiments: 6.2 against 2.4 to 2.7 on Hermes, and 13.6 against 11.3 on CrowdQueue. The single search found no parameter set that meets the stated tolerances in both experiments. Whether a better compromise exists elsewhere in the parameter space, one start cannot say.

What happened to the calibrated model

It is worth being plain about it. After Hermes we had two parameter sets that reproduced five bottleneck widths mostly within ten percent, neither passing the tolerance everywhere. At a half-metre gate one of them stalled in 46 of 63 seeds and the other was 20 to 40 % too fast. The set calibrated on the gate passes it in every seed in which everyone could be placed, with flows 12 % high on average. Separate timing parameters reproduce much of the front occupancy contrast between baseline and low motivation, in-sample and conditioned on the observed positions, but not the speed. At the unguided entrance every set tested drains the crowd two to three times faster than the people did. The one parameter that was actually measured, the free speed, was the one the first optimizer run most wanted to change.

None of this is a verdict on the Collision Free Speed model in particular. It is what cross-regime transfer testing looks like, as opposed to a fit within one regime. Other pedestrian models with a handful of scalar parameters and round agents may show similar limits, but this note tests one model on three datasets and claims no more.

Lessons

A calibration is a statement about a regime, not about a model. Parameters fitted to wide bottlenecks are parameters for wide bottlenecks. Using them for a turnstile, a narrow door or an unguided entrance is an extrapolation, and here it was off not by percent but by clogging or not clogging. Report the experiments a parameter set was fitted to and tested on; without that, “calibrated” is not information.
State the purpose, then pick the observables. Flow alone would have passed a model with the wrong density, as it did in 2014. Analyse experiment and simulation with the same code, so the observables are defined once.
Measured quantities are data, not fitting variables. Fixing the free speed at its measured value removed one unphysical value from the fit. It did not make the remaining parameters physical, and the value itself was measured on one experiment’s participants and assumed for the others.
Run the optimizer more than once. Every calibration in this note that was started twice gave two different parameter sets. Multiple starts are a robustness check, not a map of the landscape, and a set reported without one is a sample.
Screen before you calibrate. The Morris run cost 160 evaluations, removed two of seven parameters, and mapped a region where the model fails outright. Use sensitivity analysis for the grouping of parameters unless you have paid for converged indices.
Measure an observable’s noise floor before you weight it. People standing still register about 0.1 m/s with the speed definition we calibrated against, and the wide-corridor speeds are within a factor of two of that. Estimate the floor first, and do not read residuals near it as well resolved.
A fitted parameter shift is not a mechanism. The timing profiles absorb part of low motivation and the spacing slice develops a stall cliff. Both are useful diagnostics; neither says what motivation is. Keep parameter uncertainty and seed variability separate, and label as exploratory whatever was chosen after looking at the data.
Audit the pipeline before you audit the model. Reviews of earlier drafts found late participants missing from simulations, unwritten final trajectory frames, a route around rather than through the gate, empty frames averaged as zero speed, and one script that ran set D while its output was discussed as set C. Each was plausible, produced numbers, and changed a conclusion. Each was corrected and every dependent result recomputed. The drivers now record an evaluation status file for every simulation, and the pipeline corrections section of the repository README lists each error and its fix, pinned to the commit this note was built from.

Several limitations remain open: one realization per condition, correlated observables, pragmatic weights, one optimizer start for the joint search, no time-step convergence study and no quantified observation error. The next useful work is not a uniquely calibrated parameter set but independent repetitions, more spatial and temporal observables, an explicit motivation state, a pre-specified uncertainty model and verification tests independent of the calibration data. Until then the conclusion is conditional: these scripts demonstrate a reproducible calibration and validation workflow and expose regime failures. They do not establish the validity of the model.

A note on open-source software

Sargent drew the cost of validation against the confidence it buys as a curve that rises steeply near the end. Everything in this study is about pushing that curve down. The trajectories are on a public archive with a DOI and the geometry attached, the model is open source and scriptable, the analysis library is the same one used on the experiment, and the calibration toolkit has been maintained by a national laboratory for 25 years. Anyone who thinks a remaining discrepancy is a model deficiency can test it this afternoon. And a validation that lives in a text file can be reviewed, which a sentence like “parameters were chosen according to the literature” never could.

References

Adrian, J., Seyfried, A., Sieben, A. (2020). Crowds in front of bottlenecks at entrances from the perspective of physics and social psychology. Journal of the Royal Society Interface 17, 20190871. doi:10.1098/rsif.2019.0871.
Kretz, T., Eckes, L., Knappik, L., Diz, J., Lipp, D., Müller, S. (2026). Using empirical travel time distributions for calibration of a model of pedestrian dynamics. EURO Journal on Transportation and Logistics 15, 100179. doi:10.1016/j.ejtl.2026.100179.
Liao, W., Chraibi, M., Seyfried, A., Zhang, J., Zheng, X., Zhao, Y. (2014). Validation of FDS+Evac for pedestrian simulations in wide bottlenecks. IEEE ITSC 2014, 554–559. doi:10.1109/ITSC.2014.6957748.
Liao, W., Zhang, J., Zheng, X., Zhao, Y. (2017). A generalized validation procedure for pedestrian models. Simulation Modelling Practice and Theory 77, 20–31. doi:10.1016/j.simpat.2017.05.002.
Sargent, R. G. (1984). Simulation model validation. In: Simulation and Model-Based Methodologies: An Integrative View, Springer, 537–555.
Sieben, A., Schumann, J., Seyfried, A. (2017). Collective phenomena in crowds — where pedestrian dynamics need social psychology. PLOS ONE 12(6), e0177328. doi:10.1371/journal.pone.0177328.
Experiment data: Hermes bottleneck experiment, Düsseldorf 2009, doi:10.34735/ped.2009.6; CrowdQueue experiment, Wuppertal 2018, doi:10.34735/ped.2018.1; BaSiGo entrance experiment, Düsseldorf 2013, doi:10.34735/ped.2013.2. Dakota: Adams et al., Sandia National Laboratories, version 6.24.

Code

Driver scripts, Dakota input files, results and figures for every step: github.com/PedestrianDynamics/jupedsim-dakota-calibration

Software used in this study:

Dakota 6.24
JuPedSim 1.4.2
PedPy 1.4.0

By: Mohcine Chraibi

Calibrating a pedestrian model against real data, part 1: the bottleneck

Sat, 05 Sep 2026 00:00:00 +0000

Abstract. This is the first of two parts. We calibrated JuPedSim’s Collision Free Speed model with Dakota against the Hermes bottleneck experiment, three widths in the fit and two held out, with the same open analysis code on experiment and simulation. The defaults miss the measured flow by a factor of two. After calibration, flow, density and speed in front of the opening are mostly within ten percent, including at the two held-out widths, but not everywhere: both calibrated sets overshoot the flow at 3.0 m and undershoot it at 5.0 m. Two optimizer starts give two parameter sets with near-equal residual norms that differ by factors of two to four. The held-out widths are an interpolation test inside one experiment, not an independent validation. Part 2 takes the same parameters to a half-metre gate and an unguided entrance, where they fail.

Pedestrian simulations are used to size exits, plan events and argue about safety, so the question whether a model is right is not academic. The simulation community settled the vocabulary long ago. Verification asks whether the software implements the model correctly. Validation asks whether the model, within the domain it is meant for, reproduces reality with an accuracy that is good enough for its purpose (Sargent 1984, 2008; ISO 16730 as summarised by Ronchi et al. 2013). Sargent adds two points that are easy to forget: a model is never valid in the abstract but only for a purpose, and confidence in a model costs money, so one always stops somewhere short of certainty.

Two lessons from this field’s own validation work shape what follows. Liao et al. (2014) calibrated FDS+Evac on the very experiment we use here, by hand, to match the flow at one width. The flow matched. The density in front of the bottleneck was too low and the speed too high, so the model reached the right flow for the wrong reasons. Their follow-up (Liao et al. 2017) turned that into a rule: a single characteristic cannot validate a model, and the comparison has to combine several. Kurtc et al. (2018) drew the second lesson: the assessment should be automated, including the search for parameters, or it will not be done consistently. Kretz et al. (2026) reached the first lesson from the other side with a social force model: a single-parameter fit to the mean flow matched the flow either way, and only the individual travel times told a good calibration from a bad one.

This note puts the two lessons together with tools that are open and generic. The simulation is JuPedSim, the analysis of experiment and simulation alike is PedPy, and the screening, sensitivity analysis and calibration are run by Dakota, Sandia’s toolkit for exactly this kind of study. Dakota has been used in engineering for decades and almost never in pedestrian dynamics. It turns “which parameters matter, and what values fit the data” from a pile of private scripts into a text file.

What validation means here

Following Sargent, we state the purpose first. The model should reproduce the flow through wide bottlenecks in a dense, motivated crowd, and it should do so with the right density and speed in front of the bottleneck, not only the right flow. Three observables therefore enter every comparison: throughput and the two components of the upstream state. They are the complementary quantities we chose for this study, not a canonical minimum. The decision rule is a fixed relative tolerance per observable, 6 % on flow and 10 % on density and speed, and a comparison passes when the mean over simulation seeds lies within it. The same numbers weight the calibration residuals. They were set before any calibration was run, from the seed-to-seed spread of the simulation plus an allowance for measurement error, but they were chosen as calibration weights. Adopting them as the acceptance tolerance is a decision we made when the results were assessed, not a pre-registered one. Two of the five bottleneck widths are held out: the optimizer sees 2.4, 3.6 and 5.0 m and never 3.0 and 4.4 m. That is an interpolation test within one apparatus and one participant pool. It is a stronger check than a fit, and a weaker one than a second experiment.

Step 1 — the experiment

The Hermes bottleneck experiment, Düsseldorf 2009: overhead view of a run, and the setup sketch with the semicircular holding zones at 3 persons per square metre. Images: Pedestrian Dynamics Data Archive, Forschungszentrum Jülich.

The Hermes bottleneck experiment was run in May 2009 in Hall 2 of the Düsseldorf fairground with about 350 participants (Seyfried et al. 2009; Liao et al. 2014). The bottleneck was built from boards higher than two metres, 1 m long, and its width was varied from 2.4 to 5.0 m in five runs, one run per width. Participants waited in a semicircular holding area of radius 8.6 m directly in front of the bottleneck, which puts the initial density at three persons per square metre, and walked through on command. The free walking speed of 42 participants was measured separately: 1.55 m/s with a spread of 0.18 m/s.

The trajectories are open data with a DOI. Since 2024 they come as HDF5 files that carry a geometry along, so PedPy loads both with one call each. One caveat, found only by reading the paper next to the file: the polygon in the archive is a 12 m box around the camera window with the wall blocks cut off at its edges, and three of its five gaps are 0.1 m narrower than the run’s stated width. The simulation therefore takes the gap width from the run parameter and the corridor and holding area from the paper. A corrected polygon for each run is in the repository. The initial positions are not the measured ones either: the simulation places 350 agents on a jittered lattice inside the semicircle at the stated density, so the holding area is an initialization region, not a reconstructed crowd.

import pedpy, pathlib
f = pathlib.Path("ao-360-400.h5")
traj = pedpy.load_trajectory_from_ped_data_archive_hdf5(f)
area = pedpy.load_walkable_area_from_ped_data_archive_hdf5(f)

We measure three things per run: the flow through the bottleneck, from the slope of the N(t) curve at a line across the gap; and the density and mean speed in a 2.8 by 2 m area directly in front of it, averaged over the jam phase. The area is the same for all widths, so it covers more than the 2.4 m opening and only the middle of the 5.0 m one. These are local measurements of the crowd state upstream of the gap, and the calibration compares them like for like. The same three functions are applied to the simulated trajectories. Seyfried and Schadschneider’s (2008) warning that the measurement method changes the result is answered the simplest way: there is only one method, applied to both sides.

Measurement setup for the 2.4 m and 5.0 m runs. Red: the flow line across the gap. Blue: the area for density and speed. The dashed line is the simulation geometry: a 20 m wide corridor with the bottleneck boards as a wall band. The dotted lines mark the camera window of the experiment.

Step 2 — the defaults are wrong by a factor of two

We rebuilt the setup in JuPedSim: the corridor, the boards, the semicircular holding area with 350 agents at three per square metre, and an exit line behind the bottleneck. The model is the Collision Free Speed model (Tordeux et al. 2016) with its default parameters.

Measured N(t) curves (solid) against the default model (dashed), and flow against bottleneck width. The defaults underpredict the flow by a factor of two to three, worst at the widest opening.

This is not a surprise if you know the model. In the Collision Free Speed model the time gap sets the headway an agent keeps to the one in front, and the radius sets how close agents can pack side by side; together with the neighbor repulsion they fix the throughput of a jammed opening, and the desired speed hardly enters once the crowd is jammed. With the defaults, one second of headway and a 0.2 m radius, the model delivers a third to a half of the measured flow. No amount of tweaking the desired speed fixes that. The question is which of the seven parameters do, and by how much.

Step 3 — let Dakota find the parameters that matter

Dakota needs two things: an input file describing the study, and a driver script that turns a parameter vector into responses. Our driver reads Dakota’s parameter file, runs one JuPedSim simulation per bottleneck width, computes the three observables with PedPy and writes them back. That is all the coupling there is.

method
psuade_moat
samples = 160
partitions = 3
variables
continuous_design = 7
descriptors 'desired_speed' 'radius' 'time_gap' 'strength_neighbor'
'range_neighbor' 'strength_geometry' 'range_geometry'
lower_bounds 0.8 0.12 0.10 2.0 0.02 1.0 0.01
upper_bounds 1.8 0.25 1.20 15.0 0.50 10.0 0.20
interface
fork
analysis_drivers = 'python3 driver.py'
parameters_file = 'params.in'
results_file = 'results.out'
work_directory named 'runs/run' directory_tag
asynchronous evaluation_concurrency = 3
responses
response_functions = 9
no_gradients
no_hessians

The first stage is a Morris screening: twenty trajectories of eight evaluations each, every step changing one parameter at a time, 160 evaluations in all. It is cheap and it ranks the parameters per observable. It is a ranking, not a measurement, and it is conditional on the ranges it was sampled over.

Morris screening over the full parameter ranges, twenty trajectories. Each cell is the mean absolute elementary effect of a parameter on an observable, scaled to the strongest parameter in that column. Steps where either end failed were discarded, and n per row is the number of valid steps out of the twenty drawn. The time gap, the radius and the neighbor repulsion range dominate; the desired speed has a moderate effect on flow and speed over this wide range; the two wall parameters are weakest.

The screening also revealed something no fit would have: in 43 of the 160 evaluations, about a quarter, the model pushed agents through the walls and JuPedSim aborted the run. Those combinations have strong, long-range neighbor repulsion and weak wall repulsion. This is a failure of the model’s numerical domain, not an observation about pedestrians, but it is one an engineer can land in by hand. We narrowed the bounds accordingly and froze the wall parameters at their defaults.

With five parameters left, and narrower ranges, a Sobol analysis gives the variance decomposition. Dakota computes the indices itself; the input file changes the method block to sampling with variance_based_decomp, the variables become uncertain variables with uniform distributions over the narrowed ranges, and the cost is the base sample size times the number of parameters plus two. Each parameter point uses simulation seed 1; seed variability is measured separately later rather than treated as another physical input.

How many samples are enough is an empirical question, so we answered it empirically: three replicate runs with 40 base samples, then 80 and 160, 2520 evaluations in all. The magnitudes are not converged: the time gap’s total index on flow and speed still rises from about 0.4 to about 0.7 between 80 and 160 samples. What is stable across every run is the small group of influential parameters per observable: time gap and radius for flow, radius and neighbor range for density, time gap and neighbor range for speed, with the desired speed a minor contributor. Since only the grouping is used below, we did not spend further computation on the numbers.

Sobol indices for the five remaining parameters from the run with 160 base samples, 1120 evaluations, none failed. The time gap dominates flow and speed with total indices of 0.6 to 0.8; the radius dominates density with 0.5 to 0.6, followed by the neighbor repulsion range with about 0.3; the desired speed contributes 0.06 to flow and 0.12 to 0.17 to speed.

Sobol total indices against the base sample size. Bars at 40 span three replicate runs with different seeds. The group of influential parameters per observable is the same in every run; the magnitudes are not settled at 160 base samples.

Step 4 — calibration, and why we did not use gradients

Our first attempt at calibration, on a simpler two-room test case before touching the experiment, used Dakota’s gradient-based least-squares solver with finite-difference gradients. It stalled at the starting point: two simulations whose parameters differ in the sixth digit produced evacuation times that differed by seconds, and the differences used as gradients were noise. We did not repeat the test on the Hermes objective, so this is an observation about that solver on that test case, not a property of the calibration problem here.

The method we used instead is Dakota’s efficient_global, an implementation of Efficient Global Optimization (Jones et al. 1998). The objective is a weighted sum of squared residuals over the nine observables, three per calibration width,

$$ f(\theta) = \sum_{i=1}^{9} w_i \,\bigl(y_i(\theta) - \hat{y}_i\bigr)^2 , $$

where $\theta$ is the parameter vector, $y_i(\theta)$ the simulated observable, $\hat{y}_i$ the measured one, and $w_i = 1/\sigma_i^2$ with $\sigma_i$ the assumed uncertainty, 6 % of the measured flow and 10 % of the measured density and speed; the weights in the input file below are these inverse squared uncertainties. The residual norm quoted throughout is $\sqrt{f}$, in units of the assumed uncertainty. Each evaluation of $f$ costs three JuPedSim runs at one fixed seed, so the optimizer sees a single realization of a stochastic simulator; seed variability is measured afterwards, not modelled in the fit. Finite differences of such an objective are noise, so the optimizer never differentiates it. Instead it fits a Gaussian process to all evaluations made so far, which gives a prediction $\mu(\theta)$ and a standard deviation $s(\theta)$ at every untried point, and it picks the next point where the expected improvement over the best value $f_{\min}$ seen so far is largest:

$$ \mathrm{EI}(\theta) = \bigl(f_{\min} - \mu(\theta)\bigr)\,\Phi(z) + s(\theta)\,\varphi(z), \qquad z = \frac{f_{\min} - \mu(\theta)}{s(\theta)} , $$

with $\Phi$ and $\varphi$ the standard normal distribution and density functions. The first term rewards points the surrogate expects to be good, the second rewards points where it is uncertain, so the search balances exploiting the current best region against exploring the rest of the box. It stops when the largest expected improvement falls below a threshold, or when the new point is too close to an evaluated one, or at the iteration cap. Each step costs one simulation batch and a surrogate fit, so a calibration takes a few dozen evaluations rather than the hundreds a gradient method would spend on noise.

method
efficient_global
max_iterations = 60
responses
calibration_terms = 9
weights = 7.6 3.2 977 3.4 4.3 517 1.45 6.0 252
no_gradients
no_hessians

After 42 evaluations, about twenty minutes on a laptop, the optimizer converged with all nine residuals below one standard deviation of the assumed uncertainty. And the desired speed had gone to the lower bound of its range, 0.8 m/s, half the free speed measured on the participants.

This is the moment where a fit and a validation part ways. The Sobol analysis had already said that the desired speed is a minor influence on the jam observables, so the misfit is flat along it and the optimizer was free to put it wherever the residuals were a hair smaller. A model with participants walking at 0.8 m/s reproduces the bottleneck. It does not reproduce the participants. Sargent calls this data validity: parameters that were measured are data, not fitting variables. Liao et al. (2014) made the same move for FDS+Evac. So we fixed the desired speed at the measured 1.55 m/s and calibrated the remaining four parameters. That converged in 24 evaluations.

parameter	default	all five free	speed fixed, set A	speed fixed, set B
desired speed [m/s]	1.2	0.80 (at bound)	1.55 (measured)	1.55 (measured)
radius [m]	0.20	0.14	0.13	0.15
time gap [s]	1.0	0.56	0.81	0.55
neighbor repulsion strength	8	6.1	2.1	9.4
neighbor repulsion range [m]	0.10	0.15	0.25	0.10
weighted residual norm (six-seed recalculation)			2.7	2.4

Set B is the reason for the last column. A surrogate-based optimizer with a few dozen evaluations can settle in one region of a flat landscape, so we repeated the fixed-speed calibration from a different optimizer seed. Dakota reported objective norms of 1.0 and 1.5 at the selected evaluations. Recomputing the norm from the six-seed means used for the later comparison gives 2.7 for set A and 2.4 for set B; those are the values in the table. We have no uncertainty on these norms, the observables are correlated, and we did not sample the objective’s seed distribution, so we cannot say that the two sets are statistically equivalent or that the landscape has exactly two optima. What two starts do show is that the reported optimum depends on where the optimizer started. Reported without a second start, a calibrated parameter set is one sample.

Both sets reproduce flow and density at the three calibration widths to within about ten percent, and speed to within twenty, with residual norms we cannot distinguish. They are nowhere near each other. Time gap, repulsion strength and range shift by factors of two to four between them. Three observables at three widths do not pin four free parameters, and a good fit says nothing about whether the parameters mean what their names say. Fixing the one parameter that was measured removes the freedom to be wrong in that direction. It does not make the others physical: the radius that comes out, 0.13 to 0.15 m, is an effective size in a model with round bodies, and the remaining parameters still compensate for one another.

Step 5 — does it hold on the widths it never saw?

The 3.6 m run, tracked experiment inside the camera window beside one seed of the simulation with set A, cropped to the same window. Both clocks start at the first crossing of the flow line (red); the dashed box is the density and speed area. Illustration only: the simulation starts from a synthetic lattice, not from the measured positions.

Flow, density and speed against bottleneck width: experiment with the acceptance band (6 % on flow, 10 % on density and speed, the uncertainty assumed in the calibration weights), default model, and the two calibrated sets with the desired speed fixed. Error bars are the minimum and maximum over six simulation seeds; the experiment is one run per width. The dotted widths, 3.0 and 4.4 m, were not used in the calibration.

With the desired speed fixed, both calibrated sets improve every observable substantially, at the held-out widths too. Applied as a decision rule, the tolerance passes neither set everywhere. Both overshoot the flow at the held-out 3.0 m by 10 to 13 %; set A’s speed is 12 % low at the held-out 4.4 m; and at the 5.0 m calibration width both are 10 % low on flow. Where the two sets differ is what else they get wrong at the wide end. Set A follows the measured increase of speed with width up to 3.6 m and then flattens, 17 % low at 5.0 m, while it keeps the density within tolerance. Set B keeps the speed within tolerance at every width and instead lets the density fall 10 % short at 5.0 m. So at the widest opening set A preserves the local density and underpredicts the speed, set B improves the speed and underpredicts the density, and both underpredict the flow. That is a milder version of what Liao et al. (2014) saw in FDS+Evac.

Deviation of the simulation from the experiment, mean over six seeds, with the band at 6 % for flow and 10 % for density and speed:

width [m]	flow A / B	density A / B	speed A / B
2.4	+3 % / +5 %	+7 % / +1 %	−6 % / −5 %
3.0 (held out)	+10 % / +13 %	+2 % / −3 %	+5 % / +9 %
3.6	+3 % / +4 %	0 % / −6 %	−3 % / +4 %
4.4 (held out)	+2 % / +4 %	+8 % / −1 %	−12 % / +1 %
5.0	−11 % / −10 %	−1 % / −10 %	−17 % / −4 %

The seed-to-seed standard deviation is 2 to 4 % for flow and density and 3 to 6 % for speed, so the differences between A and B at the wide openings are larger than the simulation noise. The experiment, however, is one run per width, and we have not quantified its measurement error; the band is an assumption, not a measured interval.

The honest summary of Step 5 is: a substantial improvement over the defaults, and a partial transfer to the held-out widths, but not a model validated for the full 2.4 to 5.0 m range by the stated rule. At the calibration widths 2.4 and 3.6 m both sets meet the tolerance on all three observables, which is a fit, not validation. The held-out evidence is two widths: at 4.4 m set B meets the tolerance on all three observables and set A fails it on speed; at 3.0 m both fail it on flow. And both miss the widest calibration opening in a way the two calibrations move between speed and density but neither removes. Whether that is acceptable depends on the purpose, which is why the purpose has to be stated first.

What part 1 establishes, and what it does not

A reproducible calibration. Every step is a Dakota input file and a driver script in the repository. Anyone can rerun it with another observable, model or experiment from the same archive.
The defaults are not usable for this purpose. They miss the flow by a factor of two to three, and no single parameter fixes it.
Two parameter sets, not one. Two optimizer starts give sets with near-equal residual norms that differ by factors of two to four. The data do not identify the parameters; they constrain combinations of them.
Partial width transfer. Held-out widths inside one experiment are an interpolation test. The model passes it in part. This is not evidence for any other geometry, crowd or motivation.
A pipeline audit. Between the first draft and this one, reviews found participants missing from simulations, unwritten final trajectory frames and a geometry error. Each was corrected and every dependent result recomputed; the Hermes steady-phase observables changed by less than half a percent. The record is the pipeline corrections section of the repository README, pinned to the commit this note was built from.

The natural next question is what the two sets do outside this experiment. Part 2 takes both to a half-metre gate and to an unguided entrance. One of them stalls, the other is too fast, and a second calibration on the new data disagrees with the first about which parameters matter.

References

Jones, D. R., Schonlau, M., Welch, W. J. (1998). Efficient global optimization of expensive black-box functions. Journal of Global Optimization 13, 455–492. doi:10.1023/A:1008306431147.
Kretz, T., Eckes, L., Knappik, L., Diz, J., Lipp, D., Müller, S. (2026). Using empirical travel time distributions for calibration of a model of pedestrian dynamics. EURO Journal on Transportation and Logistics 15, 100179. doi:10.1016/j.ejtl.2026.100179.
Kurtc, V., Chraibi, M., Tordeux, A. (2018). Automated quality assessment of space-continuous models for pedestrian dynamics. arXiv:1809.01862.
Liao, W., Chraibi, M., Seyfried, A., Zhang, J., Zheng, X., Zhao, Y. (2014). Validation of FDS+Evac for pedestrian simulations in wide bottlenecks. IEEE ITSC 2014, 554–559. doi:10.1109/ITSC.2014.6957748.
Liao, W., Zhang, J., Zheng, X., Zhao, Y. (2017). A generalized validation procedure for pedestrian models. Simulation Modelling Practice and Theory 77, 20–31. doi:10.1016/j.simpat.2017.05.002.
Ronchi, E., Kuligowski, E. D., Reneke, P. A., Peacock, R. D., Nilsson, D. (2013). The process of verification and validation of building fire evacuation models. NIST Technical Note 1822.
Sargent, R. G. (1984). Simulation model validation. In: Simulation and Model-Based Methodologies: An Integrative View, Springer, 537–555.
Sargent, R. G. (2008). Verification and validation of simulation models. Proceedings of the Winter Simulation Conference, 157–169.
Seyfried, A., Schadschneider, A. (2008). Fundamental diagram and validation of crowd models. ACRI 2008, LNCS 5191, 563–566.
Seyfried, A., Passon, O., Steffen, B., Boltes, M., Rupprecht, T., Klingsch, W. (2009). New insights into pedestrian flow through bottlenecks. Transportation Science 43, 395–406.
Tordeux, A., Chraibi, M., Seyfried, A. (2016). Collision-free speed model for pedestrian dynamics. Traffic and Granular Flow ‘15, 225–232.
Experiment data: Hermes bottleneck experiment, Düsseldorf 2009, doi:10.34735/ped.2009.6. Dakota: Adams et al., Sandia National Laboratories, version 6.24.

Code

Driver scripts, Dakota input files, results and figures for every step: github.com/PedestrianDynamics/jupedsim-dakota-calibration

Software used in this study:

Dakota 6.24
JuPedSim 1.4.2
PedPy 1.4.0

By: Mohcine Chraibi

Why do roaming crowds drift counterclockwise — and where?

Mon, 29 Jun 2026 00:00:00 +0000

Recently a paper in Nature Communications (Echeverría-Huart et al.) reports something simple and odd. If you put people in a circular arena and let them walk around freely, the crowd slowly drifts counterclockwise. The authors saw it in experiments in Spain and Japan, so it is not a fluke of one room or one group. Their explanation is that the drift does not come from people avoiding each other. It comes from each person, on their own, having a slight tendency to turn left. (The paper even made it into The New York Times.)

This sounds intriguing, so we explored it a bit further.

We used JuPedSim. It splits the problem into two parts. One part is how people avoid bumping into each other, which JuPedSim handles with a couple of operational models. The other part is where each person wants to go next, which is up to us. We give each agent a desired direction; JuPedSim returns motion without collisions. That split is convenient here. We drive each agent ourselves with direct steering: every step, we pick a heading and place the agent’s target just ahead of it.

sim = jps.Simulation(model=jps.SocialForceModel(), geometry=disk)
steering = sim.add_direct_steering_stage()
journey = sim.add_journey(jps.JourneyDescription([steering]))
while sim.agent_count() > 0:
 for agent in sim.agents():
 heading = roam(agent, cfg)
 agent.target = ahead(agent.position, heading)
 sim.iterate()

The for loop is our little model. sim.iterate() is JuPedSim’s machinery. Here is one agent following the target we hand it each step:

One roaming agent (blue) and the target we set each step (red ring, dashed line). With a single agent there is nothing to avoid, so this shows only the steering loop.

The collision model is set on the first line, so we can later run the exact same behavior through different JuPedSim models. This also gives us a clean test of the paper’s claim. The collision handling has no built-in left or right preference. So if a crowd with symmetric rules does not rotate, that supports the idea that the bias is individual. And if adding the bias makes the crowd rotate, that supports the proposed mechanism. We model a 5 m arena with agents roaming in random directions. We measure the same quantity as the paper: M, each person’s velocity projected onto the counterclockwise direction, averaged over the crowd. M̄ > 0 means the crowd circulates counterclockwise.

Symmetric avoidance does not rotate

With no bias, the crowd does not rotate. Across crowd sizes and random seeds, M̄ stays near zero. The collision model on its own picks no direction.

This observation is model agnostic. JuPedSim ships several collision models, some velocity-based and some force-based, and they work quite differently inside. So we ran the same no-bias case through four of them: Social Force, WarpDriver, Collision-Free Speed, and Anticipation Velocity. All four stay within a few hundredths of zero, far below the experimental value. No rotation comes out of symmetric avoidance, whatever model you use. It has to be added.

One real experimental run next to the four model controls. The real crowd circulates counterclockwise; the bare models just mill around.

The no-bias case for four collision models, against the experimental M̄ (dashed line). All four sit near zero.

Adding the bias

The paper’s main claim is an always-on individual bias: a slight per-person tendency to turn left, present at every step, even when walking alone with no walls. The specific “turn left when facing a wall” rule comes from earlier work the authors build on, and is one concrete way that bias could show up. We tried both.

We first tried putting the bias in free-space curvature, a constant gentle left turn at every step, which is the closest simple stand-in for the always-on bias. On its own that did not give the confined rotation (more on why below). So we also tried the wall-only rule, and bingo, that reproduced the effect.

How strong the rotation gets depends on how many people have the bias. A single turn strength saturates, so the useful knob is the fraction of people who turn left, which also matches the paper’s mixed population. With nobody biased the crowd does not rotate. With everybody biased it rotates strongly. With about a third of the crowd turning left, we get M̄ ≈ 0.2, the experimental value.

Mean polarization against the share of left-turners. Zero with no bias, rising with the fraction; about a third matches the experiment.

The distribution of M. The control sits on zero; adding left-turners shifts it counterclockwise, as in the paper.

Three crowds side by side (0%, ~45%, 100% left-turners). The control mills without a net direction; more left-turners make the crowd circulate.

Matching the number is not matching the mechanism

The paper’s data let us go past the single number M̄. Their trajectory files include a per-agent polarization value, and recomputing it with our metric matches theirs exactly. So we can compare not just the average rotation but where in the arena it happens. This was the most useful test, because it can disagree with us instead of just confirming us. (Shout out to the authors for sharing their data on Zenodo!)

The experiment shows a rotation that fills the disk. It is positive at every radius, weak at the center, strongest in the outer-middle, then easing right at the wall. We compared this against two simple models: the turn-left-at-wall rule above, and an “intrinsic” version where every agent has a small constant left veer at every step, which is closer to a bias that is present even away from walls.

Neither matches the experiment, and they fail in different ways. The wall-turn model gets the total amount of rotation right but puts it in the wrong place: a thin spike at the rim and an almost still interior, because the bias only acts at the wall. The intrinsic veer spreads rotation into the interior, closer to where the experiment peaks, but it comes out clockwise, the wrong sign. In open space a left veer does turn counterclockwise, and we checked that directly. But once the arena is closed, the way the veer meets the wall flips the net direction.

Local rotation against distance from the center, for the experiment and the two models. The experiment is positive throughout and peaks in the outer-middle; the wall-turn model spikes at the rim; the intrinsic veer goes negative.

The same three as spatial maps (blue counterclockwise, red clockwise). The experiment is mostly blue across the disk; the wall-turn model is a blue ring; the intrinsic veer is a red band.

We read this as support for the paper. It seems that matching one number with a single knob is easy. The spatial pattern of the real rotation does not fall out of either simple rule, so more work needs to be done on the model.

What we learned, and what’s open

We set out to test the mechanism behind a real effect, and we learned a few concrete things. Symmetric collision avoidance produces no preferred rotation. A turn-left-at-wall rule recovers counterclockwise motion at roughly the reported strength. And looking at where the rotation lives shows that matching the average is the easy part, while matching the spatial structure is far from being a trivial task.

There is more to do. The experiment also reports counterclockwise motion with no walls at all, and for people walking alone, which neither of our simple models covers. Reproducing the full spatial field will need a richer individual model than one knob. That leaves the real open question: with no wall to turn at and no crowd to follow, what makes a single person drift counterclockwise? The authors checked the obvious candidates — handedness, footedness, eye dominance — and ruled them all out, calling the cause an individual bias that is “most likely biologically rooted” and leaving the origin open. So this is a known open problem, worth further research.

A last thought about the tool. None of this is a simulator telling us the right answer, or even ruling out the wrong ones. It does not work that way, and no simulator does. But what it gave us was a cheap way to ask: what if avoidance is symmetric, what if a third of people turn left, what if the bias is at the wall instead of everywhere. Each one is a few lines in the roam() rule, run, and looked at. We did not find the mechanism here, but we did sharpen the question and make a couple of stories look more or less likely — and that is most of what a simulation is good for.

Code and credits

Code, full-quality videos, and data: github.com/PedestrianDynamics/clockwise

Credit for the phenomenon and the experiments goes to Iñaki Echeverría Huarte, Claudio Feliciani, Iker Zuriguel Ballaz and colleagues. There was a charming talk at TGF26 by Iker Zuriguel Ballaz — here is the abstract.

By: Mohcine Chraibi

Why don't airlines board planes the optimal way?

Fri, 12 Jun 2026 00:00:00 +0000

In 2008, astrophysicist Jason Steffen worked out the mathematically optimal way to board an airplane. It boards passengers in a precise interleaved order — window seats first, every other row — so that many people stow their luggage at the same time and nobody ever has to climb over a seated neighbor. In ideal conditions it’s dramatically faster than how airlines actually do it. A counter-intuitive part of the result is that boarding back-to-front, which many airlines use, is no better than boarding everyone at random, and front-to-back is worse than random.

No airline uses Steffen’s method. We rebuilt the classic studies in simulation to look at why.

Step 1 — reproduce the classic ranking

We modelled a single-aisle 180-seat cabin and ran six boarding strategies. Each passenger walks the aisle to their row, holds there while they stow luggage and let neighbors shuffle past (the two things that actually cause boarding delays), then sits. Same geometry, same luggage draws, 20 repetitions per method.

Six methods boarding side by side; the Steffen variants pull ahead while front-to-back jams.

The ranking comes out as in the literature: Steffen fastest, back-to-front barely better than random, front-to-back worst.

Boarding time for each method, uniform passengers, 180-seat cabin, 20 runs per method. Each box shows the spread across runs. Steffen-Perfect is fastest and Front-to-Back slowest; Back-to-Front is close to Random.

This reproduces established results. The next question is the one that matters in practice.

Step 2 — why the optimum stays on paper

Steffen’s method assumes 180 strangers will form one perfect single-file queue and board in a strict choreographed order. Real passengers don’t. Families board together. People show up late. Not everyone follows the plan.

A very recent paper — Dong, Yanagisawa & Nishinari (2025), Physica A, on boarding the future blended-wing-body aircraft — quantified one of these effects: a compliance rate, the share of passengers who board in their assigned slot. We reproduced their compliance sweep on our single-aisle cabin.

The same method at 100% vs. 50% compliance. The red dots are passengers boarding out of their assigned slot; the 50% panel finishes later.

The result matches their paper: as compliance drops, the optimized methods lose their advantage and move toward random boarding, while random itself barely changes. At zero compliance every method gives the same time, because there is no longer any order to follow.

Mean boarding time as compliance falls from 100% (left) to 0% (right), one line per method. As fewer passengers board in their assigned slot, the optimized methods move toward random boarding; Random stays roughly flat and Front-to-Back improves.

Two completely different models — their discrete cellular automaton on a wide multi-aisle aircraft, our continuous pedestrian simulation on a narrow-body — agree on the trend.

Step 3 — the optimum is the fragile one

We pushed the same idea two more ways. We gave passengers realistic profiles (fast young travelers, heavy luggage, elderly with reduced mobility, families with kids), and we let travel groups board together instead of in perfect order.

Boarding at 80% travel groups; each color is one party boarding together, which clumps the order the Steffen method relies on.

In each case, Steffen’s “perfect” method is the most fragile. It is the most sensitive to a mixed passenger crowd, it is the only method that gets slower as more people travel in groups, and it loses its advantage when people do not comply. A coarser, more practical Steffen-style variant, which real passengers could plausibly follow, stays robust and overtakes the optimum.

Boarding time as travel groups grow; the optimum rises while the practical variant holds flat.

The practical takeaway is that the benefit of an optimized boarding order depends almost entirely on passengers following it.

Code and credits

Steffen showed the ranking in 2008 and confirmed it experimentally in 2012; Dong et al. did the robustness analysis in 2025.

Code, high-quality videos & data: github.com/PedestrianDynamics/boarding

By: Mohcine Chraibi

derRSE 2026

Tue, 03 Mar 2026 00:00:00 +0000

The sixth German Conference on Research Software Engineering (derRSE26) will be held at the University of Stuttgart from March 3rd to 5th, 2026.

Talk

JuPedSim: From Obscure Executables to a Publicly Documented Library

Speaker: Mr Kai Kratz

JuPedSim is an open-source framework for simulating how pedestrians move through buildings, public spaces, and city environments. For years, it was one of those research tools that “just worked” — if you knew how to use it. The software has been part of the pedestrian dynamics community for over 15 years, supporting research and practice from evacuation planning to urban mobility studies.

By 2020, though, JuPedSim had grown into a collection of hard-to-maintain executables driven by complex XML files: powerful, but opaque, and increasingly inaccessible to new users. At some point, we realized we couldn’t keep patching it forever. We needed to rethink what JuPedSim was for and who it was for. That’s where the real RSE story begins.

We set out to reshape JuPedSim into a modern, community-oriented library — and learned a lot along the way. The transformation was not primarily a technical one; it was a story about stakeholders, sustainability, and focus.

We rebuilt JuPedSim as a Python module with a modular C++ core, but the key change was philosophical: API over configuration, explicit over implicit. Instead of hiding behavior behind large configuration files, we encourage users to express their ideas directly in code. It’s not about convenience; it’s about clarity and reproducibility. We want people to understand what their simulation does, not just run it. This shift changed how the community interacts with the tool — making experiments more reproducible, the codebase more maintainable, and the system easier to teach and extend.

This journey wasn’t purely technical, it required full stakeholder commitment. None of this would have happened without support from domain scientists, developers, and institutions. Gaining that trust meant learning to communicate in their terms: why change was necessary, what would be gained, and how it would affect their work. Along the way, we faced the classic second-system problem — the temptation to rebuild everything, and learned that doing less was often better. We focused on the essentials, dropped non-core tools like JpsEditor, and spun out independent, focused projects such as PedPy, a Python analysis library that became a success in its own right.

Finally, we’ll reflect on governance, documentation, and long-term sustainability: how to keep a research software project alive beyond funding cycles and how to turn a codebase into a community.

This talk isn’t about code or performance. It’s about what it takes to turn a long-running research codebase into something sustainable — building trust, reducing complexity, and creating space for others to contribute.

In short: how we learned to do less, but better.

Keywords: pedestrian dynamics, research software engineering, sustainability, governance, community building, Python/C++ hybrid, documentation.

Download

IBIT Conference 2026

Mon, 26 Jan 2026 00:00:00 +0000

Overview and Topics

The 11th IBIT Conference on Crowdmanagement & Event Safety Management brought together professionals from diverse fields, including representatives from authorities, event organizers, operators, safety planners, security services, and researchers. This year’s conference featured six distinct sessions addressing various aspects of contemporary challenges in event safety:

Back to Basics: Introduction to Modern Safety Planning 2026
Hostile Vehicle Mititgation as a Guarantee for Safe Events?
Best Practice: Transferable Insights and Pragmatic Approaches
Where do the Rules come from? Regulations, Standards, and Guidelines
Operational Aspects and Procedures
Legal Aspects

Unlike previous years, the conference did not include a dedicated research session. However, research topics were integrated into several sessions, highlighting the importance of scientific input for practical applications in event safety—and vice versa.

Concert crowd model in the exhibition area

Contributions from the Pedestrian Dynamics Community

Speaker: Anna Sieben, Jette Degenhardt and Mohcine Chraibi

Transfer Workshop

The workshop “Pedestrian Dynamics: From Laboratory to Practice – From Practice to Laboratory,” presented by Anna Sieben and Jette Degenhardt, addressed the challenge of knowledge transfer between research and practice.

In research, knowledge is generated that does not always reach practical implementation. Scientific findings are relevant for all stakeholders, as they continuously influence guidelines and standards. Practitioners, on the other hand, possess extensive knowledge about crowds, often communicated differently than scientific insights and therefore not systematically incorporated into research. There is a lack of mechanisms to actively shape research based on practical input (and vice versa).

Anna and Jette provided an overview of current research questions and facilitated a discussion. The discussion revealed that practitioners desire more insights on the following topics:

Influence of lighting conditions and different floor surfaces on crowds
Research on perception of signs and wayfinding systems
Behavior under stress conditions and in extreme situations (perceptible dangers such as loud noises or fires)
Analysis of decision-making processes in dynamic crowd situations
Further development of VR tools for crowd simulation and training purposes

JuPedSim News: From Research to Practical Application

Mohcine Chraibi presented “From Research for Practice: A Simulation Tool for Pedestrian Flow Analysis,” introducing the JuPedSim Web App. The open-source software, previously used primarily within the scientific community, is now accessible through this online platform.

The tool enables the configuration and simulation of pedestrian flow scenarios based on validated models. While currently focused on basic scenarios, the long-term vision includes application to more complex situations. The JuPedSim team aims to enhance the tool’s practical applicability and seeks user feedback to further improve its functionality for broader use.

By: Jette Degenhardt

Pedestrian and Evacuation Dynamics 2025

Tue, 09 Sep 2025 00:00:00 +0000

The 12th edition of Pedestrian and Evacuation Dynamics (PED 2025) was held from 9 to 12 September 2025 at the Czech Technical University in Prague, Czechia, organized by the Faculty of Information Technology and the Faculty of Civil Engineering—the youngest and oldest faculties of the university, respectively.

Altogether, 138 researchers from 19 countries participated in the conference, discussing over 60 oral and 41 poster presentations. Forty-two contributions were transformed into proceedings papers, published on 22 June 2026 in Acta Polytechnica CTU Proceedings.

Starting with a Practical PhD Workshop

The first day, 9 September 2025, was dedicated to the PhD workshop focusing on data acquisition from pedestrian experiments. The organising team and two invited lecturers welcomed 35 participants.

After a brief introduction by Hana Najmanová to the ethical and legislative aspects of pedestrian experiments, the participants signed informed consents and became part of a bottleneck evacuation experiment on site. With the help of the Image Processing Laboratory, the experiment was recorded and the data were directly available for the afternoon practical part of the event.

Snapshot from the workshop experiment. Participants of the workshop were provided with the raw data together with pre-processed trajectories for comparison.

The crucial part of the workshop consisted of lectures on “Planning and analyzing controlled pedestrian experiments” and “Extraction of trajectories with help of PeTrack” by Alica Kandler and Maik Boltes from Forschungszentrum Jülich. These theoretical tutorials were directly linked to afternoon practical sessions in PC labs, during which the participants had the option to acquire trajectory data by means of PeTrack or simple manual techniques, to obtain calibration data from the camera, and also to discuss various bureaucratic aspects of pedestrian experiments.

The organizers received very positive feedback regarding the practical focus of the workshop, with the main complaint being: “The practical part should not be for three hours, but rather for three days!”

Main Conference

The main conference from 10 to 12 September 2025 took place in the Atrium of the Faculty of Civil Engineering and adjacent lecture and seminar rooms, providing an open and well-connected place for meetings and discussions, not only during the two poster sessions.

Group photo in the Atrium. Unfortunately, some participants are missing due to discussions after a wonderful keynote.

Every day started with a keynote session: Maik Boltes summarized the journey towards precise data acquisition from pedestrian experiments (related to the workshop); He Wang focused on AI and deep-learning methods rising in pedestrian dynamics research (two parallel sessions were dedicated to machine learning and AI); and Simo Hostikka reminded the audience of the fire safety engineering origin of the field by talking about the coupling of fire and evacuation models.

Parallel sessions were mostly divided (as is standard at PED) into a track on the analysis of pedestrian data (experimental or real-world) and a track on modelling and simulations. However, a high number of contributions combined both approaches, which is great for pursuing the common research goal. A significant portion of the contributions focused on crowd management and behavioural aspects of pedestrian dynamics, demonstrating the high interdisciplinarity of the PED conferences.

Thanks to the supporting technical staff, all presentations were streamed and recorded to allow community members who could not attend to stay partially connected.

Social Events

In addition to the opening and closing social gatherings, the conference programme was enriched by two main social events: a concert of the Academic Orchestra of the Czech Technical University in Prague in the Bethlehem Chapel, and the conference dinner in the brewery of the Strahov Monastery. The organizers deliberately let the participants get to the social event venues on their own via public transport (the best way to travel around Prague), which turned out to be a good occasion to explore the city.

Closing Words

The 13th PED conference will be hosted by Tsinghua University in Beijing, China, in 2027.

Now, three years after being given the opportunity to host PED, and almost one year after the conference, the journey for us as organizers ends with the publication of the proceedings in Acta Polytechnica CTU Proceedings — 42 papers, open access — and the writing of this report. We believe that thanks to the PED community, the valuable contributions, and the nice atmosphere among the participants, PED 2025 became a great and memorable event.

By: Pavel Hrabák, Hana Najmanová, Jana Vacková

Edited by: Mohcine Chraibi

PTED 2024

Fri, 13 Dec 2024 00:00:00 +0000

Workshop on Pedestrian Traffic and Evacuation Dynamics 2024

The 5th Workshop on Pedestrian Traffic and Evacuation Dynamics (2024) was held in Hefei, Anhui Province in China from December 13th to 15th. It was organized by Prof. Weiguo Song and Prof. Jun Zhang from the State Key Laboratory of Fire Science of University of Science and Technology of China and the National & Local Joint Engineering Research Center of Thermal Safety Technology to discuss new theories, new methods, new technologies and the trends in the research of pedestrian traffic and evacuation dynamics.

Group photo of the workshop.

Participants

110 scholars from 40 institutions including the University of Melbourne, Tsinghua University, Chongqing University, University of Science and Technology of China, Beijing Jiaotong University, Sun Yat-sen University, Wuhan University of Technology, Southwest Jiaotong University, Shanghai Institute of Technology, Shanghai Maritime University, Jiangsu University, Fuzhou University, Nanjing University of Technology, Beihang University, University of Tokyo, Hong Kong Polytechnic University, and City University of Hong Kong participated.

Participants and inaugural opening speech by Prof. Song.

Talks

Prof. Majid Sarvi from the University of Melbourne in Australia, Prof. Zhou Tiejun from Chongqing University, and Prof. Chen Tao from Tsinghua University gave Plenary speeches at the workshop. Prof. Han Xiao from Beijing Jiaotong University, Associate Prof. Cao Shuchao from Jiangsu University, Associate Prof. Xiao Yao from Sun Yat-sen University, Prof. Fu Ming from the Hefei Institute of Public Safety at Tsinghua University, and Associate Prof. Xie Qimiao from Shanghai Maritime University delivered keynote speeches. They shared the latest progress in experimental research, modeling, and application of new technologies in pedestrian traffic and emergency evacuation and provided prospects for the future development of this field.

Besides, 33 experts and scholars gave wonderful oral presentations. The attendees had in-depth exchanges and discussions on topics such as mixed traffic dynamics, pedestrian evacuation under disaster conditions, personnel safety protection technology, virtual reality technology, and personnel psychology. They shared their latest research results, experiences, and viewpoints, and actively explored new research directions and solutions. The open and interactive environment provides valuable learning and communication opportunities for participants, promoting the collision of ideas and academic progress in various research directions.

QA-session and moderation by Prof. Zhang.

During the workshop, some of the researchers visited the State Key Laboratory of Fire Science and watched the fire whirlwind demonstration, and learned about the development of the institute and its research achievements in fire science and emergency evacuation.

Looking ahead

The 6th Workshop on Pedestrian Traffic and Evacuation Dynamics (2025) is set to convene next year again at the University of Science and Technology of China, Hefei, China. The exact date has not yet been set but will be announced soon.

By: Jun Zhang

Edited by: Mohcine Chraibi

Traffic and Granular Flow 2024

Mon, 02 Dec 2024 00:00:00 +0000

The 15th edition of the Traffic and Granular Flow (TGF24) conference took place from December 2nd to 5th, 2024, in the beautiful city of Lyon, France. Held at the historic Palais Hirsch and the equally magnificent Palais de la Bourse, the conference offered an inspiring setting along the Rhône River.

Emergence of behavioral repertoires in entrance situations: Interplay between physical, collective and inter-personal factors - Tom Postmes

Since its inception in 1995, TGF has been a biennial event that brings together a diverse, interdisciplinary audience, including experts in physics, computer science, engineering, and fields such as pedestrian dynamics, granular flows, collective animal behavior, and urban mobility. This year’s conference continued the tradition of fostering collaboration and innovation by presenting the latest developments and ideas in these areas.

The scope of TGF covered a wide range of topics, including:

Internet traffic and cyber systems,
Pedestrian and evacuation dynamics,
Urban Mobility,
Trafic,
Collective animal behavior,
Swarm behavior, and
Collective dynamics of biological systems.

Researchers from across the globe—spanning the USA, UK, Europe, Australia, New Zealand, Japan, China, India, and more—gathered to share their work and form new collaborations, which were facilitated by the smooth organization offering plenty of time for discussions between the partecipants. A particular highlight was the keynote by Andreas Schadschneider, who delivered a retrospective on nearly 30 years of TGF, celebrating its role in connecting scientists across disciplines and continents. Andreas explained how the conference started from a small group of physicists in Germany interested in traffic from a physical perspective and evolved into the global community studying a variety of subjects as it is today.

TGF Conference Series: A retrospective by Andreas Schadschneider

Another memorable presentation was by Ernesto Altshuler, whose talk, From Vortices to People: An Experimental Approach to Traffic, captivated the audience. Ernesto showed, through beautiful videos and funny anecdotes, how systems of different sizes show properties related to traffic systems and concluded his presentation showing some puzzling experiments on collective motion of people.

The conference also featured impactful contributions from attendees. For instance, Paul Geoerg and his co-authors presented their research on Mapping Crowd Safety Risks at a Global Scale, which introduced a meticulously curated database of 1,780 crowd-related events from 2000 to 2025. This database tracks incidents such as crowd crushes, terror attacks, and cancellations due to extreme weather, providing valuable insights for safety planning and risk assessment.

During the conference there was also the opportunity to showcase various software that have been created to help researchers working in the community. For instance, Peter Thompson presented a tool intended to visualize results from simulation or experiment in 3D to help gaining a more multi-facet inside into collected data. Ahmed Alia presented a comparison of Deep Learning-based Computer Vision tools used to track people’s heads in videos, a kind of software that is increasingly used by a variety of researchers to collect data on crowds.

In her invited talk Marta Gonzalez discussed traffic from a macroscopic perspective focusing on mobility patterns of people through data collected on mobile phones. Marta provided an overview of her impressive career presenting a number of graphs in videos in an informative way. On a smaller scale Xiaolu Jia presented a method to obtain origin-destination data of people moving inside a train station in Tokyo (Japan).

The conference also included a large number of presentations on pedestrian dynamics, vehicular traffic, collective animal behavior, network theory, and also some touching aspects closer to psychology.

Impressions and Atmosphere

TGF is renowned not only for its academic depth but also for its welcoming and collaborative environment. The informal atmosphere encourages open dialogue, making it easy for participants to engage in meaningful discussions regardless of their gender, nationality, or professional position. First-time attendees often remark on the strong sense of community and the openness to new ideas. Although TGF is held biannually, some participants also meed during the PED conference sharing some of the topics but specifically focused on pedestrian traffic.

This inclusivity has been a hallmark of TGF, thanks to the efforts of long-standing contributors like Andreas Schadschneider and Armin Seyfried, who have nurtured this culture over the years. Once again, the conference attracted a wide range of scientists from around the world, with a remarkable balance of ages and genders. Such diversity and openness ensure that the value of the content takes precedence over status.

Although the famous festival of lights was not part of the official program of the conference, many participants took the chance to enjoy the light show in a crowded scene. Probably given the crowd researchers in the conference a more realistic feeling for high densities.

Festival of lights, Lyon 2024.

Looking Ahead

As the conference came to a close, the steering committee announced that the next edition, TGF26, will be hosted by Nikolai Bode at Bristol University. With its strong foundation of inclusivity and innovation, TGF is set to continue thriving as a platform for cutting-edge research and collaboration. TGF26 is planning to continue the multidisciplinary tradition of TGF and will further extent its coverage to include also robotics, which is increasigly playing an important role in traffic through automated vehicles and swarm robotics. Of course, content from more “classical” topics typical of the past editions of TGF will be also welcome.

By: Ahmed Alia, Thomas Chatagnon, Claudio Feliciani

Edited by: Mohcine Chraibi

IBIT Conference 2023

Wed, 08 Nov 2023 00:00:00 +0000

The 9th IBIT Conference on Crowdmanagement & Event Safety Management took place in Cologne in November. Related to the slogan “Knowledge, Network, Responsibility”, current challenges, basic knowledge, best practices and research activities were reported and discussed over two days.

Venue of the IBIT Conference: RheinEnergieSTADION in Cologne

Sabine Funk and her IBIT team put together a comprehensive program that included the following topics, among others:

Cooperation and communication during events
Inclusion of people with disabilities
Technical support for event safety (e.g. simulations, AI and drones)
Last mile / Zone Ex from different perspectives (legal, operational, …)
Best practice examples for e.g. specific safety concepts and managing threat situations

Most contributions and discussions were in German - although the English track is growing. On the 2nd day of the conference there was an international track presented by the YES Group, in which, among other things, the lessons learned from the Loveparade were discussed.

In addition to presentation and discussion panels, the exhibition area offered the opportunity to get in touch with various companies, institutes and organizations active in the field of event safety and to discuss crowd management strategies using this vivid model as an example:

New approach for modeling crowds at major events ;)

As a part of the side program, the topic “Knowledge transfer from research to practice and vice versa” was highlighted in an open discussion round. Barriers to transfer were collected, such as the lack of a common vocabulary, lack of detailed understanding of each other’s activities and different time scales of research/practice. While the list of challenges grew longer as the discussion progressed, there were only a few proposed solutions - such as the preparation of scientific publications for specific target groups from practice. The discussion showed that there is still a lot of work to be done in terms of transfer. However, the conference with its established research track is already an effective platform for advancing the exchange between research and practice.

The IBIT conference will take place again next year. The exact date and location have not yet been set but will be announced here.

By: Jette Degenhardt

IAFSS Symposium 2023

Mon, 23 Oct 2023 00:00:00 +0000

The International Association for Fire Safety Science organized a Symposium in Tsukuba, Japan during October 2023. This week-long symposium included a track on human behaviour in fire and evacuation and a dedicated workshop on the topic which took place the Sunday before the conference. The Sunday workshop included a set of presentations and discussions related to the ongoing work of the permanent group on human behaviour in fires led by Enrico Ronchi and Erica Kuligowski. This included the arrangement of a webinar series and the development of a research agenda for the field. Several important discussion points were brought up, including the need for training opportunities for early career researchers and a specific event entirely dedicated to the field.

During the symposium, three sessions were dedicated to evacuation along with two keynote talks. The conference was opened by a tribute to Rita Fahy by the chair of IAFSS, Brian Meacham, Erica Kuligowski and Daniel Nilsson. Rita was a pioneer of the field of evacuation modelling with her model Exit89 which was developed in the 80ies and has been used for the investigation of the World Trade Centre performed by the National Institute for Standards and Technology (NIST). Rita passed away during summer 2023, leaving a strong impact on the human behaviour in fire community with her work on the panic misconception, evacuation simulations and quantification of evacuation behaviour for fire engineering applications. The tribute was followed by the Howard Emmons award-winning talk by David Purser focused on understanding evacuation behaviour in presence of smoke and toxic gases.

Exploring single-line walking in immersive virtual reality by Luke de Schot.

The talk gave a comprehensive overview on the research in the domain providing valuable insights into the state-of-the-art of fire evacuation behaviour. The second keynote of the symposium was also on the topic of evacuation and was given by Yuki Akizuki. The focus of Akizuki’s talk was evacuation route design considering the impact of reduced visibility conditions due to smoke. Presentations during the symposium included three sessions on evacuation and human behaviour in fire and were chaired by Enrico Ronchi, Rosaria Ono and Daniel Nilsson. They included a wide range of topics, such as the use of novel methods to collect pedestrian movement and evacuation behaviour data (e.g., VR and online experiments), exploration of laypeople knowledge of evacuation, design of evacuation signage and way-finding systems, and evacuation of vulnerable groups such as children and people with functional limitations.

Sheldon Tieszen Student Awards (Foto by Peter Thompson).

On the Friday after the Symposium, a group of approximately 15 researchers interested in evacuation behaviour gathered in an informal meeting arranged at the University of Tokyo by Claudio Feliciani. This was a valuable opportunity to discuss collaborative efforts in the field and share experiences concerning pedestrian and evacuation data collection and modelling activities. All in all, the IAFSS Symposium and its related activities have been a fantastic opportunity for the evacuation community to meet and learn from each other research efforts in the field.

The 15th IAFSS International Symposium is set to convene in 2026 at the scenic locale of La Rochelle, France. This event promises to continue the esteemed tradition of fostering dialogue and sharing groundbreaking research in the field of fire safety science.

By: Enrico Ronchi

Edited by: Mohcine Chraibi