Next Article in Journal
RETRACTED: Gładyszewska-Fiedoruk et al. Simulation and Analysis of Various Ventilation Systems Given in an Example in the Same School of Indoor Air Quality. Energies 2019, 12, 2845
Previous Article in Journal
Optimal Daily Trading of Battery Operations Using Arbitrage Spreads
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Sizing of Electric Vehicle Charging Stations Considering Urban Traffic Flow for Smart Cities

1
Master of Electricity Program, Universidad Politécnica Salesiana, Quito 170525, Ecuador
2
Smart Grid Research Group, Universidad Politécnica Salesiana, Quito 170525, Ecuador
3
Inclusive Education Research Group, Universidad Politécnica Salesiana, Quito 170525, Ecuador
*
Author to whom correspondence should be addressed.
Current address: Postgraduate Department, Girón Campus, Av. 12 de Octubre N 23-52, Quito 170525, Ecuador.
Energies 2021, 14(16), 4933; https://doi.org/10.3390/en14164933
Submission received: 3 June 2021 / Revised: 16 July 2021 / Accepted: 3 August 2021 / Published: 12 August 2021

Abstract

:
Achieving high penetration of electric vehicles (EVs) is one of the objectives proposed by the scientific community to mitigate the negative environmental impact caused by conventional mobility. The limited autonomy and the excessive time to recharge the battery discourage the final consumer from opting for new environmentally friendly mobility alternatives. Consequently, it is essential to provide the urban road network with charging infrastructure (CI) to ensure that the demand generated by EV users is met. The types of terminals to be considered in charging stations (CS) are fast and ultra-fast. The high energy requirements in these types of terminals could stress the electrical systems, reducing the quality of service. To size and forecast the resources needed in CI, it is of great interest to model and predict the maximum number of EVs, in each hour, that each CS will have to serve according to the geographic area in which they are located. Our proposal is not based on an assumed number of vehicles to be served by each CS, but rather it is based on the analysis of vehicular traffic in geo-referenced areas, so that the load managers can design the topology of the CS. The maximum vehicular concentration is determined by some considerations such as the road system, direction of the route, length of the road segment, the intersections, and the economic zone to which it belongs. The topology of the road map is freely extracted from OpenStreetMap to know the latitude and longitude coordinates. Vehicular traffic will be modeled through the topology obtained from OpenStreetMap and other microscopic variables to understand the traffic engineering constraints. In addition, the Hungarian algorithm is used as a minimum cost decision tool to allocate demand to the CS by observing vehicular traffic as a cost variable. The multi commodity flow problem (MCFP) algorithm aims to make commodities flow through the road network with the minimum cost without exceeding the capacities of the road sections. Therefore, it is proposed to solve the transportation problem by directing the vehicular flow to the CS while minimizing the travel time. This situation will contribute significantly to the design of the topology of each CS, which will be studied in future research.

1. Introduction

The use of conventional transportation is a current topic of study due to the reduction of air quality by the emission of C O 2 noise, and visual contamination [1,2]. It possible to solve this problem; electric mobility has attracted a great deal of interest from the scientific community. However, new engineering challenges arise, such as the sizing of CS to recharge the batteries of EVs that travel on a road network. The sizing problem consists of deciding on the number and type of chargers to install in the CS. The introduction of EVs on a large scale represents a heavy load on the electrical distribution systems (EDS) since they demand high electrical energy in short periods in systems with fast and ultra-fast charging terminals. Thus, the massive introduction of EVs to the system may stress EDS, causing voltage instability. It is possible to reduce the stress on the electricity system; it is necessary to predict the demand to estimate the necessary infrastructure resources for optimal sizing of the CS. The need to limit the level of EV penetration and the optimal distribution of available energy among them to make sure that the stress on distribution grids is reduced is noted in [3]. Thus, these constraints directly impact the user experience, so that the potential EV consumer would not be able to find any incentive to buy non-conventional mobility alternatives. Therefore, it is essential to create charging infrastructures that ensure the quality of electric service and maximize the user experience by recharging batteries to promote the consumption of electric vehicles.
The essential difficulty facing this new sustainable mobility alternative is the most significant travel distance, which is limited by the battery capacity and the lack of CI for EVs [4,5,6]. As EVs represent an essential load in EDS and match the residential load profiles, the need to size the CS arises considering the different topologies in terms of trajectory and how these are related to vehicular traffic variables. In order to model the effects of the massive introduction of EVs in the electric distribution networks, it is imperative to estimate the most considerable demand generated by electric mobility in geo-referenced areas. Thus, this study will make it possible to analyze the possibility of introducing new load peaks in substations to avoid stress on the power system. If the power demanded by electric mobility exceeds the installed power of the substations, it is necessary to expand the substation’s capacity and the transmission lines of the distribution power system. Therefore, a complex problem is faced due to EV user satisfaction, battery capacity, EC occupancy, charging station topology, medium and low voltage substation chargeability, and the amount of power available from the grid. On the other hand, considering a single electrical system to cover the demand generated by electric mobility would be very costly and technically unfeasible. Therefore, by analyzing the most significant traffic density in a specific geo-referenced area, it is possible to design the topology of the CS and the primary energy sources that will give the electrical power needed for each CS.
Thus, a heuristic model is proposed using Hungarian and MCFP algorithms based on linear problems. Previous work has solved [7] the optimal place of charging stations based on linear problems; for such a reason, the present article assumes that the locations of charging stations are the optimal locations and spend efforts in directing the vehicular flow to the different CS at the lowest cost [8,9,10]. The microscopic model simulates vehicular density in a geo-referenced area and understands vehicular flow’s basic variables. The essential traffic flow variables are flow (or volume), speed, and density. Flow is defined as the number of vehicles per unit of time, distance per unit of time is speed, and several vehicles per unit of distance are defined as vehicle density [11,12]. These variables will allow us to relate the different nodes by applying graph theory as an essential tool to formulate the different equations in an optimization problem [13,14,15]. An essential aspect being considered in the microscopic analysis is that they are helpful in situations of free traffic, i.e., the movement of the vehicle on the traffic road does not take into account the behavior of the driver, who decides the characteristics of such movement at any given moment. The analysis of vehicle traffic will allow us to look at the physical characteristics of the studied area, the characteristics of the vehicles, and variables such as vehicle spacing and the average speed at which an electric vehicle travels on a road map. Spacing refers to the spatial advance between vehicles. The speed and spacing between vehicles are modeled with a random vector under the discrete Poisson distribution [16,17,18].
Sizing and allocating resources to CS minimizes the capital cost and thus allows for testing the use of existing electrical systems or offering incentives for the construction of new distribution systems, in order to maximize the user experience with the construction of CIs with fast and ultra-fast charging terminals. The reviewed literature shows studies that address CS sizing under expected demand without discussing the optimization of EC sizing solutions. This work aims to develop a strategy to decide the most significant vehicle density and how this is directed towards the multiple CS options with the shortest travel paths. Furthermore, they decide the number and types of chargers based on an assumed number of vehicles that need to be recharged. Figure 1 shows the conceptualization of the problem to be solved. The contributions of this article are as follows.
  • A theoretical sizing model is introduced to develop charging station infrastructures for electric vehicles in a heterogeneous transportation system considering the paths; this application will allow us to foresee the highest and lowest demands in periods, thus allowing the sizing of the charging infrastructure with the minimal cost.
  • The algorithm is designed to reduce vehicle paths using a geo-referenced vehicle flow routing system considering the topology of the road map.
  • The organization and management of road vehicle traffic are solved with the Hungarian algorithm and multi-commodity flow problem under linear programming considerations.
  • The proposed academic model provides industrial application potential for charging station infrastructure providers in smart cities for competitive environments.
