Evaluation of Multi-Objective Optimization Techniques for Resilience Enhancement of Electric Vehicles

: The pervasiveness of electric vehicles (EVs) has increased recently, which results in the interdependence of power and transport networks. Power outages may adversely impact the transportation sector, and the available energy may not be sufﬁcient to meet the needs of all EVs during such events. In addition, EVs will be used for diverse purposes in the future, ranging from personal usage to emergency response. Therefore, the allocation of energy to different EVs may have different degrees of societal-, community-, and individual-level beneﬁts. To capture these diverse aspects, the energy allocation problem to EVs during outages is modeled as a multiobjective optimization (MOO) problem in this study. Three indices are formulated to quantify the value of different EVs for societies, communities, and individuals during outages, and, correspondingly, three objective functions are formulated. The formulated MOO problem is solved using the ﬁve most widely used MOO solution methods, and their performance is evaluated. These methods include the weighted-sum method, lexicographic method, normal boundary intersection method, min–max method, and nondominated sorting genetic algorithm II. To compare the performance of these methods, two indices are proposed in this study, which include the demand fulﬁllment index and total demand fulﬁllment index. The former is for analyzing the demand fulﬁllment ratio of different priority EVs, while the latter is for the demand fulﬁllment analysis of the whole EV ﬂeet requiring a recharge. In addition, the computational complexity, variance, and additional constraints required by each method are also analyzed. The simulation results have shown that the lexicographic method has the best performance when the relative priorities are known, while the min–max method is the most suitable method if the priorities are not known.


Introduction
Multi-objective optimization (MOO) is generally used to solve problems involving two or more conflicting objectives, and most of the problems in the real world are multiobjective in nature [1]. In contrast to single-objective optimization problems, a single solution generally cannot minimize/maximize all the objective functions simultaneously in MOO problems [2]. Instead, several solutions can be considered as optimal considering trade-offs among the different objectives. The set of all these possible optimal solutions form a Pareto front. A solution is known as Pareto optimal if none of the objective functions can be improved without deteriorating other objective functions. All Pareto optimal solutions are considered equally good if no preference information is provided [3]. Due to the ability to provide acceptable solutions for various complex real-world problems, MOO is used in several fields, including power systems and microgrids [4][5][6][7][8].
each for the societal level, community level, and individual level benefits of allocating energy to a particular EV. Then, one index is formulated for each of the three benefit levels mentioned above. Finally, EVs are ranked based on their value for society, the community, and individuals using the three indices. Five major MOO methods are selected, and their performance is analyzed for allocating energy to EVs during power outages. At least one method is selected from the three major categories of MOO, which includes methods with prior, posterior, and no articulation of preferences. In addition, the performance of the NSGA-II is also analyzed. To compare the performance of these methods, two indices are also proposed in this study: the demand fulfillment index and total demand fulfillment index. The demand fulfillment index considers the relative priority of different EVs while evaluating the number of EVs with the demand being fulfilled. However, the total demand fulfillment index compares the total number of EVs requiring a recharge to the total number of EVs with the demand being fulfilled. The performance of all five methods is compared based on the proposed indices. In addition, several other factors, such as the variance, computational complexity, and additional required constraints, are also compared for all the methods.

