# Exact Boundary Controllability for Free Traffic Flow with Lipschitz Continuous State

Authors: Martin Gugat
Journal: Mathematical Problems in Engineering
 (2016)

---

## Abstract

We consider traffic flow governed by the LWR model. We show that a Lipschitz continuous initial density with free-flow and sufficiently small Lipschitz constant can be controlled exactly to an arbitrary constant free-flow density in finite time by a piecewise linear boundary control function that controls the density at the inflow boundary if the outflow boundary is absorbing. Moreover, this can be done in such a way that the generated state is Lipschitz continuous. Since the target states need not be close to the initial state, our result is a global exact controllability result. The Lipschitz constant of the generated state can be made arbitrarily small if the Lipschitz constant of the initial density is sufficiently small and the control time is sufficiently long. This is motivated by the idea that finite or even small Lipschitz constants are desirable in traffic flow since they might help to decrease the speed variation and lead to safer traffic.


---

## Body

## 1. Introduction

The Lighthill Whitham and Richards (LWR) model is a macroscopic model for traffic flow (see [1, 2]). In the model, the traffic flow is described by the solution of an initial boundary value problem with a hyperbolic partial differential equation (pde). The pde can be solved using the method of characteristics. In fact, since the pde is a scalar conservation law, for this model the characteristic curves are straight lines. A derivation of this macroscopic model as a limit of a microscopic model is given in [3].In the control of traffic flow, for example, for the flow through a tunnel, it is desirable to have free-flow traffic in the tunnel where the density remains below a certain critical density that corresponds to maximal throughput in the system. In this way, congested flow in the tunnel is avoided. This can be achieved by controlling the traffic inflow into the tunnel; see [4, 5].It is well known that the solutions of the LWR model can develop shocks in finite time. However, in order to decrease variations in speed among the cars and along the road solutions without shocks are desirable. In this paper we construct Lipschitz continuous controls that generate Lipschitz continuous states. This is motivated by the idea that the traffic safety is increased if the Lipschitz constant is decreased, since this might help to decrease the speed variation among vehicles on the road and as vehicles drive along the road, which in turn reduces the expected accident rate; see [6, 7]. Moreover, we hope that as the Lipschitz constant decreases, also the fuel efficiency of the traffic flow increases. Thus, the proposed controls point out possible improvements for the control of free-flow traffic that may help to improve the safety and fuel efficiency in the free-flow regime. However, in the example of tunnel traffic flow control there is a trade-off between decreasing the Lipschitz constant of the traffic flow density flow within the tunnel and possible longer waiting time of the drivers before they are allowed to enter the tunnel. The reason is that in order to generate Lipschitz continuous states we impose an upper bound on the rate of change of the controlled inflow traffic density. We want to emphasize that the method presented in this paper cannot dissolve existing congested flow within the tunnel.Exact boundary controllability has been studied intensively for classical solution of quasilinear hyperbolic partial differential equations, in particular by Li and his group; see [8]. In this paper we want to consider solutions that are less regular, namely, solutions generated by Lipschitz continuous initial and boundary data. It turns out that the solutions can still be represented using the method of characteristics. Also in [9, 10] the exact controllability has been studied in the framework of classical solutions. Both studies focus on systems where the eigenvalues can change their sign. In [9] the St.-Venant system that is a model for the flow of water through channels is considered, whereas in [10] a general class of systems with vanishing characteristic speed is studied. In the LWR model the characteristic speed vanishes exactly at the critical density where the traffic flow changes from free flow to congested flow. Since our exact controllability is not restricted to desired states that are close to the initial state, we can consider it as a global exact controllability result within the free-flow regime. A result of this kind has also been presented in [11] for the St.-Venant system. In fact, the basic idea of the proof in [11] is similar to the proof presented here. However, in [11] classical solutions are considered.In this paper we show how Lipschitz continuous control functions can be chosen in such a way that they generate Lipschitz continuous states and control the given initial density in finite time to any desired density below the critical density. While the density can be increased arbitrarily fast without generating shock, in order to decrease the density to a substantially smaller density some intermediate steps can be necessary. As an example, a piecewise linear controller is presented explicitly in Section5 of the paper.The paper has the following structure. In Section2, the boundary control system that we want to consider is defined and a well-posedness result in the framework of solutions that are defined in the sense of characteristics is given. Based upon this result, in Section 3, we study the exact controllability of the system. We present a result on global exact boundary controllability to constant states with Lipschitz continuous solutions. To prove this result, we first show two results on local exact controllability. In Section 4 we present a sufficient condition for Lipschitz continuity of the generated states and provide Lipschitz constants explicitly. In Section 5, we summarize the constructed piecewise linear control functions in explicit form and describe the Lipschitz continuous states that they generate.


## 2. The System

First we introduce some notation for the flux function.


### 2.1. The Flux Function

Let a maximal densityρ
m
a
x
>
0 be given. Let a nonzero concave function f that maps [
0
,
ρ
m
a
x
] to 0
,
∞ with f
0
=
0 and f
ρ
m
a
x
=
0 be given.The functionf defines the fundamental diagram (see [12]) and serves as the flux function in our traffic flow model. Define the critical density ρ
c as the smallest point where the function f attains its maximum value on [
0
,
ρ
m
a
x
]. Assume that f is continuously differentiable on [
0
,
ρ
c
]. The concavity of f implies that it is strictly increasing on [
0
,
ρ
c
]. Moreover, f
′ is decreasing on [
0
,
ρ
c
].Assume thatf
′ is Lipschitz continuous on [
0
,
ρ
c
] with the Lipschitz constant L
f
′
>
0; that is, for all ρ
1, ρ
2
∈
[
0
,
ρ
c
] we have the Lipschitz inequality: (1)
f
′
ρ
1
-
f
′
ρ
2
≤
L
f
′
ρ
1
-
ρ
2
.Forε
∈
0
,
ρ
c define (2)
Λ
_
ε
=
min
ρ
∈
0
,
ρ
c
-
ε
⁡
f
′
ρ
=
f
′
ρ
c
-
ε
>
0
.Remark 1.
Since in our analysis we only work in the uncongested phase, in fact only the assumptions onρ as a function defined on the interval [
0
,
ρ
c
] are relevant for our analysis.


### 2.2. The Initial Boundary Value Problem

