# Three-detector problem and Newell's method

Last updated

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 ${\displaystyle v_{f}}$, wave velocity -w and maximum density ${\displaystyle 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 ${\displaystyle v_{f}}$ that reaches P from the upstream detector. Such a wave must be emitted ${\displaystyle \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 ${\displaystyle \tau _{1}}$ time units earlier. Since${\displaystyle \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 ${\displaystyle \tau _{1}}$.

Downstream. Likewise, since the state over the downstream detector is queued, there will be a wave reaching P from a location ${\displaystyle P_{2}}$ with wave velocity ${\displaystyle -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: ${\displaystyle k_{j}(-w)}$. Therefore, in the time that it takes for the wave to reach the middle location, ${\displaystyle \tau _{2}=L_{D}/(-w)}$, the change in count is ${\displaystyle \delta =k_{j}(-w)\tau _{2}=k_{j}L_{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 ${\displaystyle \tau _{2}}$ units and up ${\displaystyle \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 (${\displaystyle x_{u}}$, ${\displaystyle x_{m}}$, ${\displaystyle x_{d}}$) as follows:

${\displaystyle N(t,x_{m})=\min\{\ N(t-L_{U}/v_{f},x_{u})\ ,\ N(t+L_{D}/w,x_{d})+k_{j}L_{D}\ \}\qquad (1)}$

where ${\displaystyle L_{U}=x_{m}-x_{u}}$ and ${\displaystyle 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 ${\displaystyle 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, ${\displaystyle 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, ${\displaystyle C_{AB}}$. The relative capacity formula tells us that it is: ${\displaystyle C_{AB}=r(v^{0})\Delta {t}}$. For TFD and using ${\displaystyle v_{AB}}$ for the slope of segment AB, ${\displaystyle C_{AB}}$ can be written as:

${\displaystyle 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 ${\displaystyle C_{AB}}$ along path L we get an upper bound for ${\displaystyle N_{p}}$. This upper bound applies to any observer that moves with speeds in the range ${\displaystyle [-w,v_{f}]}$. Thus we can write:

${\displaystyle N_{P}\leq 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 ${\displaystyle N_{P}}$ is the largest possible value, subject to the capacity constraints. Thus the VT recipe is:

${\displaystyle 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 ${\displaystyle 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, ${\displaystyle N_{U}}$  2. Find the minimum downstream count, ${\displaystyle N_{D}}$  3. Choose the lower of the two, ${\displaystyle N_{P}=\min\{N_{U}\ ,\ N_{D}\}}$

### Step 1

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

${\displaystyle 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 ${\displaystyle \Delta {t}=t_{P}-t_{Q}={\frac {(x_{M}-x_{U})}{v_{QP}}}}$ for ${\displaystyle v_{QP}\in [0,v_{f}]}$ and ${\displaystyle \Delta {x}=x_{M}-x_{U}}$

Thus we need to minimize ${\displaystyle {N_{Q}+q_{0}(t_{P}-t_{Q})-k_{0}(x_{M}-x_{U})}}$; i.e.,

${\displaystyle N_{U}=\min _{t_{Q}}\{N_{Q}+q_{0}(t_{P}-t_{Q})-k_{0}(x_{M}-x_{U})\}\qquad (6)}$

Since ${\displaystyle dN_{Q}/dt\leq q_{0}}$, we see that the objective function is non-increasing and therefore ${\displaystyle t_{Q}^{*}=t_{P_{1}}}$. So Q should be placed at ${\displaystyle P_{1}}$ and we have:

${\displaystyle C_{QP}=C_{P_{1}P}=q_{0}\left({\frac {x_{M}-x_{U}}{v_{f}}}\right)-k_{0}(x_{M}-x_{U})=0\qquad (7)}$

Thus, ${\displaystyle N_{U}=N_{P_{1}}}$

### Step 2

We have:${\displaystyle 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 ${\displaystyle N_{D}}$ is minimized when ${\displaystyle Q^{'}=P_{2}}$. And at point ${\displaystyle P_{2}}$ we get:

${\displaystyle 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, ${\displaystyle {\frac {q_{0}}{w}}+k_{0}=k_{j}}$. Therefore, (8) reduces to:

${\displaystyle 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 ${\displaystyle N_{U}}$ and ${\displaystyle N_{D}}$.

${\displaystyle 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.

The fundamental diagram of traffic flow is a diagram that gives a relation between the traffic flux (vehicles/hour) and the traffic density (vehicles/km). A macroscopic traffic model involving traffic flux, traffic density and velocity forms the basis of the fundamental diagram. It can be used to predict the capability of a road system, or its behaviour when applying inflow regulation or speed limits.

Microscopic traffic flow models are a class of scientific models of vehicular traffic dynamics.

Microsimulation is a category of computerized analytical tools that perform highly detailed analysis of activities such as highway traffic flowing through an intersection, financial transactions, or pathogens spreading disease through a population. Microsimulation is often used to evaluate the effects of proposed interventions before they are implemented in the real world. For example, a traffic microsimulation model could be used to evaluate the effectiveness of lengthening a turn lane at an intersection, and thus help decide whether it is worth spending money on actually lengthening the lane.

## Related Research Articles

The heat equation is a parabolic partial differential equation that describes the distribution of heat in a given region over time.

In mathematics, Green's theorem gives the relationship between a line integral around a simple closed curve C and a double integral over the plane region D bounded by C. It is named after George Green, though its first proof is due to Bernhard Riemann and is the two-dimensional special case of the more general Kelvin–Stokes theorem.

In calculus, the product rule is a formula used to find the derivatives of products of two or more functions. It may be stated as

In the theory of stochastic processes, the Karhunen–Loève theorem, also known as the Kosambi–Karhunen–Loève theorem is a representation of a stochastic process as an infinite linear combination of orthogonal functions, analogous to a Fourier series representation of a function on a bounded interval. The transformation is also known as Hotelling transform and eigenvector transform, and is closely related to principal component analysis (PCA) technique widely used in image processing and in data analysis in many fields.

Virtual work arises in the application of the principle of least action to the study of forces and movement of a mechanical system. The work of a force acting on a particle as it moves along a displacement be different for different displacements. Among all the possible displacements that a particle may follow, called virtual displacements, one will minimize the action. This displacement is therefore the displacement followed by the particle according to the principle of least action. The work of a force on a particle along a virtual displacement is known as the virtual work.

In continuum mechanics, the finite strain theory—also called large strain theory, or large deformation theory—deals with deformations in which strains and/or rotations are large enough to invalidate assumptions inherent in infinitesimal strain theory. In this case, the undeformed and deformed configurations of the continuum are significantly different, requiring a clear distinction between them. This is commonly the case with elastomers, plastically-deforming materials and other fluids and biological soft tissue.

In pharmacology, the clearance is a pharmacokinetic measurement of the volume of plasma from which a substance is completely removed per unit time; the usual units are mL/min. The quantity reflects the rate of drug elimination divided by plasma concentration.

In general relativity, the Gibbons–Hawking–York boundary term is a term that needs to be added to the Einstein–Hilbert action when the underlying spacetime manifold has a boundary.

In mathematics and civil engineering, traffic flow is the study of interactions between travellers and infrastructure, with the aim of understanding and developing an optimal transport network with efficient movement of traffic and minimal traffic congestion problems.

The finite element method (FEM) is a powerful technique originally developed for numerical solution of complex problems in structural mechanics, and it remains the method of choice for complex systems. In the FEM, the structural system is modeled by a set of appropriate finite elements interconnected at discrete points called nodes. Elements may have physical properties such as thickness, coefficient of thermal expansion, density, Young's modulus, shear modulus and Poisson's ratio.

In physics, Hamilton's principle is William Rowan Hamilton's formulation of the principle of stationary action. It states that the dynamics of a physical system is determined by a variational problem for a functional based on a single function, the Lagrangian, which contains all physical information concerning the system and the forces acting on it. The variational problem is equivalent to and allows for the derivation of the differential equations of motion of the physical system. Although formulated originally for classical mechanics, Hamilton's principle also applies to classical fields such as the electromagnetic and gravitational fields, and plays an important role in quantum mechanics, quantum field theory and criticality theories.

In mathematics, the class of Muckenhoupt weightsAp consists of those weights ω for which the Hardy–Littlewood maximal operator is bounded on Lp(). Specifically, we consider functions f on Rn and their associated maximal functions M( f ) defined as

In mathematics, Maass wave forms or Maass forms are studied in the theory of automorphic forms. Maass wave forms are complex-valued smooth functions of the upper half plane, which transform in a similar way under the operation of a discrete subgroup of as modular forms. They are Eigenforms of the hyperbolic Laplace Operator defined on and satisfy certain growth conditions at the cusps of a fundamental domain of . In contrast to the modular forms the Maass wave forms need not be holomorphic. They were studied first by Hans Maass in 1949.

In mathematics, a Bose–Mesner algebra is a special set of matrices which arise from a combinatorial structure known as an association scheme, together with the usual set of rules for combining those matrices, such that they form an associative algebra, or, more precisely, a unitary commutative algebra. Among these rules are:

The fractional Schrödinger equation is a fundamental equation of fractional quantum mechanics. It was discovered by Nick Laskin (1999) as a result of extending the Feynman path integral, from the Brownian-like to Lévy-like quantum mechanical paths. The term fractional Schrödinger equation was coined by Nick Laskin.

Vehicular traffic can be either free or congested. Traffic occurs in time and space, i.e., it is a spatiotemporal process. However, usually traffic can be measured only at some road locations. For efficient traffic control and other intelligent transportation systems, the reconstruction of traffic congestion is necessary at all other road locations at which traffic measurements are not available. Traffic congestion can be reconstructed in space and time based on Boris Kerner’s three-phase traffic theory with the use of the ASDA and FOTO models introduced by Kerner. Kerner's three-phase traffic theory and, respectively, the ASDA/FOTO models are based on some common spatiotemporal features of traffic congestion observed in measured traffic data.

Heat transfer physics describes the kinetics of energy storage, transport, and energy transformation by principal energy carriers: phonons, electrons, fluid particles, and photons. Heat is energy stored in temperature-dependent motion of particles including electrons, atomic nuclei, individual atoms, and molecules. Heat is transferred to and from matter by the principal energy carriers. The state of energy stored within matter, or transported by the carriers, is described by a combination of classical and quantum statistical mechanics. The energy is also transformed (converted) among various carriers. The heat transfer processes are governed by the rates at which various related physical phenomena occur, such as the rate of particle collisions in classical mechanics. These various states and kinetics determine the heat transfer, i.e., the net rate of energy storage or transport. Governing these process from the atomic level to macroscale are the laws of thermodynamics, including conservation of energy.

The acoustoelastic effect is how the sound velocities of an elastic material change if subjected to an initial static stress field. This is a non-linear effect of the constitutive relation between mechanical stress and finite strain in a material of continuous mass. In classical linear elasticity theory small deformations of most elastic materials can be described by a linear relation between the applied stress and the resulting strain. This relationship is commonly known as the generalised Hooke's law. The linear elastic theory involves second order elastic constants and yields constant longitudinal and shear sound velocities in an elastic material, not affected by an applied stress. The acoustoelastic effect on the other hand include higher order expansion of the constitutive relation between the applied stress and resulting strain, which yields longitudinal and shear sound velocities dependent of the stress state of the material. In the limit of an unstressed material the sound velocities of the linear elastic theory are reproduced.

In the ADM formulation of general relativity one splits spacetime into spatial slices and time, the basic variables are taken to be the induced metric, , on the spatial slice, and its conjugate momentum variable related to the extrinsic curvature, ,. These are the metric canonical coordinates.

Ashtekar variables, which were a new canonical formalism of general relativity, raised new hopes for the canonical quantization of general relativity and eventually led to loop quantum gravity. Smolin and others independently discovered that there exists in fact a Lagrangian formulation of the theory by considering the self-dual formulation of the Tetradic Palatini action principle of general relativity. These proofs were given in terms of spinors. A purely tensorial proof of the new variables in terms of triads was given by Goldberg and in terms of tetrads by Henneaux et al. Here we in particular fill in details of the proof of results for self-dual variables not given in text books.

## 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