Multi-Objective Optimization Methods and Energy Allocation to EVs
MOO, also known as vector optimization, is the process of simultaneously optimizing a set of objective functions [2]. Therefore, MOO techniques are known to be more realistic options for designers due to their ability to satisfy several objectives simultaneously [1]. However, even for a trivial problem, it is unlikely to find a single solution that optimizes all the functions at the same time [3]. In most cases, the objective functions are conflicting in nature, i.e., optimizing one can deteriorate the other and vice versa. Therefore, the goal of MOO techniques is to find a trade-off between these conflicting objectives. There could be several trade-off solutions for each MOO problem, and these solutions are known as Pareto optimal solutions. A solution is known to be Pareto optimal if none of the objective functions can be improved without deteriorating other objective functions. The set of all the Pareto optimal solutions is known as the Pareto front. All the solutions in the Pareto front are considered equally good and are considered optimal if there is no specific preference for any particular solution.
For the sake of visualization, an overview of a hypothetical bi-objective optimization problem is shown in Figure 1. In Figure 1, the anchor points refer to the optimal points of the individual objective functions. In addition, the anchor points refer to the worst value for the other objective function in the case of bi-objective functions [17]. The utopia point is the ideal point where both the objective functions are optimized, which is generally unattainable. However, several methods try to minimize the distance from the utopia point. Finally, all the points on the Pareto front are attainable.
ing outages. The developed MOO problem comprises three objective functions, one for the societal level, community level, and individual level benefits of allocating en to a particular EV. Then, one index is formulated for each of the three benefit levels m tioned above. Finally, EVs are ranked based on their value for society, the commu and individuals using the three indices. Five major MOO methods are selected, and t performance is analyzed for allocating energy to EVs during power outages. At least method is selected from the three major categories of MOO, which includes methods w prior, posterior, and no articulation of preferences. In addition, the performance of NSGA-II is also analyzed. To compare the performance of these methods, two indices also proposed in this study: the demand fulfillment index and total demand fulfillm index. The demand fulfillment index considers the relative priority of different EVs w evaluating the number of EVs with the demand being fulfilled. However, the total mand fulfillment index compares the total number of EVs requiring a recharge to the number of EVs with the demand being fulfilled. The performance of all five method compared based on the proposed indices. In addition, several other factors, such as variance, computational complexity, and additional required constraints, are also c pared for all the methods.

Multi-Objective Optimization Methods and Energy Allocation to EVs
MOO, also known as vector optimization, is the process of simultaneously opti ing a set of objective functions [2]. Therefore, MOO techniques are known to be mor alistic options for designers due to their ability to satisfy several objectives simultaneo [1]. However, even for a trivial problem, it is unlikely to find a single solution that o mizes all the functions at the same time [3]. In most cases, the objective functions are flicting in nature, i.e., optimizing one can deteriorate the other and vice versa. There the goal of MOO techniques is to find a trade-off between these conflicting object There could be several trade-off solutions for each MOO problem, and these solutions known as Pareto optimal solutions. A solution is known to be Pareto optimal if non the objective functions can be improved without deteriorating other objective functi The set of all the Pareto optimal solutions is known as the Pareto front. All the solut in the Pareto front are considered equally good and are considered optimal if there i specific preference for any particular solution.
For the sake of visualization, an overview of a hypothetical bi-objective optimiza problem is shown in Figure 1. In Figure 1, the anchor points refer to the optimal poin the individual objective functions. In addition, the anchor points refer to the worst v for the other objective function in the case of bi-objective functions [17]. The utopia p is the ideal point where both the objective functions are optimized, which is gener unattainable. However, several methods try to minimize the distance from the ut point. Finally, all the points on the Pareto front are attainable.

Categorization of Multi-Objective Optimization Techniques
The MOO techniques can be divided into three major groups based on the articulation of preferences. These groups comprise methods with a prior articulation of preferences, methods with a posterior articulation of preferences, and methods with no preferences. In the first case (prior articulation), the preference or priority of different objective functions is known in advance and this information is provided to the optimization algorithm. In the case of the second group (posterior articulation), the Pareto front is obtained first and then an optimal solution is selected from the Pareto front. In addition to these three categories, multi-objective genetic algorithms are also widely used for solving MOO problems. Therefore, an overview of all four of these categories and their representative examples is shown in Figure 2. The performance of all the methods shown in Figure 2 is compared for an energy allocation problem to EVs during system contingencies.

Categorization of Multi-Objective Optimization Techniques
The MOO techniques can be divided into three major groups based on the articulation of preferences. These groups comprise methods with a prior articulation of preferences, methods with a posterior articulation of preferences, and methods with no preferences. In the first case (prior articulation), the preference or priority of different objective functions is known in advance and this information is provided to the optimization algorithm. In the case of the second group (posterior articulation), the Pareto front is obtained first and then an optimal solution is selected from the Pareto front. In addition to these three categories, multi-objective genetic algorithms are also widely used for solving MOO problems. Therefore, an overview of all four of these categories and their representative examples is shown in Figure 2. The performance of all the methods shown in Figure 2 is compared for an energy allocation problem to EVs during system contingencies.