Letl
>
0 denote the length of the considered road section. In particular, this could be the length of a tunnel. Let a Lipschitz continuous initial density ρ
0
∈
C
0
,
l with values in [
0
,
ρ
c
] and the Lipschitz constant L
ρ
0 be given. Let a time T
>
0 be given. Let a Lipschitz continuous boundary density u
∈
C
0
,
T with the Lipschitz constant L
u be given. Consider the system: (3)
S
ρ
0
,
x
=
ρ
0
x
,
x
∈
0
,
l
ρ
t
,
0
=
u
t
,
t
∈
0
,
T
ρ
t
+
f
ρ
x
=
0
,
on
0
,
T
×
0
,
l
.Theorem 3 contains sufficient conditions for the existence of a solution of S in the sense of the characteristic curves on the time interval [
0
,
T
]. To avoid the generation of shocks at t
=
0, it is assumed that u
0
<
ρ
c and that u and ρ
0 are compatible in the sense that u
0
=
ρ
0
0. The sufficient conditions are stated in terms of the size of the Lipschitz constants for the initial density and the boundary density. Note that in Theorem 3 we do not control the growth of the Lipschitz constants during the evolution of the solution. The reason is that we only assume that the time T is strictly smaller than the breaking time, which can be estimated quite accurately. The corresponding result for the Burgers equation has been stated in [13]. The size of the Lipschitz constants of the generated states is studied in Section 4. In Theorem 12 explicit representations for these Lipschitz constants are presented that are increasing functions of the Lipschitz constants of the initial density and the boundary density and also of the time T. To make sure that the solution exists on a given possibly large time interval [
0
,
T
], the Lipschitz constants L
ρ
0 and L
u must be sufficiently small. Note that the values of the boundary density and hence also the initial density at x
=
0 are assumed to be strictly below the critical density ρ
c.Remark 2.
Foru
t
∈
0
,
ρ
c, the boundary condition ρ
t
,
0
=
u
(
t
) is equivalent to f
ρ
t
,
0
=
f
u
t. Thus, the boundary condition can be implemented by controlling the influx of vehicles per second as in ramp metering (see, e.g., [14–16]).For a given upper boundρ
c
-
ε for the values of the boundary density and given values of the Lipschitz constants L
ρ
0 and L
u, Theorem 3 guarantees the existence of the solution without shocks only on a possibly short time interval [
0
,
T
].In the definition ofS and thus also in Theorem 3 the boundary conditions at x
=
l are not considered. If congested flow occurs at x
=
l, a shock wave might travel backwards into the region 0
,
T
×
0
,
l, thus interfering with the regular solution that is constructed in Theorem 3. To avoid this effect, throughout the paper we assume free outflow at x
=
l on the whole time interval 0
,
T; that is, we assume that x
=
l is an absorbing boundary, where we do not impose boundary conditions.Theorem 3 (quasilinear initial boundary value problem).
LetT
>
0 and ε
∈
0
,
ρ
c be given. Assume that the initial density ρ
0 is Lipschitz continuous on [
0
,
l
] with the Lipschitz constant L
ρ
0
>
0 and that for all x
∈
[
0
,
l
] one has 0
≤
ρ
0
x
≤
ρ
c.
Assume that the boundary densityu is Lipschitz continuous on [
0
,
T
] with the Lipschitz constant L
u
>
0 and that for all t
∈
[
0
,
T
] one has 0
≤
u
t
≤
ρ
c
-
ε.
Assume that theC
0-compatibility conditions between ρ
0 and u hold: that is, (4)
ρ
0
0
=
u
0
≤
ρ
c
-
ε
.Assume that (5)
T
<
min
⁡
1
L
f
′
L
ρ
0
,
Λ
_
ε
L
f
′
L
u
.
Then systemS has a solution on [
0
,
T
]
×
[
0
,
l
] in the sense of characteristics without shocks and rarefaction fans. In particular, the solution is continuous.Proof of Theorem3.
We consider a solution ofS in the sense of characteristics. We define the characteristic curves that intersect the x-axis at the time t
=
0; we have for x
∈
[
0
,
l
] and s
≥
0
(6)
ξ
s
,
x
,
0
=
x
+
s
f
′
ρ
0
x
.For the characteristic curves that intersect the axis x
=
0 at a time t
∈
[
0
,
T
] we have for s
≥
0
(7)
ξ
s
,
0
,
t
=
s
-
t
f
′
u
t
.Then for each point s
,
z
∈
[
0
,
T
]
×
[
0
,
l
] there exists a characteristic curve such that ξ
s
,
x
,
0
=
z or such that ξ
s
,
0
,
t
=
z.
In order to obtain a well-defined solution in the sense of characteristics it is sufficient to show thatT is chosen in such a way that it is impossible that the characteristic curves intersect for s
≤
T.
Then at each points
,
z
∈
[
0
,
T
]
×
[
0
,
l
] the solution is well defined since the value of ρ is constant along the characteristic curves.
There are three possibilities of how two different characteristic curves can intersect at a points
,
z
∈
[
0
,
T
]
×
[
0
,
l
]:(1)
Two characteristic curves of type (6) intersect: the equation ξ
s
,
x
1
,
0
=
ξ
s
,
x
2
,
0 with x
1, x
2
∈
0
,
l and x
1
≠
x
2 implies f
′
ρ
0
x
1
≠
f
′
ρ
0
x
2 and (8)
s
=
-
1
f
′
ρ
0
x
2
-
f
′
ρ
0
x
1
/
x
2
-
x
1
.
Due to our regularity assumptions onf
′ and ρ
0, the function x
↦
f
′
ρ
0
x is Lipschitz continuous with the Lipschitz constant L
f
′
L
ρ
0. Hence, (8) implies the inequality s
≥
1
/
L
f
′
L
ρ
0. Hence, if T
<
1
/
L
f
′
L
ρ
0, an intersection of this type cannot occur for s
∈
[
0
,
T
].
(2)
Two characteristic curves of type (7) intersect: the equation ξ
s
,
0
,
t
1
=
ξ
s
,
0
,
t
2 with t
1, t
2
∈
[
0
,
T
] and t
1
≠
t
2 implies (9)
s
-
t
2
=
f
′
u
t
1
f
′
u
t
2
-
f
′
u
t
1
/
t
2
-
t
1
.
Due to our regularity assumptions onf
′ and u, the function t
↦
f
′
u
t is Lipschitz continuous with the Lipschitz constant L
f
′
L
u. Thus, using the definition of Λ
_
ε in (2) we arrive at the inequality (10)
s
-
t
2
≥
Λ
_
ε
L
f
′
L
u
.
Thus, ifT
<
Λ
_
ε
/
L
f
′
L
u, an intersection of this type is impossible for s
∈
0
,
T.
(3)
Characteristic curves of type (6) intersect characteristic curves of type (7): for t
∈
0
,
T and x
∈
0
,
l the equation ξ
s
,
0
,
t
=
ξ
s
,
x
,
0 implies (11)
s
=
1
f
′
u
t
-
f
′
ρ
0
x
/
t
f
′
u
t
+
x
.
Due to our compatibility assumption (4) we have (12)
f
′
u
t
-
f
′
ρ
0
x
=
f
′
u
t
-
f
′
u
0
+
f
′
ρ
0
0
-
f
′
ρ
0
x
≤
L
f
′
L
u
t
+
L
f
′
L
ρ
0
x
≤
L
f
′
L
u
t
f
′
u
t
Λ
_
ε
+
L
f
′
L
ρ
0
x
≤
m
a
x
L
f
′
L
u
Λ
_
ε
,
L
f
′
L
ρ
0
t
f
′
u
t
+
x
.
Thus, we obtain the inequality(13)
s
≥
1
max
⁡
L
f
′
L
u
/
Λ
_
ε
,
L
f
′
L
ρ
0
.
IfT satisfies (5) an intersection of this kind cannot occur for s
∈
[
0
,
T
].
Thus, we have proved Theorem3.Remark 4.
Ifρ
0 satisfies the assumptions of Theorem 3, the compatibility assumption (4) holds, ρ
0
x
≤
ρ
0
0 for all x
∈
0
,
l, and the control function u is continuous and increasing on [
0
,
T
] with values in (
0
,
ρ
c
), then the solution exists as a continuous solution in the sense of characteristics if T
<
1
/
L
f
′
L
ρ
0.
The reason is that in this case for intersections of the characteristic curves of the second type, ift
1
<
t
2 we have s
-
t
2
<
0; hence, ξ
s
,
0
,
t
2
<
0. Thus, there is no such intersection of the second type with s
∈
0
,
T and ξ
s
,
0
,
t
2
∈
0
,
l.
Moreover, sinceu is increasing, u
t
>
ρ
0
x for all t
∈
0
,
T, x
∈
0
,
l. Hence, for intersections of the characteristic curves of the third type, we have f
′
u
t
<
f
′
ρ
0
x; thus, in this case we have s
<
0, so intersections of this type do not occur with s
∈
0
,
T.Example 5.
Letρ
m
a
x
>
0 and u
max
⁡
>
0 and for ρ
∈
[
0
,
ρ
m
a
x
] define (14)
f
ρ
=
u
max
⁡
ρ
1
-
ρ
ρ
m
a
x(see [12] and the introduction in [17]). Then ρ
c
=
(
1
/
2
)
ρ
m
a
x and L
f
′
=
2
(
u
max
⁡
/
ρ
m
a
x
).
Letl
=
2 and for x
∈
[
0
,
l
] define the initial density: (15)
ρ
0
x
=
1
16
1
+
cos
⁡
π
x
+
π
+
4
-
π
20
x
.Figure 1 shows the graph of the initial density on the space interval [
0,2
].
Chooseρ
m
a
x
≥
m
a
x
x
∈
[
0
,
l
]
ρ
0
(
x
). We have L
ρ
0
=
π
/
16
+
(
4
-
π
)
/
20, ρ
0
(
0
)
=
0, and f
′
ρ
0
0
=
u
max
⁡. Hence, (16)
1
L
f
′
L
ρ
0
=
ρ
m
a
x
u
max
⁡
π
/
8
+
4
-
π
/
10
.
Forβ
∈
0
,
ρ
c choose (17)
L
u
≤
u
max
⁡
1
-
2
β
ρ
m
a
x
π
16
+
4
-
π
20
.
Fort
≥
0 define the control: (18)
u
t
=
L
u
t
if
t
∈
0
,
β
L
u
,
β
if
t
>
β
L
u
.Then the compatibility condition (4) holds and u is Lipschitz continuous with the Lipschitz constant L
u.
Chooseε
=
ρ
c
-
β. Then we have Λ
_
(
ε
)
=
f
′
(
β
) and hence due to (17) (19)
1
L
f
′
L
ρ
0
≤
Λ
_
ε
L
f
′
L
u
.Theorem 3 implies that, for all (20)
T
<
ρ
m
a
x
u
max
⁡
π
/
8
+
4
-
π
/
10
,(5) holds. In particular, the solution in the sense of characteristics exists without shocks and rarefaction fans on the time interval [
0,2
ρ
m
a
x
/
u
max
⁡
]. If we choose L
u as the right-hand side in (17), we have (21)
β
L
u
<
1
L
f
′
L
ρ
0if and only if (22)
β
<
1
4
ρ
m
a
x
.Hence, if (22) holds the solution exists at least until u
t has attained the value β. Figure 2 shows the characteristic curves corresponding to the generated state ρ
(
t
,
x
) on [
0
,
T
]
×
[
0
,
l
]
=
[
0,2
]
×
[
0,2
] for ρ
m
a
x
=
1, u
m
a
x
=
1, β
=
0.2, and L
u chosen as the right-hand side in (17). The horizontal axis is the space axis.Figure 1
The initial densityρ
0 for Example 5 on the space interval [
0,2
].Figure 2
The characteristic curves for Example5 generated with initial density ρ
0, ρ
m
a
x
=
1, and u
m
a
x
=
1 for β
=
0.2 on the time interval [
0,2
]. Since at time 2 the characteristic curves that contain the information of the initial density have already left the space interval [
0,2
], this implies that with the constant control u
t
=
β for t
>
2 the solution continues to exist without shocks.Example 6.
Choosef and ρ
0 as in Example 5. For ω
>
0 define (23)
u
t
=
L
u
ω
sin
⁡
ω
t(see Figure 3). Then the compatibility condition (4) holds and u is Lipschitz continuous with the Lipschitz constant L
u. Choose ε
=
ρ
c
-
L
u
/
ω. Then we have Λ
_
(
ε
)
=
f
′
(
L
u
/
ω
). If (24)
L
u
≤
u
m
a
x
L
ρ
0
1
+
2
u
m
a
x
L
ρ
0
/
ω
ρ
m
a
x
,inequality (19) holds. Hence, Theorem 3 implies that, for all T
>
0 that satisfy (20), a solution without shocks exists on [
0
,
T
]. In particular, such a solution exists on the time interval [
0,2
(
ρ
m
a
x
/
u
max
⁡
)
].
Figure4 shows the characteristic curves corresponding to the state for ρ
m
a
x
=
1, u
m
a
x
=
1, ω
=
6, and L
u chosen as the right-hand side in (24) on the time interval [
0,2
]. Again the horizontal axis is the space axis [
0
,
l
]. The example illustrates that also oscillating inflow can generate a Lipschitz continuous solution if the Lipschitz constant of the boundary density is sufficiently small.Figure 3
The functionu for Example 6 for ω
=
6 and L
u chosen as the right-hand side in (24) on the time interval [
0,2
]. The horizontal axis is the time axis.Figure 4
The characteristic curves for Example6 generated with the initial density ρ
0 from Example 5, ρ
m
a
x
=
1, u
m
a
x
=
1, ω
=
6, and L
u chosen as the right-hand side in (24) on the time interval [
0,2
].


