Electric Vehicle Charging Facility Planning Based on Flow Demand—A Case Study

: The location of electric vehicle charging facilities is of great signiﬁcance in promoting the use of electric vehicles. Most existing electric vehicle location models, including the ﬂow refueling location model (FRLM) and its ﬂexible reformulation (FRFRLM), are based on ﬂow demand. At present, these models cannot effectively deal with large-scale trafﬁc networks within a limited time, and there has been little comparison of their relative beneﬁts and limitations. Additionally, there have been few evaluations of the actual construction and location of charging facilities in cities. This paper describes an algorithm that can solve the large-scale transportation network problem within a reasonable time. Using this algorithm, the FRLM and FRFRLM models are compared in a case study focused on Jiading District, Shanghai, China, which provides some direction for the future development of ﬂow demand models. Finally, to evaluate the actual construction of urban charging facilities, this paper presents an algorithm that can map the actual charging facilities to the transportation network, and compares the actual construction situation with the model output. This enables a comprehensive evaluation of the actual construction of charging facilities and provides guidance for future construction.


Introduction
With continued economic development, increasing attention is turning to the problem of environmental pollution. Motor vehicle exhaust emissions are the main source of urban air pollution [1]. Related reports indicate that the petroleum resources of the earth can only be exploited for another 200 years, while other fossil fuel resources will also be fully exploited over a similar timescale [2]. With increasingly severe environmental problems and the looming energy crisis, vigorous developments in new energy sources are of great significance [3]. The development of electric vehicles is an important component of China's new energy development strategy [4]. Electric vehicles have the advantages of low vibration, low noise, and less pollution than petrol-or diesel-fueled vehicles [5]. In most countries, the cost of electricity consumption of electric vehicles is lower than that of ordinary vehicles. However, a number of factors are restricting the development of electric vehicles [6]. First, the driving range of electric vehicles is much shorter than that of traditional cars. Relevant studies have shown that owners of electric vehicles generally suffer from "range anxiety" [7]. Secondly, electric vehicles take a long time to charge. DC fast charging requires around half an hour to fill up to 80% of the battery capacity, whereas AC slow charging may take 6-8 h to fully recharge the battery [8]. Finally, electric

Literature Review
The location model based on flow demand was first proposed by Hodgson, who developed a flow-capture location model (FCLM) based on the maximum coverage [17]. In this model, the flow can be captured as long as there is a facility on the path. The model locates p facilities and maximizes the flow that passes through them. Hodgson used 25 transportation network nodes designed by Berman and Simchi-Levi to conduct a case study. On this basis, Kuby considered the driving range of the vehicle and proposed the FRLM. Different from FCLM, FRLM requires a series of facility combinations to satisfy the need for the vehicle to go back and forth between an origin-destination (OD) pair. FCLM assumes that only one facility is needed on the path for the flow to be captured. FRLM is divided into two stages. The first stage identifies feasible combinations of candidate sites so that the vehicle can successfully complete a round-trip between the OD pair. The second stage uses these combinations as the model input to determine the site locations that maximize the passing flow. Kuby used the 25 transportation network nodes designed by Berman and Simchi-Levi to conduct a case study, while Kuby and Lim further expanded the location of candidate sites, using algorithms such as mid-path segment, minimax, and maximin to disperse candidate sites on arcs [18].
On the basis of FRLM, Upchurch and Kuby restricted the service capacity of alternative facility points, and proposed a capacitated flow refueling location model (CFRLM) that considers capacity constraints [19]. In CFRLM, the service volume of each alternative facility point is different and the locations cannot exceed their service capacity. The objective is still to maximize the flow. On the basis of FRLM, Kim and Kuby considered the deviation of the vehicle path to find facility service points, and proposed the deviation-flow refueling location model (DFRLM) [20]. DFRLM reflects the driver's willingness to deviate from a pre-planned path to find service stations within a specified distance. However, in finding the shortest path to the chosen service station, DFRLM may cause the vehicle to pass too many nodes, so that the calculation time becomes too long. For this reason, Kim and Kuby developed a heuristic algorithm to solve this problem, which effectively reduces the runtime [21]. Based on the first stage of FRLM to find a feasible combination of candidate sites that cannot effectively deal with large-scale network problems in a short time, MirHassani and Ebrazi proposed their flexible reformulation model (FRFRLM) based on flow demand. They also used the 25 transportation network nodes designed by Berman and Simchi-Levi to conduct a case study, and randomly generated transportation networks of different sizes for further experiments. Wang proposed a model that minimizes the number of charging piles to be built by considering the charging time, and found that fast charging can effectively reduce the number of charging piles [22]. On this basis, Wang proposed an optimization model for the charging piles in the case of switching power considering the capacity limitations of the service facilities [23]. The results show that increasing the capacity of service facilities can effectively reduce the number of charging stations. In these two models, Wang used a seven-node transportation network to conduct case studies. In addition, Wang proposed a flow-based set cover model [24] based on vehicle refueling logic that minimizes the cost of building all charging piles while ensuring that vehicles can successfully go back and forth along a single path. A 51-node transportation network was developed for a case study evaluating this model.

