Macroscopic Modeling and Simulation of Managed LaneFreeway Networks
Abstract
To help mitigate road congestion caused by the unrelenting growth of traffic demand, many transit authorities have implemented managed lane policies. Managed lanes typically run parallel to a freeway’s standard, generalpurpose (GP) lanes, but are restricted to certain types of vehicles. It was originally thought that managed lanes would improve the use of existing infrastructure through incentivization of demandmanagement behaviors like carpooling, but implementations have often been characterized by unpredicted phenomena that is often to detrimental system performance. Development of traffic models that can capture these sorts of behaviors is a key step for helping managed lanes deliver on their promised gains.
Towards this goal, this paper presents several macroscopic traffic modeling tools we have used for study of freeways equipped with managed lanes, or “managed lanefreeway networks.” The proposed framework is based on the widelyused firstorder kinematic wave theory. In this model, the GP and the managed lanes are modeled as parallel links connected by nodes, where certain type of traffic may switch between GP and managed lane links. Two types of managed lane configuration are considered: fullaccess, where vehicles can switch between the GP and the managed lanes anywhere; and separated, where such switching is allowed only at certain locations called gates.
We incorporate two phenomena into our model that are particular to managed lanefreeway networks: the inertia effect and the friction effect. The inertia effect reflects drivers’ inclination to stay in their lane as long as possible and switch only if this would obviously improve their travel condition. The friction effect reflects the empiricallyobserved driver fear of moving fast in a managed lane while traffic in the adjacent GP links moves slowly due to congestion.
Calibration of models of large road networks is difficult, as the dynamics depend on many parameters whose numbers grow with the network’s size. We present an iterative learningbased approach to calibrating our model’s physical and driverbehavioral parameters. Finally, we validate our model and calibration methodology with case studies of simulations of two managed laneequipped California freeways.
Keywords: macroscopic first order traffic model, first order node model, multicommodity traffic, managed lanes, HOV lanes, dynamic traffic assignment, dynamic network loading, inertia effect, friction effect
1 Introduction
Traffic demand in the developed and developing worlds shows no sign of decreasing, and the resulting congestion remains a costly source of inefficiency in the built environment. One study (Lomax et al., 2015) estimated that, in 2014, delays due to congestion cost drivers 7 billion hours and $160B in the United States alone, leading to the burning of 3 billion extra gallons of fuel. The historical strategy for accommodating more demand has been construction of additional infrastructure, but in recent years planners have also developed strategies to improve the performance of existing infrastructure, both through improved road operations and demand management, which seeks to lower the number of vehicles on the road (Kurzhanskiy and Varaiya, 2015). One such strategy that has been widely adopted in the United States and other developed countries is the creation of socalled managed lanes (Obenberger, 2004). Managed lanes are implemented on freeways by restricting the use of one or more lanes to certain vehicles. As an example, highoccupancyvehicle (HOV) lanes are intended to incentivize carpooling, which reduces the total number of cars on the road as a demand management outcome (Chang et al., 2008).
In addition to demand management, managed lanes provide an opportunity for improved road operations through realtime, responsive traffic control. For example, tolled express lanes give drivers the opportunity to pay a toll to drive parallel to the generalpurpose lanes on a (presumably lesscongested) express lane. Traffic management authorities here have an opportunity to adjust the toll amount in response to the realtime state of traffic on the network. The potential for managed lanes as instruments for reactive, realtime traffic operations–in addition to their demandmanagement purpose–has made them popular among transportation authorities (Kurzhanskiy and Varaiya, 2015).
However, the trafficoperational effects of managed lanes are not always straightforward or as rehabilitative as expected, as their presence can create complex traffic dynamics (Jang and Cassidy, 2012). Even in a freeway with simple geometry, the dynamics of traffic flow are complex and not fully understood, and adding managed lanes alongside the nonmanaged, generalpurpose (GP) lanes only exacerbates this. In effect, adding a managed lane creates two parallel and distinct, but coupled, traffic flows on the same physical structure. When used as intended, managed lanes carry flows with different densityvelocity characteristics and vehicletype (e.g., strictly HOVs) compositions than the freeway. When vehicles move between the two lane flows, complex phenomena that are unobserved in GPonly freeways can emerge (see e.g. Daganzo and Cassidy (2008); Liu et al. (2011); Jang and Cassidy (2012); Cassidy et al. (2015), and others).
Making better use of managed lanes requires an understanding of the macroscopic behavior they induce. One widelyused tool for understanding macroscopic traffic flow behavior is the macroscopic traffic flow model. A rich literature exists on macroscopic models for flows on long roads, and at junctions where those roads meet, but an extension to the parallelflows situation created by placing a managed lane in parallel with a freeway (a “managed lanefreeway network”) is not straightforward. Such a model should describe when vehicles enter or exit the managed lane, and capture the interactions caused by the two flows’ proximity, such as the socalled “friction effect” (Liu et al., 2011; Jang and Cassidy, 2012).
This paper presents macroscopic flow modeling tools we have used for simulation of managed lanefreeway networks. We begin in Section 2 with a discussion of relevant modeling tools from the literature, and how we make use of them. Section 3 describes network structures for the two common managed lane configurations: gatedaccess and fullaccess (ungated) lanes. This Section also describes models for friction and inertia effects that arise in managed lanes. Section 4 outlines how one may take these managed lanespecific constructions and add them to a general macroscopic traffic simulation process. Section 5 describes calibration methodologies for both network configurations. Section 6 presents case studies of two networks, one each of gated and fullaccess, and typical macroscopic simulation results.
2 Managed Lane Modeling
The modeling techniques presented in this paper are based on the firstorder “kinematic wave” macroscopic traffic flow model. These models describe aggregate traffic flows as fluids following a onedimensional conservation law. We briefly introduce our notation here, but do not discuss the basics of this class of models. Detailed reviews are available in many references.
2.1 Modeling basics
In this simulation framework, a road is divided into discrete cells, which we refer to as links. Links are drawn between nodes: a link begins at one node and ends at another. Many links may begin and end at each node. Each link is characterized by density , the number of cars in the link. In a firstorder model, the traffic flows are fully prescribed by the density. From timestep to , link ’s density updates according to the equation {linenomath}
(2.1) 
where is the length of link , is the flow (number of vehicles) leaving link and entering link at time , is the flow leaving link and entering link at time , is the number of links that end at link ’s beginning node, and is the number of links that begin at link ’s ending node.
Computing the interlink flows requires the use of two intermediate quantities for each link. These are the link demand , which is the number of vehicles that wish to exit link at timestep ; and the link supply, , which is the number of vehicles link can accept at time . Both and are functions of the density . The model that computes and from is often called the “fundamental diagram” or “link model,” and the model that computes the flows from all links’ supplies and demands is often called the “node model.”
A brief outline of how firstorder macroscopic simulation of a road network (sometimes called a dynamic network loading simulation) is performed could be:

