Abstract
A three-lane traffic flow model is proposed to investigate the effect of car accidents on the traffic flow. The model is an extension of the full velocity difference (FVD) model by taking into account the lane changing. The extended lane-changing rules are presented to model the lane-changing behaviour. The cases that the car accidents occupy the exterior or interior lane, the medium lane, and two lanes are studied by numerical simulations. The time-space diagrams and the current diagrams are presented, and the traffic jams are investigated. The results show that the car accident has a different effect on the traffic flow when it occupies different lanes. The car accidents have a more serious effect on the whole road when they occupy two lanes. The larger the density is, the greater the influence on the traffic flow becomes.
1. Introduction
Traffic jam is a complicated traffic phenomenon and has been investigated by various traffic models [1–18]. Traffic accidents are one of the important sources of traffic jams. The accident usually leads to vehicle delay and reduces the capacity of the roads.
Recently, the studies of the influence of the car accident on traffic flow have attracted considerable attention. Nagatani [19] applied a cellular automaton (CA) model to study the effect of a traffic accident on the dynamical jamming. Zhu et al. [20] investigated the physics of congested traffic pattern induced by a car accident using a two-lane CA model. Kerner et al. [21] proposed a simple two-lane CA model to explain the physics of traffic breakdown, highway capacity, and synchronized flow. Qian et al. [22] presented an improved CA model to evaluate the effect of an accident. Sheng et al. [23] proposed a new CA model to investigate the temporary bottleneck by a special accident.
Kurata and Nagatani [24] extended the optimal velocity (OV) model [25] to two-lane traffic by introducing the lane-changing rules and used this model to study traffic jams induced by a car accident. Sugiyama and Nagatani [26] presented an extended OV model to simulate the multiple-vehicle collision when a vehicle decelerates suddenly in single-lane traffic flow. Based on the full velocity difference (FVD) model [27], Tang et al. [28] suggested a new car-following model with consideration of the effect of the traffic interruption probability.
As for the study of the macroscopic model, Tang et al. [29] proposed a new macromodel taking into account the traffic interruption probability. This model was derived from the microscopic model [28] by using the relationship between the microscopic and macroscopic variables. Zhao [30] adopted this macromodel to explore the effect of the accident position on shock and rarefaction. Tian and Sun [31] extended this model to two-lane traffic flow. To explore the effects of static bottleneck on traffic flow, Tang et al. [32] proposed a new macromodel. Tang et al. [33] further developed an extended macromodel considering effects of multistatic bottlenecks on traffic flow.
The abovementioned models concentrate on modeling the effects of car accidents on traffic flow at a single-lane or two-lane roadway. However, there exist a large number of three-lane roads in real traffic. In comparison with two-lane roads, the greater traffic flow and more frequent lane changing in three-lane roads make the accident easier to happen. It is thus necessary to further investigate the effect of the car accidents on traffic flow in three-lane roads.
In this paper, we extend the FVD model to the three-lane traffic by taking into account the lane changing. We investigate the traffic jams caused by the car accidents which occupy the exterior or interior lane, the medium lane, and two lanes. We derive the time-space diagrams and the current diagrams to study the effect of the car accidents on traffic flow (Figure 1).

(a)

(b)

(c)
2. Model
We consider traffic flow on a highway which consists of three lanes: the exterior lane (lane 1), the medium lane (lane 2), and the interior lane (lane 3). We discuss the cases where a car accident is located at lane 1 or lane 2, two car accidents occupy lane 1 and lane 2, and two car accidents occupy lane 1 and lane 3, as shown in Figure 2, where the car accidents are marked as black. Cars can overtake the car accident by changing lanes if the lane-changing criteria are fulfilled.