Flow Refueling Location Model
Kuby considered the driving range of the vehicle in developing the original FRLM. The model can be written as where: q = index of OD pairs (and, by implication, the shortest paths for each pair); Q = set of all OD pairs; f q = flow volume on the shortest path between OD pair q; y q = 1 if f q is captured, 0 otherwise; k = potential facility location; K = set of all potential facility locations; x k = 1 if a facility is located at k, 0 otherwise; p = number of facilities to be located; The objective function in Equation (1) is intended to capture as great a flow as possible. Constraint (2) means that when y q is 1, a set of facility locations that can capture the flow must be built, and this set of facilities must be able to allow the car to pass back and forth on the path. Constraint (3) means that if a facility combination is to be built, then all facilities in that combination must be built. Constraint (4) means that we must construct p facilities.

Flexible Reformulation of the Flow Refueling Location Model
MirHassani and Ebrazi considered the flaws of the FRLM. When the network scale becomes larger, it takes too long to determine the combination of facilities that would allow the car to go back and forth on a certain path successfully. Based on the new vehicle round-trip logic, they developed the FRFRLM. The model can be written as where: q = index of OD pairs (and, by implication, the shortest paths for each pair); f q = flow volume on the shortest path between OD pair q; y i = 1 if a station is located at node i, 0 otherwise; x q ij = flow on an arc (i, j) belonging to the expanded network arcs of path q (i.e.,Â q ); Q = set of all candidate paths to be covered (usually the shortest path); N = set of all path-related nodes (N = U q∈Q N q ); Q i = subset of Q that contains all the paths passing node i.
The objective function in Equation (5) is intended to capture as much flow as possible. Constraint (6) means that the inflow at a point is equal to its outflow. Constraint (7) means that if the path selected by the car passes through a point, then a charging station must be Sustainability 2021, 13, 4952 5 of 23 built at this point. Constraint (8) means that the inflow at a point must be greater than or equal to 0. Constraint (9) means that we must construct p facilities.

Algorithm for Solving Large-Scale Transportation Network Problems
We have developed an algorithm for solving the large-scale network problem based on the flow demand model. The algorithm is divided into two stages. In the first stage, the transportation network is reconstructed by clustering, and the site types are obtained by solving the model. The second stage restores the traffic network and deletes un-sited classes from the transportation network, but does not change the traffic network structure. Algorithm 1 describes the proposed approach.

Input
Flow between nodes of the transportation network: f low; distance between nodes of the transportation network: dist; number of clusters: n; number of charging facilities to be built in the transportation network: p Output Location of charging facilities in the transportation network Step 5 modelFunction ← addConstraints(locationSet i , 1) Step 6 location ← modelFunction(p, netWork2)