At time , use the link model for each link to compute the link’s demand and supply as a function of its density .

Use the node model for each node to compute the interlink flows for all incoming links and outgoing links as functions of , , and information about vehicles’ desired movements .

Update the state of each link using (2.1).

Increment and repeat until the desired simulation end time is reached.
The tools described in this paper are compatible with any such link model. We make use of a particular node model that we have studied in Wright et al. (2016a, b). An aspect of this node model of particular relevance to managed lane modeling is our “relaxed firstinfirstout (FIFO) rule” construction (Wright et al., 2016a). This is necessary for modeling the flows between the GP and managed lanes. Without a FIFO relaxation, congestion in one of the two lane groups could block traffic in the other when that may be unrealistic (see Wright et al. (2016a, Section 2.2) for a detailed discussion).
So far, we have presented ingredients for a model that, while able to express many simple network topologies by joining links and nodes, does not capture several important behaviors in managed lanefreeway networks. The next three Sections briefly overview the additions to the standard model that will be explained in greater detail in the remainder of the paper.
2.2 Multiple classes of vehicles and drivers
In (2.1), we describe the number of vehicles in a link as a single number, . In this formulation, all vehicles are treated the same. However, for simulation in a managed lanefreeway network, it makes sense to break into different classes of vehicles and/or drivers. For example, for a freeway with an HOV lane facility, we might consider two classes: HOVs and nonHOVs. To this end, (2.1) can be rewritten as {linenomath}
(2.2) 
where indexes vehicle classes (often called “commodities” in the traffic literature).
Extending the density update equation to multiple classes means that the link and node models must also be extended to produce perclass flows . In this paper, we will not specify a particular link model, but assume use of one that produces perclass demands and overall supplies (the node model, in computing the , is responsible for splitting the available supply among the different demanding vehicle classes). Examples of this type of link model include those considered in Wong and Wong (2002), Daganzo (2002), and van Lint et al. (2008) (examples of multiclass link models of second or higherorder include those of Hoogendoorn and Bovy (2000). These types of models have higherorder analogs of supply and demand).
2.3 Toplogical expression of managed lanefreeway networks
In both (2.1) and (2.2), we describe a link in terms of its total density and its breakdown into percommodity portions, . By discretizing the road into these onedimensional links, we lose information about differences between vehicle proportions across lanes, as well as interlane and lanechanging behavior. This becomes a problem if such unmodeled behavior is of interest. In our setting, this means that modeling a freeway with a managed lane should not be done with a single link following (2.1) or (2.2), as it would be impossible to study the managed lanefreeway network behavior of interest.
Modeling differences in vehicle density across lanes is natural in microscopic and mesoscopic (see, for example, Treiber et al. (1999); Hoogendoorn and Bovy (1999); Ngoduy (2006), and others) models, but macroscopic models, in their simplicity, have less readilyaccessible avenues for including these differences. One straightforward method is to model each lane as a separate link, as in, e.g. Bliemer (2007) or Shiomi et al. (2015). However, this method has a few drawbacks. First, it requires the addition of some laneassignment method to prescribe the proportions of each vehicle class for each lane (such as a logit model as used in Farhi et al. (2013) and Shiomi et al. (2015)), which requires notalwaysaccessible data for calibration. Second, drastically increasing the number of links in a macroscopic model will necessarily increase the size of the state space and model complexity, which is, in a sense, incompatible with the overall goal of selecting a macroscopic model over a micro or mesoscopic model: some of the “macro” in the macroscopic model is lost.
Instead, in this paper we choose to model the GP lanes (or “GP lane group”) as one link and the parallel managed lanes (or “managed lane group”) as another link. Applied to an entire length of road, this creates a network topology of two “parallel chains” of links  one GP and one managed. The two chains will share nodes, but crossflows between the chains are permitted only in locations where there is physical access (i.e., no physical barriers) and policy access (i.e., no double solid lines under U.S. traffic markings). Where crossflows are possible, we do not use a logit model, but instead a driver behavior model first introduced in Wright et al. (2016a). This twochain model is similar to the one described in Liu et al. (2012), though in this reference, GPmanaged lane crossflows were not considered.
2.4 Observed phenomena in managed lanefreeway networks
The abovementioned friction effect is one of several emergent behaviors thought to be important for understanding traffic phenomena present on managed lanefreeway networks, and has been blamed for some of the “underperformance” of managed lanes (Jang and Cassidy, 2012). To expand on our earlier description, the friction effect describes a tendency of vehicles in an HOV lane to reduce their speed when the vehicles in the adjacent GP lane(s) congest to the point where they slow down. It has been hypothesized that this occurs due to HOV drivers being uncomfortable when traveling at a drastically higher speed than the adjacent GP vehicles, and slowing down to reduce the speed differential (Jang and Cassidy, 2012).
In Liu et al. (2012), the authors created a macroscopic model of a managed lanefreeway network with a friction effect model. They proposed modeling the friction effect by having two separate link models for the managed lanes – one each for when the GP lane is above or below some threshold density value.
In Section 3, we propose a more expressive model of the friction effect by replacing this piecewise implementation with a linear implementation, where the reduction in demands in the managed lane is proportional to the speed differential between the managed lane(s) and GP lane(s). A linear relationship like this was hypothesized in, e.g., Jang and Cassidy (2012, Fig. 5(c)). The magnitude of the friction effect (i.e., the degree to which managed lane drivers slow down to match GP lane drivers’ speed) is thought to be dependent on the type of separation between the managed lane(s) and GP lane(s) (e.g., painted lines vs. a concrete barrier) (Jang et al., 2012), which can be encoded by selecting the linear coefficient to this speed differential.
3 Full and GatedAccess Managed LaneFreeway Network Topologies
We will consider two types of managed lanefreeway network configurations: full access and separated with gated access. In a a fullaccess configuration, the managed lane(s) are not physically separated from the GP lane(s), and eligible vehicles may switch between the two lane groups at any location. Often, fullaccess managed lane(s) are specialuse only during certain periods of the day, and at other times they serve as GP lane(s) (e.g., HOV lanes are often accessible to nonHOVs outside of rush hour). On the other hand, in a gatedaccess configuration, traffic may switch between the managed lane(s) and GP lane(s) only at certain locations, called gates; at nongate locations, the two lane groups are separated by road markings (i.e., a double solid line in the U.S.) or a physical barrier. Usually, gatedaccess managed lanes are specialuse at all times. The implemented managed lane access scheme depends on jurisdiction. For example, fullaccess lanes are common in Northern California, and separated lanes are common in Southern California.
The differences in physical geometry and access points between the two access types requires two different types of topology in constructing a network for a macroscopic model.
3.1 Note on link and node models used in this section
As discussed in the previous Section, we attempt to be agnostic with regards to the particular link model (e.g., firstorder fundamental diagram) used in our implementations. However, in our model of the friction effect in Section 3.3, we parameterize friction being in effect on a particular link for a particular vehicle class at time by adjusting that link’s demand . For that discussion only, we specify a particular link model. This link model is reviewed in Appendix A. This link model is also the one used in the simulations presented in Section 6.
The node model used here and for the remainder of this paper, when a particular form is necessary, is the one discussed in Wright et al. (2016a) and Wright et al. (2016b). We use this node model because it handles multicommodity traffic, optimizes the utilization of downstream supply, makes use of input link priorities and has a relaxation of the “conservation of turning ratios” or “firstinfirstout” (FIFO) constraint of most node models. This last feature allows us to describe a set of GP lanes just upstream of an offramp with one link, and handle a condition of a congested offramp by having the congestion spill back onto only the offrampserving lanes of the GP link (as opposed to the entirety of the GP link). See Wright et al. (2016a, Section 3) and Wright et al. (2016b) for more discussion.
3.2 Fullaccess managed lanes
A fullaccess managed lane configuration is presented in Figure 1: GP and managed links (recall, as discussed above, all GP lanes and all managed lanes are collapsed into one link each) are parallel with the same geometry and share the same beginning and ending node pairs; traffic flow exchange between GP and managed lanes can happen at every node. Note that in Figure 1, we use a slightly irregular numbering scheme so that it is clear whether a link is a GP link, managed lane link, or ramp link. Parallel links in the graph have numbers made up of the digit of their terminating node, with GP links having one digit (i.e., link 1), managed lane links having two (i.e., link 11), and ramp links having three (i.e., link 111). Note also that we use a U.S.style, drivingontheright convention here, with the managed lane(s) on the left of the GP lanes and the ramps on the far right.
Links that are too long for modeling purposes (i.e., that create too lowresolution a model) may be broken up into smaller ones by creating more nodes, such as nodes 2 and 3 in Figure 1. Fundamental diagrams for parallel GP and managed lane links may be different (Liu et al., 2011).
We introduce two vehicle classes (): corresponds to the GPonly traffic and corresponds to the special traffic. When the managed lane(s) is (are) active, traffic is confined to the GP link, whereas traffic can use both the GP and managed lane links. We denote the portion of vehicles of class in link that will attempt to enter link as . This quantity is called the split ratio.
Split ratios for fullaccess managed lanes
We make an assumption that both vehicle classes take offramps at the same rate. For example, for node 1 in Figure 1, we might say that . Strictly speaking, it is not necessary to assume that the are equal for all . However, in practice the offramp split ratios are typically estimated from flow count data taken from detectors on the offramp and freeway. Generally speaking, these detectors cannot identify vehicle type, so the only quantity estimable is a flowweighted average of the quantities . To estimate the classspecific split ratios, one needs some extra knowledge of the tendency of each class to take each offramp (for example, that GPonly vehicles are half as likely as special traffic to take a certain offramp). Assuming that each class exits the freeway at the same rate is a simple and reasonableseeming assumption.
This same problem of unidentifiability from typical data appears in several other split ratios. First, it may not be possible to tell how many vehicles taking an offramp link come from the upstream GP, managed lane, or (if present) onramp link. In this case, some assumptions must then be made. For example, three different assumptions that may be reasonable are (1) that vehicles in each link take the offramp at the same rate; or (2) that no vehicles from the managed lane(s) are able to cross the GP lanes to take the offramp at this node, and that no vehicles entering via the onramp, if one is present, exit via the offramp at the same node; or (3) that vehicles in GP and managed links take the offramp at the same rate, while no vehicles coming from the onramp are directed to the offramp. Looking back at node 1 in Figure 1 again, assumption (1) would say the are equal for all ; assumption (2) would say ; and assumption (3) would say and . The best assumption for each node will depend on the road geometry for that particular part of the road (how near the offramp is to any onramps, how many GP lanes a vehicle in a managed lane would have to cross, etc.).
Second, the crossflows between the GP and managed lane links are not observable. Even if there exist detectors immediately upstream and downstream of the node where traffic can switch between GP and managed lanes, it is impossible to uniquely identify the crossflows. In a simulation, these crossflows must be governed by some driver choice model.
Putting together these assumptions and the specialtrafficonly policy for the managed lane, we can summarize most of the necessary split ratios needed for computing flows in a node model. For example, for node 1 in Figure 1, {linenomath}
where “n/a” means that the split ratios are not applicable, as there should be no vehicles of class in the managed lane. The split ratios marked with a dash are those abovementioned flows that are unobservable and come from some driver choice model. Of course, whatever method is chosen to compute these unknown split ratios, we must have .
Similarly, for node 2, which does not have an onramp or an offramp, {linenomath}
with “n/a” and the dash meaning the same as above.
As previously mentioned, fullaccess managed lanes often have certain time periods during which nonspecial () traffic is allowed into the managed lane. We can model this change in policy simply by changing the split ratios at the nodes. For node 1, for example, the nonrestrictive policy is encoded as {linenomath}
and for node 2 as
In other words, the managed lane link is treated as additional GP lane(s), and the split ratios governing the crossflows between the two links should be found from the driver choice model for both vehicle classes.
Node model for fullaccess managed lanefreeway networks
As mentioned in Section 3.1, we make use of the node model discussed in Wright et al. (2016a) and Wright et al. (2016b) to describe a freeway network with managed lanes. This node model differentiates itself from others in that it deals with multicommodity traffic flow, optimally utilizes the available supply, makes use of input link priorities, and has a relaxation of the common FIFO constraint. By default, link priorities can be taken proportional to link capacities. To explain relaxed FIFO, say that some link has vehicles that wish to enter both links and . If link is jammed and cannot accept any more vehicles, a strict FIFO constraint would say that the vehicles in that wish to enter will queue at ’s exit, and block the vehicles that wish to enter . In a multilane road, however, only certain lanes may queue, and traffic to may still pass through other lanes. The relaxation is encoded in socalled “mutual restriction intervals” . This interval partly describes the overlapping regions of link ’s exit that serve both links and . For , a portion of ’s lanes that serve also serve , and will be blocked by the cars queueing to enter when is congested. For example, if is served by three lanes of , and of those three lanes, the leftmost also serves , we would have .
As an example, we consider again node 1 in Figure 1. Say that the GP links (1 and 2) have four lanes, that the managed lane links (11 and 22) have two lanes, and that the onramp merges into and the offramp diverges from the rightmost GP lane. Further, we say that when the managed lane link is congested, vehicles in the GP lanes that wish to enter the managed lanes will queue only in the leftmost GP lane. On the other hand, when the GP link is congested, vehicles in the managed lanes that wish to enter the GP lanes will queue only in the rightmost managed lane. Finally, we suppose that jammed offramp (222) will cause vehicles to queue only in the rightmost GP lane. Taking together all of these statements, our mutual restriction intervals for this example are: {linenomath}
To read the above tables, recall that as written, is the congested, restricting link, and is the restricted link.
These chosen restriction intervals allow for expected behavior in this network, such as a congested GP link causing possible queueing in the right managed lane (if some drivers are trying to enter the GP link), but no spillback into the left managed lane.
3.3 Friction effect
The friction effect is an empiricallyobserved phenomenon in situations where managed lanes are relatively uncongested, but the managedlane traffic will still slow down when the adjacent GP lanes congest and slow down (see Daganzo and Cassidy (2008), Liu et al. (2011), Cassidy et al. (2015), etc.). It has been hypothesized (Jang and Cassidy, 2012) that this phenomenon arises from the managedlane drivers’ fear that slowermoving vehicles will suddenly and dangerously change into the managed lane ahead of them.
We suggest modeling the friction effect based on a feedback mechanism that uses the difference of speeds in the parallel GP and managed lane links to scale down the flow (and therefore the speed) out of the managed lane link if necessary.
To explain the concept, we again refer to Figure 1 and consider parallel links 1 (GP) and 11 (managed lane). Recall that, under a firstorder model (2.2), the speed of traffic in link at time is {linenomath}
(3.1) 
where is the theoretical free flow speed of link at time .
We say that the friction effect is present in managed lane link 11 (following the notation of Figure 1) at time if {linenomath}
(3.2) 
which means that (1) the GP link is in congestion (its speed is below its current free flow speed), and (2) the speed in the GP link is less than the speed in the managed lane link. We denote this speed differential as: {linenomath}
(3.3) 
It has been observed (Jang et al., 2012) that the magnitude of the friction effect — the degree to which managedlane drivers slow down towards the GP lane’s traffic speed — depends on the physical configuration of the road. For example, less of a friction effect will be present on managed lanes that are separated from the GP lanes by a buffer zone than those that are contiguous with the GP lanes (Jang et al., 2012), and the presence of a concrete barrier would practically eliminate the friction effect. Other factors that may affect this magnitude include, for example, whether there is more than one managed lane, or whether there is a shoulder lane to the left of the managed lane that drivers could swerve into if necessary.
To encode this variability in the magnitude of the friction effect in managed lane link 11, we introduce the friction coefficient of this link. The friction coefficient reflects the strength of the friction. Its value depends on the particular managed configuration and is chosen by the modeler. A value of means there is no friction (which may be appropriate if, perhaps, the managed lane(s) are separated from the GP lanes by a concrete barrier), and means that the managed lane link speed tracks the GP link speed exactly.
When the friction effect is active (i.e., when (3.2) is true), we adjust the fundamental diagram of the managed lane link by scaling down its theoretical free flow speed , and propagate that change through the rest of the fundamental diagram parameters. The exact mathematical changes will of course be different for every different form of fundamental diagram. For the particular fundamental diagram discussed in Appendix A, this means adjusting the free flow speed and capacity as follows: {linenomath}
(3.4)  
(3.5) 
where is the high critical density (see Appendix A for its definition), and using these adjusted values in the calculation of the sending function (A.1), {linenomath}
(3.6) 
For the fundamental diagram of Appendix A, we must also check whether {linenomath}
(3.7) 
If not, then applying friction will lead to the managed lane link speed falling below the GP link speed, and the unadjusted sending function should be used (the possibility of this happening is due to the use of two critical densities to create Appendix A’s fundamental diagram’s “backwards lambda” shape).
3.4 Separated managed lanes with gated access
A separated, gatedaccess managed lane configuration is presented in Figure 2. Unlike the fullaccess configuration, the GP and managed lane link chains do not necessarily meet at every node. Instead, they need only meet at a few locations, where vehicles can move into and out of the managed lane(s). Note that, unlike the fullaccess configuration, there is no need for GP and managed lane links to be aligned.
As labeled in Figure 2, the nodes where the two link chains meet are called gates (as an aside, one way to describe the fullaccess managed lane configuration would be that every node is a gate). Similar to our construction of excluding GPonly traffic from the managed lane in the fullaccess case, we can disable flow exchange at a given gate by fixing split ratios so that they keep traffic in their lanes. For example, to disable the gate (the flow exchange between the two lanes) at node 2 in Figure 1, we set and (this means that and ), . Thus, the fullaccess managed lane can be easily converted into the separated managed lane by setting nonexchanging split ratios everywhere but designated gatenodes.
In practice, a gate is stretch of freeway that may be a few hundreds of meters long (Cassidy et al., 2015), and, potentially, we can designate two or three sequential nodes as gates. In this paper, however, we model a gate as a single node.
For the gatedaccess configuration, we suggest setting mutual restriction coefficients in the same manner as fullaccess managed lanes, in Section 3.2.2.
Compared to the fullaccess managed lane configuration, the gatedaccess configuration has a smaller friction effect (Jang et al., 2012): drivers in the separated managed lane feel somewhat protected by the buffer, whether it is virtual (double solid line) or real (concrete), from vehicles changing abruptly from the slow moving GP lane and, therefore, do not drop speed as dramatically. The degree to which the friction effect is mitigated is disputed (e.g., see Footnote 3 in Cassidy et al. (2015)), but overall the bottlenecks created by the gates are much greater instigators of congestion (Cassidy et al., 2015). Inclusion of the friction effect in modeling separated managed lane configurations is thus not as essential as in modeling the fullaccess case.
Modeling a flow of vehicles from the managed lanes to the offramps
Recall from Section 3.2.1 that, in the fullaccess managed lane model, we can model vehicles moving from the managed lane link to offramps in a straightforward manner, by setting corresponding split ratios (for example, , , for node 1 in configuration from Figure 1). For the gatedaccess configuration, however, modeling traffic as moving from the managed lanes to offramps is more complicated: generally, gates do not coincide with offramp locations. In fact, there are typically between two and five offramps between two gate locations. We denote the offramps in the GP road segment connecting two gates as exits (see Figure 2). These offramps cannot be accessed directly from the managed lane. Instead, vehicles traveling in the managed lane link that intend to take one of the exits , must switch from the managed lane link to the GP link at gatenode 1 and then be directed to the correct offramp. Creating a simulation where vehicles behave like this requires a more involved modeling construction.
To resolve this challenge, our gatedaccess model introduces new vehicle classes in addition to the (GPonly) and (special) traffic used in the fullaccess model of Section 3.2. These additional classes will be used to distinguish subsets of the special traffic population by its destination offramp. If is the largest number of offramps between two adjacent gates, then altogether we have vehicle classes: , where indicates the class of vehicles that will exit through the th offramp after leaving the managed lane through the gate. By definition, traffic of type may exist in the GP lane segment between gate 1 and offramp , but there is no traffic of this type either in the GP link segment between offramp and gate 2 or in the managed lane link. This movement pattern is ensured by setting constant split ratios: {linenomath}
(3.8) 
where , and denotes the input GP link for the node that has the output link (see Figure 2).
Vehicles of class do not enter the network via onramps or the upstream boundary, but instead are converted from special traffic as it leaves the managed lane(s) through the gate. We perform this conversion as part of the link model computation, such that the total demand for the switching link, remains the same before and after the switch. The exact link on which this switching takes place is the managed lane link immediately upstream of the gate (e.g., link 11 in Figure 2).
We say that the amount of traffic that should change from vehicle class to vehicle class at time is (using link 11 as an example) , where is in units of vehicles per simulation timestep (this factor is included so that the switching done is proportional to link 11’s outflow, rather than its density), and is the split ratio from , the GP link immediately upstream of exit , and the exit at time . This statement is based on the assumption that the vehicles in the managed lane link will exit the freeway through exit at the same rate as special () vehicles that happened to stay in the GP lanes. That is, if a portion of vehicles intend to leave the GP lanes through exit , then a portion of the vehicles in the managed lane(s) will leave the managed lane link at the closest upstream gate and leave the network at exit when they reach it. Note that if there are offramps between two particular gates, then no vehicles should switch to type at the upstream gate, as they would have no ramp to exit through.
Using Figure 2 as a reference, we can now formally describe the procedure for destination assignment to traffic in the managed lane link.

