# The Burgers Equation for a New Continuum Model with Consideration of Driver’s Forecast Effect

Authors: Lei Yu; Bingchang Zhou
Journal: Journal of Applied Mathematics
 (2014)

---

## Abstract

A new continuum model with consideration of driver’s forecast effect is obtained to study the density wave problem and the stop-and-go phenomena. The stability condition of the new model is derived by using linear analysis. The triangular shock wave, one type of density wave, which is determined by Burgers equation in the stable region, is discussed in great detail with reductive perturbation method. The local cluster appears when we perform the numerical simulations for the new model. It also proves that the driver’s forecast effect has the positive effect of reducing the local cluster.


---

## Body

## 1. Introduction

Traffic jams, the typical signature of the complex behavior of vehicular traffic, have been studied by various traffic models [1–8]. From different theoretical basis, there are microscopic and macroscopic models to describe traffic flow. The dynamical aspects of microscopic models are based on the description of the individual vehicles’ situation. This process is largely determined by the drivers’ behavior and the physical performance of vehicles. The macroscopic models describe traffic streams as a compressible fluid obeying global rules. This coarse-grained means is needed for understanding the collective behavior of traffic, designing efficient control strategies, developing macroscopic traffic simulation, and so forth.Bando et al. [9] propose the optimal velocity (OV) model to characterize the car-following behavior. Although the OV model is shown to have the universal structure in describing many properties of traffic flow, many approaches to extending the model toward more realistic traffic model have been pursued. Helbing and Tilch [10] develop a generalized force model with a velocity difference term added into the OV model. Xue et al. [11] extend the OV model to take into account the effect of the relative velocity.The pioneer work of continuum traffic flow models is the LWR model [12, 13]. Although the LWR model can reproduce most basic traffic flow phenomena such as traffic congestion formation and dissipation in heavy traffic, this model can not describe nonequilibrium traffic flow dynamics and does not have the ability to explain the amplification of small disturbances in heavy traffic. To overcome the deficiencies in the LWR model, various macroscopic traffic models have been proposed. Payne [14] introduces a high-order continuum model which can describe the amplification of small disturbances in heavy traffic and allow fluctuations of speed around the equilibrium value. Thus, the Payne model is suitable to describe nonequilibrium situations such as stop-and-go traffic. However, a fundamental principle of traffic flow, that is, a car being influenced only by the motion of cars ahead of it, not by the motion of cars behind it, is violated in the Payne model because one characteristic speed of this model is greater than the macroscopic flow speed. Notice that cars are anisotropic particles and respond only to frontal stimuli. Jiang et al. [15] develop a macroscopic continuum model based on a car-following theory. This model overcomes the characteristic speed problem that exists in many high-order continuum models such as the Payne model.However, the consistency between microscopic and macroscopic models has been proved [16]. There are also some papers to derive macroscopic continuum models from micromodels (mainly the car-following models), such as [14, 15, 17].There are many studies to reveal nonlinear phenomena of vehicular traffic, such as stop-and-go, phase transition, self-organized, and the nonlinear waves. Kurtze and Hong [18] derive the Korteweg-de Vries (KdV) equation from one continuum traffic flow model. Zhou et al. [19] obtain the KdV equation and the modified Korteweg-de Vries (mKdV) equation from the continuum traffic flow model derived from a car-following model. Ou [17] obtains the Burgers equation and the KdV equation from a continuum version of the full velocity difference car-following model. Yu et al. [20] investigate density waves in an optimal velocity model with reaction-time delay of drivers and derive the Burgers, KdV, and mKdV equations. However, the nonlinear waves results of the continuum traffic flow models are far less than those of the car-following models. The reason is probably the complex partial differential forms of the continuum traffic flow models.In fact, the future traffic situation can be forecasted by the intelligent transportation system (ITS) based on the current traffic status, so the driver may be guided by the forecast information to adjust his/her current acceleration. However, few models consider the drivers forecast effect. Recently, a new car-following model with the driver’s forecast effect is proposed by Tang et al. [21]. Similar to the other models which consider the information of ITS, the model presented by Tang et al. [21] can improve the stability of traffic flow and reduce traffic jams.In this paper, a new macroversion is obtained based on the anisotropic continuum model proposed by Tang et al. [21]. The density wave problem and the stop-and-go phenomena are studied. In Section 2, the stability condition of the model is derived. In Section 3, the triangular shock wave, which is determined by the Burgers equation in the stable region, is discussed in great detail by using nonlinear analysis. In Section 4, the simulation results are given. These results prove that the driver’s forecast effect has the positive effect of reducing the local cluster. Finally, a summary is given.