Power Contingencies and EVs
The penetration of EVs is increasing across the globe in order to reduce the carbon emissions from the transport sector. It is predicted in the Global EV Outlook 2021 report that the size of the EV fleet will reach between 145 and 230 million by 2030 [18]. In the near future, EVs will not only be used for the personal commute but also emergency services, for example, by firefighters and hospitals. In addition, EVs will also be used for evacuation during major outages. During power outages, the locally available energy may not be sufficient to fulfill the needs of all the EVs. Therefore, different factors need to be considered to prioritize EVs and allocate the available energy. Recently, the allocation of power to EVs during major outages has been considered in several studies. For example, a non-cooperative game model is proposed in one case [16], and the energy fulfillment of higher priority EVs is ensured first by including a priority factor in the utility function of the EVs. In another study [14], three objectives are formulated for allocating energy to EVs. These objectives are social welfare, community wellbeing, and individual satisfaction. In another study, the EVs are prioritized and energy allocation is scheduled [19] during both normal congestion and power outages.
Meanwhile, the integration of stationary energy storage and renewables in fastcharging stations is also being studied [20][21][22]. These local resources can reduce the grid impacts of EVs during normal operation and can enhance the resilience of EVs during

Power Contingencies and EVs
The penetration of EVs is increasing across the globe in order to reduce the carbon emissions from the transport sector. It is predicted in the Global EV Outlook 2021 report that the size of the EV fleet will reach between 145 and 230 million by 2030 [18]. In the near future, EVs will not only be used for the personal commute but also emergency services, for example, by firefighters and hospitals. In addition, EVs will also be used for evacuation during major outages. During power outages, the locally available energy may not be sufficient to fulfill the needs of all the EVs. Therefore, different factors need to be considered to prioritize EVs and allocate the available energy. Recently, the allocation of power to EVs during major outages has been considered in several studies. For example, a non-cooperative game model is proposed in one case [16], and the energy fulfillment of higher priority EVs is ensured first by including a priority factor in the utility function of the EVs. In another study [14], three objectives are formulated for allocating energy to EVs. These objectives are social welfare, community wellbeing, and individual satisfaction. In another study, the EVs are prioritized and energy allocation is scheduled [19] during both normal congestion and power outages.
Meanwhile, the integration of stationary energy storage and renewables in fastcharging stations is also being studied [20][21][22]. These local resources can reduce the grid impacts of EVs during normal operation and can enhance the resilience of EVs during outages. However, due to the absence of connection with the grid and scarcity of locally available energy, additional factors need to be considered to prioritize EVs during outages. These factors could be the required energy, occupants of EVs, trip purpose, energy level, and time by which recharge is required. Based on these factors and the available energy, different objectives could be formulated. For example, the allocation of energy to all EVs provides individual satisfaction, but the community wellbeing and social welfare depend on the trip purpose and the occupants of the EVs since it is predicted that EVs will be the major source of transportation in the future and will be used for diverse purposes, ranging from personal use to public transport and emergency response. Therefore, these objectives need to be considered in allocating energy to EVs during power outages.

MOO for Power Allocation to EVs
To allocate energy to EVs during contingencies, a MOO problem is formulated in this section, similar to [14]. The overall objective of the MOO problem is to maximize the satisfaction of EV owners/communities using three objective functions as given by (1). The problem is formulated in such a way that the first objective function (2) has the highest precedence; it is followed by the second objective function (3), and the third objective function (4) has the least precedence. The first objective refers to societal-level benefits that can be achieved by allocating energy to a certain EV. The second objective refers to community-level benefits, and, finally, the third objective refers to individual-level benefits, which can be achieved by allocating energy to a certain EV. I 1 n , I 2 n , I 3 n represent the rank of n th EV corresponding to the first, second, and third objectives. These indices are formulated based on the trip purpose, occupants of EVs, departure time, and the difference between available and required energy. The vector x in Equations (2)-(4) is a binary vector that indicates whether the required energy of EVs is fulfilled or not. As indicated by (5), the binary variable takes a value of 1 if the energy allocated to the n th EV (E alloc n ) is greater than the required energy (E req n ). Otherwise, the binary variable corresponding to the n th EV takes a value of zero. Similarly, E de f n represents the energy deficit for n th EV and is computed as in (6). Equation (7) implies that energy allocated to each EV should be less than or equal to the amount of energy required by that EV. Finally, Equation (8) implies that the total amount of allocated energy to all EVs should be less than or equal to the available energy (E avail ). where,