Given are vehicle counts per commodity , ; free flow speed ; and offramp split ratios and , .
^{2} 
Initialize:

Assign type traffic:
(3.9) (3.10) (3.11) 
If , then set and return to step 3.

Update the state:
After the switches to class traffic have been done, there may be unresolved split ratios for both classes and at the gates (similar to the dashed split ratios in the tables in Section 3.2.1). These split ratios should be filled in with the same tools as those in Section 3.2.1: some sort of driver lane choice behavior.
3.5 Inertia effect
At the end of the previous Section, we mentioned that some split ratios will likely be undefined after switching class vehicles to the classes. In particular, class vehicles in the upstream GP link (e.g., link 1 in Figure 2) and remaining vehicles in the upstream managed lane link (e.g. link 11 in Figure 2) will need to decide whether to pass through the gate or remain in their current lane. Sample tools for modeling driver lane choices such as these include the class of “logit” logistic regression models (McFadden, 1973), or dynamic split ratio solvers such as the one presented in Wright et al. (2016a) and reviewed in Appendix B.
When applied to lane choice (e.g., as in Farhi et al. (2013)), logit models produce a set of portions in , one for each lane, that sum to one. Each lane’s value is the equilibrium portion of vehicles that will travel on that lane. In a dynamic simulation context, such as considered in this article, the differences between the logit model’s equilibrium portions and the current distribution of vehicles across lanes at time are used to select split ratios at time such that the actual distribution approaches the logit equilibrium. We argue, though, that use of this sort of split ratio solver in unmodified form might be inappropriate for computing gate split ratios in the gatedaccess managed lane configuration. Unmodified, a value function for either the GP or managed lane link might include terms such as the link’s speed of traffic, density, etc. However, often a gatedaccess managed lane is separated from the GP lanes by some buffer zone or visibilityobstructing barrier that makes switching between the two links more hazardous than switching between two contiguous lanes. Therefore, if one uses a logitbased model, it would be appropriate to modify the logit model’s value function such that staying in the current link (e.g., the movements (1,2) and (11,22) in Figure 3) has some positive value for drivers, and the gains (e.g., in travel time) for ingress/egress movements (e.g., (1,22) and (11,2) in Figure 3) must be of more value than the stayinginthelane value. We refer to this model as the inertia effect.
We may also incorporate the inertia effect into dynamic split ratio solvers, such as the one introduced in Wright et al. (2016a) and reviewed in Appendix B. At time , this particular algorithm selects split ratios in an attempt to balance the density ratio at the next timestep , , where is the jam density, or the maximum number of vehicles that link can hold. For example, if applied to the node in Figure 3, this algorithm would attempt to make and as equal as possible. For more details, see Appendix B and Wright et al. (2016a).
Here, we modify several steps of the solver so that this equalityseeking goal is balanced with a goal towards enforcing the inertia effect. We illustrate these changes with the particular example of the node in Figure 3. Ensuring that the split ratio assignment algorithm gives preferences to movements (1,2) and (11,22) over (1,22) and (11,2) can be done in step 5 of the original algorithm, setting of oriented priorities. Specifically, we modify (B.4). For this particular example, the original formula gives us: {linenomath}
Since for , , , we get and , which, according to (B.3), yields and .
For each input link that forms one lane with an output link and class , such that , (B.4) can be modified as follows: {linenomath}
(3.12) 
where the parameter is called the inertia coefficient, and indicates how strong the inertia effect is. With , (3.12) reduces to the original formula, (B.4). With , all the a priori unassigned traffic from link must stay in its lane — be directed to output link . The choice of lies with the modeler. In the case of example from Figure 3, the modified formula (3.12) yields: {linenomath}
(3.13) 
where , and picking () would give preference to movement (1,2) over (1,22) (and (11,22) over (11,2)).
The way of choosing for multiple input links is not obvious and an arbitrary choice may result in an unbalanced flow distribution among output links. Therefore, we suggest picking just one inputoutput pair , and for that input link setting , while for other input links setting , . The input link must be from the lane that is expected to have a positive net inflow of vehicles as a result of the split ratio assignment and flows computed by the node model. So, {linenomath}
(3.14) 
where {linenomath}
(3.15) 
and denotes the output link that is in the same lane as input link .
For the particular example node in Figure 3, we need to determine whether flow from link 1 will proceed to link 2 or flow from link 11 to link 22, while other a priori undefined split ratios will be computed according to the split ratio assignment algorithm. If , then , , and a priori unassigned traffic in the managed lane will stay in the managed lane (), while the a priori unassigned traffic coming from links 1 and 111, will be distributed between links 2 and 22 according to the dynamic split ratio solver.
On the other hand, if , then , , and the a priori unassigned traffic in the GP lane will stay in the GP lane (), while the a priori unassigned traffic coming from links 11 and 111, will be distributed between links 2 and 22 according to the dynamic split ratio solver.
4 Putting Together the Pieces: A Managed LaneFreeway Network Simulation Model
In this Section, we present a unified simulation algorithm for managed lanefreeway networks. This can be considered a fleshingout of the simplified firstorder macroscopic simulation method briefly outlined in Section 2.1, with extensions made by incorporating the additional items we have described in the sections between that and this one.
4.1 Definitions

