Abstract

In order to prevent the occurrence of traffic accidents, drivers always focus on the running conditions of the preceding and rear vehicles to change their driving behavior. By taking into the “backward-looking” effect and the driver’s anticipation effect of flux difference consideration at the same time, a novel two-lane lattice hydrodynamic model is proposed to reveal driving characteristics. The corresponding stability conditions are derived through a linear stability analysis. Then, the nonlinear theory is also applied to derive the equation describing traffic congestion near the critical point. Linear and nonlinear analyses of the proposed model show that how the “backward-looking” effect and the driver’s anticipation behavior comprehensively affect the traffic flow stability. The results show that the positive constant , the driver’s anticipation time , and the sensitivity coefficient play significant roles in the improvement of traffic flow stability and the alleviation of the traffic congestion. Furthermore, the effectiveness of linear stability analysis and nonlinear analysis results is demonstrated by numerical simulations.

1. Introduction

With the rapid development of wireless technologies and intelligent vehicles in the past years, the vehicle-to-vehicle (V2V) communication [1, 2] can not only alleviate the congested traffic flow but also provide pleasurable driving experiences. To reveal the realistic traffic phenomena, there have emerged various traffic models such as car-following models [310], cellular automation models [11, 12], macrotraffic models [1315], lattice hydrodynamic models [1621], continuum models [2224], and gas kinetic models [25].

By incorporating the ideas of car-following as well as continuum models, the first lattice hydrodynamic model [26] with a conservation equation (1) and a flow equation (2) was proposed in 1998.where and denote the optimal flow and the actual flow, respectively. is the average density, and refers to the driver’s sensitivity. Then, the stability condition and equation can be, respectively, obtained with the linear analysis and nonlinear analysis methods. Subsequently, many extended works have been developed by taking different factors into accounts, such as driver’s memory [2731], driver’s anticipation effect [3234], density difference [35], traffic interruption probability [3638], and “backward-looking” effect [3941].

In addition to the information from the preceding vehicle, drivers always focus on the running conditions of the rear vehicles from the rear-view mirror to change their driving behavior. By considering the “backward-looking” effect, Ge and Cheng [39] presented the BL-LV model as follows:where and represent the Heaviside function and the safety density, respectively. denotes a parameter close to 1. is the optimal velocity function for “forward-looking.” While is the optimal velocity function for “backward-looking,” which works if and only if holds. Functions and have exactly the opposite effects and are defined aswhere is a positive constant that represents the relative role of the “backward-looking.”

For the sake of safe and comfortable driving, drivers always predict the velocity of the preceding vehicles so that they can adjust their driving velocity at the next moment. This driving behavior is known as the anticipation effect [42] or predictive effect [43, 44]. Wang et al. [43] added the driver’s anticipation effect into the optimal velocity (OV) function in a single-lane lattice model. Later, Kaur and Sharma [44] extended the single-lane model [43] to the two-lane case along with the optimal current difference (OCD) effect. By considering the anticipation effects in the OCD term as well as the OV function, Sharma [45] constructed a new lattice hydrodynamic model with an anticipation coefficient , where the positive represents the anticipation effect while the negative represents the delay effect. Then, Li et al. [32] adopted this anticipation effect to the multiple optimal current differences’ anticipation version. Meanwhile, numerical results indicated that only the information of three preceding lattices [46] was enough for the improvement of traffic stability. In view of the anticipation effect of flux difference, Chang and Cheng [47] proposed a new lattice model as follows:where represents the driver’s anticipation time and and denote the deviation degree of speed and the feedback gain, respectively. Some previous studies have considered the “backward-looking” effect and the anticipation effect of flux difference separately. In a real traffic system, the velocity of the current vehicles is not only affected by the preceding vehicles but also affected by the rear vehicles. Generally, the information of the rear vehicles could be received by the rear-view mirror. According to the information, drivers could change their driving behavior to avoid the traffic accidents. Meanwhile, the velocity of the preceding vehicles could be predicted by drivers on the basis of the vehicle distance information, so that they can adjust their driving velocity in a period of anticipation time. However, the above two factors have never been considered simultaneously in the existing lattice models. This motivates us to develop a new lattice model by incorporating the “backward-looking” effect with the anticipation effect of flux difference.