Finally, this article is organized as follows: Section 2 discusses the requirements for charging infrastructure. Section 3 discusses the relevance of geo-referenced path analysis and then details the problem formulation. Section 4 presents the results. Finally, Section 5 concludes this article.

2. EV Charging Station Infrastructure

The importance of building charging infrastructure for electric mobility, considering vehicular traffic in geo-referenced areas, depends on the sizing of the charging station capacity and how the traffic is distributed on the road map. The problem of sizing a recharging station usually consists of minimizing the cost of the capital investment in recharging systems, provided that CS can satisfy a fixed number of recharging requests [19]. A queuing theory-based CS sizing algorithm that benefits EV users, cost savings for CI, better station utilization, and reduced impact on power systems and the environment is proposed in [20]. Furthermore, it explains that grid energy loss and traffic congestion are vital factors that affect the location and sizing of CS. Thus, the primary considerations for the building of CIs in this article are vehicle density in geo-referenced scenarios and the routing of traffic to the various charging stations to optimally size and allocate resources to ensure that the demand generated by electrical mobility is satisfied. Therefore, three areas of research are related to the development of this article: demand assignment to charging stations, shortest trip trajectories, and sizing of operational capacity in EV charging stations. There are microscopic and macroscopic models that relate different traffic variables such as volume, speed, interval density, and spacing [21,22,23]. Consequently, to understand the different elements of vehicular traffic and how they relate to each other in a geo-referenced area, the microscopic model is used. The results generated from a microscopic analysis are valuable to characterize the most significant vehicle density in urban areas. As a result, transportation evaluation strategies are a severe problem to be solved because they involve geographic elements, vehicle distribution, and traffic variability subject to human behavior that make it impossible to predict with precision the path of a vehicle’s patterns. The trajectory and demand allocation have been considered decision variables for the least-cost routing of vehicular traffic and the optimal sizing of the charging station. In [24], the path of vehicles is considered based on realistic trip patterns; however, considering a rigid trip pattern reduces much of the size of the analysis space, i.e., a starting and ending point is considered, reducing the sample size where multiple users participating in the road network can recharge their batteries. Therefore, this article studies the most significant vehicular density that a zone could experience to give the necessary resources for CS distributed in the study and to guarantee the possibility of recharging users who use the road network.
The sizing of the charging stations relates to the technology of the terminals required to charge the battery of the electric vehicles on the road networks. Likewise, these technologies are approved and refer to the capacity of the batteries’ charging times: in fact, the range of travel of an electric vehicle on a road network with 100% of the battery charged. The average range of an EV for the class by European homologation criteria L7e-C is 123 km with an average battery capacity of 15 kW. In addition, the automotive mass must be less than or equal to 400 kg, excluding the mass of the batteries for electric vehicles. One more detail is that L7e-C corresponds to four-wheeled micro cars. Table 1 shows some characteristics of vehicles designed for urban use of type-approval class L7e-C.
Therefore, the decision problem is to direct the flow of vehicles to charge stations, reduce trip times, allocate users to each charging station, and generate efficient dispatch policies [25,26]. Moreover, the increased penetration of unconventional vehicles in electrical distribution systems needs particular attention because the use of high power in short periods could result in voltage instability, causing the reduction of the quality, safety, and continuity of the electrical service [27]. Therefore, the charging infrastructure for batteries in electric vehicles plays an essential role in motivating consumption and developing this new unconventional mobility alternative. The importance of sizing CS with optimizing criteria for electric mobility lies in the fact that the electricity generated cannot be stored, as is the case with conventional service stations, until users have access to the electricity available at the charging stations.
In [28,29,30], the authors propose modeling mixed vehicular traffic where conventional and unconventional mobility participates; furthermore, it is suggested that multiple charging terminals be considered according to homologation according to the requirements needed by each user. Several variables such as limited capacity, path chance, trip time, flow conservation, charging time, spatio-temporal collaboration, terminal types, and traffic events are considered in [31,32,33], and the main aspects to obtain robust solutions that respond to the different situations where unconventional mobility practices exist are introduced. There are three types of terminals that have a direct relationship with EV charging time. Type I and II terminals are applicable in residential areas for slow and semi-fast charging, respectively, while type III terminals, corresponding to ultra-fast charging, apply exclusively to charging infrastructures in the public sector. Each type of terminal to recharge batteries in electric vehicles has its topologies and specific technical features. For example, in ultra-fast charging terminals, the input voltage is DC, while in slow charging terminals, AC is used. Therefore, the optimal dimensions of the charging stations are critical so that they can meet the hourly demand for battery recharging. In [34,35], methods for clustering and the construction of minimum spanning trees using graph theory as optimization tools under a linear programming formulation are proposed.
Current contributions in the literature have been focused on multi-objective models to solve the charging station planning problem to incentivize the mass introduction of electric vehicles. However, microscopic analysis has not been considered to study the density of acceptable vehicular traffic in an area under study. Neither has it considered the multiple options that an electric vehicle can use to recharge its battery. These EV battery recharging options will be limited by distances and the availability of roads to design paths that involve reduced vehicle density, thus optimizing trip times to the charging stations. Another fundamental contribution is that by using free data downloaded from the OpenStreetMap platform, the case study can be recreated to add to the analysis of vehicular traffic using the geographic characteristics of the road network, topology, and lengths of the different road cross-sections.
Consequently, the interest of this article is to build a linear programming model that incorporates a study of free traffic situations that allows us to test the most significant vehicle density concerning paths in geo-referenced areas in order to direct the vehicular flow to the available charging stations. Predicting the most significant number of vehicles that could circulate in a geo-referenced area and understanding how they are distributed in trajectories, the designer of the charging infrastructure could use the proposed model to make the best decision about the topological design of the CS. Therefore, the study developed in this article has a high potential for the development of service providers for the creation of charging infrastructures, and the analysis is based on the study of geo-referenced traffic in terms of trajectories, looking at the features of the road network to decide the input variables for the linear optimization model. Finally, Table 2 shows the taxonomy of the related studies, where one can see a summary of the problem statement, the objective function, and the main constraints analyzed with which they discussed the problem.

