# Dynamic Congested Traffic States of Density Difference Lattice Hydrodynamic Model with On-Ramp

Authors: Jun-fang Tian; Zhen-zhou Yuan; Bin Jia; Wang Tao
Journal: Discrete Dynamics in Nature and Society
 (2013)

---

## Abstract

The density difference lattice hydrodynamic model with on-ramp is proposed. Stochastic and deterministic rules for the on-ramp flow entering into the main road are designed. Under the stochastic rule, various empirical spatiotemporal patterns can be well reproduced, such as the pinned localized cluster (PLC) and homogeneous synchronized traffic (HST). Under the deterministic rule, four new types of congested traffic states, such as the stationary oscillatory congested traffic (OCT) upstream of on-ramp and the stationary and localized OCT emerging in HST, are identified for the first time. Comparisons with the macroscopic speed gradient model are carried out, and some advantages of our model are clarified. To our knowledge, it is the first research showing that the lattice hydrodynamic model could reproduce various congested patterns induced by the on-ramp system.


---

## Body

## 1. Introduction

With the rapid development of urbanization, traffic congestion becomes one of the most serious problems that undermine the operation efficiency of modern cities. In order to understand the mechanism of traffic congestion, many models and analyses have been carried out to explain the empirical findings [1–7], such as the hysteresis, stability, and phase transitions. Usually, there are two types of traffic models, that is, microscopic and macroscopic models. Microscopic models, including cellular automata [8] and car-following models [9–11], focus on the dynamics of individual vehicles. Macroscopic models, including gas-kinetic-based models [12], continuum models [13, 14] and lattice hydrodynamic models [15–29], regard the whole traffic flow as a flow of continuous medium composed of vehicles. Individual vehicles do not appear explicitly in macroscopic models.Recently, phase transitions induced by the on-ramp have been observed and studied intensively. Helbing et al. [30–36] found that traffic states in the open systems with bottlenecks are mainly determined by the instability diagram. The preconditions for the possible occurrence of different traffic states can be illustrated by the phase diagram. These empirical traffic states include the pinned localized cluster (PLC), moving localized cluster (MLC), triggered stop-and-go traffic (TSG), oscillating congested traffic (OCT), homogeneous synchronized traffic (HST), and homogeneous congested traffic (HCT); see Figure 1.Examples of congested traffic states measured on the German freeway A5 close to Frankfurt. (a) Pinned localized cluster (PLC), (b) Moving localized cluster (MLC), (c) Triggered stop-and-go waves (TSG), (d) Oscillating congested traffic (OCT), (e) Homogeneous synchronized flow (HST), and (f) Homogeneous congested traffic (HCT). This figure is taken fromhttp//:www.traffic-states.com/.
(a)
Pinned localized cluster (PLC) A5 South, September 19, 2001
(b)
Moving localized cluster (MLC) A5 South, June 06, 2001
(c)
Triggered stop-and-go waves (TSG) A5 South, May 07, 2001
(d)
Oscillating congested traffic (OCT) A5 South, July 31, 2001
(e)
Homogeneous synchronized traffic (HST) A5 South, April 09, 2001
(f)
Homogenous congested traffic (HCT) A5 South, April 25, 2001In previous works, both microscopic and macroscopic models except the lattice hydrodynamic models were investigated under the on-ramp system. The lattice hydrodynamic models were firstly proposed by Nagatani [15, 16] and later extended by many other researchers [17–29], due to their simple structure, ease of theoretical analysis and simulation, and deep connection with the macroscopic continuum models and microscopic car-following models. Maybe the lacking of the lattice hydrodynamic model with on-ramp is the reason why the phase transitions and congested states induced by the on-ramp have never been investigated by these models.In this paper, we firstly propose the density difference lattice hydrodynamic model with on-ramp. Later two types of on-ramp system are designed. We show that not only all observed congested traffic states can be simulated by the model mentioned before but also some new congested states are found. The remainder of the present paper is organized as follows. In the following section, DDLM is introduced. In Section3, DDLM with on-ramp is proposed and the on-ramp systems for lattice hydrodynamic models are designed. In Section 4, numerical simulations are carried out. Finally, conclusions are given in Section 5.


## 2. Density Difference Lattice Hydrodynamic Model