1.
Clustering. We cluster the traffic network nodes into n classes according to the latitude and longitude coordinates of the traffic network nodes, so that the network can be effectively processed in a short time; 2.
Transportation network reconstruction. After clustering, it is necessary to reconstruct the distance between traffic network classes and the OD flow. The distance between classes is the farthest distance between points in the two classes, and the OD flow between classes is taken from all points in one class to all points in another class. For example, in Figure 1, the distance from class 1 to class 2 is the maximum distance from points 1, 3, and 6 in class 1 to points 25, 12, and 63 in class 2. If this does not exist, the distance between class 1 and class 2 is infinite. The OD flow of class 1 to class 2 is the sum of the flows from points 1, 3, and 6 to points 25, 12, and 63. Of course, in this way, part of the flow from origin to destination will be lost in the network. For example, the flow from point 1 to point 3 is lost; 3.
We use the model to solve the problem. We thereby obtain the clustered charging facility location class and the classes where charging facilities will not be located; 4.
Restore the transportation network. On the premise of not damaging the transportation network, the locations that will not have charging facilities are deleted from the transportation network. For example, as shown in Figure 1, assume that charging facilities are built in classes 1, 2, and 4. Class 3 no longer has a charging service facility. For the classes that do not have charging facilities, the length of the road section is retained but the node data are removed, that is, points 15, 2, and 31 are removed from the network. If the path from point 6 to point 199 passes through point 15, however, the distance is not changed by removing point 15; 5.
Add a constraint. For the charging facility location class, add the constraint that a charging facility must be located in this class. LocationSet denotes the clustered charging facility location class ∑ k∈locationSet i x k = 1, ∀ locationSet i ∈ locationSet;

6.
By the above steps, we obtain a transportation network that gets rid of the classes where charging facilities will not be located. We use the distance data between points and the OD flow data to bring into the model to obtain the layout of the charging facilities; 5. Add a constraint. For the charging facility location class, add the constraint that a charging facility must be located in this class. LocationSet denotes the clustered charging facility location class ∑ ∈ 1, ∀ ∈ ; 6. By the above steps, we obtain a transportation network that gets rid of the classes where charging facilities will not be located. We use the distance data between points and the OD flow data to bring into the model to obtain the layout of the charging facilities; Figure 1. Traffic network.

Algorithm for Mapping Real Charging Facilities to the Transportation Network
Currently, the construction of electric vehicle charging facilities cannot be evaluated. Thus, we have developed an algorithm that maps real charging facilities to the transportation network. We can then compare the results with the model output to evaluate the construction of regional electric vehicle charging facilities. Algorithm 2 describes the proposed approach.

Input
Parking and charging latitude and longitude data: 1; transportation network node latitude and longitude data: 2; distance between charging points in the same charging station: 1; number of clusters: ; distance: 2 Output Construction of charging facilities in the transportation network Step 1 1 , 1 1, 1 Step 2 3 ← Step 3 2 , 3 2 2 1 Note 3 denotes the cluster center points after clustering The data required for this algorithm are as follows: (1) latitude and longitude coordinate data of vehicle parking and charging on a working day; (2) latitude and longitude coordinate data of the traffic network nodes. The following parameters need to be set: (i) the distance between charging points in the same charging station is denoted as 1; (ii) the number of clustering categories is denoted as ; (iii) the distance between a traffic network node and the cluster center point is denoted as 2.
1. The vehicle charging data are screened. One piece of data from each charging station is retained, as this is convenient for later clustering of the charging stations. One piece of data for which the distance between the vehicle charging latitude and longitude coordinates is less than 1 is also retained;

Algorithm for Mapping Real Charging Facilities to the Transportation Network
Currently, the construction of electric vehicle charging facilities cannot be evaluated. Thus, we have developed an algorithm that maps real charging facilities to the transportation network. We can then compare the results with the model output to evaluate the construction of regional electric vehicle charging facilities. Algorithm 2 describes the proposed approach. Algorithm 2. Algorithm for mapping real charging facilities to the transportation network.

Input
Parking and charging latitude and longitude data: L1; transportation network node latitude and longitude data: L2; distance between charging points in the same charging station: d1; number of clusters: n; distance: d2 Output Construction of charging facilities in the transportation network L3 denotes the n cluster center points after clustering The data required for this algorithm are as follows: (1) latitude and longitude coordinate data of vehicle parking and charging on a working day; (2) latitude and longitude coordinate data of the traffic network nodes. The following parameters need to be set: (i) the distance between charging points in the same charging station is denoted as d1; (ii) the number of clustering categories is denoted as n; (iii) the distance between a traffic network node and the cluster center point is denoted as d2.