3. Framework and Model

This article proposes a theoretical model to give resources and solve the transportation problem using linear programming in heterogeneous scenarios for CI implementation. The main aim is to allocate resources and guide the flow of demand to the charging station using the most extraordinary transit capacity to cut travel times from a node where demand originates to a supply station.

3.1. Path Analysis with Graph Theory Using Free Data

The study of traffic is highly complex since it involves people’s behavior and their particular mobilization needs. However, in this article, the purpose is to test the most critical scenario, which is when there is a high traffic density, for the optimal sizing of CS. Sizing for charging stations is directly proportional to the number of vehicles served by each charging point during peak demand periods; this data allows for building and designing with technical criteria the best topology required for each charging station distributed in the study area. The CS topology refers to the required number of fast chargers and ultra-fast chargers, and even battery replacement could be considered as an alternative to cut the times needed in a supply station for electrical mobility. Consequently, the purpose of this article is to build a linear programming model that includes a study of free traffic scenarios that allows us to test the most significant vehicle density in terms of trajectories in geo-referenced areas, in order to be able to direct the vehicular flow to the available charging stations in the case study.
Another essential detail in this article is that this model works with data available on a freeway as obtained from OpenStreetMap. In addition, OpenStreetMap will allow us to download essential data of latitude and longitude coordinates together with their mapping to use them in real places, allowing us to execute dimension processes in charging infrastructures from free data obtained from the web. Therefore, since the maximum number of vehicles that could enter a CS is predicted, IC designers would not have to start from an assumed number of vehicles, as has been done so far, to design the CS topology.
A non-convex combinatorial problem is used to solve the multi-commodity flow problem, including the Hungarian algorithm. The traffic rates and traffic directed to CS. Multiple commodities are generated with the optimal allocation. Thus, the following necessary constraints must be considered: (i) commodities, (ii) travel demand, and (iii) capacity of crossing between two intersections. The architecture of the proposed method is shown in Figure 2.

3.2. Target Location Problem for Multi-Commodity Flow

