Three-detector problem and Newell's method

The Three-detector problem[1] is a problem in traffic flow theory. Given is a homogeneous freeway and the vehicle counts at two detector stations. We seek the vehicle counts at some intermediate location. The method can be applied to incident detection and diagnosis by comparing the observed and predicted data, so a realistic solution to this problem is important. Newell G.F.[2][3][4] proposed a simple method to solve this problem. In Newell's method, one gets the cumulative count curve (N-curve) of any intermediate location just by shifting the N-curves of the upstream and downstream detectors. Newell's method was developed before the variational theory of traffic flow was proposed to deal systematically with vehicle counts.[5][6][7] This article shows how Newell's method fits in the context of variational theory.

A special case to demonstrate Newell's method

Assumption. In this special case, we use the Triangular Fundamental Diagram (TFD) with three parameters: free flow speed v_f, wave velocity -w and maximum density k_j (see Figure 1). Additionally, we will consider a long study period where traffic past upstream detector (U) is unrestricted and traffic past downstream detector (D) is restricted so that waves from both boundaries point into the (t,x) solution space (see Figure 2).

The goal of three-detector problem is calculating the vehicle at a generic point (P) on the "world line" of detector M (See Figure 2). Upstream. Since the upstream state is uncongested, there must be a characteristic with slope v_f that reaches P from the upstream detector. Such a wave must be emitted \tau_1=L_U/v_f times unit earlier, at point P' on the figure. Since the vehicle number does not change along this characteristic, we see that the vehicle number at the M-detector calculated from conditions upstream is the same as that observed at the upstream detector \tau_1 time units earlier. Since \tau_1 is independent of the traffic state (it is a constant), this result is equivalent to shifting the smoothed N-curve of the upstream detector (curve U of Figure 3) to the right by an amount \tau_1.

Downstream. Likewise, since the state over the downstream detector is queued, there will be a wave reaching P from a location P_2 with wave velocity -w<0. The change in vehicular label along this characteristic can be obtained from the moving observer construction of Figure 4, for an observer moving with the wave. In our particular case, the slanted line corresponding to the observer is parallel to the congested part of TFD. This means that the observer flow is independent of the traffic state and takes on the value: k_j(-w). Therefore, in the time that it takes for the wave to reach the middle location, \tau_2=L_D/(-w), the change in count is \delta=k_j(-w)\tau_2=k_jL_D; i.e., the change in count equals the number of vehicles that fit between M and D at jam density. This result is equivalent to shifting the D-curve to the right \tau_2 units and up \delta units.

Actual count at M. In view of the Newell-Luke Minimum Principle, we see that the actual count at M should be the lower envelope of the U'- and D'-curves. This is the dark curves, M(t). The intersections of the U'- and D'- curves denote the shock's passages over the detector; i.e., the times when transitions between queued and unqueued states take place as the queue advances and recedes over the middle detector. The area between the U'- and M-curves is the delay experienced upstream of location M, trip times are the horizontal separation between curves U(t), M(t) and D(t), accumulation is given by vertical separations, etc.

Mathematical expression. In terms of the function N(t,x) and the detector location (x_u, x_m, x_d) as follows:

 
    N(t,x_m)= \min\{\ N(t-L_U/v_f,x_u)\ ,\ N(t+L_D/w,x_d)+k_jL_D\ \}  \qquad (1)

where L_U=x_m-x_u and L_D=x_d-x_m.

Basic principles of variational theory (VT)

Goal. Suppose we know the number of vehicles (N) along a boundary in a time-space region and we are looking for the number of vehicles at a generic point P (denoted as N_p) beyond that boundary in the direction of increasing time(see Figure 5).[8]

Suppose, again, that an observer starts moving from the boundary to point P along path L. We know the vehicle number the observer sees, N_L. We then break the path of the observer into small sections (such as the one show between A and B) and note that we also know the maximum number of vehicles that can pass the observer along that small section is, C_{AB}. The relative capacity formula tells us that it is: C_{AB}=r(v^0)\Delta{t}. For TFD and using v_{AB} for the slope of segment AB, C_{AB} can be written as:


    C_{AB}=r(v_{AB})\Delta{t}=q_0\Delta{t}-k_0\Delta{x}=q_0(t_B-t_A)-k_0(x_B-x_A); for\ v_{AB}\in[-w,v_f] \qquad (2)

