next up previous
Next: About this document

First approach -- track motion of each car Let tex2html_wrap_inline716 be the location of the j'th car at time t, measured in miles along the road, say. The primary assumption that allows us to track the motion of the cars is that the velocity of the j'th car depends on the distance between this car and the next. In other words, drivers adjust their speed depending on how crowded the road is just in front of them. The distance to the next car is tex2html_wrap_inline724 .

The velocity of the j'th car depends on tex2html_wrap_inline728 in some way, but it is more usual to re-express it in terms of the density of traffic, which might be measured in cars/mile. The symbol tex2html_wrap_inline730 is usually used for density, and can vary between tex2html_wrap_inline732 (empty road) and tex2html_wrap_inline734 where L is the length of a car (in miles!) so that tex2html_wrap_inline734 corresponds to bumper-to-bumper traffic.

Then the density seen by the j'th driver is

displaymath700

Now we will suppose the velocity of traffic is given by some function of the density, tex2html_wrap_inline742 is the velocity in miles/hour. A simple example might be the linear relation shown below. The velocity is 0 when tex2html_wrap_inline734 (bumper-to-bumper) and increases to some maximum value, say 55 miles/hour, as the density decreases.

figure363

The motion of the j'th car is then determined by a differential equation,

displaymath701

This can be approximated by taking discrete time steps of length tex2html_wrap_inline748 and moving each car a little bit over each time step based on the density seen at the corresponding time. Let tex2html_wrap_inline750 denote the approximate location of the j'th car at time tex2html_wrap_inline754 . Then we can update the car positions by

displaymath702

where the density seen by the j'th car at time tex2html_wrap_inline758 is tex2html_wrap_inline760 . This is the approach used to generate the animation seen on the web page.

Second approach -- A finite volume method for density In modeling the flow of air it is not possible to keep track of how every molecule moves, since there are on the order of tex2html_wrap_inline762 molecules in any region of interest. The approach usually used is to split the region up into small pieces (finite volumes) and keep track of the density in each region.

In the case of traffic flow, this amounts to subdividing the highway into short stretches of highway as indicated below:

figure373

Each stretch (grid cell) has length tex2html_wrap_inline778 (in miles) and if tex2html_wrap_inline780 is an estimate of the density in the i'th cell at time tex2html_wrap_inline758 (in cars/mile), then the number of cars in the i'th cell is roughly tex2html_wrap_inline788 .

How many cars will be on this stretch of highway at time tex2html_wrap_inline790 ? Assuming cars don't vanish or appear spontaneously, the number will change only becuase of cars entering the cell at tex2html_wrap_inline792 or leaving at tex2html_wrap_inline794 . The flow of cars past a given point per unit time is called the flux of cars. This can be computed as the product of the density (cars/mile) times the speed (miles/hour) giving a flux in cars/hour:

displaymath703

The flux of cars passing the point tex2html_wrap_inline792 is approximately

displaymath704

based on the density in the cell to the left, from which cars are arriving at tex2html_wrap_inline792 .

Over a short time period tex2html_wrap_inline748 , the number of cars passing tex2html_wrap_inline792 will be approximately tex2html_wrap_inline804 , so the total number of cars in this cell at time tex2html_wrap_inline790 can be approximated by

eqnarray692

This gives a formula for updating the density:

displaymath705

In computing this, we use the relation tex2html_wrap_inline808 to calculate the velocity based on the density in each cell, e.g., tex2html_wrap_inline810 .

Note that we can rewrite this formula as

displaymath706

If we let tex2html_wrap_inline778 and tex2html_wrap_inline748 go to zero, this suggests a partial differential equation:

displaymath707

Partial differential equations of this sort are the basic tool for studying fluid dynamics.

Here's how the density evolves with this sort of finite volume method. (Actually a more complicated method is used which gives ``sharper'' results.) Each circle represents one value tex2html_wrap_inline780 . In this calculation there were 100 cells. Here the density function was rescaled so that tex2html_wrap_inline818 corresponds to bumper-to-bumper traffic, i.e., the density is measured in ``cars per car length''.

Can you figure out how cars are moving from these density plots?

figure410




next up previous
Next: About this document

Randy LeVeque
Sun Apr 19 23:35:52 PDT 1998