Programming in the analysis of geo-distributed data is proposed as a multi-commodity flow target location problem. The commodities can be sent from its resources; MCFP aims to locate its targets so that the multi-commodity flow is optimized in some direction.
Let G ( V , E ) be a directed graph where V represents the vertices and E the edges. The vertex represents the intersections where the demand for the commodity d and the CS are generated; moreover, V is the form of length n + m . The scalar n indicates the number of crosses and m the number of CS. The edges represent the different directions a path will take from a vertex u i to a vertex v j . The set of edges is e i , j or e j , i , where e i n E and i , j indicate the direction of flow that originates the optimal trajectories for each commodity set ( K ) composed of origin ( u k ) , destination ( v k ), and demand ( d k ). Therefore, V ( u , v ) , E ( e i , e j ) , and k = ( u k , v k , d k ) K . At each cross section of length ( l ) , there is a number of vehicles ( α ) that generates a traffic rate ( q ) . The vector of concentration of vehicles ( X ) will be given by the relation between q and the average speed V E L at which the vehicles circulate in a time interval t s . The length of the goods vector K will be directly proportional to the number of crossings it generates with x 0 , where x X and represents the partial concentration of vehicles at each node. The equations that define what has been discussed so far are presented below.
α i , j = l i , j + s p s p λ [ v e h ]
q i , j = α i , j t s [ v e h h ]
x = q i , j V E L [ v e h k m ]
With (1), calculate the largest number of vehicles in a cross-section of the road ( l ) , considering the spacing between vehicles ( s p ) . λ is the annual increase rate of the vehicle park. The rate or frequency of circulation of a vehicle in a given time interval is calculated with (2). Finally, (3) calculates the concentration of vehicles that will be received as a flow of goods over each node. A discreet Poisson probability distribution is used to model the stochastic of vehicle traffic. Therefore, the MCFP can be formulated as follows:
min X f T X
A X b
A e q X = b e q
b l X b u
The objective function presents with (4) and the constraint set with (5)–(7). f, b, and b are all constant vectors. A and A e q are constant matrices. The limits b l and b u can include constants or infinity. Finally, X is the vector of vehicle concentration variables. The MCFP variables are expressed as (8).
x k , e f l o w = Flow of Commodity k at Edge e , k K , e E
The bounds define the feasibility region for the MCFP variables, as expressed in (9) and (10).
x k , e f l o w 0 ; k K , e E
z 0
Let us say that percent injected flow (PIF) describes the relationship between the flow of a commodity injected into the network and the demand for the commodity as expressed in (11). The objective of the MCFP can be described as maximizing the lowest PIF of all commodities as expressed in (12).
z k = v V x k , ( s k , v ) f l o w d k
Objective : Maximize min { z k : k K }
The capacity constraint can be summarized as the sum of the flow on any edge; therefore, it must be less than or equal to the capacity of the edge as presented in (13). The net flow for commodity k of edge e is represented in (14).
c e k K x k , e f l o w ; e E
A e q ( k , u v ) = w V x k , ( u v , w ) f l o w w V x k , ( w , u v ) f l o w
b e q ( k , u v ) = 0 , k K , u v { w V : w u k , v k } b e q ( k , u v ) = d k , k K , u v { w V : w = u k } b e q ( k , u v ) = d k , k K , u v { w V : w = v k }
With (15), the conservation constraints are defined. The flow network describes a commodity through the sum of flows leaving and entering a node. The conservation constraint can be defined as a flow network that would be equal to zero if the node is not a source of a commodity; however, if the node is a source, the network flows would be equal to d k . The network flow objective must be equal to d k . Finally, Algorithm 1 is used to solve the problem defined in this article. In step 1, demand is assigned to each CE considering vehicular traffic as a cost variable. Furthermore, in step 2, the commodities are each represented as a 3-tuple containing an origin node, a destination node, and a demand. The objective function is defined in step 3. The constraint and conservation equations are defined with steps 4 and 5. Finally, in steps 6 and 7, the problem is solved, and the results are interpreted, respectively.
Algorithm 1 Optimal Sizing for Charging Stations Infrastructure (OSCSI).
Input: Case study
connectivity matrix = get ( . o s m )
G 0 = digraph(connectivity matrix)
Output: Global solution
G s o l v e ;
Variable: u k , v k , d k , k , q , w , z , x K ( u k , v k , d k ) f l o w , c o s t q , W = g e t (connectivity matrix);
Step 1: Allocation
[ u k , v k , x ] = Hungarian ( c o s t q , x ) ;
Step 2: Commodities
K = [ u k , v k , d k ] ; K ( f i n d ( K ( 3 ) = 0 ) ) = [ ] ;
Step 3: Objective function
z k = v V x k , ( u k , v ) f l o w d k
A ( k , e ) = z k ; A ( k , z + 1 ) = 1
f ( z + 1 ) = 1 ; Define the objective function
Step 4: Capacity constraints
A ( e , k ) f l o w = 1 ; b ( e , k ) f l o w = G 0 .Edges.Weight
Step 5: Conservation constraints
t e m p = w V x k , ( u v , w ) f l o w w V x k , ( w , u v ) f l o w
A e q ( k , u v ) = t e m p ; A e q ( u v , k ) = t e m p
b e q ( k , u v ) = 0 , k K , u v { w V : w u k , v k }
b e q ( k , u v ) = d k , k K , u v { w V : w = u k }
b e q ( k , u v ) = d k , k K , u v { w V : w = v k }
Step 6: Dual Simplex solver
[ x ] = l i n p r o g ( f, A, b, Aeq, beq)
Step 7: Return
[ G s o l v e ] = digraph( x , u k , v k )

4. Analysis of Results