We have a network consisting of set of links and a set of nodes .

A node always has at least one incoming link and one outgoing link.

A link may have an upstream node, a downstream node, or both.


We have different vehicle classes traveling in the network, with classes indexed by .

Let denote the simulation timestep.

In addition, while we have not covered them here, the modeler may optionally choose to define control inputs to the simulated system that modify system parameters or the system state. Such control inputs may represent operational traffic control schemes such as ramp metering, changeable message signs, a variable managedlane policy, etc. In the context of this paper, we suggest including the parameterization of the friction effect and the classswitching construction of the separatedaccess managed lane model as control actions.
Link model definitions

For each link , let there be a timevarying dimensional state vector , which denotes the density of the link of each of the vehicle classes at time .

Each element of this vector, , updates between timestep and timestep according to (4.1).

Also define for all , the initial condition of the system.


We define three types of links:

Ordinary links are those links that have both beginning and ending nodes.

Origin links are those links that have only an ending node. These links represent the roads that vehicles use to enter the network.

Destination links are those links that have only a beginning node. These links represent the roads that vehicles use to exit the network.


For each origin link , define a dimensional timevarying vector, , where the th element denotes the exogenous demand of class into the network at link .
Node model definitions

For each node , let denote the incoming links and denote the outgoing links.

For each node, define the timevarying split ratios for each triplet . Each split ratio may be fully defined, partially defined, or fully undefined. If some split ratios for a node are undefined, then also define for that node some split ratio solver (e.g., logit or dynamic) to fill in undefined split ratios at time .

