Location Design of Electrification Road in Transportation Networks for On-Way Charging

Electric vehicles tend to be a great mobility option for the potential benefits in energy consumption and emission reduction. Onway charging (OWC) has been recognized to be a promising solution to extend driving range for electric vehicles. Location of the electrification road (ER) is a critical issue for future urban traffic management to accommodate the new mobility option. -is paper proposes a mathematical program with equilibrium constraint (MPEC) approach to solve this problem, which minimizes the total travel time with a limited construction budget. To describe the drivers’ routing choice, a path-constrained network equilibrium model is proposed to minimize their travel time and prevent running out of charge. We develop a modified active set algorithm to solve the MPEC model. Numerical experiments are presented to demonstrate the performance of the model and the solution algorithm and analyze the impact of charging efficiency, battery size, and comfortable range.


Introduction
Electrification systems based on renewable energy power sources are introduced in the urban transportation system for positive environmental effect and carbon reduction [1,2]. Recently, with the increasing concern about sustainable transportation [3], electric vehicles (EVs) are widely adopted in urban travel. Electrification road (ER) is considered as a promising domain for sustainable electrical energy harvesting to support on-way charging (OWC) for EVs, thanks to the advancement of application of photovoltaic, piezoelectric, and pyroelectric materials converting solar energy, physical pressure, and thermal energy into electrical energy, e.g., [4][5][6][7][8]. OWC provides a new mode of charging for EV drivers to extend driving range, who are suffering from range anxiety of running out of power on the way [9,10] and long charging time ranging from 0.5 to 2 hours for a full charge [11]. OWC has raised the interests of organizations and is being tested all over the world. ER deployed on a network has been tested in South Korea [12]. UK has conducted a series of OWC recharging tested on highways [13]. In Sweden, the Volvo Group and the Swedish Transport Association planned to build 500-meter ER for wireless power transfer [14]. Such tests and experiments show that OWC would come true in future. Owing to the high investment in the ER, which reaches $4 million per lane mile [15], it is imperative to find the optimal location.
Most previous studies related to the location design of charging facilities are based on user equilibrium (UE) problems with EVs which investigate how the limited driving range of EVs and location of charging facilities affect the routing choice and subsequently the flow equilibrium distribution on road networks. Among these studies, Jiang et al. [16] first proposed a path-constrained assignment model where lengths of selected routes are within the driving ranges of EVs. e model is further extended to consider mixed gasoline and electric vehicle flows and their combined routing choices under different problem settings [17,18]. He et al. investigated the optimal prices of electricity and the integrated prices of electricity and roads based on the multiclass spatial distribution of electric vehicles across the transportation network [19,20]. All these studies assumed the energy consumption of EVs is flow-independent. e aforementioned studies do not consider the location of charging facilities in networks. Problems related to the optimal location of charging stations or battery swapping stations have received much attention. He et al. [21] first proposed a UE-based framework for locating charging stations to maximize social welfare, considering route choices of electric vehicles, and further developed three models of different flow dependencies and energy consumption [22,23]. Xu et al. [24] defined a network with battery swapping stations where swapping time is independent with charging capacity, considering the impact of queuing to find the best charging strategy. ese existing relevant studies are intended for charging stations or battery swapping stations, where EVs reach the station, stop, and get replenished. Such models cannot address issues for OWC directly.
OWC, unlike charging stations or battery swapping stations, can charge EVs during the ride. OWC technology has been implemented in a network in South Korea [12]. A limited number of research studies have been done on the OWC problem to select the optimal ER locations so far. Riemann et al. [25] proposed a mixed-integer nonlinear program to describe the flow-capturing location model considering driver's routing behavior. Chen et al. [26] investigated the optimal deployment of charging lanes based on the speed control of EV drivers on the charging lanes. Chen et al. [27] proposed the integrated deployment of charging lanes and charging stations along a long traffic corridor, which cannot be applied in the transportation network directly. Liu and Song [28] adopted a robust optimization methodology to provide robust location optimization of wireless charging facilities for the electric bus system to reduce total cost of batteries and charging facilities with different uncertainty levels. Liu et al. [29] developed a mixed-integer linear program for OWC location and battery size optimization and battery size for an electric bus system with overlapping bus lines. Ahmed et al. [30] introduced a method to find the best combination of battery capacity and wireless charger characteristics to solve the tradeoff between maximizing charge sustaining, minimum battery capacity, and minimizing the initial investment. Bi et al. [31] adopted a genetic algorithm to optimize the rollout of OWC infrastructure both spatially and temporally in order to minimize life cycle costs and energy burdens. Zhao et al. [32] proposed a biobjective optimization problem for integrated EV location and on-board battery size design of electric bus systems to minimize deployment cost and reduce energy consumption of electrified systems. e above research studies can solve the ER location optimization to some extent, but they tend to assume that charging amount is only related to the charging time. Under such assumption, drivers tend to deliberately slow down to extend the charging time, and even EV may revisit the same charging lanes. All the vehicles behind the EV that is moving slowly to get charged are delayed. A road contains more than two lanes, one of which is the ER. Vehicles without charging intention would perform lane-changing, which reduces the road capacity. Relationship between the cost of building the ER and width is seldom considered.
In this paper, charging amount is assumed to be related with the driving distance, while charging and the candidates of the ER have to satisfy the minimum energy supply. e width of the ER is considered to affect ER building cost. e link performance function is set as the BPR function in numerical examples. erefore, the route choice of drivers will be complicated, and the UE flows of EVs in the transportation network will be different.
is paper investigates the location design problem of the ER based on the analysis of the routing behavior of EV drivers, which is different from the previous studies with the charging station or swapping station. Our model inherits from the paper written by He et al. [21], but charging behaviors and routing behaviors are different tracking back to the specific situation of OWC in contrast to charging at station, which leads to a totally different location design. Assuming that the energy consumption rate and recharging rate are flow-independent, this paper proposes a path-constrained network equilibrium model (PCNE) considering the routing and recharging behavior of EV drivers as well as the constraint of driving range. en, we develop an iterative solution framework to solve the PCNE problem. Lastly, we investigate a mathematical program with equilibrium constraint (MPEC) model to optimize the ER location to minimize the system total travel time under a given construction budget limit. With the proposed model, discussions about the effect of different recharging rates, comfortable ranges, and battery sizes on the location design and system performance are conducted. e remainder of the paper is organized as follows. e notation, assumption, formulation of the PCNE problem, and a solution framework are elaborated in Section 2. Section 3 formulates the location design problem as a MPEC program. Section 3.1 proposes a modified active set method to solve the problem. Section 4 presents numerical cases based on the Sioux Falls network, and Section 5 offers the conclusions of this study and implications for future research on the transportation network of the ER. Assumption 3. Given an ER location design, a physical path from the origin to the destination may not be feasible for EV drivers due to the limit for the EV battery.
Assumption 4. All vehicles in networks are EVs. is assumption can be relaxed as the model can be extended to accommodate both EVs and regular vehicles.

