Energy efficient route planning for electric vehicles with special consideration of the topography and battery lifetime

In contrast to conventional routing systems, which determine the shortest distance or the fastest path to a destination, this work designs a route planning specifically for electric vehicles by finding an energy-optimal solution while simultaneously considering stress on the battery. After finding a physical model of the energy consumption of the electric vehicle including heating, air conditioning, and other additional loads, the street network is modeled as a network with nodes and weighted edges in order to apply a shortest path algorithm that finds the route with the smallest edge costs. A variation of the Bellman-Ford algorithm, the Yen algorithm, is modified such that battery constraints can be included. Thus, the modified Yen algorithm helps solving a multi-objective optimization problem with three optimization variables representing the energy consumption with (vehicle reaching the destination with the highest state of charge possible), the journey time, and the cyclic lifetime of the battery (minimizing the number of charging/discharging cycles by minimizing the amount of energy consumed or regenerated). For the optimization problem, weights are assigned to each variable in order to put emphasis on one or the other. The route planning system is tested for a suburban area in Austria and for the city of San Francisco, CA. Topography has a strong influence on energy consumption and battery operation and therefore the choice of route. The algorithm finds different results considering different preferences, putting weights on the decision variable of the multi-objective optimization. Also, the tests are conducted for different outside temperatures and weather conditions, as well as for different vehicle types.


Introduction
Combustion engine driven cars have been dominating our world for more than a century. With global warming concerns (Gota et al. 2019) and increasingly stringent environmental policies ahead, vehicles with electric motors could be part of the solution, a precondition being that electrical power comes from renewable sources. In addition, further challenges come up with this new technology option in transportation. Despite higher efficiency than combustion engines, the task to store energy in the vehicle is still challenging. The battery offers a much more limited range compared with conventional cars. Moreover, recharging an electric vehicle is more time-consuming than refilling a tank. This leads to a more complex trip planning with an electric vehicle. A good advice is to identify the location of charging stations beforehand and to plan the route accordingly. Energy consumption can vary significantly depending on the path chosen. This work elaborates on optimal route planning to a desired destination while considering the special characteristics of electric vehicles. The main focus lies on improving the battery lifetime, as well as minimizing energy consumption and journey time while taking into account impacts of topography.
The first objective of this work is to find a route from a starting to a destination point with the least amount of energy used. This task is expanded such that it is possible to calculate the shortest journey time as well as the best route to maximize battery lifetime. Those route planning options are combined into a weighted multi-objective optimization problem. Shortest path algorithms are used in networks consisting of nodes, edges, and edge costs to solve the optimization problem. With the help of these algorithms, it is possible to find the path from one node to another with the smallest edge costs. Based on a model of the electric vehicle, which describes the physics of different driving modes, the energy consumption and travel time are approximated. It includes the energy required for driving as well as additional loads for air conditioning and other accessories. The model also considers engine operation as either a motor or as a generator and thus respects regenerative energy.
The following "State of the art" section gives an overview of relevant work covering this topic. The "Methodology" section explains the model, which is used to calculate the vehicle's energy consumption, and the impact additional loads and efficiencies have on the vehicle. Furthermore, the street network model, the shortest path algorithm, and the multiobjective optimization are explained. The "Results" section presents the results of the optimization algorithm tested for an urban as well as a suburban area. In addition, a sensitivity analysis investigates the effects of temperature and different vehicle models from different manufacturers. The "Conclusion and outlook" section presents conclusions and elaborates on still open questions in this field of research.