2.1. Basic Model
The car motion is divided into two parts: the forward motion and the sideways motion. We employ the FVD model [27] to describe the forward motion. This model has the form where is the position of car at time , is the headway of car , denotes the optimal velocity, is the sensitivity, represents the velocity difference between car and its adjacent front car , and is the sensitivity coefficient of a driver to the velocity difference.
The optimal velocity function can take many forms. We adopt the form suggested by Bando et al. [25]: where is the maximal velocity and is the safety distance.
2.2. Lane-Changing Rules
In general, a driver has an incentive to change lane if the adjacent lane can provide better driving conditions. Lane-changing decision of a driver needs to satisfy two criteria: incentive criterion (i.e., the adjacent lane has a better driving condition) and the security criterion (i.e., a driver must ensure the safety in the process of lane changing). As an example, Figure 2 shows a schematic illustration of lane changing for car on lane 2. In Figure 2, , and denote, respectively, the distances between car and its adjacent front car in lane 1, lane 2, and lane 3, while , and represent, respectively, the distances between car and its following car in lane 1, lane 2, and lane 3. In this paper, we adopt the symmetric lane-changing rules. That is, cars on lane 2 can change to lane 1 or lane 3 while cars on lane 1 and lane 3 can also change to lane 2.
Based on the lane-changing rules [24, 34], we propose the following extended lane-changing rules at time .(a)cars on lane 1 change to lane 2.
We adopt the rules of Kurata and Nagatani [24]: or or (b)Change from lane 3 to lane 2, if or or (c)Change from lane 2 to lane 1, if or or or or or or (d)Cars on lane 2 change to lane 3.
The rules are similar to those of cars on lane 2 changing to lane 1 (consult (9)–(15)).(e)Change from lane 2 to lane 1 and lane 3 with equal probability, if or or
In (3)–(18), , and denote the velocity of car on lane 1, lane 2, and lane 3, respectively.
The explanation of lane-changing rules is as follows. The rule described by (3) means that a driver has a motivation to change lane when the velocity of his car is 1.02 times larger than that of the front car and the distance between them is less than , and his car can change lane successfully if the distance between his car and its adjacent front car in lane 2 (target lane) is larger than and the distance between his car and its following car in lane 2 is larger than . In comparison with (3), in the rule given by (4) the incentive criterion remains, and the security criterion is that the velocity of adjacent front car in lane 2 (target lane) is larger than that of car , and the distance between car and its adjacent front car and the distance between car and its following car in lane 2 are both larger than . The rule given by (5) means that a driver has an incentive to change lane when the headway is less than , and his car can change lane successfully if the distance between his car and its adjacent front car in lane 2 (target lane) is larger than his headway and the distance between his car and its following car in lane 2 is larger than . The rules given by (6)–(8) are similar to those described by (3)–(5). The rules given by (9)-(10) are an extension of two-lane rule described by (3) for three-lane. A car in lane 2 can change to lane 1 because two criteria of lane changing are fulfilled, but it is not able to change to lane 3 because the security criterion or is not satisfied. Similarly, the rules given by (11)–(13) are an extension of the rule described by (4), and the rules given by (14)-(15) are an extension of the rule described by (5). The rules given by (16)–(18) mean that when a car in lane 2 can change to either lane 1 or lane 3, then it changes two lanes with equal probability.
3. Simulation and Results
In this section, we use the above model to explore the effects of car accidents on traffic flow. There are a certain number of cars in each lane at the beginning of simulation. The cars’ motion is split into two parts: one is the forward motion and the other is the lane changing. The forward motion is described by (1). The judging criteria of the lane changing are determined by (3)–(18). In the simulation, model (1) is solved by using the Euler method, where the time interval is , the sensitivity , the maximal velocity , the safety distance , and the sensitivity coefficient .
The simulation is performed by varying the initial density (or initial headway ) for road length . All cars in each lane are numbered in order from 1 to , where is the total number of cars on the lane. We adopt the open boundary condition for each lane. At each time of interval , we check whether the position of the leading car is longer than the length of the road , and if so we will remove the leading car and set the car behind as a new leading car. If the distance between the last car and the entrance is greater than the initial headway , a new car is added as the last car.
3.1. The Car Accident in the Exterior or Interior Lane
We first investigate the case that the accident happens on the exterior lane (i.e., lane 1). The situation that the car accident is in the interior lane (i.e., lane 3) can be discussed similarly. Assume that a car accident is located at in lane 1. Figures 3(a)–3(c), respectively, show the time-space diagrams of lane 1, lane 2, and lane 3 at a low density, . The car accident has a little effect on the traffic flow. This is because lane 2 and lane 3 can ease the effect of the car accident at this initial low density. Not all cars in lane 1 which pass through the car accident change lane back to lane 1 since the traffic conditions of lane 2 and lane 3 are also good. Figures 4(a)–4(c), respectively, show the time-space diagrams of lane 1, lane 2, and lane 3 at a high density, , using the same parameters as those in Figure 3. The jam appears just behind the car accident in lane 1. This is caused by the cars that are hindered on lane 1 and cannot change lane to the right in time. Corresponding to the local jam of lane 1, a region with an increase in density appears in lane 2. The result is due to the lane changing. The car behind the car accident in lane 1 has to change to lane 2. It will prevent the following car in lane 2 from getting ahead. Then the following car in lane 2 slows down and a localized high density region eventually is formed. When the traffic circumstance of lane 2 becomes worse, more cars will change to lane 3. The density of lane 3 increases.