Assumption 5.
Recharging rate is not affected by speed of the EV.
When EV drivers travel from their origins to destinations, they are assumed to select routes to minimize their travelling cost. Travelling cost includes electricity cost and travel time cost. Electricity cost is much smaller than travel time cost [23]. erefore, we simply assumed that EV drivers select routes of the least travelling time. While running on the ER, EV drivers can decide whether to charge and charging amount. Charging amount is not greater than the maximum charging amount that equals charging rate multiplied by the distance of the ER.
Consider a transportation network containing ER running EVs, denoted by G � (N, A), where N and A are sets of nodes and link, respectively. e link (i, j) ∈ A is directional and emanates from the starting node i ∈ N to the ending node j ∈ N with the length D ij and the road grade G ij . Travel demands are between a set of O-D pairs W. q od represents the travel demand between the O-D pairs, od ∈ W, and R od represents the set of paths between the O-D pairs.

Formulation.
Following definitions are introduced here to describe the PCNE problem.

Definition 1.
A path is usable if an EV is able to complete the path with the remaining charge state never falling below the comfortable range.
EV drivers are assumed to select routes of the least travel time among all the usable paths. At least one usable path exists for each O-D pair here for model generality. As the energy consumption rate and recharging rate are assumed to be flow-independent, the usability of a path is independent of the drivers' routing choice, which can be predetermined before flow distribution given the location design of the ER. e following path-constrained network equilibrium proposed by He et al. [23] is adopted here.

