A Node Model Capturing Turning Lane Capacity and Physical Queuing for the Dynamic Network Loading Problem

An analytical dynamic node-based model is proposed to represent flows on a traffic network and to be utilized as an integral part of a dynamic network loading DNL process by solving a continuous DNL problem. The proposed model formulation has an integrate base to be structured with a link load computing component, where physical queuing and its influence were explicitly taken into account by dividing a link into two parts: running part and queuing part. The solution to the model is obtained by a hybridization algorithm of simulation and analytical approach, where an iteration process is conducted to update time-dependent network flow conditions after a reasonable discretization of the problem. The performance of the proposed model, as a DNL model, is tested on a sample network. It is seen that the proposed model provides consistent approximations to link flow dynamics. The dynamic node model proposed in this paper is unique in that it explicitly models directional queue in each turning lane and the First-In-First-Out FIFO rule at lane level rather than link level is pursued.


Introduction
The dynamic network loading DNL problem is the reproduction of variable link performances and network flow conditions by considering nodal and exiting flow characteristics 1-4 or by explicitly considering the flow propagation and the provision of real time pathlink information 5-7 , given path flows and link-performance functions.In the past two decades, DNL has been extensively studied owing to the needs of simulating urban traffic and solving dynamic traffic assignment DTA problems.The variation in model structure is heavily dependent on both assumptions made to obtain a solution for the problem, that is, the discretization dimension, queuing, and the criteria that affect the computation of link loads and path-link traveling times.One way to categorize the different approaches is discrete-flow models and continuous-flow models.Discrete-flow models, also referred to as simulations, are usually classified into microscopic simulations and macroscopic simulations.In this paper we concentrate on continuous-flow models which are also referred as "analytic" or "macroscopic" models, and we do so for the usual reasons: low number of parameters to be calibrated, good computational performance, mathematical tractability.Microsimulations are time-based, meaning that individual vehicles are moved in short-time intervals 0.1-1 s according to car following and other traffic behavior models 8 .The main assumption that is made during the construction of a model to solve DNL problem is on queuing and can be divided into point queuing and physical queuing.Considering the queuing assumption, the models of DNL can be further classified as link based, which vary with the adopted performance function such as link exit function and travel-time-function and node based.Although the performance-function-adopted-link-based modeling framework has been widely used for the analytical formulation of the DNL problem, little research has been conducted to explore their capabilities of reproducing realistic flow dynamics 9-11 .
The term node based is generally used for models that explicitly consider the flow splitting rates which are the proportions of traffic leaving a node, assigned to each exiting link, and has been originated from the pioneering proposition of the cell-based traffic flow model in 12 : the cell transmission model CTM .Node models can be classified as follows according to how link interactions are modeled: 1 competition-free nodes: only the flow conservation law is obeyed at such nodes.The competition-free node model is often seen in the analytical DTA research 13 .2 Uncontrolled competition nodes: traffic from different incoming links and/or heading to different outgoing links would have to compete against each other for the limit capacities 14 .A typical example is freeway junctions on-and offramps without metering facilities.3 Controlled competition nodes: the competition among different traffic streams is managed by a controlled logic, such as signalized intersections.No matter if controlled or uncontrolled competition nodes, the local demand and supply flow concept are utilized to put up an unifying framework for the modeling of intersections simultaneously with the imposition of boundary conditions to various network flow models such as first-order wave models in 12, 15 , second-order wave models in 16 , models that explicitly incorporate the spatial queuing effect 9 , and link performance models 17 .Most of the aforementioned node-based models are proposed to overcome the deficiencies of linkbased models and to be utilized in network traffic control and management applications, where merges and diverges excessively break down the stability of traffic flow.
However, some details that are really important and strongly influence loading results have not got satisfactory representation and enough attention.For example, most node models deal with flow propagation through nodes by distributing local demand of upstream links to downstream links constrained by prevailing supply, where the prevailing demand of a link depends on flows waiting to exit and the capacity of the link rather than the capacity of turning lanes.Queuing on an upstream link is treated as a whole rather than capturing turning directions of vehicles or directional lanes they belonging to.To satisfy FIFO at link level, such models are liable to cause unrealistic description of flow dynamics.For example, rightturning vehicles have to wait even if its downstream link has enough supply just because left-turning vehicles cannot cross the intersection for lacking of downstream supply, whether these vehicles right turning and left turning are in the same directional turning lane.In real traffic conditions, as pointed in HCM, flows propagating through an intersection are limited by the capacity of directional entry lanes, and drivers generally choose a turning lane according to their turning directions when approaching an intersection, which leads to several directional queues with different flow dynamics even on the same link.
Being cognizant of the insufficiency of link-based models in reproducing actual traffic dynamics and the motivation for node-based models, we have proposed an analytical dynamic node-based model that explicitly considers the influence of directional turning lane and physical queuing in different lanes.In our model, queues in each turning lane are explicitly modeled and FIFO at lane level is pursued.We have chosen to solve the continuous DNL problem analytically with our proposed methodology rather than formulating a theoretical high-order wave model for a node.The proposed node model has an integrate base to be structured with a link load computing component.The constraints of link dynamics, flow conservation, flow propagation, and boundary conditions are considered both in nodal rules and the link model.We obtained the solution of the model by coding a simulation-based hybridization algorithm after designing a discrete version of the problem.
This paper is organized as follows.In Section 2, model description, including the link model and the node model formulation is given.The solution procedure is explained in Section 3. Section 4 holds the numerical results of the solution method that is employed to solve the DNL problem on a sample network.Findings and discussions based on the results obtained conclude this paper as the final section.