For each node, define a “node model” that, at each time , takes its incoming links’ demands and split ratios , its outgoing links’ supplies , and other nodal parameters, and computes the flows . In this paper, we refer to a specific node model with a relaxed firstinfirstout (FIFO) construction that has additional parameters (mutual restriction intervals) and (the incoming links’ priorities).
State update equation definitions

All links update their states according to the equation {linenomath}
(4.1) which is a slightly generalized form of (2.2).

For all ordinary and destination links, {linenomath}
(4.2) where is the beginning node of link .

For all origin links, {linenomath}
(4.3) 
For all ordinary and origin links, {linenomath}
(4.4) where is the ending node of link .

For all destination links, {linenomath}
(4.5)
4.2 Simulation algorithm

Initialization: {linenomath}
for all , .

Perform all control inputs that have been (optionally) specified by the modeler. For our purposes, this includes:

For each managed lane link whose downstream node is a gate node in a gatedaccess managed lane configuration, perform class switching as detailed in Section 3.4.1

For each link and commodity , compute the demand, using the link’s link model.

For each ordinary and destination link , compute the supply using the link model. For origin links, the supply is not used.

For each node that has one or more undefined split ratios , use the node’s split ratio solver to complete a fullydefined set of split ratios. Note that if an inertia effect model is being used, the modified split ratio solver, e.g. the one described in Section 3.5, should be used where appropriate.