1.
The vehicle charging data are screened. One piece of data from each charging station is retained, as this is convenient for later clustering of the charging stations. One piece of data for which the distance between the vehicle charging latitude and longitude coordinates is less than d1 is also retained; 2.
Use the K-means clustering method to cluster the filtered data. Setting the number of clusters to n, we obtain n cluster center points after clustering; 3.
Map the cluster center points to the traffic network nodes. According to the distance between the latitude and longitude coordinates of the transportation network node and the cluster center point, if there is a cluster center point around d2 from a transportation network node, then the transportation network node can be considered as a charging facility construction point.

Case Study
We conducted a case study using data from Jiading District in Shanghai, China. The transportation network in Jiading District is shown in Figure 2. The OD flow of electric  The OD  flow table is presented in Table 1, and the distances between nodes in the traffic network of Jiading District are given in Table 2.
a charging facility construction point.

Case Study
We conducted a case study using data from Jiading District in Shangha transportation network in Jiading District is shown in Figure 2. The OD fl vehicles and the distance between the OD points in Jiading District during a are provided. There are 271 charging points and 25,259 OD pairs in a work OD flow table is presented in Table 1, and the distances between nodes in t work of Jiading District are given in Table 2.

Transportation Network Node
Transportation Network Node Flo   175  218  1  181  143  3  153  130  2  35  10  66  129  223  19  25  10  15  143  185  65  203  157  2  92  89  2  13 106 28   175  218  1  181  143  3  153  130  2  35  10  66  129  223  19  25  10  15  143  185  65  203  157  2  92  89  2  13  106  28  62  110  11  129  117  42  39  14  4  244 160 2 Under the assumption that the driving range of an electric vehicle is 150 km, we used the proposed algorithm and FRLM to solve this large-scale transportation network problem with the number of clustering categories set to 120. The transportation network problem was solved with 10, 20, 30, 40, 50, and 60 given facilities, respectively; the results of the first stage of solving the flow capture are presented in Table 3, and the situation of the location classes is shown in Figure 3. We used a 2.3 GHz Intel®Core TM i7-10875H processor and 16 GB RAM to solve this problem. The results of the second stage are presented in Table 4, and the situation of the sites is shown in Figure 4.  The results of the second stage are presented in Table 4, and the situation of the sites is shown in Figure 4.

Using FRFRLM and the Proposed Algorithm to Solve the Problem
Under the same conditions and assumptions as described in Section 4.1, FRFRLM and the proposed algorithm were used to solve the large-scale transportation network problem. The results from the first stage of solving the flow capture problem are presented in Table 5, and the situation of the location classes is shown in Figure 5. The results of the second stage are presented in Table 6, and the situation of the facilities is shown in Figure 6.

Using FRFRLM and the Proposed Algorithm to Solve the Problem
Under the same conditions and assumptions as described in Section 4.1, FRFRLM and the proposed algorithm were used to solve the large-scale transportation network problem. The results from the first stage of solving the flow capture problem are presented in Table 5, and the situation of the location classes is shown in Figure 5.

Comparative Analysis of Model Solution Coverage
The FRLM and FRFRLM flow capture coverage in the first stage is compared in Figure 7. Figure 8 compares the FRLM and FRFRLM flow capture coverage in the second stage.
From Figure 7, it can be seen that, given the different numbers of location points, the solution coverage of FRLM and FRFRLM is exactly the same in the first stage. Figure 8 shows that, for different numbers of locations, FRLM and FRRSLM also have essentially the same solution coverage rate in the second stage. In summary, there is only a slight difference in the results when the two models are used to solve the flow capture problem.

