pith. sign in

arxiv: 2211.16652 · v1 · submitted 2022-11-30 · 🧮 math.DS · math.AP

Canards in a bottleneck

Annalisa Iuorio , Gaspard Jankowiak , Peter Szmolyan , Marie-Therese Wolfram This is my paper

Pith reviewed 2026-05-24 10:26 UTC · model grok-4.3

classification 🧮 math.DS math.AP
keywords canardsbottleneckFokker-Planck equationstationary profilesbifurcation diagramsingular perturbationdensity transitionnonlinear boundary conditions
0
0 comments X

The pith

Canard solutions generate density transitions exactly at the minimum width of a corridor bottleneck in the small-diffusion limit.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper constructs the stationary profiles of a nonlinear Fokker-Planck equation on corridors whose width reaches a global minimum inside the domain. In the small-diffusion limit three families appear: high-density states with a boundary layer only at the entrance, low-density states with a layer only at the exit, and mixed states whose high-to-low jump occurs inside the bottleneck. The mixed states are realized by canard solutions that are born at the narrowest cross-section; a bifurcation diagram in the plane of inflow and outflow rates determines which family is selected. The constructions are verified by direct numerical solution on corridors of varying shape.

Core claim

In the small diffusion limit the profiles are obtained constructively by using methods from geometric singular perturbation theory. Three main types appear: high density with an entrance boundary layer, low density with an exit boundary layer, and transitions from high to low density inside the bottleneck. The last family uses canard solutions generated at the narrowest point; their occurrence is organized by a detailed bifurcation diagram in the in- and outflow rates. The analytic picture is confirmed by computational experiments on corridors with variable width.

What carries the argument

Canard solutions generated at the narrowest point of the bottleneck, which carry the high-to-low density transition in the singular limit.

If this is right

  • High-density profiles with an entrance layer exist when inflow exceeds a critical threshold set by the bottleneck geometry.
  • Low-density profiles with an exit layer exist when outflow is sufficiently strong relative to inflow.
  • Mixed profiles with an interior canard transition exist in an open region of the inflow-outflow parameter plane.
  • The bifurcation diagram partitions the rate plane into regions corresponding to each of the three profile families.
  • The same constructions hold for any corridor whose width function satisfies the nondegenerate-minimum condition.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same canard mechanism is likely to organize time-dependent or stochastic versions of the same equation when diffusion remains small.
  • The organizing role of the width minimum suggests analogous interior transitions may appear in other transport models on domains with a single narrowest cross-section.
  • Numerical continuation in the diffusion parameter could be used to track how the canard persists for finite diffusion values.

Load-bearing premise

The corridor width has a single global nondegenerate minimum inside the domain.

What would settle it

Direct numerical solution of the stationary equation for successively smaller diffusion coefficients showing the density jump displaced from the width minimum would falsify the location of the canard.

Figures

Figures reproduced from arXiv: 2211.16652 by Annalisa Iuorio, Gaspard Jankowiak, Marie-Therese Wolfram, Peter Szmolyan.

