REVIEW 4 major objections 5 minor 38 references
Numerical investigation of the effect of macro control measures on epidemics transport via a coupled PDE crowd flow - epidemics spreading dynamics model
T0 review · 4 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read This paper builds a coupled PDE model of crowd flow, airborne contagion, and ventilation, and uses it to rank macro control measures by how many people become exposed during an evacuation.
desk verdict Worth a careful look for the ventilation-direction reversal, but Eq. (19) has a units and boundary-condition problem that undermines the quantitative exposure comparisons until fixed. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The carrying object is the infection-coefficient field $\beta(x,t)$, a transportable quantity that converts infected presence into new exposures: the source term for exposed density is $\beta_I\rho_S$. The field is advected by the ventilation velocity $U_G$, diffuses with effective turbulent viscosity $\sigma$, decays by aerosol settling $\nu\beta$, and is replenished proportionally to infected density $\rho_I/\rho$. Computing $U_G$ as the gradient of a Laplace potential $\Psi$ with Neumann duct boundary conditions is what lets ventilation rate and direction enter the model without resolving indoor turbulence. The crowd-flow component supplies both the advection velocity for the epidemic densities and the Eikonal-based desired direction toward exits, so the spatial density profile—not just a uniform mixing assumption—controls where exposure happens.
What would settle it
Replace the potential-flow air field with a measured or turbulence-resolving airflow in the same geometries and compare predicted exposed totals; if the direction of airflow that minimizes exposure changes, the ventilation conclusions are artifacts of the ideal airflow assumption.
Extended reading notes
Core claim
This paper's central claim is that epidemic transport in a closed, evacuating space can be described by a single coupled system in which the infection field is not imposed by hand but emerges from crowd motion, ventilation, and infected density. The force of infection is $\beta_I = i_0\beta(x,t)$, where $\beta$ solves a drift-reaction-diffusion equation $\beta_t + \nabla\cdot(\beta U_G) = \nabla\cdot(\sigma\nabla\beta) - \nu\beta + \rho_I/\rho$, and $U_G = \nabla\Psi$ is a steady potential-flow velocity field obtained from Laplace's equation with duct-like boundary conditions. Advecting the SEISV densities with the crowd velocity from a second-order macroscopic flow model closes the loop. In the numerical tests, the predicted exposed total falls when ventilation rate rises, when air moves against pedestrian flow, when maximum speed rises, when pressure coefficient $C_0$ rises (larger average distances), and when masked or vaccinated pedestrians are present; the direction effect reverses in large crowds near an exit, where a clean-air pocket of the size seen in smaller rooms no longer reaches most of the crowd.
Load-bearing premise
The ventilation air is treated as a steady, smooth, incompressible potential flow; real indoor airflow is turbulent and altered by the crowd itself, so the exposure reductions attributed to ventilation rate and direction rest on this idealization.
Editorial extensions
If this is right
- If the model is right, ventilation rate can be treated as a tunable exposure-reduction knob, but its direction should be chosen per geometry; airflow against walking direction helps small single-exit rooms, while a pattern that sweeps clean air across a large waiting area helps big crowds.
- Faster pedestrian movement reduces exposure mainly by shortening evacuation time; the benefit saturates or reverses when a narrow exit turns higher speed into stronger clogging.
- Enforcing larger average interpersonal distance (higher $C_0$) lowers exposure and evacuation time in single-exit congestion, yet contributes little when multiple or wide exits already keep densities low.
- A fraction of masked or vaccinated pedestrians reduces the exposed total additively with the other levers, bringing the lowest observed exposure in the tests to about 4% of the crowd.
- Real-time spatial density and infection-field profiles, rather than room-averaged totals, are the natural signals for deciding when to change ventilation or issue spacing and speed recommendations.
Reading between the lines
- Beyond the paper, the potential-flow airflow is the linchpin of the ventilation conclusions; replacing it with a turbulence-resolving airflow solver in the same geometries would directly test whether the direction rankings survive real indoor air motion.
- Because the model outputs full spatiotemporal density and infection fields, it invites a closed-loop control formulation—ventilation rate, direction, and speed or spacing recommendations updated from measured crowd state—which the paper mentions as a possibility but does not implement.
- The simulations fix recovery, latency, and vaccination rates at zero; extending the SEISV equations with nonzero recovery, latency, and vaccination rates would determine whether the same control rankings hold for longer exposures or partially immune populations.
- A tracer-gas experiment in a ventilated room with a controlled walking crowd could probe the predicted clean-air pocket at the exit and the reversal of direction preference in large crowds.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper presents a coupled PDE model for crowd flow and epidemic spreading in confined spaces. The crowd component is a second-order Payne-Whitham-style model with an Eikonal equation for the desired walking direction; the epidemic component is a macroscopic SEISV model; infection transmission is represented through a drift-diffusion equation for the infection coefficient, with air flow obtained from a potential-flow solution of Laplace's equation. The authors implement this model with finite volume schemes (Roe for the crowd model and Rusanov for the advection equations, with MUSCL reconstruction and a fast sweeping method for the Eikonal equation) and run numerical experiments for one-exit, two-exit, and corridor geometries, varying ventilation rate and direction, maximum pedestrian speed, pressure coefficient, exit width, and the fraction of masked/vaccinated individuals. The main reported findings are that total exposure decreases with higher ventilation rate, with air flow against pedestrian motion, with higher maximum speed, with higher pressure coefficient, and with masked/vaccinated individuals, with a noted reversal for large crowds near an exit.
Significance. If taken as a deliberately idealized modeling study, the paper is useful: it constructs a complete coupled PDE system, states all parameter values, and performs a broad sweep of control scenarios without fitting parameters to reproduce the exposure percentages. The numerical infrastructure is standard and the parameter choices are mostly traceable to the cited literature. The main value is qualitative hypothesis generation for macro-control measures in evacuating crowds. However, the quantitative exposure percentages and the rankings of ventilation scenarios are not yet fully grounded, because the infection-coefficient equation has an apparent dimensional inconsistency and missing boundary conditions, and because no grid-convergence or sensitivity study is reported. The ventilation conclusions also rest entirely on an inviscid potential-flow idealization of indoor air movement.
major comments (4)
- [Section 2.2, Eq. (19)] Equation (19) is not dimensionally consistent as written. Since beta_I = i0 beta in Eqs. (15)-(17) is an infection rate with units 1/time, beta itself has units 1/time when i0 is dimensionless. Then beta_t, the advection term, and the diffusion term each carry units 1/time^2, while the source term rho_I/rho is dimensionless. If beta is instead interpreted as dimensionless, the differential terms carry 1/time and the source term remains dimensionless. Either way, one term in Eq. (19) has mismatched units. This directly affects the quantitative exposure percentages in Section 4 and the interpretation of i0. Please introduce an explicit rate constant in the source term, state the units of i0, sigma, nu, and beta, and rerun or justify the simulations with a consistent formulation.
- [Section 3, boundary conditions paragraph] No boundary conditions are given for the infection-coefficient equation (19). The ghost-cell conditions listed in Section 3 apply to the crowd flow model only, and Section 2.2 supplies only beta(x,0)=0. Boundary conditions matter most at the inflow and outflow ducts, where the clean-air pockets and ventilation-direction effects in Figs. 7-12 and 22 are controlled by how beta is transported into and out of the domain. A Dirichlet condition beta=0 at supply ducts versus homogeneous Neumann conditions at walls and exhausts can change both the magnitude of the exposed fraction and the ranking of airflow directions. Please specify the boundary conditions used in the simulations and test the sensitivity of the reported trends to that choice.
- [Section 4] No grid-convergence or sensitivity study is reported. All one-exit and two-exit runs use Delta x = 0.05 m and Delta t = 2e-3 s, while the corridor runs use Delta x = 0.1 m and Delta t = 5e-3 s, with no refinement study. Many of the reported exposure differences are only a few percentage points, for example the reductions of up to about 6% in Section 4.1.2, so numerical dissipation and dispersion could be of the same order as the effects being claimed. Please add at least one mesh-refinement test for each geometry and report the resulting exposure percentages to confirm the conclusions are not grid-dependent.
- [Section 2.2, Eqs. (20)-(21)] The ventilation velocity field U_G is obtained as a potential flow of an inviscid, irrotational, incompressible fluid. This ignores turbulence, crowd-air interactions, and duct geometry, and the spatial structure of U_G is exactly what determines the 'pocket of clean air' near the exit that drives the ventilation-direction conclusions in Figs. 7-8 and Sections 4.1-4.3. The paper should state this limitation more prominently and, for at least one scenario, compare the potential-flow velocity field with a more realistic field, such as a RANS solution or measured data, to show that the qualitative trends are robust.
minor comments (5)
- [Section 2.2, Eq. (20)] Equation (20) contains a typo: 'where Psi(x is' should be 'where Psi(x) is'.
- [Section 2.2, after Eq. (18)] The sentence defining the initial condition states rho0 = rhoS0 + rhoE0 + rhoS0 + rhoV0; the third term should presumably be rhoI0.
- [Section 3, CFL condition] In the CFL condition (29), 'therm' should be 'term'.
- [Appendix] The limiter is named 'van Albaba' but should be 'van Albada'.
- [Figure 1 caption] The caption says 'Equilibrium speed-density relations for flow (left) and velocity (right)', but both panels appear to be speed-density relations; please clarify what distinguishes the two panels.
Circularity Check
No significant circularity: the reported simulation outcomes are direct outputs of the stated coupled PDE model, with no fitted-input predictions, self-citation chains, or imported uniqueness claims.
full rationale
The paper does not fit parameters to reproduce the reported exposure percentages; all results are computed from the stated PDEs with declared parameter values (e.g., i0 = 0.04, ν = 0.5, σ = 1.2e-3), and no data or target outcomes are used to tune the model. The ventilation, speed, spacing, and vaccination findings are consequences of the model equations and of the chosen simulation scenarios, not of parameters reverse-engineered from the results. The model equation for β, Eq. (19), is introduced following external references [27] and [26], which are not the present authors' own works, so no load-bearing self-citation is present. No uniqueness theorem is invoked, and no ansatz is smuggled in through a self-citation chain. The potential lack of boundary conditions and the apparent unit inconsistency in Eq. (19) are validity or correctness concerns that could undermine confidence in the quantitative predictions, but they do not make the predictions equivalent to the model inputs by construction; an ill-posed equation is a different failure mode from circularity. Overall, the derivation chain is self-contained as a numerical modeling study, and no circular step can be exhibited from the paper's own equations.
Assumptions & free parameters
free parameters (9)
- alpha in speed-density relation =
7.5
- C0 (anticipation factor) =
0.5, 0.8, 1.2
- tau (relaxation time) =
0.6 s
- i0 (infectivity scaling) =
0.04
- nu (aerosol settling rate) =
0.5
- sigma (effective turbulent viscosity) =
1.2e-3
- umax (maximum pedestrian speed) =
1.4 and 2 m/s
- uin (ventilation speed) =
5 and 10 m/s
- initial infected fraction =
0.25 and 0.1 of initial density
assumptions (6)
- domain assumption Pedestrian flow is described by a second-order macroscopic Payne-Whitham-type model with isotropic pressure and the specified speed-density relation.
- domain assumption Pedestrians choose their desired direction from the negative gradient of the Eikonal solution, i.e., they have global knowledge of the density field and seek the shortest travel time.
- domain assumption The ventilation-induced air flow is inviscid, incompressible, and irrotational, so it is fully described by a potential flow from Laplace's equation.
- domain assumption The infection coefficient beta evolves according to the drift-reaction-diffusion equation with source rho_I/rho and constant parameters, independent of the crowd dynamics except through the source.
- domain assumption During the time horizon, recovery, progression from exposed to infected, and vaccination are negligible (kappa=theta=xi=0).
- standard math The numerical schemes (Roe with entropy fix, Rusanov, MUSCL reconstruction, fast sweeping method) converge to the exact solution of the PDE system for the chosen grids and time steps.
Cite this review
Pith. "Pith review of Numerical investigation of the effect of macro control measures on epidemics transport via a coupled PDE crowd flow - epidemics spreading dynamics model." pith.science (2026). https://pith.science/paper/UD6MG4NM
@misc{pith2026241116223,
author = {Pith},
title = {Pith review of: Numerical investigation of the effect of macro control measures on epidemics transport via a coupled PDE crowd flow - epidemics spreading dynamics model},
year = {2026},
howpublished = {\url{https://pith.science/paper/UD6MG4NM}},
note = {Machine review of arXiv:2411.16223}
}
abstract
This work aims to provide an approach to the macroscopic modeling and simulation of pedestrian flow, coupled with contagion spreading, towards numerical investigation of the effect of certain, macro-control measures on epidemics transport dynamics. To model the dynamics of the pedestrians, a second-order macroscopic model, coupled with an Eikonal equation, is used. This model is coupled with a macroscopic Susceptible-Exposed-Infected-Susceptible-Vaccinated (SEISV) contagion model, where the force-of-infection $\beta$ coefficient is modeled via a drift-diffusion equation, which is affected by the air-flow dynamics due to the ventilation. The air-flow dynamics are obtained assuming a potential flow that can imitate the existence of ventilation in the computational domain. Numerical approximations are considered for the coupled model along with numerical tests and results. In particular, we investigate the effect of employment of different, epidemics transport control measures, which may be implemented through real-time manipulation of i) ventilation rate and direction, ii) maximum speed of pedestrians, and iii) average distances between pedestrians, and through iv) incorporation in the crowd of masked or vaccinated individuals. Such simulations of disease spreading in a moving crowd can potentially provide valuable information about the risks of infection in relevant situations and support the design of systematic intervention/control measures.
Figures
Figures from the paper (21 more)
Reference graph
Works this paper leans on
-
[1]
J.P. Agnelli, F. Colasuonno, and D. Knopoff. A kinetic theory appro ach to the dynamics of crowd evacuation from bounded domains. Mathematical Models and Methods in Applied Sciences , 25(1):109 – 129, 2015. 27 0 5 10 15 20 25 30 35 40 45 50 0 5 10 15 20 0 5 10 15 20 25 30 35 40 45 50 0 5 10 15 20 Figure 21: Corridor with a 4 m exit: Velocity field UG vector...
work page 2015
-
[2]
A spatial kinetic model of crowd evacuation dynamics with infectious d isease contagion
Juan Pablo Agnelli, Bruno Buffa, Dami´ an Knopoff, and Germ´ an Torres. A spatial kinetic model of crowd evacuation dynamics with infectious d isease contagion. Bulletin of Mathematical Biology , 85(4), 2023
work page 2023
-
[3]
R. M. Anderson, H. Heesterbeek, D. Klinkenberg, and T. D. Hollin gsworth. How will country-based mitigation measures influence the course of the covid-19 epidemic? Lancet, 395:931–934, 2020
work page 2020
-
[4]
J. Arino and P. Van den Driessche. A multi-city epidemic model. Mathe- matical Population Studies , 10:175–193, 2003
work page 2003
-
[5]
Ventilation for acceptable indoor air quality, 2015
ASHRAE. Ventilation for acceptable indoor air quality, 2015
work page 2015
-
[6]
G. Bertaglia and L. Pareschi. Hyperbolic models for the spread of epidemics on networks: Kinetic description and numerical methods. ESAIM Math. Model. & Numerical Analysis , 55:381–407, 2021. 28 0 20 40 60 80 100 120 0 5 10 15 20 25 30 35 No ventilation Airflow pattern I Airflow pattern II 0 20 40 60 80 100 0 5 10 15 20 25 30 35 No ventilation Airflow pat...
work page 2021
-
[7]
S. Buchm¨ uller and U. Weidmann. Parameters of pedestrians, pe destrian traffic and walking facilities. IVT Schriftenreihe. Institut f¨ ur Verkehrspla- nung und Transportsysteme (IVT), ETH Z¨ urich, 132, 2006
work page 2006
-
[8]
M. Delmastro and G. Zamariola. Depressive symptoms in response to covid- 19 and lockdown: a cross-sectional study on the italian population. Nature Scientific Reports , 10, 2020. article no 22457
work page 2020
Show all 38 references
-
[9]
Riemann Solvers and Numerical Methods for Fluid Dynamics
E.F.Toro. Riemann Solvers and Numerical Methods for Fluid Dynamics . Springer Berlin, Heidelberg, 2009
2009
-
[10]
L. Guan, C. Prieur, L. Zhang, C. Prieur, D. Georges, and P. Be llemain. Transport effect of covid-19 pandemic in france. Annual Reviews in Control, 50:394–408, 2020
2020
-
[11]
M. Haghani. Crowd dynamics research in the era of covid-19 pan demic: Challenges and opportunities. Safety Science, 153, 2022. paper no. 105818. 29 Figure 23: Corridor with a 4 m exit: Total density profiles (top), infection coefficient profile (middle), and exposed pedestrians’ p...
2022
-
[12]
F. S. Hanseler and S. Hoogendoorn. Optimal crowd managemen t for con- gested metro stations. In In hEART 2018: 7th Symposium of the European Association for Research in Transportation, 5-7 September, Athens, Greece, 2018
2018
-
[13]
Harten and J
A. Harten and J. M. Hyman. Self adjusting grid methods for one - dimensional hyperbolic conservation laws,. J. Comput. Phys. , 50:253–269, 1983
1983
-
[14]
Hosseinloo, S
A.H. Hosseinloo, S. Nabi, A. Hosoi, and M.A. Dahleh. Data-driven c ontrol of covid-19 in buildings: a reinforcement-learning approach. IEEE Trans. on Automation Science and Engineering , 2023
2023
-
[15]
R. L. Hughes. A continuum theory for the flow of pedestrians. Transp. Res. Part B: Methodological , 36(6):507, 2002
2002
-
[16]
Kachroo, S.A
P. Kachroo, S.A. Wadoo, S. J. Al-nasur, and A. Shende. Pedestrian Dy- namics: Feedback Control of Crowd Evacuation . Springer, Berlin, 2008
2008
-
[17]
M. Kaushik. Potential flow theory. In Theoretical and Experimental Aero- dynamics. Springer, Singapore, 2019. 30 Figure 24: Corridor with a 4 m exit: Total density profiles (top), infection coefficient profile (middle), and exposed pedestrians’ profile at diffe rent time instances wi...
2019
-
[18]
Kim and A
D. Kim and A. Quaini. Coupling kinetic theory approaches for pede strian dynamics and disease contagion in a confined environment. Math. Models and Meth. in Applied Sciences , 30:1893–1915, 2020
1915
-
[19]
R. J. Leveque. Finite Volume Methods for Hyperbolic Problems . Cambridge University Press, U.K., 2002
2002
-
[20]
R. J. Leveque. Finite Difference Methods for Ordinary and Partial Dif- ferential Equations: Steady-State and Time-Dependent Pro blems. SIAM, Philadelphia, 2007
2007
-
[21]
M. J. Lighthill and G.B. Whitham. On kinematic waves, i:flow movement in long rivers. ii:a theory of traffic on long crowded roods. Proc. Royal Soc. A, A229:281–316, 1955
1955
-
[22]
Sundar, and J¨ org Kuhnert
Somnath Maity, S. Sundar, and J¨ org Kuhnert. A high-resolut ion meshfree particle method for numerical investigation of second-order macr oscopic pedestrian flow models. Applied Mathematical Modelling , 131:205 – 232, 2024. 31
2024
-
[23]
Canudas-de-Wit, A
M.U.B Niazi, C. Canudas-de-Wit, A. Kibangou, and P.-A. Bliman. Opt imal control of urban human mobility for epidemic mitigation. In 2021 60th IEEE Conference on Decision and Control (CDC) , pages 6958–6963, TX, USA, 2021
2021
-
[24]
H. J. Payne. Models of freeway traffic and control. Math. Models Publ. Sys. Simul. Council Proc. , 28:51–61, 1971
1971
-
[25]
Qian and X
H. Qian and X. Zheng. Ventilation control for airborne transmis sion of human exhaled bio-aerosols in buildings. Journal of Thoracic Disease , 10:S2295–S2304, 2018
2018
-
[26]
P. S. A. Salam, W. Bock, A. Klar, and S. Tiwari. Disease contagion models coupled to crowd motion and mesh-free simulation. Math. Models and Meth. in Applied Sciences , 31:1277–1295, 2021
2021
-
[27]
P. S. A. Salam, W. Bock, A. Klar, and S. Tiwari. Coupling pedestria n flow and disease contagion models. In N. Bellomo and L. Gibelli, editors, Crowd Dynamics, Volume 4: Analytics and Human Factors in Crowd Mod eling, pages 223–246. Springer International Publishing, Cham, 2023
2023
-
[28]
Sattenspiel and K
L. Sattenspiel and K. Dietz. A structured epidemic model incor porating geographic mobility among regions. Mathematical Biosciences, 128:71–91, 1995
1995
-
[29]
Tizzoni, P
M. Tizzoni, P. Bajardi, A. Decuyper, G. Kon Kam King, C. M. Schn eider, V. Blondel, Z. Smoreda, M. C. Gonzalez, and V. Colizza. On the use of human mobility proxies for modeling epidemics. Plos Computational Biology, 10, 2014
2014
-
[30]
Treiber and A
M. Treiber and A. Kesting. Traffic flow dynamics: Data, Models and Sim- ulation. Springer Berlin, Heidelberg, 2013
2013
-
[31]
Twarogowska, P
M. Twarogowska, P. Goatin, and R. Duvigneau. Comparative st udy of macroscopic pedestrian models. Transportation Research Procedia, 2:477– 485, 2014
2014
-
[32]
Twarogowska, P
M. Twarogowska, P. Goatin, and R. Duvigneau. Macroscopic mo deling and simulations of room evacuation. Applied Mathematical Modelling , 38(24):5781–5795, 2014
2014
-
[33]
Duives, and Serge P
Arco van Beek, Yan Feng, Dorine C. Duives, and Serge P. Hooge ndoorn. Studying the impact of lighting on the pedestrian route choice using v irtual reality. Safety Science , 174, 2024. Art. no. 106467
2024
-
[34]
Vynnycky and R
E. Vynnycky and R. White. An Introduction to Infectious Disease Mod- elling. Oxford Univ. Press, 2010
2010
-
[35]
S. A. Wadoo and P. Kachroo. Feedback control of crowd evac uation in one dimension. IEEE Transactions on Intelligent Transportation Systems , 11:182–193, 2010. 32
2010
-
[36]
Willem, F
L. Willem, F. Verelst, J. Bilcke, N. Hens, and P. Beutels. Lessons from a decade of individual-based models for infectious disease transmiss ion: a systematic review (2006-2015). BMC Infectious Diseases , 17, 2017. paper no. 612
2006
-
[37]
B. Zhao, Z. Zhang, X. Li, and D. Huang. Comparison of diffusion c har- acteristics of aerosol particles in different ventilated rooms by num erical method. In ASHRAE Transactions , volume 110 PART 1, page 88 – 95. 2004
2004
-
[38]
H. Zhao. A fast sweeping method for eikonal equations. Mathematics of Computation, 74(250):603–627, 2005. Appendix The idea behind the MUSCL scheme is to replace the piecewise constan t ap- proximation of the conserved variables of the first-order finite vo lume scheme by recons...
2005
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.