For each node , use the node model to compute throughflows for all .

If , then stop. Otherwise, increment and return to step 2.
5 Calibrating the Managed LaneFreeway Network Model
Typically, a traffic modeler will have some set of data collected from traffic detectors (e.g., velocity and flow readings), and will create a network topology with parameter values that allow the model to reproduce these values in simulation. Then, the parameters can be tweaked to perform prediction and analysis. For our managed lanefreeway networks, the parameters of interest are:

Fundamental diagram parameters for each link. Calibration of a fundamental diagram is typically agnostic to the node model and network topology, and there exists an abundant literature on this topic. For the purposes of the simulations in the following Section, we used the method of Dervisoglu et al. (2009), but any other method is appropriate.

Percentage of special (that is, able to access the managed lane) vehicles in the traffic flow entering the system. This parameter depends on, e.g., the time of day and location as well as on the type of managed lane. It could be roughly estimated as a ratio of the managed lane vehicle count to the total freeway vehicle count during periods of congestion at any given location.

Inertia coefficients. These parameters affect only how traffic of different classes mixes in different links, but they have no effect on the total vehicle counts produced by the simulation.

Friction coefficients. How to tune these parameters is an open question. In Jang and Cassidy (2012) the dependency of a managed lane’s speed on the GP lane speed was investigated under different densities of the managed lane, and the presented data suggests that although the correlation between the two speeds exists, it is not overwhelmingly strong, below 0.4. Therefore, we suggest setting friction coefficients to values not exceeding .