Definition 2.
At the path-constrained network equilibrium, all utilized paths are usable, and travel time of all utilized paths of a given O-D pair is the same, which is not more than that of any unutilized usable paths of the same O-D pair. e formulation of the path-constrained network equilibrium is as follows: PCNE: where R od represents the set of all usable paths between a given O-D pair, od ∈ W. Energy consumption rate and charging rate depend on the travel time, which are flowindependent. erefore, the usability of any path is flowindependent, which can be calculated. f � (. . . , f od r , . . .) is a vector of path flow. f od r represents the traffic flow on path r ∈ R od of O-D pair, od ∈ W. Constraint (2) demonstrates the flow balance for each O-D pair. Constraint (3) ensures nonnegativity of flow for each usable path.
Here, χ od ij,r is a binary variable to represent whether path r traverses link (i, j) ∈ A, which equals 1 if path r traverses link (i, j) ∈ A, and 0 otherwise. v ij represents the traffic flow Here, assume that link travel time takes the following form of Bureau of Public Roads (BPR): where t 0 ij represents the free-flow travel time of link (i, j) and c ij represents the capacity of link (i, j).

Solution Procedure.
If we can enumerate all usable paths beforehand, the proposed PCNE problem is a regular nonlinear optimization which can be solved easily by commercial nonlinear solvers such as CONOPT. However, enumeration of usable paths is a time-consuming work. Here, we adopt the solution procedure proposed by He et al. [22] to obtain the solution.
e procedure starts with a subset of R od , od ∈ W, and solves a restricted version of PCNE defined upon the subset. en, the feasibility of the solution to the restricted PCNE is tested. If not, a new usable path will be generated and added to the subset. Iteration proceeds until termination. Some new variables are introduced in the subproblem. Here, l max and l 0 are the battery size and the initial state of charge. For an EV travelling between the O-D pair, od ∈ W, the state of charge at node i is l od i . δ ij is a binary parameter representing whether link (i, j) ∈ A is equipped with the ER to support OWC. e variable equals 1 if the street (i, j) is equipped with the ER, and 0 otherwise. A p ⊆A denotes the set of all links equipped with the ER. Let ω ij denote the energy consumption rate on link (i, j) ∈ A and ϖ denote the recharging rate of the ER. erefore, link energy consumption is denoted as c ij � (ω ij − δ ij ϖ)D ij . m od is the minimum charge state within the comfortable range for drivers of the O-D pair, od ∈ W. N t (n) is the set of tails of those links heading to node n, and N h (n) is the set of heads of those links emanating from node n. K and M are sufficient large constants. x od ij is a binary variable indicating utilization incidence, which equals 1 if link (i, j) is utilized for travel demands between the O-D pair, od ∈ W, and 0 otherwise. Accordingly, ρ od ij is a variable that equals 0 if link (i, j) is utilized for travel demands between the O-D pair, od ∈ W , and is unrestricted otherwise. ξ n d is a binary variable, which equals 1 when n � d, and 0 otherwise; ξ on is a binary variable, which equals 1 when n � o, and 0 otherwise; and e od ij represents the recharging amount at link (i, j) ∈ A.
For an O-D pair, od ∈ W, the link flow obtained (. . . , v [k] ij , . . .) at the kth iteration of PCNE is adopted at the shortest usable path solution framework (denoted as SUPFk). e formulation of SUPF-k is shown as follows: SUPF-k: where Ω consists of the following conditions: In the above, the objective function is to minimize the travel time. Constraint (6) guarantees the balance of traffic flow; constraint (7) specifies the relation between the states of the charge of EV batteries travelling from node i to node j, link (i, j) ∈ A. Constraint (8) ensures the recharging amount of electricity on the ER at link (i, j) ∈ A P does not exceed the maximum charging quantity that equals charging rate multiplied by the distance of the link. Constraint (9) sets the upper and lower bounds of the states of the charge of the EV. Constraint (10) suggests that the EV driver can only recharge when travelling on the ER. Constraint (11) sets the lower bounds of the states of charge of the EV during the trip to satisfy comfortable range for drivers. Constraint (12) specifies the initial state of charge. Constraint (13) is a given electricity consumption function with respect to the road grade. Constraint (14) requires x od ij to be binary. e formulation assumes that EV drivers can control the recharging process on the ER. ey can determine the quantity of recharging on the ER. SUPF-k is a mixed-integer linear program, which can be easily solved by commercial solvers like CPLEX 12.8 for small-or medium-sized problems. e solution to SUPF-k for each O-D pair, denoted as ij , . . .), is adopted in the set construction of the shortest usable path. e iterative steps of PCNE is shown as follows: Step 0: set the iteration variable k � 0 and traffic flow Step 1: solve the restricted NE upon R od . Denote ij , . . .) and (. . . , μ od , . . .) as the optimal solutions and multipliers associated with constraint (2).
Step 2: for each O-D pair, od ∈ W, solve SP-(k + 1) and obtain the optimal solution, (. . . , , and go to Step 1, k � k + 1. ij , . . .) is the equilibrium link flow distribution. e above procedure terminates in a finite number of steps.