So, if we now add the vehicle number on the boundary to the sum of all C_{AB} along path L we get an upper bound for N_p. This upper bound applies to any observer that moves with speeds in the range [-w,v_f]. Thus we can write:


    N_P \le N_L + \sum_L(C_{AB}),\    v_{AB} \in [-w,v_f] \qquad (3)

Equations (1) and (2) are based on the relative capacity constraint which itself follows from the conservation law.

Maximum principle. It states that N_P is the largest possible value, subject to the capacity constraints. Thus the VT recipe is:


     N_P =\min_L\{N_L+\sum_L(C_{AB})\}\qquad (4)

Equation (4) is a shortest path(i.e., calculus of variations) problem with C_{AB}=r(v_{AB}) as the cost function. It turns out that it produces the same solution as Kinematic wave theory.

Generalized solution

 Three steps:
 1. Find the minimum upstream count, N_U
 2. Find the minimum downstream count, N_D
 3. Choose the lower of the two, N_P=\min\{N_U\ ,\ N_D\}

Step 1

All possible observer straight lines between the upstream boundary and point P have to constructed with observer speeds smaller than free flow speed:


   C_{QP}=q_0\Delta{t}-k_0\Delta{x}=q_0(t_P-t_Q)-k_0(x_P-x_Q) \qquad (5)

where \Delta{t}=t_P-t_Q=\frac{(x_M-x_U)}{v_{QP}} for v_{QP}\in[0,v_f] and \Delta{x}=x_M-x_U

Thus we need to minimize {N_Q+q_0(t_P-t_Q)-k_0(x_M-x_U)}; i.e.,


   N_U=\min_{t_Q}\{N_Q+q_0(t_P-t_Q)-k_0(x_M-x_U)\}  \qquad (6)

Since dN_Q/dt \le q_0, we see that the objective function is non-increasing and therefore t_Q^*=t_{P_1}. So Q should be placed at P_1 and we have:


C_{QP}=C_{P_1P}=q_0\left(\frac{x_M-x_U}{v_f}\right)-k_0(x_M-x_U)=0   \qquad (7)

Thus, N_U=N_{P_1}

Step 2

We have:N_D=\min_{Q^'}\{N_{Q^'}+C_{Q^'P}\}=\min_{Q^'}\{N_{Q^'}+q_0\Delta{t}-k_0\Delta{x}\} So repeat the same steps we find that N_D is minimized when Q^'=P_2. And at point P_2 we get:


  N_D=N_{P_2}+q_0(\frac{x_D-x_M}{w})-k_0(x_D-x_M)  \qquad (8)

Since the FD is triangular, \frac{q_0}{w}+k_0=k_j. Therefore (8) reduces to:


   N_D=N_{P_2}+(x_D-x_M)k_j   \qquad (9)

Step 3

To get the solution we now choose the lower of N_U and N_D.


N_P=\min\{N_U\ ,\ N_D\}=\min\{N_{P_1}\ ,\ N_{P_2}+(x_D-x_M)k_j\}  \qquad (10)

This is Newell's the recipe for the 3-detector problem.

See also

References

  1. Daganzo, Carlos. 1997. Fundamentals of transportation and traffic operations. Oxford: Pergamon.
  2. Newell, G. F. 1993. "A simplified theory of kinematic waves in highway traffic. Part I, General theory". Transportation Research. Part B, Methodological. 27B (4).
  3. Newell, G. F. 1993. "A simplified theory of kinematic waves in highway traffic. Part II. Queuing at freeway bottlenecks". Transportation Research. Part B, Methodological. 27B (4).
  4. Newell, G. F. 1993. "A simplified theory of kinematic waves in highway traffic. Part III. Multi-destination flows". Transportation Research. Part B, Methodological. 27B (4).
  5. Daganzo, Carlos F. 2005. "A variational formulation of kinematic waves: solution methods". Transportation Research. Part B, Methodological. 39B (10).
  6. Daganzo, Carlos F. 2005. "A variational formulation of kinematic waves: basic theory and complex boundary conditions". Transportation Research. Part B, Methodological. 39B (2).
  7. Daganzo, Carlos F. 2006. "On the variational theory of traffic flow: well-posedness, duality and applications". Networks and Heterogeneous Media. 1 (4).
  8. Daganzo, Carlos F. Lecture notes: Operation of transportation facilities. Compiled by Offer Grembek