Mutual restriction intervals. It is also an open question how to estimate mutual restriction intervals from the measurement data. See the discussion in Section 3.2.2 for some guidelines.

Offramp split ratios.
Calibrating a traffic model, or identifying the best values of its parameters to match realworld data, is typically an involved process for all but the simplest network topologies. In particular, once we consider more than a single, unbroken stretch of freeway, the nonlinear nature and network effects of these systems means that estimating each parameter in isolation might lead to unpredictable behavior. Instead, nonlinear and/or nonconvex optimization techniques such as genetic algorithms (Poole and Kotsialos, 2012), particle swarm methods (Poole and Kotsialos, 2016), and others (Ngoduy and Maher (2012); Fransson and Sandin (2012), etc.) are employed.
In the managed lanefreeway networks we have discussed, the key unknown parameters we have introduced are the offramp split ratios, item 6, which may be particularly hard to estimate as they are typically timevarying and explicitly represent driver behavior, rather than physical parameters of the road. The other items are not too difficult to identify using methods from the literature. The remainder of this section describes iterative methods for identification of the offramp split ratios for both the fullaccess and gatedaccess configurations.
5.1 Split ratios for a fullaccess managed lane
Consider a node, one of whose output links is an offramp, as depicted in Figure 3. We shall make the following assumptions.