Location Design of the ER
is section investigates the location design under a given budget limit to maximize social welfare (i.e., minimize the total travel time of all EVs in our paper). e following assumptions are introduced here.
Assumption 5. Construction cost of the ER is related to length, width of the street. Assumption 6. When a street is designed to be equipped with the ER, vehicles driving on any lane of the street can get charged. at is, streets in the network are not mixed with charging lanes and regular lanes.
When the ER location is determined, the equilibrium traffic flow can be derived, as well as the recharging and routing behavior. e problem of the location design for electrification road (LDER) is formulated as follows: LDER: where Ψ(δ) � f | r∈R od f od r � q od , ∀od ∈ W; f od r ≥ 0, ∀od ∈ W, r ∈ R od (δ)}; W ij is the width of link (i, j) ∈ A; B is the construction cost of the unit-area ER, which is the function of the recharging rate ϖ; Θ is the budget limit; and G is a sufficiently large constant.
In the above, the objective function represents the total driving time. Constraint (16) presents the shortest length of the street chosen as the ER, guaranteeing that Assumption 1 always holds. Constraint (17) specifies the budget limit for the ER deployment. Constraint (18) requires δ ij to be binary. Constraint (19) states the traffic flow following the result predicted by the PCNE model.

Active Set Method.
LDER is a typical bilevel discrete network design problem with NE at lower level. Many existing solutions from the literature can be explored, for example, branch and bound [33], simulated annealing [34], and SO relaxation [35] with consensus of reformulating an equivalent single-level model with equivalent optimality conditions of the lower-level problem like KKT conditions to find the solution of the bilevel program. Problems with the nonconvex feasible region generally contain many local optimal solutions and typically require a time-consuming branch-and-bound scheme to search for a globally optimal solution. Such algorithms are not readily applicable considering the nonlinear and complementary constraints, adding enormous amount of calculation.
To solve the LDEP, we utilize the active set method [36], which terminates after a finite number of iterations ending with a strongly stationary solution. e active set method has been proved to have the potential for solving larger network design problems. e active set method assigns nonnegative variables in each pair as zero to construct the location design problem as a regular nonlinear program. Adopting multipliers associated with the constraints forcing nonnegative variables to be zero, the proposed algorithm formulates the binary knapsack problem to solve the zero-value assignments to decrease the total travel time.
e active set algorithm starts by setting an initial feasible solution U � δ 0 ij |(i, j) ∈ A and considering two sets: where Ω 0 ∪Ω 1 � A and Ω 0 ∩Ω 1 � ϕ. Under the location design for the ER corresponding to the initial feasible solution U 0 , the NE problem can be solved easily by the pathbased iteration procedure proposed in Section 2. Let V * � v * ij , (i, j) ∈ A be the solution set of NE, which is unique because of the monotonically increasing BPR function. e basic principle of the active set algorithm is to exchange elements between Ω 0 and Ω 1 in order to reduce the overall travel time in each iteration until the system optimal solution is found. In each iteration, the active sets are adjusted to meet the following theorems proved by Zhang et al. [36]. Theorem 1. Given that Ω 0 and Ω 1 are feasible solution sets, if λ ij < 0 for some (i, j) ∈ Ω 0 , switching (i, j) from Ω 0 to Ω 1 yields less system travel time. If μ ij > 0 for some (i, j) ∈ Ω 1 , switching (i, j) from Ω 1 to Ω 0 yields less system travel time.

