An Equilibrium Dynamic Traffic Assignment Model with Linear Programming Formulation

Static traffic assignment. Static models assume that network flows do not depend on the time variable and thus can only capture the average characteristics of the network’s loading for a given period of time. The Beckmann model [6] with Bureau of Public Roads (BPR) functions representing the dependency between flow density and travel time on a road segment is probably the most popular model [7, 27]. Paper [25] presented the alternative Stable Dynamics static model, where the link flows are strictly limited by the capacities, motivated by the logical inconsistencies of the Beckmann model.

A significant advantage of static models is their computational efficiency. Static equilibrium models can be formulated as convex optimization problems [6] or variational inequalities if flow interactions are considered [10, 28]. Many highly performant algorithms were designed specifically for such problems [24, 5, 18, 15]. This allows for the development of combined travel demand models based on static models, which can still be solved within a reasonable time limit for decent accuracy in large networks [13, 19].

Dynamic traffic assignment. First dynamic traffic assignment models were suggested in [22, 23]. However, they were limited to searching only system optimum, which means that each single driver wants to reduce not his personal cost, but total expenses of all drivers in a system, according to Wardrop’s second principle. These first models were limited to single destination and single commodity due to computational complexity of dynamic models. In the work of Carey [9] it was shown that when models consider multicommodities and multiple destinations then we should also think about first-in-first-out (FIFO) condition. Independently of what type of equilibrium we search – user or system – there is non-convexity arise due to FIFO constraint. Author of the paper considers FIFO condition on average. This means that in reality some particular vehicles can pass each other and travel at different speed, but on average traffic which enter a road first will exit from it first.

Similar idea is highlighted in [32]. The authors conclude that FIFO is not a necessary constraint for building dynamic traffic assignment problem, however FIFO assumption can simplify building link and traffic dynamic models.

The study and analysis of traffic congestion during rush hours have a long-standing history, beginning with the pioneering work of Vickrey in 1969 [30]. In Vickrey’s fundamental model, it is assumed that a set of commuters, aim to reach a single destination using a single route that features a bottleneck with limited capacity at the desired time, which is evenly distributed over a closed interval. Each commuter decides on their departure time from home to minimize their total travel cost, which includes travel time, waiting time at the bottleneck, and the costs of arriving early or late at their destination. Several authors have extended this model. For example, [29] investigated the existence of equilibrium for a wider class of arrival cost functions. The follow-up work [11] described additional properties of the equilibrium in the model of [29]. Papers [4, 21] consider heterogeneities in the cost of travel time or delay in the schedule. There are also works (for example [11]) in which different types of user heterogeneities are assumed (work start times and value of times), but with a very stringent assumption about the piecewise linearity of the schedule function. However, in [26], the authors showed that there is no guaranteed convergence in the straightforward formulation of the equilibrium conditions. As noted earlier, all these works consider the simplest network design — the case of a single origin-destination pair connected by one route with a single bottleneck. In the articles [20] and [3] the authors decided to systematize the models of departure-time choice equilibrium (DTCE) for this type of road network; in the second work the authors additionally proposed a novel approach to modeling departure time equilibrium. They utilize continuous-time linear programming and optimal transport theory to achieve analytical solutions. Special attention is given to the Monge property, which determines sorting patterns of equilibrium, allowing models to account for multiple types of user heterogeneity.

There are also many studies focusing on other types of road networks. The paper by [2] considered a road system with several bottlenecks, while in [16], the authors examined a network with multiple origin-destination pairs and several roads between them, each containing a single bottleneck.

In Section 2 we formulate an equilibrium departure time and route choice problem in continuous time, and provide a toy example which can be solved analytically. Then, in Section 3 we derive a formulation of a discrete variant of the problem as a minimum cost multicommodity flow problem. We show that the solution to the obtained optimization problem is a user equilibrium network loading. Last, in Section 4 we discuss numeric algorithms suitable for solving the discretized problem.

Transportation Network. We consider the road network’s representation as a simple directed graph, where nodes correspond to junctions and centroids (artificial origin and destination nodes) and links correspond to road segments and centroid connectors.

