Optimum Location of Autonomous Vehicle Lanes: A Model Considering Capacity Variation

,is paper proposes a model to find the optimal location of autonomous vehicle lanes in a transportation network consisting of both Autonomous Vehicles (AVs) and Human-Driven Vehicles (HDVs) while accounting for the roadway capacity variation.,e main contribution of the model is considering a generalized definition of capacity as a function of AV proportion on a link and incorporating it into the network design problem. A bilevel optimization model is proposed with total travel time as the objective function to be minimized. At the upper-level problem, the optimal locations of AV lanes are determined, and at the lower level which is a multiclass equilibrium assignment, road users including both AVs and HDVs seek to minimize their individual travel times. It is shown that if capacity variation is ignored, the effect of AV lane deployment can be misleading. Since there will be a long transition period during which both AVs and HDVs will coexist in the network, this model can help the networkmanagers to optimally reallocate the valuable road space and better understand the effects of AV lane deployment at the planning horizon as well as during the transition period. Employing this model as a planning tool presents how the proposed AV lane deployment plan could consider the AV market penetration growth during the transition period. Numerical analysis based on the Sioux Falls network is presented in two cases with and without variable capacity to illustrate the application of this model. At the 60% penetration rate of AVs, the improvement in total travel time was 3.85%with a fix capacity while this improvement was 9.88%with a variable capacity.


