# The Car Following Model with Relative Speed in Front on the Three-Lane Road.

1. Introduction

In this paper, the new LRVD model is proposed, which studies the cars behavior on the middle lane of one road with three lanes. In Section 2, the LRVD model is elaborated in detail. In Section 3, the linear analysis theory of LRVD model, including the neutral curve, is presented. In Section 4, by studying the nonlinear theory, we get the coexisting curve of the model, which divides the traffic flow into stable, metastable, and unstable regions. In Section 5, simulation results also prove LRVD model indeed enhances the stability of the traffic flow. The last section concludes some remarks of LRVD model.

2. Model

At the very beginning, Chowdhury et al. [4] proposed car following model on one single road just as follows:

d[x.sup.2.sub.n]/d[t.sup.2] = [f.sub.sti]([v.sub.n], [DELTA][x.sub.n], [DELTA][v.sub.n]) , (1)

where the function [f.sub.sti](x) represents the drivers' response to the stimulus, which is related to the speed of nth car [v.sub.n], the relative speed between the cars in front of the nth car and the nth car, and the distance between them.

Newell [6] proposed a simple single road car following model.

d[x.sub.n] (x + [tau])/dt = V([DELTA][x.sub.n](t)), (2)

where V([DELTA][x.sub.n](t)) is the optimal velocity function, which just contains the distance between the nth car and the car in front. [tau] is the delay time; the driver of the nth car can respond to the stimulus of the (n + 1)th car in front after the time lag [tau].

Yu [8] considered that the optimal velocity may be determined not only by the successive space but also by the response to the stimulus of the relative velocity. In this case, d[x.sub.n](t + [tau])/dt = V([DELTA][x.sub.n](t), [DELTA][v.sub.n] (t)), which further describes the car following process as follows, namely, RV model:

[d.sup.2][x.sub.n](t)/d[t.sup.2] = [alpha][V([DELTA][x.sub.n](t)) - [v.sub.n](t)] + [k.sub.1][DELTA][v.sub.n](t). (3)

Yu et al. [17] proposed the LVD model based on the RV model, which considers the impact from the relative velocity on the main and adjacent lanes. By expanding (2) into d[x.sub.n](t + [tau])/dt = V([DELTA][x.sub.n](t), [DELTA][v.sub.n](t), [DELTA][v.sub.n.sup.LVD](t)), and assuming V([DELTA][x.sub.n](t), [DELTA][v.sub.n](t), [DELTA][v.sub.n][sup.LVD](t)) = V([DELTA][x.sub.n](t)) + [[lambda].sub.1][DELTA][v.sub.n](t) + [[lambda].sub.2][tau][DELTA][v.sub.n][sup.LVD](t), the LVD model is obtained as follows:

[d.sup.2][x.sub.n](t)/d[t.sup.2] = [alpha][V([DELTA][x.sub.n](t)) - [v.sub.n](t)] + [k.sub.1][DELTA][v.sub.n](t) + [[lambda].sub.2][DELTA][v.sub.n.sup.LVD](t). (4)

The scholars have improved the car following model in many aspects before; however no paper was taken with regard to the cars behavior on the middle lane of the three-lane road. In fact, one road with three same lanes, as Figure 1 has shown, is more common than one single road or two-lane road in real life. Considering the complex phenomena above, the LRVD model is proposed in this paper. Equation (2) can be expanded into such differential equation:

d[x.sub.n](t + [tau])/dt = V([DELTA][x.sub.n](t), [DELTA][v.sub.n](t), [DELTA][v.sup.LVD.sub.n](t), [DELTA][v.sup.RVD.sub.n](t)), (5)