We would like to represent the entire transportation network as a simple directed graph, where links are associated with capacities $\hat{x}_{e}$ and free-flow times $\bar{t}_{e}$, without offloading any additional information, such as allowed turn directions, to other data structures. This can be achieved by splitting each node corresponding to a junction in the following way. For each input link $(v_{1},v_{2})$ adjacent to the intersection node $v_{2}$, we replace this node in the link with an artificial input node $v_{12}$. Similarly, we replace the first node in an output link $(v_{2},v_{3})$ with a node $v_{23}$. The splitting procedure is illustrated at Fig. 1. If the turn $(v_{1},v_{2},v_{3})$ is allowed at the junction, we add a link $(v_{12},v_{23})$ with capacity and free-flow time representing the turn capacity and the time required to make the turn. This results in that we have three disjoint subsets of links: $E=E_{R}\cup E_{J}\cup E_{C}$, where $E_{R}$ consists of links corresponding to road segments such as $(v_{1},v_{12})$, $E_{J}$ consists of links corresponding to junctions such as $(v_{12},v_{23})$ and $E_{C}$ are centroid connectors. The resulting set of all nodes is denoted by $V$, and the whole transportation graph is denoted by $G=(V,E)$. Sets of origin and destination centroids are denoted by $I\subset V$ and $J\subset V$ respectively.

Travel Cost. We assume that the travel cost, that each user strives to selfishly minimize, sums up from the travel time and a cost $\tau(t,T_{0})$ which is a function of the actual arrival time $t$ parametrized by the desired arrival time $T_{0}$ (desired arrival times an be different for different users). The arrival time cost represents the fact that most users aims to appear at some place at some specific time, and arriving late or too early is undesirable. The travel time sums from free-flow travel times of traversed links $e\in E$ denoted by $\bar{t}_{e}$ and time spent in queues.

This paper adopts the point queue link model to represent the queuing effect of congestion. In the point queue model, queue hold no physical space: vehicles that cannot live the link just accumulate at the end of the link and do not interfere with the other vehicles, which travels freely to the end of the link.

We assume that all road segments have constant flow limit for all of its sections, and all capacity constraints are posed only for junctions. Any change in link’s capacity can be modeled by introducing an artificial junction, possibly with single input and single output. We denote the capacity (in standardized vehicles per unit of time) of link $e\in E$ by $\hat{x}_{e}$, setting it to infinity for edges in $E_{R}$.

Demands. Traffic between origins and destinations is represented by the time-dependent matrix $\{d_{ij}\}_{i\in I,j\in J}$, where $d_{ij}(t)$ is the total number/density of vehicles in zone $i$ that aim to arrive to zone $j$ at time $t$.

To provide intuition about dynamic models in continuous time here we derive an equilibrium solution for a simple single origin-single destination network. This example was shared with us by Yurii Nesterov. A similar model differing in assumption about the distribution of arrival time preferences was also proposed in [30].

In the morning, $N$ cars must leave a residential area to reach a work area. Each driver wants to arrive to the work area precisely at $T_{0}=$ 9 a.m.. If a driver arrives after $T_{0}$, he may be late for work. Conversely, arriving before $T_{0}$ results in lost sleep and waiting time before starting the workday. We represent the costs of arriving later or earlier in the same units as travel costs. Costs of arriving later or earlier are proportional to the time difference with coefficients $\alpha$ and $\beta$ respectively: $\tau(t,T_{0})=\alpha\max\left\{0,t-T_{0}\right\}+\beta\max\left\{0,T_{0}-t\right\}$.

There is a single road between the residential and the work areas. The usual travel time on a free road is $\bar{t}$ minutes. However, there is a narrow section midway along the road where the throughput is limited to $\hat{x}$ cars per hour.

We need to find the equilibrium distribution of cars depending on departure times from the residential area $n(t)$.

It is impossible for all cars to arrive to the destination area at the same time. The driver arriving at exactly $T_{0}$ will be called a punctual driver. Let us denote his departure time as $t_{2}$. Let also denote the time when the first car leaves the residential area by $t_{1}$. The last car will depart at $t_{3}$.

If we want to thoroughly describe the car distribution we need to find the densities of our two flows: before and after the punctual driver $n_{1}(t)$ and $n_{2}(t)$, as well as the time borders limiting those distributions ($t_{1}$, $t_{2}$, $t_{3}$).

First, let’s find the departure times of the first driver $t_{1}$ and last driver $t_{3}$. The first and the last driver don’t need to to be stuck at the traffic jam, because they can move in the very beginning of the traffic flow when the jam had not formed yet or move in the end when the jam has already finished. According to our task the travellers costs depend on the time difference between the desired and actual arrival times to the destination area. The arrival time for the first driver can be expressed as ($t_{1}+\bar{t}$), hence for the last driver it will be ($t_{3}+\bar{t}$). The time difference between the desired and actual arrival time for the first driver ($T_{0}-(t_{1}+\bar{t})$), and for the last driver ($(t_{3}+\bar{t})-T_{0}$).