## 2.1. The Flux Function

Let a maximal densityρ
m
a
x
>
0 be given. Let a nonzero concave function f that maps [
0
,
ρ
m
a
x
] to 0
,
∞ with f
0
=
0 and f
ρ
m
a
x
=
0 be given.The functionf defines the fundamental diagram (see [12]) and serves as the flux function in our traffic flow model. Define the critical density ρ
c as the smallest point where the function f attains its maximum value on [
0
,
ρ
m
a
x
]. Assume that f is continuously differentiable on [
0
,
ρ
c
]. The concavity of f implies that it is strictly increasing on [
0
,
ρ
c
]. Moreover, f
′ is decreasing on [
0
,
ρ
c
].Assume thatf
′ is Lipschitz continuous on [
0
,
ρ
c
] with the Lipschitz constant L
f
′
>
0; that is, for all ρ
1, ρ
2
∈
[
0
,
ρ
c
] we have the Lipschitz inequality: (1)
f
′
ρ
1
-
f
′
ρ
2
≤
L
f
′
ρ
1
-
ρ
2
.Forε
∈
0
,
ρ
c define (2)
Λ
_
ε
=
min
ρ
∈
0
,
ρ
c
-
ε
⁡
f
′
ρ
=
f
′
ρ
c
-
ε
>
0
.Remark 1.
Since in our analysis we only work in the uncongested phase, in fact only the assumptions onρ as a function defined on the interval [
0
,
ρ
c
] are relevant for our analysis.


## 2.2. The Initial Boundary Value Problem