where [DELTA][v.sup.LVD.sub.n](t) = mean([v.sup.L.sub.m], [v.sup.L.sub.m+1],..., [v.sup.L.sub.m+s-1]) - [v.sub.n], and [DELTA][v.sup.RVD.sub.n](t) = mean([v.sup.R.sub.k], [V.sup.R.sub.k+1],..., [v.sup.R.sub.k+s-1]) - [v.sub.n]. [v.sub.n] is the speed of the nth car on the main lane. s is the car number on the left lane and the right lane. The mean([v.sup.L.sub.m], [v.sup.L.sub.m+1],..., [v.sup.L.sub.m+s-1]) is the average speed of all the cars on the left lane. The mean([v.sup.R.sub.k], [v.sup.R.sub.k+1],..., [v.sup.R.sub.k+s-1]) is the average speed of all the cars on the right lane.

Here we suppose the linear combination of V([DELTA][x.sub.n](t), [DELTA][v.sub.n](t), [DELTA][v.sup.LVD.sub.n](t), [DELTA][v.sup.RVD.sub.n](t)):

[mathematical expression not reproducible]. (6)

Then after applying Talor expanding to (5), we get the following model:

[d.sup.2][x.sub.n](t)/d[t.sup.2] = [alpha][V([DELTA][x.sub.n](t)) - [v.sub.n](t)] + [k.sub.1][DELTA][v.sub.n](t) + [[lambda].sub.2][DELTA][v.sub.n.sup.LVD](t) + [[lambda].sub.3][DELTA][v.sup.RVD.sub.n](t), (7)

where [k.sub.1] = [[lambda].sub.1]/[tau] and [[lambda].sub.1] is the response coefficient of the relative speed between the nth car and the (n + 1)th car in front. [[lambda].sub.2] is the response coefficient of the relative speed between cars on the main and left adjacent lane, which are just in front of the nth car. [[lambda].sub.3] is the response coefficient of the relative speed between cars on the main and right adjacent lanes, which are just in front of the nth car.

V([DELTA][x.sub.n](t)) is the optimal velocity of nth car at time t, which is defined as follows:

V([DELTA][x.sub.n]) = [v.sub.max]/2 (tanh ([DELTA][x.sub.n] - [h.sub.c]) + tanh ([h.sub.c])), (8)

where [h.sub.c] is the safe distance of successive cars and [v.sub.max] is the maximum velocity. Note that when we ignore the relative speed on the right adjacent lane, namely, [[lambda].sub.3] = 0, the LRVD model proposed can be degenerated to the LVD model. Then if we also ignore the relative speed on the left adjacent lane and the right adjacent lane, [[lambda].sub.2] = [[lambda].sub.3] = 0, the LRVD model proposed can be degenerated to the RV model. Further, if we continue to ignore the relative speed on the main lane, namely, [[lambda].sub.1] = [[lambda].sub.2] = [[lambda].sub.3] = 0, the new LRVD model can be degenerated into the original OV model. So the theory continuity of the OV, RV, LVD and LRVD models indicates the LRVD model proposed is reasonable and it can indeed explain some complicated phenomenon on three-lane road that scholars have not studied before. In a word, LRVD model is the basement of the three-lane road traffic phenomenon.

The dynamic equation (7) can be discretized by the asymmetric forward difference, so we get

[mathematical expression not reproducible]. (9)

Here, [[lambda].sub.1], [[lambda].sub.2], and [[lambda].sub.3] are defined as follows:

[mathematical expression not reproducible], (10)

where [pi], [eta], [phi] are the given response coefficients and [s.sub.c] is a critical successive space. If the distance between the two successive cars on the main lane is more than [DELTA][x.sub.n], the speed difference of these two cars would not affect the following car's behavior at all. If the distance between the car on the main lane and the car in front of it on the left adjacent lane is beyond the critical successive space, these two cars' relative speed can be ignored. If the distance between the car on the main lane and the car in front of it on the right adjacent lane is more than the critical successive space [s.sub.c], these two cars' relative speed can be ignored together.

Eq. (9) can be rewritten as

[mathematical expression not reproducible]. (11)

3. Linear Stability Analysis

The hypotheses of LRVD model are (a) cars cannot change their lane, (b) the three lanes on one road are the same, (c) neither overtaking nor ramp exists in this model, and (d) all the cars on the right and left adjacent lanes can be affected by the cars in front on the middle lane and also on their own respective lanes.