Overview of Selected MOO Techniques
The performance of five MOO techniques is analyzed in this study for the energy allocation problem to EVs as described above. At least one method is selected from the four categories (prior articulation, posterior articulation, no articulation, and genetic algorithms). In addition, two indices are proposed in this study to compare the performance of these techniques and determine the most suitable method for the said problem. Therefore, in the following sections, an overview of all five techniques along with the proposed indices are presented.

Weighted-Sum Method
The WS method belongs to the methods with a prior articulation of preferences, and it is known to be the most widely used MOO method [17]. In this method, all the objectives are combined to form a single-objective function, as given by (9). Different weights (w i ) are assigned to each objective function. Generally, the weights need to fulfill certain conditions, as shown in (10). If the weights are positive and the problem is convex, this method can find the global optimal point.

Lexicographic Method
The lexicographic method (LEX) also belongs to the group of methods with a prior articulation of preferences. In this method, the objective functions are arranged in the order of their priority. Then, the most important objective function is solved first, as shown in (11). Then, the second objective function is solved and an additional constraint is added, as given by (12). The additional constraint implies that the first objective function is not allowed to deteriorate from its optimal point (x j * ). This process is repeated till all the objective functions are solved or a unique solution is obtained [2]. It can be observed that, in this method, i − 1 constraints are added after solving each objective function, where "i" refers to the number of the current objective function.

Normal Boundary Intersection Method
The normal boundary intersection (NBI) method belongs to the group of methods with a posterior articulation of preferences [23]. A generalized formulation for the normal boundary intersection method is shown in Equations (13)- (16). Φ is the payoff matrix of size n x n. Each entry in the payoff matrix can be computed as outlined in (14), where the diagonal elements will be zero. w is a vector of positive scalers, and the sum of all elements in each row is equal to 1. The quasi-normal vector n can be computed using (16), where e is a column vector of ones. In the case of NBI, the solution point is independent of how the objective functions are scaled [17]. max x∈S β Subject to Φ · w + β · n = F(x) − F 0 (14) Equations (4)-(7)

Min-Max Method
This min-max (MM) method belongs to the category of MOO algorithm with no articulation of preferences. In some cases, it is difficult for the designers to articulate preferences or provide a preference for any particular objective function. This method provides a way to solve such MOO problems without any articulation of preferences. The MM method can be solved by defining an additional variable β as given by (17). Then, all the objective functions are constrained to be not below the variable β, as given by (18).

Nondominated Sorting Genetic Algorithm II
The nondominated sorting genetic algorithm II (NSGA-II) [24] is the revised version of the original NSGA [25]. The off-spring values (c1, c1) can be generated from parent values (p1, p2) using Equations (19) and (20). The value β is computed using (21), where u is a random number between 0 and 1 while n c is the user-defined crossover parameter [17]. Finally, the crowding distance is computed using (22).

Performance Evaluation Indices
To compare the performance of all five MOO techniques, discussed in the previous sections, two indices are proposed in this study. The objective of the energy allocation problem to EVs is to maximize the number of EVs with fulfilled energy demands. Therefore, the first index is formulated to evaluate the number of EVs under priority level i with demand being fulfilled, as given by (23), where N s i is the number of EVs in priority level i with energy demand fulfilled under MOO scheme s. Higher values of demand fulfillment index (DI i ) represent a higher number of EVs with demand being fulfilled and vice versa. The second index determines the total number of EVs with demand being fulfilled and is named as total demand fulfillment index (TI). In (24), N s is the total number of EVs with energy demand fulfilled under MOO scheme s.