State of the art
Creating an optimal route planning system for electric vehicles is multi-disciplinary and requires profound knowledge of electric vehicles, batteries, route planning algorithms, and dynamic optimization.
Modeling electric vehicles Energy efficiency and battery conservation are the main goals of the proposed optimal route planning system. Therefore, the electric vehicle is described based on a physical model in order to enable calculation of the energy consumption. In this context, one special characteristic of electric vehicles is the ability to regenerate energy from braking or driving downhill. In Lv et al. (2015) the energy flow from regenerative braking is modeled in detail. Another detailed model of the vehicle's energy consumption is presented in Yi and Bauer (2017). The goal is to find a high-resolution powertrain efficiency estimation. The energy consumption is described by the physical forces acting on the moving vehicle. A similar concept can be found in Maia et al. (2011), where an energy consumption simulator is presented. An alternative approach to estimate the energy consumption of an electric vehicle is found in Hayes et al. (2011), where a range estimator is created by using different drive cycles, which are originally based on combustion engines and then adapted for electric vehicles. However, the approach in Hayes et al. (2011) avoided modeling the physics of an electric vehicle.
Battery lifetime in electric vehicles Having batteries as an energy storage in vehicles comes along with many new challenges (see, e.g., Rizoug et al. (2018)) compared with cars with conventional combustion engines. Apart from their costly production, batteries have a limited lifetime. According to Li et al. (2017), by 2020 nearly one million electric vehicle battery packs will be "retired" in the USA alone, meaning that they cannot be used in electric vehicles anymore. This happens when (i) the battery only provides less than 80% of the original capacity and maximum power, or (ii) there are functional failures occurring. Several models to describe the battery degradation can be found in Pelletier et al. (2017), where cycle aging as well as calendar aging mechanisms are considered. Fernández et al. (2013) propose an optimal charging energy management that minimizes the degradation of the battery as a function of temperature and depth of discharge. In Wang et al. (2016) a battery degradation model of the cycle losses is used, which is linear dependent in cumulative current throughput, but quadratic in temperature. One interesting result is found in Peterson et al. (2010), where capacity fade of commercial Li-ion batteries caused by driving and vehicle-to-grid usage is studied. Battery degradation was found to be related to the amount of energy processed in the battery, while the depth of discharge has less effects on the battery lifetime. The importance of methods to increase the battery lifetime corresponds with the work in Hawkins et al. (2012), where it is shown that replacing a battery frequently may harm the environment and increases the risk of toxicity (e.g., for humans and freshwater) and metal depletion impacts.
Route planning algorithms Route planning problems are typically solved by applying shortest path algorithms. The Dijkstra algorithm (Dijkstra 1959) is a very efficient method to work with weighted graphs and is widely used in network theory. A Dijkstrabased speed-up technique for constrained shortest path problems is applied in Storandt (2012), where route planning is conducted for bicycles. The Dijkstra algorithm is also applied in Song et al. (2015). In that work, route planning in urban areas with influences of traffic is considered. The Bellman-Ford (Bellman 1958) algorithm has a higher complexity than Dijkstra, but can be used in networks with negative edge costs, which is useful for route planning for electric vehicles. Despite its high complexity, the Bellman-Ford algorithm works well with street networks. The Yen algorithm (Yen 1970) is based on Bellman-Ford, but it is generally able to solve the shortest path problem faster than Bellman (1958).
Route planning for electric vehicles It is mentioned in Neaimeh et al. (2013) that many people, who tested electric vehicles, experience so-called range anxiety. It turns out that many drivers would change their driving behavior and especially their route-choices to the destination. In order to find an optimal route that can extend the range of the vehicle, Dijkstra's shortest path algorithm is applied. The goal is not to simply find the energy-minimal route, but also to help with range anxiety, making drivers feel more comfortable with e-mobility. In Nunzio and Thibault (2017) a range estimation for online use is created, which is based on calculating the energy-optimal route. The vehicle's energy consumption is modeled including influences of traffic and with a shortest path algorithm (Bellman-Ford) the range can be estimated. In Storandt and Funke (2012) battery switch station are included in the network and a modification of Dijkstra's algorithm is used. This modification is done by using Johnson's shifting technique (see Johnson 1977) in order to include regenerated energy of the electric vehicle. In practice, battery switch stations are not established, but charging stations become increasingly visible. In Storandt et al. (2013), the approach is similar to Storandt and Funke (2012), but including charging stations instead of battery switch stations. It is also proposed to work on a multi-criteria optimization that includes the journey time and a maximum number of recharging events.
The present work implements a multi-objective optimization approach determining a path that seeks to meet best the drivers' requirements. In detail, the optimization algorithm uses weighting factors on these three optimization variables: The goal of this paper is to include the drivers' preferences in a very flexible way. For example, the drivers can choose a single-objective optimization of either of the above-mentioned variables, or give equal weights to all three. Charging stations for electric vehicles are not included in this work.