In 1998, Nagatani [15, 16] proposed the first lattice hydrodynamic model. Nagatani’s model was derived by discretizing the following continuum model:(1)∂tρ+∂x(ρv)=0,(2)∂t(ρv)=aρ0V(ρ(x+δ))-aρv,
where a is the sensitivity of the driver and δ represents the average headway, δ=1/ρ0. The ρ(x+δ) is the local density at position x+δ at time t, which is related to the inverse of headway h(x,t); that is, ρ(x+δ)=1/h(x,t). V(ρ(x+δ)) is the optimal velocity.In Nagatani’s original paper [15, 16], ρ0 is defined as the average density. Then (2) can describe correctly steady-state traffic flow (subscript e): ρ0=ρ, ρ(x+δ)=ρ(x), and d(ρv)/dt(ρv)=0, so ve(ρ)=V(ρ) which is consistent. However, in this interpretation, a definition of the averaging procedure to obtain ρ0 from ρ has not been presented. Moreover the real traffic system is an open system, which leads to difficulties formulating this procedure. Therefore, we think ρ0 should be considered as a constant parameter. Although it makes the steady-state relationship not consistent since from (2) it follows that ve(ρ)=ρ0/ρV(ρ) which diverges for ρ to zero rather than tending to V(0), but it is definite that the traffic is in the very low density free flow state if ρ is smaller than ρ due to the small value of ρ0. We think this state is not the focus of traffic flow, which is the congested traffic. Thus, this deficiency will have little influence on the application of the model.The idea of this model is that a driver adjusts the vehicle velocity according to the observed headwayh(x,t) or density ahead ρ(x+δ). This is similar to the idea in the car-following models with the optimal velocity V(h(x,t)). Then, for dimensionless space x (let x~=x/δ, and x~ is indicated as x hereafter), the lattice model is expressed as:
(3)∂tρj+ρ0(ρjvj-ρj-1vj-1)=0,∂t(ρjvj)=aρ0V(ρj+1)-aρjvj,
where j denotes the site on the one-dimensional lattice and ρj and vj denote the local density and the local average velocity on site j at time t, respectively.Subsequently, many extended models [17–29] have been proposed to suppress traffic jams. In particular, the density difference lattice hydrodynamic model (DDLM) [28, 29], by incorporating the density difference ahead into Nagatani’s model, shows a close relationship with the gas-kinetic-based model [30, 31] and the information that density difference plays a significant role in traffic flow. DDLM consists of continuity equation (4) and dynamic equation (5):
(4)∂tρj+ρ0(ρjvj-ρj-1vj-1)=0,(5)∂t(ρjvj)=aρ0V(ρj+1)-aρjvj+λ(ρj-ρj+1)ρ0,
where λ is the reaction coefficient to the density difference. By eliminating v in (4) and (5), the following equation is obtained:
(6)∂t2ρj+aρ02(V(ρj+1)-V(ρj))+a∂tρj0000+λ(2ρj-ρj+1-ρj-1)=0.
One can see that DDLM has considered the following lattice j-1. For the convenience of simulation, we rewrite (6) into the difference form:
(7)ρj(t+2τ)=2ρj(t+τ)-ρj(t)-aτ2ρ02(V(ρj+1(t))-V(ρj(t)))-aτ(ρj(t+τ)-ρj(t))-λτ2(2ρj(t)-ρj+1(t)-ρj-1(t)).
The optimal velocity is chosen as follows [16]:
(8)V(ρ)=vmax2[tanh(2ρ0-ρρ02-1ρc)+tanh(1ρc)],
where vmax is the maximum speed and ρc is the critical density. This function has the turning point at ρ=ρ0=ρc. Typical values of these parameters are ρc=ρ0=0.25 and  vmax=2 m/s [15, 16]. Figure 2 is the schematic illustration of the speed and flow as a function of the density. According to the maximum value of flow on the flow-density curve, free traffic and congested traffic can be classified, which is applied to distinguish the traffic state in the following simulations.Schematic illustration of speed and flow as a function of the density forρc=ρ0=0.25,  vmax=2 m/s. (a) The speed-density curve. (b) The flow-density curve.
(a)
(b)


## 3. Density Difference Lattice Hydrodynamic Model with On-Ramp