Model Description
Given path flow and link performance functions, the CDNLP consists of determining timedependant network flow conditions such as link travel times, link inflows, and link outflows.In the loading procedure, the modeling of traffic flow has two major facets: the representation of traffic dynamics on a link homogeneous road segment and on a node boundary of several links, intersection .We propose an approach that evolves a link model with a node model, where a link is divided into two parts: the running part and the queuing part.The partial length of the running part is dependent on the prevailing traffic condition.Flow dynamics on the running part are described by a Travel Time Function TTF-based link model, and flow propagation through intersections is represented with a node model that describes directional queuing in different turning lanes explicitly.Without loss of generality, all links are assumed to have exclusive turning lanes for each turning direction at intersections.Time-dependent exit flow of a directional lane rather than link exit flow is modeled and calculated as the minimum value of prevailing demand and supply of the lane.For links holding mixed turning lanes, exit flow from mixed lanes involves distributing lane demand to downstream links, which can be dealt with using methods adopted in previous node models 18, 19 by taking a mixed lane as a link.
Some basic notations and variables are given first.The physical traffic network can be represented by a directed network G N, A , where N is the set of nodes, and A is the set of links.In the following, the index r denotes an origin node, the index s denotes a destination node, the index a denotes a link, the index p denotes a path between the origin O and the destination D , and the index aa denotes a directional turning lane which connects link a and its downstream link a .An O-D pair r, s is designated, where O ⊂ S, D ⊂ S, r ∈ O, s ∈ D, rs ∈ R, and R ⊂ O * D .The ⊂ of paths between an O-D pair r, s is denoted by p rs .