Optimization problem
The objective of this work is to design a route planning system for electric vehicles in order to optimize the battery lifetime, energy consumption, and journey time. A weighted multi-objective optimization approach is used to prioritize one aspect or the other.

Flowchart
To achieve the objective of this work-finding the weighted optimal path to a desired destination with an electric vehicle-it is important to have a model describing the energy consumption of the electric vehicle and the street network. Then, shortest path algorithms can be used for optimization. A flow chart (see Fig. 1) presents an overview of this optimization problem. At first, a starting and a destination point have to be defined. Then, it is necessary to have a road network that includes topography information.
Each road section has a defined length s, a velocity v, and a slope q in %, calculated from the topography information. The next step is to define the vehicle parameters. The initial state of charge of the battery SoC init and its maximum capacity SoC max are necessary, as well as the mass m, the drag coefficient c w , and the cross-sectional area A of the vehicle. When having information on the outside temperature T , all parameters are available for the calculation of the energy consumption model of the electric vehicle and the journey time on all the possible roads of the network. As this work formulates an optimization problem with constraints, the boundary conditions have to be defined. Naturally, the state of charge of the battery SoC cannot be negative or above the maximum capacity. In the interest of the vehicle owner, the SoC should not fall below a certain minimum capacity, because deep discharge can decrease battery lifetime. A factor a, 0 ≤ a < 1, is added to the optimization problem, such that the constraint equation becomes The task will be extended to a multi-objective optimization with the variables energy, time, and battery lifetime. The factors γ and δ, with γ + δ ≤ 1, are used to give weights to the optimization variables. The last step is to apply a shortest path algorithm.

Energy consumption for driving
We start with the calculation of the energy consumption of an electric vehicle for driving. The basic principles of the calculations are similar to conventional vehicles. The driving force F drive is calculated as the sum of the rolling resistance F roll , the air resistance F air , and the gradient resistance F grad , while other factors are neglected for simplicity: The next step is the calculation of the power resulting from the driving force and then finally for the energy required over a certain period of time t 1 ≤ t ≤ t 2 .

Total energy consumption of the electric vehicle
Apart from regenerative energy from braking or driving downhill, the battery is the only source of energy in electric vehicles. Other loads than driving must Other so-called accessory loads are headlights, fan, windshield wipers, rear window heating, and radio. These need low voltage (12 V) and therefore a DC-DC converter. Figures 2 and 3 show the energy flow from the battery to the components of the vehicle with the engine operating as a motor or as a generator, respectively. E out is the energy the battery provides. E acc and E hc represent the low voltage accessories, and the heating and cooling, respectively. E drive is used for driving, as in (5). The amount of energy E that decreases the state of charge of the battery considering efficiency factors is if the engine operates as a motor (Fig. 2). With electric vehicles, recharging the battery with regenerated energy from braking or driving downhill is possible. In this case, the driving power P drive in (4) is negative and the engine operates as a generator (see Fig. 3). There are two cases to distinguish: (i) the battery charges, only if the regenerated energy is sufficient to cover the accessory loads; (ii) the battery discharges, if the regenerated energy cannot cover all the accessory loads. The energy E out = −E in can be calculated as and the total energy including charging and discharging efficiencies is If the efficiency of the engine operating as a generator is the same as for a motor (η m = η g ), Eqs. (6) and (7) become If the charging and discharging efficiencies of the battery are equal (η cha = η dis ), Eqs. (6) and (8) merge to one new equation with E out as in Eq. (9).
Modeling the street network The next task of the route planning is to model the real-life street network with its properties such as slope, distances, speed limits, and road type. In order to solve a shortest path problem, a specific structure of the street network is needed.

