Microscopic Models under a Macroscopic Perspective

Transcription

1 Mathematical Models of Traffic Flow, Luminy, France, October 27 Microscopic Models under a Macroscopic Perspective Ingenuin Gasser Department of Mathematics University of Hamburg, Germany Tilman Seidel, Gabriele Sirito, Bodo Werner

2 Outline: General Dynamics of the microscopic model (homogeneous case) Dynamics of the microscopic model (non-homogeneous case and roadworks) Macroscopic view Fundamental diagrams Future Basic concept: Take a very simple microscopic model (Bando), study the full dynamics, take a macroscopic view on the results.

3 Microscopic Bando model on a circular road (scaled) N cars on a circular road of lenght L: Behaviour: x j position of the j-th car ẍ j (t) = { V ( x j+ (t) x j (t) ) ẋ j (t) }, j =,..., N, x N+ = x +L V = V (x) optimal velocity function: V () =, V strictly monoton increasing, lim x V (x) = V max

4 System for the headways: y j = x j+ x j ẏ j = z j ż j = { V (y j+ ) V (y j ) ż j }, j =,..., N, yn+ = y Additional condition: N j= y j = L quasistationary solutions: y s;j = L N, z s;j =, j =,..., N. Linear stability-analysis around this solution gives for the Eigenvalues λ: (λ 2 + λ + β) N β N =, β = V ( L N ) Result (Huijberts ( 2)): For For +cos 2π N +cos 2π N > β max = max x V (x) asymptotic stability = V ( N L ) loss of stability

5 What kind of loss of stability? (I.G., G.Sirito, B. Werner 4): Eigenvalues as functions of β = V ( L N ) k=4 k=3 k=2 B k= Imaginary axis M C A A O Beta 2 3 k= k=2 k=3 B k= Real Axis L/N Bifurcation analysis gives a Hopfbifurcation. Therefore we have locally periodic solutions. Are these solutions stable? (i.e. is the bifurcation sub- or supercritical?)

6 Criterion: Sign of the first Ljapunov-coefficient l Theorem: ( ( ) l = c 2 V L V ( L N )) 2 N V ( N L ) Conclusion: For the mostly used (Bando et al (95)) V (x) = V maxtanh(a(x )) + tanh a + tanh a the bifurction is supercritical (i.e. stable periodic orbits). But: Similar functions V give also subcritical bifurcations.

7 Problem: It seems to be very sensitive with respect to V Global bifurcation analysis: numerical tool (AUTO2) Norm of solution H 3 H Norm of solution Norm of solution H 2 H H L L L Conclusion: Globally similar functions V give similar behaviour. The bifurcation is macroscopically subcritical Conclusion for the application: the critical parameters from the linear analysis are not relevant

8 Example Example 2 Example Norm of solution H Norm of solution H Norm of solution H L L L Relative Velocity Relative Velocity.5.5 Relative Velocity Headway Headway Headway

9 More bifurcations: Eigenvalues as functions of β = V ( L N ) k=4 k=3 k=2 B k= Imaginary axis M C A A O 2 3 k= k=2 k=3 B k= Real Axis Conclusion: There are many other (weakly unstable) periodic solutions

10 (J.Greenberg 4, 7) Solutions with many oscillations finally tend to a solution with one oscillation (G. Oroz, R.E. Wilson, B.Krauskopf 4, 5) Qualitatively the same global bifurcation diagram for the model with delay

11 Extension to standart microscopic model Every driver is aggressive with weight α ẍ j (t) = α { ( V xj+ (t) x j (t) ) ẋ j (t) } +α { ẋ j+ (t) ẋ j (t) }, τ j =,..., N, x N+ = x + L optimal velocity-function: loss of stability similar Aggressive drivers stabilize the fraffic flow! unfortunately also the number of accidents increases! (Olmos & Munos, Condensed matter 24)

12

13 V max ǫ Road works Symmetry breaking, the above theory is not easily applicable A solution is called ponies on a Merry-Go-Round solution (short POM), if there is a T R, such that (i) x i (t + T) = x i (t) + L (i =,..., N) (ii) x i (t) = x i ( t + T N ) (i =,..., N) hold (Aronson, Golubitsky, Mallet-Paret 9). We call T rotation number and N T the phase (phase shift).

14 Theorem:The above model has POM solutions for small ǫ > ẋ x Velocity of the quasistationary solution (no roadwork) versus roadwork solution (The red line indicates maximum velocity).

15 Technique: Poincare maps Π(η) = Φ T(η) (η) Λ, where Φ is the induced flow and Λ reduces the spacial components by L. Π(ξ) Λ ξ x x + Λ Σ Φ T(ǫ,ξ) Σ + Λ T Study fixed points of the corresponding Poincare and reduced Poincare maps. Roadworks are (regular) perturbations. x ǫ

16 Bifurcation diagram in (L, ǫ)-plane: I ε II L A curve of Neimark-Sacker bifurcations in the (L, ǫ)-plane for N = 5.

17 Four different attractors.: ǫ I II ǫ = trivial POM x Hopf periodic solution ǫ > POM x ǫ quasi-pom i.e. here POM s are typically perturbed quasistationary solutions quasi-pom s are perturbed (Hopf) periodic solutions

18 Invariant curves:.8.8 velocity v(ρ).6.4 velocity v(ρ) headway /ρ headway /ρ Two closed invariant curves (ǫ = and ǫ > ) of the reduced Poincaré map π. On the left also the optimal velocity function V is given in gray.

19 The 4 different scenarios: v.5 v x mod L L x L v.5 v x mod L L x L above: no roadworks, below: with roadworks

20 Macroscopic view of the 4 different scenarios: t t x mod L 2 6 x mod L t t x mod L x mod L 8 2 above: no roadworks, below: with roadworks

21 Macroscopic view of density and velocity: t t x mod L x mod L strong road work influence (ǫ =.32) 25

22 Fundamental diagrams: A real world point of view on the reduced Poincaré map π for N =, ǫ =.

23 Fundamental diagrams I: ρ v(ρ) ρ ρ v(ρ) ρ velocity v(ρ) headway /ρ Overlapped fundamental diagrams for N =, L = 5,..., 4 measuring at a fixed point.

24 Fundamental diagrams II: ρ v(ρ) ρ v(ρ) N/L N/L Fundamental diagram of time-averaged flow versus average density for N =, L =

25 Fundamental diagrams III (with roadworks): traffic flow ρ i v i density ρ i Overlapped fundamental diagrams for N =, L = 5,...,4, ǫ =. measuring at a fixed point.

26 Current and future work: Is this dynamics contained in macroscopic models? Which marcoscopic model has the same (rich) dynamics than the basic Bando model Micro-macro link (Aw, Klar, Materne, Rascle 22)