A new control method based on the lattice hydrodynamic model considering the double flux difference
Qin Shunda1, 2, 3, Ge Hongxia1, 2, 3, Cheng Rongjun1, 2, 3, †
Faculty of Maritime and Transportation, Ningbo University, Ningbo 315211, China
Jiangsu Province Collaborative Innovation Center for Modern Urban Traffic Technologies, Nanjing 210096, China
National Traffic Management Engineering and Technology Research Center Ningbo University Sub-centre, Ningbo 315211, China

 

† Corresponding author. E-mail: chengrongjun@nbu.edu.cn

Abstract

A new feedback control method is derived based on the lattice hydrodynamic model in a single lane. A signal based on the double flux difference is designed in the lattice hydrodynamic model to suppress the traffic jam. The stability of the model is analyzed by using the new control method. The advantage of the new model with and without the effect of double flux difference is explored by the numerical simulation. The numerical simulations demonstrate that the traffic jam can be alleviated by the control signal.

1. Introduction

With the development of cities, peopleʼs daily lives are deeply affected by traffic problems, especially traffic congestionm which is one of the most serious problems. Many experts and scholars have paid more and more attention to alleviating traffic jams. A variety of traffic flow models[118] have been studied, such as the lattice hydrodynamic model[47] and the car-following model,[812] to work out the formation of the traffic jam and reduce traffic congestion. All models can be categorized as either a macroscopic model,[4,5] microscopic model,[9,10] or mesoscopic model.[19,20]

In 1961, a famous car-following model was proposed by Newell[8,9] with an optimal velocity (OV) function, which has a profound influence on later research. Scholars began to build models by considering the OV function to describe the traffic problems. In 1995, Bando et al.[21,22] put forward a classic car-following model called the optimal velocity model (OVM). The OVM shows that the driver adjusts the velocity based on the headway distance and exposes the essence of the traffic congestion. Later, many car-following models were developed based on the OVM. In 2006, Zhao et al. proposed a coupled-map (CM) car-following model,[23,24] and the velocity feedback control[25,26] was taken into account. According to the CM model, it is helpful to suppress the traffic congestion considering the feedback control signal in the models.

The lattice hydrodynamic model is a macroscopic model, which was put forward by Nagatani[27] based on the theory of Kerner and Konhauser.[28] Since then, more and more scholars have utilized the lattice hydrodynamic method to analyze the macroscopic traffic flow. The lattice hydrodynamic model draws lessons from the car-following model with a similar description of the density waves and stability analysis. In 2015, Ge[29] proposed a new lattice hydrodynamic model considering the feedback signal called flux difference to suppress the traffic jam in view of modern control theory.[30] In the same year, Redhu[31] presented a delayed-feedback control (DFC) model with the consideration of the driving behavior. Subsequently, the effect of the density change rate difference was considered as a feedback signal by Li.[32] In 2016, Zhu[33] proposed a new control signal called the variation rate of the optimal velocity in the lattice hydrodynamic model.

In the lattice hydrodynamic model, the traffic flow is divided into lots of lattices by the discretization method. Obviously, during the actual driving condition, the influence of multi-lattices is significant. The drivers adjust their behavior based on the flux difference of the preceding and local lattices. In 2012, Wang[34] considered the multiple flux difference effect in his new model. In Wangʼs research, the density wave analysis was adopted to explain the effect of the multiple flux difference. In addition, the flux of the two lattices in front of the local lattice has an important effect on driving, and the effect has diminished since the third lattice. As far as we know, the double flux difference which has an important effect on the real traffic, has not been studied. Up to date, most feedback signal methods were applied to the car-following model to alleviate traffic congestion. The feedback signal was seldom considered in the lattice hydrodynamic model. For these reasons, a new lattice hydrodynamic model with the consideration of the double flux difference is proposed. With the feedback signal of multi-lattices, a new optimal velocity function is proposed to keep the traffic system stable. And we use the control method to explore the effect of the double flux difference on the traffic flow to explore its impact on the traffic flow.

The outline of the paper is as follows. In Section 2, the basic lattice hydrodynamic model is derived, and we introduce the control signal for the flux difference model. In Section 3, a control signal is added to the basic lattice hydrodynamic model and the feedback control theory is used to analyze the stability conditions. In Section 4, several numerical simulations are carried out to verify the theoretical results. Conclusions are given in Section 5.

2. Basic model

The basic lattice hydrodynamic model has the advantages of macroscopic models and car-following model. The equations are organized as follows:[28]

where ρ is the current density at time t and is the average density. δ is the average space headway and x is the position of the local lattice, so is the current density at position . Especially, and δ can be expressed as . a represents the sensitivity coefficient of the driver. We set and . Nagatani creatively put forward the discretization method to explain the equations as follows:
where j represents the location of the lattice. In the above equations, and vj are the current density and the current velocity, respectively. By simplifying the equation, making , and considering the feedback signal, we have
where the control signal uj is given as
with k being the flux difference parameter.