Simulation Results
In this section, the performance of the five MOO techniques, discussed in the previous section, is analyzed individually for the energy allocation problem. Specifically, the Pareto front of the different MOO methods is discussed. In the next section, several performance metrics are used to compare the performance of all five techniques. NSGA-II was implemented in Java [26] and CPLEX [27] was integrated with Java for carrying out the optimization in the remaining four methods.

Input Data
The data required to solve the MOO problem are the required energy by each EV and the values of three indices for each EV, as shown in Equations (2)-(4). For the sake of visualization, a fleet of 15 EVs is considered, similar to [14]. The rank factors of all the EVs and required energy by each EV are taken from [14]. An overview of the EV-wise rank and the required energy are shown in Table 1. The zero value in the table implies that those EVs are not contributing to that particular index due to the nature of the task and occupants of the EV. The total energy required by all the EVs is 439.5 kWh, and the available energy is 220 kWh.

Weighted-Sum Method
In the case of the WS method, the weights of all three functions are varied and five extreme cases are simulated, as shown in Table 2. The values of the objective functions are recorded for each case and the results are shown in Figure 3. For the sake of visualization, three pairs of the three objective functions (F1, F2, and F3) are made and the results are shown in Figure 3a-c. In addition, Figure 3d gives an overview of all three of the functions together. It can be observed from the result that the outcome changes with the change in weights, which is a drawback of the WS method. In the case of the WS method, determining the optimal weights for each method is another optimization problem. If the weights are not optimal, the outcome may not be optimal and may also not be in accordance with the preferences of the designers and policymakers.

Weighted-Sum Method
In the case of the WS method, the weights of all three functions are varied and five extreme cases are simulated, as shown in Table 2. The values of the objective functions are recorded for each case and the results are shown in Figure 3. For the sake of visualization, three pairs of the three objective functions (F1, F2, and F3) are made and the results are shown in Figure 3a-c. In addition, Figure 3d gives an overview of all three of the functions together. It can be observed from the result that the outcome changes with the change in weights, which is a drawback of the WS method. In the case of the WS method, determining the optimal weights for each method is another optimization problem. If the weights are not optimal, the outcome may not be optimal and may also not be in accordance with the preferences of the designers and policymakers.

Lexicographic Method
In the case of the LEX method, the priority of all three of the objective functions needs to be defined. The LEX method is not prone to the weights of the individual functions, and it gives the same outcome irrespective of the magnitude of the weights; it only compares the relative values of each function. To cover all possible cases, four cases are simulated by considering the different priorities of the three objective functions (F1, F2, F3), as shown in Table 2. In the first case, all three objectives are given the same priority, while, in the other three cases, each objective is given the highest priority, respectively. The results of all four cases are shown in Figure 4. In the case of the LEX method, only these four unique outcomes were recorded. All other combinations have given the same results due to the consideration of the relative priority of different objectives in the LEX method.

Normal Boundary Intersection Method
The NBI is the method with a posterior articulation of preferences. Therefore, all possible combinations of different weights need to be evaluated first and the policymakers will choose one of the outcomes. Therefore, the weights of all three of the functions (F1, F2, and F3) are varied in the range of 0.8 to 0.1 and the corresponding value of the optimization parameter β is monitored. The values of the weights and the corresponding value of β for each run are shown in Table 3. Some of the combinations have resulted in infeasible outcomes and those combinations are omitted from the results. lated by considering the different priorities of the three objective functions (F1, F2, F3), as shown in Table 2. In the first case, all three objectives are given the same priority, while, in the other three cases, each objective is given the highest priority, respectively. The results of all four cases are shown in Figure 4. In the case of the LEX method, only these four unique outcomes were recorded. All other combinations have given the same results due to the consideration of the relative priority of different objectives in the LEX method.

