Optimal Planning of Electric Vehicle Charging Stations Considering Trafﬁc Load for Smart Cities

: The massive introduction of electric vehicles as a mobility alternative requires deploying an infrastructure of charging stations for electric cars (ICSEC). This new mobility concept will mitigate the environmental harm caused by the emission of CO2 generated by conventional internal combustion mobility methods. The sustainability of the ICSEC depends not only on the capacity to meet the demand for charging batteries for electric vehicles (EV) but also on an adequate number of public/private charging stations (CS) distributed in a geolocalized area. It is noted that the distribution of CS must respond to a set of real mobility constraints, such as vehicular ﬂow capacity, road capacity, and trajectories. The planning, intelligent location of public charging stations (PCS), and contingency analysis will enable us to study the increase in demand for electrical distribution substations (EDS). Therefore, the present model explains the need to plan the massive introduction of EVs by observing the user’s conditions at the trajectory level through ﬁnite resource allocation processes, segmentation, and minimum spanning trees, by observing heterogeneous vehicular ﬂow criteria through microscopic analysis, to understand the space–time relationship of vehicular ﬂow in georeferenced scenarios. Consequently, the computational complexity of the model is of the combinatorial type, and it is deﬁned as NP-Hard given the multiple variables and constraints within the transportation problem.


Introduction
The energy supply for charging EVs plays an essential role in the sustainability of the mass introduction of EVs.This new mobility concept contributes directly to reducing greenhouse gases and decreasing noise and visual pollution compared to conventional mobility.In Europe, it is estimated that approximately 80% of noise pollution caused on the streets is generated by vehicles that use fossil fuels as an energy source.Furthermore, it is kept in mind that human mobility has a fundamental role in the development and economic growth of towns; however, with conventional mobility, the environmental cost becomes critical and necessary for future generations due to the current radical climate change.
The most challenging problem is managing the supply chain for EV recharging; therefore, it is essential to plan the deployment of CS to EVs [1].Consequently, two optimizing variables are identified, (i) operating costs and (ii) heterogeneous charge balance [2].Therefore, studying and simulating the EV flow that should serve each CS in different zones within a geolocalized area is essential to plan, predict demand, and guarantee electric service continuity.Forecasting the peak demand in the CS is necessary because it would help us decide whether an electricity grid can take on new charges.If not, creating new supply circuits for each electric mobility CS is necessary.It would significantly contribute to maintaining the voltage stability in EDS and guaranteeing service continuity.Therefore, the planning can be divided into three stages: (i) planning to satisfy demand, (ii) allocation and location of CS, and (iii) EV routing to an energy source to recharge the vehicle's battery.
The article considers features in the model that bring the transport problem closer to different real-world scenarios.Microscopic examination in traffic models considers fundamental aspects such as speed and spacing between vehicles; it lets us know vehicle flow or density in geolocated scenarios.In addition, consider microscopic analysis to understand the space-time relationship in vehicular networks [3].
Therefore, a multi-objective approach to address the problems of routing, siting, and sizing of CS that are required for electric mobility charging infrastructures (CI) is noticed in this article.Capacities are allocated to each deployed CS to satisfy the demand, and the sum of the needed parts must equal the resource allocated in the planning problem.Candidate sites for siting CS could be parks, shopping malls, green spaces, and car parks.
Finally, the best topology for EV routing decisions is achieved by looking at road capacity, vehicle flow, and spacing variables.Hence, the aim is to solve multiple resource allocation problems in a georeferenced area of NP-hard computational complexity, which is supported by optimization tools, graph theory, and routing algorithms; it has been possible to reduce the computational complexity to a problem of NP-type [4].
Figure 1 presents the architecture proposed in this article.It shows a spatial approach on two levels: (I) route-based demands and (II) equipment for charging infrastructure.These demands must be covered by spreading a set of CS in permitted areas.Moreover, the equipment for a charging infrastructure should consider the forecast of vehicle demand as a function of the traffic density generated in georeferenced scenarios.Consequently, the heuristic proposed in this paper will solve the resource allocation by considering the required charging infrastructure equipment for EVs using real transport constraints at the level of trajectories under a spatial approach.Level I analyzes the flow generated on each road section as a function of vehicular density.Vehicular flow helps to understand traffic characteristics and behavior, which are essential requirements for planning, road and street design, and operation within the land transport system.There are microscopic and macroscopic models that relate different variables such as volume, speed, density, interval, and spacing [5] for understanding traffic density in geolocated scenarios.
Capacity and levels of service concepts applied to different road elements can be understood using the above models.However, it is not a trivial problem because these transport evaluation strategies are a great challenge to solve as they involve geographical elements (roads, orography, flyovers), vehicle arrangements are unknown, and due to their high traffic variability, it is not possible to predict with certainty trajectory patterns [6].Then, to develop optimal routes at the level of trajectories, it is essential to pay special attention to the aspects related to the variables of vehicular flow, the probabilistic description of vehicular traffic, the distribution of vehicles in viability, and the statistical distribution of vehicular density.
Consequently, level I studies and assigns a weight to each road section by observing actual mobility criteria for vehicle flow analysis with microscopic analysis.Level II examines the different topologies generated to achieve connectivity or supply of the multiple resources developed in the other CS to satisfy the EV demand in vehicle flow management processes [7][8][9][10].Vertices are defined as each track section's CS, EV, and intersections.Edges constitute the geolocated topology [11,12] of connection and the connectivity relationship between the multiple vertices that are part of a directed graph.
Henceforth, the present article is organized as follows: Section 2 briefly reviews related work.Section 3 presents the traditional problem formulation and the methodology to solve it.Section 4 contains the results analysis and the proposed model's validation.Finally, in Section 5, conclusions are presented.