Theorem 2.
e active set method converges after a finite number of iterations. λ ij and μ ij are the Lagrangian multipliers associated with equations (20) and (21), respectively.
To ensure that budget limit and shortest-length limit are satisfied in each iteration according to constraints (16) and (17), an embedded program is proposed as follows: minimize where g ij � 1 records a shift in (i, j) ∈ Ω 0 to Ω 1 and h ij � 1 records a shift in (i, j) ∈ Ω 1 to Ω 0 . e objective of programming (22) is to minimize the estimated decrease, the negative value of which implies a potential reduction of the objective value. e process terminates when equation (22) is zero. Constraint (23) is consistent with constraint (16). Constraint (24) is consistent with constraint (17). Parameter θ is set to guarantee a decrease in the objective function. In the first iteration, θ � − ∞. In the next iteration, θ can be calculated by the following equation: where ε is a sufficiently small constant. Equation (27) prevents the solution from degenerating back to the previous iteration by increasing θ by ε.

Journal of Advanced Transportation
Steps of the algorithm are shown in the following. Here, SO denotes the objective value of programming (15), and TD denotes the objective value of programming (22). (Algorithm 1).

Basic Settings. In this section, numerical examples on the Sioux Falls network in South
Dakota are presented to demonstrate the performance of the proposed model as shown in Figure 1, which consists of 24 nodes, 76 links, and 14 OD pairs. All links are available to be equipped with the ER. Table 1 reports the link characteristics, including freeflow travel time (FFTT), width, length, and capacity. e OD demand is given in Table 2.
We assume that the battery size l max is 25 kWh. e initial state of the charge l 0 is set as 25% l max . For simplicity, it is assumed that energy consumption rate ω at all links in the network is set as 0.3(kWh/mile), and the recharging rate ϖ is set as 2.5(kWh/mile). e comfortable range for all drivers m od is set as zero. Assume that the cost of building unit-area ER is (($2 × 10 5 )/(mile × meter)) at the recharging rate ϖ � 2.5 (kWh/mile). e budget limit Θ is set as $200, 000, 000. All the above values are chosen for illustrative purpose. We adopt CPLEX12.8 to solve the PCNE, SP, and LDER.

Base Case Result.
e LDER problem is solved under the aforementioned basic setting, as shown in Figure 2(a). Figure 2(b) depicts the convergence performance of the SP-k algorithm, indicating the potential of the SP-k algorithm in solving the PCNE problem. e equilibrium link flows are reported in Table 3, with the total travel time of 9.20 × 10 4 h · veh and the average travel time 0.575 h per vehicle.

Strategy Comparison under Different Budget Limits.
In this section, we consider different location designs of the proposed model under different budget limits to demonstrate the system performance, as shown in Figure 3(a). Figure 3 compares the total travel time, cost of the ER, and the number of ERs under different budget limits. As expected, total travel time of the network decreases with an increase in the budget. More ERs are built with an increasing budget when the budget ranges from $1.3 × 10 8 to $2 × 10 8 , generating more usable route choice for the EV to complete the tour, which promotes average distribution of traffic flow leading to reduction in the system travel time. Nevertheless, more ERs are built under the budget of $1.2 × 10 8 , comparing with the budget of $1.3 × 10 8 and $1.4 × 10 8 . It can be observed from Figure 3(a) that when the budget exceeds $1.4 × 10 8 , link (20,18) is always chosen to be equipped with ER. e usability of link (20,18) has a significant impact on EV drivers' routing behavior due to its special location. However, ER on link (20,18) of the 15 m width is an expensive project, which has to be weighed according to the budget limit. Figure 3 If not, derive (g ij , h ij )) and the objective value TD η and go to step 4. Obtain a new solution to (Ω η+1