Link Model
Link models for the DNL problems enable the specification of flow dynamics on a link in three ways: bottleneck models, whole link function models, and hydrodynamic models.In Mathematical Problems in Engineering bottleneck-type models, vehicles always move along a link at free-flow speed before they arrive at the exit node, where they form a FIFO queue if the outflow rate they induce exceeds the maximum discharge rate bottleneck capacity of the link 19 .In whole-link function models, a performance function Travel-Time-based Function or Link Exit Function is adopted to describe flow dynamics taking the whole link as homogeneous 20 .While hydrodynamic models view traffic as a continuous fluid represented by density, speed, and flowrate; they are also known as kinematic wave KW models because their solutions can be categorized by combinations of kinematic waves in any of the three quantities 21 .We follow the approach of bottleneck models because flow dynamics of vehicles running on homogeneous segments and of those waiting in a queue for exiting should be treated differently.Also, some improvements are made to model congestions on a link segment and to capture the effect of physical queuing.In our approach, a link is divided into two parts: the running part and the queuing part.The prevailing partial length of the running part of link a denoted with L a t depends on the number of vehicles on the link and is calculated as below: where v a t denotes the travel speed for vehicles entering link a at time t and is calculated with a modified Green-Shields equation 22 as given by 2.3 : where v a, min , v a, f , k a,j denote the minimum speed, free-flow speed, and jam density of link a, respectively.The prevailing density of the running part of link a at time t denoted with k a t is determined by the following equation: where lan a denotes the lane number of link a.
The path-specific exit flow of the running part of link a is determined by 2.5 :

Node Model
We only consider nodes with competition.Most node models presented in previous research deal with flow propagation through nodes by distributing demand of upstream links to downstream links.The demand of a link, D, is the maximum possible exit flow rate, that is, D min{C, Q}, 2.9 the supply of a link, S, is the maximum possible receiving flow rate, that is, where C is the flow exit capacity; Q is the rate of flow ready to exit; R is the maximum entry flow rate to the downstream link permitted by the current traffic condition.The partial flow exit from an upstream link entering a downstream link is calculated according to the turning proportion and constrained by prevailing supply.Taking Figure 1 as an example, the queue on link a or link b is modeled as a whole, time-dependent partial flow exiting from link a to link c or d is calculated as follows: exit flow from link a is determined first with the formula mentioned above, and the partial flow is calculated with the total exit flow from link a multiplying by a turning proportion.To satisfy FIFO rule at link level, such models may lead to unnecessary queuing of vehicles.For example, vehicles on link a ready to enter link c will have to wait even if link c has enough supply, just for the reason that link d got a traffic jam and cannot allow any entering flow at current time.In other words, vehicles with different turning directions will have to travel synchronously to satisfy FIFO at link level.
In real traffic conditions, as pointed in HCM, there exist directional turning lanes at intersections especially those with signal controlled, and the capacity of turning lanes rather than that of the whole link actually plays role in limiting flow exiting a link.Drivers generally change or choose turning lanes according to their turning directions when approaching an intersection, which leads to different queues with different flow dynamics in turning lanes.To overcome this shortcoming in describing real dynamics of traffic flow, queuing in different turning lanes should be specifically described, and FIFO at lane level rather than link level should be pursued.
In such situation, the input to the node model is time-dependent entering flow to the queuing part of a link in each turning lane, which can be determined by the exit flow of the running part of the same link according to path flows as given below: where uq aa rsp t is the partial entering flow to the queuing part of link a in turning lane aa of rsO-D pair on path p at time t, and uq aa t is the total entering flow to the queuing part of link a in turning lane aa at time t.
Based on the entering flow, the exit flow of the directional queuing part can be calculated with prevailing demand and supply as following: The splitting rate λ p aa for the flow that is exiting link a in turning lane aa along path p is calculated as given by 2.18 .The constraint associated with λ p aa is given by 2.19 , which expresses that at turning lane level, the FIFO rule holds.