The feedback control signal is similar to the velocity difference feedback signal indicating the flux difference of lattices j and j+1. Nagatani[27] also proposed the new optimal velocity function of the lattice hydrodynamic model

where vmax represents the maximum velocity of a vehicle on this road and is the safety density.

3. Control scheme

In 2015, the control method was first introduced by Ge[29] to consider the flux difference. In reality, the traffic flow is a complex system which consists of many factors. In order to be closer to the actual traffic flow, we propose a lattice hydrodynamic model considering the double flux difference

where , , k is the ratio coefficient of the flux difference of lattices j, j+1, and j+2. Due to the influence of the multi-lattice feedback signal, the new optimal velocity function is given as

Based on the control theory, we compare the expected and real flux and density. and represent the theoretical data of flux and density. The following lattices have the steady state described as

We apply the linear stability criterion to analyze the system considering a small perturbation. The following equation can be derived by the control method:

where , , , and .

Taking the Laplace transformation of Eq. (5), we can easily obtain

where and . L indicates the Laplace transform and s represents the transform function variable. Equation (10) can be written in the matrix form
By simplifying Eq. (11) with the control theory, the transfer function is written as follows:
where the characteristic polynomial . From the deduction of Eq. (12), we can easily obtain

By the Hurwitz stability criterion, we can easily draw a conclusion that when the traffic flow is smooth, characteristic function is greater than 0. To make the system stable, it can be confirmed that is satisfied by the judgement of the control theory and inequalities , , . The following derivation provides an effective measure to work out the multi-variable system stable condition. must be smaller than 1 for all to ensure the stability of the system

The sufficient condition can be derived as

From the above analysis, the sensitivity coefficient of driver a in the traffic system should satisfy

4. Numerical simulations

In this section, several simulations are provided to verify the stability of the feedback control method for the traffic flow system. The original traffic flow is also investigated in our simulation for comparison. We assume that there are N lattice sites on a signal-way road and the specific parameters for the lattice hydrodynamic model are set as follows: N = 140, , , , .

We suppose 280 vehicles in the traffic flow system at the steady state introduced by Eq. (6). Then we set the initial density of sites 50–55 as 0.5 and sites 56–60 as 0.2. The densities of the 5th, 25th, 55th and 80th lattices are depicted in detail.

Figure 1 presents the original lattice hydrodynamic model and Figure 2 is an extended lattice hydrodynamic model with the consideration of the double flux difference. Comparing Fig. 1 with Fig. 2, it is clear that when the double flux difference is considered in the lattice hydrodynamic model, the traffic congestion is suppressed a great deal. From Figs. 25, with flux parameter λ increasing, the density–time curve becomes smoother. According to the patterns of the density in Figs. 58, we can conclude that when ratio coefficient k grows to 0.85, the traffic flow can become more stable. The traffic system will be most stable with k = 0.85 and λ = 0.9.

Fig. 1. (color online) The density–time curve of the 5th, 25th, 55th and 80th lattices without the control signal.
Fig. 2. (color online) The density–time curve of the 5th, 25th, 55th and 80th lattices with λ = 0.6 and k = 0.85.
Fig. 3. (color online) The density–time curve of the 5th, 25th, 55th and 80th lattices with λ = 0.7 and k = 0.85.
Fig. 4. (color online) The density–time curve of the 5th, 25th, 55th and 80th lattices with λ = 0.8 and k = 0.85.
Fig. 5. (color online) The density–time curve of the 5th, 25th, 55th and 80th lattices with λ = 0.9 and k = 0.85.
Fig. 6. (color online) The density–time curve of the 5th, 25th, 55th and 80th lattices with λ = 0.9 and k = 0.75.
Fig. 7. (color online) The density–time curve of the 5th, 25th, 55th and 80th lattices with λ = 0.9 and k = 0.8.
Fig. 8. (color online) The density–time curve of the 5th, 25th, 55th and 80th lattices with λ = 0.9 and k = 0.9.

Figures 912 show the relationship of the density difference and the densities at the 5th, 25th, 55th and 80th lattices. The Y axis indicates the density difference between time t and from t = 1 to t = 500, and the X axis represents the current density. Figures 912 are sets of dispersed points in the phase space plot which show the chaos of the traffic flow. The figures demonstrate that the densities of the 5th, 25th, 55th and 80th lattices will change with the feedback gain. Comparing patterns (a) with (b) of Figs. 912, we can clearly see that the amplitude of the density wave of every lattice approaches 0 under the influence of the positive feedback signal. To describe the stability of the system intuitively, local enlarged images are depicted in patterns (b) of the above figures.

