The Riemann Solver for Traffic Flow at an Intersection with Buffer of Vanishing Size

The paper examines the model of traffic flow at an intersection introduced in [2], containing a buffer with limited size. As the size of the buffer approach zero, it is proved that the solution of the Riemann problem with buffer converges to a self-similar solution described by a specific Limit Riemann Solver (LRS). Remarkably, this new Riemann Solver depends Lipschitz continuously on all parameters.


Introduction
Starting with the seminal papers by Lighthill, Witham, and Richards [12,13], traffic flow on a single road has been modeled in terms of a scalar conservation law: (1.1) Here ρ is the density of cars, while v(ρ) is their velocity, which we assume depends of the density alone. To describe traffic flow on an entire network of roads, one needs to further introduce a set of boundary conditions at road junctions [7]. These conditions should relate the traffic densities on incoming roads i ∈ I and outgoing roads j ∈ O, depending on two main parameters: (i) Driver's turning preferences. For every i, j, one should specify the fraction θ ij ∈ [0, 1] of drivers arriving to from the i-th road, who wish to turn into the j-th road.
(ii) Relative priorities assigned to different incoming roads. If the intersection is congested, these describe the maximum influx of cars arriving from each road i ∈ I, allowed to cross the intersection.
Various junction models of have been proposed in the literature [4,5,7,9]. See also [1] for a survey. A convenient approach is to introduce a Riemann Solver, i.e. a rule that specifies how to construct the solution in the special case where the initial data is constant on every incoming and each outgoing road. As shown in [4], as soon as a Riemann Solver is given, the general Cauchy problem for traffic flow near a junction can be uniquely solved (under suitable assumptions).
The recent counterexamples in [3] show that, on a network of roads, in general the Cauchy problem can be ill posed. Indeed, two distinct solutions can be constructed for the same measurable initial data. On a network with several nodes, nonuniqueness can occur even if the initial data have small total variation. To readdress this situation, in [2] an alternative intersection model was proposed, introducing a buffer of limited capacity at each road junction. For this new model, given any L ∞ initial data, the Cauchy problem has a unique solution, which is robust w.r.t. perturbations of the data. Indeed, one has continuous dependence even w.r.t. the topology of weak convergence.
A natural question, addressed in the present paper, is what happens in the limit as the size of the buffer approaches zero. For Riemann initial data, constant along each incoming and outgoing road, we show that this limit is described by a Limit Riemann Solver (LRS) which can be explicitly determined. See (2.15)-(2.17) in Section 2.
We recall that, in a model without buffer, the initial conditions consist of the constant densities ρ ♦ k on all incoming and outgoing roads k ∈ I ∪O, together with the drivers' turning preferences θ ♦ ij . On the other hand, in the model with buffer, these initial conditions comprise also the length of the queues q ♦ j , j ∈ O, inside the buffer. One can think of q ♦ j as the number of cars already inside the intersection (say, a traffic circle) at time t = 0, waiting to access the outgoing road j. Our main results can be summarized as follows.
(i) For any given Riemann data ρ ♦ k , θ ♦ ij , one can choose initial queue sizes q ♦ j such that, for all t > 0 the solution of the problem with buffer is exactly the same as the solution determined by the Riemann Solver (LRS).
(ii) For any Riemann data ρ ♦ k , θ ♦ ij , and any initial queue sizes q ♦ j , as t → ∞ the solution of the problem with buffer approaches asymptotically the solution determined by the Riemann Solver (LRS).
Using the fact that the conservation laws (1.1) are invariant under space and time rescalings, from (ii) we obtain a convergence result as the size of the buffer approaches zero.
Our present results apply only to solutions of the Riemann problem, i.e. with traffic density which is initially constant along each road. Indeed, for a general Cauchy problem the counterexamples in [3] remain valid also for the Riemann Solver (LRS), showing that the initial-value problem with measurable initial data can be ill posed. Hence no convergence result can be expected. This should not appear as a paradox: for every positive size of the buffer, the Cauchy problem has a unique solution, depending continuously on the initial data. However, as the size of the buffer approaches zero, the solution can become more and more sensitive to small changes in the initial conditions. In the limit, uniqueness is lost.
An extension of our results may be possible in the case of initial data with bounded variation, for a network containing one single node. In view of the results in [4,7], we conjecture that in this case the solution to the Cauchy problem with buffer converges to the solution determined by the Riemann Solver (LRS).

