Linear Stability and Nonlinear Analysis of an Extended Optimal Velocity Model Considering the Speed Limit ()
1. Introduction
With the development of social economy and the increasing number of motor vehicles, traffic jam has become the common bottleneck of urban development. The problem of traffic congestion can be attributed to the stability and solitary waves of traffic flow models. Generally speaking, traffic flow models have been divided into macroscopic models, microscopic models, and mesoscopic models according to the aggregation level. There are many microscopic models in traffic flow, and the most commonly studied model is the car-following model. The car-following model is a microscopic traffic flow model used to describe the behavior of a single driver. In order to study the complex characteristics of traffic flow, scholars have proposed various car-following models to understand the formation and transmission mechanism of traffic congestion.
The early research work on the car-following theory was originated from Reuschel [1] and Pipes [2]. Since then, scholars have carried out many related research works. In 1961, Newell [3] proposed a car-following model of velocity-governing equation by assuming that the velocity adjustment of the rear car depends on the optimized velocity of the headway. In 1995, Bando et al. [4] proposed the optimal velocity model (OVM), in which the optimal velocity was determined only by the headway with the forward vehicle. Despite the simplicity and few parameters, the OVM can be used to describe various characteristics of the actual traffic flow, such as stop-and-go phenomenon, the instability of traffic flow and the evolution of traffic jams. In 1998, Helbing and Tilch [5] verified OVM by using empirical data and proposed a generalized force model (GFM), which overcomes the problems of unrealistic deceleration and excessive acceleration of OVM. In 2001, Jiang et al. [6] found that the GFM could not be used well to describe the delay time and the kinematic wave speed at jam density and proposed a full velocity difference model (FVDM) by considering the influence of positive and negative speed differences on the car-following models. In 2008, in order to overcome the high deceleration of FVDM, Ge et al. [7] proposed a two-velocity difference model (TVDM) based on the application of intelligent transportation system (ITS). Recently, some new car-following models were submitted to describe traffic nature more realistically [8] [9] [10] [11] [12]. Undoubtedly, these scholars have made various contributions to the establishment and development of traffic flow. For the above car-following models, there have been many researches focusing on the stability and density waves in the past time.
In the past, more research has focused on the density waves for the car-following models. However, it is still enormous significance to derive thermodynamic theory of traffic flow model and to link traffic meta-stability with spinodal decomposition in first-order phase transition [13]. The TDGL equation can refer to non-equilibrium phase transition [14]. Nagatani [15] firstly proposed the thermodynamic theory of traffic flow, and derived the TDGL equation of two simple lattice hydrodynamic models near the critical point and the thermodynamic potential of traffic flow by applying the reduced perturbation method. So far, there has been little research on deriving TDGL equation from car-following models [16] [17] [18] [19]. Since the TDGL equation of the traffic flow model can be used to represent the thermodynamics theory of jamming transition, it is very important to establish the traffic flow model and derive the corresponding TDGL equation.
In recent years, traffic jams are still intense as the number of cars in big cities increases and road infrastructure becomes more complex. Exploitation of intelligent transport system is a part of the best ways for settling transport problems. Considering the influence of an intelligent transportation system, many researchers have established some car-following models to describe the real nature of traffic flow [20] [21] [22] [23] [24]. In fact, intelligent vehicle navigation system can provide service for the vehicles, when the actual speed of the vehicle deviates from the speed limit. If the vehicle overspeeds, it is easy to cause traffic congestion and traffic accidents, and the driver can learn the speed limit information in advance, he can adjust in advance to achieve a stable speed. At this time, the traffic can develop smoothly. In other words, drivers’ knowledge of speed limit information in advance plays an important role in stabilizing the traffic flow. However, most of the car-following models have not considered the effect of the speed limit.
In order to reveal the effect of the speed limit on traffic flow, an extended optimal velocity model is introduced by considering the speed limit in this paper. The content of this paper is organized as follows. In Section 2, an extended optimal velocity model with consideration of the speed limit is introduced. In Section 3, the linear stability condition is obtained from the linear stability analysis method. In Section 4, the TDGL equation and its corresponding soliton solution are obtained from the reductive perturbation method. In Section 5, the mKdV equation is derived in the unstable region. In Section 6, conclusions are given.
2. An Extended OVM with Speed Limit
In 1995, Bando et al. [4] proposed the OVM based on the analysis of the characteristics of traffic flow. The equation of motion is as follows:
(1)
where a is the driver’s sensitivity coefficient,
is the position of the vehicle at time t,
is the optimal velocity function, and
represents the space headway between the car ahead and the car following at time t.
As people make full use of the latest information of the intelligent transportation system, some traffic problems are gradually solved. For example, intelligent vehicle navigation system uses computer and communication technology to provide information relevant to the cars on the road and generate dynamic traffic information, such as the speed limit [25]. It can determine at what appropriate velocity the driver should drive under the current conditions according to the current traffic conditions and signal status. When the actual velocity of the vehicle deviates from the speed limit by the navigation system, the driver will adjust the actual velocity according to the stable speed limit to improve the traffic capacity and relieve traffic congestion. Therefore, we propose the following extended optimal velocity model by considering the speed limit:
(2)
where
is the speed limit provided by the navigation system, and
denotes the reactive coefficient related to the speed limit
.
We chose the optimal velocity function
as proposing in reference [4].
(3)
where
is the maximum velocity of the vehicles,
is the safety distance of the vehicles. The optimal velocity function
is a function of
, and it’s a monotonically increasing function with an upper bound.
For the convenience of analysis in the following, we discretize Equation (2) with asymmetric forward difference [24], and rewrite it in terms of headway as follows:
(4)
where
is time step and
.
3. Linear Stability Analysis
In this section, we study the extended optimal velocity model of Equation (4) by the method of linear stability analysis.
First, we assume that the initial state is stable, the headway of the vehicle is h and the corresponding optimal velocity is
. At this time, the vehicle position of the steady-state traffic flow can be expressed as
(5)
where N is the number of cars, L is the road length, and h is the average headway. Let
be a small deviation from the steady-state solution
, then the perturbed solution is given as
(6)
Substituting (6) into Equation (4), we obtain the following formulation:
(7)
where
and
.
Let
, Equation (7) can be rewritten as follows:
(8)
where
. Since
as
, z can be expressed by a long wave as
.
Substituting it into (8) and neglecting the higher order items, we obtain the coefficients of the first-and second-order term of
, as follows:
(9)
Clearly, if
is negative, the initial steady uniform flow will become unstable; if
is positive, the original steady flow state remains unchanged. Thus, the neutral stability condition of an extended optimal velocity model considering the speed limit is obtained as follows:
(10)
We are concerned with
relating to
, the reactive coefficient relating to
. Thus, the uniform traffic flow remains stable if the following condition holds:
(11)
As
, we have
. This shows that the stable region of (11) is larger than that of (12) the OVM (where
) given as follows:
(12)
From (11), we can see that the parameter
has an important effect on the stabilization of traffic flow. Notice that
,
, at this time,
is increases monotonically with respect to
, so it holds that
,
. As a result, the stability condition (11) is weaker than (12), and the stable region is larger than OVM.
Figure 1 shows the phase diagram in the
-phase where h (meter) is the headway and a (1/second) is sensitivity which corresponds to the inverse of the delay time. It shows the stable neutral lines of (10) with different values of
, where
is taken as the maximal velocity and
is the safe distance. In Figure 1, above the neutral stability curve is a stable region, which represents a free phase without traffic congestion, and below the neutral stability curve is an unstable region, which represents the traffic jam phase that evolves backward as the density waves stop-and-go. The peak of each curve represents the critical point
. Figure 1 also shows that as the sensitivity coefficient
to the speed limit increases, the stable region also gradually expands. This shows that if the driver can obtain the speed limit information through the intelligent vehicle navigation system in advance, he can adjust the vehicle speed earlier, effectively enhance the stability of the traffic flow, and avoid traffic congestion.
4. TDGL Equation
Now we use the long-wavelength modes to derive the TDGL equation for describing the pedestrian flow on a coarse-grained scale [26]. The long-wavelength
Figure 1. Phase diagram in headway-sensitivity space
.
expansion is the simplest way to describe the behavior of the long-wavelength models.
First, we analyze the slow-varying behavior of long waves near critical point
. Introduce the slow scales
of space variable n and time variable t [27], and define the slow variables X and T as follows:
(13)
where b is a constant. The headway
is defined by
(14)
Next, by expanding Equation (4) to the fifth-order of
with the use of (13) and (14). We obtain the expression:
(15)
Here, the coefficients
are given in Table 1. Where
, and
.
Now, we study the traffic flow near critical point
. By taking
, we eliminate the second- and third-order terms of
from Equation (15). The simplified equation of Equation (15) is as follows:
(16)
Here, the coefficients
are given in Table 2.
By transforming variables X and T into variables
and
, and taking
, Equation (16) can be rewritten as follows:
Table 1. The coefficients
of (15).
Table 2. The coefficients
of (16).
(17)
By adding term
to both of the left-and right-hand sides of (17) and transforming
and
in (17), we obtain
(18)
We define the thermodynamic potentials:
(19)
By taking (19) into Equation (18), the TDGL equation becomes
(20)
with
(21)
where
is given by Equation (19), and
is the derivative of function. The TDGL Equation (20) has two steady-state solutions except for a trivial solution
. One is the uniform solution
(22)
and the other is the kink solution
(23)
where
is a constant. Equation (23) represents the coexisting phase.
By the condition
(24)
Hence, substituting Equation (19) into Equation (24), we can obtain the coexisting curve related to the original parameters
(25)
The spinodal line is given by the following condition
(26)
From Equation (18), we obtain the spinodal line described by the following equation
(27)
The critical point is given by the condition
(28)
Substituting Equation (19) into Equation (28), we have the critical point related to the original parameters as follows:
(29)
5. mKdV Equation and Its Connection to TDGL
In this section, based on the stability condition (11) in Section 3, we use the reductive perturbation method to derive the mKdV equation for the model (4) in the unstable region of traffic flow. Similar to the derivation of the TDGL equation in Section 4, we study the slowly varying behavior at long wavelengths near the critical point. We also give slow scales for space variable n and time variable t.
By inserting
and
into Equation (15), making
, and near the critical point, we obtain:
(30)
Here, the coefficients
are given in Table 3.
In Table 3,
,
. We make the following transformation to Equation (30):
(31)
so the standard mKdV equation as follows:
(32)
If we ignore the term of
, Equation (32) is the modified KdV equation. Its kink solution is given as
(33)
Now, assuming that
, we take into account the
correction. For the purpose of determining the selected value of the velocity c for the kink solution, it is necessary to satisfy the solvability condition
as
, where
. We get the general velocity c [28],
(34)
Therefore, the general kink-antikink soliton solution is obtained as follows:
(35)
Table 3. The coefficients
of (30).
Since the kink soliton solution represents the coexisting phase in the space
, and the kink solution (35) is agreed with the solution (23) obtained from the TDGL Equation (20). We see that the jamming transition can be described by both the TDGL equation with a nontravelling solution and the mKdV equation with a propagating solution [29].
6. Conclusion
We propose an extended optimal velocity model (2) by considering the speed limit with intelligent prompts. Using the linear stability theory, the neutral stability line and the critical point of the new model are derived. The stable condition (11) shows that the effect of intelligent prompt limiting speed has a positive effect on expanding the stable area of traffic flow and easing traffic congestion. In addition, the TDGL equation is derived to describe the traffic behavior near the model critical point by using the reduced perturbation method, and the corresponding two steady-state solutions are obtained. From the TDGL, the spinor line and the critical point equation are calculated. At the same time, we derive the mKdV equation in the unstable region and obtain the relation between TDGL and mKdV equation.
Acknowledgements
We thank the Editor and the referee for their comments. This study was supported by the National Natural Science Foundation of China (Grant No.11162020).