Comparative Analysis of Coincidence Rate of Model Site Selection
The coincidence rate of FRLM and FRRSLM site selection in the first stage is shown in Figure 9. Figure 10 shows the coincidence rate of FRLM and FRFRLM site selection in the second stage. Figure 9 indicates that, when the number of given location points is small, the location points output by FRLM and FRRSLM in the first stage completely coincide. When the number of given location points is large, the location points output by FRLM and FRFRLM are more diverse, but the majority of sites are the same. The reason for this may be the difference in the shortest path selection between FRLM and FRRSLM. It can be seen from Figure 10 that the results in the second stage are essentially the same as those in the first stage. In summary, when the number of given location points is small, FRLM and FRFRLM output the same location points. When the number of given location points is large, the output location points of FRLM and FRFRLM are slightly different, but the majority of locations are the same.  The results of the second stage are presented in Table 6, and the situation of the facilities is shown in Figure 6.

Comparative Analysis of Program Runtimes
The relationship between the number of FRLM and FRFRLM site selection points and the program runtime in the first stage is shown in Figure 11. The relationship between the number of FRLM and FRFRLM site selection points and the program runtime in the second stage is shown in Figure 12.
In Figure 11, the total number of alternative point categories for FRLM charging facilities in the first stage are 120, when the given facility categories are 65, 60, 55, 50, 45, 40, 35, 30, 25, 20, 15, 10, and 5, respectively, we can find that the FRLM runtime is always less than that of FRFRLM.
In the first stage, the alternative point categories are 120. When we give 65 facility categories, according to the FRLM and the solution of the first stage of the algorithm, the "65 facility categories" consist of 166 charging facilities. Therefore, we need to choose 65 charging facility from the 166 charging facilities in the second stage. The alternative point categories are 120 in the first stage, when the given facility categories are 65, 60, 55, 50, 45, 40, 35, 30, 25, 20, 15, 10, and 5, respectively. According to the FRLM and the solution of the first stage of the algorithm, the total number of alternative points for FRLM charging facilities in the second stage, from largest to smallest, are 166, 152, 151, 143, 130, 120, 109 Figure 12 that when the total number of alternative points exceeds 130, the runtime of FRLM is always greater than that of FRFRLM. The reason is that when the total number of site selection points is large, FRLM takes a long time to filter out the combinations that allow vehicles to travel between OD pairs. FRFRLM does not perform this process, but there are too many variables in the establishment of the model, resulting in a long runtime. In sum, when solving small-scale problems, the solution efficiency of FRLM is better than that of FRFRLM. When solving large-scale problems, FRFRLM is more efficient than FRLM.

Comparative Analysis of Model Solution Coverage
The

Comparative Analysis of Model Solution Coverage
The FRLM and FRFRLM flow capture coverage in the first stage is compared in Figure 7. Figure 8 compares the FRLM and FRFRLM flow capture coverage in the second stage.   From Figure 7, it can be seen that, given the different numbers of location points, the solution coverage of FRLM and FRFRLM is exactly the same in the first stage. Figure 8 shows that, for different numbers of locations, FRLM and FRRSLM also have essentially the same solution coverage rate in the second stage. In summary, there is only a slight difference in the results when the two models are used to solve the flow capture problem.

Comparative Analysis of Coincidence Rate of Model Site Selection
The coincidence rate of FRLM and FRRSLM site selection in the first stage is shown in Figure 9. Figure 10 shows the coincidence rate of FRLM and FRFRLM site selection in the second stage.     From Figure 7, it can be seen that, given the different numbers of location points, the solution coverage of FRLM and FRFRLM is exactly the same in the first stage. Figure 8 shows that, for different numbers of locations, FRLM and FRRSLM also have essentially the same solution coverage rate in the second stage. In summary, there is only a slight difference in the results when the two models are used to solve the flow capture problem.