Normal Boundary Intersection Method
The NBI is the method with a posterior articulation of preferences. Therefore, all possible combinations of different weights need to be evaluated first and the policymakers will choose one of the outcomes. Therefore, the weights of all three of the functions (F1, F2, and F3) are varied in the range of 0.8 to 0.1 and the corresponding value of the optimization parameter β is monitored. The values of the weights and the corresponding value of β for each run are shown in Table 3. Some of the combinations have resulted in infeasible outcomes and those combinations are omitted from the results.
The results of all three of the objective functions are shown in Figure 5 for the NBI method. It can be observed that the number of outcomes is higher than the previous two methods due to the evaluation of all feasible combinations of the three objective functions. However, it can be noticed that each combination does not produce a unique result for the objective functions. Some of the combinations give the same results, which is evident from the number of points in Figure 5 and the total number of runs, listed in Table 3. However,  The results of all three of the objective functions are shown in Figure 5 for the NBI method. It can be observed that the number of outcomes is higher than the previous two methods due to the evaluation of all feasible combinations of the three objective functions. However, it can be noticed that each combination does not produce a unique result for the objective functions. Some of the combinations give the same results, which is evident from the number of points in Figure 5 and the total number of runs, listed in Table 3. However, it provides a better picture to the policymakers for making decisions based on their preferences while analyzing all possible combinations.

Min-Max Method
The MM method produces only one unique solution due to the absence of any preferences. The energy allocated to each EV by the MM method is shown in Table 4. The binary variable ( n x ) indicated whether the energy demand of the particular EV is fulfilled or not, while def n E is the deficit energy amount. It can be observed that the MM method has allocated energy to 8 EVs, while 7 EVs have a deficit of energy. However, due to the

Min-Max Method
The MM method produces only one unique solution due to the absence of any preferences. The energy allocated to each EV by the MM method is shown in Table 4. The binary variable (x n ) indicated whether the energy demand of the particular EV is fulfilled or not, while E de f n is the deficit energy amount. It can be observed that the MM method has allocated energy to 8 EVs, while 7 EVs have a deficit of energy. However, due to the absence of any preference, the selected EV is not based on its importance, which was the case for the previous methods. The value of β was 110.2 for the MM method.

Nondominated Sorting Genetic Algorithm II
In the case of NSGA-II, the Pareto front is obtained by running the algorithm several times and filtering the best Pareto solutions. The obtained results for the NSGA-II are shown in Figure 6. It can be observed that the NSGA-II also provides several unique solutions for the energy allocation problem to EVs. The range of objective functions for NSGA-II is slightly different from the first three methods (WS, LEX, and NBI). However, these solutions need to be analyzed for their quality considering different metrics/indices. The comparison of all five of the methods is presented in the following section.

Nondominated Sorting Genetic Algorithm II
In the case of NSGA-II, the Pareto front is obtained by running the algorithm several times and filtering the best Pareto solutions. The obtained results for the NSGA-II are shown in Figure 6. It can be observed that the NSGA-II also provides several unique solutions for the energy allocation problem to EVs. The range of objective functions for NSGA-II is slightly different from the first three methods (WS, LEX, and NBI). However, these solutions need to be analyzed for their quality considering different metrics/indices. The comparison of all five of the methods is presented in the following section.

Discussion and Analysis
In this section, the performance of the five MOO techniques discussed in the previous sections is compared. In the first part, the performance is analyzed using the proposed indices. In the second part, the complexity and the variance in the results are analyzed.