Solution Procedures
The solution procedures sought for the DNL problem can be clustered into three approaches: 1 an analytical approach, 2 a simulation-based approach, and 3 a hybridization of these two approaches.The DNL problem solved with an analytical approach possesses the desired properties of both the solution and the link performance function 23 .The simulationbased approach, that is, INTEGRATION, DYSMART, and other propositions, provides more flexibility and consequently enables the reproduction of detailed flow dynamics both in dealing with more complex traffic patterns and in obtaining more realistic loading results lacking the desired solution properties 24 .The hybrid approach set out by the incorporation of an analytical formulation with simulation and utilized in several studies, including 12, 22, 25 and this paper, has become an efficient alternative in network performance modeling.The feature that it simultaneously possesses both ability of analytical models in obtaining more accurate solutions and the ability of simulation models in obtaining computationally efficient solutions in network traffic modeling context has motivated us to employ a hybrid approach in this paper.The analysis period is divided into small time intervals with same length σ.The interval σ is treated as unit of time, and time period that is not an exact multiple of σ involved in loading procedure is processed with linear approximation.Given path flows during O-D traffic period, we can get entering flow to the running part of each link by 2.6 .
The following algorithm summarizes the determination of flows and the number of vehicles on all links by the proposed node-based modeling approach when a network is considered.9 Determine the disaggregate variable vq aa rsp t and xq aa rsp t with 2.17 and 2.20 , respectively.
10 If the demand is discharged from all the paths during the analysis period, stop; otherwise, set t following time interval and go to step 2 .

Sample Network
The performance of the proposed model is tested on a sample network shown in Figure 2 with several paths between the given O-D pair.The link characteristics assigned are given in Table 1, and the characteristics of turning lane s at the entrance of each intersection are given in Table 2.The traffic flow interval, σ, is 1 min, and the departure period, T , is 30 min and is divided into 5 intervals with same length.There is only one O-D pair and seven routes used in this network.Route departure flows, f rsp t , have a constant value during each departure interval, given in Table 3 and f rsp t 0 at other times.

Test Results
The performance of the proposed model is evaluated with a number of critical terms from the simulation.The outflow diagram of flow propagation at node 5 is given in Figure 3.The dynamic outflow and queue length diagrams in each turning lane on link d are given in turn in Figures 4 and 5.
In Figure 3, it is seen that the loadings to link k required longer time, due to a longer flow profile of paths involving this link.It can be seen from Figure 4 that outflows of different turning lanes on link d appear to vary in similar manners, which reach saturation in a few minutes and keep the state for some time 41 minutes in the right turn lane, 52 minutes in the left turn lane, and 69 minutes in the through lane before a sharp decline.From Figure 5,   we can see that the number of queuing vehicles in each turning lane increases to a maximum value after the entrance lanes got saturation.

Comparison of Different Models
In this part, link d was implemented using other models for comparison.Travel-Time-Function-based link models TTF that has been used in 22 , node model pursuing FIFO at link level for simplicity named NM-1 adopted in 18 , and node model presented in this paper for simplicity named NM-2 are tested.In TTF-based link model, a link was treated as a whole, and the travel speed for a vehicle entering link a at time t was calculated with the modified Green-Shields model shown in 2.6 .The difference is that the traversal time was calculated with the total length of link a, not the partial length of the running part, divided by the speed at time t.NM-1 has a similar framework with NM-2, in which a link is divided into two parts and the flow dynamics on running part were described with a travel time function.It is the description of flow dynamics on queuing parts that makes NM-2 different from NM-1.In NM-1, time-dependent exit flow of a whole link but not a turning lane is calculated with  its prevailing demand and supply; queuing vehicles with different directions are described as mixed together; the prevailing demand waiting to exit is limited by the capacity of the link, rather than specified to turning lanes.
Figure 6 compares the time-dependent link traversal times on link d produced by the three models mentioned above.As shown, NM-1 overestimated the congestion on link d, while TTF led to a slight underestimation.The reason lies in that the TTF-based link model ignores queuing delay before intersections, and NM-1, pursuing FIFO at link level, may lead to unnecessary queuing of some vehicles.