Introduction
In recent years, Autonomous Vehicles (AVs) have received considerable attention from different research groups all over the world. ey are expected to have potential benefits on safety, mobility, traffic efficiency, driver productivity, road capacity, and environment [1][2][3][4].
Despite significant advances in autonomous vehicle technology, there are still many stages and many years left to make driverless cars ready for broad adoption [5]. Considering the patterns of previous vehicle technology deployments such as automatic transmissions, airbags, and hybrid vehicles, it takes decades for autonomous vehicle technology to reach market saturation [5]. Reviewing studies about AV adoption [5][6][7][8] leads to the conclusion that there will be a long transition period in which AVs and human-driven vehicles (HDVs) share the road. erefore, policymakers and practitioners should prepare strategies for the time that AVs and HDVs coexist and share the same transportation network. ese strategies need to recognize different sources of uncertainty in transportation planning for AVs such as technology, adoption rates, roadway capacity, land use changes, and travel demand [9].
Due to the high uncertainty in the penetration rate of AVs, flexible and low-cost strategies are required. Among available strategies, providing dedicated lanes to autonomous vehicles using existing infrastructure is a flexible and low-cost one. It is a tactical decision which can benefit the network through existing resources [10]. However, conversion of a general purpose lane to a dedicated lane may lead to public dissatisfaction. A dedicated AV lane is one of the managed lane strategies with a similar concept to High Occupancy Vehicle (HOV) lanes, High Occupancy Toll (HOT) lanes, and express lanes. Previous experiences with these strategies show if they are designed and operated properly, they can improve safety, efficiency, and the environment quality [11]. Despite the similarities between dedicated AV lanes and other managed lane strategies, there are yet differences such as interactions between AVs and HDVs, which make it necessary to specifically examine the adaptation of dedicated AV lanes.
Some studies suggest that many potential AV benefits require platooning, which can be achieved by dedicating lanes to AVs [5]. On the other hand, effective platooning and consequently successful implementation of dedicated AV lanes highly depends on their market penetration rate; for example, at low market penetration rates benefits may not be significant [12]. Several studies have used simulation models to investigate the effect of a dedicated AV lane on traffic flow. Microscopic simulation models are one of the tools to consider the interactions between AVs and HDVs. Talebpour et al. [12] developed a microscopic simulation model to explore how congestion and travel time reliability are affected by reserving a lane for AVs on a hypothetical twolane segment and a four-lane highway segment. eir results indicate that optional use of dedicated lanes by AVs can improve link performance both in terms of congestion and travel time reliability while limiting AVs to dedicated lanes decreases the overall performance. Vander Laan and Sadabadi [13] used macroscopic traffic simulation to address effects of a dedicated AV lane on the performance of a congested corridor. ey applied it to a four-lane corridor with one AV lane. eir results showed that the best system performance is obtained at medium market penetration rates (i.e., 30%, 40%, and 50%). Ye and Yamamoto [14] reached a similar conclusion. ey proposed a microscopic simulation model to evaluate the performance of overall traffic throughput under different AV market penetration rates while allocating a lane to connected and autonomous vehicles (CAVs) in a heterogeneous traffic flow. ey found that the positive effects of CAV dedicated lanes emerge at medium market penetration rates. Ghiasi et al. [15] formulated the capacity of a one-lane highway shared by CAVs and HDVs considering stochastic headways between different vehicle types. ey also presented a lane management model that determines the optimal number of CAV dedicated lanes on a highway segment while maximizing throughput.
While most studies have explored the effect of a dedicated AV lane in a highway segment, there are only a limited number of studies investigating these effects at the network level. To understand the effects of exclusive AV lanes at the network level, Ivanchev et al. [16] designed an agent-based simulation for the city of Singapore. ey converted one regular lane of each highway segment to an exclusive AV lane. e results indicate that the introduction of AV lanes to all highway segments worsens the entire network performance compared to the network without AV lanes for all AV penetration rates. e reason may be that introducing AV lanes to all highway segments limits the access of HDVs to a large part of the network. erefore, the optimal location of AV lanes should be determined in order to obtain benefits at the network level. Wu et al. [17] examined network flow pattern under user equilibrium (UE) and system optimum (SO) in a linear traffic corridor with dedicated AV expressways and nonautonomous local streets. ey also developed a problem of AV expressway exit as an optimization bilevel problem and explored the effect of exit locations on UE and SO patterns. Chen et al. [18] developed a mathematical framework which determined the optimal location, the time to deploy, and the number of AV lanes in a transportation network. ey considered the capacity increase associated with AV shorter headways in an exclusive AV lane; however, capacity impacts of the presence of AVs in the traffic flow have been ignored. In fact they did not consider that the capacity of a roadway shared by AVs and HDVs is different from the capacity of the same roadway with homogeneous traffic flow of AVs or HDVs. e approach used in previous studies [18] may lead to inaccurate conclusions as we show in this study. To avoid such inexact results, we have introduced the capacity function of mixed traffic to the lower level of the network design problem investigated in this paper.
To bridge this research gap, this study aims to develop a mathematical model that determines the optimal location of AV lanes in a network accounting for capacity changes of a roadway shared by AVs and HDVs. We considered capacity of each link as a function of link autonomy level which is the proportion of AVs on a link to total link flow. We call this approach variable capacity case, which can be compared against the fixed capacity case (i.e., the problem solved by [18]). is problem is designed as a network design problem in which at the upper level, the location of AV lanes are determined. e lower-level problem is a multiclass equilibrium assignment in which both AVs and HDVs follow the Wardrop equilibrium principle to minimize their individual travel costs. No study to date has formulated a network design problem for a network with mixed traffic of AVs and HDVs while examining the effects of capacity variations associated with presence of AVs. e main advantage of the proposed model is considering a more realistic capacity function which is already developed in previous link level studies [15,19] to the network level study by applying it in traffic assignment and the network design problem. e proposed change improves the estimated travel time accuracy as well as link flows in a network. To better explain the concept and contribution of this study, an analysis on a small network of Figure 1 is provided below.
Consider the network of Figure 1. Let w be the set of origin-destination (OD) pairs. ere are two OD pairs, where t a is the link travel time, t 0 a is the link free-flow travel time which is assumed to be 6 minutes for all links of this example, and V a and C a are link flow and capacity, respectively. We computed total travel time at equilibrium for two cases: the fixed capacity case and the variable capacity case (i.e., the approach of this study). Per-lane capacity for the fixed capacity case is assumed to be 2400 vehicles and for the variable capacity case is computed by equation (2) [16], where h AV and h HDV are headway of AVs and HDVs to a leading vehicle, respectively, and p a is proportion of AVs on a link: Let h AV � 0.85 and h HDV � 1.5, the total travel time is equal to 15,231 vehicle minutes for the first case and 61,341 for the second case.
e results indicate that integrating capacity variations into the link performance function improves the accuracy of travel time estimation. Since traffic assignment is the lower-level problem of a network design problem, the travel time estimation affects the optimal solution too. e remainder of this paper is divided into five sections. Section 2 provides a brief review of studies related to our work. Section 3 proposes the mathematical framework, and Section 4 explains the solution method used to solve the problem. In Section 5, results of numerical analysis are presented. Section 6 concludes the paper.

