Source-linked AI summary
Self-sustained nonlinear waves in traffic flow
Morris R. Flynn, Aslan R. Kasimov, Jean-Christophe Nave, Rodolfo R. Rosales, Benjamin Seibold
TL;DR
Because viscous regularization of traffic shocks has an ambiguous physical rationale, the paper examines whether inviscid continuum traffic equations support self-sustained waves. Its Lax-condition analysis identifies allowable shock speeds and finds analytic jamiton structures that agree with simulations; jamitons can include sharp density rises relevant to roadway design.
Problem
Because the physical rationale for viscous smoothing of traffic shocks is ambiguous, the paper asks whether inviscid traffic equations support self-sustained waves that help explain jams without bottlenecks.
Source: §1 · p. 3Method
The paper analyzes self-sustained traveling-wave solutions of inviscid continuum traffic equations, treating shocks as discontinuities and applying Lax entropy conditions.
Source: §1 · p. 3, §6 · p. 22Results
Lax conditions permit only self-sustained shocks overtaken by vehicles; analytic jamitons for the stated speed and pressure forms agree excellently with simulations and feature sharp, compact density rises.
Source: §6 · p. 22Takeaways & Limitations
Jamitons can produce downstream densities above ρM beyond a critical average density, informing choices of speed limits and roadway carrying capacity.
Source: §6 · p. 22Their spatiotemporal stability remains a major objective for future research.
Source: §6 · p. 22
Abstract
from arXiv · showhide
In analogy to gas-dynamical detonation waves, which consist of a shock with an attached exothermic reaction zone, we consider herein nonlinear traveling wave solutions, termed "jamitons," to the hyperbolic ("inviscid") continuum traffic equations. Generic existence criteria are examined in the context of the Lax entropy conditions. Our analysis naturally precludes traveling wave solutions for which the shocks travel downstream more rapidly than individual vehicles. Consistent with recent experimental observations from a periodic roadway (Sugiyama et al., New Journal of Physics, 10, 2008), our numerical calculations show that, under appropriate road conditions, jamitons are attracting solutions, with the time evolution of the system converging towards a jamiton-dominated configuration. Jamitons are characterized by a sharp increase in density over a relatively compact section of the roadway. Applications of our analysis to traffic modeling and control are examined by way of a detailed example.
1 Introduction and problem formulation
Because the physical rationale for viscous smoothing of traffic shocks is ambiguous, the paper investigates self-sustained traveling waves in inviscid continuum traffic equations. It connects these waves to detonation structure and proposes that they can represent the saturation of phantom jams.
1 Introduction and problem formulation
- The study uses a one-dimensional Payne–Whitham continuum model on a periodic road, with density, speed, pressure, and a desired-speed function.The pressure is assumed convex in specific volume, while desired speed decreases with density.Source: §1 · p. 2
- The analysis adopts ũ = ũ0(1 − ρ/ρM)^n for n close to 1, while deferring exact jamiton existence conditions on ũ to later work.The desired speed ũ0 is the desired travel speed on an otherwise empty road, and ρM is maximum density.Source: §1 · pp. 2–3
- The paper seeks self-sustained traveling waves in inviscid traffic equations, whose shock structure resembles self-sustained detonation waves.The motivation is that the physical rationale for the viscous term often used to smooth shocks is ambiguous.Source: §1 · p. 3
- Numerical experiments suggest that phantom jams arising from linear instabilities can ultimately saturate as jamitons, which may be analytically described under selected conditions.The paper relates jamiton formation from small perturbations to ignition and detonation in a metastable explosive medium.Source: §1 · pp. 3–4
2 Traveling wave solutions – jamitons
The section derives conditions for periodic jamitons by joining monotone traveling-wave profiles across admissible shocks. It rules out right-shock waves and gives a sonic-point condition and construction procedure for the physically consistent waves.
2 Traveling wave solutions – jamitons
- Except for pathological or unphysical choices, the first-order traveling-wave ODE has no smooth periodic solution, so periodic waves must join monotone profiles with shocks.The simplest case has one shock per period, as in the Sugiyama et al. experiment.Source: §2 · p. 5
- Admissible shocks must conserve mass and momentum and satisfy Lax conditions, which require one characteristic family to converge on the shock while the other passes through.These conditions enforce dynamical stability and causal shock evolution.Source: §2 · pp. 5–6
- For a left shock, m > 0 and the shock is compressive: density rises while vehicles decelerate as they overtake it.The continuous profile between shocks has decreasing density and increasing velocity.Source: §2 · p. 7
- Negative-mass-flux traveling waves are mathematically inconsistent, excluding right-shock jamitons from the admissible solutions.For m < 0, the smooth-profile requirement conflicts with the assumed decreasing ˜u(ρ) and the implied behavior of ˜u − u.Source: §2 · pp. 6, 8
- For traveling waves with m > 0, the desired solution has du/dη > 0 and dρ/dη < 0; smoothness requires a same-order numerator zero to cancel the sonic-point denominator zero, yielding a regularized ODE.Source: §2 · p. 8
- The regularization condition is analogous to the Chapman–Jouguet condition in detonation theory.Source: §2 · p. 9
- The construction fixes wave speed from the downstream state, applies Rankine–Hugoniot conditions for the upstream state, then integrates the profile to determine its period.The resulting period is the circumference of the periodic roadway; vehicle number is evaluated with the roadway’s fixed total vehicle count.Source: §2 · p. 10
- Jamitons smoothly accelerate traffic and reduce density downstream, then produce an abrupt speed drop and density increase at the shock, consistent with Sugiyama et al..The algorithm parameterizes periodic jamitons by downstream state; for examined cases, upstream density maps one-to-one to vehicle number.Source: §2 · p. 10
3 Alternative description of the momentum equation
The section tests whether self-sustained shock waves depend on the Payne–Whitham momentum equation. It outlines related Aw–Rascle and Helbing formulations and shows how their traveling-wave analysis connects to the earlier one.
3 Alternative description of the momentum equation
- The existence of self-sustained shock solutions is presented as a generic feature of hyperbolic conservation laws with forcing, not a detail unique to one continuum traffic model.Source: §3 · p. 11
- Under Aw–Rascle, the traveling-wave equation matches the earlier form, but its sonic point occurs when u − s = m^1/2c.The sonic singularity is regularized by aligning it with a root of ˜u − u; the analysis otherwise proceeds similarly with modified shock conditions.Source: §3 · p. 11
- Helbing’s generalized formulation recovers Payne–Whitham and Aw–Rascle through different choices of its two traffic pressures.Payne–Whitham uses p1 = p, p2 = 0; Aw–Rascle uses p1 = 0, p2 = −uρ^2 dp/dρ.Source: §3 · pp. 11–12
4 Phase plane analysis
The phase-plane analysis uses the sonic point's linearization to constrain jamiton profiles and links a zero eigenvalue to instability onset; numerical experiments show unstable uniform traffic saturating into jamitons.
4 Phase plane analysis
- Introducing phase-plane variable ξ yields two ODEs with the sonic point as a critical point, whose Jacobian provides eigenvalues r1 and r2.Source: §4 · pp. 12–13
- At the sonic point, the sign of r2 determines the sign of du/dη, while the Lax entropy conditions require u to increase monotonically away from a shock.Source: §4 · p. 13
- Because r1 = 0 makes the critical point linearly degenerate, fully analyzing nearby solution behavior requires examining leading-order nonlinearities.Source: §4 · p. 13
- The condition r2 = 0 coincides with the boundary where a constant uniform base state becomes unstable to infinitesimal disturbances.Source: §4 · p. 13
- When a uniform traffic state is linearly unstable, the instability consistently saturates into a jamiton-dominated state.Source: §4 · p. 13
- In the Helbing analysis, modifying F2,u does not alter the Jacobian eigenvalues r1 and r2.Source: §4 · pp. 13–14
5 An example
The example develops exact jamiton solutions for selected traffic pressure and equilibrium-speed laws, then compares those predictions with a mass-conservative numerical method. The simulations approach the predicted waves, while the results connect jamiton formation and stability to uniform-flow stability and identify model-specific breakdown conditions.
5.1 Preliminaries
- The example selects relatively simple pressure and equilibrium-speed functions; the resulting regularized equation has an exact implicit solution, with comparable solutions also available for p ∝ ρ.The shock-speed relation can require numerical solution; for the chosen pressure law, the related polynomial has at most two physically relevant roots.Source: §5.1 · pp. 14–16
5.2 Numerical method
- The Lagrangian particle method conserves mass by construction and remains macroscopic, with typically two orders of magnitude more particles than vehicles.A small amount of numerical viscosity smooths bounded but unphysical oscillations where density or speed changes abruptly.Source: §5.2 · pp. 16–17
5.3 Results
- For λ=500 m, asymptotic numerical profiles agree favorably with theoretical jamitons at matching vehicle counts, except that the scheme smears the shock.Figure 3 compares profiles for N=27 and N=38 after evolution from a small perturbation of an unstable uniform state.Source: §5.3 · p. 18, Figure 3 · p. 19
- The passages predict model breakdown when ρ+ > ρM for ρavg/ρM > 0.277, and when u+ < 0 for ρavg/ρM > 0.391.Source: §5.3 · p. 19
- The example uses a 500 m periodic track, selected pressure and speed laws, and a deliberately low β; larger β may delay, but not necessarily avoid, collisions.For the parameters germane to Figure 4, collisions are anticipated when mean vehicle spacing falls below approximately 3.6ℓ.Source: §5.3 · pp. 18–19
- Multimodal waves often coalesce into one; unimodal waves appear stable to small perturbations except on circuits longer than approximately 2 km, where new waves may grow.On longer circuits, the downstream state can exceed the linear-stability density requirement, allowing further waves to develop.Source: §5.3 · p. 20
- At the stability boundary, jamiton amplitude vanishes as the uniform state becomes less unstable, while coalescing u1 and u2 precludes a self-sustained shock.The results identify jamitons as possible when the corresponding uniform state becomes unstable.Source: §5.3 · p. 21
6 Conclusions
The conclusions establish jamitons as self-sustained traffic shocks linked to detonation-wave structure, and identify which shocks are allowable and when density spikes can reach collision conditions. The appendix also gives exact solutions for specified pressure laws.
6 Conclusions
- The detonation analogy hinges on a sonic point, which permits self-sustained waves unlike inert gas-dynamical shocks.The paper interprets jamitons as saturated phantom jams arising from small instabilities in relatively dense traffic.Source: §6 · pp. 21–22
- For realistic pressure and equilibrium-speed functions, Lax entropy conditions allow only self-sustained shocks overtaken by vehicles; analytical structures for γ = 1, 2 agree excellently with simulations.The analytical cases use ũ = ũ0(1 − ρ/ρM) and p ∝ ργ.Source: §6 · p. 22
- Above a critical average density, the downstream density can exceed ρM, corresponding to collisions; roadway design should avoid triggering such densities anywhere in the domain.Jamitons are described as sharp density increases over short distances and as collision hot spots.Source: §6 · p. 22
- Analytical spatio-temporal stability remains a major future question, with transonic solutions making it difficult to describe.The authors identify stability analysis as potentially useful for examining multiple lower-density waves versus a single density spike.Source: §6 · p. 22
A Exact solutions when p ∝ρ
- For p ∝ ρ, the appendix gives corresponding exact expressions, while traffic speed remains an implicit function of η.The preceding equations provide exact analytical solutions for u, λ, and N when p ∝ ρ².Source: A Exact solutions when p ∝ρ · p. 23
B Model behavior for ρ > ρM
The model’s limiting behavior links negative forcing to densities above the maximum, while its collision condition is reached before the condition for nonsensical speeds.
B Model behavior for ρ > ρM
- The collision condition ρ+ = ρM is necessarily reached before the nonsensical-speed condition u+ < 0.As u+ approaches zero, the mass-flux relation implies s < 0.Source: B Model behavior for ρ > ρM · p. 23
- The jamiton existence condition requires uM, the speed at maximum density, to remain positive.As u+ approaches zero, the shock speed s is negative.Source: B Model behavior for ρ > ρM · p. 23
- For the Lighthill-Whitham-Richards forcing, ũ < 0 requires ρ > ρM, a regime already physically unrealistic.Setting ũ = 0 above ρM would correct the model equations in that regime, but does not alter the physically relevant behavior described.Source: B Model behavior for ρ > ρM · p. 23
C Linear stability of the Payne-Whitham model considered in § 5
The appendix linearizes the Payne–Whitham base state using normal-mode perturbations, then derives a growth-rate condition for linear stability.
C Linear stability of the Payne-Whitham model considered in § 5
- The analysis starts from the constant base state for forcing (5.1) and introduces hatted perturbations expressed as normal modes.The horizontal wave number is k and the corresponding growth rate is σ.Source: C Linear stability of the… · p. 24
- Setting the perturbation matrix determinant to zero and eliminating Λ2 yields a polynomial expression.The derivation introduces Υ2 and an identity for Λ2² along the way.Source: C Linear stability of the… · pp. 24–25
- Linear stability requires a non-positive growth rate, and the stated condition completes the derivation of (5.22).The condition follows from (C.10).Source: C Linear stability of the… · p. 25