Fig. 9. (color online) Scatter plot of the 5th lattice at time t = 1–500 s (a) without control signal (b) and with λ = 0.9, k = 0.85.
Fig. 10. (color online) Scatter plot of the 25th lattice at time t = 1–500 s (a) without control signal (b) and with λ = 0.9, k = 0.85.
Fig. 11. (color online) Scatter plot of the 55th lattice at time t = 1–500 s (a) without the control signal (b) and with λ = 0.9, k = 0.85.
Fig. 12. (color online) Scatter plot of the 80th lattice at time t = 1–500 s (a) without the control signal (b) and with λ = 0.9, k = 0.85.

In brief, the simulation manifests the feasibility of taking into account the double flux difference to alleviate traffic congestion, and the control signal also plays a vital role in the lattice hydrodynamic model. Thus, it can be concluded that the proposed model is useful for suppressing the increasingly serious traffic jams.

5. Conclusion

A new feedback signal considering the double flux difference and the new optimal velocity function is proposed based on the control method for the lattice hydrodynamic model. The stability condition of the double flux difference model is derived by the control theory, and the transform function and the judgement of its norm expansion are deduced in the paper. The control coefficients k and λ can be computed by the numerical simulation of flux curve and density curve. In addition, the analysis of chaotic dispersed points validates the consistency between the theoretical analysis and the simulations. The result of simulations show that the double flux difference plays a vital role in alleviating traffic congestion. Therefore, the feedback control is helpful to suppress traffic jams and this factor should be considered in the traffic flow model. In reality, drivers adjust their driving behavior with a time delay. So research on the effect of the double flux difference model considering delay-time feedback is our next task.

Reference
[1] Xu Y Jia L S Teng W Z Lu W Z 2015 Commun. Nonlinear. Sci. Numer. Simul. 22 285
[2] Jiang R chen J Y Ding Z J Ao D C et al. 2016 Trans. Res. Part C- Emering Tech. 69 527
[3] Tang T Q Huang H J Shang H Y 2017 Physica 468 322
[4] Peng G H Liu C Q Tuo M 2015 Physica 436 952
[5] Peng G H 2015 Mod. Phys. Lett. 29 1550174
[6] Li Z P Xu X Xu S Z 2017 Commun. Nonlinear Sci. Numer. Simul. 42 132
[7] Gupta A K Redhu P 2014 Nonlinear Dyn. 76 1001
[8] Tang T Q Huang H J Wong S C 2009 Chin. Phys. 18 975
[9] Tang T Q He J Yang S C 2014 Physica 413 583
[10] Li Z P Liu Y C 2006 Chin. Phys. 15 1570
[11] Ge H X Cheng R J Li Z P 2008 Physica 387 5239
[12] Zaki M Tarek S Khaled S 2014 Transport. Res. Rec. 2434 103
[13] Helbing D Tilch B 1998 Phys. Rev. 58 133
[14] Wilson R E Ward J A 2011 A Trans. Plan Tech. 34 3
[15] Yu L Li Tong Shi Z K 2010 Physica 389 2607
[16] Jiang R Wu Q S Zhu Z J 2001 Phys. Rev. 64 017101
[17] Zhou J Shi Z K Cao J L 2014 Nonlinear Dyn. 78 1765
[18] Kerner Boris S 2016 Physica 450 700
[19] Helbing D Hennecke A Shvetsov V Treiber M 2001 Transp. Res. 21 183
[20] Nagatani T 1996 J. Phys. Soc. Japan 65 3150
[21] Bando M Hasebe K Nakanishi K 1998 Phys. Rev. 58 5429
[22] Zeng Y Z Zhang N Liu L J 2014 Acta Phys. Sin. 63 068901 in Chinese
[23] Zhao X M Gao Z Y 2006 Physica 366 513
[24] Zhou T Sun D H Li H M 2014 Chin. Phys. 23 050203
[25] Zhao X M Xie D F 2013 Phys. lett. 377 3161
[26] Konishi K Kokame H Hirata K 1999 Phys. Rev. 60 4000
[27] Nagatani T 1998 Physica 261 599
[28] Zhu H B Dai S Q 2008 Physica 387 4367
[29] Ge H X Cui Y Zhu K Q Cheng R J 2015 Commun. Non. Sci. Numer. Simul. 22 903
[30] Cheng R J Han X L Lo S M 2014 Chin. Phys. 23 030507
[31] Redhu P Gupta A K 2015 Commun. Non. Sci. Numer. Simul. 27 263
[32] Li Y Zhang L Zheng T et al. 2015 Commun. Non. Sci. Numer. Simul. 29 224
[33] Zhu C Zhong S Li G et al. 2017 Physica 468 445
[34] Wang T Gao Z Y Zhao X M 2012 Chin. Phys. 21 020512