Literature Review
In this section, studies related to this work is summarized in three main topics: (a) network design problems solved for AVs so far, (b) equilibrium assignment problem in transportation networks with AVs and HDVs, and (c) capacity impacts of AVs in a heterogeneous traffic flow (mix of AVs and HDVs). Little research has been conducted on the first and second topics. However, the last issue has drawn the main attention from researchers.

Network Design Problem (NDP) in the Context of Autonomous
Vehicles. NDP is a well-known optimization problem in transportation planning literature, which aims to find the optimal value of specific objectives (e.g., total travel time) under limited resources. e ultimate goal of the NDP is to determine optimal projects among a set of alternatives while accounting for route choice behavior of road users. NDP can be considered as a two-stage leaderfollower Stackelberg game which is usually formulated as a bilevel problem [20]. At the upper level, the network manager (leader) aims to optimize network performance under resource constraints, and at the lower level, the network users (followers) choose a route with the minimum travel cost.
Research on a network design problem for different transportation modes with various objective functions, decision variables, and a wide range of solution methods has a long history in transportation [21,22]. However, a few network design problems have considered the context of AVs.
Chen et al. [18] proposed an optimization model for AV lane deployment on a transportation network in a mixed AV and HDV environment. e upper-level problem determined the deployment plan of AV lanes to minimize the social cost while considering a diffusion model to forecast the AV adoption rate. e lower-level problem was a multiclass user equilibrium which followed Wardrop's first equilibrium principle [23]. e objective function was the total social cost, which included the value of travel time for both modes and the cost related to the unsafety factor for using HDVs in the planning horizon. ey used the US Bureau of Public Roads (BPR) formula as a volume-delay function while assuming a fix per-lane capacity for regular links in the network. e per-lane capacity for AV links was assumed fixed and greater than the per-lane capacity of the adjacent regular link, given that AVs can achieve shorter headways compared to HDVs. Applying the proposed model to the South Florida network showed that there was no considerable gain from deploying AV lanes at low market penetrations (e.g., less than 20%). Furthermore, the optimal deployment plan was at the 25th year of a 40 year horizon where the social cost was reduced by 3.91%.
Another related network design problem is the optimal design of AV zones in transportation networks [24]. e problem was expressed as a mixed integer bilevel model that determined the optimal AV zones in a general network while both AVs and HDVs followed the user equilibrium outside the AV zones. AVs assumed to follow the system optimum within the AV zones. In their study, the AV zone is defined as a zone with links dedicated to AVs. Similar to their previous study, it was assumed that the capacity of links within an AV zone was more than mixed AV and HDV links and the capacity of links with mixed AV and HDV and HDV only links were equal. Numerical analyses showed that designing AV zones in a transportation network reduced the total travel cost. Mathematical Problems in Engineering