Statement of the main results
Consider a family of n + m roads, joining at a node. Indices i ∈ {1, . . . , m} = I denote incoming roads, while indices j ∈ {m + 1, . . . , m + n} = O denote outgoing roads. On the k-th road, the density of cars ρ k (t, x) is governed by the scalar conservation law Here t ≥ 0, while x ∈ ] − ∞, 0] for incoming roads and x ∈ [0, +∞[ for outgoing roads. The flux function is is the speed of cars on the k-th road. We assume that each flux function f k satisfies where ρ jam k is the maximum possible density of cars on the k-th road. Intuitively, this can be thought as a bumper-to-bumper packing, so that the speed of cars is zero. For a given road k ∈ {1, . . . , m + n}, we denote by  Moreover, we say that Given initial data on each road in order to determine a unique solution to the Cauchy problem we must supplement the conservation laws (2.1) with a suitable set of boundary conditions. These provide additional constraints on the limiting values of the vehicle densities near the intersection. In a realistic model, these boundary conditions should depend on: (i) Relative priority given to incoming roads. For example, if the intersection is regulated by a crosslight, the flow will depend on the fraction η i ∈ ]0, 1[ of time when cars arriving from the i-th road get a green light.
(ii) Drivers' choices. For every i ∈ I, j ∈ O, these are modeled by assigning the fraction θ ij ∈ [0, 1] of drivers arriving from the i-th road who choose to turn into the j-th road. Obvious modeling considerations imply Since we are only interested in the Riemann problem, throughout the following we shall assume that the θ ij are given constants, satisfying (2.6).
In [2] a model of traffic flow at an intersection was introduced, including a buffer of limited capacity. The incoming fluxes of cars toward the intersection are constrained by the current degree of occupancy of the buffer. More precisely, consider an intersection with m incoming and n outgoing roads. The state of the buffer at the intersection is described by an n-vector Here q j (t) is the number of cars at the intersection waiting to enter road j ∈ O (in other words, the length of the queue in front of road j). Boundary values at the junction will be denoted by Conservation of the total number of cars implieṡ at a.e. time t ≥ 0. Here and in the sequel, the upper dot denotes a derivative w.r.t. time. Following [7], we define the maximum possible flux at the end of an incoming road as This is the largest flux f i (ρ) among all states ρ that can be connected toρ i with a wave of negative speed. Notice that the two right hand sides in (2.9) coincide ifρ i = ρ max i .
Similarly, we define the maximum possible flux at the beginning of an outgoing road as As in [2], we assume that the junction contains a buffer of size M . Incoming cars are admitted at a rate depending of the amount of free space left in the buffer, regardless of their destination. Once they are within the intersection, cars flow out at the maximum rate allowed by the outgoing road of their choice.
Single Buffer Junction (SBJ). Consider a constant M > 0, describing the maximum number of cars that can occupy the intersection at any given time, and constants c i > 0, i ∈ I, accounting for priorities given to different incoming roads.
We then require that the incoming fluxesf i satisfȳ In addition, the outgoing fluxesf j should satisfy Here ω i = ω i (ρ i ) and ω j = ω j (ρ j ) are the maximum fluxes defined at (2.9)-(2.10). Notice that (SBJ) prescribes all the boundary fluxesf k , k ∈ I ∪ O, depending on the boundary densities ρ k . It is natural to assume that the constants c i satisfy the inequalities These conditions imply that, when the buffer is empty, cars from all incoming roads can access the intersection with the maximum possible flux (2.9). The analysis in [2] shows that, with the above boundary conditions, the Cauchy problem on a network of roads has a unique solution, continuously depending on the initial data.
The main goal of this paper is to understand what happens when the size of the buffer approaches zero. More precisely, assume that (2.11) is replaced bȳ (2.14) Notice that (2.14) models a buffer with size M ε. When j q j = M ε, the buffer is full and no more cars are admitted to the intersection.
We will show that, as ε → 0, the limit of solutions to the Riemann problem with buffer of vanishing size can be described by a specific Limit Riemann Solver.
be the corresponding maximum possible fluxes at the boundary of the incoming and outgoing roads, as in (2.9)-(2.10). Consider the oneparameter curve Then for t > 0 the Riemann problem is solved by the incoming fluxes In turn, by the conservation of the number of drivers, the outgoing fluxes arē By specifying all the incoming and outgoing fluxesf i ,f j at the intersection, the entire solution of the Riemann problem is uniquely determined. Indeed: (i) For an incoming road i ∈ I, there exists a unique boundary state ρ 0 In this case the density of cars on the i-th road remains constant: , then the solution to the Riemann problem contains only waves with negative speed. In this case, the density of cars on the i-th road coincides with the solution of (2.18), for x < 0.
(ii) For an outgoing road j ∈ O, there exists a unique boundary state ρ 0 In this case the density of cars on the j-th road remains constant: , then the solution to the Riemann problem  Remark 1. For the Riemann Solver constructed in [3], the fluxesf k are locally Hölder continuous functions of the data ρ ♦ k , θ ij , on the domain where θ ij > 0, ω j > 0 for all j ∈ O. The Riemann Solver (LRS) has even better regularity properties. Namely, the fluxesf k defined at (2.15)-(2.17) are locally Lipschitz continuous functions of ρ ♦ k , θ ij . Unfortunately, as remarked earlier, this additional regularity is still not sufficient to guarantee the well-posedness of the Cauchy problem, for general measurable initial data.
Our first result refers to "well prepared" initial data, where the initial lengths of the queues are suitably chosen.

20)
be assigned along each road, together with drivers' turning preferences θ ij .
Then one can choose initial values q ♦ j , j ∈ O for the queues inside the buffer in such a way that the solution to the Riemann problem with buffer coincides with the self-similar solution determined by the Limit Riemann Solver (LRS).
Our second result covers the general case, where the initial sizes of the queues are given arbitrarily, and the solution of the initial value problem with buffer is not self-similar.
Theorem 2. Let the assumptions (2.2), (2.13) hold. Let Riemann data (2.20) be assigned along each road, together with drivers' turning preferences θ ij > 0 and initial queue sizes Then, as t → +∞, the solution (ρ k (t, x)) k∈I∪O to the Riemann problem with buffer asymptotically converges to the self-similar solution (ρ k (t, x)) k∈I∪O determined by the Limit Riemann Solver (LRS). More precisely: A proof of the above theorems will be given in Sections 4 and 5, respectively. By an asymptotic rescaling of time and space, using Theorem 2 we can describe the behavior of the solution to a Riemann problem, as the size of the buffer approaches zero.
Corollary (limit behavior for a buffer of vanishing size). Let f k , θ ij , c i , M be as in Theorem 2. Let Riemann data (2.20) be assigned along each road, together with drivers' turning preferences θ ij > 0 and initial queue sizes as in (2.21).
For ε > 0, let (ρ ε k (t, x)) k∈I∪O be the solution to the initial value problem with a buffer of size M ε, obtained by replacing (2.11) with (2.14) and choosing q ε j (0) = εq ♦ j as initial sizes of the queues.
Callingρ k the self-similar solution determined by the Limit Riemann Solver (LRS) with the same initial data (2.20), for every τ > 0 we have Proof of the Corollary. Let (ρ k (t, x)) k∈I∪O be the solution constructed in Theorem 2. Then, for every i ∈ I, by a rescaling of coordinates we obtain In the last step we used Theorem 2 in connection with the variable change t = τ /ε. For j ∈ O, the difference |ρ ε j −ρ j | is estimated in an entirely similar way.

The Riemann problem with buffer
We consider here an initial value problem with Riemann data, so that the initial density is constant on every incoming and outgoing road.
We decompose the sets of indices as depending on whether these roads are initially free or congested. More precisely: Observe that • If i ∈ I c , then the i-th incoming road will always be congested, i.e. ρ i (t, x) ≥ ρ max i for all t, x.
• If j ∈ O f , then the j-th outgoing road will always be free, i.e. ρ j (t, x) ≤ ρ max j for all t, x.
• If i ∈ I f , then part of the i-th road can become congested (Fig. 3, left).
• If j ∈ O c , then part of the j-th road can become free (Fig. 3, right).
x 0 t 0 x 1 t 2 t 3 Figure 3: Left: an incoming road which is initially free. For t 1 < t < t 2 part of the road is congested (shaded area). Right: an outgoing road which is initially congested. For 0 < t < t 3 part of the road is free (shaded area). In both cases, a shock marks the boundary between the free and the congested region.
The next lemma plays a key role in the proof of Theorem 2. It shows that, for any t > 0, the maximum possible flux at the boundary of any incoming or outgoing road is greater or equal to the maximum flux computed at t = 0.
x), k ∈ I ∪ O be the solution of the Riemann problem with initial data (3.1). As in (2.9)-(2.10) call ω ♦ k = ω k (ρ ♦ k ) the maximum possible fluxes. Similarly, for t > 0 call ω k (t) = ω k (ρ k (t)) the corresponding maximum fluxes. Then Proof. 1. We first consider an incoming road i ∈ I. . For a given t > 0, two subcases may occur.
(i) There exists a characteristic with positive speed, reaching the point (t, 0). Since this characteristic must start at a point x 0 < 0, we conclude that (ii) There exists a neighborhood of (t, 0) covered with characteristics having negative speed.