According to the first Wardrop principle, which is an application of the Nash’s rule to the transportation task, the costs should be the same for all drivers or someone will change his strategy in order to optimise his costs. In our case the strategy means the departure time that the drivers choose to start their journeys. The equilibrium requires the same costs for the drivers. We can write an equation balancing the costs per minute of arriving early $\beta$, and cost per minute of arriving late $\alpha$

Obviously, the time from the first to the last car will be spread among all the other cars to pass through the bottleneck.

Now lets calculate $t_{3}$. The term $\frac{N}{\hat{x}}$ denotes the total time taken by $N$ cars to pass through the bottleneck at a rate of $\hat{x}$ cars per hour

The previous two equations are resolved as a system to obtain $t_{1}$ and $t_{3}$:

As drivers travel, they either arrive late, early, or on time. Let’s find the departure time $t_{2}$ of the driver who arrived on time. We will call this driver punctual. The time for the punctual driver balancing out the costs adjusted for travel costs of the other drivers. As before, we will use the first Wardrop principle and consider that the benefit to each driver is the same

Using previously found $t_{1}$ the solution for $t_{2}$ is:

Now we can determine the distribution of cars that arrive before, after the punctual driver, and at the same time. All the cars that arrive before the punctual driver will form a car flow that value work time higher then rest time, on the other hand the cars that arrive after value their personal time higher than being at work at a certain time mark. This analysis focuses on the specific scenario when the road has a bottleneck with limited capacity.

The model is formulated based on the law of conservation of mass and Kirchhoff’s first law, which asserts that traffic entering a road link must equal the traffic exiting, given no sources or sinks in between. We can find $n_{1}$ and $n_{2}$ which are the traffic flow rates in two successive time intervals before and after the punctual driver and $\hat{x}$ is the capacity of the bottleneck. The number of cars leaving the departure area, moving in the first part of the flow before the punctual driver can be found by multiplying the first traffic density $n_{1}$ by the departure time difference between the first and the punctual driver. The number of cars of the same flow arriving to the destination area can be found by taking the multiplication of the bottleneck throughput by the time taken by the first part of the flow to come to the destination area from the first driver to the punctual driver. When we add the number of cars of our two flows we must get the total number of cars moving from the departure to the destination area. Lets resolve those equations as a system:

Lets find the time which takes to the first part of the flow to leave the departure area ($t_{2}-t_{1}$)

To apply the conservation of mass law lets find the number of cars that will arrive to the destination area as the time needed for all the cars of the first part of the flow to arrive to the destination area times the capacity of the bottle neck

Now we can find the density of the first part of the flow

Lets find the time taken by the second part of the flow to leave the departure area

Adding together our two flows we should end up with a total number of cars. Using that principle allows find the density of the second part of the flow:

Now we can describe our car flow as a set of inequalities for time borders and the densities of each part

Complementarity. By the first Wardrop’s principle, the equilibrium time and route choice problem formulates as the problem of assigning a departure time and a path in the spatial transportation network for each user in such a way that no user can change his strategy (departure time and spatial path) to strictly decrease travel cost. If we embed the departure time in the definition of path, we can say that in equilibrium each user is assigned to a shortest path.

For single spatial path, time spent in queue may depend on the network loading. But in time-expanded graph, each path has fixed travel cost. The way, in which users affect each other choices, from the time-expanded graph point of view, is that some paths may be prohibited to use, because they are already occupied by other agents so that their capacity limit is reached. Therefore, we need to distinguish available and fully occupied paths.

Mathematically, the equilibrium problem can be stated as a nonlinear complementarity problem:

where $P_{ij}$ is the set of all paths from origin node $i\in I$ to destination node $j\in J$ of the $ij$-th demand in the time-expanded graph; $P^{\text{avail}}_{ij}(x)$ is the set of paths $p\in P_{ij}$ such that flow $x_{p}$ can be increased by a $\delta_{p}>0$ without violating capacity constraints (18); $P=\bigcup_{i\in I,j\in J}P_{ij}$; $x$ is the path flows vector with components $x_{p},\,p\in P$; $T_{p}$ is the travel cost for path $p$; The meaning of the inequalities is as follows: (17) ensures that flows are non-negative; (16) stands for either $x_{p}=0$ or the travel cost for path $p$ is less or equal than the cost of the shortest available path for $ij$-th demand (thus the problem is called complementarity problem); and (18),(19) express capacity and demand constraints respectively.

Variational Inequality. Let $X$ denote the set of all feasible path flows vectors, i.e. path flows vectors satisfying (17)-(19). Then, the nonlinear complementarity problem can be reformulated as the following variational inequality problem (VI) on $x$:

Linear Programming. One can note that (16’) coincides with the first-order optimality criteria for a problem of minimizing a convex function $f(x)$ on a convex set $X$ if $\partial f(x)/\partial x_{p}=T_{p}$ [8, Propositon 1.3]. The set $X$ defined above is indeed convex, and the convex function $f(x)=\sum_{p\in P}x_{p}T_{p}$, which is equal to total travel cost, satisfies the equality for partial derivatives. Therefore, the following convex optimization problem

which, in fact, is a linear programming problem, is equivalent to VI (16’). Let us assume that $X$ is nonempty. Due to (17) and (18), $X$ is bounded. This implies that problem (20) has a solution. Together with Lemma 1 this gives us the following theorem.

Problem (20) is defined for the path flows vector variable. The problem can also be reformulated for link flows variables [1]:

where $x_{ei}$ denotes flow on edge $e$ originating from origin $i\in I$; $x_{e}$ denotes total flow on edge $e$; $\text{In}(v)\subset{\mathcal{E}}$ and $\text{Out}(v)\subset{\mathcal{E}}$ denote sets of links coming in and out of node $v\in{\mathcal{V}}$ respectively. Equation (24) imposes flow continuity constraint in form of Kirchhoff’s current law.

In this form MCMF problem has $O(mT|I|)$ variables and $O(nT|I|)$ constraints. Constraint (24) can be written in a matrix form using the incidence matrix of the time-expanded graph. This allows to apply existing linear programming solvers [14, 12], but for larger networks or finer time discretization this approach may not work due to the problem’s size.

The column generation approach [1, Section 17.5] can be applied to solve the MCMF problem in the path-flows form (20). Basically, it is a specific variant of simplex method, which tackle the problem of huge dimension of $P$ by operating only on a limited subset of paths from $P$ which carry nonzero flow. Column generation simplex method at each iteration solves the shortest paths problem for each demand and then solves a linear system of equations to distribute path flows between previously used paths and newly found shortest paths. Dual variable’s values are updated to satisfy a condition that all active paths for given demand have the same travel costs (with respect to the additional costs introduced by the dual variables). The primal variables in the linear system, also called basic variables, are nonzero path flows and nonzero slack variables for constraint (23). The number of basic variables is the number of paths with nonzero flow plus the number of not fully loaded link, and it is equal to $|K|+|{\mathcal{E}}|$. Since $|{\mathcal{E}}|=O(mT)$, the system can be prohibitively large.

However, it is likely that only a small subset of links will carry nonzero flow at each iteration: for every $t\in{\mathcal{T}},~{}u,v\in V$ links $\left((u,t),(v,t^{\prime})\right)$ can be expected to be used only for $t^{\prime}\geq t$ in a small subset of ${\mathcal{T}}$, corresponding to limited variation of waiting time between users in the same queue. After solving the shortest path subproblem, unused links can be removed from consideration for the linear system solving step, greatly reducing number of variables in the system. With this optimization in mind, we can expect each step of the column generation simplex method require solving linear system on $O(m+|K|)$ variables and solving the shortest path problem on network with $O(Tn)$ nodes and $O(Tm)$ links for each demand, what can be done by $|I|$ invocations of Dijkstra’s algorithm.

By using Lagrangian duality to handle the capacity constraints (18) in problem (20), one can obtain the following primal-dual problem [1, Section 17.4]:

where dual variables $y_{e}$ penalize usage of congested links. As in column generation approach, one iteration of a subgradient method for solving the dual problem (by variable $y$) require finding shortest paths for all demands with respect to additional costs introduced with dual variables. Demands are then distributed over the shortest paths what induces the so-called all-or-nothing traffic assignment on links. Comparing to column generation, the linear system solving step is replaced with a much simpler operation: the flow distribution from the previous iteration is averaged with the obtained all-or-nothing assignment. Dual variables values $y_{e}$ are increased for congested links and decreased for links with available capacity (proportionally to the difference $\sum_{p\ni e}x_{p}-\bar{x}_{e}$ in case of classical subgradient method). Again, only variables, corresponding to active (used) links are updated at each iteration.

In [18] a detailed discussion of this approach is given for general MCMF problem (referred to by the alias of the stable dynamic model) with convergence analysis of adaptive primal-dual subgradient methods.

Since practically efficient simplex method has unsatisfying worst-case complexity [17], we cannot make informative analytical comparison of complexities of above mentioned approaches. However, one can expect link-based approach 4.1 only work for small networks, while column generation 4.2 and dual subgradient 4.3 methods can be applied for large networks with the latter have computationally more affordable but less effective iteration.