Solving a Multi-Class Trafﬁc Assignment Model with Mixed Modes

: In comparison to conventional human-driven vehicles (HVs), connected and automated vehicles (CAVs) provide beneﬁts (e


Introduction
In the transportation system of the future, connected and automated vehicles (CAVs) might replace human-driven vehicles (HVs) [1], and drivers' travel behaviors will change dramatically in the future when CAVs become more common [2]. As CAVs will reduce human driving errors, travel safety could be significantly enhanced [3] and travel time will be improved due to an increased capacity by reducing headway. Tientrakool et al. [4] suggested that highway capacity might be doubled if 60% of traffic is CAVs and improved four to five times if all vehicles were changed to CAVs. Tran and Bae [5] demonstrated that CAVs can increase the average speed by 1.27 times. Besides these two primary benefits, there are countless more, such as improved trip quality and reduced fuel use. However, before the period of fully CAVs appears, it will encounter a situation in which both HVs and CAVs are present, and the traffic flow pattern of the mixed modes will be different when compared to that of only a single class (e.g., HV or CAV). In the literature, in terms of mixed flows, Gkartzonikas and Gkritza [6] focused on the value of time when people drive using CAVs, whereas Levin and Boyles [7] and Xie et al. [8] investigated the optimal control of the mixed traffic flows based on the road level. Chen et al. [9,10] studied the network design problem for allocating dedicated lanes and zones for CAVs. Recently, Li et al. [11] and Zhang and Nie [12] used mixed behavior equilibrium theories to simulate CAV flows. Levin and Boyles [7] demonstrated that even if a small part of a HV is replaced by a CAV, the network flows may significantly lower travel time.
To summarize, estimating network flows is critical for transportation planning during the period when both CAVs and HVs are present. As CAVs can follow the car in front of • User equilibrium model [13][14][15][16]; • Nonadditive traffic equilibrium problem [17]; • Stochastic user equilibrium problem [18,19]; • Elastic demand traffic assignment problem [20].
However, in real networks, the use of GP to solve the multi-class TAP with asymmetric vehicle interactions is quite limited. In this paper, the GP technique was extended to solve a multi-class TAP that explicitly considered mixed modes (e.g., HV and CAV). The remainder of the paper is organized as follows. After the introduction, the next section provides the modeling issues, which were vehicle interactions between HVs and CAVs and the multi-class TAP. Section 3 describes a GP method intended specifically to solve the proposed multi-class TAP. Two numerical experiments were adopted to demonstrate the proposed solution algorithm and model formulation in Section 4. Section 5 presents concluding remarks.

Multi-Class Traffic Assignment
This section discusses the travel cost function for HVs and CAVs and the model formulation for the multiclass TAP with asymmetric vehicle interactions. A list of variables is provided first for convenience. Vehicle interactions occur when multiple vehicle classes interact in the same roadway space. These interactions typically occur between vehicles with varied maneuvering and operational characteristics. Due to the lower reaction times described in the introduction, CAVs have a smaller headway than HVs. With this assumption, headway in mixed flows can be described as a function of the proportion of HVs and CAVs [7]. The function of headway based on the proportion of mixed flows is shown in Equation (1).

Travel Cost Function
The capacity of a road refers to its maximum traffic density, which can be expressed as follows in terms of headway.
If only a single class (e.g., HV or CAV) is present, we can obtain Equation (3), and Equation (1) can be written as Equation (4) using Equation (3).
With the derived headway in Equation (4), the capacity for mixed flows can be expressed as Figure 1 with Equation (2).
Vehicle interactions occur when multiple vehicle classes interact in the same roadw space. These interactions typically occur between vehicles with varied maneuvering a operational characteristics. Due to the lower reaction times described in the introducti CAVs have a smaller headway than HVs. With this assumption, headway in mixed flo can be described as a function of the proportion of HVs and CAVs [7]. The function headway based on the proportion of mixed flows is shown in Equation (1).
The capacity of a road refers to its maximum traffic density, which can be express as follows in terms of headway. If only a single class (e.g., HV or CAV) is present, we can obtain Equation (3), a Equation (1) can be written as Equation (4) using Equation (3).
With the derived headway in Equation (4), the capacity for mixed flows can expressed as Figure 1 with Equation (2).     Figure 2a illustrates a typical link travel time function (i.e., a Bureau of Public Roads (BPR) function) for the TAP and the used capacity for HVs and CAVs. Using the input data in Figure 2a, Figure 2b illustrates travel time change and capacity change as a function of the proportion of CAVs. As CAVs have a smaller headway than HVs, the total capacity for mixed flows increased as the proportion of CAVs increased. However, the increasing rate was not a linear pattern. With increased capacity, travel time was intuitively decreased, and this was also not a linear trend. Hence, the travel time function was asymmetric with regard to HVs and CAVs because the travel time of each mode was affected by the other modes. The following link travel time function in Equation (5) was used to model the asymmetric interactions between various modes in a general situation involving a variety of vehicle types.    Figure 2a, Figure 2b illustrates travel time change and capacity change as a function of the proportion of CAVs. As CAVs have a smaller headway than HVs, the total capacity for mixed flows increased as the proportion of CAVs increased. However, the increasing rate was not a linear pattern. With increased capacity, travel time was intuitively decreased, and this was also not a linear trend. Hence, the travel time function was asymmetric with regard to HVs and CAVs because the travel time of each mode was affected by the other modes. The following link travel time function in Equation (5) was used to model the asymmetric interactions between various modes in a general situation involving a variety of vehicle types.

Model Formulation
Let c and f be the vectors of route costs . . . , c rs km , . . . T and route flows . . . , f rs km , . . . T .
Then, the multi-class TAP with asymmetric vehicle interactions can be formulated. A variational inequality (VI) problem is shown as follows: Find a route flow solution f * ∈ Ω, such that where Ω represents the feasible set defined by Equations (7)-(9): ∑ k∈K m rs f rs km = q rs m , ∀rs ∈ RS, m ∈ M Appl. Sci. 2022, 12, 3678 Equation (7) is the travel demand conservation constraint; Equation (8) is a definitional constraint that sums up all route flows that pass through a given link a; and Equation (9) is a non-negativity constraint on the route flows. The VI formulation meets Wardrop's equilibrium conditions [21] based on the Karush-Kuhn-Tucker (KKT) conditions of the VI formulation [22].
where o rs m is the minimum O-D cost of class m between O-D pair rs.

Solution Algorithms
The single-class traffic assignment using the GP algorithm was initially applied by Jayakrishnan et al. [14]. The GP algorithm is based on Bertsekas's [23] GoldsteinLevitin-Polyak (GLP) projection approach, which was designed to address the nonlinear multicommodity problem in communication networks. In each iteration, the algorithm solves the shortest route problem to find a descending search direction, and then scales the result using the second derivative information. The following are the updated equations.
Bertsekas et al. [20] recommended utilizing α(n) = 1 as the step size and using the second derivative information to avoid determining an appropriate step size. This approach significantly reduced the computational time required to determine the step size and has been demonstrated to be quite robust in previous experiments using the "one-at-a-time" flow updating strategy on networks of various sizes and topologies [13,14]. Furthermore, demand conservation constraints were incorporated. As could be seen in the formulation, the TAP has only two constraints, which are the travel demand conservation constraint in Equation (7) and a non-negativity constraint for the decision variable in Equation (9), while the route travel time consisted of the sum of the link travel time. From this special structure, it was possible that demand conservation constraints in Equation (7) were embedded into the objective function, as in Equation (12). As a result, the simple projection just needed to verify that the route flows were not negative, which was readily accomplished by making the route flow zero if it results were negative. Figure 3 illustrates the solution procedure strategy for the multi-class TAP with asymmetric vehicle interactions.
12, x FOR PEER REVIEW 6 of 12

Numerical Experiments
This section conducts numerical experiments on a small network and a real network to demonstrate the capabilities of the multi-class TAP in mixed modes and the efficiency of the route-based GP algorithm.

Small Network
The network's structure and characteristics are presented in Figure 4. It consisted of seven links, five nodes, and two O-D pairs (O-D pairings (1)(2)(3)(4) and (1)(2)(3)(4)(5)). To compute the travel time, the BPR link performance function ( a α = 0.15 and a β = 4.00) was used, and the travel time weight parameter ( m p ) was considered to be 1.0.
• Update link costs: • Find the shortest path based on : • Update path set: • Determine an appropriate step size α • Update non-shortest path flows: • Update shortest path flows • Set n=0, rs=0, m=1 • Initialize path flows and link flows: • Initialize path set:    Figure 3. Solution procedure strategy for the multi-class TAP with asymmetric vehicle interactions.

Numerical Experiments
This section conducts numerical experiments on a small network and a real network to demonstrate the capabilities of the multi-class TAP in mixed modes and the efficiency of the route-based GP algorithm.

Small Network
The network's structure and characteristics are presented in   To begin, we assessed the network's overall travel time as a result of multi-class TAP in Figure 5. As expected, when the percentage of CAVs increased, overall travel time decreased due to increased capacity. However, overall travel time was reduced nonlinearly, although the percentage of CAVs was raised linearly. Increases in the proportion of CAVs at an early stage (0 percent to 10% of CAVs) were more effective than increases at a later stage (90 percent to 100 percent of CAVs). Additionally, when the percentage of CAVs was 100%, capacity improved by 50%, as shown in Figure 2b, while overall travel time improved by 60%. As a result of the findings, we could infer that the network was highly congested when only HV demand was used.    To begin, we assessed the network's overall travel time as a result of multi-class TAP in Figure 5. As expected, when the percentage of CAVs increased, overall travel time decreased due to increased capacity. However, overall travel time was reduced nonlinearly, although the percentage of CAVs was raised linearly. Increases in the proportion of CAVs at an early stage (0 percent to 10% of CAVs) were more effective than increases at a later stage (90 percent to 100 percent of CAVs). Additionally, when the percentage of CAVs was 100%, capacity improved by 50%, as shown in Figure 2b, while overall travel time improved by 60%. As a result of the findings, we could infer that the network was highly congested when only HV demand was used.   To begin, we assessed the network's overall travel time as a result of multi-class TAP in Figure 5. As expected, when the percentage of CAVs increased, overall travel time decreased due to increased capacity. However, overall travel time was reduced nonlinearly, although the percentage of CAVs was raised linearly. Increases in the proportion of CAVs at an early stage (0 percent to 10% of CAVs) were more effective than increases at a later stage (90 percent to 100 percent of CAVs). Additionally, when the percentage of CAVs was 100%, capacity improved by 50%, as shown in Figure 2b, while overall travel time improved by 60%. As a result of the findings, we could infer that the network was highly congested when only HV demand was used.    Second, by increasing the percentage of CAVs, the trajectories of the link capacity were explored and presented in Figure 6. When the percentage of CAVs was increased from 0% to 60%, the capacities of links 2, 3, 6, and 7 increased, while the capacities of the other links (i.e., links 1, 4, and 5) increased above 60% of CAVs. This result indicated that more CAVs were allocated to routes 3 and 6 until 60% of CAVs was reached, while CAV flows were increased in routes 1 and 4 after 60% of CAVs was reached. This was because the links in routes 3 and 6 were approaching their maximum capacity, and the excess CAV flows were diverted to other underutilized links to increase capacity and decrease travel time accordingly. As could be seen, capacity expansion also did not follow a linear trend. Second, by increasing the percentage of CAVs, the trajectories of the link capacity were explored and presented in Figure 6. When the percentage of CAVs was increased from 0% to 60%, the capacities of links 2, 3, 6, and 7 increased, while the capacities of the other links (i.e., links 1, 4, and 5) increased above 60% of CAVs. This result indicated that more CAVs were allocated to routes 3 and 6 until 60% of CAVs was reached, while CAV flows were increased in routes 1 and 4 after 60% of CAVs was reached. This was because the links in routes 3 and 6 were approaching their maximum capacity, and the excess CAV flows were diverted to other underutilized links to increase capacity and decrease travel time accordingly. As could be seen, capacity expansion also did not follow a linear trend.

Korea Network
This section conducts numerical experiments on the South Korea network. These experiments were performed to investigate (a) the convergence characteristics of the route-based GP algorithm for solving the multi-class TAP and (b) the effect of flow allocation. Figure 7 depicts the Korea network, which consists of 247 zones, 45,434 nodes, 96,621 links, and 51,541 O-D pairs. The Korea Transportation Institution provided the network structure, O-D trip tables, and link travel time function. Four scenarios were created to investigate the effect of the multi-class TAP with mixed modes on different roadway classes. Scenario 1 assumed that CAVs were not permitted to drive on any roads, which meant they drove manually on all roadways like HVs. On the other hand, Scenario 2 assumed that CAV travel was permitted exclusively on freeways. It was assumed that CAVs were permitted on freeways and arterials in scenario 3. Finally, scenario 4 assumed that CAVs were capable of traveling on all available roadways. The aggregated demands were assumed to be disaggregated by 50% and 50% to obtain the O-D trips for HVs and CAVs, respectively. Maximum capacity was expected to be 1.5 times more when driving only CAVs than when driving only HVs. The GP algorithm was coded in Intel Visual FORTRAN XE 2021 and ran on a 3.60 GHz processor with 16.00 GB of RAM.

Korea Network
This section conducts numerical experiments on the South Korea network. These experiments were performed to investigate (a) the convergence characteristics of the routebased GP algorithm for solving the multi-class TAP and (b) the effect of flow allocation. Figure 7 depicts the Korea network, which consists of 247 zones, 45,434 nodes, 96,621 links, and 51,541 O-D pairs. The Korea Transportation Institution provided the network structure, O-D trip tables, and link travel time function. Four scenarios were created to investigate the effect of the multi-class TAP with mixed modes on different roadway classes. Scenario 1 assumed that CAVs were not permitted to drive on any roads, which meant they drove manually on all roadways like HVs. On the other hand, Scenario 2 assumed that CAV travel was permitted exclusively on freeways. It was assumed that CAVs were permitted on freeways and arterials in scenario 3. Finally, scenario 4 assumed that CAVs were capable of traveling on all available roadways. The aggregated demands were assumed to be disaggregated by 50% and 50% to obtain the O-D trips for HVs and CAVs, respectively. Maximum capacity was expected to be 1.5 times more when driving only CAVs than when driving only HVs. The GP algorithm was coded in Intel Visual FORTRAN XE 2021 and ran on a 3.60 GHz processor with 16.00 GB of RAM.   Figure 8 illustrates the convergence characteristics of the route-based GP method used to solve the multi-class TAP with mixed modes. As can be seen, the method was capable of achieving convergence with a relative gap (RG) of 1 × 10 −6 . The computational effort required to achieve a relative gap of 1 × 10 −6 was around 3600 s. Based on the results, computational CPU times were reduced as the available roads for CAVs increased. This was expected since adding available roads for CAVs to the network resulted in a reduced congestion level with increased capacity and the algorithm required less computational efforts to equilibrate the path flows of all vehicle classes.  Figure 8 illustrates the convergence characteristics of the route-based GP method used to solve the multi-class TAP with mixed modes. As can be seen, the method was capable of achieving convergence with a relative gap (RG) of 1 × 10 −6 . The computational effort required to achieve a relative gap of 1 × 10 −6 was around 3600 s. Based on the results, computational CPU times were reduced as the available roads for CAVs increased. This was expected since adding available roads for CAVs to the network resulted in a reduced congestion level with increased capacity and the algorithm required less computational efforts to equilibrate the path flows of all vehicle classes.   Figure 8 illustrates the convergence characteristics of the route-based GP method used to solve the multi-class TAP with mixed modes. As can be seen, the method was capable of achieving convergence with a relative gap (RG) of 1 × 10 −6 . The computational effort required to achieve a relative gap of 1 × 10 −6 was around 3600 s. Based on the results, computational CPU times were reduced as the available roads for CAVs increased. This was expected since adding available roads for CAVs to the network resulted in a reduced congestion level with increased capacity and the algorithm required less computational efforts to equilibrate the path flows of all vehicle classes.

Assigned Flow Comparison
To illustrate the assigned flow difference and effect of travel time alleviation when CAVs drove on different roadway classes, Figure 9 illustrates the link flow differences between scenario 1 and other scenarios on a Geographic Information System (GIS) map. It should be noted that the scenarios were based on the specifications presented in Figure 7c. Links were categorized by thickness and color to emphasize the degree of the difference. The green color showed that the assigned flow in scenario 1 was higher than the assigned flows in other scenarios, and the red color indicated the reverse. The results in Figure 9 clearly show that the multi-class TAP considering mixed modes had a significant impact on flow allocations. The majority of the differences occurred on the freeway, where the freeway generally had a larger capacity; hence capacity was increased more when CAVs traveled more on the freeway compared to other roadways. Using vehicle kilometers traveled (VKT) and vehicle hours traveled (VHT) to evaluate the effect of allowing different rode classes for CAVs, VKT was not significantly affected by CAV driving (the values changed by 0.4%, 0.2%, and −0.5% in scenarios 2, 3, and 4, respectively, compared with scenario 1, whereas VHT was significantly decreased when CAVs drove on more roadways (the values changed by −2.3%, −4.1%, and −8.3% in scenarios 2, 3, and 4, respectively, compared with scenario 1). The reasons for the results were that more CAVs driving on the roadways could promote an increase in the overall network capacity, and this increased capacity influenced travel time reduction. Specially, when CAVs traveled more on the freeway, capacity on the freeway was increased compared to other roadways and travelers changed their route passing through the freeway to minimize their travel time, despite the travel distance increase. As a result, VKT was slightly increased in the comparison between scenario 1 and scenario 2 and between scenario 1 and scenario 3.

Assigned Flow Comparison
To illustrate the assigned flow difference and effect of travel time alleviation when CAVs drove on different roadway classes, Figure 9 illustrates the link flow differences between scenario 1 and other scenarios on a Geographic Information System (GIS) map. It should be noted that the scenarios were based on the specifications presented in Figure  7c. Links were categorized by thickness and color to emphasize the degree of the difference. The green color showed that the assigned flow in scenario 1 was higher than the assigned flows in other scenarios, and the red color indicated the reverse. The results in Figure 9 clearly show that the multi-class TAP considering mixed modes had a significant impact on flow allocations. The majority of the differences occurred on the freeway, where the freeway generally had a larger capacity; hence capacity was increased more when CAVs traveled more on the freeway compared to other roadways. Using vehicle kilometers traveled (VKT) and vehicle hours traveled (VHT) to evaluate the effect of allowing different rode classes for CAVs, VKT was not significantly affected by CAV driving (the values changed by 0.4%, 0.2%, and −0.5% in scenarios 2, 3, and 4, respectively, compared with scenario 1, whereas VHT was significantly decreased when CAVs drove on more roadways (the values changed by −2.3%, −4.1%, and −8.3% in scenarios 2, 3, and 4, respectively, compared with scenario 1). The reasons for the results were that more CAVs driving on the roadways could promote an increase in the overall network capacity, and this increased capacity influenced travel time reduction. Specially, when CAVs traveled more on the freeway, capacity on the freeway was increased compared to other roadways and travelers changed their route passing through the freeway to minimize their travel time, despite the travel distance increase. As a result, VKT was slightly increased in the comparison between scenario 1 and scenario 2 and between scenario 1 and scenario 3

Conclusions
The GP algorithm for solving the multi-class TAP with mixed modes was developed in this paper. It proposed the use of asymmetric travel time functions to represent the effect of vehicle interaction between HVs and CAVs. Two numerical experiments were performed to investigate the proposed model's features. The numerical results indicated the following: (1) when the proportion of CAVs was increased, travel time decreased significantly; (2) travel time was reduced nonlinearly although the percentage of CAVs was raised linearly; (3) there was a considerable change in the assigned traffic pattern with and without CAV demand; (4) the flow pattern on freeways was more affected by CAV driving; and (5) computational time was more efficient when CAV driving was allowed on more roads. These findings have the potential to have a considerable impact on various fields of research, including road construction, safety, and environmental impact studies. In the future, the issue of CAV performance in terms of travel time function and route choice behavior should be investigated. Another possibility is to consider other types of CAV, such as trucks (medium and heavy trucks) to improve the realism as each type of CAV has unique characteristics. Additionally, computational experiments should be conducted on a variety of network topologies, congestion levels, and zonal densities to determine the route-based GP algorithm's effectiveness and efficiency.

Conclusions
The GP algorithm for solving the multi-class TAP with mixed modes was developed in this paper. It proposed the use of asymmetric travel time functions to represent the effect of vehicle interaction between HVs and CAVs. Two numerical experiments were performed to investigate the proposed model's features. The numerical results indicated the following: (1) when the proportion of CAVs was increased, travel time decreased significantly; (2) travel time was reduced nonlinearly although the percentage of CAVs was raised linearly; (3) there was a considerable change in the assigned traffic pattern with and without CAV demand; (4) the flow pattern on freeways was more affected by CAV driving; and (5) computational time was more efficient when CAV driving was allowed on more roads. These findings have the potential to have a considerable impact on various fields of research, including road construction, safety, and environmental impact studies. In the future, the issue of CAV performance in terms of travel time function and route choice behavior should be investigated. Another possibility is to consider other types of CAV, such as trucks (medium and heavy trucks) to improve the realism as each type of CAV has unique characteristics. Additionally, computational experiments should be conducted on a variety of network topologies, congestion levels, and zonal densities to determine the route-based GP algorithm's effectiveness and efficiency.

Data Availability Statement:
The data presented in this study are available, upon request, from the corresponding author.