The microscopic study presented with Figure 3 allows for understanding three fundamental variables of vehicular traffic, speed, flow (or volume), and vehicular concentration, considering the topology of the road map, the speed at which a vehicle travels, and frequency with which a vehicle passes through a point. Thus, the microscopic analysis will make it possible to describe vehicular traffic in terms of the road map in geo-referenced areas, bringing the problem closer to reality. Data obtained from the microscopic model concerning vehicular concentration are entered into the MCFP mathematical model proposed in this article.
The analysis in Figure 3 shows the predicted vehicle concentration over 20 years with a theoretical annual rate of increase of λ = 5.4 % . During the first 10 years, the growth of the vehicle park is modeled with a constant average speed of 15 km/h, from which about 4816 vehicles are obtained that will circulate at most during one hour in the most critical condition. The critical condition occurs when s p is minimum, given that a section of track could accommodate the largest capacity of vehicles given its length. A random vector with Poisson distribution will be used in this article, where 0.03 s p 2.5 km will be used. Figure 3 in year 11 includes microscopic analysis with three types of average traffic speeds 10 , 15 and 25 km/h to understand vehicle behavior as a function of traffic speed. Speed has an inverse effect on vehicular concentration, i.e., the less vehicular freedom of movement in a geo-referenced area indicates that there is a high vehicular concentration, and the greater the vehicular freedom of movement, the lower the vehicular concentration. This is because speed is inversely proportional to the concentration of vehicles in a time range. Thus, high vehicular concentration in a geo-referenced area is modeled when the travel speed is the least, and s p is the least. Therefore, the functionality of the microscopic model, implemented in the heuristic, is demonstrated for the analysis of vehicular traffic considering features in geo-referenced areas, allowing the problem to approach reality.
Figure 4 shows the first conditions where the topology of the road map is defined considering the high vehicle stochastic and the most significant number of vehicles under Poisson random distribution that can circulate in a cross-section of the road map. Furthermore, Figure 4 shows the CS labeled with nodes 55 to 61 in green color, and nodes 1–54 are the points where the demand originates. Remember that the aim is to allocate resources to charging stations for electric vehicles at the least cost to meet the demand generated as a function of the most significant vehicle density in a geo-referenced area. Thus, each intersection is allocated a demand based on the concentration of vehicles, which is analyzed with the microscopic model mentioned in previous paragraphs. The main aim is to obtain the smallest path for the supply of electric power from a demand node to a supply node by looking at the most significant and most negligible limits of vehicular flow present in a cross-section of the road map. The demand is generated by the electric vehicles that need to recharge their batteries, and the power supply to the batteries will be provided by the charging stations distributed in the study area.
As there are multiple demands and a restricted number of CS, the commodities originate from the charging stations to the demand nodes, and the path approached by the algorithm will be the optimal one, considered the existing flow of vehicles in each time interval analyzed. The objective function of the problem is directly related to the flow of vehicles X; this will allow for cutting, as much as possible, the trip time from a demand node to an energy supply center for battery recharging in electric mobility. The demand generated at each intersection is related to K and will be the variable that will allow sizing the CS based on the vehicular flow directed to each charging station.
The hourly distribution of vehicles that will circulate in the case study is presented in Figure 5. In Figure 5, there is a total of vehicles without classifying what type of vehicle each is according to its primary energy, that is, we only know that at the hour 07:00 H, about 384 vehicles are circulating, without identifying how many of them are combustion or electric, buses, or trucks. Thus, an essential piece of advice for CI sizing is to test the percentage penetration of electric vehicles in the different regions where a case study is conducted. In the present case study, we assume that 100% of the number of vehicles distributed on an hourly basis in Figure 5 are electric. This consideration will allow us to size the necessary resources to cover the demand for battery recharging considering the 100% accessibility of electric vehicles in the vehicle park that participates in the study area. For this case study, it has been assumed that only light vehicles move in this area, and we consider that all of them are electric vehicles of class L7e-C with an average power of 15 kW for the sizing of their engine and that they are exclusively for urban use with an average range of 100 km and a battery capacity of 15 kWh. Figure 5 reveals another vital detail: the empirical criteria used to assume the most significant number of vehicles that will use the road map, μ + 2 σ , where μ corresponds to the average value of the demands over 24 h, and σ corresponds to the standard deviation. This criterion can be helpful to foresee the most significant power of the electric substation that will give the electric resource to the charging infrastructure for electric vehicles ensuring that the demand for battery recharging is covered in high demand hours from 06:00 H to 20:00 H. Therefore, the least energy demand considered in the worst scenario is 6.4 MWh.
The vehicular traffic concentration and the different trajectories are shown in Figure 6. It has been considered that the hour of highest demand occurs at 07:00 H (see Figure 5) in executing the algorithm, thus resulting in assigning a CS to each vehicular flow from the demand node to the consumption node in periods of high vehicular traffic concentration. Linear programming is based on the observance of constraints that limit the problem, and globally optimal solutions are obtained. An interesting detail that can be verified with Figure 6 is the observance of the conservation constraints, where it is satisfied that the supply provided by the CS satisfies the demand, achieving equality between the source (charging station) and the demand (need for battery recharging).
Below, with Table 3, a summary from Figure 6 is shown for the reader’s convenience.
Table 3 shows the number of vehicles to be served by each charging station. The standard deviation of the number of vehicles to be covered by each CS at 07:00 H is five vehicles, which means that the variability of the number of vehicles between charging stations corresponds to only 10.38%. As a result, given that there is no significant variability, it can be concluded that the mathematical model proposed in this article not only assigns the resource with a lower cost but is also efficient in balancing the demand in the CI.
Therefore, in Figure 6, the planner can be supported to plan the sizing of each CS by analyzing the vehicular flow paths. Thus, the designer may be able to make engineering decisions considering real geo-located scenarios. On the other hand, the length of the accumulated roadway of each cross-section is another consideration: in this case, 13.98 km of the road network. This cumulative length parameter, when multiplied by the traffic concentration with units Veh km , will allow us to calculate the number of vehicles that a CS should be able to serve (see Table 3). As a result, a linear model can give globally optimal solutions by observing the problem’s constraints. Future work will study the management of recharging at CS to increase customer satisfaction and control the power delivered to the EV batteries to guarantee the protection of the electrical distribution networks.

5. Conclusions

This article presents theoretical models for sizing and allocating resources to CS infrastructures for EVs in a heterogeneous transport system, considering trajectories and restrictions of capacity in each section of the road. This application will allow us to predict the most significant and most minor demands in hourly periods, thus allowing us to size the charging infrastructure with the least cost. Considering actual vehicular traffic variables in geo-referenced scenarios has allowed us to bring the problem closer to a real scenario. Microscopic analysis has made it possible to consider traffic variables, obtaining results that show the most significant vehicle density generated according to the topology of the scenario. Thus, the results obtained guarantee reliability for decision making in the different stages of engineering for the sizing of CI. Furthermore, with the microscopic analysis, it has been possible to understand vehicular traffic behavior and how it affects the resource allocation processes in vehicle systems in smart cities. Moreover, this algorithm is designed to reduce vehicle trip length with a geo-referenced vehicular flow routing system; the topology of the road network and its physical characteristics obtained free of charge from the OpenStreetMap platform are considered. Therefore, a novel approach is proposed in this article, considering the spatial study of the area under analysis, i.e., how the vehicle’s speed, the physical dimensions of the vehicle, the lengths from intersection i to j, and the spacing between vehicles interact. On the other hand, the organization and management of vehicular traffic on the road to allocate resources to charging stations is solved using the Hungarian algorithm and the multi-commodity flow problem. Linear programming is used to solve the problem, where a scalable model is proposed to deal with the increase of electric vehicles and possible changes in hourly patterns of mobility. As a result, the proposed academic model has a potential industrial application for infrastructure providers of charging stations in smart cities in developing countries. Finally, the study of vehicular traffic should be extended to macroscopic models designed to describe vehicles’ general behavior when the flow has constant interruptions. A further detail that is not considered is the efficient use of CS to manage CS use.

Author Contributions