2.
For an outgoing road j ∈ O, the analysis is similar.
In this case the j-th road always remains free and we have ω j (t) = ω ♦ j = f max j , for every t ≥ 0.

CASE 2:
The road is initially congested, namely ρ ♦ j > ρ max j . For a given t > 0, two subcases may occur.
(i) There exists a characteristic with negative speed, reaching the point (t, 0). Since this characteristic must start at a point x 0 > 0, we conclude that ρ j (t, 0+) = ρ j (0, (ii) There exists a neighborhood of (t, 0) covered with characteristics having positive speed. In this case ρ j (t, 0+) ≥ ρ max j , hence ω j (t) = f max j .

Proof of Theorem 1
Let ρ ♦ k , k ∈ I ∪ O be the initial densities of cars on the incoming and outgoing roads, and let θ ij be the drivers' turning preferences, as in (2.6). Call ω ♦ i , ω ♦ j the maximum possible boundary fluxes on the incoming and outgoing roads, and defines as in (2.16). Two cases will be considered, shown in Fig. 4.
. This is the case where none of the incoming roads remains congested, and all the drivers arriving at the intersection can immediately proceed to the outgoing road of their choice.
In this case we choose the initial queues With these choices, the solution of the Cauchy problem with buffer coincides with the selfsimilar solution determined by the Limit Riemann Solver (LRS). The buffer remains always empty: q j (t) = 0 for all t ≥ 0 and j ∈ O.
CASE 2:s < M . In this case there exists an index j * ∈ O such that When this happens, the entire flow through the intersection is restricted by the number of cars that can exit toward the single congested road j * . We then define and choose the initial queues to be Then the corresponding solution coincides with the self-similar solution determined by the Limit Riemann Solver (LRS). Indeed, by the definition of γ(s), for every j ∈ O we have with equality holding when j = j * . By (4.4), all queues remain constant in time, namely q j * (t) = q * and q j (t) = 0 for j = j * . Remark 2. In the proof of Theorem 1, the queue sizes q ♦ j may not be uniquely determined. Indeed, in Case 2 there may exist two distinct indices j * 1 , j * 2 ∈ O such that When this happens, we can choose the queue sizes to be

Proof of Theorem 2
In this section we prove that, for any initial data, as t → +∞ the solution to the Riemann problem with buffer converges as to the self-similar function determined by the Limit Riemann Solver ( Given the densities ρ ♦ i on the incoming roads i ∈ I, call ω ♦ i the corresponding maximal flows, as in (2.9). Callq i the value of the queue such that Without loss of generality, we can assume 0 ≤q m ≤ · · · ≤q 2 ≤q 1 .
At an intuitive level, we have • If the queue inside the buffer is small, i.e. q <q i , then all drivers arriving from the i-th road can access the intersection, and the i-th road will become free.
• If the queue inside the buffer is large, i.e. q >q i , then not all drivers coming from the i-th road can immediately access the intersection, and the i-th road will become congested.
This can be formulated in a more precise way as follows. By the definition (2.16), if q > M −s one has On the other hand, if q < M −s, let j * ∈ O be an index such that (4.1) holds. We then have The proof is achieved in several steps.

1.
We first study the case where, in the solution determined by the Limit Riemann Solver, at least one of the outgoing roads is congested (Fig. 4, right), so that (4.1) holds. Lets be as in (2.16). As in (4.2), define the asymptotic size of the queue to be q * = M −s > 0. To fix the ideas, assume 0 ≤q m ≤ · · · ≤q ν+1 ≤ q * <q ν ≤ · · · ≤q 2 ≤q 1 .
(i) If t ≥ t ℓ−1 , then we have the implication For all times t ≥ t ℓ the incoming road ℓ is free. Hence its flux near the intersection satisfies The above claim is proved by induction on ℓ = 1, . . . , ν.
Next, for t > τ 1 , the flux of cars arriving to the intersection from road 1 is If road 1 is congested near the intersection, then for t > τ 1 the outgoing flux is for some δ ′ 1 > 0. As a consequence, road 1 must become free within time This proves (iii), in the case ℓ = 1.
The general inductive step is very similar. Assume that the statements (i)-(iii) have been proved for ℓ − 1. Then for t ≥ t ℓ−1 the incoming roads i = 1, . . . , ℓ − 1 are free. The flux of cars reaching the intersection from these roads is f i (t) = ω ♦ i . Now assume that t > t ℓ−1 and q(t) ≥q ℓ . In this case, q(t) ≥q i for all i ∈ I, i ≥ ℓ. Using Lemma 1, for any j ∈ O we thus obtaiṅ Therefore, if q j (t) > 0, thenq j (t) ≤ − 2ε ℓj < 0 for some constants ε ℓj . By continuity, there exists δ ℓ > 0 such that Setting ε ℓ . = min j ε ℓj , we obtain (5.6).
From the implication it follows q(t) ≤q ℓ − δ ℓ for all t ≥ τ ℓ sufficiently large. This yields (ii).
Finally, for t > τ ℓ , the flux of cars arriving to the intersection from road ℓ is If road ℓ is congested near the intersection, then for t > τ ℓ the outgoing flux is for some δ ′ ℓ > 0. As a consequence, road ℓ must become free within time This proves (iii). By induction on ℓ, our claim is proved.

2.
We now prove that, for any ε > 0, there exists a time t ε > t ν large enough so that Indeed, if t > t ℓ , then the same arguments used before yield the implication 3. In this step we prove a lower bound on the queue. Namely, for any ε ♯ > 0 there exists a time t ♯ > 0 such that Toward this goal, choose j * ∈ O such that (4.1) holds. Theṅ Now assume q(t) < q * − ε ♯ . If road j * is initially free, then it remains free for all times t ≥ 0.
In this case we havė for some δ ♯ > 0. Therefore (5.11) holds with Next, we consider the more difficult case where the outgoing road j * is initially congested. We claim that, if q(τ ) ≤ q * at some time τ ≥ t ν , then q(t) ≤ q * for all t ≥ τ . Indeed, for any j ∈ O, if q j (t) > 0, thenq Observing that the right hand side of (5.12) is nonpositive when q(t) = q * , our claim is proved.
Finally, assume that q(τ ) < q * − ε ♯ , for some τ ≥ t ν . Then (5.13) and (5.15) imply (5.11), with 4. Denote by ρ k (t, x), k ∈ I ∪ O, the solution to the Riemann problem with buffer, and σ k (t, x) the self-similar solution determined by the Limit Riemann Solver (LRS). From the limit lim t→∞ q(t) = q * proved in the previous steps, it follows that all boundary fluxes f k (t) converge to the corresponding boundary fluxesf k in the self-similar solution determined by (LRS). Now consider an incoming road i ∈ I. Since the initial data coincide for every t > 0 by [11] we have the estimate For outgoing roads j ∈ O, the estimates are entirely similar. This achieves a proof of Theorem 2 in the case where (4.1) holds for some j * ∈ O.

5.
It remains to consider the case (Fig. 4, left) where i ω ♦ i θ ij < ω ♦ j (5.17) for every j ∈ O. In this case, the arguments in step 1 show that, for all t ≥ t m sufficiently large, all incoming roads become free. In this case, for all times t ≥ t ♯ sufficiently large the incoming fluxes are f i (t) = ω ♦ i =f i . i ∈ I.
Moreover, for t large all queue sizes become q j (t) = 0, and the outgoing fluxes are Inserting these identities in (5.16), we conclude the proof as in the previous case.