In previous studies [15–29], theoretical analysis and numerical simulations for lattice hydrodynamic models focus on suppressing traffic jams and linear and nonlinear stability analysis. Traffic states induced by the on-ramp systems have not been investigated.In order to simulate the empirical phenomena observed by Helbing et al. [30–36], two simple methods are designed to model the on-ramp system for the lattice model. In Figure 3(a), the position of the on-ramp is located at xon and the region [xon,xon+Lramp] is selected as the inserting area of the flux from on-ramp, where Lramp is the length of the on-ramp. At each time step, a lattice s in this region will be selected randomly. However, the length of the inserting area is set as Lramp=1 in Figure 3(b); that is, the flux from on-ramp could enter the main road by only one lattice. According to the different inserting rules of the on-ramp flux, Figures 3(a) and 3(b) are simply named as the “stochastic on-ramp” and “deterministic on-ramp,” respectively.The sketches of the on-ramp systems. (a) The stochastic on-ramp, (b) the deterministic on-ramp.
(a)
(b)The on-ramp fluxQon(t) will be inserted into the selected lattice. The density difference lattice model with the on-ramp could be described as follows:
(9)∂tρj+ρ0(ρjvj-ρj-1vj-1)=ρ0Q¯j,∂t(ρjvj)=aρ0V(ρj+1)-aρjvj+λ(ρj-ρj+1)ρ0,
where
(10)Q¯j(t)={Qon(t):ifj=s,0:ifj≠s.By eliminatingv in (9) and writing the results into the difference form, the following model is obtained:
(11)ρj(t+2τ)=2ρj(t+τ)-ρj(t)-aτ2ρ02(V(ρj+1(t))-V(ρj(t)))-aτ(ρj(t+τ)-ρj(t))-λτ2(2ρj(t)-ρj+1(t)-ρj-1(t))+aτ2ρ0Qj(t).


## 4. Numerical Simulations

The initial conditions for numerical simulations are set as follows. The main road is composed byN=400 lattices and the open boundary condition is applied; that is,
(12)ρ1(t)=ρ2(t),ρN-1(t)=ρN(t).
The main road is prepared to start with a homogeneous traffic flow, whose density is ρh. The on-ramp region is located at xon=300. The length of the stochastic on-ramp is Lramp=10 lattices. The sensitivity a is set as a piecewise constant function: a=0.5 for ρj≤0.19 and a=1 for ρj>0.19. Other parameters are set as follows: τ=0.1,  ρc=ρ0=0.25,  vmax=2 m/s, and  λ=0.2.


### 4.1. Stochastic On-Ramp Simulation

In Figure4, the congested traffic is localized at the on-ramp, which is characterized by the localized reduction of speed and increase of density. Thus, it is the pinned localized cluster (PLC). The downstream front of the PLC is fixed at the bottleneck and its upstream front does not propagate upstream continuously over time. Upstream and downstream PLC, there is free traffic flow. Figure 5 shows the moving localized cluster (MLC). The traffic is congested with a single density wave. This density wave is a compression wave and propagates upstream with a constant speed. In contrast to the PLC, both the downstream and upstream fronts of MLC propagate on the road. Figure 6 shows the triggered stop-and-go traffic (TSG). Compared with MLC, there are several moving clusters on the road. Figure 7 exhibits the oscillatory congested traffic (OCT). Density waves emerge at the downstream and propagate upstream continuously. Figure 8 shows the homogeneous synchronized traffic (HST). The downstream front of HST is fixed at the downstream boundary and the upstream front of HST propagates upstream continuously over time. The density profile in Figure 9 is homogeneous over space except for the neighborhood of the downstream boundary. It is the homogeneous congested state (HCT). Differently from HST, vehicles’ speed is very low due to the high density in HCT.The spatiotemporal evolutions of PLC forρh=0.08,  Qon=0.13.
(a)
(b)The spatiotemporal evolutions of MLC forρh=0.05,  Qon=0.1.
(a)
(b)The spatiotemporal evolutions of TSG forρh=0.07,  Qon=0.11.
(a)
(b)The spatiotemporal evolutions of OCT forρh=0.25,  Qon=0.12.
(a)
(b)The spatiotemporal evolutions of HST forρh=0.09,  Qon=0.09.
(a)
(b)The spatiotemporal evolutions of HCT forρh=0.25,  Qon=0.18.
(a)
(b)Furthermore, it should be noted that the congested clusters simulated by DDLM emerge in the congested traffic flow rather than free traffic flow. Comparing with Figure1, one can see that the empirical congested traffic states can be well reproduced by DDLM.


### 4.2. Deterministic On-Ramp Simulation

In this section, we investigate the congested traffic states induced by the deterministic on-ramp. We found that HST and HCT can be reproduced, but the PLC, TSG, and MLC cannot be simulated. Moreover, the following new types of congested traffic states have been found.Figure10 shows the moving congested states. Its upstream front is composed by several clusters, that is, MLCs, and there are two clusters in Figure 10. Oscillations make up its downstream front, that is, moving OCT. Free flow exits between MLCs and OCT. In the other regions of the road, the HST is maintained. Thus, MLCs, OCT, free flow, and HST coexist in this congested state. Figure 11 shows another type of moving congested state. MLCs firstly occur in HST in this state. Then, free flow appears in the downstream of MCLs. Later, the moving OCT emerges in free flow. Usually, several regions of moving OCT and free flow alternately arise in this congested state. Figure 12 shows the free flow, OCT. and HST coexistence state. Before this state is formed, MLCs occur in HST. Then the OCT emerges far from the on-ramp. In the downstream front of OCT, HST is maintained. In the upstream front of OCT, free flow is preserved. The OCT is stationary and localized on the road. Figure 13 shows the HST and OCT coexistence state. The downstream front of OCT is away from the on-ramp; that is, this OCT is a stationary OCT. Downstream of the front, HST is maintained.The spatiotemporal evolutions forρh=0.05,  Qon=0.15.
(a)
(b)The spatiotemporal evolutions forρh=0.05,  Qon=0.18.
(a)
(b)The spatiotemporal evolutions forρh=0.12,  Qon=0.14.
(a)
(b)The spatiotemporal evolutions forρh=0.18,  Qon=0.18.
(a)
(b)Moreover, we compare these results with those of other models, such as the macroscopic speed gradient model (SGM) [37]. The following conclusions can be made. Firstly, all empirical congested states can be simulated by both models. However, the simulated states of DDLM are more in accordance with the empirical data, such as the simulated PLC; TSG by SGM is only composed by two wave fronts. No congested regions exist between these fronts, which is inconsistent with the empirical PLC and TSG (see Figure 1). Secondly, the congested states (Figures 10–12) of DDLM induced by the deterministic on-ramp cannot be reproduced by SGM. Thirdly, similar state of Figure 13 is reproduced by SGM, which is made up by the stationary OCT and HCT.


## 4.1. Stochastic On-Ramp Simulation

In Figure4, the congested traffic is localized at the on-ramp, which is characterized by the localized reduction of speed and increase of density. Thus, it is the pinned localized cluster (PLC). The downstream front of the PLC is fixed at the bottleneck and its upstream front does not propagate upstream continuously over time. Upstream and downstream PLC, there is free traffic flow. Figure 5 shows the moving localized cluster (MLC). The traffic is congested with a single density wave. This density wave is a compression wave and propagates upstream with a constant speed. In contrast to the PLC, both the downstream and upstream fronts of MLC propagate on the road. Figure 6 shows the triggered stop-and-go traffic (TSG). Compared with MLC, there are several moving clusters on the road. Figure 7 exhibits the oscillatory congested traffic (OCT). Density waves emerge at the downstream and propagate upstream continuously. Figure 8 shows the homogeneous synchronized traffic (HST). The downstream front of HST is fixed at the downstream boundary and the upstream front of HST propagates upstream continuously over time. The density profile in Figure 9 is homogeneous over space except for the neighborhood of the downstream boundary. It is the homogeneous congested state (HCT). Differently from HST, vehicles’ speed is very low due to the high density in HCT.The spatiotemporal evolutions of PLC forρh=0.08,  Qon=0.13.
(a)
(b)The spatiotemporal evolutions of MLC forρh=0.05,  Qon=0.1.
(a)
(b)The spatiotemporal evolutions of TSG forρh=0.07,  Qon=0.11.
(a)
(b)The spatiotemporal evolutions of OCT forρh=0.25,  Qon=0.12.
(a)
(b)The spatiotemporal evolutions of HST forρh=0.09,  Qon=0.09.
(a)
(b)The spatiotemporal evolutions of HCT forρh=0.25,  Qon=0.18.
(a)
(b)Furthermore, it should be noted that the congested clusters simulated by DDLM emerge in the congested traffic flow rather than free traffic flow. Comparing with Figure1, one can see that the empirical congested traffic states can be well reproduced by DDLM.


## 4.2. Deterministic On-Ramp Simulation

In this section, we investigate the congested traffic states induced by the deterministic on-ramp. We found that HST and HCT can be reproduced, but the PLC, TSG, and MLC cannot be simulated. Moreover, the following new types of congested traffic states have been found.Figure10 shows the moving congested states. Its upstream front is composed by several clusters, that is, MLCs, and there are two clusters in Figure 10. Oscillations make up its downstream front, that is, moving OCT. Free flow exits between MLCs and OCT. In the other regions of the road, the HST is maintained. Thus, MLCs, OCT, free flow, and HST coexist in this congested state. Figure 11 shows another type of moving congested state. MLCs firstly occur in HST in this state. Then, free flow appears in the downstream of MCLs. Later, the moving OCT emerges in free flow. Usually, several regions of moving OCT and free flow alternately arise in this congested state. Figure 12 shows the free flow, OCT. and HST coexistence state. Before this state is formed, MLCs occur in HST. Then the OCT emerges far from the on-ramp. In the downstream front of OCT, HST is maintained. In the upstream front of OCT, free flow is preserved. The OCT is stationary and localized on the road. Figure 13 shows the HST and OCT coexistence state. The downstream front of OCT is away from the on-ramp; that is, this OCT is a stationary OCT. Downstream of the front, HST is maintained.The spatiotemporal evolutions forρh=0.05,  Qon=0.15.
(a)
(b)The spatiotemporal evolutions forρh=0.05,  Qon=0.18.
(a)
(b)The spatiotemporal evolutions forρh=0.12,  Qon=0.14.
(a)
(b)The spatiotemporal evolutions forρh=0.18,  Qon=0.18.
(a)
(b)Moreover, we compare these results with those of other models, such as the macroscopic speed gradient model (SGM) [37]. The following conclusions can be made. Firstly, all empirical congested states can be simulated by both models. However, the simulated states of DDLM are more in accordance with the empirical data, such as the simulated PLC; TSG by SGM is only composed by two wave fronts. No congested regions exist between these fronts, which is inconsistent with the empirical PLC and TSG (see Figure 1). Secondly, the congested states (Figures 10–12) of DDLM induced by the deterministic on-ramp cannot be reproduced by SGM. Thirdly, similar state of Figure 13 is reproduced by SGM, which is made up by the stationary OCT and HCT.


## 5. Conclusions

Congested traffic states induced by the on-ramp have been observed and investigated widely by the macroscopic and microscopic traffic flow models. However, these have never been studied by the lattice hydrodynamic models.In this paper, we have proposed the density difference lattice hydrodynamic model with on-ramp, in order to study the congested traffic states induced by the on-ramp. Stochastic and deterministic on-ramps are designed. Under the stochastic on-ramp, various empirical congested states are reproduced, such as the pinned localized cluster (PLC), moving localized cluster (MLC), triggered stop-and-go traffic (TSG), oscillating congested traffic (OCT), homogeneous synchronized traffic (HST), and homogeneous congested traffic (HCT). Under the deterministic on-ramp, four new types of congested traffic states are identified, such as the moving congested states whose upstream and downstream fronts are composed by MLCs and moving OCT, respectively, the moving congested state composed by several regions of moving OCT and free flow, the congested state that the free flow, the stationary and localized OCT and HST coexistence state, and the HST and OCT coexistence state. It is the first time that these new congested states have been found. Empirical data should be explored to demonstrate whether they exist in the real traffic. At last, comparisons with the macroscopic speed gradient model are carried out. We found that the simulation results of DDLM are more in accordance with the real traffic.At last, it should be noted that the essential difference between the stochastic and deterministic on-ramps is the form of perturbations imposed on the main road. According to results in Section4, we can conclude that different form of perturbations, not the amplitudes of perturbations, can lead to different types of congested traffic states. Thus, we infer that the empirical findings of three-phase theory may be reproduced if a suitable form of perturbations is applied.



---

*Source: 941238-2013-10-03.xml*