User Equilibrium Assignment Problem in Mixed Traffic.
Static equilibrium assignment is the most common problem used as the lower level of the NDP. Although many researchers have studied the problem properties, formulation, extensions, and solution methods since the introduction of Wardrop principles in 1952 [23], there is little literature on the static equilibrium assignment with mixed autonomy to the best of the authors' knowledge. Of course, the main reason is the emerging nature of the subject. It will be of particular importance to formulate the equilibrium assignment problem appropriately when solving an NDP where AVs and HDVs share same roads, and all network users have a selfish route behavior. e mixed autonomy traffic assignment is a multiclass equilibrium model with special properties. e multiclass user model, which was formulated by Dafermos [25] is more complex than a single user model. e single user model can be formulated as a convex optimization model with a unique solution under the assumption of separable and increasing link performance function [26]. However, there exists no convex formulation for a general network with multiple classes of vehicles. In addition, class-specific link flow solution is not unique for a general network. However, the total flow of each link is unique in certain cases. For instance, with a BPR volume-delay function, fixed capacity for all links, and same capacity for all vehicle classes, uniqueness of total link flows is guaranteed [27]. Levin and Boyles [28] is the first study that accounted for roadway capacity changes in the presence of AVs in a traffic assignment model. ey developed a four-stage travel demand model to investigate the effects of AV ownership on travel behavior. ey used a static assignment model to identify the route choice behavior of multiple classes of road users, which are divided according to their value of time. e travel time of their model was based on the BPR function with capacity specified as a linear function of the jam density as defined in Greenshield's model (equation (3). In this study, jam density is defined by equation (4) as a function of the proportion of AVs on the link: where C a (x a ) is the capacity, k a (x a ) is the jam density of link a in units of vehicles per mile, ρ is the calibration constant, and l H and l A are the maintained vehicle spacing for HDVs and AVs, respectively. x y a is the flow of class y and Γ y AV is a Boolean variable indicating if class y uses AVs entirely or not. Under the assumption 2l A > l H , the resulting volume-delay function satisfies the monotonicity condition with respect to each single class which is the necessary condition (but not sufficient) for the convexity of a multiclass problem. e Frank Wolfe algorithm was used to solve the equilibrium problem applied to Austin, Texas, and downtown network. Although the results of the numerical analysis showed the convergence of the algorithm, Levin and Boyles mentioned that the algorithm may not converge, and multiple equilibria are certainly possible.
A recent study of Mehr and Horowitz [27] on the equilibrium state of traffic network in a mixed AV and HDV environment showed that introducing capacity impacts of autonomous vehicles in the mixed flow makes the multiclass equilibrium problem more complex. ey found that the equilibrium problem does not have a unique solution in a general network with both AVs and HDVs. In this problem, neither class-specific link flows nor total link flows are unique. e results are based on BPR volume-delay function with the following form: where f a is the total link flow of link a, t 0 a is the free flow travel time, c a and β a are nonnegative parameters, and C a (α a ) is the capacity of link a which is a function of the proportion of AVs on link a. Despite other available capacity functions for estimating link capacity in a mixed condition, which will be discussed later in Section 2.3, capacity was defined as follows: where m a and M a are the capacity of link a when all vehicles on the link are HDVs and AVs, respectively (m a < M a ). μ a is defined as the ratio of m a to M a and is conceptually the degree of capacity asymmetry of link a.
Mehr and Horowitz proved that using the volume-delay function of form (5) and the capacity function of form (6), as well as assuming a homogeneous degree of capacity asymmetry (i.e., same μ a for all links) in a general network with a fixed AV market penetration, the uniqueness of total travel time is guaranteed for all Wardrop equilibrium flow vectors.
Wang et al. [29] developed a static multiclass traffic assignment model for a mixed AV and HDV environment. eir proposed model offers a different type of route choice behavior for each group of road users. AV users are assumed to follow the user equilibrium principle while HDV users are assumed to choose their route according to a crossnested logit model. ey considered the fact that AVs could have lower reaction time than HDVs which can lead to an increased link capacity. To estimate link travel time, they used the BPR function with the link capacity formulated as follows: where C a is the capacity of link a, C a,A and C a,H are the capacity of link a when all vehicles on the link are AVs and HDVs, respectively, and variables v a,A and v a,H are the AV and HDV flows on link a. ey have applied a new routeswapping-based solution algorithm to solve the problem. e suggested algorithm is path-based and can only find one local optimal depending on the initial path flow pattern.

Capacity Impacts of Autonomous Vehicles in the Heterogeneous Traffic Flow.
Several studies suggest that AVs can improve roadway capacity through the possibility of platooning and shorter safe-intervehicle distance compared to HDVs. Tientrakool et al. [30] compared the highway capacity of three vehicle technologies: manual vehicles, vehicles with sensors, and communicating vehicles. e results indicated that when all vehicles were equipped with sensors, the capacity was about 1.5 times of when all vehicles were manual. Furthermore, communicating vehicles could lead to a 237% increase in the highway capacity.
Chen et al. [31] proposed an analytical formulation for capacity of a single-lane and a two-lane highway with heterogeneous traffic of both AVs and HDVs. ey extend their formulation to a general case of a multilane highway. According to their study, the AV penetration rate, platoon size, and spacing characteristics were the parameters that affect capacity. ey also proposed capacity functions under different lane policies including exclusive AV or HDV lanes and mixed-used lanes.
ey found that policies that segregate AVs and HDVs (i.e., have exclusive HDV and AV lanes), led to a lower capacity than policies including mixedused lanes.
Ghiasi et al. [15] formulated a mixed traffic capacity as a function of the AV market penetration, platooning intensity, and mixed traffic headway settings based on a Markov chain model. ey showed that unlike the generally accepted assumption, higher market penetration and platooning intensity did not necessarily lead to higher capacity, since headway setting is dependent on the AV technology and plays an important role in determining the capacity.
Lazar et al. [19] proposed two capacity models under two different scenarios. e main assumption in the first scenario was that AVs are capable of maintaining reduced headways regardless of the preceding vehicle. e second scenario assumed that the headway depends on the preceding vehicle, and a short headway is possible when an AV follows another AV. In this case, the distribution of vehicles on the road is important. is study presented a capacity formulation assuming that vehicle types are determined as a result of the Bernoulli process. In both scenarios, capacity was the function of proportion of AVs on the road. e capacity calculated in the first scenario served as an upper bound for the capacity in the second scenario [32].

Problem Statement and Formulation
In this study, we consider a network with two types of lanes: regular lanes and AV lanes. Two modes use the transportation network: AVs and HDVs. AVs can use all lanes in the network and HDVs can only use regular lanes. e demand for each origin-destination (O-D) pair is assumed fixed for each travel mode. e route choice behavior of both modes is based on Wardrop's first principle.
We formulated AV lane location problem as an integer bilevel programming (IBP). e IBP consists of optimally selected lanes from the existing network to be dedicated to AVs. Below we summarize the notations used in the model. x w,m a : flow of travel mode m ∈ M on link a ∈ A between OD pair w ∈ W x a : aggregate flow on link a ∈ A y k : if one lane of the kth link pair is converted to an AV lane, y k � 1; otherwise y k � 0 C a : total capacity of link a ∈ A c a : per-lane capacity of link a ∈ A p a : proportion of AVs on link a ∈ A

Network Representation. Let G(N, A) represent a general transportation network with node set N and link set A. Define A as a set of candidate links to introduce an AV lane (A ⊂ A).
Each candidate AV link is paired with a regular link. Candidate links are links with more than two lanes because there should be at least one regular lane between each OD pair to maintain network connectivity. Figure 2 illustrates a sample network topology. In Figure 2, the solid line represents regular links, and the dotted line denotes the candidate AV links. e set of paired links is defined as K. Specifically, in Figure 2

Multiclass User Equilibrium
Problem. An equilibrium model for a general multiclass transportation network can be Mathematical Problems in Engineering written as a Variational Inequality (VI) problem of equation (8) [33,34]. e solution of the following VI formulation with the feasible region in Ω is x * . It can be shown that x * satisfies the multiclass user equilibrium condition as defined by Dafermos [25]: e feasible region Ω in a multiclass transportation network with exclusive AV lanes is mathematically described by equations (9)-(13). Constraint (9) simply ensures the conservation of flow in the network. Constraint (10) is defined on sets of regular links (i.e., A\A). Nonnegativity of flows is expressed by (10) and (12). By (10) and (11), we allow HDVs on regular links, and by (12) x w,1 a ≥ 0, ∀a ∈ A\A, ∀w ∈ W, x w,1 a � 0, ∀a ∈ A, ∀w ∈ W, x w,2 a ≥ 0, ∀a ∈ A, ∀w ∈ W, x a � m∈M w∈W e link performance function, which calculates travel time, is the BPR function in the form of equation (14), where t 0 a is the free flow travel time of link a and C a is the capacity of link a. In this study, we assume that the capacity of each link is the per-lane capacity times the number of lanes (equation (15). e per-lane capacity is considered similar to that of Mehr and Horowitz [27] (equation (16) with homogeneous degree of capacity asymmetry (i.e., for all links (h AV /h HDV ) is constant).
is function gives the lane capacity based on the share of AVs in the flow and also ensures that a unique total travel time is resulted despite the fact that multiple Wardrop equilibria might exist. According to Lazar et al. [19], the proposed capacity function gives an upper bound for capacity; thus, we consider ξ as a capacity reduction factor in mixed traffic. It is a realistic assumption since it is unlikely that shared use lanes reach their maximum capacity due to complex interactions between AVs and HDVs. ξ is a positive coefficient less than 1.0 for shared use lanes and is equal to 1.0 for exclusive AV lanes. Link autonomy level which is defined as the share of AVs on a link is calculated by equation (17): C a p a � c a p a L a ,

AV Lane Location
Problem. e problem of determining the optimal location of AV lanes in the network is a bilevel program in which the lower level is a multiclass user equilibrium expressed by equations (8) Objective (18) is to minimize the total travel time of the transportation network where x a is determined by the lower level user equilibrium model. A binary variable y k is introduced to indicate whether an AV candidate link is converted to AV link or not. d k a is defined as a parameter which can take three values as below:  e capacity of link a is calculated by equation (19). In brief, the problem is to find the optimal location of AV lanes in a network with mixed traffic flow of AVs and HDVs while minimizing the total travel time as the objective function and accounting for capacity variation. us, given the transportation network with a specified demand matrix, the output of the model is the optimal location of AV lanes. Below, we summarize basic assumptions of our model:

Solution Method
e transportation network design problem (TNDP) is NPhard even in the simplest form [35]. Several solution methods, such as exact methods, heuristics, and metaheuristics have been developed to solve the TNDP [22,36,37]. Hybrid heuristics and metaheuristics approaches are other methods that can be used to solve network design problems.
e exact methods such as the branch and bound algorithm proposed by Leblanc [38] will become computationally inefficient in case of large-size networks due to its high accuracy. On the contrary, heuristics and metaheuristics are more efficient but less accurate methods. Hybrid methods which combine the two methods are more efficient than one heuristic or metaheuristic algorithm, as shown by some studies [39][40][41] To solve the bilevel problem proposed in this study, we applied a heuristic method proposed by Bagloee et al. [42] which is a hybrid machine learning and optimization. e method employs a multivariate linear regression and an integer linear program to solve the bilevel problem iteratively. e algorithm starts with a feasible solution. en, the traffic assignment problem (i.e., the lower-level problem) is solved based on this solution, and the objective function of the upper-level problem is computed. e objective function is regressed against the decision variables, and the coefficients of decision variables are calibrated. e next step is to solve an integer linear program (ILP), which finds a new feasible decision variable for the next iteration. e lowerlevel problem is solved based on the new solution, and a new regression is applied to the updated data. e termination criterion is the number of iterations, which is a user-specified parameter. e algorithm is applicable to large scale networks and can quickly converge to a local optimum [42]. Step 0: specify i max ; set i � 0; set y k 0 � 0 as the initial feasible solution.
Step 1: solve the traffic assignment problem for the given y k i and compute the objective function Z i . Step 2: estimate new regression coefficients of equation (18) using y k i and the corresponding Z i : Step 3: update the ILP of equations (19)- (25) according to the current binary solution y k i and estimated coefficients: Step 4: if the termination criterion is met (i.e., i � i max ), stop and report the best objective function (Z i ) and corresponding (y k i ) as the final solution, else, solve the updated ILP in step 3, introduce y k i+1 as a new feasible solution, set i � i + 1, and go to Step 1.

Basic Settings.
e performance of the proposed model on the Sioux Falls network is presented in this section. e network has 76 links and 24 nodes and its data is available at a repository maintained by Bar-Gera [43] including node coordinates, start node, and end node of each link, capacity, length, and free flow travel time of each link. e given OD flows are 0.1 of the original daily flows in LeBlanc et al. [43].
ere is no explicit data on the number of lanes of each link. erefore, we have converted the given capacity to the number of lanes assuming that capacity is defined for the same period as OD flows and per-lane capacity of each link is 2400 vehicles per hour. Table 1 shows the results of the calculation. Note that the values calculated for number of lanes are rounded up to the nearest integer. Since the only constraint for a link to be a candidate of introducing an AV lane is to have more than one lane, there are 16 two way links that make our alternative set. Dotted links in Figure 3 are the candidates to get an AV lane. In Figure 3, both directions of a link are considered as one alternative (e.g., (1,2) and (2, 1) is considered as one alternative to get an exclusive AV lane or not). Note that each candidate for an AV lane is paired with a regular link, as shown in Figure 2. e capacity reduction factor (ξ) in mixed traffic is set to 0.8 for a mixed-used lane and 1 for an exclusive AV lane. e headway of AVs and HDVs from the preceding vehicle is 0.85 and 1.5 seconds, respectively.   To solve the proposed bilevel problem, a code in python is developed which is synchronized to GAMS (the General Algebraic Modeling System), which is used to formulate the ILP problem as one part of the heuristic algorithm explained in Section 4. e rest part of the algorithm is coded in python.

Optimal Solution.
During the transition period, market penetration rate grows from very low values to relatively high values until it reaches market saturation. Within the first stages of AV implementation, not much benefit would be obtained from the introduction of AV lanes to the network; however, more benefits are expected at higher market penetration rates. To approve this hypothesis we have conducted a sensitivity analysis on market penetration rates, the results of which are shown in Table 2. e benefits began to be significant at the 40% penetration rate, while improvements at lower rates (i.e., 10%, 20%, and 30%) are so marginal that a network manager may prefer not to implement any AV lanes. According to the analysis, an improvement in the objective function value was maximized at the market penetration rate of 60%. e benefits of deploying AV lanes were reduced at both lower and higher rates. e optimal location plan at a market penetration rate of 50% is shown in Figure 4 as an example. is solution contains a selection of 12 alternatives out of 16 possible alternatives. e total travel cost at the optimal scenario was 89371 vehicle hours, which was 3.26 percent better than the base scenario (i.e., do nothing).
In addition to the optimal plan, a practical plan is also proposed. If we want to present a practical AV lane deployment plan to the network manager while accounting for growth in the AV market penetration, we should consider practical hassle of adding or removing an AV lane every year during the transition period. Investigating the results shows that there are common alternatives in the optimum plan of most market penetration rates. For instance, alternatives 1 and 10 are selected in all penetration rates above 20%. Using this observation and taking into account the fact that little benefit is gained when the market penetration rate was below 40%, it can be recommended to start the implementation of AV lanes when 40% of all vehicles on the network are autonomous. e optimal solution at 40% includes alternatives 1, 7, 8, 10, 11, 12, and 14. We add alternatives 3,9,14, and 16 at 50%, and for 60%, 70%, and 80% alternative 2 is added to the previous set of alternatives. When the AV market size grows to 90%, fewer lanes need to be converted to exclusive AV lanes; therefore, alternatives 3, 7, and 8 were no longer selected to have an AV lane. Some alternatives that were present in the optimal plan are not considered in the practical plan to decrease the practical hassle of adding and removing AV lanes every year during the transition period. Running the traffic assignment and calculating the objective function for this deployment plan (see Table 3) indicates that there is a little difference between this plan and the optimal plan of Table 2; therefore, the network benefits from implementation of AV lanes are maximized with minimum disturbances to the network structure.
In order to show the computational efficacy of the applied algorithm, the number of iterations and the CPU time to reach the optimal solution is presented in Table 4.   We have used a laptop computer with Intel Corei7 of 2.70 GHz and 12 GB of RAM. It is worth noting that the optimal solution in all cases is the global optimal solution which is already identified by the exhaustive enumeration. e results show that the algorithm converges to the optimal solution in early iterations with a reasonable CPU time.

e Variable Capacity Case Versus the Fixed Capacity
Case. Table 5 shows the optimal plan for different market penetration rates assuming a fixed per-lane capacity. In this case, the capacity of each regular lane is set to 2400 vehicles per hour and the capacity of each AV lane is set to 4200 vehicles per hour, which equals to 1.75 ((h AV /h HDV ) � (1.5/0.85) � 1.75) times the per-lane capacity of a regular lane. e comparison of the fixed capacity with variable capacity cases indicates that not only the optimal plan varies at each market penetration rate but also the trend of estimated benefits obtained by the AV lane deployment is significantly overestimated. As shown in Table 5, the reduction in the total travel time caused by an AV lane deployment increased up to market penetration rate of 60% and remained constant with the increase of AV market size. Compared to its corresponding result in the variable capacity case, the benefits are overestimated in the fixed capacity case. e reason is that, in the variable capacity case, the base scenario total travel time (TTT bs ) decreases as the AV market size grows. is decrease is due to the dependence of capacity to market penetration rate, while the fixed capacity case ignores the change of capacity and thus, takes the base scenario total travel time as constant. It can be concluded that the introduction of AVs to the network improved the TTT bs itself, as shown in Figure 5 and, therefore, the estimated benefit in the fixed capacity case is an overestimated value.
is point can have significant implications in network planning and should be considered in the process of decision making by the network manager to have a reasonable expectation from the network improvement obtained by an AV lane deployment.

Conclusion
is paper proposes a mathematical model that determines the best scheme of road space reallocation between AVs and HDVs by locating AV lanes in a transportation network. e bilevel optimization problem is formulated as a discrete network design problem which is an NP-hard problem and a hybrid machine-learning and optimization method is employed to solve it. e location of AV lanes is specified at the upper-level problem. e multiclass assignment problem at the lower level is formulated as a variational inequality since it is not a convex problem. e main contribution of the model is investigating the capacity impacts of autonomous vehicles at the network level. To achieve this goal, we used the BPR volume-delay function as the travel time function in the route choice procedure. e link capacity is one of the parameters of BPR function which is usually assumed as a constant value due to the link specifications. In this study the link capacity is considered as a function of AV proportion in a link, and, therefore, it can take different values according to the class-specific link flows. e model was applied on the Sioux Falls network and the numerical results were presented at various market penetration rates. It was found that the benefits (i.e., reduction in total travel time) gained from AV lane deployment is highly affected by the market penetration rate. To derive benefits from the implementation of AV lanes requires relatively wide adaptation of AVs (i.e., more than 30%). e results indicate that the benefit is maximum at medium ranges (i.e., 60%) of the market penetration rate. At low market penetration rates (i.e., 10%, 20%, and 30%), AV lane deployment does not lead to significant network improvements. e results also point out that ignoring the capacity variation would lead to different optimal plans as well as overestimation in benefits. A practical deployment plan during transition period, when market penetration rate grows from low levels to the saturation rate, is presented in this study. e practical plan considers the disturbances that adding or removing an AV lane may cause in the network. is study can be followed up in different directions. In the proposed model, the total travel time serves as the objective function, while it is possible to add other social costs and account for different values of time of each user class. It is also interesting to examine the results using different capacity functions developed in various studies. e model can be extended to a multistage problem in which time is one of the dimensions. Since the introduction of AVs can affect travel demand and urban development, another field of research is to integrate the optimal AV lane location problem with AV travel impacts such as spatial distribution [44,45], ridesharing [46], parking demand [47], and mode choice [48], which are examined separately in some previous studies.

Data Availability
e data used to support the findings of this study are available from Leblanc et al. [43].

Conflicts of Interest
e authors declare that they have no conflicts of interest.