Abstract

A nonlinear car-following model with driver’s reaction time is studied from the synchronization transition viewpoint. We investigate the traffic congestion from the view of chaos system synchronization transition. Our result shows that the uniform flow corresponds to the complete synchronization and the stop-and-go congested state corresponds to the lag synchronization of the vehicles. An analytical criterion for synchronization manifolds stability is obtained; the analytical result and the numerical result are consistent. The synchronization transition is also trigged by the driver’s reaction time. We analyze the car-following model by the use of the nonlinear analysis method and derive the modified KdV equation describing the kink density wave.

1. Introduction

Traffic jams are very annoying in our life and have been studied by many physicists [1]; however, the precise mechanism for generation and propagation of traffic jams is still not clearly understood. Some traffic models have been established to study the dynamics of traffic flow, such as car-following models, cellular automaton models, gas kinetic models, and hydrodynamic models [2].

Car-following models or microscopic models describe the behaviors of individual vehicles, which are described by ordinary differential equation or delay differential equation. Recently, one more widely used car-following model is the optimal velocity model [3, 4]; by this kind of model, the effects of fluctuations of traffic jams are analyzed, the jamming transitions and density waves have been invested [59], and the bifurcation phenomena of the oscillating solution are also explored by using numerical continuations techniques [1012]. Another popular car-following model is the intelligent driver model that was introduced by Treiber et al. [13, 14]; in this model, all parameters have a clear physical meaning; congested traffic states have been empirically observed and microscopically simulated. Many models are able to explain uniform flow as well as stop-and-go waves; however, the transition [1528] between the two qualitatively different solutions is still not clarified.

According to the car-following model theory, for each individual vehicle, an equation of motion is the analogue of Newton’s equation for each individual particle in a system of interacting classical particles. In Newtonian mechanics, the acceleration may be regarded as the response of the particle to the stimulus, and it receives in the form of force which includes both the external force and those arising from its interaction with the other particles in the system. In traffic system, each driver can respond to the surrounding traffic conditions only by accelerating or decelerating the vehicle; so, the basic philosophy of the car-following theories can be summarized by the equation for the th vehicle ().

We know, for a chaotic system, the appearances and robustness of chaotic synchronization states have been established by means of different coupling schemes [29], one of which is Pecora and Carroll method (unidirectional coupling or drive-response coupling).

In this paper, we regard the car-following system as a drive-response coupling chaotic system; we make use of chaotic systems synchronization transition [3034] method to study the synchronization transition of microscopic movement of the vehicles and further reveal the relationship between synchronization transition and traffic congestion. First, we use the long wave expansion method to give an analytical criterion for synchronization manifolds stability. In order to verify the analytical result, we use DDE-BIFTOOL to perform a two-parameter bifurcation analysis of the model; due to the demands on the CPU time and the memory, the investigation was restricted to the setting of vehicles, and the analytical result is consistent with the numerical result. Second, we find different transition region in two-parameter plane and the vehicles display different motion. Third, we consider how the driver’s reaction time impacts on the synchronization transition. Finally, we analyze the car-following model by the use of the nonlinear analysis method and we derive the modified KdV equation describing the kink density wave.

The layout of this paper is organized as follows. We introduce the car-following model in Section 2 and derive analytical criteria for synchronization manifolds stability in Section 3. Numerical simulation is in Section 4. Nonlinear analysis is in Section 5. In Section 6, we present the main results.

2. Car-Following Model

Here, we consider a single lane of traffic with identical vehicles; displacements and velocities are denoted as and , respectively; the spacing of adjacent vehicle is called headway (). For the sake of simplicity, we suppose that vehicles are placed on a circular road of length (). The headway consists of the condition The acceleration of the th vehicle is given by The is called the sensitivity of the vehicles and the is driver’s reaction time. The function is optimal velocity function; it has the following properties.(a) is a nonnegative, continuous, and monotone increasing function.(b) as .(c)There exists a jam headway such that for .

The dimensional parameter OV function [12] is We introduce the rescaled variables , ; the OV function becomes Model (1) and model (2) are transformed into Here , ; in the remainder of this paper we take the rescaled OV function and model. In order to express simplicity, the OV function and model are expressed:

3. Synchronization Transition

The synchronization manifold of the system (7) and (8) is The system ((7), (8)) possesses uniform flow equilibrium All vehicles reach complete synchronization. The stability of the synchronization manifolds will change when the parameter is varied. To see whether the synchronization manifold of the system ((7), (8)) is stable or not, we add a small perturbation

According to (7) and (8), we can calculate a linearized equation with respect to the uniform flow (10). The linearized equation is

Indeed, we have , , and .

Orosz et al. have used the dynamical system approach and numerical continuation technique to give out the stable curves, but the numerical continuation method (DDE-BIFTOOL) needs a lot of CPU time and the memory; below, we can find a simple method to analyze the linear stability of the synchronization manifold [10].

Physical Approach. Let , , ; the characteristic equation is given by We know that the leading term of is of order . When , . We can derive the long wave expansion of , which is determined order by order around . By expanding the first- and second-order terms of are obtained: If is a positive value, the synchronization manifold is stable; if is a negative value, the synchronization manifold is unstable.

By , we get Equation (16) is the analytical criterion for synchronization manifold stability.