The total flow entering the offramp, , at any given time is known (from measurements) and is not restricted by the offramp supply: .

The portions of traffic sent to the offramp from the managed lane and from the GP lane at any given time are equal: , .

None of the flow coming from the onramp (link 111), if such flow exists, is directed toward the offramp. In other words, , .

The distribution of flow portions not directed to the offramp between the managed lane and the GP output links is known. This can be written as: , where , as well as , , , , are known.

The demand , , , and supply , , are given.
At any given time, is unknown and is to be found.
If were known, the node model would compute the inputoutput flows, in particular, , . Define {linenomath}
(5.1) 
Our goal is to find from the equation
(5.2) 
such that , where . Obviously, if , the solution does not exist, and the best we can do in this case to match is to set , directing all traffic from links 1 and 11 to the offramp.
Suppose now that . For any given , we assume is a monotonically increasing function of (this assumption is true for the particular node model of Wright et al. (2016a)). Moreover, , while . Thus, the solution of (5.2) within the given interval exists and can be obtained using the bisection method.
The algorithm for finding follows.

Initialize:

If , then are not enough vehicles to satisfy the offramp demand. Set and stop.

Use the node model with and evaluate . If , then set and stop.

Use the node model with and evaluate . If , then set and stop.

If , then update:
Else, update:

Set and return to step 4.
Here, represents the lower bound of the search interval at iteration and the upper bound.
5.2 Split ratios for the separated managed lane with gated access
The configuration of a node with an offramp as one of the output links is simpler in the case of a separated HOV lane, as shown in Figure 4. Here, traffic cannot directly go from the HOV lane to link 222, and, thus, we have to deal only with the 2input2output node. There is a caveat, however. Recall from Section 3.4 that in the separated managed lane case we have destinationbased traffic classes, and split ratios for destinationbased traffic are fixed.
We shall make the following assumptions:

The total flow entering the offramp, , at any given time is known (from measurements) and is not restricted by the offramp supply: .

All the flow coming from the onramp (link 111), if such flow exists, is directed toward the GP link 2. In other words, and , .

The demand , , , and supply are given.

We denote the set of destinationbased classes as . The split ratios for are known. Let the split ratios for , where is to be determined (i.e., we assume all nondestinationbased classes exit at the same rate).
The first three assumptions here reproduce assumptions 1, 3 and 5 made for the fullaccess managed lane case. Assumption 4 is a reminder that there is a portion of traffic flow that we cannot direct to or away from the offramp, but we have to account for it.
Similarly to the fullaccess managed lane case, we define the function : {linenomath}
(5.3) 
where , are determined by the node model. The first term of the righthand