Related Works
The work proposed in this article represents a significant advance in electro-mobility infrastructure planning.By considering actual transport variables at the trajectory level, the current job has developed a heuristic model that can plan and allocate resources to deploy electric vehicle charging stations.The model considers the complex interplay between vehicular flow, road capacity, and trajectories and uses advanced optimization algorithms such as multi-commodity flow problems and graph theory to provide finite and scalable solutions over time.
One of the key contributions of this work is its potential to become a computational model and infrastructure manager for planning and resource allocation in electro-mobility.The scalable model has been validated through iterative processes that deal with contingency concepts to predict the maximum density of vehicular traffic in road networks.Moreover, the model provides hourly analysis of vehicle flow and allows for the identification of maximum and minimum demands experienced for each charging station, which contributes to coordinating and managing battery-charging processes in electric vehicles.Overall, this work represents a significant step forward in developing sustainable transportation systems and has important implications for the future of electro-mobility.
This paper studies the optimal deployment and allocation of CS resources to satisfy demand at the trajectory level by considering vehicular flow characteristics with heterogeneous users.This paper is closely related to three lines of research: optimal CS location, Resource allocation for EV charging at the trajectory level, and ICSEC planning.This section summarizes the most current and relevant work on the proposed study.

Assignment and Location Problem
The location of energy supply centers for EVs was developed in [13,14] by considering the routing problem and EV autonomy as the decision variables.Table 1 sets out brief characteristics of designed vehicles for urban usage of the L7e-C category, in which the average range is observed to be 123 km.The decision problem is determining the optimal supply for charging infrastructures, the locations and users' allocation to each CS, and the optimal dispatch policy.In one paper [15], the problem is expanded to a dynamic planning model by considering itineraries by EV usage.Recently, in [16], the location problem was studied by analyzing supply and demand, including the driver's psychological component, and in [17], in addition to deploying CS, it is suggested to involve the integration of distributed generation sources.
The authors showed that integrating distributed generation helps relieve unexpected charges on existing distribution electrical grids, which have generally been operated in cities for several years.However, it has not been possible to demonstrate a rigorous analysis considering the road capacity in each section and vehicle flow analysis that defines city traffic in decision variables to allocate and locate CS in electric mobility infrastructures.

Resource Allocation under Optimal Criteria
Both exact and approximate solutions have been applied to solve the resource allocation problem.Accurate solutions involve integer linear programming problems, while approximate solutions involve heuristics and metaheuristics development.Table 2 provides a taxonomy of related studies with a summary of the problem statement, objective function, main constraints, solution methods, and other fundamental characteristics to address the problem of allocation and CS display.
Recent developments focused on solving planning problems based on a convex multiobjective model, which will depend on the number of variables and constraints considered.Among the issues that need to be adequately considered are those features that bring the problem closer to real-world situations, such as the heterogeneous vehicular fleet and the multiple CS options, which are used to route the vehicular flow to ensure continuity of coverage for EV battery recharging.

Planning of Charging Infrastructure for Batteries in Electric Mobility
The technological advancement and the constant effort to introduce new concepts of electric mobility on a large scale have encouraged the scientific community to develop products that enable the reduction of CO 2 emissions by conventional private and public transport.Today there are multiple electric mobility alternatives, such as bicycles, cars, and public transport buses.However, it has not been possible to introduce EVs massively into the land transport system because of variables that do not make large-scale purchases attractive to potential users in urban and rural areas.These unattractive variables to consumers could be limited autonomy, long charging times, battery life, high costs, and lack of charging infrastructure for EVs.This article focuses on the problem of planning and deployment of CS to support the mass introduction of electric vehicles into the vehicle fleet.The spatial-temporal collaboration between CS, allocation, and optimal routing based on bidirectional flows, undoubtedly, gives place to ICSEC.In one paper [27], the authors proposed tackling the problem by modeling the hybrid vehicle traffic with different charge terminals and observing EV driving range constraints.Furthermore, the mathematical formulation and solution are addressed by linear programming.
In another paper [28][29][30][31][32], the problem of allocation and deployment of CS is addressed by partially considering several aspects which guide the achievement of robust results; these aspects are limited capacity, probability of trajectories, travel time, flow conservation, charge time, spatial, temporal collaboration, terminal type, and traffic events.
Planning, CS deployment, and selecting the charging terminal are important because EVs present technical constraints such as maximum driving ranges and limited payload.Airports have been classified into three types, depending on the charge current which EVs can withstand.
Table 3 briefly overviews the different types of existing terminals and their various applications.It is well-known that as the charging current increases, the charging time decreases, and the types of terminals differ.Type I and II application are in residential areas, and type III is exclusively a public/private service.Consequently, this article is essential because it covers the planning problem in its integral form at both levels, as explained in the previous paragraphs.Remember that to dimension CS, it is necessary to model the vehicular traffic at the trajectory level and the user's need for battery recharging.Therefore, it has to consider ICSEC an essential actor during the planning process.
To look for integrated solutions that safeguard the operability of existing electricity grids, according to their lifetime and the fact that many of them have reached their nominal capacities, they could be affected by the massive introduction of EVs in the traditional electricity distribution system.
It will allow us to decide whether or not the existing networks can support new charges, thus making it possible to choose to use the existing equipment or whether it is necessary to install new distribution circuits to satisfy the demand for battery recharging in electric mobility, to guarantee the quality, continuity, and safety of the electricity service.In [33,34], a method of clustering and constructing minimum spanning trees is proposed as an indispensable tool in electrical equipment's optimization and planning processes.

Problem Description
This section describes the problem through the stochastic analysis of vehicular traffic at the trajectory level and allocating resources for ICSEC.Furthermore, the strategy and methodology proposed are presented in this paper, followed by the formal problem statement.