If the number of the vehicles is large, to obtain the stable boundary curve by the numerical continuations method requires a lot of time. In order to compare with numerical method, we now focus on the case of vehicles. For , the stable boundary curve is shown in Figure 1. We use the numerical continuation techniques (DDEBIFTOOL) to get the blue curve and we get the analytical curve (red curve) according to (16). We can find that the blue curve and the red curve completely overlap, and this shows that our analytical approach is correct. In Section 4 numerical simulations also verify our conclusion.

When the stability of the synchronization manifold is lost, we can get , , , and , and is the period. For two adjacent vehicles, the speed of the two vehicles has a certain time delay; the vehicles will get lag synchronization.

Figure 1(a) shows that the synchronization manifold is stable above the sable boundary curve and the vehicles can achieve complete synchronization; the synchronization manifold is unstable below the stable boundary curve and the vehicles can achieve lag synchronization. The synchronization transition of the vehicles arises when the sensitivity () is fixed and the headway () varies from the right to the left through the critical curve. Figure 1(b) shows that the unstable region is unbounded for .

Below we focus on the synchronization transition of the vehicles. If , according to Figure 2, by adding a vertical line across the top point () of the critical curve, we can divide the whole region into three regions, region I, region II, and region III. Region I and region II are stable; region III is unstable. If , the region is divided analogous to the situation when , as is shown in Figure 1(b).

4. Numerical Simulation

In Figure 3, we show the space-time evolutions of the vehicles’ velocity () and location () plotted against time (). We consider the velocity of the first and the second vehicles, the location of fifteen vehicles. First, we take a point in region III where , , , , and . The corresponding velocity-time plot and the space-time plot are shown in Figures 3(a) and 3(b). In region III, the headway is larger and the length of the road is longer; all vehicles travel with the same high speed. This is a freely moving phase. Second, we take a point in region II where , , , , and . The corresponding velocity-time history plot and the space-time plot are shown in Figures 3(c) and 3(d). The vehicles travel with stop-and-go state and get lag synchronization; the velocity of the adjacent vehicles has a certain time delay. Last, we take a point in region I where , , , , and . Figures 3(e) and 3(f) show the velocity-time history plot and the space-time relation; the vehicles travel with a low speed close to zero and get complete synchronization, that is, a uniformly congested phase.

In Figure 4 the stability under different time delay influence is compared. We know that the unstable region is bounded for ; when , the unstable region is unbounded. We find that the unstable region becomes larger as the time delay increases.

In Figure 5 we consider the time delay effect on the synchronization transition. For , the vehicles travel with the same speed and they get complete synchronization. For , the vehicles get lag synchronization and travel with moderate amplitude. If the time delay is increased to , the vehicles travel with large amplitude. If the time delay continues to increase, the traffic is completely congested.

According to the analytical stability condition (23), the stability of the synchronization manifold is strictly dependent on the three parameters .

5. Nonlinear Analysis

We analyze the car-following model by the use of the nonlinear analysis method. We derive the modified KdV equation describing the kink density wave. We now consider long-wavelength models in the traffic flow on coarse-grained scales. The simplest way to describe the long-wavelength models is the long wave expansion. We consider the slowly varying behavior at long wavelengths near the critical point . We assume that the value of is determined adiabatically by . This statement is expressed by the relation

We define the headway as and and are a scaled position and time defined by Then, from (18), is expressed as Since is a small parameter, a Taylor expansion can be applied to terms like in (20). This leads to the expression Here, remains unknown; we determine it by assuming the form From the form of (22), we find the terms are expressed as where , and so forth are constants, which are calculated in the following way. First, we substitute (22) (using (23)) into (8), where is expressed as a function of , According to (18), we get the expression about By expanding (25) to the fifth order of , we get Let where We expand the optimal velocity function at the turning point:

Then where Then, collecting terms of equal order on both sides and comparing coefficients of , and so forth, we obtain : Substituting derived in this way into (7), we get a reduced equation for : By taking , , the second-order term of and the third-order term of are eliminated from (33); then we get where Let if we ignore the terms in (34), then (34) is transformed into a Kdv equation In order to derive the regularized equation, we make the following transformation for (37): then (37) is transformed into a regularized Kdv equation If we ignore the terms in (39), this is just the modified Kdv equation with a kink solution as the desired solution The selected value of propagation velocity for the kink solution is determined from the term.

Next, assuming that , we take into account the correction: where In order to determine the selected value of the propagation velocity for the kink solution (37), it is necessary to satisfy the solvability condition where

Then, we can get

6. Conclusion and Discussion

In this paper, we study a car-following model with driver’s reaction time. First, we investigate the traffic congestion from the view of chaos system synchronization transition. Our result shows that the uniform flow corresponds to the complete synchronization and the stop-and-go congested state corresponds to the lag synchronization of the vehicles. Second, we derive analytical criteria for synchronization manifolds stability; the analytical result and the numerical result are consistent. Third, the synchronization manifold stability regions are further classified; we find that the vehicles exhibit different states of motion in different regions, and the synchronization transition can also be trigged by the driver’s reaction time. Last, we analyze the car-following model by the use of the nonlinear analysis method and we derive the modified KdV equation describing the kink density wave. Our result also shows that it is crucial to pay attention to the behavior of individual drivers in order to understand the emergent behavior of traffic. The driver’s reaction time is influenced by individual differences, personality, environment, and other factors. In further work, we need to study how to better characterize the reaction time of the drivers.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This research was supported by the National Natural Science Foundation of China (Grant no. 11032009) and National Natural Science Foundation of China (Grant no. 11272236).