M.C.: Conceptualization, Methodology, Validation, Writing—review and editing. E.I.: Conceptualization, Methodology, Software, Writing—original draft. M.C.: Data curation, Formal analysis. E.I.: Supervision. J.C.: Writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Universidad Politécnica Salesiana and GIREI—Smart Grid Research Group under the project Charging of electric vehicles in buildings and its impact on the sizing and planning of electrical distribution networks.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This work was financially funded by the Smart Grid Research Group (GIREI) of Universidad Politécnica Salesiana (Project: Chargeability of the electrical distribution network considering the massive inclusion of electric vehicles) and the Power Grids and Smart Cities (RECI), Salesian Institutions of Higher Education (IUS).

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclatures

G ( V , E ) Directed graph
V E L Average speed
VSet of nodes
ESet of edges
nNumber of intersections
mNumber of charging stations
u i Origin node
v j Destination node
e i , j Edge origin-destination node
e j , i Edge destination-origin node
i , j Path direction
KCommodity set
u k Source node of k K
v k Target node of k K
d k Demand of k K
l i , j Length cross-sectional area
α i , j Number of vehicles
q i , j Traffic rate
f T Transposed vector of linear objective function f
t s Time range
s p Spacing between vehicles
λ The annual rate of increase in the vehicle park
AInequality matrix A
bInequality vector b
A e q Equality matrix A
b e q Equality vector b
b l Vector limit lower
b u Vector limit upper
zObjective variable
z k PIF of commodity k K
c e Edge capacity constraint
Δ k , u v Net flow for commodity k of edge e
P I F Percentage of injected flow
xScalar vehicle concentration
XVector of variables for the flow of commodities