## 2. The Model and Its Stability Analysis

The new car-following model with the consideration of the driver’s forecast effect can be written as follows [21]:
(1)
d
v
n
(
t
)
d
t
=
κ
(
V
(
Δ
x
n
(
t
)
)
-
v
n
(
t
)
)
+
β
κ
(
V
(
Δ
x
n
(
t
+
τ
)
)
-
v
n
(
t
+
τ
)
)
,
where κ is the reactive coefficient, β is the coefficient of the driver’s forecast effect, and τ is the time-step of the driver forecast. Using the transformation between microvariables and macrovariables, the corresponding anisotropic macrocontinuum model of (1) is obtained [21], that is, the following equations:
(2)
∂
ρ
∂
t
+
∂
(
ρ
v
)
∂
x
=
0
,
∂
v
∂
t
+
v
∂
v
∂
x
=
(
1
+
β
)
(
V
e
(
ρ
)
-
v
)
T
+
β
τ
+
β
τ
u
e
′
(
h
)
c
0
∂
v
∂
x
,
where T is the reactive time given by the inverse of the coefficient κ, c
0
=
ϵ
/
(
β
τ
+
T
)
>
0 is the propagating velocity of the small perturbation, ϵ is the distance between the following and leading vehicles in micromodel, and
(3)
u
e
(
h
)
=
V
e
(
ρ
)
,
h
=
1
/
ρ
.To derive the stability condition, the Burgers equation, and its shock solution of the macroversion of (1), the following new macromodel is derived from (2) according to (3); that is,
(4)
∂
ρ
∂
t
+
∂
(
ρ
v
)
∂
x
=
0
,
∂
v
∂
t
+
v
∂
v
∂
x
=
γ
(
V
e
(
ρ
)
-
v
)
-
β
τ
c
0
ρ
2
V
e
′
(
ρ
)
∂
v
∂
x
,
where γ
=
(
1
+
β
)
/
(
T
+
β
τ
). Obviously, model (4) is anisotropy.The linear stability theory will be applied to derive the linear stability condition of model (4). Assume traffic to be initially in a state differing infinitesimally from the uniform steady flow. Similar to the decomposition of the flow of (4) into a linear combination of Fourier modes in [18], we have
(5)
ρ
(
x
,
t
)
=
ρ
0
+
∑
k
‍
ρ
^
k
exp
⁡
(
i
k
x
+
σ
k
t
)
,
v
(
x
,
t
)
=
v
0
+
∑
k
‍
v
^
k
exp
⁡
(
i
k
x
+
σ
k
t
)
,
where ρ
0 and v
0 are the uniform steady states of (4). Substituting (5) into (4), linearizing and neglecting the higher-order terms of the small perturbations ρ
^
k and v
^
k, we have
(6)
(
σ
k
+
i
k
v
0
)
ρ
^
k
+
i
k
ρ
0
v
^
k
=
0
,
(7)
σ
k
v
^
k
+
i
k
v
0
v
^
k
=
γ
(
V
e
′
(
ρ
0
)
ρ
^
k
-
v
^
k
)
-
β
τ
c
0
ρ
0
2
V
e
′
(
ρ
0
)
i
k
v
^
k
.
From (6), we have
(8)
ρ
^
k
=
-
i
k
ρ
0
v
^
k
σ
k
+
i
k
v
0
.
Substituting (8) into (7), the following quadratic equation is obtained:
(9)
(
σ
k
+
i
k
v
0
)
2
+
(
γ
+
β
τ
c
0
ρ
0
2
V
e
′
(
ρ
0
)
i
k
)
(
σ
k
+
i
k
v
0
)
+
γ
ρ
0
V
e
′
(
ρ
0
)
i
k
=
0
.
Consider the long wave expansion of σ
k in (9), which is determined order by order around i
k
≈
0 [8]. By expanding σ
k
=
σ
1
(
i
k
)
+
σ
2
(
i
k
)
2
+
⋯ and separating the real part and the imaginary part, the coefficients of i
k and (
i
k
)
2 are derived as follows:
(10)
σ
1
=
-
v
0
-
ρ
0
V
e
′
(
ρ
0
)
,
σ
2
=
1
γ
ρ
0
2
V
e
′
2
(
ρ
0
)
(
ϖ
ρ
0
-
1
)
,
where ϖ
=
β
τ
c
0.The linear stability condition of (4) is decided by the real part of σ
k. In fact, the uniform traffic flow described by (4) is stable against all infinitesimal disturbances when σ
2
>
0; that is,
(11)
ϖ
ρ
0
>
1
.
It means that the traffic flow described by (4) is stable when (11) holds.