Conclusion
In this paper, an analytical node-based model has been proposed for the continuous dynamic network loading problem CDNLP .Taking account of a set of analytical rules, an algorithm has been derived.The new dynamic node model proposed in this paper is unique in that it explicitly captures the turning lane capacity and directional physical queuing in the description of traffic propagation through nodes.
In the proposed node model, directional queue in each turning lane is explicitly considered and the FIFO rule at lane level rather than link level is pursued.A directional exit flow function is presented to update the exit flow from each turning lane with the constraints of lane capacity, downstream link capacity, and flow conservation, and so forth.A travel time function TTF-based link model has been evolved to a node model in the modeling structure.To capture queue spillback caused by the capacity of turning lanes as well as the downstream link, the concept of prevailing partial length is adopted, which tells the length of partial link that can be used by running vehicles and is determined by the number of running vehicles as well as queuing vehicles.The time to travel the running part of a link equals to the partial length divided by the prevailing speed calculated with a modified Green-shields formulation and is common to vehicles entering the link at the same time.An algorithm is presented for the CDNL problem with our node-based modeling approach, which is based on simulation and updates the performance of running parts as well as queuing parts in order of time step.
The drawback is that the node model integrates a priori mean effects on flows of a traffic signal without explicitly representing the alternation of green and red stages, and traffic wave can not be tracked, for the running part of a link is treated as a homogeneous road segment.
The proposed node-based model can be easily integrated as a flow modeling component of a dynamic traffic-assignment process, enabling its utilization in a wide range of intelligent transportation system applications.For example, the realistic representation of traffic flow dynamics enables the proposed model to be easily utilized in advanced traveler information systems ATISs .The prediction on link performances can be obtained in terms of realtime flow volume data inputs.These predictions can be basic inputs to ATIS applications, such as variable message signs for route guidance.Dynamic signal optimization and ramp metering are other possible topics that the proposed model's extensions can study for capacity management and speed regulation.

Figure 2 :
Figure 2: Sample network used for testing trials.

Figure 3 :
Figure 3: Outflows at node 5 of sample network.

Figure 4 :
Figure 4: Outflows from turning lane s of link d.
Queuing vehicles in right turn lane Queuing vehicles in left turn lane Queuing vehicles in through lane

Figure 5 :Figure 6 :
Figure 5: Queue length in turning lane s of link d.
a t and xq a t denote the number of vehicles on the running part and the queuing part of link a at time t.The traversal time on the running part of link a experienced by vehicles entering at time t denoted with dr a t is a function of link volume as given in 2.2 : * aa t is the number of vehicles in a turning lane of link a with longest queue at time t; xr a rsp t and vr a rsp t denote entering flow and exit flow of the running part of link a of rsO-D pair on path p at time t.ur a rsp t can be determined beforehand as below: is the exit flow from the queuing part of link a − in turning lane a − a of rsO-D pair on path p at time t; f rsp t is the departure flow rate on path p between an rsO-D pair departing at time t.Equation 2.6 ensures the FIFO behavior on the running part by forcing vehicles that enter the link at t to be pushed out at t dr a t .The number of vehicles existing on the running part of link a is updated accordingly as given by 2.7 and 2.8 : a t − dr a t − 1 , ∀r, s, t, a, ∀p ∈ p rs , 2.5where ur rsp t denotes the number of vehicles on the running part of link a of rsO-D pair on path p at time t.
denotes the number of vehicles on the queuing part of link a in turning lane aa at time t; c aa denotes the capacity of turning lane aa on link a.The partial supply of link a allocated to link a at time t denoted with R aa t can be calculated by 2.15 , as shown below: aa is the proportionality coefficient which depends on the lane number and control mode of an intersection and treated as a constant; h a is the maximum number of vehicles that link a can accommodate; s a is the entering capacity of link a ; a h is the head node of link a; a t is the tail node of link a .With partial exit flow of each turning lane known, the total exit flow of the queuing part of link a denoted with vq a t can be determined by 2.16 as given below:

1
Set t 0 first time interval , clear the network, that is, Calculate the aggregate inflow uq aa t with 2.12 and calculate the number of vehicles in each turning queue xq aa t with 2.21 .

Table 1 :
Assigned link characteristics of sample network.

Table 2 :
Capacities of turning lanes at intersection entrance.

Table 3 :
Route departure flow rates in 30 minutes.