Networks
Networks consist of so-called nodes and edges. The nodes are connected with each other by edges, which can have assigned values called "edge weights" or "edge costs." Nodes Nodes are the decision points of the network. Each node is connected to different nodes via edges. Some nodes are directly connected, others are only reachable by passing other nodes on the way. In the real-world street network, nodes are intersections of roads.
Edges Edges are the roads connecting the nodes. In this network model, the edges are split into straight sections with constant slope and constant speed limit.
Edge costs Networks have the purpose to express relations between the nodes. There are different ways to characterize such relations. The most elementary ones are a binary relations, describing whether there is a direct connection (an edge) or not. Edge costs can have physical meanings, such as the distance between nodes. A shortest path algorithm finds the path that cumulates the least amount of edge costs on its way. When the physical meaning of the edge costs is the energy consumption of an electric vehicle, the network becomes bidirectional. The main reason is the topography, as two nodes can be on different altitudes and, therefore, the required energy is different for both directions.

Topography data
Topography data is obtained from the "USGS Earth Explorer" website (https://earthexplorer.usgs.gov/), where geographical data is available for download. The file obtained from the USGS is a TIFF-file with 3601 × 3601 pixel, containing a digital elevation map from the NASA Shuttle Radar Topography Mission (SRTM). It has a high spatial resolution of one arc-second for longitude (east-west) as well as for latitude (north-south). The height resolution is one meter.
Since two adjacent pixels are one arc-second apart, the TIFF-file covers an area of one degree in longitude and one degree in latitude. For the calculations of the energy consumption of the electric vehicle, it is necessary to convert the distances given in degrees into distances in meters. The conversion is achieved by an approximation of the earth as a perfect sphere 1 with radius of R earth = 6371 km = 6371000 m. φ 1 and ψ 1 are the coordinates of a point in degree longitude and degree latitude, respectively. φ 2 and ψ 2 are the coordinates of another point, and φ = φ 1 − φ 2 and ψ = ψ 1 − ψ 2 are the distances between point 1 and point 2 (in degrees). The distances x for longitude and y for latitude in meters are converted using and x = 2πR earth 360 cos(ψ 1 π/180) φ.
It can be noticed that Eq. (12) depends on the coordinates of the latitude. At the equator, we have cos(ψ 1 π/180) = 1. The coordinate lines of the longitude converge when moving further away from the equator. Two places with a φ of 1 • at the equatorial line have a higher x than two places with the same φ at another degree of latitude. The resolution of one arc-second of the map obtained from NASA means having a resolution of about 31 m north-south. It has to be mentioned that Eqs. (11) and (12) are only approximations of the real distances. The next task is to calculate the slope of a road section that is straight and has a constant gradient. The difference in altitude z between the start and the end point of the section, which is derived from the TIFF-file mentioned in the beginning of this section, is divided by the Euclidean distance between those points and multiplied by 100 in order to get the slope q in %:

Street network
After obtaining topographical data, the next step is to get an adequate network of roads containing the region of interest. One option is using "OpenStreetMaps," a big database for streets, roads, and paths. The use of "OpenStreetMaps" is problematic because many roads are missing labels with no possibility of distinguishing if they are footpaths, waterways, or streets. The same goes for data obtained from official government websites. Thus, a more simple approach is chosen. First, it is decided which roads should be part of the network.
Only the main roads in the area of interest are included by taking coordinates of road sections. Some of the roads (edges) are a few kilometers long and split up into straight consecutive sections with a length varying from 50 m to a few hundred meters. For the purpose of testing the algorithm, this resolution is satisfying.
For an advanced route planning system, a higher resolution would be needed. The coordinates of the start and end points of the road segments and the matching topography information obtained, the distances and height differences are calculated as explained in the "Topography data" section. This ensures each segment is straight and has a constant slope.
Optimization with a shortest path algorithm The multi-objective optimization is solved by applying a shortest path algorithm. Shortest path algorithms use networks made of nodes and edges. The goal of shortest path algorithms is to find the path with the least amount of edge costs. In this case the edge costs represent the energy consumption, journey time, or the battery degradation. The "Yen algorithm" section explains the shortest path algorithm for one optimization variable. Then, the multi-objective optimization based on the shortest path algorithm is described in the "Multi-objective optimization" section.

Yen algorithm
A well-known shortest path algorithm is the Bellman-Ford algorithm (see Bellman 1958). An improved version of Bellman-Ford is the algorithm proposed by Jin Y. Yen (see Yen 1970). Both algorithms find the shortest paths from all nodes of the network to one specific destination node. The iterative search goes backwards, always starting from the destination node. In contrast, the state of charge of the battery (SoC) is calculated forward, starting from the starting point. The constraints are violated if the SoC exceeds the boundaries (see Eq. (1)). In addition, the goal for energy-optimal route planning is to arrive at the destination with the highest state of charge possible. The backwards iterations of the Yen algorithm are, therefore, not suitable for the optimization problem with SoC constraints. Hence, the following solution is introduced. Instead of computing all minimum paths to the destination, all minimum paths from the start are calculated. The principles of the Yen algorithm remain the same, but there are a few modifications. The algorithm is applied to find the energy-optimal path in a network with N nodes in the following manner, starting with iteration 0, where d 1i are the edge costs from node 1 to node i = 2, . . . , N, representing the energy consumption required to travel from node 1 to i. The edge costs are calculated adding up the energy consumption E l (derived from Eq. (10)) of M adjacent road sections: f (k) i are the costs from node 1, the start node, to node i at iteration k. A distinction is made between odd and even numbered iterations. For odd iterations, the minimum is found using with i = 2, 3, . . . , N. For even iterations, the minimization is with i = N − 1, N − 1, . . . , 2. For all iterations, the battery constraints are checked. If SoC > SoC max , the path is still feasible, but the current state of charge is set to the maximum capacity of the battery (SoC = SoC max ), because no further charging is possible. The energy consumption f (k) j +d ji is adapted accordingly. If SoC < a SoC max , that path is considered unfeasible and f (k) i is set to infinity, such that this path is neglected in all further calculations. If all minimum paths f i of iteration k remain the same compared to the previous iteration k − 1, the optimum is found: For time-optimal route planning, the optimum g (k) i of each iteration k is calculated having the edge costs t ji (see Eq. (A1.1))), the journey time from node j to node i (details in Appendix 1). Similarly, the optimum for the battery lifetime h (k) i works with the edge costs a ji (see Appendix 2). a ji represents the absolute energy flow to and from the battery. The main difference between d ij and a ij concerns regenerative energy, which decreases the value f (k) i , but increases h (k) i , and therefore the Yen algorithm finds different minima.