Here we study the uniform traffic flow stability on each lane, which is given as the solution of the uniformly steady state of (9), which is shown as

[mathematical expression not reproducible], (12)

where [N.sub.1], [N.sub.2] and [N.sub.3] represent the number of cars on the main lane, the left adjacent lane, and right adjacent lane; h, [h.sub.l], and [h.sub.r] are the successive space on the main lane, the left adjacent lane, and the right adjacent lane; L is the length of the road; V(h), [V.sub.l]([h.sub.l]), and [V.sub.r]([h.sub.r]) indicate the optimal speed of the main lane, the left adjacent lane, and the right adjacent lane; [x.sup.(0).sub.n](t), [x.sup.L(0).sub.m](t), and [x.sup.R(0).sub.k](t) are the position of the nth car on the main lane, the mth car on the left adjacent lane, and the kth car on the right adjacent lane at t = 0 in the stable state.

Considering the small perturbation into the steady state solutions in (12), [x.sub.n](t) = [x.sup.(0).sub.n](t) + [y.sub.n](t), [x.sup.L.sub.m](t) = [x.sup.L(0).sub.m](t) + [y.sup.L.sub.m](t), and the linear stability equation can be obtained:

[mathematical expression not reproducible], (13)

where V'(h) = dV([DELTA]x)/d[DELTA]x. Expanding [DELTA]y into the (13) with [DELTA][y.sub.n](t) = Aexp{ikn + zt},

[DELTA][y.sup.L.sub.m](t) [congruent to] [gamma]A exp {ikn + zt}, [DELTA][y.sup.R.sub.k](t) [congruent to] [gamma]A exp {ikn + zt}, (14)

where [gamma] represents the average strength coefficient which describes the perturbations between the main lane and the adjacent lanes. Then we get (13) as follows:

[mathematical expression not reproducible]. (15)

Inserting z = [z.sub.1](ik) + [z.sub.2][(ik).sup.2] + ... into (15), we obtain

[z.sub.1] = V' + ([[lambda].sub.2] + [[lambda].sub.3]) ([gamma] - 1), [z.sub.2] = -3/2 [z.sup.2.sub.1][tau] + V'/2 + [[lambda].sub.1][z.sub.1] + ([[lambda].sub.2] + [[lambda].sub.3]) [gamma]s - 1/2. (16)

The flow is unstable if [z.sub.2] < 0, and stable if [z.sub.2] > 0. So the neutral stability condition is given by [z.sub.2] = 0:

[tau] = V' + 2[[lambda].sub.1][z.sub.1] + ([[lambda].sub.2] + [[lambda].sub.3])([gamma]s - 1)/3[z.sup.2.sub.1]. (17)

Considering the inverse of delay time [tau] which is the critical sensitivity [[alpha].sub.c], the neutral stability curve is shown as follows:

[[alpha].sub.c] = 3 [[V' + ([[lambda].sub.2] + [[lambda].sub.3])([gamma] - 1)].sup.2]/ V' + 2[[lambda].sub.1] [V' + ([[lambda].sub.2] + [[lambda].sub.3]) ([gamma] - 1)] + ([[lambda].sub.2] + [[lambda].sub.3]) ([gamma]s - 1). (18)

The unstable and stable areas are

[mathematical expression not reproducible]. (19)

In Figure 2, suppose that the OV model with [[lambda].sub.1] = [[lambda].sub.2] = [[lambda].sub.3] = 0, RV model with [[lambda].sub.1] = 0.2, [[lambda].sub.2] = [[lambda].sub.3] = 0, the LVD model with [[lambda].sub.1] = [[lambda].sub.2] = 0.2, [[lambda].sub.3] = 0, [gamma] = 0.4, the LRVD model with [[lambda].sub.1] = [[lambda].sub.2] = [[lambda].sub.3] = 0.2, [gamma] = 0.4, s = 1, and [v.sub.max] = 3 m/s. Through Figure 2, we know the stable area of LRVD model is the biggest compared to others, which means considering the cars' relative speed on the third adjacent lane indeed improves the traffic flow stability. So it is necessary to propose the new LRVD model in this paper, which also can describe the real traffic phenomenon on the three-lane road in real life.

As Figure 3 has shown, if drivers on the middle lane pay more attention to the cars behavior on the two adjacent lanes, the stable area is much bigger. The more attention is paid, the more stable the area is. So if the drivers on the middle lane pay more attention to the cars' speed on the two side lanes, the traffic flow on the middle lane will be more stable.

As Figure 4 has shown, if the drivers on the middle lane pay the same attention to the two other lanes, the more cars on the two adjacent lanes are considered by the drivers on the middle lane and the more stable the traffic flow on the middle lane will be.

Through Figures 3 and 4, we realize the cars' driving truth on the three-lane road. Firstly, if the drivers on the middle lane pay more attention to the other lanes, rather than only see the cars on middle lane or one adjacent lane just like the existing model only on the single road or two-lane road, the traffic flow on the middle lane must be more stable. Secondly, in the same condition of the middle lane, if more cars in front on the two adjacent lanes are noticed by the driver on the middle lane, the middle lane traffic flow is certain to be more stable.

4. Nonlinear Stability Analysis

To prove that the relative speed on the three lanes on one road indeed affects the traffic flow on the middle lane effectively, we conduct a nonlinear stability of the LRVD model. Here the nonlinear method is used to study the slowly varying behavior of long waves near the critical point ([h.sub.c], [a.sub.c]), and the slow scales for space variable n and time variable t.The slow variables X and T are defined as follows:

X = [epsilon](n + bt), T = [[epsilon].sup.3]t, 0 < [epsilon] [much less than] 1, (20)

where b is a constant to be decided later. Considering a small disturbance [epsilon]R(X, T), the distances between the two successive cars can be rewritten as follows:

[DELTA][x.sub.n](t) = [h.sub.c] + [epsilon]R(X, T). (21)

Then expanding (11) to the [epsilon] fifth order by substituting (20) and (21),we can get the nonlinear partial differential equation as follows:

[mathematical expression not reproducible]. (22)

Near the critical point ([h.sub.c], [a.sub.c]), (22) is written as

[mathematical expression not reproducible], (23)

where

[mathematical expression not reproducible]. (24)

Then a transformation for (23) is shown as follows:

T' = [g.sub.1]T, R = [square root of ([g.sub.1]/[g.sub.2])] R'. (25)

Here, the regularized equation of (23) is obtained:

[mathematical expression not reproducible]. (26)

Ignoring the O([epsilon]) term in (26), we get the modified KdV equation

[R'.sub.0](X, T') = [square root of (c)] tanh [square root of (c/2)](X - cT'), (27)

where c is the propagation velocity of the kink wave, which is determined by O([epsilon]) term. In order to decide the value of c, the solvability condition must be satisfied:

([R'.sub.0], M[[R'.sub.0]]) = [[integral].sup.[infinity].sub.-[infinity]] dX[R'.sub.0](X, T') M[[R'.sub.0](X, T')] = 0, (28)

where

[mathematical expression not reproducible]. (29)

After integration, we get the selected speed c = 5[g.sub.2][g.sub.3]/ (2[g.sub.2][g.sub.5] - 3[g.sub.1][g.sub.4]). Then, the solution of the modified KdV equation is obtained

R(X, T) = [square root of ([g.sub.1]c/[g.sub.2])] tanh [square root of (c/2)](X - c[g.sub.1]T). (30)

From (8), we know V' = [v.sub.max]/2, V'" = -[v.sub.max]. The amplitude of the kink solution is given by A = [[([g.sub.1]c/[g.sub.2])([a.sub.c]/a - 1)].sup.1/2] with

[mathematical expression not reproducible]. (31)

Through the mKdV equation, the coexisting curves are obtained. The coexisting phase consists of congested phase with high density and freely moving phase with low density. The successive space in freely moving phase is [DELTA]x = [h.sub.c] + A, and that in congested phase is [DELTA]x = [h.sub.c] - A. The coexisting curve and neutral stability curve divide the whole area into stable, metastable, and unstable areas, just as Figure 2 has shown.

5. Simulation

In this section, we use simulation results to describe the LRVD model. The initial parameters are given as follows:

On the main lane [v.sub.max] = 4 m/s, [h.sub.c] = 4m, a = 2[s.sup.-1], [v.sub.n](0) = 2 m/s, N = 50, and T = 0.1 s, where N is the car number and T is the time step.

On the left and right side lanes: [v.sup.L.sub.m](0) = [v.sup.R.sub.k](0) = 2 m/s, [N.sub.L] = [N.sub.R] = 50, and [h.sub.m] = [h.sub.k] = 4m; all the cars on the two side lanes and the main lane keep the same characteristics during the whole simulation. All the cars are driving at the speed of 2 m/s until 100 s on their own lanes. From 100 s to 103 s, the leading car on the main lane stops for three seconds as a disruption to the traffic flow, and then the leading car runs later at the recovery speed of 2 m/s.

Here in the SLVD model, [[lambda].sub.1] = [[lambda].sub.2] = [[lambda].sub.3] = 0.2. In the LVD model, [[lambda].sub.1] = [[lambda].sub.2] = 0.2, [[lambda].sub.3] = 0. Figure 5 shows the speed fluctuation of the 3rd, 25th, and 50th cars after 100 s, when [[lambda].sub.2] = [[lambda].sub.3] = 0.2 in LRVD model and [[lambda].sub.2] = 0.2, [[lambda].sub.3] = 0 in the LVD model.

When perturbation is set in the stable traffic flow, it will be amplified with time.The traffic jam, go-and-stop wave and speed fluctuation will all appear in the simulation. In Figure 5, the speed fluctuation in LRVD model is smaller than that in LVD model, which means the stability of LRVD model has been improved compared to the LVD model. The cars run with less perturbation effect when considering the relative speed on the other adjacent lane in LRVD model.

To explore the improved stability of LRVD model, we pay attention to the distances between the first car and every following car behind it on the main lane. Figure 6 shows the distance-time patterns of LRVD model and LVD model, from which we realize the leading car' perturbation propagates backward after it stops at 100 s. From Figure 6, we know the LRVD model recovers its stable state much earlier than LVD model. The distance between the 50th car and the leading car in LVD model recovers the distance in stable state at about 145 s; but the distance between the 50th car and the leading car in LRVD model recovers its stable distances at 141 s. In fact, the traffic flow will recover its stable state with different time, just because of the different stability of the different model. The LVD model needs more time, but LRVD model needs less time. This simulation results mean, under the same traffic situation in real life, the LRVDismore stable than LVD model.

Results of the simulation shown in Figure 6 are described below (reaction time [tau] is 0.1 s).

(1) When t = 100 s, the first car stops; then all the distances between the first car and each car that follows start decreasing. At that time, the driver of the second car cannot react from the first car's stopping. So all the following cars, from the second car to the last one, go on running at the uniform linear motion on the main lane with the speed of 2 m/s.

(2) When t = 100.1 s, the second car reacts to the first car stopping and starts to decelerate. Because of the time delay [tau], the second driver decides its acceleration value according to the acceleration formula with the relative speed, relative distance, and the second car's speed information and so on at the last time. Here the second driver gets its deceleration value from the relative speed, its own speed, and relative distance with the leading car at t = 100 s. Then, based on the acceleration value, the second driver moves according to the [v.sub.100.1 s] = [v.sub.100 s] + at, [x.sub.100.1 s] = [x.sub.100 s] + [v.sub.100 s] t + (1/2)a[t.sup.2] from t = 100 s to t = 100.1 s, by which we get all the relevant parameters value at t = 100.1 s. At t = 100.1 s, [v.sub.1] = 0, [v.sub.2] < [v.sub.3], [v.sub.3] = 2 m/s, and [DELTA]x(2) < [DELTA]x(3) < 4 m. ([DELTA]x(i) = x(i) - x(i - 1)). The remaining cars, from the third car to the last one, go on running at the uniform linear motion with a speed of 2 m/s.

(3) When t = 100.2 s, the third car realizes the second car has decelerated, so it starts to decelerate. At this time, [DELTA]x(2) < [DELTA]x(3) < [DELTA]x(4) < 4 m, [DELTA]x(i) = x(i) - x(i - 1), [a.sub.1] = 0, [a.sub.2] < [a.sub.3] < 0. The remaining cars, from the fourth car to the last one, go on running at the uniform linear motion with a speed of 2 m/s.

(4) The following ith car decelerates according to the (i - 1)th car ahead deceleration so that [DELTA]x(2) < [DELTA]x(3) < ... < [DELTA]x(i + 1) < 4 m, and [a.sub.1] = 0, [a.sub.2] < [a.sub.3] < ... <[a.sub.i] < 0, [a.sub.i+1] = 0, [v.sub.1] = 0, [v.sub.2] < [v.sub.3] < ... < [v.sub.i] < 2 m/s, [v.sub.i+1] = 2 m/s. During the time interval from 100 s to 103 s, we may expect 30 cars to be affected by this disturbance. However, the disturbance does not spread so fast: in fact, the cars in front, from the second one to the fifth car, are affected every [tau] until 100.4 s because of [absolute value of [a.sub.i]] > 0. But the speed of the cars, which are behind jth car (j > i, here j = 6 in our simulation), will change much less every [tau] time, because [absolute value of [a.sub.j] [approximately equal to] 0, [absolute value of [a.sub.j]t [approximately equal to] 0, [v.sub.j] = [v.sub.j,0] + [a.sub.j]t. In this case, the successive distance after jth car even seems to be not changed. But as time passes by, during the 100 s to 103 s, the cars in front keep decelerating, so the following cars will repeat the same decelerating process above; in fact, the decelerating wave spreads backward.

(5) At t = 103 s, only 11 cars in front decelerate, while all the following cars keep with the same speed of 2 m/s and the successive space of 4 m in the stable state. All the distances between the first car and each following car must decrease during the whole three seconds so that [DELTA]x(2) < [DELTA]x(3) < ... < [DELTA]x(12) < 4 m, and [DELTA]x(13) = ... = [DELTA]x(50) = 4 m. At t = 103 s, we assume that the first car instantaneously recovers its stable speed of 2 m/s; soon before this instant, all the distances between the first car and each following car have reached their minimum value.

(6) After the first car has recovered the stable speed (2 m/s), the second car reacts and begins to accelerate at the 103.1 s, but its speed being lower than 2 m/s, [DELTA]x(2) keeps increasing until [DELTA]v(2) = 0. After this time, [DELTA]x(2) recovers the stable value of 4 m with some fluctuation in the end.

(7) Later, at t = 103.2 s, the third one realizes the second car begins to increase its speed, but according to the adopted numerical values of the parameters, the simulation shows initially that v(2) < v(3), despite the fact that the third car keeps decelerating while the second car is accelerating until v(2) = v(3). At that time, the 3rd car acceleration is zero, which has reached its minimum value; then, the third car begins to accelerate because of v(2) > v(3), and [DELTA]v(3) keeps increasing until the third car recovers its stable speed of 2 m/s. Then, after some fluctuation, [DELTA]x(3) will recover the stable distance of 4 m, but this phenomenon will happen after [DELTA]x(2) finishes recovering the same stable distance. From 103 s, the distance between the first car and third car keeps increasing from the minimum value up to reach the stable state with some fluctuation in the end.

(8) The 11 cars in front, which have already decreased their speed, will repeat similar process as shown in part 7. Take the ith car as an example; the ith car decelerates firstly because of v(i - 1) < v(i) until v(i - 1) = v(i). At this time, [a.sub.i] = 0, the successive space [DELTA]x(i) decreases at the lowest distance. Then, with the (i - 1)th car accelerating, ith car keeps accelerating, v(i - 1) > v(i). Finally [DELTA]x(i) keeps increasing up to the stable value of 4 m, which occurs after the (i-1)th car has recovered the same distance with respect to the (i - 2)th car.

(9) Other following cars behind the 11th car have not decelerated yet. So, each car has to repeat the same process as 11 cars ahead, in the deceleration phase firstly during 100 s to 103 s, and after 103 s keeps decelerating at first till the relative speed increases to zero and begins to accelerate later until reaching the stable state, as parts 7 and 8 have described above. For each following car the same process occurs, so the following car indeed needs more time to react to the behavior of the car ahead, which is in compliance with the kinematic wave theory.

(10) After 103 s, when each following car begins to decelerate, the distance between the first car and this following car starts to increase at once up to reach the stable state distance with some fluctuation in the end. So the recovery wave in Figure 6 propagates obliquely up and back until the last car recovers the stable state.

6. Conclusions

In this paper, we presented a new LRVD model through considering the relative speed in front on the three-lane road. The original OV model has been developed into RV, LVD, and especially LRVD model in this paper, which can portray more complex realistic traffic phenomenon of the cars driving behavior on the middle lane of the three-lane road. The linear analysis is studied to get the neutral stability curve. The nonlinear analysis is studied to obtain the mKdV equation and the coexisting curve, by which the area is divided into three regions of stable, metastable, and unstable. The simulation results about distance-time profile are accordant with the analysis. The stability of traffic flow of LRVD model has been proved to be improved compared to the LVD model on two lanes. The LRVD model shows if the drivers on the middle lane pay more attention to more cars on the two adjacent lanes in real life, the middle lane traffic flow is certain to be more stable. All the methods, including the linear, nonlinear analyses, and simulation prove that LRVD model with three lanes indeed stabilize traffic flows much better and capture traffic flow characteristic on the three-lane road in real life more effectively. On the complex three-lane road, if intelligent traffic management system based on the huge traffic data for drivers is applied in the real life, it is very helpful to ensure traffic safety, which is also the trend of transportation development in future.

https://doi.org/10.1155/2018/7560493

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

References

[1] M. Bando, K. Hasebe, A. Nakayama, A. Shibata, and Y. Sugiyama, "Dynamical model of traffic congestion and numerical simulation," Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 51, no. 2, pp. 1035-1042, 1995.

[2] M. Bando, K. Hasebe, K. Nakanishi, and A. Nakayama, "Analysis of optimal velocity model with explicit delay," Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 58, pp. 5429-5435, 1998.

[3] D. Helbing and B. Tilch, "Generalized force model of traffic dynamics," Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 58, no. 1, pp. 133-138, 1998.

[4] D. Chowhury, L. Santen, and A. Schadschneider, "Statistical physics of vehicular traffic and some related systems," Physics Reports, vol. 329, pp. 199-329, 2000.

[5] R. Jiang, Q. Wu, and Z. Zhu, "Full velocity difference model for a car-following theory," Physical Review E, vol. 64, no. 1, pp. 1-4, 2001.

[6] G. F. Newell, "Nonlinear effects in the dynamics of car following," Operations Research, vol. 9, pp. 209-229, 1961.

[7] Z. Liang, S. Zhong-ke, and L. Tong, "An improved FVD model for intersection," Journal of Transportation Systems Engineering and Information Technology, vol. 13, no. 6, pp. 80-85, 2013.

[8] X. Yu, "Analysis of the stability and density waves for traffic flow," Chinese Physics, vol. 11, no. 11, pp. 1128-1134, 2002.

[9] T. Q. Tang, C. Y. Li, and H. J. Huang, "A new car-following model with the consideration of the driver's forecast effect," Physics Letters A, vol. 374, no. 38, pp. 3951-3956, 2010.

[10] M. Treiber, A. Hennecke, and D. Helbing, "Derivation, properties, and simulation of a gas-kinetic-based, nonlocal traffic model," Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 59, no. 1, pp. 239-253, 1999.

[11] M. Treiber, A. Hennecke, and D. Helbing, "Congested traffic states in empirical observations and microscopic simulations," Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 62, no. 2, pp. 1805-1824, 2000.

[12] T. Q. Tang, Y. P. Wang, X. B. Yang, and Y. H. Wu, "A new car-following model accounting for varying road condition," Nonlinear Dynamics, vol. 70, no. 2, pp. 1397-1405, 2012.

[13] T. Q. Tang, W. F. Shi, X. B. Yang, Y. P. Wang, and G. Q. Lu, "A macro traffic flow model accounting for road capacity and reliability analysis," Physica A: Statistical Mechanics and its Applications, vol. 392, no. 24, pp. 6300-6306, 2013.

[14] X.-B. Yang, B.-F. Si, and M. Huan, "Mixed traffic flow modeling near Chinese bus stops and its applications," Journal of Central South University, vol. 19, no. 9, pp. 2697-2704, 2012.

[15] S. Li, Z. Shi, and T. Liu, "An improved model based on FVD model and experimental comparison," Journal of Highway and Transportation Research and Development, vol. 31, no. 2, pp. 109-113, 2014.

[16] T. Nagatani, "Stabilization and enhancement of traffic flow by the next-nearest-neighbor interaction," Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 60, no. 6, pp. 6395-6401, 1999.

[17] G. Yu, P. Wang, X. Wu, and Y. Wang, "Linear and nonlinear stability analysis of a car-following model considering velocity difference of two adjacent lanes," Nonlinear Dynamics, vol. 84, no. 1, pp. 387-397, 2016.

[18] J. Chen, Z. Peng, and Y. Fang, "Effects of car accidents on three-lane traffic flow," Mathematical Problems in Engineering, vol. 2014, Article ID413852, 11 pages, 2014.

Zichu Gao (iD), (1) Ning Zhang, (1) Livia Mannini, (2) and Ernesto Cipriani (2)

(1) School of Economics and Management, Beihang University, Beijing 100083, China

(2) Department of Engineering, Roma Tre University, 00185 Rome, Italy

Correspondence should be addressed to Zichu Gao; zhuzhusenlin@163.com

Received 3 November 2017; Accepted 19 February 2018; Published 10 April 2018

Caption: Figure 1: Three-lane road traffic situation.

Caption: Figure 2: Neutral stability curves and coexisting curves of OV, RV, LVD, and LRVD models.

Caption: Figure 3: Different stable curves with different [[lambda].sub.2] and [[lambda].sub.3] in LVRD model.

Caption: Figure 4: Different stable curves with different car numbers on the two adjacent lanes in LVRD model.

Caption: Figure 5: (Color online) speed-time profile of the LRVD model (a) and LVD model (b) after the first car stops from the 100 s to 103 s. (a) [[lambda].sub.1] = [[lambda].sub.2] = [[lambda].sub.3] = 0.2 (b) [[lambda].sub.3] = 0, [[lambda].sub.1] = [[lambda].sub.2] = 0.2. As time passes by, the red line represents the speed fluctuation of the 3rd car behind the first car, the blue line represents the speed fluctuation of the 25th car, and the green line represents the speed fluctuation of the 50th car.

Caption: Figure 6: (Color online) distance-time profile of the LRVD model (a) and LVD model (b). The first car stops from 100 s to 103 s with (a) LRVD model [[lambda].sub.1] = [[lambda].sub.2] = [[lambda].sub.3] = 0.2 (b) [[lambda].sub.3] = 0, [[lambda].sub.1] = [[lambda].sub.2] = 0.2. The red line shows the end of the distance change. The LRVD needs 141 s, and LVD model needs 145 s.