Letl
>
0 denote the length of the considered road section. In particular, this could be the length of a tunnel. Let a Lipschitz continuous initial density ρ
0
∈
C
0
,
l with values in [
0
,
ρ
c
] and the Lipschitz constant L
ρ
0 be given. Let a time T
>
0 be given. Let a Lipschitz continuous boundary density u
∈
C
0
,
T with the Lipschitz constant L
u be given. Consider the system: (3)
S
ρ
0
,
x
=
ρ
0
x
,
x
∈
0
,
l
ρ
t
,
0
=
u
t
,
t
∈
0
,
T
ρ
t
+
f
ρ
x
=
0
,
on
0
,
T
×
0
,
l
.Theorem 3 contains sufficient conditions for the existence of a solution of S in the sense of the characteristic curves on the time interval [
0
,
T
]. To avoid the generation of shocks at t
=
0, it is assumed that u
0
<
ρ
c and that u and ρ
0 are compatible in the sense that u
0
=
ρ
0
0. The sufficient conditions are stated in terms of the size of the Lipschitz constants for the initial density and the boundary density. Note that in Theorem 3 we do not control the growth of the Lipschitz constants during the evolution of the solution. The reason is that we only assume that the time T is strictly smaller than the breaking time, which can be estimated quite accurately. The corresponding result for the Burgers equation has been stated in [13]. The size of the Lipschitz constants of the generated states is studied in Section 4. In Theorem 12 explicit representations for these Lipschitz constants are presented that are increasing functions of the Lipschitz constants of the initial density and the boundary density and also of the time T. To make sure that the solution exists on a given possibly large time interval [
0
,
T
], the Lipschitz constants L
ρ
0 and L
u must be sufficiently small. Note that the values of the boundary density and hence also the initial density at x
=
0 are assumed to be strictly below the critical density ρ
c.Remark 2.
Foru
t
∈
0
,
ρ
c, the boundary condition ρ
t
,
0
=
u
(
t
) is equivalent to f
ρ
t
,
0
=
f
u
t. Thus, the boundary condition can be implemented by controlling the influx of vehicles per second as in ramp metering (see, e.g., [14–16]).For a given upper boundρ
c
-
ε for the values of the boundary density and given values of the Lipschitz constants L
ρ
0 and L
u, Theorem 3 guarantees the existence of the solution without shocks only on a possibly short time interval [
0
,
T
].In the definition ofS and thus also in Theorem 3 the boundary conditions at x
=
l are not considered. If congested flow occurs at x
=
l, a shock wave might travel backwards into the region 0
,
T
×
0
,
l, thus interfering with the regular solution that is constructed in Theorem 3. To avoid this effect, throughout the paper we assume free outflow at x
=
l on the whole time interval 0
,
T; that is, we assume that x
=
l is an absorbing boundary, where we do not impose boundary conditions.Theorem 3 (quasilinear initial boundary value problem).
LetT
>
0 and ε
∈
0
,
ρ
c be given. Assume that the initial density ρ
0 is Lipschitz continuous on [
0
,
l
] with the Lipschitz constant L
ρ
0
>
0 and that for all x
∈
[
0
,
l
] one has 0
≤
ρ
0
x
≤
ρ
c.
Assume that the boundary densityu is Lipschitz continuous on [
0
,
T
] with the Lipschitz constant L
u
>
0 and that for all t
∈
[
0
,
T
] one has 0
≤
u
t
≤
ρ
c
-
ε.
Assume that theC
0-compatibility conditions between ρ
0 and u hold: that is, (4)
ρ
0
0
=
u
0
≤
ρ
c
-
ε
.Assume that (5)
T
<
min
⁡
1
L
f
′
L
ρ
0
,
Λ
_
ε
L
f
′
L
u
.
Then systemS has a solution on [
0
,
T
]
×
[
0
,
l
] in the sense of characteristics without shocks and rarefaction fans. In particular, the solution is continuous.Proof of Theorem3.
We consider a solution ofS in the sense of characteristics. We define the characteristic curves that intersect the x-axis at the time t
=
0; we have for x
∈
[
0
,
l
] and s
≥
0
(6)
ξ
s
,
x
,
0
=
x
+
s
f
′
ρ
0
x
.For the characteristic curves that intersect the axis x
=
0 at a time t
∈
[
0
,
T
] we have for s
≥
0
(7)
ξ
s
,
0
,
t
=
s
-
t
f
′
u
t
.Then for each point s
,
z
∈
[
0
,
T
]
×
[
0
,
l
] there exists a characteristic curve such that ξ
s
,
x
,
0
=
z or such that ξ
s
,
0
,
t
=
z.
In order to obtain a well-defined solution in the sense of characteristics it is sufficient to show thatT is chosen in such a way that it is impossible that the characteristic curves intersect for s
≤
T.
Then at each points
,
z
∈
[
0
,
T
]
×
[
0
,
l
] the solution is well defined since the value of ρ is constant along the characteristic curves.
There are three possibilities of how two different characteristic curves can intersect at a points
,
z
∈
[
0
,
T
]
×
[
0
,
l
]:(1)
Two characteristic curves of type (6) intersect: the equation ξ
s
,
x
1
,
0
=
ξ
s
,
x
2
,
0 with x
1, x
2
∈
0
,
l and x
1
≠
x
2 implies f
′
ρ
0
x
1
≠
f
′
ρ
0
x
2 and (8)
s
=
-
1
f
′
ρ
0
x
2
-
f
′
ρ
0
x
1
/
x
2
-
x
1
.
Due to our regularity assumptions onf
′ and ρ
0, the function x
↦
f
′
ρ
0
x is Lipschitz continuous with the Lipschitz constant L
f
′
L
ρ
0. Hence, (8) implies the inequality s
≥
1
/
L
f
′
L
ρ
0. Hence, if T
<
1
/
L
f
′
L
ρ
0, an intersection of this type cannot occur for s
∈
[
0
,
T
].
(2)
Two characteristic curves of type (7) intersect: the equation ξ
s
,
0
,
t
1
=
ξ
s
,
0
,
t
2 with t
1, t
2
∈
[
0
,
T
] and t
1
≠
t
2 implies (9)
s
-
t
2
=
f
′
u
t
1
f
′
u
t
2
-
f
′
u
t
1
/
t
2
-
t
1
.
Due to our regularity assumptions onf
′ and u, the function t
↦
f
′
u
t is Lipschitz continuous with the Lipschitz constant L
f
′
L
u. Thus, using the definition of Λ
_
ε in (2) we arrive at the inequality (10)
s
-
t
2
≥
Λ
_
ε
L
f
′
L
u
.
Thus, ifT
<
Λ
_
ε
/
L
f
′
L
u, an intersection of this type is impossible for s
∈
0
,
T.
(3)
Characteristic curves of type (6) intersect characteristic curves of type (7): for t
∈
0
,
T and x
∈
0
,
l the equation ξ
s
,
0
,
t
=
ξ
s
,
x
,
0 implies (11)
s
=
1
f
′
u
t
-
f
′
ρ
0
x
/
t
f
′
u
t
+
x
.
Due to our compatibility assumption (4) we have (12)
f
′
u
t
-
f
′
ρ
0
x
=
f
′
u
t
-
f
′
u
0
+
f
′
ρ
0
0
-
f
′
ρ
0
x
≤
L
f
′
L
u
t
+
L
f
′
L
ρ
0
x
≤
L
f
′
L
u
t
f
′
u
t
Λ
_
ε
+
L
f
′
L
ρ
0
x
≤
m
a
x
L
f
′
L
u
Λ
_
ε
,
L
f
′
L
ρ
0
t
f
′
u
t
+
x
.
Thus, we obtain the inequality(13)
s
≥
1
max
⁡
L
f
′
L
u
/
Λ
_
ε
,
L
f
′
L
ρ
0
.
IfT satisfies (5) an intersection of this kind cannot occur for s
∈
[
0
,
T
].
Thus, we have proved Theorem3.Remark 4.
Ifρ
0 satisfies the assumptions of Theorem 3, the compatibility assumption (4) holds, ρ
0
x
≤
ρ
0
0 for all x
∈
0
,
l, and the control function u is continuous and increasing on [
0
,
T
] with values in (
0
,
ρ
c
), then the solution exists as a continuous solution in the sense of characteristics if T
<
1
/
L
f
′
L
ρ
0.
The reason is that in this case for intersections of the characteristic curves of the second type, ift
1
<
t
2 we have s
-
t
2
<
0; hence, ξ
s
,
0
,
t
2
<
0. Thus, there is no such intersection of the second type with s
∈
0
,
T and ξ
s
,
0
,
t
2
∈
0
,
l.
Moreover, sinceu is increasing, u
t
>
ρ
0
x for all t
∈
0
,
T, x
∈
0
,
l. Hence, for intersections of the characteristic curves of the third type, we have f
′
u
t
<
f
′
ρ
0
x; thus, in this case we have s
<
0, so intersections of this type do not occur with s
∈
0
,
T.Example 5.
Letρ
m
a
x
>
0 and u
max
⁡
>
0 and for ρ
∈
[
0
,
ρ
m
a
x
] define (14)
f
ρ
=
u
max
⁡
ρ
1
-
ρ
ρ
m
a
x(see [12] and the introduction in [17]). Then ρ
c
=
(
1
/
2
)
ρ
m
a
x and L
f
′
=
2
(
u
max
⁡
/
ρ
m
a
x
).
Letl
=
2 and for x
∈
[
0
,
l
] define the initial density: (15)
ρ
0
x
=
1
16
1
+
cos
⁡
π
x
+
π
+
4
-
π
20
x
.Figure 1 shows the graph of the initial density on the space interval [
0,2
].
Chooseρ
m
a
x
≥
m
a
x
x
∈
[
0
,
l
]
ρ
0
(
x
). We have L
ρ
0
=
π
/
16
+
(
4
-
π
)
/
20, ρ
0
(
0
)
=
0, and f
′
ρ
0
0
=
u
max
⁡. Hence, (16)
1
L
f
′
L
ρ
0
=
ρ
m
a
x
u
max
⁡
π
/
8
+
4
-
π
/
10
.
Forβ
∈
0
,
ρ
c choose (17)
L
u
≤
u
max
⁡
1
-
2
β
ρ
m
a
x
π
16
+
4
-
π
20
.
Fort
≥
0 define the control: (18)
u
t
=
L
u
t
if
t
∈
0
,
β
L
u
,
β
if
t
>
β
L
u
.Then the compatibility condition (4) holds and u is Lipschitz continuous with the Lipschitz constant L
u.
Chooseε
=
ρ
c
-
β. Then we have Λ
_
(
ε
)
=
f
′
(
β
) and hence due to (17) (19)
1
L
f
′
L
ρ
0
≤
Λ
_
ε
L
f
′
L
u
.Theorem 3 implies that, for all (20)
T
<
ρ
m
a
x
u
max
⁡
π
/
8
+
4
-
π
/
10
,(5) holds. In particular, the solution in the sense of characteristics exists without shocks and rarefaction fans on the time interval [
0,2
ρ
m
a
x
/
u
max
⁡
]. If we choose L
u as the right-hand side in (17), we have (21)
β
L
u
<
1
L
f
′
L
ρ
0if and only if (22)
β
<
1
4
ρ
m
a
x
.Hence, if (22) holds the solution exists at least until u
t has attained the value β. Figure 2 shows the characteristic curves corresponding to the generated state ρ
(
t
,
x
) on [
0
,
T
]
×
[
0
,
l
]
=
[
0,2
]
×
[
0,2
] for ρ
m
a
x
=
1, u
m
a
x
=
1, β
=
0.2, and L
u chosen as the right-hand side in (17). The horizontal axis is the space axis.Figure 1
The initial densityρ
0 for Example 5 on the space interval [
0,2
].Figure 2
The characteristic curves for Example5 generated with initial density ρ
0, ρ
m
a
x
=
1, and u
m
a
x
=
1 for β
=
0.2 on the time interval [
0,2
]. Since at time 2 the characteristic curves that contain the information of the initial density have already left the space interval [
0,2
], this implies that with the constant control u
t
=
β for t
>
2 the solution continues to exist without shocks.Example 6.
Choosef and ρ
0 as in Example 5. For ω
>
0 define (23)
u
t
=
L
u
ω
sin
⁡
ω
t(see Figure 3). Then the compatibility condition (4) holds and u is Lipschitz continuous with the Lipschitz constant L
u. Choose ε
=
ρ
c
-
L
u
/
ω. Then we have Λ
_
(
ε
)
=
f
′
(
L
u
/
ω
). If (24)
L
u
≤
u
m
a
x
L
ρ
0
1
+
2
u
m
a
x
L
ρ
0
/
ω
ρ
m
a
x
,inequality (19) holds. Hence, Theorem 3 implies that, for all T
>
0 that satisfy (20), a solution without shocks exists on [
0
,
T
]. In particular, such a solution exists on the time interval [
0,2
(
ρ
m
a
x
/
u
max
⁡
)
].
Figure4 shows the characteristic curves corresponding to the state for ρ
m
a
x
=
1, u
m
a
x
=
1, ω
=
6, and L
u chosen as the right-hand side in (24) on the time interval [
0,2
]. Again the horizontal axis is the space axis [
0
,
l
]. The example illustrates that also oscillating inflow can generate a Lipschitz continuous solution if the Lipschitz constant of the boundary density is sufficiently small.Figure 3
The functionu for Example 6 for ω
=
6 and L
u chosen as the right-hand side in (24) on the time interval [
0,2
]. The horizontal axis is the time axis.Figure 4
The characteristic curves for Example6 generated with the initial density ρ
0 from Example 5, ρ
m
a
x
=
1, u
m
a
x
=
1, ω
=
6, and L
u chosen as the right-hand side in (24) on the time interval [
0,2
].