Strategy Comparison with Different Recharging Rates.
In order to investigate the impact of the charging rate on location design and equilibrium flow (see Figure 4), we solve the DLEP and NE with charging rate ϖ set as 2kWh/mile (LEVEL1), 2.5kWh/mile (LEVEL2), and 3.5kWh/mile (LEVEL3). Cost of the ER is set as $1 × 10 5 , $2 × 10 5 , and $3 × 10 5 , respectively. Budget limit is set as $1.9 × 10 8 . Interestingly, the minimum total travel time can be obtained under this budget in LEVEL 1 and LEVEL 2 cases. at is, even when the budget limit increases, the total travel time will not be reduced. To explore the minimum travel time in LEVEL 3, we repeat to solve the DLEP and NE in the LEVEL 3 case with the budget limit increasing gradually and find the minimum travel time of 9.13 × 10 4 h · veh at the budget of $2.73 × 10 8 , which is less than the travel time of 9.20 × 10 4 h · veh in the LEVEL 2 case. Figure 4(b) compares the total   1  20  16  20  1  16  1  13  8  13  1  8  1  7  12  7  1  12  2  13  10  13  2  10  2  20  8  20  2  8  7  13  14  13  7  14  7 20 12 20 7 12 Journal of Advanced Transportation 7 travel time and cost in these four cases. It can be observed that, under the budget limit of $1.9 × 10 8 , the ER of LEVEL 2 gets the best performance, which is much better than the other two cases. ree explanations are offered here. First, due to constraint (16), shortest length of the ER is different as recharging rate changes. In the case of LEVEL 1, only 26       the ER location design. In this section, the impact of comfortable range on location design and system performance is tested with DLEP of different comfortable ranges m od � 1.5 kWh and m od � 2kWh under the budget limit of $2.1 × 10 8 and $2.7 × 10 8 solved, respectively, comparing with m od � 0. Location designs of the ER are shown in Figure 5(a). Figure 5(b) compares the total travel time and cost in these six cases. Under the budget limit of $2.1 × 10 8 , change of m od from 0 to 2kWh increases the total travel time by 78%. To further explore the impacts of comfortable range, we raise the budget limit to $2.7 × 10 8 and find that the variations of the total travel time with the change in comfortable range are not obvious. at is to say, under a given budget of $2.1 × 10 8 , increased comfortable range leads to the increment in the total travel time. When the budget limit is set as $2.7 × 10 8 , the location design of m od � 0 is the same as that under the budget limit of $2.1 × 10 8 , with the construction cost of $1.865 × 10 8 . e location designs of m od � 1.5 kWh and m od � 2kWh change to yield the minimum travel time, which is close to the number when m od � 0, with the construction cost increasing to $2.63 × 10 8 and $2.5 × 10 8 , respectively.

Impacts of the Battery Size on Location Design and System
Performance. Different battery sizes affect the routing choice and location design in two ways. First, the shortest length of the ER is different according to constraint (16). Second, the state of charge and charging behavior are different according to constraints (7) and (9). In this section, budget limit is set as $4 × 10 8 to explore the location design and system performance under different battery sizes l max � 20, 25, and 30 kWh. Initial state of charge is set as l 0 � 25% l max . As revealed in Figure 6, location design with the battery size of l max � 25 kWh has the least total travel time and construction cost of the ER. On the one hand, as commonly considered, a larger battery size allows more charging amount on the ER, which may generate more usable path and reduce the total travel time. On the other hand, constraint (16)  minimum length of the ER which is related to the battery size. When l max � 3 kWh, only 26 links whose length is more than 12 miles can be the candidate of the ER, which affects the set of usable paths and the total travel time.
From the above discussion, we can find that (1) System total travel time of the network decreases with an increase in the budget. More ERs are built with an increasing budget, generating more usable route choice for the EV to complete the tour, which promotes average distribution of traffic flow leading to reduction in the system travel time. (2) Recharging rate sets the bottom boundary for the length of the candidate of the ER, which affects the EV drivers' routing choice. Besides, recharging rate is related to the construction cost, which should be balanced with the budget limit. (3) Increase in comfortable range changes the sets of usable path. erefore, system travel time increases. (4) Battery sizes affect the shortest length of the ER as well as the initial state of charge and charging behavior, which should be balanced under different budgets.

Conclusion
To enhance the transportation network performance as ER is equipped, this paper proposes a LDER model to optimize the location of the ER, considering the routing and charging behavior of EV drivers by a PCNE model. EV drivers are assumed to decide their routing and recharging plan to minimize their travel time and prevent running out of charge. We develop a modified active set algorithm to obtain the optimal location of the ER. Numerical examples are conducted to demonstrate the performance of the model in comparison with the location designs in different budget limits, recharging rates, comfortable range, and battery sizes. e energy consumption rate and the recharging rate of the proposed path-constrained network equilibrium are assumed to be flow-independent. Our future study will extend our model to investigate the effect of the traffic congestion on the energy consumption rate and the recharging rate and try to find the solution algorithm. Besides, the transportation network can be extended to accommodate both EVs and regular vehicles, with impact of the mixed traffic flow investigated.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request. Disclosure e authors are responsible for all views and opinions expressed in this paper.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.