Comparative Analysis of Coincidence Rate of Model Site Selection
The coincidence rate of FRLM and FRRSLM site selection in the first stage is shown in Figure 9. Figure 10 shows the coincidence rate of FRLM and FRFRLM site selection in the second stage.    Figure 9 indicates that, when the number of given location points is small, the lo tion points output by FRLM and FRRSLM in the first stage completely coincide. When number of given location points is large, the location points output by FRLM and FRFRL are more diverse, but the majority of sites are the same. The reason for this may be difference in the shortest path selection between FRLM and FRRSLM. It can be seen fro Figure 10 that the results in the second stage are essentially the same as those in the fi stage. In summary, when the number of given location points is small, FRLM a FRFRLM output the same location points. When the number of given location points large, the output location points of FRLM and FRFRLM are slightly different, but the m jority of locations are the same.

Comparative Analysis of Program Runtimes
The relationship between the number of FRLM and FRFRLM site selection points a the program runtime in the first stage is shown in Figure 11. The relationship between number of FRLM and FRFRLM site selection points and the program runtime in the s ond stage is shown in Figure 12. Running time of the program(s) Figure 10. Coincidence rate of FRLM and FRFRLM location points in the second stage for various numbers of site selection points.
ncidence rate of FRLM and FRFRLM location points in the second stage for various numbers of site selection Figure 9 indicates that, when the number of given location points is small, the location points output by FRLM and FRRSLM in the first stage completely coincide. When the number of given location points is large, the location points output by FRLM and FRFRLM are more diverse, but the majority of sites are the same. The reason for this may be the difference in the shortest path selection between FRLM and FRRSLM. It can be seen from Figure 10 that the results in the second stage are essentially the same as those in the first stage. In summary, when the number of given location points is small, FRLM and FRFRLM output the same location points. When the number of given location points is large, the output location points of FRLM and FRFRLM are slightly different, but the majority of locations are the same.

Comparative Analysis of Program Runtimes
The relationship between the number of FRLM and FRFRLM site selection points and the program runtime in the first stage is shown in Figure 11. The relationship between the number of FRLM and FRFRLM site selection points and the program runtime in the second stage is shown in Figure 12.    In Figure 11, the total number of alternative point categories for FRLM charging cilities in the first stage are 120, when the given facility categories are 65, 60, 55, 50, 45, 35, 30, 25, 20, 15, 10, and 5, respectively, we can find that the FRLM runtime is always le than that of FRFRLM.
In the first stage, the alternative point categories are 120. When we give 65 facil categories, according to the FRLM and the solution of the first stage of the algorithm, t "65 facility categories" consist of 166 charging facilities. Therefore, we need to choose charging facility from the 166 charging facilities in the second stage. The alternative po categories are 120 in the first stage, when the given facility categories are 65, 60, 55, 50, 40, 35, 30, 25, 20, 15, 10, and 5, respectively Figure 12 that when the total number of alternative points exceeds 1 the runtime of FRLM is always greater than that of FRFRLM. The reason is that when t total number of site selection points is large, FRLM takes a long time to filter out the co binations that allow vehicles to travel between OD pairs. FRFRLM does not perform t process, but there are too many variables in the establishment of the model, resulting in long runtime. In sum, when solving small-scale problems, the solution efficiency of FRL is better than that of FRFRLM. When solving large-scale problems, FRFRLM is more e cient than FRLM.

Summary of the Model Comparison
According to the above analysis, the models are summarized in Table 7.

Summary of the Model Comparison
According to the above analysis, the models are summarized in Table 7. Table 7. Comparison and analysis between FRLM and FRFRLM.

Comparing the Program Runtime between FRLM and FRFRLM
Modeling is based on flow demand and considers the limitations of vehicle driving range.
When the number of given site selection points is the same, the flow captured by the two models is almost the same.
When the number of given site selection points is small, the site selection points of the two models are exactly the same; when the number of given site selection points is large, the site selection points of the two models are mostly the same.
When the total number of traffic network nodes is less than 120, the program runtime of FRFRLM is always greater than that of FRLM. When the total number of traffic network nodes is greater than 130, the program runtime of FRFRLM is always less than that of FRLM. As the scale of the transportation network increases, the program runtime of FRLM increases significantly compared with that of FRFRLM.