## 3. Exact Controllability to Constant Densities

Based on the solutions ofS that we have obtained in Theorem 3, we can show that starting from a given initial density ρ
0 with values in [
0
,
ρ
c
] we can control the traffic flow to any constant density b
<
ρ
c in finite time. As before, this requires the assumption that at x
=
l free traffic outflow occurs.Theorem 7 (exact controllability to constant states).
Assume that the initial densityρ
0 is Lipschitz continuous on [
0
,
l
] with the Lipschitz constant L
ρ
0 and that for all x
∈
[
0
,
l
] one has 0
≤
ρ
0
x
≤
ρ
c and ρ
0
0
<
ρ
c. Moreover, assume that (25)
L
ρ
0
<
f
′
ρ
0
0
l
L
f
′
.
Then for allβ
∈
0
,
ρ
c there exist a time T
>
0 and a piecewise linear Lipschitz continuous boundary control function u
t such that system S has a continuous solution on 0
,
T
×
0
,
l without shocks in the sense of characteristics which satisfies for all x
∈
0
,
l the end condition: (26)
ρ
T
,
x
=
β
.Thus, the traffic flow can be controlled arbitrarily close to a state of maximal throughput with a shock-free state.To prove Theorem7, we show first that any initial state with a sufficiently small Lipschitz constant and ρ
0
0
<
ρ
c can be controlled to the constant state ρ
(
x
)
=
ρ
0
(
0
) in finite time.Theorem 8 (exact controllability toρ
0
0).
Assume that the initial densityρ
0 is Lipschitz continuous on 0
,
l with the Lipschitz constant L
ρ
0 and that for all x
∈
0
,
l one has 0
≤
ρ
0
x
≤
ρ
c and ρ
0
(
0
)
<
ρ
c. Moreover, assume that (27)
L
ρ
0
<
f
′
ρ
0
0
l
L
f
′
.Let T
=
l
/
f
′
ρ
0
0. Then the constant control u
(
t
)
=
ρ
0
(
0
) generates a continuous solution on 0
,
T
×
0
,
l without shocks in the sense of characteristics for system (
S
) that satisfies for all x
∈
0
,
l the end condition: (28)
ρ
T
,
x
=
ρ
0
0
.Thus, the traffic flow can be controlled to a constant state without the generation of shocks.Proof of Theorem8.
Letε
=
ρ
c
-
ρ
0
0
>
0. Define the control: (29)
u
t
=
ρ
0
0
.Then u is Lipschitz continuous for any Lipschitz constant L
u
>
0. Moreover, it satisfies the compatibility assumption (4) and 0
≤
u
(
t
)
≤
ρ
c
-
ε. Hence, Theorem 3 implies that for all T
>
0 that satisfy (5), system (
S
) has a continuous solution on [
0
,
T
]
×
[
0
,
l
] without shocks in the sense of characteristics
The definition ofu implies that the slope of the characteristic curve that goes through x
=
0 at t
=
0 is f
′
(
ρ
0
(
0
)
). Hence, the travel time until a characteristic curve that transports the values of the boundary density u
(
t
) reaches x
=
l is (30)
τ
=
l
f
′
ρ
0
0
.Hence, if there exists a number T
>
0 that satisfies (5) and (31)
τ
≤
T
,end condition (26) holds. Due to (30), this is equivalent to the inequality (32)
l
f
′
ρ
0
0
<
min
⁡
1
L
f
′
L
ρ
0
,
f
′
ρ
0
0
L
f
′
L
u
.Due to assumption (27) we have (33)
l
f
′
ρ
0
0
<
1
L
f
′
L
ρ
0
.Moreover, we can choose L
u
>
0 such that (34)
L
u
<
f
′
ρ
0
0
2
L
f
′
.Then (32) holds. Thus, we have proved Theorem 8.Theorem8 states that we can steer initial densities with sufficiently small Lipschitz constants to a constant state. Starting from this constant state, we can steer the system to any other constant state in 0
,
ρ
c that is sufficiently close to the initial state in finite time.Theorem 9 (local exact controllability between constant states).
Assume that the initial density is constant; that is,ρ
0
(
x
)
=
α on [
0
,
l
] with α
∈
[
0
,
ρ
c
). Let a number β
∈
[
0
,
ρ
c
) be given. Assume that (35)
β
-
α
<
min
⁡
f
′
α
,
f
′
β
L
f
′
.
Then there exist a timeT
>
0 and piecewise linear Lipschitz continuous control u
(
t
) that generates a continuous solution on [
0
,
T
]
×
[
0
,
l
] without shocks in the sense of characteristics for system (
S
) that satisfies for all x
∈
[
0
,
l
] the end condition: (36)
ρ
T
,
x
=
β
.
Thus, the traffic flow can be controlled between two free-flow constant states without the generation of shocks.Proof of Theorem9.
Letε
=
min
⁡
ρ
c
-
α
,
ρ
c
-
β
>
0. For L
u
>
0 define the control: (37)
u
t
=
α
+
sign
⁡
β
-
α
L
u
t
if
t
∈
0
,
β
-
α
L
u
,
β
if
t
>
β
-
α
L
u
.Then u is Lipschitz continuous with the Lipschitz constant L
u. Moreover, it satisfies the compatibility assumption (4) and 0
≤
u
(
t
)
≤
ρ
c
-
ε. Hence, Theorem 3 implies that (
S
) has a continuous solution on [
0
,
T
]
×
[
0
,
l
] without shocks in the sense of characteristics if T satisfies (5). Moreover, the minimal slope of the characteristic curves that go through x
=
0 is Λ
_
(
ε
); hence, the maximal travel time until a characteristic curve that transports the values of the boundary density u
(
t
) reaches x
=
l is (38)
τ
m
a
x
=
l
Λ
_
ε
.The definition of u implies that, after the time β
-
α
/
L
u, the travel time until a characteristic curve that transports the values of the boundary density u
(
t
)
=
β reaches x
=
l is l
/
f
′
(
β
)
. Hence, if (39)
β
-
α
L
u
+
l
f
′
β
≤
T
,end condition (36) holds. Thus, it remains to find values of β such that there exists a time T that satisfies (39) and (5) at the same time. This is the case if (40)
β
-
α
L
u
+
l
f
′
β
<
min
⁡
1
L
f
′
L
ρ
0
,
Λ
_
ε
L
f
′
L
u
.Since we can choose L
ρ
0 arbitrarily small, we only have to show that (41)
β
-
α
L
u
+
l
f
′
β
<
Λ
_
ε
L
f
′
L
u
.Assumption (35) implies that (42)
β
-
α
<
Λ
_
ε
L
f
′
.Thus, we can choose a number L
u
>
0 such that (43)
L
u
<
f
′
β
l
Λ
_
ε
L
f
′
-
β
-
α
.Then (41) holds. Thus, we have proved Theorem 9.Example 10.
In this example we study a decreasing control function for a constant initial state. Letρ
0
(
x
)
=
α
∈
(
0
,
ρ
c
) and β
∈
[
0
,
α
). Choose L
u
>
0.
Fort
≥
0 define the control: (44)
u
t
=
α
-
L
u
t
if
t
∈
0
,
α
-
β
L
u
,
β
if
t
>
α
-
β
L
u
.Then compatibility condition (4) holds and u is Lipschitz continuous with the Lipschitz constant L
u. Theorem 9 implies that if (45)
α
-
β
<
f
′
α
L
f
′
,we can choose L
u sufficiently small such that a solution in the sense of characteristics is generated. More precisely this is the case if L
u satisfies (43): that is, (46)
L
u
<
f
′
β
l
f
′
α
L
f
′
-
α
-
βand at the time t
∗
=
(
α
-
β
)
/
L
u
+
l
/
f
′
(
β
), the system has reached the state ρ
(
x
,
t
∗
)
=
β.
Figure5 shows the characteristic curves corresponding to the state for f as in Example 5: ρ
m
a
x
=
1, u
m
a
x
=
1, α
=
0.4, β
=
0.35, and L
u
=
0.007 on the time interval [
0,14
]. Again the horizontal axis is the space axis [
0
,
l
].Figure 5
The characteristic curves for Example10 generated with the initial density ρ
0
=
0.4 and β
=
0.35, L
u
=
0.007, ρ
m
a
x
=
1, and u
m
a
x
=
1.Example 11.
In this example we study an increasing control function for a constant initial state. Letρ
0
(
x
)
=
α
∈
(
0
,
ρ
c
) and β
∈
(
α
,
ρ
c
). Choose L
u
>
0. For t
≥
0 define the control: (47)
u
t
=
α
+
L
u
t
if
t
∈
0
,
β
-
α
L
u
,
β
if
t
>
β
-
α
L
u
.Then compatibility condition (4) holds and u is Lipschitz continuous with the Lipschitz constant L
u. As pointed out in Remark 4, the control generates a solution without shocks and rarefaction fans for arbitrarily large T.
At the timet
∗
=
(
β
-
α
)
/
L
u
+
l
/
f
′
(
β
), the system has reached the state ρ
(
x
,
t
∗
)
=
β.
The generated solution is Lipschitz continuous. In fact for all(
t
1
,
x
), (
t
2
,
x
)
∈
[
0
,
T
]
×
[
0
,
l
], we have the inequality (48)
ρ
t
1
,
x
-
ρ
t
2
,
x
≤
L
u
t
1
-
t
2
.For all (
t
,
ξ
1
), (
t
,
ξ
2
)
∈
[
0
,
T
]
×
[
0
,
l
], we have the inequality (49)
ρ
t
,
ξ
1
-
ρ
t
,
ξ
2
≤
L
u
1
f
′
u
T
ξ
1
-
ξ
2
.This can be seen as follows. Let ξ
1
=
ξ
(
t
,
0
,
τ
1
) and ξ
2
=
ξ
(
t
,
0
,
τ
2
) with 0
≤
τ
1
<
τ
2
≤
t. Then we have ξ
2
<
ξ
1 and (50)
τ
2
-
τ
1
=
τ
1
-
t
f
′
u
τ
1
-
f
′
u
τ
2
+
ξ
1
-
ξ
2
f
′
u
τ
2
.Since u is increasing, we have f
′
(
u
(
τ
1
)
)
-
f
′
(
u
(
τ
2
)
)
≥
0; thus, (50) implies (51)
0
≤
τ
2
-
τ
1
≤
1
f
′
u
τ
2
ξ
1
-
ξ
2
.Thus, (52)
ρ
t
,
ξ
2
-
ρ
t
,
ξ
1
=
u
τ
2
-
u
τ
1
≤
L
u
τ
2
-
τ
1
≤
L
u
f
′
u
T
ξ
1
-
ξ
2
.If one or both characteristic curves through (
t
,
ξ
1
), (
t
,
ξ
2
) start at τ
=
0 with x
>
0, where the initial density is prescribed, (49) also follows since ρ
0 is constant.
Due to (48) and (49), for all (
t
1
,
ξ
1
), (
t
2
,
ξ
2
)
∈
[
0
,
T
]
×
[
0
,
l
], we have the inequality (53)
ρ
t
1
,
ξ
1
-
ρ
t
2
,
ξ
2
≤
L
u
t
1
-
t
2
+
1
f
′
u
T
ξ
1
-
ξ
2
.Hence, the generated state is Lipschitz continuous and the size of the Lipschitz constant on [
0
,
T
]
×
[
0
,
l
] is proportional to L
u.
Figure6 shows the characteristic curves corresponding to the state for f as in Example 5: α
=
0, ρ
m
a
x
=
1, u
m
a
x
=
1, L
u
=
6, and β
=
0.4 on the time interval [
0,10
]. Again the horizontal axis is the space axis [
0
,
l
].Figure 6
The characteristic curves for Example11 generated with the initial density ρ
0
=
0 and β
=
0.4, L
u
=
6, ρ
m
a
x
=
1, and u
m
a
x
=
1.Using Theorems8 and 9, we show that it is possible to steer initial densities with sufficiently small Lipschitz constants to any free-flow constant state.Proof of Theorem7.
Starting with the constant control from Theorem8, we can control the initial state in finite time to a constant state a
=
ρ
0
(
0
). On account of Theorem 9, it suffices that there exist a finite number of states z
0
,
z
1
,
…
,
z
N
∈
[
0
,
m
a
x
{
a
,
β
}
] with z
0
=
a and z
N
=
β such that for all i
∈
{
1
,
…
,
N
} it is possible to steer the system from z
i
-
1 to z
i in finite time with a piecewise linear Lipschitz continuous control. Inequality (35) in Theorem 9 implies that this is possible if for all i
∈
{
1
,
…
,
N
} we have the inequality (54)
z
i
-
z
i
-
1
≤
min
⁡
f
′
z
i
,
f
′
z
i
-
1
L
f
′
.Since for all j
∈
{
0
,
…
,
N
} we have (55)
f
′
z
j
≥
min
⁡
f
′
a
,
f
′
β
, (54) holds if z
i are chosen such that for all i
∈
{
1
,
…
,
N
} we have (56)
z
i
-
z
i
-
1
≤
min
⁡
f
′
a
,
f
′
β
L
f
′
.This is possible if N is sufficiently large, for example, if the points z
i are distributed equidistantly. Then we can reach the desired system state b after N intermediate control steps of the type described in Theorem 9. Thus, we have proved Theorem 7.


