Title: Cellular Automata Modelling of Traffic in Human and Biological Systems
1Cellular Automata Modelling of Traffic in Human
and Biological Systems
Andreas Schadschneider Institute for Theoretical
Physics University of Cologne
www.thp.uni-koeln.de/as
www.thp.uni-koeln.de/ant-traffic
2Introduction
- Modelling of transport problems
- space, time, states can be discrete or continuous
- various model classes
3Overview
- Highway traffic
- Traffic on ant trails
- Pedestrian dynamics
- Intracellular transport
- Unified description!?!
4Cellular Automata
- Cellular automata (CA) are discrete in
- space
- time
- state variable (e.g. occupancy, velocity)
- Advantage very efficient implementation for
large-scale computer simulations - often stochastic dynamics
5Asymmetric Simple Exclusion Process
6Asymmetric Simple Exclusion Process
Caricature of traffic
- Asymmetric Simple Exclusion Process (ASEP)
- directed motion
- exclusion (1 particle per site)
For applications different modifications
necessary
7Influence of Boundary Conditions
Applications Protein synthesis
Surface growth
Boundary induced phase transitions
exactly solvable!
8Phase Diagram
Maximal current phase JJ(p)
Low-density phase JJ(p,?)
High-density phase JJ(p,?)
9 Highway Traffic
10Cellular Automata Models
- Discrete in
- Space
- Time
- State variables (velocity)
velocity
11Update Rules
- Rules (Nagel-Schreckenberg 1992)
- Acceleration vj ! min (vj 1, vmax)
- Braking vj ! min ( vj , dj)
- Randomization vj ! vj 1 (with
probability p) - Motion xj ! xj vj
(dj empty cells in front of car j)
12Example
Configuration at time t Acceleration (vmax
2) Braking Randomization (p 1/3) Motion
(state at time t1)
13Interpretation of the Rules
- Acceleration Drivers want to move as fast as
possible (or allowed) - Braking no accidents
- Randomization
- a) overreactions at braking
- b) delayed acceleration
- c) psychological effects (fluctuations in
driving) - d) road conditions
- 4) Driving Motion of cars
14Simulation of NaSch Model
- Reproduces structure of traffic on highways
- - Fundamental diagram
- - Spontaneous jam formation
- Minimal model all 4 rules are needed
- Order of rules important
- Simple as traffic model, but rather complex as
stochastic model
15Fundamental Diagram
Relation current (flow) density
16Metastable States
- Empirical results Existence of
- metastable high-flow states
- hysteresis
17VDR Model
- Modified NaSch model
- VDR model
(velocity-dependent randomization) - Step 0 determine randomization pp(v(t))
- p0 if v 0
- p(v) with
p0 gt p - p if v gt 0
- Slow-to-start rule
18Simulation of VDR Model
NaSch model
VDR model
VDR-model phase separation Jam stabilized by
Jout lt Jmax
19 Dynamics on Ant Trails
20Ant trails
ants build road networks trail system
21Chemotaxis
- Ants can communicate on a chemical basis
- chemotaxis
- Ants create a chemical trace of pheromones
- trace can be smelled by other
- ants follow trace to food source etc.
22Ant trail model
- motion of ants
- pheromone update (creation evaporation)
Dynamics
q
q
Q
f f f
parameters q lt Q, f
23Fundamental diagram of ant trails
velocity vs. density
non-monotonicity at small evaporation rates!!
Experiments Burd et al. (2002, 2005)
different from highway traffic no egoism
24Spatio-temporal organization
- formation of loose clusters
early times
steady state
coarsening dynamics
25 Pedestrian Dynamics
26Collective Effects
- jamming/clogging at exits
- lane formation
- flow oscillations at bottlenecks
- structures in intersecting flows ( D.
Helbing)
27Pedestrian Dynamics
- More complex than highway traffic
- motion is 2-dimensional
- counterflow
- interaction longer-ranged (not only nearest
neighbours) -
28Pedestrian model
idea Virtual chemotaxis chemical trace
long-ranged interactions are translated into
local interactions with memory
- Modifications of ant trail model necessary since
- motion 2-dimensional
- diffusion of pheromones
- strength of trace
29Floor field cellular automaton
- Floor field CA stochastic model, defined by
transition probabilities, only local interactions - reproduces known collective effects (e.g. lane
formation)
Interaction virtual chemotaxis (not
measurable!)
dynamic static floor fields interaction with
pedestrians and infrastructure
30Transition Probabilities
- Stochastic motion, defined by
- transition probabilities
- 3 contributions
- Desired direction of motion
- Reaction to motion of other pedestrians
- Reaction to geometry (walls, exits etc.)
- Unified description of these 3 components
31Transition Probabilities
- Total transition probability pij in direction
(i,j) - pij N Mij exp(kDDij)
exp(kSSij)(1-nij) - Mij matrix of preferences (preferred
direction) - Dij dynamic floor field
(interaction between pedestrians) - Sij static floor field
(interaction with geometry) - kD, kS coupling strength
- N normalization (? pij 1)
32Lane Formation
velocity profile
33 Intracellular Transport
34Intracellular Transport
- Transport in cells
- microtubule highway
- molecular motor (proteins) trucks
- ATP fuel
35Kinesin and Dynein Cytoskeletal motors
Fuel ATP
ATP ADP P
Kinesin
Dynein
- Several motors running on same track
simultaneously - Size of the cargo gtgt Size of the motor
- Collective spatio-temporal organization ?
36Practical importance in bio-medical research
Goldstein, Aridor, Hannan, Hirokawa,
Takemura,.
37ASEP-like Model of Molecular Motor-Traffic
ASEP Langmuir-like adsorption-desorption
Parmeggiani, Franosch and Frey, Phys. Rev. Lett.
90, 086601 (2003)
D
A
q
a
b
Also, Evans, Juhasz and Santen, Phys. Rev.E. 68,
026117 (2003)
38Spatial organization of KIF1A motors experiment
MT (Green)
10 pM
KIF1A (Red)
100 pM
1000pM
2 mM of ATP
2 mm
position of domain wall can be measured as a
function of controllable parameters.
Nishinari, Okada, Schadschneider, Chowdhury,
Phys. Rev. Lett. (2005)
39Summary
- Various very different transport and traffic
problems can be described by similar models - Variants of the Asymmetric Simple
Exclusion Process - Highway traffic larger velocities
- Ant trails state-dependent hopping rates
- Pedestrian dynamics 2d motion, virtual
chemotaxis - Intracellular transport adsorption desorption
40Collaborators
Thanx to
- Duisburg
- Michael Schreckenberg
- Robert Barlovic
- Wolfgang Knospe
- Hubert Klüpfel
Rest of the World Debashish Chowdhury
(Kanpur) Ambarish Kunwar (Kanpur) Katsuhiro
Nishinari (Tokyo) T. Okada (Tokyo)
- Cologne
- Ludger Santen
- Alireza Namazi
- Alexander John
- Philip Greulich
many others