The outline of this paper is organized as follows. In Section 2, we introduce a new lattice hydrodynamic model considering the flux anticipation difference effect and the “backward-looking” effect and use the linear stability theory to derive the stability condition of the traffic flow. In Section 3, the propagation behavior of traffic density waves is presented by means of the nonlinear analysis method. In Section 4, numerical simulations are carried out to validate the results of linear and nonlinear stability analyses. At last, the conclusions are drawn in Section 5.

2. The Extended Lattice Hydrodynamic Model and Linear Stability Analysis

In light of the aforementioned facts, a novel lattice model is presented by taking the flux anticipation difference effect and the “backward-looking” effect (FADBE) into account as follows:where is the sensitivity coefficient of the anticipation flux difference. The greater is, the stronger the driver’s anticipation effect is. means the anticipation time, and represents the anticipation flux difference between the preceding lattice and the current lattice . The optimal velocity functions and the function are the same as those in Ge and Cheng’s model [39].

After eliminating the velocity term in equations (1) and (6), the following density equation can be derived:

When , the model reduces to Ge and Cheng’s model [39]. In addition, the lattice model coincides with Nagatani’s model [26] when and .

The influence of the “backward-looking” effect and the anticipation effect has been studied by linear stability analysis. The uniform traffic stream with the fixed density and the optimal velocity is defined as the stable state. One can obtain the solution of the steady-state vehicular system as follows:

Then, is supposed as a small perturbation added into the steady-state solution at lattice , and the perturbed solution is

Equation (9) is substituted into equation (7), and the derived equation is linearized. The equation is obtained as

Let the small perturbation , and inserting it into equation (10), we gain

Assuming and putting it into equation (11), we reserve the first-order and second-order terms of by neglecting higher-order terms. Then, one can get

The uniform steady-state flow remains stable as is positive. On the contrary, the uniform flow tends to be unstable if the value of is negative. When is equal to zero, the formula of the neutral stability condition is acquired aswhere . The stability condition for the uniform traffic flow is obtained as

It is obvious that the stable condition (14) is reduced to that of Ge and Cheng’s model [39] when . Equation (14) shows that the parameters , , and play important roles in the performance of the traffic system.

The phase diagram of is shown in Figure 1, where is density and denotes sensitivity. Note that denotes a positive constant considering the “backward-looking” effect [40, 41]. Three solid curves in each pattern of Figure 1 represent the neutral stability curves under . The phase diagram is divided into two regions: the stable and unstable regions are above and below the solid neutral stability curve, respectively. In the stable region, the traffic flow will remain stable with a perturbation. In the unstable region, the congested traffic will happen even with a small perturbation. Figure 1 depicts the neutral stability curves of the FADBE model for different values of driver’s anticipation time and the sensitivity coefficient of the anticipation flux difference when . As shown in Figure 1, the corresponding stable region expands gradually with the increasing values of and , which means that the driver’s anticipation effect can contribute to the stabilization of traffic flow obviously.

Figure 2 demonstrates the neutral stability curves of the FADBE model for different values of and when . As shown in Figure 2, the corresponding neutral stability curves and critical points decrease gradually with the increasing values of and . When the value of is fixed in each pattern of Figure 2, the stability of the traffic flow can be improved efficiently with the increasing value of anticipation time .

Figure 3 exhibits that the phase diagram of of the FADBE model for different values of the positive constant and the sensitivity coefficient when . With the increasing value of the positive constant in each pattern of Figure 3, the amplitudes of neutral stability curves all fall down gradually when the value of the parameter is fixed. For different values of , it can be seen from each pattern of Figure 3 that the corresponding peak value of solid curves lowers down with the increasing value of . These phenomena demonstrate that when the value of or increases, the area of the stable region becomes larger. Thus, it can be concluded that the traffic flow stability can be improved by incorporating the positive constant and the sensitivity coefficient .

3. Nonlinear Stability Analysis

To investigate the influence of the “backward-looking” effect and the anticipation effect of flux difference, we carried out the reduction perturbation method to obtain the mKdV equation. The slow variables and are defined as follows:where is a small positive scaling parameter.

Substituting equation (15) into equation (7), each term of equation (7) is expanded up to the fifth-order of with the Taylor expansion technique as follows:where and . Near the critical point , the value of is defined as . The squared and cubic terms of are eliminated in equation (16) with the consideration of . One can get the simplified equation as follows:where

To derive the standard equation, the following transformations are applied:

Thus, equation (17) is converted into the standard equation as follows:

When the term is ignored, the kink-antikink soliton solution of the equation can be rewritten aswhere represents the propagation velocity. The specific value of can be provided only when the following solvability condition is satisfied:where . By solving equation (22), the general solution of propagation velocity can be calculated as

Hence, the corresponding kink-antikink solution is derived aswhere , and the amplitude of the density wave is defined as . According to the kink-antikink solution equation (24), one can get the coexisting phases which include the freely moving phase with low density and the congested phase with high density .

Based on the linear and nonlinear theory analyses, numerical simulations are carried out with the comprehensive consideration of the anticipation effect of flux difference and the “backward-looking” effect.

4. Numerical Simulations

As a positive constant considering the “backward-looking” effect has been studied fully in some previous studies [36, 40, 41], this paper will focus on the effects of parameters and in this part. Periodic boundary conditions are adopted by a step function as follows:where is the initial perturbation and other parameters are set as , , and .

4.1. The Effect of Anticipation Time

In order to analyze the anticipation time’s influence on the stability of the traffic system, the temporal and spatial evolution of density waves between time is depicted in Figure 4 for different anticipation time with fixed parameters and . Since the given parameters in Figures 4(a)4(c) do not satisfy the linear stable condition equation (14), kink-antikink density waves appear and the initial stable traffic flow evolves into nonuniform traffic flow with the added perturbation. To reveal the amplitude of the density wave more clearly, Figure 5 describes the density profile of each subfigure shown in Figure 4 at t = 10200. Moreover, Figures 5(a)5(c) reveal that the amplitudes of density waves corresponding to Figures 4(a)4(c) weaken as the anticipation time increases. When , i.e., the parameters in Figure 4(d) satisfy the linear stability condition, the added perturbation will disappear, and the traffic flow will eventually return to a uniform state. Figure 5(d) can get results similar to Figure 4(d). This phenomenon indicates that driver’s anticipation effect of flux difference could significantly promote the stabilization of traffic flow.

4.2. The Effect of Sensitivity Coefficient

To expose the sensitivity coefficient’s influence on the stability of the traffic system, the temporal and spatial evolution of density waves between time is exhibited in Figure 6 for different values of with fixed and . When , the FADBE model degrades into Ge and Cheng’s model [39]. Note that the parameters in Figures 6(a)6(c) do not satisfy the linear stable condition equation (14). Therefore, the stop-go traffic density waves occur with the small added perturbation in Figures 6(a)6(c). Meanwhile, the amplitude of the density wave decreases gradually as the value of increases. Especially in Figure 6(d), the congested state entirely disappears when . The relationship between Figures 6 and 7 is similar to the relationship between Figures 4 and 5. Figure 7 describes the density profile of each subfigure shown in Figure 6 at t = 10200. As parameter increases, the traffic flow becomes more gentle and stable. This indicates that the sensitivity coefficient contributes to suppressing the traffic congestion.

Figures 8(a) and 8(b) describe the flux against the density with different parameters corresponding to the traffic flow in Figures 4 and 6, respectively. In Figure 8, the amplitude of the hysteresis loops reduced with the increasing value of or . When the corresponding parameter or , the hysteresis loop turns into a stable point, which means that the traffic flow eventually becomes stable even with the added small perturbation. To sum up, as the anticipation time or the sensitivity coefficient increases, the stability of the system can be enhanced by incorporating the driver’s anticipation effect and the “backward-looking” effect.

5. Conclusion

In this paper, a novel lattice hydrodynamic model with the anticipation effect of flux difference and the “backward-looking” effect is constructed to explore the traffic driving behavior for a single-lane traffic system. The stability conditions and the solution of the equation are deduced, respectively, based on the linear stability analysis and nonlinear theory. The results of the theoretical analysis display that the increasing values of both the driver’s anticipation time and the sensitivity coefficient of the anticipation flux difference can contribute to enlarging the stable region. Numerical simulations reveal that the FADBE model is effective in stabilizing traffic flow even with perturbations. By considering the anticipation effect of flux difference and the “backward-looking” effect simultaneously, the traffic congestion can be alleviated effectively, which is in accordance with the conclusions of the theoretical analysis. Future work will focus on solving the traffic congestion on multilane road problem, the automation for longitudinal driving, and the vehicular emergency warning system.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

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

Acknowledgments

This work was supported in part by the National Natural Science Foundation of China (Grant no. 61903203) and the student’s innovation and entrepreneurship training program from Jiangsu Provincial Department of Education (Grant no. 202110304055Y).