Evaluation of Charging Facility Construction
Sample electric vehicle charging data for a working day in Jiading District, Shanghai, are presented in Table 8; the full dataset contains 8636 items. Longitude and latitude coordinate data of transportation network in Jiading District, Shanghai, are presented in Table 9; the full dataset contains 271 items. Using the algorithm that maps the real charging facilities to the transportation network, with d1 = 50 m, n = 120, and d2 = 500 m, we obtained the construction points of charging facilities on the network, as shown in Figure 13.  In Figure 14, according to the first stage output from FRLM, there are currently 34 real construction points in Jiading District in the model output construction category. The final FRLM output indicates that there are currently 15 real construction points that overlap with the model output. The first stage of the FRFRLM output shows 36 real construction points that are currently in the model output construction category. The final output d1 means the maximum distance between charging points in the same charging station; we think that 50 m is appropriate in Jiading District, because of the shortage of land resources in Jiading District. n means the number of clusters; there is no special requirement for this, usually we make it twice as many charging facilities as given in FRLM and FRFRLM. d2 means the distance between a cluster center point and a transportation network node. We think that if the distance between a cluster center point and a transportation network node is less than 500 m, then the transportation network node should be considered as having a charging facility. Of course, for different cities, there can be different d1 and d2.
In Figure 14, according to the first stage output from FRLM, there are currently 34 real construction points in Jiading District in the model output construction category. The final FRLM output indicates that there are currently 15 real construction points that overlap with the model output. The first stage of the FRFRLM output shows 36 real construction points that are currently in the model output construction category. The final output from the FRFRLM model indicates that there are currently 15 real construction points that overlap with the model output.

Conclusions and Future Research
Existing flow demand models cannot solve large-scale network problems quickly and effectively. This paper describes an algorithm that solves large-scale network problems through the clustering reconstruction of the transportation network. We compare FRLM with FRFRLM in terms of the flow capture results, site selection, and program runtime using data from Jiading District, Shanghai. When the given number of charging facilities are the same, we find that the flow capture results and site selection are essentially the same, while FRFRLM has an obvious advantage in program runtime when we solve large-scale transportation network problems. Although a series of extended models, such as CFRLM and DFRLM, have been developed based on FRLM, there are few extensions of FRFRLM. In the future, we can expand FRFRLM considering the capacity of charging facilities and the deviation of vehicles from the original shortest path to find available charging sites. Through our analysis of the current construction of charging facilities in Jiading District, Shanghai, we can conclude that there is considerable room for the optimization of such facilities. Different cities can get the current status of the city's charging facilities construction by adjusting the parameters that suit them.
Author Contributions: Data acquisition, curation, and analysis, C.W. and Z.G.; conceptualization, Z.G. and C.W.; funding acquisition, P.Y. and Z.L.; methodology, Z.G. and P.Y.; resources: C.W. and From the above data, it can be concluded that the overlap rate between electric vehicle charging facilities in Jiading District and the optimal solution of the model is about 25%, whereas the overlap rate between the construction type output in the first stage of the model is about 58%. From the perspective of transportation network nodes, there is currently a lack of charging facilities in the lower-left area of Jiading District. Generally speaking, according to the output of the model, there is considerable room for optimization in the construction of charging facilities in Jiading District, Shanghai.

Conclusions and Future Research
Existing flow demand models cannot solve large-scale network problems quickly and effectively. This paper describes an algorithm that solves large-scale network problems through the clustering reconstruction of the transportation network. We compare FRLM with FRFRLM in terms of the flow capture results, site selection, and program runtime using data from Jiading District, Shanghai. When the given number of charging facilities are the same, we find that the flow capture results and site selection are essentially the same, while FRFRLM has an obvious advantage in program runtime when we solve large-scale transportation network problems. Although a series of extended models, such as CFRLM and DFRLM, have been developed based on FRLM, there are few extensions of FRFRLM. In the future, we can expand FRFRLM considering the capacity of charging facilities and the deviation of vehicles from the original shortest path to find available charging sites. Through our analysis of the current construction of charging facilities in Jiading District, Shanghai, we can conclude that there is considerable room for the optimization of such facilities. Different cities can get the current status of the city's charging facilities construction by adjusting the parameters that suit them.