(a)

(b)

(c)

(a)

(b)

(c)
Figure 5 shows the plots of traffic currents against the density. In Figure 5, the solid curve represents the theoretical current curve which is computed as , and other lines indicate the traffic currents on three lanes obtained by simulation. When the density is low, the current increases with density. Moreover, the currents on three lanes are similar to each other and nearly agree with the theoretical current curve. With the increase of density, the currents of lane 2 and lane 3 approach a saturated state while the current on lane 1 drops and then keeps almost constant value. When the density increases even further, the current on three lanes decreases with increasing density.

3.2. The Car Accident in the Middle Lane
For a comparison with the case that the car accident is in lane 1, we next study the case that a car accident occupies the middle lane. The car accident is located at in lane 2. Figures 6(a)–6(c) show the time-space diagrams of lane 1, lane 2, and lane 3 at a low density, , respectively. It can be seen that the effect of the car accident is eased by lane 2 and lane 3 at this initial low density. Figures 7(a)–7(c) show the time-space diagrams of lane 1, lane 2, and lane 3 at a high density, , respectively. A localized jam appears behind the car accident in lane 2. Comparing with the results of the car accident located at lane 1, we can find that the length of jam behind the car accident is shorter when the car accident is located at lane 2. This is because we adopt the symmetric lane-changing rules and the cars behind the car accident in lane 2 can directly change to lane 1 or lane 3.

(a)

(b)

(c)

(a)

(b)

(c)
Figure 8 shows the plots of traffic currents against density when the car accident is in lane 2. From Figure 8, it can be found that the current of each lane increases with increasing density when the density is low. The transition point from free to congestion state is larger than the case that the car accident occupies lane 1, because there are two lanes to ease the effect of the car accident at the same time. When the density is high, the currents of lane 1 and lane 3 approach a saturated state, while the current of lane 2 decreases and then keeps almost constant. As the density is higher, the current of each lane decreases with increasing density, and the currents of lane 1 and lane 3 are higher than lane 2. Moreover, we can see that the currents of lane 1 and lane 3 are similar to each other since the symmetric lane-changing rules are used.