Figure 1
Figure 1. Figure 1: Fast dynamics in (j, ρ)-space for a fixed value of ξ. The blue curve represents C0 consisting of the two branches C a 0 (attracting), C r 0 (repelling), and the fold line F. The green lines indicate orbits of the layer problem (11), while the blue dot represents the line of fold points F. The phase space for the reduced problem is [0, 1/4] × [0, 1], where j = 1/4 corresponds to the fold line. It follows fr… view at source ↗
Figure 2
Figure 2. Figure 2: Illustration of the reduced flow associated to Equations (8)-(9) described in Lemma 1 for [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Schematic representation of L (orange line) and L+ (orange curve) for (a) 0 < α < 1 2 and (b) 1 2 < α < 1. The orange dot corresponds to l, the blue curve represents C0, and the green lines correspond to the orbits of the layer problem. 10 [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Schematic representation of R (purple line) and R− (purple curve) for (a) 0 < β < 1 2 and (b) 1 2 < β < 1. The purple dot corresponds to r, the blue curve represents C0, and the green lines correspond to the orbits of the layer problem. Based on this geometric interpretation of the boundary conditions, we proceed with the con￾struction of the singular orbits by connecting L + ∪l and R− ∪r by means of the r… view at source ↗
Figure 5
Figure 5. Figure 5: Schematic illustration of the special orbits [PITH_FULL_IMAGE:figures/full_fig_p012_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Representation of the (α, β) bifurcation diagram for ε = 0. Here k(ξ) = 1 + a cos  2πξ b  with a = 0.3, b = 1.5. Red regions correspond to high density, blue regions to low density, and green regions to transitions from high to low density regimes. In the insets, the density ρ is shown as a function of ξ. The blue parts correspond to solutions of the reduced problem (16), whereas the green parts indicate… view at source ↗
Figure 7
Figure 7. Figure 7: Schematic representation of a singular solution of (8)-(9) with [PITH_FULL_IMAGE:figures/full_fig_p017_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Schematic representation of singular solutions of type 1-4 (rows 1-4, respectively). [PITH_FULL_IMAGE:figures/full_fig_p018_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Schematic representation of singular solutions of type 5-8 (rows 1-4, respectively). [PITH_FULL_IMAGE:figures/full_fig_p019_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Schematic representation in (ξ, ρ)-space of the slow portions (blue curves) of the possible singular orbits for (a) (α, β) ∈ γ17, (b) (α, β) ∈ γ13 ∪ γ24, and (c) (α, β) ∈ γ37 ∪ γ58. The orange and purple curves correspond to the projection of L+ and R−, respectively, on the (ξ, ρ)-space. Fast jumps from the slow solution in C r 0 to the slow solution in C a 0 are possible (a) at each ξ ∈ [0, 1], (b) for ξ… view at source ↗
Figure 11
Figure 11. Figure 11: Schematic representation of a solution (continuous red curve) to the full problem (5a)-(5b) obtained [PITH_FULL_IMAGE:figures/full_fig_p021_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: (α, β) bifurcation diagram for ε = 10−3 . Here k(ξ) = 1 + a cos  2πξ b  with a = 0.3, b = 1.5. The insets show the density ρ as a function of ξ. low density stationary states with ρ < 1 2 and a boundary layer on the right boundary. • Small β (which corresponds to low outflow as in the red regions G7 and G8), which leads to high density stationary states with ρ > 1 2 and a boundary layer on the left boun… view at source ↗
Figure 13
Figure 13. Figure 13: Representation of the 2D domain associated with a supergaussian [PITH_FULL_IMAGE:figures/full_fig_p023_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Phase diagrams in the (α, β) parameter space from the GSPT analysis for the supergaussian k along with some typical non-singular solutions for wm = 0.9 (top) and wm = 0.5 (bottom). faster as α and β increase. The maximum of Jk will also decrease linearly with minξ k(ξ). Nu￾merical experiments with two narrow sections of varying width suggest that this applies also for functions k with several (nondegenera… view at source ↗
Figure 15
Figure 15. Figure 15: Illustration of the flow J1 (solid colors on the left, wireframe on the right) and Jk (solid colors on the right) as a function of α and β. Recall that ρ 0,1 c → 1 2 as k(ξ ∗) → 1, so that from the green rectangles, only the darker one (top right) remains in the limit. Conclusion In this work, we investigate the steady-states of a 1D area averaged model describing pedes￾trian dynamics for unidirectional f… view at source ↗
read the original abstract

In this paper we investigate the stationary profiles of a nonlinear Fokker-Planck equation with small diffusion and nonlinear in- and outflow boundary conditions. We consider corridors with a bottleneck whose width has a global nondegenerate minimum in the interior. In the small diffusion limit the profiles are obtained constructively by using methods from geometric singular perturbation theory (GSPT). We identify three main types of profiles corresponding to: (i) high density in the domain and a boundary layer at the entrance, (ii) low density in the domain and a boundary layer at the exit, and (iii) transitions from high density to low density inside the bottleneck with boundary layers at the entrance and exit. Interestingly, solutions of the last type involve canard solutions generated at the narrowest point of the bottleneck. We obtain a detailed bifurcation diagram of these solutions in terms of the in- and outflow rates. The analytic results based on GSPT are further corroborated by computational experiments investigating corridors with bottlenecks of variable width.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 3 minor

Summary. The paper applies geometric singular perturbation theory (GSPT) to the stationary nonlinear Fokker-Planck equation with small diffusion and nonlinear inflow/outflow boundary conditions on a corridor whose width attains a global nondegenerate interior minimum. It constructs three families of limiting profiles: (i) high-density solutions with an entrance boundary layer, (ii) low-density solutions with an exit boundary layer, and (iii) transitional profiles that cross from high to low density via canard solutions generated at the bottleneck minimum. A bifurcation diagram in the in- and outflow rate parameters is obtained for the transitional case, and the analytic results are stated to be corroborated by numerical experiments on corridors of variable width.

Significance. If the GSPT construction and matching are complete, the manuscript supplies an explicit, geometrically organized bifurcation diagram for canard-mediated transitions in a concrete PDE model; this is a concrete advance in the application of slow-fast methods to diffusion problems with nonlinear boundary conditions. The numerical corroboration and the explicit geometric hypotheses on the corridor width are positive features.

minor comments (3)
  1. The statement of the critical manifold and its loss of normal hyperbolicity (presumably in §2 or §3) should include the explicit reduced flow on the manifold and the precise transversality condition at the bottleneck minimum to make the canard existence transparent.
  2. The bifurcation diagram in the rate parameters would benefit from an explicit statement of the curves separating the three regimes (e.g., the critical values of the inflow/outflow rates at which the canard appears).
  3. A brief remark on the error estimates between the constructed GSPT profiles and the true stationary solutions of the original PDE would strengthen the claim that the three types are exhaustive in the small-diffusion limit.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive summary of our manuscript, for highlighting its significance in applying GSPT to nonlinear Fokker-Planck equations with nonlinear boundary conditions, and for recommending minor revision. We are pleased that the geometric construction, bifurcation diagram, and numerical corroboration were viewed favorably.

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper applies standard geometric singular perturbation theory (GSPT) to the given stationary nonlinear Fokker-Planck equation with the stated nonlinear boundary conditions and the explicit geometric assumption that the corridor width has a global nondegenerate interior minimum. The three profile types, including the canard transitions generated at the bottleneck, and the bifurcation diagram in the rate parameters are obtained directly from the slow-fast structure and loss of normal hyperbolicity on the critical manifold; none of these steps reduce by construction to fitted parameters, self-definitions, or load-bearing self-citations. The analytic construction is corroborated by independent numerical experiments, confirming that the central claims remain self-contained against external GSPT results.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the applicability of geometric singular perturbation theory to construct stationary profiles in the small-diffusion limit for the given nonlinear boundary-value problem; the nondegenerate minimum of the width function is an additional structural assumption required for canard location.

axioms (1)
  • domain assumption Geometric singular perturbation theory applies to the stationary nonlinear Fokker-Planck equation with the given nonlinear boundary conditions in the small-diffusion limit.
    The abstract states that profiles are obtained constructively by using methods from GSPT.

pith-pipeline@v0.9.0 · 5706 in / 1372 out tokens · 25115 ms · 2026-05-24T10:26:04.983176+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

17 extracted references · 17 canonical work pages

  1. [1]

    Burger, S

    M. Burger, S. Hittmeir, H. Ranetbauer, and M.-T. Wolfram. Lane formation by side-stepping. SIAM Journal on Mathematical Analysis , 48(2):981–1005, 2016

    work page 2016

  2. [2]

    Burger, P

    M. Burger, P. A. Markowich, and J.-F. Pietschmann. Continuous limit of a crowd motion and herding model: analysis and numerical simulations. Kinetic & Related Models, 4(4):1025, 2011. 25

    work page 2011

  3. [3]

    Burger and J.-F

    M. Burger and J.-F. Pietschmann. Flow characteristics in a crowded transport model. Non- linearity, 29(11):3528–3550, 2016

    work page 2016

  4. [4]

    Cristiani, B

    E. Cristiani, B. Piccoli, and A. Tosin. Multiscale modeling of pedestrian dynamics, volume 12. Springer, 2014

    work page 2014

  5. [5]

    Derrida, E

    B. Derrida, E. Domany, and D. Mukamel. An exact solution of a one-dimensional asymmetric exclusion model with open boundaries. J. Stat. Phys. , 69(3-4):667–687, 1992

    work page 1992

  6. [6]

    Fenichel

    N. Fenichel. Geometric singular perturbation theory for ordinary differential equations. Jour- nal of Differential Equations , 31(1):53–98, 1979

    work page 1979

  7. [7]

    Iuorio, G

    A. Iuorio, G. Jankowiak, P. Szmolyan, and M.-T. Wolfram. A PDE model for unidirectional flows: stationary profiles and asymptotic behaviour. Journal of Mathematical Analysis and Applications, page 126018, 2022

    work page 2022

  8. [8]

    Iuorio, N

    A. Iuorio, N. Popovi´ c, and P. Szmolyan. Singular perturbation analysis of a regularized MEMS model. SIAM Journal on Applied Dynamical Systems , 18(2):661–708, 2019

    work page 2019

  9. [9]

    Iuorio and F

    A. Iuorio and F. Veerman. The influence of autotoxicity on the dynamics of vegetation spots. Physica D: Nonlinear Phenomena , 427:133015, 2021

    work page 2021

  10. [10]

    C. K. R. T. Jones. Geometric singular perturbation theory. In Dynamical Systems , pages 44–118. Springer Berlin Heidelberg, 1995

    work page 1995

  11. [11]

    C. Kuehn. Multiple Time Scale Dynamics . Springer International Publishing, 2015

    work page 2015

  12. [12]

    M. J. Lighthill and G. B. Whitham. On kinematic waves II. A theory of traffic flow on long crowded roads. Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences, 229(1178):317–345, 1955

    work page 1955

  13. [13]

    Maury and S

    B. Maury and S. Faure. Crowds in Equations: An Introduction to the Microscopic Modeling of Crowds. 09 2018

    work page 2018

  14. [14]

    P. I. Richards. Shock waves on the highway. Operations research, 4(1):42–51, 1956

    work page 1956

  15. [15]

    Szmolyan and M

    P. Szmolyan and M. Wechselberger. Canards in R3. Journal of Differential Equations , 177(2):419–453, 2001

    work page 2001

  16. [16]

    Szmolyan and M

    P. Szmolyan and M. Wechselberger. Relaxation oscillations in R3. Journal of Differential Equations, 200(1):69–104, 2004

    work page 2004

  17. [17]

    A. J. Wood. A totally asymmetric exclusion process with stochastically mediated entrance and exit. Journal of Physics. A. Mathematical and Theoretical , 42(44):445002, 10, 2009. 26

    work page 2009