Evaluation via Performance Indices
The two indices proposed in this study are used in this section to analyze the performance of the five MOO methods. The first index (demand fulfillment index, DI) indicates the ratio of EVs with a particular rank to the total number of EVs having a non-zero value for that rank. Higher values imply more EVs with demand being fulfilled and vice versa. It can be observed from Figure 7a that DI 1 is the highest in the case of LEX and NBI, while it is the lowest for the NSGA-II. The value on DI 2 is the highest for WS and lowest for MM, while the value of DI 3 is the highest for WS and lowest for LEX. Since the priority of F1 is the highest and that of F3 is the lowest, the results of LEX and NBI will be considered better. The WS is increasing the values of DI 2 and DI 3 at the cost of DI 1 , which is generally not preferred.
Similarly, the second index (total fulfillment index, TI) shows the ratio of the total number of EVs with demand being fulfilled to the total number of EVs requiring a recharge. It can be observed from Figure 7d that the value of TI is the highest for the WS method, while it is the lowest for the NSGA-II and LEX methods. The highest value in the WS is due to the consideration of the accumulated indices in the case of the WS, while, in the case of LEX, higher priority ranks are preferred irrespective of their accumulated indices. In the case of NSGA-II, there is no specific preference for any particular rank.
If other factors are not considered, the performance of LEX and NBI is suitable for objectives with clear priority/preferences due to their ability to maximize the allocation based on the priority. The performance of NSGA-II was the lowest in most of the cases, thus making it unsuitable for this problem. Finally, the WS method is better for maximizing the total number of EVs having their demands fulfilled at the cost of some degradation in the performance of the higher rank EVs. However, other factors also need to be considered to determine the best suitable candidate method, which is discussed in the following section.
the case of LEX, higher priority ranks are preferred irrespective of their accumulated indices. In the case of NSGA-II, there is no specific preference for any particular rank.
If other factors are not considered, the performance of LEX and NBI is suitable for objectives with clear priority/preferences due to their ability to maximize the allocation based on the priority. The performance of NSGA-II was the lowest in most of the cases, thus making it unsuitable for this problem. Finally, the WS method is better for maximizing the total number of EVs having their demands fulfilled at the cost of some degradation in the performance of the higher rank EVs. However, other factors also need to be considered to determine the best suitable candidate method, which is discussed in the following section.

Complexity and Variance Analysis
In the previous section, the best-performing case was selected for each method. However, in the case of the WS, LEX, and NBI methods, the outcome changes with a change in the priority or weights. Therefore, the maximum and minimum number of EVs survived by different methods and their variance is tabulated in Table 5. It can be observed that the variance is maximum for the LEX while it is minimum for the NBI method (the MM and NSGA give single unique solutions). The higher variance in the LEX is due to a significant change in the outcome with the change in the priorities of the objective functions. In terms of additional constraints, only one additional constraint is required for WS, while NBI and MM require I number of additional constraints, where I is the number of objective functions. Similarly, in the case of LEX, I-1 additional constraints are required, while NSGA-II does not require any additional constraints. An increase in the number of objective functions could increase the complexity of the problem where additional constraints are required.
The computational complexity of WS is the lowest due to the requirement of scalers only and a single additional constraint, while the complexity of NBI is higher due to the requirement of several runs and having I number of additional constraints. Similarly, the computational complexity of NSGA-II is also higher due to the requirement of more time in finding the optimal solution. MM and LEX have a medium level of complexity, where MM requires I additional constraints but provides a single solution. LEX requires I-1 additional constraints and provides different solutions based on the priority.
The analysis in this section and the previous section shows that NSGA-II is not suitable for this problem due to lower performance despite having higher computational complexity. LEX is the most suitable method if the priority of different EVs is known due to its ability to allocate energy based on the priority, similar to the findings of [14]. In addition, it has a medium level of computational complexity. MM is suitable for cases where the priority is not known, and it can provide a single unique solution with a medium level of computational complexity.

Conclusions
A multi-objective optimization problem was developed for allocating energy to electric vehicles during power outages. Three objective functions were developed based on three indices to quantify the benefits of different electric vehicles for societies, communities, and individuals. The performance of five multi-objective optimization methods was analyzed for solving the formulated problem. At least one method was selected from the three major categories of the multi-objective optimization, which include methods with a prior articulation, methods with a posterior articulation, and methods with no articulation of preferences. In addition, the well-known non-dominated sorting genetic algorithm was also considered. Two indices were proposed and used to analyze the performance of all five methods. The comparison showed that the lexicographic method has the best performance if the relative preference/priority of the different objectives is known in advance. In addition, the computational complexity of this method is not too high. The min-max method also showed a reasonable performance and is suitable if the priority of different objective functions is not known. This method also has a medium level of computational complexity. The performance of the non-dominated sorting genetic algorithm was the lowest, and it had higher computational complexity as well. Therefore, it is not deemed suitable for the energy allocation problem to electric vehicles.