3.3. The Car Accidents in Two Lanes
In real traffic, many traffic accidents happen when cars are changing lanes. Under this situation, the car accidents usually occupy two lanes at the same time. Thus, we finally explore the cases that there are two car accidents.
The first situation is that two car accidents are, respectively, located at in lane 1 and in lane 2. Figures 9(a)–9(c) illustrate the time-space diagrams of lane 1, lane 2, and lane 3 at a low density, , respectively. From Figure 9, we can see that the jam is formed behind the car accident in lane 1, and a region with a slight increase in density appears in lane 3 because all cars behind the accident car in lane 1 and lane 2 need to change to lane 3 and go through the accident car. Figures 10(a)–10(c) show the time-space diagrams of lane 1, lane 2, and lane 3 at a high density, , respectively. The high density region appears and extends behind the car accident in lane 1 and lane 2. Correspondingly, a region with a very large increase in density appears in lane 3. At this initial high density, all cars behind the car accident in lane 1 and lane 2 need to change to lane 3 while the capacity of lane 3 is limited. Thus more cars in lane 1 and lane 2 are not able to change lane with time and the extended jams are formed behind the car accidents. On the other hand, the cars changing from lane 1 to lane 2 gradually lead to the extended high density region in lane 3. In addition, from Figures 9 and 10 we can see that the density in the downstream of the car accident in lane 2 and lane 3 is lower than that in lane 1. The reason is as follows. After the cars go through the car accident, many of them attempt to change from lane 3 to lane 1 where the driving condition is better. Thus the accident not only has a significant effect on the upstream of the car accident, but also considerably affects the downstream of the car accident.

(a)

(b)

(c)

(a)

(b)

(c)
Figure 11 shows the currents against the density when two car accidents occupy lane 1 and lane 2. When the density is low, vehicles move freely. With the increase of density, the current of lane 3 increases and then keeps a high value, while the current of lane 1 decreases and then keeps a low value, and the current of lane 2 almost saturates and keeps a constant value. The current of lane 2 is higher than of lane 1. When the density is higher, the currents decrease with increasing density.

The second situation is that two car accidents are located at in lane 1 and in lane 3, respectively. Figures 12(a)–12(c) illustrate the time-space diagrams of lane 1, lane 2, and lane 3 at a low density, , respectively. It can be seen that the localized jam is formed behind the car accident in lane 1 and lane 3, and a region with an increase in density arises in lane 2. Figures 13(a)–13(c) show the time-space diagrams of lane 1, lane 2, and lane 3 at a high density, , respectively. The high density region appears and extends behind the car accidents in lane 1 and lane 3, because there are few cars changing from lane 1 and lane 3 to lane 2 at this initial high density. Correspondingly, a region with a very large increase in density appears in lane 2. Comparing Figure 10 with Figure 13, we can see that the effect of the car accidents located at lane 1 and lane 3 on the traffic flow is greater than the case of the car accidents located at lane 1 and lane 2 because the cars behind the accident cars in lane 1 and lane 3 all try to change lane to lane 2, but this causes more serious traffic confusion.

(a)

(b)

(c)

(a)

(b)

(c)
Figure 14 shows the currents against the density when two car accidents occupy lane 1 and lane 3. The current increases with density at a low density. With the increase of density, the currents of lane 1 and lane 3 decrease and then keep at a very low value because of the existence of the car accident, while the current of lane 2 almost saturates and keeps a high value. When the density is higher, the currents of each lane decrease. The current on lane 3 is consistent with that on lane 1 since the symmetric lane-changing rules are adopted.

4. Conclusions
The car accidents have an important influence on traffic flow. To address this problem, we have presented an extended three-lane traffic model to investigate the traffic jam caused by the car accidents. We obtain the time-space diagrams and the current diagrams for the cases that the car accidents occupy lane 1, lane 2, lanes 1 and 2, or lane 1 and lane 3. Results indicate that the effect of the car accident located at lane 1 on the traffic flow is greater than the case of the car accident located at lane 2, while the car accidents have a more serious effect on the whole road when they occupy two lanes. These results are in coincidence with our daily experiences in real traffic. In this paper, the proposed model is used to explore the effect of the car accidents upon the traffic flow on a highway. In our future work, this model will be extended to city traffic and to study the traffic behaviors discussed in the works of Toledo et al. [35], Varas et al. [36], Pastén et al. [37], and Toledo et al. [38].
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This paper is supported by the National Natural Science Foundation of China (Grant no. 11102165), the Natural Science Basis Research Plan in Shaanxi Province of China (Grant nos. 2012JM1001 and 2013JQ7014), and NPU Foundation for Fundamental Research (Grant no. NPU-FFR-JC201254). The authors would like to thank two anonymous referees for their valuable suggestions on improving the paper.