## 4. Lipschitz Constants for the State

Up to now, we have shown that, for a givenε
>
0 that is sufficiently small, a continuous state is generated with suitable Lipschitz initial and boundary data. In this section we provide upper bounds for the Lipschitz constants of the generated states as functions of ε.Theorem 12 (Lipschitz continuous states).
LetT
>
0 and ε
∈
(
0
,
ρ
c
) be given. Assume that the conditions stated in Theorem 3 hold.
Define the numbers(57)
L
a
=
L
u
Λ
_
ε
-
T
L
f
′
L
u
,
L
b
=
L
ρ
0
1
-
T
L
f
′
L
ρ
0
.
Then the solution of(
S
) on [
0
,
T
]
×
[
0
,
l
] is Lipschitz continuous in the sense that for all (
t
1
,
ξ
1
), (
t
2
,
ξ
2
)
∈
[
0
,
T
]
×
0
,
l we have the inequality (58)
ρ
t
1
,
ξ
1
-
ρ
t
2
,
ξ
2
≤
L
a
+
L
b
ξ
1
-
ξ
2
+
Λ
_
ε
t
1
-
t
2
.Proof of Theorem12.
Define the sets(59)
A
a
=
s
,
x
∈
0
,
T
×
0
,
l
:
x
≥
ξ
s
,
0,0
,
A
b
=
s
,
x
∈
0
,
T
×
0
,
l
:
x
≤
ξ
s
,
0,0
.Then A
a
∪
A
b
=
[
0
,
T
]
×
[
0
,
l
] and S
=
A
a
∩
A
b
=
{
(
s
,
x
)
∈
[
0
,
T
]
×
[
0
,
l
]
:
x
=
ξ
(
s
,
0,0
)
}. We show that the solution ρ of (
S
) is Lipschitz continuous on A
a with the Lipschitz constant L
a and ρ is Lipschitz continuous on A
b with the Lipschitz constant L
b in the sense that for j
∈
{
a
,
b
} and for all (
t
1
,
ξ
1
), (
t
2
,
ξ
2
)
∈
A
j we have the inequality (60)
ρ
t
1
,
ξ
1
-
ρ
t
2
,
ξ
2
≤
L
j
ξ
1
-
ξ
2
+
Λ
_
ε
t
1
-
t
2
.
Let(
t
1
,
ξ
1
), (
t
2
,
ξ
2
)
∈
A
a be given. Then there exist points x
1, x
2
∈
[
0
,
l
] such that ξ
1
=
ξ
(
t
1
,
x
1
,
0
) and ξ
2
=
ξ
(
t
2
,
x
2
,
0
). We have (61)
ξ
1
-
ξ
2
=
x
1
-
x
2
+
t
1
f
′
ρ
0
x
1
-
f
′
ρ
0
x
2
+
t
1
-
t
2
f
′
ρ
0
x
2
.Hence, we have (62)
ξ
1
-
ξ
2
≥
x
1
-
x
2
-
T
L
f
′
L
ρ
0
x
1
-
x
2
-
Λ
_
ε
t
1
-
t
2
.Thus, we have (63)
x
1
-
x
2
≤
1
1
-
T
L
f
′
L
ρ
0
ξ
1
-
ξ
2
+
Λ
_
ε
t
1
-
t
2
.This yields (64)
ρ
t
1
,
ξ
1
-
ρ
t
2
,
ξ
2
=
ρ
0
x
1
-
ρ
x
2
≤
L
ρ
0
x
1
-
x
2
≤
L
ρ
0
1
-
T
L
f
′
L
ρ
0
ξ
1
-
ξ
2
+
Λ
_
ε
t
1
-
t
2
.Thus, we have shown (60) for j
=
a.
Let(
t
1
,
ξ
1
), (
t
2
,
ξ
2
)
∈
A
b be given. Then there exist points τ
1, τ
2
∈
[
0
,
T
] such that ξ
1
=
ξ
(
t
1
,
0
,
τ
1
) and ξ
2
=
ξ
(
t
2
,
0
,
τ
2
). We have (65)
ξ
1
-
ξ
2
=
t
1
-
t
2
f
′
u
τ
1
+
t
2
-
τ
1
f
′
u
τ
1
-
f
′
u
τ
2
+
τ
2
-
τ
1
f
′
u
τ
2
.Hence, we have (66)
ξ
1
-
ξ
2
≥
Λ
_
ε
-
t
2
-
τ
1
L
f
′
L
u
τ
1
-
τ
2
-
Λ
_
ε
t
1
-
t
2
.Thus, we obtain (67)
τ
1
-
τ
2
≤
1
Λ
_
ε
-
T
L
f
′
L
u
ξ
1
-
ξ
2
+
Λ
_
ε
t
1
-
t
2
.This yields (68)
ρ
t
1
,
ξ
1
-
ρ
t
2
,
ξ
2
=
u
τ
1
-
u
τ
2
≤
L
u
τ
1
-
τ
2
≤
L
u
Λ
_
ε
-
T
L
f
′
L
u
ξ
1
-
ξ
2
+
Λ
_
ε
t
1
-
t
2
.Thus, we have shown (60) for j
=
b.
For points(
t
1
,
ξ
1
)
∈
A
a, (
t
2
,
ξ
2
)
∈
A
b define the point (
t
~
,
ξ
~
)
∈
S as the point where the line that connects (
t
1
,
ξ
1
) and (
t
2
,
ξ
2
) intersects S. Then we have (69)
ρ
t
1
,
ξ
1
-
ρ
t
~
,
ξ
~
≤
L
a
ξ
1
-
ξ
~
+
Λ
_
ε
t
1
-
t
~
,
ρ
t
2
,
ξ
2
-
ρ
t
~
,
ξ
~
≤
L
b
ξ
2
-
ξ
~
+
Λ
_
ε
t
2
-
t
~
.This yields (70)
ρ
t
1
,
ξ
1
-
ρ
t
2
,
x
2
≤
ρ
t
1
,
ξ
1
-
ρ
t
~
,
ξ
~
+
ρ
t
2
,
ξ
2
-
ρ
t
~
,
ξ
~
≤
L
a
+
L
b
ξ
1
-
x
2
+
Λ
_
ε
t
1
-
t
2
.Thus, we have shown (58) and Theorem 12 is proved.Remark 13.
Theorem12 shows that the size of the Lipschitz constant L
a
+
L
b of the state ρ can be made arbitrarily small by the choice of ε, T, L
u, and L
ρ
0.
IfL
ρ
0 is sufficiently small, L
b
>
0 can be made arbitrarily small.
Moreover, ifL
u is sufficiently small, L
a
>
0 can be made arbitrarily small.
Hence, in the situation of Theorem8, if L
ρ
0 is sufficiently small the Lipschitz constant of the generated state ρ can be made arbitrarily small.
In the situation of Theorem9, L
b
>
0 can be chosen arbitrarily small. Due to (43) we can choose L
u arbitrarily small. Hence, we can make L
a arbitrarily small. Thus, we can control the system from a constant state α to a constant state β in such a way that the Lipschitz constant L
a
+
L
b
>
0 of the generated state remains arbitrarily small.
Using the method described in Theorem7, this implies that if L
ρ
0 is sufficiently small, it is possible to control the system from the initial state ρ
0 with ρ
0
(
0
)
<
ρ
c to a constant state β
∈
(
0
,
ρ
c
) in such a way that the generated state is Lipschitz continuous and that the Lipschitz constant L
a
+
L
b
>
0 of the generated state remains arbitrarily small.