## 3. Nonlinear Analysis

To investigate the system behavior in the stable region, we consider the long wavelength modes on coarse-grained scales. The reductive perturbation method as in [8] is applied to (4). Introduce the slow scales for space variable x and time variable t and define slow variables X and T as follows:
(12)
X
=
ɛ
(
x
-
b
t
)
,
T
=
ε
2
t
,
where b is a constant to be determined and 0
<
ɛ
≪
1. We set the density and velocity as
(13)
ρ
(
x
,
t
)
=
ρ
0
+
ɛ
ρ
^
(
X
,
T
)
,
v
(
x
,
t
)
=
v
0
+
ɛ
v
^
(
X
,
T
)
.
Substituting (12) and (13) into (4) and making the Taylor expansions to ε
3, the following nonlinear partial differential equations are derived:
(14)
ε
2
(
ρ
0
∂
X
v
^
+
v
0
∂
X
ρ
^
-
b
∂
X
ρ
^
)
+
ε
3
(
∂
T
ρ
^
+
ρ
^
∂
X
v
^
+
v
^
∂
X
ρ
^
)
=
0
,
(15)
ɛ
γ
(
v
^
-
V
e
′
(
ρ
0
)
ρ
^
)
+
ε
2
[
ρ
^
v
^
+
(
v
0
-
b
+
ϖ
ρ
0
2
V
e
′
(
ρ
0
)
)
∂
X
v
^
1
2
-
V
e
′
(
ρ
0
)
ρ
^
2
-
1
2
γ
V
e
′′
(
ρ
0
)
ρ
^
2
]
+
ε
3
[
∂
T
v
^
+
v
^
∂
X
v
^
-
1
6
γ
V
e
′′′
(
ρ
0
)
ρ
^
3
1
6
+
(
ϖ
ρ
0
2
V
e
′′
(
ρ
0
)
+
2
ϖ
ρ
0
V
e
′
(
ρ
0
)
)
ρ
^
∂
X
v
^
]
=
0
.
From the coefficient of ɛ in (15), we have
(16)
v
^
=
V
e
′
(
ρ
0
)
ρ
^
+
O
(
ε
2
)
.
The relationship between the perturbation of density and velocity is given by (16) which is the basis of Burgers equation’s derivation. According to (16), we have the value of b from the second term of ɛ in (14); that is,
(17)
b
=
ρ
0
V
e
′
(
ρ
0
)
+
v
0
.
Then we have
(18)
q
^
=
ρ
^
v
^
=
(
V
e
′
(
ρ
0
)
+
1
2
γ
V
e
′
′
(
ρ
0
)
)
ρ
^
2
+
(
1
-
ϖ
ρ
0
)
ρ
0
V
e
′
2
(
ρ
0
)
∂
X
ρ
^
.
Combining (14) and (18), the following equation is derived:
(19)
∂
T
ρ
^
+
(
2
V
e
′
(
ρ
0
)
+
γ
V
e
′
′
(
ρ
0
)
)
ρ
^
∂
X
ρ
^
=
(
ϖ
ρ
0
-
1
)
ρ
0
V
e
′
2
(
ρ
0
)
∂
X
2
ρ
^
.
In accordance with stability condition (11), the coefficient of the second derivative term on the right hand side of (19) is positive in the stable region. Thus, in the stable region, (19) is just the Burgers equation.Ifρ
^
(
X
,
0
) is of compact support, the solution ρ
^
(
X
,
T
) of Burgers equation, (19), behaves like a train of N-triangular shock waves; that is,
(20)
ρ
^
(
X
,
T
)
=
1
|
c
1
|
T
(
X
-
η
n
+
η
n
+
1
2
)
-
η
n
+
1
-
η
n
2
|
c
1
|
T
*
tanh
(
c
2
4
|
c
1
|
T
(
η
n
+
1
-
η
n
)
(
X
-
ξ
n
)
)
,
where c
1
=
2
V
e
′
(
ρ
0
)
+
γ
V
e
′
′
(
ρ
0
), c
2
=
(
ϖ
ρ
0
-
1
)
ρ
0
V
e
′
2
(
ρ
0
), ξ
n are the coordinates of the shock fronts, and η
n are the coordinates of the intersections of the slopes with the x-axis (n
=
1,2
,
…
,
N). ρ
^
(
X
,
T
) decays to 0 like O
(
1
/
T
) when T
→
+
∞. That is to say, any shock wave expressed by (19) in stable traffic flow region will evolve to a uniform flow in the course of time. We see this phenomenon in Figures 1(a) and 1(e) of Section 4.Temporal evolution of traffic flow before 40 minutes for differentρ
0 with β
=
0.2.
(a)
ρ
0
=
0.03
veh
/
m
(b)
ρ
0
=
0.042
veh
/
m
(c)
ρ
0
=
0.05
veh
/
m
(d)
ρ
0
=
0.06
veh
/
m
(e)
ρ
0
=
0.08
veh
/
m