Stochastic Analysis of Vehicular Traffic and Resource Allocation at ICSEC
Service stations in the field of public-private transport have an important role because they supply primary energy to provide autonomy to the electric mobility vehicle fleet.The massive introduction of EVs in daily traffic has to be planned, which implies having heterogeneous vehicle networks where the main actors are conventional vehicles and EVs from the public and private sectors.
Most research focuses on users and reducing travel times by reorganizing vehicle traffic according to road capacity and considering heterogeneous traffic networks to allocate resources.A fundamental detail is that the demand for vehicle charge is not focused on the nodes (electric or conventional vehicles) but on the flow of the vehicular traffic network.Consequently, the research objective is to solve the problem of location and allocation of resources to CS by considering heterogeneous vehicle density.
Therefore, the vehicular traffic flow can simulate the demand generated by EVs displayed on a heterogeneous road network.It is assumed that installing a service center for EV battery charging cannot satisfy the induced demand due to the limited capacity of the CS on a road section.This limitation will cause critical traffic flow because of high coincidence factors due to battery recharging requirements.Consequently, battery capacity is fundamental to CS deployment and resource allocation.Therefore, it will obtain the following considerations: (i) if there is a CS between a start and end point of a trajectory, the EV's autonomy must be greater than the length of the route, and (ii) if there is only one CS in the route, the EV's autonomy must be more than double the route length.
An important detail already exposed is that the algorithm must obtain knowledge of the multiple options at the trajectory level that an EV can take in terms of spatial and temporal constraints (distance, trajectory, and flow variability).In addition, it considers that the maximum demand for battery recharging in a given area will be provided by the maximum vehicle density traveling in a given time interval.This maximum demand or traffic density occurs when the vehicle spacing and the speed with which the vehicle is traveling are the minimum.
Based on microscopic models for vehicle density analysis, it studied the maximum allowable density in a georeferenced area by considering vehicle spacing and speed variables.Using a random vector under Poisson distribution, it simulates the vehicular traffic density within a time interval considering maximum traffic density at minimum speeds and spacing.The objective is to select and allocate resources to CS to maximize traffic flow and provide users with coverage in a given area.
The deployment planning of CS can consider the location of the different charging infrastructures and their influence on the power quality, security, and economics of the system's operation.Consequently, it is defined as a typical multi-objective optimization problem.This article will solve the planning problem by considering a multigraph in levels: (i) traffic and (ii) distribution networks.These categories can overlap, meaning two nodes of each type can be located in the same area.In level 1 of the traffic network, the vehicle flow is characterized.Therefore, in level 2, the capacity of the charging station is studied according to the number of EVs that can be supported.

Proposed Strategy and Methodology
This paper proposes a theoretical model of location and assignment of CS with traffic flow restrictions in heterogeneous scenarios.The main objective is to optimize resources using the maximum vehicle flow capacity on each road section in georeferenced areas.A non-convex combinatorial problem is proposed under a mixed integer linear programming structure.Consequently, it will use the column generation algorithm to solve the multiproduct flow problem frequently used in MILP.
The column generation algorithm divides the problem into two stages (i) restricted master problem (RMP) and (ii) sub-problem.Solving RMP with the minimum number of variables obtains dual costs for each restriction master problem; this solution is used in the objective function of the sub-problem and solved.If the accurate value of the sub-problem is a negative reduced cost variable, this is added to RMP and is again solved.This process is iterative, and the stopping criterion happens when the objective value of the sub-problem is greater than or equal to zero.When this happens, the restricted master problem can be considered optimal.
At the first level, CS display and location are performed by assuming that a driver will choose the shortest path from the starting point to the destination by observing distance and traffic flow constraints.These restrictions will reduce to an optimal number the vector length of candidate CS sites located at strategic locations, such as shopping centers, parks, car parks, customer service centers, and conventional charging stations.
At the second level, the capacity of the CS must satisfy the daily EV charging demand, which is determined by the charging time of the simulated vehicles and the number of EVs that it may be able to handle in a specific time interval.Different locations and capacities of the CS will cause other impacts on the distribution grid.An additional detail is that the Dijkstra algorithm is used to find the shortest path between the multiple ways generated in the function of the vehicle density.
Four main groups are proposed to guide the methodological development of this article: (i) scenario analysis, (ii) process, (iii) simulation, and (iv) solution.Figure 2 shows a solution to the transport and resource allocation problem flowchart.In scenario analysis, the model suggests having full knowledge of the case study, which means that the model knows the coordinates of the electrical equipment (charging stations) and, based on the positions of each vehicle, the maximum vehicle density at the trajectory level is known.Once the scenario has been characterized, the experimental procedure defined by the process and simulation detailed in Figure 2 is carried out.This process explains the actions, considerations, and variables that enter the simulation.After the iterative process, the proposed heuristic returns the hourly local optimum in the function of the incoming vehicular traffic density variables.It used Matlab software for the simulation.
Finally, this article considers basic research because of its experimental nature by implementing simulation processes.It has a quantitative approach to determining the allocation of resources and minimization of economic and social impacts.This model is exploratory; therefore, it is essential to analyze different study cases.