## 5. The Lipschitz Controller

For the convenience of the reader, we present the constructed Lipschitz continuous controller in explicit form. Let an initial densityρ
0 that satisfies the assumptions of Theorem 7 be given. Let a desired density β
∈
[
0
,
ρ
c
) be given. Define a
=
ρ
0
(
0
).We consider three cases.(1)
Ifβ
=
ρ
0
(
0
), we can directly apply Theorem 8 where the constant control u
(
t
)
=
ρ
0
(
0
) is used.
(2)
The second case isβ
∈
(
ρ
0
(
0
)
,
ρ
c
). Define the density z
0
=
a and the time t
~
0
=
l
/
f
′
(
a
) and let (71)
t
^
1
=
t
~
0
+
β
-
a
L
u
,
t
~
1
=
t
^
1
+
l
f
′
β
.
Similarly as in Example11, for L
u
>
0 define the increasing control: (72)
u
t
=
a
if
t
∈
0
,
t
~
0
,
a
+
L
u
t
-
t
~
0
if
t
∈
t
~
0
,
t
^
1
,
β
if
t
∈
t
^
1
,
t
~
1
.
Thenu is Lipschitz continuous with the Lipschitz constant L
u. Theorem 8 implies that, at the time t
~
0, the system is in the constant state a. As pointed out in Example 11, at the time t
~
1 the system has reached the constant state β. Moreover, since the control is constant on the time interval [
0
,
t
~
0
], Theorem 12 implies that the generated state is Lipschitz continuous on [
0
,
t
~
0
] with the Lipschitz constant: (73)
L
∗
,
0
=
L
b
=
L
ρ
0
1
-
t
~
0
L
f
′
L
ρ
0
.
Note that due to (25), we have L
∗
,
0
>
0. In Example 11 we have stated that on the time interval [
t
~
0
,
t
~
1
] the generated state is Lipschitz continuous in the sense of (53). In particular, the Lipschitz constant of the state can be made arbitrarily small by choosing L
u arbitrarily small.
(3)
Now we come to the third case whereβ
∈
[
0
,
ρ
0
(
0
)
). Here if the distance between the initial state and the target state is not sufficiently small, we have to steer the state through a certain number of intermediate states. In this way we construct a decreasing piecewise linear control that is constant on the intermediate states.
For this purpose, choosez
0
,
z
1
,
z
2
,
…
,
z
N such that a
=
z
0
>
z
1
>
⋯
>
z
N
=
β and a number L
u
>
0 such that for all i
∈
{
1
,
…
,
N
} we have (74)
L
u
<
f
′
z
i
l
f
′
z
i
-
1
L
f
′
-
z
i
-
1
-
z
i
.
This can always be achieved with intermediate states on an equidistant grid, where for alli
∈
{
1
,
…
,
N
} we have z
i
=
z
0
+
i
(
(
β
-
a
)
/
N
) if N is sufficiently large. However, in order to reduce the number of intermediate states often it makes sense to increase the step size (
z
i
-
1
-
z
i
) as i increases, since f
′
(
z
i
-
1
) also increases with i.
Define the timet
~
0
=
l
/
f
′
(
a
) and for all i
∈
{
1
,
…
,
N
} let (75)
t
^
i
=
t
~
i
-
1
+
z
i
-
1
-
z
i
L
u
,
t
~
i
=
t
^
i
+
l
f
′
z
i
.
ForL
u
>
0 define the control: (76)
u
t
=
a
if
t
∈
0
,
t
~
0
,
z
i
-
1
-
L
u
t
-
t
~
i
-
1
if
t
∈
t
~
i
-
1
,
t
^
i
,
i
∈
1
,
…
,
N
,
z
i
if
t
∈
t
^
i
,
t
~
i
,
i
∈
1
,
…
,
N
.
Thenu is Lipschitz continuous with the Lipschitz constant L
u. Theorem 8 implies that, at the time t
~
0, the system has reached the constant state z
0. For all i
∈
{
1,2
,
…
,
N
}, for the time interval [
t
~
i
-
1
,
t
~
i
] due to (74) the initial state z
i
-
1 and the terminal state z
i satisfy (35). Moreover, condition (43) holds. Hence, the proof of Theorem 9 shows that, for all i
∈
{
1,2
,
…
,
N
}, at the times t
~
i the system has reached the constant state z
i. Thus, at the terminal time (77)
T
=
t
~
N
=
a
-
β
L
u
+
∑
j
=
0
N
l
f
′
z
i
end condition (26) holds for all x
∈
[
0
,
l
]. Moreover, since the control is constant on the time interval [
0
,
t
~
0
], Theorem 12 implies that the generated state is Lipschitz continuous on [
0
,
t
~
0
] with the Lipschitz constant: (78)
L
∗
,
0
=
L
b
=
L
ρ
0
1
-
t
~
0
L
f
′
L
ρ
0
.
Note that due to (25), we have L
∗
,
0
>
0. For all i
∈
{
1,2
,
…
,
N
}, define T
~
i
=
t
~
i
-
t
~
i
-
1
. For the time interval [
t
~
i
-
1
,
t
~
i
] that has the length T
~
i, the initial state is constant. Due to the construction in the proof of Theorem 9, T
~
i satisfies (5). Hence, Theorem 12 implies that the generated state is Lipschitz continuous on [
t
~
i
-
1
,
t
~
i
] in the sense of (58) with the Lipschitz constant: (79)
L
∗
,
i
=
L
u
f
′
z
i
-
1
-
T
~
i
L
f
′
L
u
>
0
.


## 6. Conclusions

We have shown an exact controllability result for traffic flow described by the LWR model for free flow. We have shown that initial densities with sufficiently small Lipschitz constants can be steered to an arbitrary constant free-flow density in such a way that the Lipschitz constant of the generated state remains arbitrarily small. We have given the corresponding controllers explicitly. The density can be increased arbitrarily fast and still generate a Lipschitz continuous solution. However, in order to decrease the density and still generate a Lipschitz continuous solution some intermediate constant states are necessary.The extension of this result to the case of road networks, for example, to study ramp metering, is a subject of future studies. Our results indicate that it makes sense to include upper bounds on the Lipschitz constant of the control functions as constraints in the corresponding optimal control problems. For the case of states where shocks are admitted, these problems have been studied, for example, in [14–16].Recently, there has been a renewed interest in second-order models for traffic flow; see [17, 18]. In the two-phase models considered in [18], the uncongested phase is modelled by the LWR model, so this case is still covered by our results. It would be interesting to study whether the exact controllability results can also be adapted to the one-phase second-order models.



---

*Source: 2743251-2016-04-12.xml*