## 4. The Numerical Simulation

To check the theoretical results, we carry out numerical simulations for model (4) by using the numerical scheme in [21]. The difference equations are as follows:
(21)
ρ
i
j
+
1
=
ρ
i
j
+
Δ
t
Δ
x
ρ
i
j
(
v
i
j
-
v
i
+
1
j
)
+
Δ
t
Δ
x
v
i
j
(
ρ
i
-
1
j
-
ρ
i
j
)
,(a) ifv
i
j
<
-
ϖ
(
ρ
i
j
)
2
V
e
′
(
ρ
i
j
),
(22)
v
i
j
+
1
=
v
i
j
+
Δ
t
Δ
x
(
-
ϖ
(
ρ
i
j
)
2
V
e
′
(
ρ
i
j
)
-
v
i
j
)
(
v
i
+
1
j
-
v
i
j
)
+
γ
Δ
t
(
V
e
(
ρ
i
j
)
-
v
i
j
)
,(b) ifv
i
j
≥
-
ϖ
(
ρ
i
j
)
2
V
e
′
(
ρ
i
j
),
(23)
v
i
j
+
1
=
v
i
j
+
Δ
t
Δ
x
(
-
ϖ
(
ρ
i
j
)
2
V
e
′
(
ρ
i
j
)
-
v
i
j
)
(
v
i
j
-
v
i
-
1
j
)
+
γ
Δ
t
(
V
e
(
ρ
i
j
)
-
v
i
j
)
,
where i, j, Δ
x, and Δ
t represent the road section, time, spatial step, and time step, respectively, ρ
i
j
≈
ρ
(
x
i
,
t
j
), v
i
j
≈
v
(
x
i
,
t
j
), x
i
=
i
Δ
x, t
j
=
j
Δ
t, and i
,
j are integers.To study the local cluster effect of (4), we perform numerical simulations over a system of 32.2 km long highway using (21)–(23). The local cluster effect corresponds to the stop-and-go wave observed in the traffic flow due to the amplification of a small disturbance. We simulate the traffic flow under the periodic boundary conditions. The following variation of the initial density ρ
0 is used as in [15]:
(24)
ρ
(
x
,
0
)
=
ρ
0
+
Δ
ρ
0
{
cosh
⁡
-
2
[
160
L
(
x
-
5
L
16
)
]
-
1
4
cosh
⁡
-
2
[
40
L
(
x
-
11
L
32
)
]
}
,
where L
=
32.2 km is the length of the road section under consideration. The periodic boundary conditions are given as follows:
(25)
ρ
(
L
,
t
)
=
ρ
(
0
,
t
)
,
v
(
L
,
t
)
=
v
(
0
,
t
)
.
Here we use the equilibrium speed-density relationship proposed in [22]:
(26)
V
e
(
ρ
)
=
v
f
{
[
1
+
exp
⁡
(
ρ
/
ρ
m
-
0.25
0.06
)
]
-
1
-
3.72
×
1
0
-
6
}
.
Assume the initial flow to be in local steady state; that is, v
(
x
,
0
)
=
V
e
(
ρ
(
x
,
0
)
). Let Δ
ρ
0
=
0.01veh/m, let the space interval Δ
x be 100 m, and let the time interval Δ
t be 1 s. The choice of Δ
x and Δ
t satisfies the Courant-Friedrichs-Levy (CFL) stability condition. The other parameter values we take are as follows:
(27)
τ
=
5
s
,
T
=
10
s
,
v
f
=
30
m
/
s
,
c
0
=
11
m
/
s
,
ρ
m
=
0.2
veh
/
m
.Figures1 and 2 show the evolution of initial uniform traffic flow under the small disturbance.Temporal evolution of traffic flow before 40 minutes for differentρ
0 with β
=
0.4.
(a)
ρ
0
=
0.03
veh
/
m
(b)
ρ
0
=
0.042
veh
/
m
(c)
ρ
0
=
0.05
veh
/
m
(d)
ρ
0
=
0.06
veh
/
m
(e)
ρ
0
=
0.08
veh
/
mFigure1 shows the temporal evolution of traffic flow before 40 minutes for different ρ
0 with β
=
0.2. In pattern (a) of Figure 1, the initial density ρ
0
=
0.03 of traffic flow is so low that the disturbance dies out without any amplification with time. With the initial density ρ
0 increasing, small disturbance is amplified and so leads to the instability of traffic flow. The pattern (b) in Figure 1 shows that several local clusters form for the initial density ρ
0
=
0.042. In pattern (c) of Figure 1, the stop-and-go phenomenon, which is a complex local structure consisting of multiclusters, can be observed. Continuing to increase the initial density ρ
0, we can see a dipole-like structure which is illustrated by Figure 1(d). In pattern (e) of Figure 1, a stable traffic flow is reached again with the initial density ρ
0
=
0.08.Figure2 shows the temporal evolution of traffic flow before 40 minutes for different ρ
0 with β
=
0.4. In patterns of Figure 2, there are no clusters and the disturbances propagate backward without any amplification with different initial densities ρ
0.Comparing every subfigure of Figures1 and 2, the value of β is the only different parameter. It is obvious that the increase of β leads to the stabilization of traffic flow. So the driver’s forecast effect has the positive effect of reducing the local cluster. This result is coincident with the one in [21].The above numerical results are not so perfect. The propagation speeds of the perturbation waves we get from our figures are about 40–60 km/h against the traffic direction, which are not very consistent with the ones from the real data. References [23, 24] point out that the propagation velocity which is obtained from the real data should be about -
12–20 km/h. Many factors may lead to this result such as the option of equilibrium function, the difference scheme, or the parameters. The control and optimization problems of the parameters in our model can also be studied [25–30]. We will continue to study these problems in future.


## 5. Summary

In this paper, a new macroversion is proposed with the consideration of the driver’s forecast effect proposed. The stability condition of the model is derived by using linear analysis. The triangular shock wave which is determined by Burgers equation in the stable region is discussed with reductive perturbation method. The stop-and-go phenomena appear when we carry out the numerical simulation for the model. The driver’s forecast effect which has the positive effect of reducing the local cluster is proved.



---

*Source: 539459-2014-08-24.xml*