Multi-objective optimization
The multi-objective optimization has three optimization variables: energy consumption of the vehicle, journey time, and cyclic lifetime of the battery. The algorithm applied for this task is based on the modified version of the Yen algorithm from the "Yen algorithm" section. The optimal path from the start node number 1 to all nodes i is calculated in each iteration, taking into account all three variables and their weights, as well as the battery constraints. We start with iteration 0, setting all three variables in the same manner: The next step is to go into the details of the odd iterations. At first, are set. Then, for the other nodes i, the energy-optimal path is calculated following Eq. (16). If at least one path is feasible, E min = |f (2k−1) i | is set and T min , the minimum in time, and A min , the smallest absolute energy consumption, are calculated according to Eqs. (A1.4) and (A2.3), respectively. All three results are combined to multi-objective optimization with the weights γ for energy, δ for time, and 1 − (γ + δ) for battery lifetime. The optimization problem solved for node i = 2, 3, . . . , N in iteration 2k − 1 is: All energy costs are normalized by the factor E min , the journey time by T min , and the absolute energy by A min . The normalized values equal to one at the optimum and greater than one for all the other paths. The deviation from the optimum is weighed with γ , δ and 1 − (γ + δ). After the minimum x are set according to the resulting multi-objective optimal path. For even iterations the situation is very similar. We start with Then, for node i the optimal values E min , T min , and A min are computed and multi-objective optimization is performed, iterating from i = N − 1 to i = 2:

Reference values and assumptions
The proposed multi-objective optimization algorithm is implemented in MATLAB 2 (MATLAB 2019) and two vehicles selected for testing are the Nissan Leaf, see data sheet (https://www.nissan.co.uk), and the Mitsubishi i-MiEV, see data sheet (https://www. mitsubishi-motors.com/en/products/#search-models). All parameters necessary for the calculations are summarized in Table 1. For the mass m, the curb weight of the vehicle, is combined with the weight of one person, the driver. The driver's mass is assumed to be 75 kg. To calculate the air resistance F air according to Eq. (3), the cross-sectional area A and the drag coefficient c w are necessary. If A is not specified in the official data sheet, it may be estimated according to Haken (2013), knowing the width w and height h of the vehicle: Accessory loads and efficiencies are obtained from the measurements in Geringer and Tober (2012). Some efficiency factors cannot be found in Geringer and Tober (2012), therefore they are assumed to be 100%, except for the engine efficiency, which is assumed to be 90%. The main differences between the two vehicle types are mass and the dimensions of height, weight, and drag coefficient.

Suburban area: near Vienna, Austria
The street network covers the main roads of the Wienerwald area West of Vienna. The city of Vienna is excluded. The network has a total of 31 nodes and 41 edges. Different routes were tested with the multiobjective optimization algorithm. Planning a trip from Passauerhof (a in Fig. 4) to Maria Gugging (b) provides interesting results (see Fig. 4). There are two different solutions of the multi-objective optimization. The routes are: (a) The first path, where the electric vehicle needs the least amount of energy, starts at Passauerhof, and passes through Katzlsdorf, Königstetten, and St. Andrä before it reaches the destination Maria Gugging (see Fig. 4). (b) The other route (b) is both the fastest and the shortest route, but presents more variation in topography than (a), passing small villages such as Unterkirchen and Hintersdorf before reaching Maria Gugging (see Fig. 4). Table 2 presents the solutions of the multi-objective optimization for various combinations of the weights.
It can be noticed that combinations with small values of γ and δ result in option (a). This route has less variation in topography than (b), as it can be seen in Figs. 5 and 6. If γ and δ are small, the focus is on increasing the battery lifetime. In contrast to the optimization variable for energy consumption, where the regenerated energy has a negative sign and, therefore, helps minimizing the costs, the optimization variable representing the battery lifetime takes the absolute value of the energy. In this case, the regenerated energy counts the same way as the consumed energy. Therefore, driving downhill or braking increases the costs and is avoided by the optimization algorithm. A route with little elevation up and down is preferred. It can be seen in Fig. 6 that the absolute energy flow to and from the battery for route (a) is smaller than for route (b). Combinations with a high value of δ lead to the fastest route (b). The difference in time is almost four minutes, which is a quarter of the total journey time of route (b). The difference in energy consumption is less significant (see Fig. 5). Choosing route (a), only

Comparing under different temperatures
The route planning system is also tested for different outside temperatures. During winter, the temperatures frequently fall below zero in Austria. Therefore, heating is necessary for the passengers in the vehicle. With the help of the measurements in Geringer and Tober (2012), it is possible to find a linear function of the Nissan Leaf's power demand for heating, which is P heat = 90 W/ • C. The energy consumption is expected to rise with lower temperatures, which may influence the results of the multi-objective optimiza-tion problem. Assuming an outside temperature of T = −10 • C, we have with T 0 = 20 • C. For the relatively short trip of around 20 min, 0.6 kWh are used only for heating. In addition, battery (dis-)charging efficiency η dis is decreasing due to low temperatures (see Table 1). Table 3 shows the results for varying combinations of the optimization weights. These results are not surprising considering that the energy consumption of all accessories, including heating and cooling, is time-dependent. This factor makes fast routes more attractive also when looking at energy-optimal route planning. In this case, route (a), which has less energy consumption at 20 • C, now requires more energy than route (b), making route (b) the most energy and time efficient route at temperatures of − 10 • C outside (see Fig. 7). When focusing on the cyclic lifetime of the battery only (low γ and δ), route (a) is still chosen because of the topography, as explained previously. Figure 8 shows that the cumulated absolute energy consumption and regeneration of route (a) is still smaller than that of route (b), but with a narrower margin at − 10 • C compared with 20 • C.

Comparing two different types of electric vehicles
So far, the tests are performed exclusively with the Nissan Leaf. In the following part, the Mitsubishi i-MiEV vehicle is added, assuming its characteristics as given in Table 1. The Mitsubishi vehicle is almost 30% lighter in weight than the Nissan, but its c w -value is higher and its efficiencies are lower. In Table 4, the results for the Nissan Leaf and Mitsubishi i-MiEV are compared. Figure 9 shows that the energy consumption on both routes, (a) and (b), is higher for the Mitsubishi despite having less weight. Another finding is that on both routes the cumulated absolute value (Fig. 10) of the consumed and regenerated energy of the Mitsubishi is smaller compared to Nissan. This result shows that the Mitsubishi regenerates less energy. Only once does the result change compared to the Nissan, which is for γ = 0.2 and δ = 0.4. The absolute energy has the weight of 1 − (γ + δ) = 0.4. In this case, the objective to be minimized in (27) is almost equal for both routes, but with route (b) being slightly more efficient.
Urban area: San Francisco, CA Further tests are conducted for the city of San Francisco, CA, because of its interesting topography. The street network of San Francisco is fundamentally different from the street network of Wienerwald, Austria. It is an urban area with a high road density. The network is a grid of perpendicular streets. Each intersection represents a node of the network. Because of the high road density, there is also a high number of nodes. Therefore, only a small part of downtown San Francisco is selected in order to keep the network more compact. Modeling an urban street network may be challenging because there are many factors that are hard to predict. Traffic flow is on top of the agenda in this context. When searching for the fastest route, the current traffic situation and the timing of traffic lights are main concerns, other factors among many others may different routes highlighted in bold be being zebra crossings, stop signs, and bus stops. They influence not only the journey time but also the energy consumption. Nevertheless, the impacts of traffic and road signs are neglected for simplification in the present work. The focus is on the energy consumption and battery lifetime, which are expected to be considerably impacted by the topography of the city. Multi-objective optimization is performed with δ = 0.

Route: from Pier 39 to Russian Hill
Different starting and destination points are tested. One journey starts close to Pier 39, at the corner Beach Street/Grant Avenue. The destination point is Union Street/Hyde Street, close to Russian Hill. All combinations of the multi-objective optimization weights result in the same route (see Fig. 11). Going in the opposite direction, from Russian Hill to Pier 39, gives the result seen in Fig. 12, for all values of γ . Figure 13 (dashed lines) shows the topography of the two scenarios: Going from Pier 39 to Russian Hill (Union/Hyde Street) and the return. The energy consumption on both trips (see Fig. 13) strongly depends on the topography. Since Pier 39 is almost at sea level, going to Russian Hill takes a lot of energy, while the electric vehicle regenerates energy by going in the Further results of San Francisco show that focusing on the battery lifetime often results in routes differing from those focusing on energy efficiency. Switching the destination point to another one close to the original results in the three routes represented in Fig. 15. Figures 16 and 17 show the topography as well as the energy consumption and absolute energy consumption of the three results. Route (a) is the result for γ = 0 and γ = 0.2. It is longer and not as energy efficient as route (b) and (c), but shows better performance in terms of battery lifetime due to a different topography. Route (c) is the most energy efficient route (γ = 0.8 and γ = 1).

Conclusion and outlook
This work shows that route planning specifically designed for electric vehicles has more to offer than conventional routing systems, which only consider time or distance to the destination. The multiobjective optimization has three aspects : energy consumption, journey time, and battery lifetime. With different weights on the optimization variables, different solutions are obtained. Tests are performed for existing street networks. In some cases, there is only one optimal route for all variables, while other cases show very distinct results for one of the optimization variables.
The results demonstrate the influence of the topography on the routes. The energy consumption depends significantly on the topography. Especially with electric vehicles, which may regenerate energy by driving downhill, there is a strong correlation between energy consumption and topography. The third aspect of the multi-objective optimization is to increase the battery lifetime. The approach in this paper evaluates battery degradation by minimizing the total energy flow of the battery and therefore the number of charging cycles over the battery's lifetime, which is different to arriving at the destination with the highest state of charge possible, which is essentially the first aspect of the multi-objective optimization. Routes with high topographical variation increase stress on the battery and are therefore avoided with focus on the battery lifetime. The battery is an expensive part and key element of the vehicle and, as such, should be protected from degradation as much as possible. The influence of additional loads such as heating, air condition, lights, and fan is significant. Those loads are powered by the battery and therefore increase the total energy consumption, with an important factor being the duration of the journey, as more energy is used the longer these additional loads operate. Accordingly, the results of the multi-objective optimization may change, due to fast routes becoming more energy efficient.
Using a weighted multi-objective optimization may be more reasonable than using a single-objective optimization. If the second fastest route were much more energy efficient, but only slightly slower than the fastest route, a time-only optimization leads to an unreasonable result, but when different aspects are considered at the same time, the most practical solution is achieved. In order to use this type of route planning for real-world scenarios, more accurate street data is necessary as well as some additional information: Traffic lights and stop signs should be included, as well as traffic flow, as a theoretically optimal route may in practice prove to be infeasible due to traffic holdups. Integration of the multi-objective optimization algorithm (knowing vehicle type and ambient conditions) into existing route planning systems with better street data could be an addition appreciated by EV owners. Future work on this topic should also include charging stations such that route planning is possible for a journey farther than the battery range of the vehicle. Lastly, another important aspect to be included in future work is driving behavior. Speed and acceleration influence the energy consumption of the vehicle. In urban areas, the journey time is mainly dependent on the traffic situation rather than the speed, while in non-urban areas the speed may have a significant impact. Assuming a certain speed profile could help achieving realistic results. Either a statistical approach to obtain a speed profile, or learning from previously obtained data of the driver's preferences could be utilized.
Funding Open access funding provided by TU Wien (TUW).

Compliance with ethical standards
Conflict of interest The authors declare that they have no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creativecommonshorg/licenses/by/4.0/.

Appendix 1. Time-optimal route planning
Compared with the search for an energy-optimal path, finding a time-optimal path is a lot less complicated. On the one hand, there are no constraints assumed concerning the journey time. It is just about finding the fastest itinerary. On the other hand, there are no negative edge costs. In fact, the graph is undirected if it is assumed that the same road with the same speed limit is available for the way back. This optimization problem could be solved with a less complex and faster algorithm than Bellman-Ford or Yen, because there are no negative edge costs. Although, when energy and time are combined to multi-objective optimization, it is easier to just use the same algorithm, because the computation time is not crucial in this work. The procedure is the same as for finding an energy-optimal path. In this case, the edge costs represent the journey time. The journey time of an edge connecting node i and j is the sum of M adjacent road sections t ij = M l=1 t l , (A1.1) with each road section having a designated speed level v l and the length s l : Iteration 0 is given as where t 1i are the edges costs from node 1 to node i = 2, . . . , N. For odd iterations, the minimum is found using Appendix 2. Energy-optimal route planning in order to increase battery lifetime The third criterion, which may be used for route planning is increasing battery lifetime. Battery wear-off should be avoided for as long as possible, as it can decrease the capacity as well as limit the range of the vehicle. The battery is also a very expensive part of the electric vehicle to replace. There are some actions to be taken in order to increase the battery lifetime. First, avoiding deep discharge is helpful. This can be included in the factor a, which can be tuned to a value that is safe for the battery. A lithium-ion battery, which is the power source for most electric vehicles, usually has a specific cyclic lifetime. Decreasing the total number of cycles of the battery helps to increase its lifetime. This is the third criterion of the multi-objective optimization. It is performed in a similar way as the energy-optimal planning, but takes the absolute value of the energy as the edge costs. Power provided by the battery for driving or for accessories is treated the same as regenerated energy from braking and driving downhill, thus increasing the costs instead of decreasing them. The edge costs from node i to j is a ij , the sum of the absolute energy values of all road sections belonging to this edge. Similar to (A1.1), a ij is calculated a ij = M l=1 a l , (A2.1) with a l = |E l |. h i are the costs from node 1 to node i. Iteration 0 is given as where a 1i are the edges costs from node 1 to node i = 2, . . . , N. For odd iterations, the minimum is found using