Problem Formulation
The product flow problem considers two necessary restrictions: (i) travel demand, which means that all products or EVs must be transported to their destinations or CS, and (ii) the capacity of each intersection of each road section; this means that the traffic flow at each intersection does not have to exceed its capacity.Therefore, the multi-product problem would appear to be a combination of several single-product problems, given both the limited roadway capacity and the diversification of vehicle density.Some considerations and the formal description of the problem statement are presented in this section.
Furthermore, the interaction between EVs denotes multiple complexities to solve the single-product flow problem independently.Consequently, the demand generated by each EV is related to various supply options for this demand, and its origin-destination node pair, which refers to the position of the EV and the CS, respectively, can identify it.The located CSS provides the battery recharging supply.
The optimal collection of CS is selected based on spatial-temporal collaboration, which involves minimum displacements by observing traffic flow constraints.Moreover, it assumes that the shortest or available route depending on the nodal traffic flow and capacity at the trajectory level can acquire the essential product supply.The variables are detailed below in Table 4. Consider G(ϕ, V, A) a directed graph where ϕ corresponds to the autonomy constraint of an EV towards a CS; the distance computation is performed with the haversine equation and it is given in km.The vertices set V represent the multiple nodes located in the study area, such as electric vehicles (ν), charging stations (ε), the intersections (ζ); therefore, ν, ε, ζ ∈ V and the displacement is between two nodes origin-destination ν u,v , ε u,v , ζ u,v ∈ V and the set of edges u, v ∈ A represent the origin and destination, respectively.A allows the construction of a specific topology; that is, it relates the multiple nodes existing in G leading to the creation of spanning trees at each level.On each track section (ι), there is a certain number of vehicles (α).Each track section is limited in capacity ( ), and it is formed between two vertices in a set of corners (β).At each β, it reflects a demand or cumulative cost (W β ).Consequently, the multiple paths (S 1,2,3... ), where S ∈ V, is formed from a set of intermediate nodes β, which, dimension the total flow f β by assigning a weight W β on each path.A track section can be either unidirectional or bidirectional β i,j = β j,i .In other words, {β ∪ S} ∈ V. Consider κ ⊆ V a subset of candidate sites where κ / ∈ β, ν because the selected positions to locate CS are parks, shopping centers, public car park.Moreover, u = {i ∈ V : ∃β ∈ S| i = 1, 2, 3 • • • } denotes the subset of source nodes in a trajectory and v = {j ∈ V : ∃β ∈ S| j = 1, 2, 3 • • • } it denotes the subset of destination nodes, where m is the length of the vector {ν ∪ ε ∪ ζ ∪ β}.In other words, the different trajectories (products) are constructed only from the charge nodes β to ζ.This is possible, since ∑ m i=1 ∑ m j=1 β i,j , ∀(i = j), ∈ Z + .Finally, consider δ κ the cost to install a node ζ j where j ∈ κ and it is n j the maximum number of vehicles to cover the capacity ψ of each CS.
Consider M ⊆ A subsets of assigned edges that define the enabled roads in the whole studied scenario.Therefore, it isM = β i,j , β i ∈ u y j ∈ κ.Consider L ⊆ A a subset of edges of a path L i , L j where i, j ∈ {κ ∪ β} with this it guarantees that there are no arcs from other families.Next, it defines the set of binary variables of resource selection.
• φ j location variable ε, is 1 when it is found any node that is part of the set j, j ∈ κ, 0 in any other way.• θ k subset of assigned arcs, if the arc k ∈ M it is used, otherwise 0. • τ j e path enabled, is 1 if the edge e ∈ L and j ∈ κ, otherwise 0.
Equations ( 1) and ( 2) are used to calculate vehicles' total and partial concentration in a specific area.The Equation (3) expresses the current flow rate on a road section and represents the frequency at which a certain number of vehicles pass through at a given time.
Equation (4) calculates the average number of vehicles that may be on a lane in the function of the cross-sectional length of the road and the spacing s between cars and Θ it is an annual rate of increase under a discrete Poisson probability distribution.
Finally, with Equation ( 5), the number of vehicles as a function of cross-sectional length and spacing is calculated.It assumes that the spacing between cars is uniform and that the maximum capacity of a road-cross section occurs at the minimum distance s between vehicles and minimum speed (V).It is essential to mention that a microscopic model has been used for the traffic flow analysis, which considers the spacing s and the individual speeds (V e ).Therefore, the multi-product transport problem can be formulated as follows: Objective function: Subject to: where Equation ( 6) corresponds to the objective function; consequently, Equations ( 7) and ( 8) guarantee that the origin and destination of the route are from service nodes to demand response.Equations ( 9)-( 11) correspond to the corresponding routing's vehicle concentration and electric vehicle autonomy restrictions.Equation (12) ensures that dual variables are searched and aggregated.Finally, the binary variable is declared with Equation (13).Furthermore, Q k = [q 1 , q 2 , . . ., q m ] T where k = 1, 2, . . ., t and represents the flow rate of each basic product.
The accumulated capacity and autonomy vector of the EV is represented by The adjacency matrix between vertices and edges is represented by G = [G ij ] m×m .The length of m × m is a function of the total number of vertices in the scenario.Finally, . ., t is defined as follows.
any other case ( 14) Column Generator This article used the column generation method.This method allows us to solve extensive linear programs.The idea is to take advantage by generating variables that have the potential to improve the objective function.The algorithm process is as follows.We first consider the following assumptions: (i) the variables are non-basic, (ii) The variables will assume a value higher or equal to zero in the optimal solution, (iii) an initial subset of variables must be considered to solve an initial problem, and (iv) the problem can be divided into the primal and dual problem.
The primal problem is solved with a subset of variables considered; with this solution, it obtains dual costs for each restriction of the primal problem.It assumes that χ k * is the optimal value of the dual objective function, which was solved by the revised simplex method; if the real double deal is a negative reduced cost such that χ k * < 0, ∀k = 1, 2, . . ., t this variable is added as a column to the primal problem, and it is again solved iteratively.
When χ k * ≥ 0, ∀k = k 1 , k 2 , . . ., k p , the algorithm stops, and it can conclude that the primal problem is optimal.
For each linear problem, there is a problem that is solved in parallel; the last one is known as a dual problem.We consider the following issue to solve the primal and dual problems.
Subject to: x ≥ 0 The linear problem is defined according to Equations ( 16) and ( 17) with a dual counterpart.
wA ≤ c (20) Equations ( 19)-( 21) have the same terms, except for the term w, which is a variable that is assigned to Equation ( 17) of the first formulation, where w = C BV B − 1, which corresponds to an iteration of a revised simplex.It is known as a dual variable arrangement.The variables C BV y B − 1 correspond to primary variable coefficients and connectivity matrix that relates vertices and edges, respectively.Consequently, the problem illustrated in Equations ( 6)-( 12) can be separated into two primal and dual issues, respectively.The equations defining each primal and dual problem are presented as follows.
Primal-Primary Problem Objective function: Subject to: The new set P is used in each iteration to identify the generated routes, and Ψ p is a partial solution to the primal problem.Equation ( 22) is the objective function of the main problem, and Equations ( 23)-( 25) are vehicular concentration and autonomy restrictions, respectively, that permit restricting the solution space to routes that only fulfill a little value.With Equation ( 26), a non-zero solution to the second problem is guaranteed.
Dual-secondary problem Objective function: Min: Subject to: Equation ( 27) is the objective function of the secondary problem, so χ k * is the reduced cost variable, which, as long as χ k * = 0 could find dual variables which are inserted into the main problem.Equations ( 28) and ( 29) are restrictions that enable source and demand nodes so the model can search for the best route.In Equation (30), the coefficient r is selected from the set S that will allow verifying the cost of the negative reduced variable.Finally, the binary selection variable is declared with Equation (31). Figure 3 presents the flowchart for resource location and allocation required to ICSEC.

Analysis of Results
This section presents the results obtained by the heuristic model MP-LRAC for creating ICSEC.Table 5 shows the simulation parameters, where the need to consider input variables through heterogeneous and scalable vehicular traffic analysis in georeferenced areas is noted.
The microscopic model was used to simulate the input variables, which allows us to study the relationship between speed, length, and spacing.It has tried to recreate scenarios close to reality to determine the average number of vehicles circulating in a given area.Therefore, it will be possible to deploy charging stations for electric cars and allocate resources for the energy supply by recharging batteries in electric mobility.Information, such as connectivity matrix, track direction, and weights assigned for each section in the optimization model, are essential.Consequently, the model's longitudinal relationships are well-known from an intersection i to a meeting j.Case study A is illustrated in Figure 4, and here it can find the topology analyzed in this article.
Considering the actual data of case A will allow us to approach reliable solutions.Moreover, it will permit an understanding of the maximum load capacity admitted in a given time interval in this study case.It means the maximum traffic capacity on each road section is determined according to the topology (longitudinal characteristics of each road section, road direction, spacing, and vehicular speed).
The maximum number of vehicles circulating in study case A is determined according to the complete vehicle density.Figure 5 presents the variation in vehicular concentration in the function of different incremental rates during a 20-year analysis period.The total rate corresponds to a theoretical value, which the vehicular network designer can adjust.Its value will depend on studies and estimated projections in each area of its different central governments.
This analysis of the increase in the vehicle fleet was performed over 20 years with four annual rates Θ as shown in Figure 5.There is a technically linear trend in the first ten years of analysis, which does not occur from the 11th year onwards.The literature suggests that discrete Poisson probability has recreated the scenario's randomness.
The slopes cease to be continuous as the annual rate of the vehicle fleet increases, and they adopt exponential trends in their addiction.Therefore, a first reflection is that it is necessary to generate control policies that mitigate for this unfavorable impact by not obtaining control over the increased vehicle fleet.
-120.625 -120.The number of heterogeneous vehicles deployed in each case in Figure 5 responds to the maximum capacity at the trajectory level that involves variables explained in previous paragraphs, such as cross-section, spacing, and speed.Consequently, it can be inferred from Figure 5 that as the average speed with which vehicles circulate on each section of the road increases, the inverse happens with the vehicle concentration.It means that as the speed increases, the vehicle density decreases.
This phenomenon is because the flow rate is inversely proportional to the speed.The purple color represents the empirical values that indicate the number of values above and below the average value of the three speeds analyzed from 11 to 20 years.It means that approximately 95% of data are above and below average, and their values are visualized in Figure 5. Therefore, speed, geometric road layout, spacing, and annual growth rates of the vehicle fleet are considered in this article as fundamental variables to achieve results that satisfy near-optimal solutions in real scenarios.
Figure 6 illustrates the number of vehicles distributed over 24 h for case study A. The trend in Figure 6 corresponds to the weekly average vehicle traffic behavior under standard working day conditions.The proposed model in this article is evaluated over four years from year zero, and an incremental rate of 3% will be considered.The total rate is taken as a reference from the literature to adapt this variable to any required reality or case study.
The hourly variability of the number of vehicles that circulate in a given area responds to a normalized random vector by considering the original trend (year zero) and the annual rate of increase.Moreover, Figure 6 shows that the highest vehicle density is distributed from 07:00 to 20:00 h.Consequently, the charging stations should be designed to meet the demand during peak traffic hours in a given area.1.9 2.9 3.9 4.9   Consequently, in Figure 7, we illustrate the estimated number of vehicles counted in a specific hour from 0:00 to 24:00 h through the frequency graph.It is correct to think that the number of cars that circulate in a particular time will follow the base mobility pattern curve, Year: 0. Finally, Figure 7 presents the zones of highest vehicle concentration illustrated by the edges, which in turn are randomly distributed in case study A.

Latitude Longitude
Frequency (q) (veh/hrs) vehicular density PDF The adjacency matrix allows us to understand the data distribution within the case study.Figure 8 presents two techniques for reordering the sparse adjacency matrix: (i) Cuthill-McKee inverse algorithm and (ii) Column Permutations.In the solid blue color in Figure 8, the inverse Cuthill-McKee algorithm is illustrated, which has a clear tendency to group its non-zero elements closer to the diagonal while preserving the special relationship from each edge (i, j) to its nodes i = j and j = i.The column permutation algorithm illustrated in Figure 8 using a heat graph permits reordering columns of the sparse-type matrix A in non-decreasing order of non-zero count.This metric allows us to visualize the different magnitudes of vehicle density that will be considered to evaluate the model and the relationship between each node in the scenario.
Therefore, Figure 8 shows that the input data are stochastic and that their origin is the result of a microscopic analysis that allows us to understand the interaction of the different variables involved in vehicular mobility, such as spacing, speed, the topology of the study area, and lengths of the vehicles.It has enabled us to bring the problem closer to the stochastic reality of land transport mobility by considering its possible trajectories within a geolocated area.Hence, it uses the microscopic model to solve PM-UARC as a directed graph to model the traffic rate and the maximum vehicle concentration at the trajectory level.
The discrete Poisson distribution generates different case studies based on a typical curve of vehicle traffic operation during working hours from Monday to Friday.This base curve is considered to be year 0. The scalability analysis of the vehicle traffic evaluated in the proposed model's performance is based on it, and its results are shown in Figure 9.      9 illustrates the algorithm performance evaluated in a scalable scenario from year 0 to year 4. Figure 9a shows the machine time to solve the primary problem.It is easily observed that it depends on the time of day; the machine time spent to solve the primal problem is variant; this response to the stochasticity of the vehicular flow in each time interval.
The linear problem is constructed as the primary concern.In Figure 9a, the time trend of the machine is a function of the annual growth rate of the vehicle fleet.It is increased with a minimum of 3 min at year 0 and a maximum of 4.7 min at year 4. In Figure 9c, the observer's machine time is presented, which has the function of deciding the growth of the column generator; this means they are the one who verifies whether the cost in terms of objective function decreases to be able to add it as a solution.Next, if a non-zero vector exists, it is added to the solution.Once added, the algorithm repeats the iterations from the primary and secondary problems, and the observer makes decisions based on the negative variables.Therefore, the time of the machine is higher than the time spent on the secondary problem, as illustrated by comparing Figure 9a,b.
A curious thing to note in Figure 9d is that the cost of the objective function in each scenario tested maintains the same trend; this means there is no data deviation from year 0 to year 4. The movement in Figure 9d responds to the immutability of the scenario topology.Finally, Figure 9d shows that the objective function is inversely proportional to the hourly characteristic curve that defines the vehicle density.
This inverse behavior of the objective function with the hourly distribution of vehicular density is linked to the operating costs; this means that the design and sizing of charging stations for electric vehicles must be according to the maximum demand that occurs from 07:00 to 20:00 h on ordinary days from Monday to Friday.Therefore, the objective function has a minimum cost in high-demand hours because it manages to justify the charging station sizing.The accurate function cost is maximum, whereas traffic demand is minimum.It is because, in those low-demand hours, the charging infrastructure is not used, so there is an oversizing in low-demand hours.
Table 6 presents the resulting metric for case study A. In each column, respectively, it can be seen the ID of the charging station, the geographical coordinates in latitudelongitude, the vehicle concentration, and the number of vehicles to be served in the hour of maximum demand.The charging station has ID 4 experiments with a complete vehicle charge of 48.71 per 1 h.
According to what has been analyzed in the preceding paragraphs, the demand in year 0 occurs at 19:00 h, with 1151.34 vehicles distributed in the geographical area under study.The distribution of vehicle charges to be handled by each charging station during the hour 19:00 in year 0 is presented in Table 6.A percentage pattern allocates vehicle charges for the day's remaining hours (see Figure 6).However, in this case, it has considered the maximum vehicle demand in a time interval to estimate the maximum power installed at the different charging stations.
The installed power of a CS should be directly proportional to the number of vehicles estimated to be met in the peak hour of vehicle charging.For example, if it considers the time 19:00 in Figure 6 and looks at Table 6, which illustrates the vehicle charge distribution, it observes that the charging station with ID 4 should serve 48.71 vehicles.If it assumes that all users want to charge their cars simultaneously and that the CS topology has ultra-fast charging terminals with 350 kW power, for each terminal, it should assign a design power for the EC with ID 4 of 17.05 MW.
However, it is essential to remember that the model performs the analysis for 1 h; it will determine where the maximum vehicle demand occurs by considering the topology of the geolocated area involved; this means that the 48.71 vehicles will be distributed over 60 min.Here, a fundamental detail must be analyzed in the design of the CS topology, which will necessarily involve the battery charger technology.It means that the charging times and terminals are studied to satisfy the 1-h high demand for charging batteries for electric vehicles.
The model's adaptability to any scenario will be verified by analyzing case study B. This study case is in Arroyo Grande County, San Luis Obispo County, California, United States.It will simulate scenarios where the CSS goes out of operation for any condition and analyze the redistribution of the vehicle flow.Furthermore, it will explore what happens with the different vehicle flows in each hourly interval with a growth rate of 3.0%, as well as verify contingencies of up to N-2.
In addition, how it affects the sizing of the charging stations in EV charging infrastructures.Figure 10 presents the geographical area in which Case B will be analyzed, where nine EV charging stations are identified with ID tags.The vehicular density analysis has been studied with microscopic examination.An additional detail visualized in dark green is the bi-directionality of each road section and how this is related to each intersection.Moreover, Figure 10 identifies the costs directly associated with the vehicle density on each path.Consequently, a case study involving traffic variables such as vehicle spacing and travel speed in vehicular networks is presented in Figure 10.These variables are entered through a random vector with the Poisson distribution.
Figure 11 presents the graphical solution to the transport problem.The magnitude of the vehicle density generated at each link is visualized in black; in red, at each intersection, the resulting vehicle density can be easily identified, which is the sum of the flows that enter each node minus the flows that leave.A fundamental condition for the convergence of the model is that the supply centers with a negative sign cover the demand with a positive signal at each intersection.This condition is critical to maintaining the nodal equilibrium criterion, which satisfies the equations in the mathematical problem statement.The red negative sign indicates the flows that will have to be served by each charging station; let us remember that the negative sign indicates that it is a supply node (charging station).
A fundamental criterion for equipment sizing in planning algorithms is considering the worst-case scenario.This consideration will allow us to approach the optimal sizing of the resources required in the planning processes.Consequently, Table 7 presents the result of the analysis case in Figure 10, where the case of the highest conflict is considered.
The highest conflict occurs at maximum vehicle density moving through the geolocated area.These extreme conditions occur at minimum travel speeds with minimum vehicle spacing distances.It means that vehicle density in a geolocated area is inversely proportional to vehicle speed and spacing.Column one in Table 7 presents the charging station index; column two illustrates the resulting vehicle concentration to be served by each charging station.
The sign in column two is negative because the model considers it a node contributing to the transport system to form the balance equations.The supply nodes must cover the sum of the partial flows observed at each intersection.Consequently, with a positive sign, it denotes those consumption nodes, and with a negative character, it identifies the nodes that satisfy the equilibrium equations of the multiple-product problem.
In addition, in Table 7, the third column shows the number of vehicles each charging station will have to serve in critical conditions of high vehicle density, without considering contingencies, in a time interval where the speed and vehicle spacing variables are minimal.
Finally, the last row of Table 7 identifies the maximum number of vehicles in an instant of time according to the scenario topology, which in this case is 348.34 vehicles as the maximum admissible density in the area considered for the analysis.Consequently, with the number of cars each charging station serves and the maximum traffic density at the trajectory level, it can size and foresee the necessary resources for constructing charging infrastructure in geolocalized areas.Moreover, Table 8 presents the optimal sizing of charging stations by considering events or contingencies (N-1 and N-2) generated in a geolocated area from normal operation conditions (N-0).These contingencies respond to removing one or two charging stations at maximum vehicle density.The out-of-service charging stations are randomly selected.Column two in Table 8 shows the rates of charging stations that go out of service, and columns 3 to 11 present the flux density at each charging station.
Considering the contingency analysis for out-of-service charging stations, the traffic density to satisfy the battery-charging requirement must be routed or rearranged with the charging stations enabled.The experience with charging station ID 1 manifests the increased traffic density generated by operating one or two charging stations.If one looks at column three in Table 8 in normal operating conditions, the vehicle density to be served by the charging station ID 1 is 29.95 ≈ 30 vehicles in a time interval when the vehicle density is maximum.When the model experiences the operation exit of the charging station with ID 4, the number of vehicles to be served given contingency N-1 is 48.22 × 48.If the charging infrastructure experiences the operation exit of two charging stations with ID 3 and 7 of contingency N-2, the traffic density to be served is 53.81 × 54 vehicles in a given time interval.
The increase represents approximately 61% and 79.7% of the different traffic densities that will have to be served when the charging infrastructure disconnects charging stations in N-1 and N-2 contingencies, respectively.Likewise, it will perform the analysis for the charging station with ID 5, which presents an increase in vehicle density of 95.4% in contingency N-1 and an increase of 167.6% in contingency N-2 that it will have to serve if one or two charging stations are taken out of operation.
If one pays attention to detail, the charging station with ID 5 in contingency N-2 experiences some increases in vehicle density that exceeds 100% in the regular operation in contingency N-0.It is understood that the algorithm reorganizes the vehicle flow in the function of a cost and capacity vector for each road section, satisfying the mathematical variables of the linear equation; however, the balance equations are fulfilled for each iteration.The charging stations as a set will give rise to the charging station infrastructure.Therefore, if the designer of the charging infrastructure knows which percentage would increase or reduce the vehicle density to be served by each charging station, they will be able to make intelligent decisions to foresee the electrical, civil, and logistical equipment required for the location of charging infrastructures in electric mobility to guarantee coverage in the road network.To determine the sizing of the charging stations, consider the empiric rule to assess the magnitude of values within a band around the average in a normal distribution with a width of two times the standard deviation.Consequently, the green color in Table 8 indicates the optimal sizing for each charging station for electric vehicles.
By multiplying the cumulative length in kilometers of the road network times the vehicle density considering the empirical sizing criteria, one can obtain the number of vehicles each charging station will have to serve in each event studied.Furthermore, Table 9 allows us to verify how much the sizing of the charging station should increase given a circumstance.Consequently, the information presented makes planning the charging infrastructure for electric mobility in smart cities possible.Once the number of electric vehicles that will circulate through each scenario studied has been identified in detail, the energy required to implement the charging infrastructure is determined.For the load forecast, it is assumed that the input power of each charger is 52.63 kVA with a power factor of 0.95, so the active power consumed by each charger is 50 kW.In addition, the charging time is assumed to be 0.5 h. Figure 12 shows the energy demand required by the charging infrastructure for electric vehicles in cases A and B studied.It is noticeable that as the service area increases, the energy needed to start up the infrastructure rises proportionally.Figure 12 follows the theoretical trend of the hourly distribution of vehicular traffic.As can be seen, there are at least four peaks of maximum demand at 07:00, 14:00, 16:00, and 18:00 h.During these hours, the vehicle density increases and thus increases the probability that an electric vehicle needs to charge the battery of the electric car.The expected number of electric vehicles and the minimum total power required in the charging infrastructure is illustrated in Figure 13.The average number of charging terminals needed for each service station for cases A and B are 6 and 7, respectively.This slight difference responds to the geographical size of each scenario, i.e., the larger the coverage area with charging terminals for electric mobility, the greater the vehicle density and, consequently, the greater the probability that electric vehicle users will require charging their batteries.The ID of each station is represented in the x-axis of Figure 13.In addition, as the number of vehicles increases, the need for energy demand at each charging station increases.This is because more charging terminals will be needed to satisfy the demand for battery recharging in considerable waiting times to guarantee greater comfort to the users of the charging stations.The charging station with ID 4 has six charging terminals and requires an installed power of 297.71 MVA.The following equation was used to calculate the minimum number of charging stations required.ET = VT * CT PO (32) ET is the number of charging terminals for electric vehicles, VT is the vehicle traffic, CT is the terminal charging time, and PO is the charging station power output.For example, from Figure 13, we extract the information of the number of vehicles in case study A. Let us take as a reference a charging station with ID 4. The maximum number of electric vehicles traveling is approximately 600, then f rac 600 Veh * 0.5 h 50 kW = 6 by requiring at least six 50 kW charging stations 6 * 50 = 300 kW or 315.9 kVA if 0.95 is considered to be power factor.Therefore, the proposed model has demonstrated the ability to adapt to different case studies that allow us to analyze contingencies, which are very useful at the time of sizing because it will permit us to consider the electrical, civil, and logistical equipment to lead to the construction of adequately planned charge infrastructures.

Discussion
The relationship between the deployment of public and private charging stations with economic conditions and policy guidance presents itself as a complex challenge, as the adoption of electric vehicle charging technologies depends on a combination of factors, including tax incentives, vehicle acquisition costs, charging infrastructure, and policy guidance towards the energy transition.Government policies and regulations are essential in encouraging investment in charging infrastructure.Still, it is also vital to consider charging stations' long-term profitability and financial sustainability, especially in the context of economic and political fluctuations.Therefore, an integrated approach involving key industry, political, and societal stakeholders is needed to achieve effective and sustainable deployment of EV charging stations worldwide.
The growing demand for electric vehicles is increasing research and development of more efficient batteries and faster and more accessible recharging technologies.However, it is necessary to remember that the mass adoption of electric vehicles poses complex technical and economic challenges that require coordinated and collaborative solutions from industry, governments, and society.It is essential to promote a transition to sustainable electric mobility in a strategic, planned, and sustainable manner over time to ensure an effective reduction in greenhouse gas emissions and a positive impact on air quality and public health.

Figure 1 .
Figure 1.Proposed architecture for the deployment of electric charging stations.

Figure 2 Figure 2 .
Figure 2 reveals that the georeferenced data are extracted from a free database published in OpenStreetMap.

Figure 3 .
Figure 3. Flowchart for solving the Multiple-Product Problem in Location and Resource Allocation by considering Capacity (MP-LRAC).

Figure 4 .
Figure 4. Case A, Grover Beach with an area of 2.13 km 2 , a city in San Luis Obispo County, US state of California.

Figure 5 .
Figure 5. Case A: Microscopic traffic flow analysis using discrete Poisson probability.

Figure 6 .
Figure 6.Vehicular traffic scalability with the incremental rate 3% for case study A.

Figure 7 .
Figure 7. Vehicle density distribution located in geolocated scenario case A.

Figure 8 .
Figure 8. Case A: Distribution of connectivity matrix and its nodal relationship.
Machine time to solve the primary problem.
Machine time to solve secondary problem.
Cost of the objective function.

Figure 9 .
Figure 9. Case A: Performance of the algorithm considering scalability of 4 years with a growth rate of 3%.

Figure
Figure9illustrates the algorithm performance evaluated in a scalable scenario from year 0 to year 4. Figure9ashows the machine time to solve the primary problem.It is easily observed that it depends on the time of day; the machine time spent to solve the primal problem is variant; this response to the stochasticity of the vehicular flow in each time interval.The linear problem is constructed as the primary concern.In Figure9a, the time trend of the machine is a function of the annual growth rate of the vehicle fleet.It is increased with a minimum of 3 min at year 0 and a maximum of 4.7 min at year 4. In Figure9c, the observer's machine time is presented, which has the function of deciding the growth of the column generator; this means they are the one who verifies whether the cost in terms of objective function decreases to be able to add it as a solution.Next, if a non-zero vector exists, it is added to the solution.Once added, the algorithm repeats the iterations from the primary and secondary problems, and the observer makes decisions based on the negative variables.Therefore, the time of the machine is higher than the time spent on the secondary problem, as illustrated by comparing Figure9a,b.

Figure 10 .
Figure 10.Case B: Weights assigned to each track section with microscopic analysis in an area of 0.9887 km 2 .

Figure 11 .
Figure 11.Case B: Graphical solution to the multiple-product problem.

Figure 12 .
Figure 12.Energy consumption from the massive introduction of electric vehicles.

Figure 13 .
Figure 13.Case A. Energy consumption and density of electric vehicles per expected charging station.

Table 2 .
Taxonomy of related studies.

Table 4 .
Notation used in this article.

Table 6 .
Case A: Maximum partial demand in each CS at hour 19:00 with 32.97 km of road network in contingency N-0 corresponding to year 0.

Table 7 .
Case B: Maximum partial demand in each CS with 9984 km of road network in contingency N-0.

Table 8 .
Case B: Optimal sizing of charging stations for EVs based on multi-product flow problem.

Table 9 .
Maximum partial demand at each charging station with 9984 km of road network given a contingency.