References

  1. Tang, X.; Bi, S.; Zhang, Y.J.A. Distributed Routing and Charging Scheduling Optimization for Internet of Electric Vehicles. IEEE Internet Things J. 2018, 6, 136–148. [Google Scholar] [CrossRef]
  2. Yin, X.; Zhao, X. Planning of Electric Vehicle Charging Station Based on Real Time Traffic Flow. In Proceedings of the 2016 IEEE Vehicle Power and Propulsion Conference, VPPC 2016-Proceedings, Hangzhou, China, 17–20 October 2016; pp. 1–4. [Google Scholar] [CrossRef]
  3. Hussain, S.; Ahmed, M.A.; Lee, K.B.; Kim, Y.C. Fuzzy logic weight based charging scheme for optimal distribution of charging power among electric vehicles in a parking lot. Energies 2020, 13, 3119. [Google Scholar] [CrossRef]
  4. Hacardiaux, T.; Tancrez, J.S. Assessing the environmental benefits of horizontal cooperation using a location-inventory model. Cent. Eur. J. Oper. Res. 2019, 28, 1363–1387. [Google Scholar] [CrossRef]
  5. Ferro, G.; Laureri, F.; Minciardi, R.; Robba, M. An optimization model for electrical vehicles scheduling in a smart grid. Sustain. Energy Grids Netw. 2018, 14, 62–70. [Google Scholar] [CrossRef]
  6. Croce, A.I.; Musolino, G.; Rindone, C.; Vitetta, A. Sustainable mobility and energy resources: A quantitative assessment of transport services with electrical vehicles. Renew. Sustain. Energy Rev. 2019, 113, 109236. [Google Scholar] [CrossRef]
  7. Campaña, M.; Inga, E. Optimal Allocation of Public Charging Stations based on Traffic Density in Smart Cities. In Proceedings of the 2019 IEEE Colombian Conference on Applications in Computational Intelligence, ColCACI 2019-Proceedings, Barranquilla, Colombia, 5–7 June 2019. [Google Scholar] [CrossRef]
  8. Wang, X.; Yuen, C.; Hassan, N.U.; An, N.; Wu, W. Electric Vehicle Charging Station Placement for Urban Public Bus Systems. IEEE Trans. Intell. Transp. Syst. 2017, 18, 128–139. [Google Scholar] [CrossRef]
  9. Ferro, G.; Paolucci, M.; Robba, M. An Optimization Model for Electrical Vehicles Routing with time of use energy pricing And partial Recharging. IFAC-PapersOnLine 2018, 51, 212–217. [Google Scholar] [CrossRef]
  10. Arafat, Z.; Fehling, T.; Kleiss, G.; Nacke, B. Application of rechargeable batteries of electrical vehicles as time dependent storage resource for the public electricity grid. Energy Procedia 2018, 155, 478–491. [Google Scholar] [CrossRef]
  11. Hidas, P.; Wagner, P. Review of Data Collection Methods for Microscopic Traffic Simulation. In Proceedings of the World Conference on Transport Research WCTR, Istanbul, Turkey, 4–8 July 2004; Volume 2, pp. 1–19. [Google Scholar]
  12. Khondaker, B.; Kattan, L. Variable speed limit: A microscopic analysis in a connected vehicle environment. Transp. Res. Part C Emerg. Technol. 2015, 58, 146–159. [Google Scholar] [CrossRef] [Green Version]
  13. Akbari, M.; Brenna, M.; Longo, M. Optimal locating of electric vehicle charging stations by application of Genetic Algorithm. Sustainability 2018, 10, 1076. [Google Scholar] [CrossRef] [Green Version]
  14. Li, Q.; Xiao, Y.; Zhao, S.; Zhu, X.; Wang, Z.; Liu, Z.; Wang, L.; Yan, X.; Wang, Y. Performance status evaluation of an electric vehicle charging infrastructure based on the fuzzy comprehensive evaluation method. World Electr. Veh. J. 2019, 10, 35. [Google Scholar] [CrossRef] [Green Version]
  15. Zhang, J.; Zhou, H.; Li, H.; Liu, H.; Li, B.; Liu, C.; Yan, J. Multi-objective planning of charging stations considering vehicle arrival hot map. In Proceedings of the 2017 IEEE Conference on Energy Internet and Energy System Integration, EI2 2017-Proceedings, Beijing, China, 26–28 November 2017; pp. 1–6. [Google Scholar] [CrossRef]
  16. Demir, E.; Huckle, K.; Syntetos, A.; Lahy, A.; Wilson, M. Vehicle Routing Problem: Past and Future; Springer International Publishing: Berlin/Heidelberg, Germany, 2019; pp. 97–117. [Google Scholar] [CrossRef]
  17. Kijewska, K.; Iwan, S.; Małecki, K. Applying multi-criteria analysis of electrically powered vehicles implementation in urban freight transport. Procedia Comput. Sci. 2019, 159, 1558–1567. [Google Scholar] [CrossRef]
  18. Liu, H.; Zhang, B.; Cui, M.; Xiao, P.; Li, L.; Wu, J. The output planning model of the electric vehicle charging stations based on game theory. In Proceedings of the Asia-Pacific Power and Energy Engineering Conference, APPEEC, Xi’an, China, 25–28 October 2016; pp. 1704–1708. [Google Scholar] [CrossRef]
  19. Khaksari, A.; Tsaousoglou, G.; Makris, P.; Steriotis, K.; Efthymiopoulos, N.; Varvarigos, E. Sizing of electric vehicle charging stations with smart charging capabilities and quality of service requirements. Sustain. Cities Soc. 2021, 70, 102872. [Google Scholar] [CrossRef]
  20. Asna, M.; Shareef, H.; Achikkulath, P.; Mokhlis, H.; Errouissi, R.; Wahyudie, A. Analysis of an Optimal Planning Model for Electric Vehicle Fast-Charging Stations in Al Ain City, United Arab Emirates. IEEE Access 2021, 9, 73678–73694. [Google Scholar] [CrossRef]
  21. Zhang, Y.; Ioannou, P.A. Integrated control of highway traffic flow. J. Control Decis. 2018, 5, 19–41. [Google Scholar] [CrossRef]
  22. Jiang, J.; Yue, Y.; He, S. Multimodal-transport collaborative evacuation strategies for urban serious emergency incidents based on multi-sources spatiotemporal data. Leibniz Int. Proc. Inform. LIPIcs 2018, 114, 1–8. [Google Scholar] [CrossRef]
  23. Liu, G.; Zeng, C. Location planning of charging station for electric vehicle based on urban traffic flow. In Proceedings of the China International Conference on Electricity Distribution, CICED, Hangzhou, China, 17–20 October 2016; pp. 10–13. [Google Scholar] [CrossRef]
  24. Hussain, S.; Lee, K.B.; Ahmed, M.A.; Hayes, B.; Kim, Y.C. Two-stage fuzzy logic inference algorithm for maximizing the quality of performance under the operational constraints of power grid in electric vehicle parking lots. Energies 2020, 13, 4634. [Google Scholar] [CrossRef]
  25. Luo, C.; Huang, Y.F.; Gupta, V. Stochastic dynamic pricing for EV charging stations with renewable integration and energy storage. IEEE Trans. Smart Grid 2018, 9, 1494–1505. [Google Scholar] [CrossRef]
  26. Mao, T.; Zhang, X.; Zhou, B. Intelligent energy management algorithms for EV-charging scheduling with consideration of multiple EV charging modes. Energies 2019, 12, 265. [Google Scholar] [CrossRef] [Green Version]
  27. Luo, L.; Wu, Z.; Gu, W.; Huang, H.; Gao, S.; Han, J. Coordinated allocation of distributed generation resources and electric vehicle charging stations in distribution systems with vehicle-to-grid interaction. Energy 2020, 192, 116631. [Google Scholar] [CrossRef]
  28. Zhang, X.; Rey, D.; Waller, S.T.; Chen, N. Range-Constrained Traffic Assignment with Multi-Modal Recharge for Electric Vehicles. Netw. Spat. Econ. 2019, 19, 633–668. [Google Scholar] [CrossRef]
  29. Zhao, Y.; Huang, H.; Chen, X.; Zhang, B.; Zhang, Y.; Jin, Y.; Zhang, Q.; Cheng, L.; Chen, Y. Charging load allocation strategy of EV charging station considering charging mode. World Electr. Veh. J. 2019, 10, 47. [Google Scholar] [CrossRef] [Green Version]
  30. Zhang, Y.; Ma, Q.; Zhang, X.; Gao, M.; Yang, P.; He, H.; Hu, X.; Aliya, B. Integrated Route and Charging Planning for Electric Vehicles Considering Nonlinear Charging Functions. In Proceedings of the 2018 IEEE 22nd International Conference on Computer Supported Cooperative Work in Design, CSCWD 2018, Nanjing, China, 9–11 May 2018; pp. 207–212. [Google Scholar] [CrossRef]
  31. Abedinia, O.; Lu, M.; Bagheri, M. An Improved Multicriteria Optimization Method for Solving the Electric Vehicles Planning Issue in Smart Grids via Green Energy Sources. IEEE Access 2020, 8, 3465–3481. [Google Scholar] [CrossRef]
  32. Vaziri, S.; Etebari, F.; Vahdani, B. Development and optimization of a horizontal carrier collaboration vehicle routing model with multi-commodity request allocation. J. Clean. Prod. 2019, 224, 492–505. [Google Scholar] [CrossRef]
  33. Tang, D.; Wang, P. Nodal Impact Assessment and Alleviation of Moving Electric Vehicle Loads: From Traffic Flow to Power Flow. IEEE Trans. Power Syst. 2016, 31, 4231–4242. [Google Scholar] [CrossRef]
  34. Inga, E.; Cespedes, S.; Hincapie, R.; Cardenas, C.A. Scalable Route Map for Advanced Metering Infrastructure Based on Optimal Routing of Wireless Heterogeneous Networks. IEEE Wirel. Commun. 2017, 24, 26–33. [Google Scholar] [CrossRef]
  35. Campaña, M.; Inga, E.; Hincapié, R. Optimal placement of universal data aggregation points for smart electric metering based on hybrid wireless. In Proceedings of the CEUR Workshop Proceedings, Pucon, Chile, 19–20 October 2017; Volume 1950, pp. 6–9. [Google Scholar]
  36. Gan, X.; Zhang, H.; Hang, G.; Qin, Z.; Jin, H. Fast Charging Station Deployment Considering Elastic Demand. IEEE Trans. Transp. Electrif. 2020, 6, 158–169. [Google Scholar] [CrossRef]
  37. Esteban, I.; Roberto, H.; Sandra, C. Capacitated Multicommodity Flow Problem for Heterogeneous Smart Electricity Metering Communications Using Column Generation. Energies 2019, 13, 97. [Google Scholar] [CrossRef] [Green Version]
  38. Azadi Moghaddam Arani, A.; Jolai, F.; Nasiri, M.M. A multi-commodity network flow model for railway capacity optimization in case of line blockage. Int. J. Rail Transp. 2019, 7, 297–320. [Google Scholar] [CrossRef]
  39. Bevrani, B.; Burdett, R.L.; Bhaskar, A.; Yarlagadda, P.K. A Multi Commodity Flow Model Incorporating Flow Reduction Functions; Springer: New York, NY, USA, 2020. [Google Scholar] [CrossRef]
  40. Liu, G.; Kang, L.; Luan, Z.; Qiu, J.; Zheng, F. Charging station and power network planning for integrated electric vehicles (EVs). Energies 2019, 12, 2595. [Google Scholar] [CrossRef] [Green Version]
  41. Ghasemi, P.; Khalili-Damghani, K.; Hafezolkotob, A.; Raissi, S. Uncertain multi-objective multi-commodity multi-period multi-vehicle location-allocation model for earthquake evacuation planning. Appl. Math. Comput. 2019, 350, 105–132. [Google Scholar] [CrossRef]
Figure 1. Applied network theory to solve the transport problem in MCFP.
Figure 1. Applied network theory to solve the transport problem in MCFP.
Energies 14 04933 g001
Figure 2. Architecture proposed OSCSI.
Figure 2. Architecture proposed OSCSI.
Energies 14 04933 g002
Figure 3. Heterogeneous time zone vehicle traffic average of 5 working days.
Figure 3. Heterogeneous time zone vehicle traffic average of 5 working days.
Energies 14 04933 g003
Figure 4. A study case at the level of trajectories in an area of 523.4 m 2 in Quito, Ecuador.
Figure 4. A study case at the level of trajectories in an area of 523.4 m 2 in Quito, Ecuador.
Energies 14 04933 g004
Figure 5. Heterogeneous hourly vehicular demand with V E L = 51.1 km h , with total transverse length of 13.98 km.
Figure 5. Heterogeneous hourly vehicular demand with V E L = 51.1 km h , with total transverse length of 13.98 km.
Energies 14 04933 g005
Figure 6. Optimal overall solution for battery charging in electric mobility.
Figure 6. Optimal overall solution for battery charging in electric mobility.
Energies 14 04933 g006
Table 1. European homologation class L7e-C features.
Table 1. European homologation class L7e-C features.
ElectricalAutonomyTerminalMotorMaximumRegenerativeCategory Average
Vehicle(km)Voltage (V)(kW)Speed (km/h)SystemDimensions (m)
Tazzari Zero20023015–2590
Renault Twizy1002301780Height < 2.5
Audi Urban73230–40015100-Width < 1.5
Peugeot BB11202301560-Length < 3.7
Table 2. Summary: comparison of the results with other methods.
Table 2. Summary: comparison of the results with other methods.
Author, YearAllocation SizingDensity TrafficUrban StudyHeterogeneous VehiclesScalableMicroscopic Traffic SimulationOther Characteristics
Gan et al., 2020 [36]Max profit service, nonlinear integer problem (NLIP)
Inga et al., 2019 [37]Min cost flow Multicommodity network
Demir et al., 2019 [16]Min distance travel Deterministic, Integer Linear Programming (ILP)
Azadi et al., 2019 [38]Routes, Mixed-Integer Linear Programming (MILP) Multicommodity network
Bevrani et al., 2019 [39]Nonlinear programming (NLP) Multicommodity network
Campaña et al., 2019 [7]Min congestion road Segmentation ILP Neighborhood
Liu et al., 2019 [40]Max flow Fuzzy multi-objective
Ghasemi et al., 2019 [41]Min shortages Multi-objective particle swarm optimization (PSO)
Proposed work (OSCSI)Multicommodity Flow, Revised simplex, MILP, NP-Hard, Hungarian method
Table 3. Partial demand for each charging station at hour 07:00 H.
Table 3. Partial demand for each charging station at hour 07:00 H.
IDVehicular# Vehicles
Concentration
55−2.8740
56−2.9541
57−3.6750
58−3.3847
59−3.7252
60−3.2646
61−3.7152
Total23.56328
Notes. # equals number.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Campaña, M.; Inga, E.; Cárdenas, J. Optimal Sizing of Electric Vehicle Charging Stations Considering Urban Traffic Flow for Smart Cities. Energies 2021, 14, 4933. https://doi.org/10.3390/en14164933

AMA Style

Campaña M, Inga E, Cárdenas J. Optimal Sizing of Electric Vehicle Charging Stations Considering Urban Traffic Flow for Smart Cities. Energies. 2021; 14(16):4933. https://doi.org/10.3390/en14164933

Chicago/Turabian Style

Campaña, Miguel, Esteban Inga, and Jorge Cárdenas. 2021. "Optimal Sizing of Electric Vehicle Charging Stations Considering Urban Traffic Flow for Smart Cities" Energies 14, no. 16: 4933. https://doi.org/10.3390/en14164933

APA Style

Campaña, M., Inga, E., & Cárdenas, J. (2021). Optimal Sizing of Electric Vehicle Charging Stations Considering Urban Traffic Flow for Smart Cities. Energies, 14(16), 4933. https://doi.org/10.3390/en14164933

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop