Application of Hexagonal Fuzzy MCDM Methodology for Site Selection of Electric Vehicle Charging Station

: In this paper, the application of hexagonal fuzzy multiple-criteria decision-making (MCDM) methodology for the site selection of electric vehicle charging stations is considered. In this regard, four factors and thirteen sub-factors have been taken into consideration for E-vehicle charging site selection. In this research, the geographic information system (GIS) has been incorporated with MCDM techniques. The fuzzy analytic hierarchy process (FAHP) is used to obtain a fuzzy weight of factors and sub-factors. MCDM tools fuzzy technique for order of preference by similarity to ideal solution (FTOPSIS) and fuzzy complex proportional assessment (FCOPRAS) have been used to rank the selected sites. A centroid-based method for defuzziﬁcation and distance measure between two hexagonal fuzzy numbers (HFN) has been developed for this paper. A practical example in Howrah, India, is considered to show the applicability and usefulness of the model. The results depict the suitability of the proposed research. Comparative and sensitivity analyses have been demonstrated to check the reliability, robustness and effectiveness of the proposed method.


Introduction
Electric vehicles play a very momentous role in addressing fossil fuel pollution and they are capable of making a paradigm shift in the entire transportation sector. Transportation is a significant contributor to urban air pollution and the reduction of urban city emission is the need of the hour. Electric vehicles make the world more liveable and provide a pollution-free mode of transportation in urban areas. The high level of pollution is degrading the environment and it has made the concept of sustainable development a fairy-tale phenomenon [1,2]. Sustainable consumption needs to be adopted by conducting timely environmental and sustainability assessments to prevent any large-scale ecological disaster from happening. This is where electrical vehicles come into the picture, with sustainability, environment-protecting and pocket-friendly being a few of their rewards. Electric vehicles use a minimum of one electric motor or traction motor for propulsion. They maybe self-contained with a generator or battery for converting the fuel into electricity or they may be power-drivenvia a collector scheme by using electricity from off-vehicle sources [3]. The problem of the energy crisis in the world can be tackled in the future by using this option. Ever-rising gas prices force people to look for alternative modes of

Literature Review
The construction of an electric vehicle charging station is the turning point in its life cycle. Selecting an appropriate place for setting up the charging station and determining its capacity will benefit all the stakeholders and endorse sustainable development of the entire industry. Some studies have been shown this with respect to economy and technology. Lee et al. [28] examined the price competition among EVCSs using renewable power generators by employing game theory with relevant physical constraints. Rivera et al. [29] put forward a novel architecture for plug-in electric vehicles: a DC charging station using a grid-tied neutral point clamped converter. Wang et al. [30] utilized threephase uncontrolled rectification chargers to study the harmonic amplification of EVCSs. Ding, Hu and Song [31] studied the energy storage system as a potential supplement of an electric bus fast charging station by employing mixed integer nonlinear programming for valuing the energy storage system. Fan, Sainbayar and Ren [32] calculated the effect of limiting electric vehicles' full state of charge to total charged energy and the revenue of EVCS. Capasso and Veneri [33] built a DC fast charging architecture for plug-in hybrid vehicles as well as fully electric vehicles by integrating the fleets of hybrid/road electric vehicles with renewable energy sources. Li et al. [34] studied the control of EVCSs and the management of energy and fit them in a dynamic price framework by developing a realtime simulation system for evaluating how the EVCS meets the charging and discharging requirements for grid-to-vehicle, vehicle-to-grid and vehicle-to-vehicle. Nansai et al. [35] conducted the life-cycle analysis on EVCSs in three phases of installation, transportation and production of electric vehicle charging equipment and then compared the carbon dioxide, carbon monoxide, sulphur oxides and nitrogen oxide emissions of electric vehicles and gasoline vehicles.
Research also focuses on the electric vehicle's size, its placement and the fields in which it is used. Khalkhali et al. [36] determined the optimal location and size of plug-in hybrid EVCSs using data envelopment analysis by maximizing the benefit of the distribution system management. Frade et al. [15] studied the location of EVCSs in Portugal by employing a maximal covering model for defining the capacity and number of EVCSs. Liu [37] investigated the nascent electric vehicle market in Beijing and then formulated an assignment model for different charging infrastructure assignment strategies. Wirges, Linder and Kessler [38] predicted several scenarios for charging infrastructure development until 2020 in Germany and also formulated a dynamic spatial model for the allocation of the EVCSs in the country. Wang et al. [39] developed an EVCS location model considering the electricity consumed along the roads and the oil sales, after which the EVCS quantity and layout were calculated. He et al. [40] developed an equilibrium modeling framework for deciding the optimal allocation of the charging stations in a metropolitan area. Liu et al. [41] proposed a two-step screening method for determining the optimal EVCS site, taking environmental factors and service radius into consideration.
There have been a large number of researches which have employed a type of decisionmaking method. Pashajavid and Golkar [42] put forward a scenario optimization algorithm by using multivariate stochastic modeling methodology for load demand by allocating the charging station of plug-in electric vehicles and also used a particle swarm optimization algorithm to minimize voltage deviation and energy loss in the distribution system. Chen, Kockelman and Khan [43] determined the optimal charging station location assignment in Seattle by developing a mixed-integer programming model which minimized the station access cost of the electric vehicle users and took trip attributes, parking demand, population density and local job as constraints. Sathaye and Kelley [44] used a continuous facility location model considering the demand uncertainty for finding the optimal location of the public-funded electric vehicle infrastructure on the highway corridors. Wang et al. [45] determined the optimal size and location of EVCSs with respect to smart grids by proposing a multi-objective EVCS planning method, maximizing the electric vehicle traffic flow under the constraint of battery capacity and the final optimal solution data-envelopment analysis was obtained and employed. Dong, Liu and Lin [46] formulated an activity-based assessment method for the evaluation of the feasibility of electric vehicles taking the heterogeneous travelling population and subsequently applying a genetic algorithm for determining the sub optimal location. Xu et al. [47] identified the candidate centralized charging station using a proposed mathematical model with the minimum total transportation distance. You and Hsieh [48] used the mixed-integer programming model to determinine the best location of the EVCS which would maximize the number of people who could complete round-trip itineraries, along with developing an efficient hybrid genetic algorithm which would obtain a compromised solution in a reasonable time. Lee et al. [49] first collected users' charging and traveling behaviors along with the batteries' state of charge and then proposed a location model of the rapid EVCS by using a probabilistic distribution function for the remaining fuel range. Baouche et al. [50] formulated an integer linear optimization model, taking the electric vehicle's input consumption as the optimal model. Sadeghi-Barzani [16] postulated a mixed-integer non-linear optimization approach for the determination of the optimal size and place of fast charging stations. Yao et al. [51] devised a multi-objective collaborative planning strategy to tackle the planning issue in electric vehicle charging systems and integrated power distributions by using a decomposition-based multi-objective evolutionary algorithm model and equilibrium-based traffic assignment model. Ma et al. [52] used an agent-based model to optimize the sketch of the initial EVCS. Chung and Kwon [53] devised a multi-period optimization model to perform the EVCS planning, basing it on the real traffic flow data of the Korean Expressway network in 2011. Lam, Leung and Chu [54] used four solution methods i.e., the greedy approach, chemical reaction optimization, iterative mixed-integer linear programming and effective mixed-integer linear programming for the EVCS place problem by checking them against practical and artificial cases. Cai et al. [55] found the relation between public charging infrastructure development and the travel patterns mined from big-data taken using large-scale trajectory data in Beijing.
Multi-criteria decision-making was used in transhipment site selection [56], shopping mall site selection [57], railway station site selection [58], waste management site selection [59] and macro-site selection of solar/wind hybrid power station [60].
The major drawback of the aforementioned studies which use multi-criteria decisionmaking approaches is the lack of incorporation between evaluation criteria and spatial data, since selecting a suitable site for an EVCS is a spatial decision predicament which involves conflicting, incommensurate and multiple evaluation criteria and a large set of evaluation criteria. Two unique parts of research i.e., multi-criteria decision-making and GIS can aid each other to overcome the intersections between evaluation criteria and spatial data. GIS integrates the spatially-referenced data in a problem-solving environment whereas multi-criteria decision-making refers to procedures and techniques used in structural decision problems by evaluating, prioritizing and designing the alternative decisions [61]. GIS-based multi-criteria decision-making techniques are usually used for spatial decision problems such as freight village site selection [62] wind power plant site selection [63] and refugee camp site selection [64]. Even though there have been multiple studies on the GIS-based multi-criteria decision-making approach on spatial decision problems as mentioned above, still there is a gap in the selection of EVCS sites using GIS.

Objectives of the Study
The present research has the following objectives: (a) Identification of the most preferred site for the construction of an e-vehicle charging station. (b) Application of hexagonal fuzzy numbers (HFN) in AHP-TOPSIS and AHP-COPRAS to obtain the ranking of the selected sites.

Novelties of the Study
Several researchers have explored fuzzy numbers with MCDM techniques AHP, TOPSIS, and COPRAS. Hardly any research has been done using hexagonal fuzzy numbers under MCDM methodology. The HFN defuzzification formulae has been developed and utilized. Distance measured between two HFN is also defined. The formulae has been developed to calculate the hexagonal fuzzy weight of factors and sub-factors. A technique has been developed to incorporate more than one decision maker's opinion into a single comprehensive value in terms of HFN. Two different ranking methods, fuzzy AHP-TOPSIS and fuzzy AHP-COPRAS has been used in this research. GIS software has been used for distance measurement and graphical presentation of the selected sites.

Structure of the Paper
The remainder of the paper is organized in the following way: Section 2 depicts the concept of fuzzy numbers, HFN, its properties, distance measure, and defuzzification formulae. It also includes the MCDM technique AHP, fuzzy TOPSIS and fuzzy COPRAS. Section 3 contains the numerical application and description of factors and sub-factors. Section 4 represents the comparison and sensitivity analysis, respectively. The results and discussion covered in Sections 5 and 6 describes the future scope and conclusion.
The structured framework of the study represented in Figure 1 describes the sequential steps followed in this paper. After the initial selection of factors and sub-factors, their weights are calculated by AHP and FAHP. Subsequently, FTOPSIS and FCOPRAS MCDM tools have been applied for the selection of preferred locations followed by comparative and sensitivity analysis, overall results and discussions.

Novelties of the Study
Several researchers have explored fuzzy numbers with MCDM techniques AHP, TOPSIS, and COPRAS. Hardly any research has been done using hexagonal fuzzy numbers under MCDM methodology. The HFN defuzzification formulae has been developed and utilized. Distance measured between two HFN is also defined. The formulae has been developed to calculate the hexagonal fuzzy weight of factors and sub-factors. A technique has been developed to incorporate more than one decision maker's opinion into a single comprehensive value in terms of HFN. Two different ranking methods, fuzzy AHP-TOPSIS and fuzzy AHP-COPRAS has been used in this research. GIS software has been used for distance measurement and graphical presentation of the selected sites.

Structure of the Paper
The remainder of the paper is organized in the following way: Section 2depicts the concept of fuzzy numbers, HFN, its properties, distance measure, and defuzzification formulae. It also includes the MCDM technique AHP, fuzzy TOPSIS and fuzzy COPRAS. Section 3 contains the numerical application and description of factors and sub-factors. Section 4 represents the comparison and sensitivity analysis, respectively. The results and discussion covered in Sections5 and 6 describes the future scope and conclusion.
The structured framework of the study represented in Figure 1 describes the sequential steps followed in this paper. After the initial selection of factors and sub-factors, their weights are calculated by AHP and FAHP. Subsequently, FTOPSIS and FCOPRAS MCDM tools have been applied for the selection of preferred locations followed by comparative and sensitivity analysis, overall results and discussions.
. Figure 1. Structural framework of the study. FTOPSIS: fuzzy technique for order of preference by similarity to ideal solution; FCOPRAS: fuzzy complex proportional assessment; MCDM: multiple-criteria decision-making.
represents the membership function ofT which takes value from zero to one. In real life situations, where the information is vague and uncertain, fuzzy logic can be efficiently used to deal with these problems.  Figure 2 represents the membership function of symmetric HFN. The following Figure 2 represents the membership function of symmetric HFN.  Where h 1 , h 2 , h 3 , h 4 , h 5 and h 6 are real numbers such that Note 2.1: Hexagonal fuzzy numbers (HFN) capture the hesitancy and uncertainty in broader aspects compared to triangular fuzzy numbers (TFN), trapezoidal fuzzy numbers (TrFN), and pentagonal fuzzy numbers (PFN) as the latter undertakes three, four, and five numbers, respectively, to represent the impreciseness and ambiguity of the decision maker (DM). If we consider the linguistic term corresponding to TFN then it is represented numerically as (low, medium, high) where the medium value corresponds to the best possible chances of the quantity. For TrFN, the numerical presentation is expressed as (very low, low, high, very high) and the maximum possibility lies in a range of (low, high). Further, if we take PFN, it considers the numerical behavior as (very low, low, medium, high, very high), here the middle value denotes the maximum possibility, which is 1. HFN undertakes the highest level of distribution in the numerical form as (very very low, very low, low, high, very high, very very high), i.e., the maximum possible spread can be accommodated in HFN. With respect to a developing country in an unplanned city, many of the attributes under consideration have a wider range of linguistic representation which cannot be captured using TFN, TrFN and PFN. For example, the consumption level in a city with a heterogeneous earning pattern can be better captured using HFN with a wider range. The hesitancy of decision makers in such an environment requires a broader range of values to depict uncertainty; HFN enables this.
The hesitancy of decision makers (DMs) can be better captured in HFN. Different levels of hesitancy in linguistic terms corresponding to different HFN. It is easy to understand and analyze that the linguistic terms "weakly important" represented in HFN as (1.1, 1.2, 1.3, 1.4, 1.5, 1.6), "absolutely important" represented in HFN as (4.6, 4.8, 5, 5.2, 5.4, 5.7) and so on. These range from lower values of HFN to higher values. Similarly, the inverse of these HFN, specifically, the opposites of the linguistic terms can be obtained by using Equation (7).

Arithmetic Operations of Linear Symmetric HFN
be two HFN, then their general arithmetic operations can be defined in the following way:

Centroid-Based Method for the Defuzzification of Hexagonal Fuzzy Numbers
A HFN can be considered as the union of two triangles and two trapeziums, e.g., ∆ADH, ∆BCG, ABCD and CDEF together form a HFN.
Further, the trapezium is a union of two triangles and one rectangle. Applying a centroid-based method to triangles and rectangles, and finally summing them, we obtain the centroid of the HFN. Since, the defuzzified value should remain within the range of a HFN, the formulae given below provides the required defuzzified value. Derivation of the HFN is executed in the following way: (iii). Centroid of Trapezium ABCD: (a) Centroid of (iv). Centroid of Trapezium CDEF is calculated in the similar order and we obtain: The defuzzified value is determined by summing Equations (9), (10), (14) and (15) and dividing the denominator by the sum of the quantities of the numerator.

Determination of Hexagonal Fuzzy Weights of Factors and Sub-Factors
We have extended the methodology developed by Buckley [74] for TFNs in the context of determining hexagonal fuzzy weight.
Step 1. The geometric mean value of the HFN is obtained using: Step 2. Summation of each k c Step 3. To calculate the inverse of each k c and arrange it in increasing order.
Step 4. To find the hexagonal fuzzy weight of factors and sub-factors using the following equation: Step 5. The global hexagonal fuzzy weight of sub-factors are computed by the product of factor weight with the respective sub-factor fuzzy weight.

Fuzzy Analytic Hierarchy Process (FAHP)
The AHP was introduced by Satty [75]. It is used widely for the evaluation of factor and sub-factor weights. AHP helps in solving the real-life situations with a scientific approach. The comparison of factors and sub-factors, thereby giving preference in linguistic terms can be considered as a hesitant task for DMs, thus HFN with AHP methodology captures the vagueness of the problem. The determination of factors' and sub-factors' weights are important for ranking the electric vehicle charging station. AHP works with a problem hierarchy where a comparison matrix is constructed to represent subjective judgments regarding criteria and sub-criteria. In this work, FAHP is taken instead of AHP, keeping in mind the fuzzy setting represents the uncertainties of the decision experts. The FAHP concept with fuzzy logic allows the DMs in the evaluation of reliable results. The steps of FAHP are given below.
Step 1. Construction of a comparison matrix in terms of HFN by a group of decision experts.
Let a group of 'H' decision-makers assigned for the comparison of factors and subfactors. Let each DM express their preference in the pairwise comparison of factors and sub-factors. Thus, 'h' set of matrices are obtained, T h = {t cdh }.
Where t cdh = ( m cdh , n cdh , o cdh , p cdh , q cdh ) denotes the HFN of c factor to d factor as expressed by the 'h' DM.
Step 2. Defuzzification of HFN: A HFN can be defuzzified by using the centroid-based method used in this paper. Thus using Equation (16), the HFN is transformed to a crisp value.
Step 3. Normalization of the defuzzified matrix: Step 4. Estimation of factors' and sub-factors' weights: Step 5. To test the Consistence Index (C.I) of the matrix: where j denotes the size of the matrix.
Step 6. Determination of Consistence Ratio (C.R): where R.I is stand for Random Index, and its value differs with the size of the matrix "n".
The assessment of C.R ≤ 0.1 is acceptable and indicates that the weights obtained are justified. Thus further evaluation is not essential.

Technique for Order Preference by Similarity to Ideal Solution (TOPSIS) and (FTOPSIS)
The TOPSIS MCDM tool is an extensively used technique, developed by Hwang and Yoon [76] to rank the alternatives, thus giving an idea as to which choice is most preferred. The TOPSIS method is considered to be a distance measure method in which the optimal alternative obtained is farthest away from the negative ideal solution (NIS) and nearest to the positive ideal solution (PIS). The linguistic human decisions can be reflected suitably with Fuzzy TOPSIS (FTOPSIS). The approach is useful in handling the complexity of the situation involving several factors and their sub-factors. In this research, for the selection of the best site to construct an electric vehicle charging station, it is dependent on multiple conflicting factors and sub-factors, thus the MCDM method FTOPSIS introduced by Sodhi and Prabhakar [77] is one of the most helpful and reliable methods. The fuzzy logic extends our goal to obtain more sensitive results in this regard. The steps of FTOPSIS are described below.
Step 1: Construction of the decision matrix by the help of decision experts in terms of linguistic terms. The linguistic terms are then converted to a HFN.
Step 2: To evaluate the normalized HFN fuzzy decision matrix: where B.A and N.B.A signifies the benefit attributes and non-benefit attributes, respectively.
Step 3: To evaluate the weighted fuzzy normalized matrix, the sub-factors' fuzzy weights are multiplied with the normalized fuzzy value: where Step 4: Calculate the fuzzy positive ideal solution (FPIS) (PIS + ) and fuzzy negative ideal solution (FNIS) (N IS − ), where h + g denotes the maximum value of h gh and h − g denotes the minimum value of h gh : where M B denotes the benefit attributes and M NB denotes the non-benefit attributes.
Step 5: Calculation of the distance measure of all alternatives from the PIS and NIS. The two Euclidean distances for individual alternatives can be calculated as follows: where d(., .) denotes the Euclidean distance between two fuzzy numbers.
Step 6: Determination of the relative closeness to the ideal alternatives: Step 7: Rank the alternatives: The alternatives are ranked based on the score obtained by R g . The larger value of R g signifies the better alternatives.

Fuzzy COPRAS Methodology
The complex proportional assessment (COPRAS) method was first introduced by Zavadskas, Kaklauskas and Sarka [78]. Fuzzy COPRAS is an extended approach of the COPRAS technique, widely used for decision-making problems. It uses the stepwise ranking and evaluation procedure for the alternatives with reference to significance and utility degree. A few applications of the COPRAS method are in economics, construction and property management. An extension of the COPRAS method is Fuzzy COPRAS, which is frequently used in decision-making problems. Ghose et al. [79] used a hybrid fuzzy COPRAS method for selecting the optimal material to be used for a solar car. They took into consideration 19 materials which had 14 different properties. Using a sensitivity analysis, the robustness of the model was checked. The reason for using the fuzzy-based MCDM technique was that it helps decision makers to get over the problems of ambiguous data. Tolga and Durak [80] used the fuzzy COPRAS technique in the air cargo sector for evaluating the potential capability. 18 criteria were chosen for selecting the best out of the present six alternatives. The steps of the COPRAS method are illustrated below: Step 1. Decision matrix is constructed in terms of HFN, the alternatives are given linguistic terms by the decision experts with respect to the criteria.
Step 2. Normalized decision matrix is formulated using Equation (1), in the similar way, we constructed for TOPSIS normalized matrix.
Step 3. Weighted normalized matrix is constructed by multiplying the criteria weights with fuzzy normalized matrix using Equation (19).
Step 4. Aggregation of beneficial B + g and non-beneficial indices NB − g for each alternative are evaluated.
where h = 1, 2, . . . , m represents the benefit attribute of the alternatives and h = m + 1, m + 2, . . . , t represents the non-benefit attributes of the alternative.
Step 5. Finally, the aggregated beneficial and non-beneficial indices are defuzzified using the Equation (9) and R +g and R −g are determined.
Step 6. Calculation of R g using the following formulae: Step 7. Ranking of the alternatives are done using the formulae: where R g represents the g-thdefuzzified value and R max represents the maximum defuzzified value of individual alternative. The prosperity of a nation is comprehended by the state of its economy. The assessment of the economic factors reveals the feasibility of the undertaken study. The factors considered are:

•
Land cost (C 11 ): It is the crux of the entire planning for the optimum charging station location. Land costs are based on their use, i.e., non-agricultural urban land is more costly than agricultural land. Since the purpose is to build a charging station, we can minimize the land cost by utilizing an already existing utility station. If money is saved on the cost of procuring the land, then it can be utilized for setting up the station.

•
Operating and management cost (OMC)(C 12 ): Yao, Bai and Xu [81] stated that a significant part of the budgeted amount is cut out for dealing with operating and management costs which arise in the day to day working of the project. The minimum management cost is essential since it helps in the smooth flow of the information from one department to the other. The charging station should be executed sothat the operations can be systematically planned, which will reduce the in-between costs. Electric vehicles will reduce the air pollution, hence initial operating costs are understandable since the long-term implications outweigh the costs.

•
Consumption level (C 13 ): Modrak and Soltysova [82] studied the operational complexity measure. The measure of consumption level denotes how affluent the people in a particular locality are. In the case that consumption level in an area is high, it can be expected that people will be more willing to go to further distances in search of more options. A charging station can be built in a high consumption area since the throng of people will have more ways of traveling and procuring their wants.

•
Construction cost (C 14 ): Manerba, Mansini and Perboli [83] researched the capacitated supplier selection problem considering total quantity discount policy and activation costs under uncertainty. Construction cost varies with the location, and to make the charging station a success, the initial fixed cost should be minimized as much as possible. If the location is well-connected by various transportation facilities, then the cost of transferring the construction materials will decrease, which will decrease the construction cost and the overall profitability of the charging station will increase initially.
• Public facilities (C 15 ): Kinay et al. [84] studied multi-criteria chance-constrained capacitated single source discrete facility location problems. Public facilities refer to schools, colleges, grocery stores, shopping malls, bus stops and the other everyday amenities which are used by commoners on a mass scale. In the case of a charging station being built near a location with a large density of public facilities, it will act as a boon since money will frequently change hands and thereby develop the area.

Environmental Factors (C 2 )
This refers to those factors which will influence the immediate surroundings of the charging station. A clean and green environment helps in resonating the theme of electric vehicles and thereby makes a charging station built in such a location a success.

•
Generation of noise and air pollution (GNAP)(C 21 ): In the current scenario, noise and air pollution are considered bigger hindrances than other forms of pollution. This is because they are experienced daily, which results in greater damage due to them. Electric vehicles will help in reducing both, since the batteries of the vehicles will not cause air pollution or make noise while being on the road.

•
Petrol stations (C 22 ): The availability of petrol stations nearby signifies a greater number of vehicles in the area. Building a charging station near a petrol station will cause an increase in the number of vehicles and people visiting the area and thereby turn out to be more profitable since the cost will be less.

•
Transportation stations (C 23 ): These refer to the various bus stops, bicycle charging stations, railway stations, car stops, and tram stops. The greater the frequency of the transportation vehicles, the more will be its impact on the environment since every form of transportation will bring more people, eventually harming the environment. Electric vehicles, on the other hand, help in the easy transportation of people without harming the environment.

Traffic Factors (C 3 )
This refers to those sets of factors which are only noticeable when there is a huge population in the area. A charging station which takes into consideration the traffic factors is one which will be able to have huge implications on the state of travel in the area.

•
Number of roads (C 31 ): Yao, Bai and Xu [81] studied the impact of the number of roads on China's thermal power industry. The number of roads stands for the various options which the vehicles can take in the case of congestion or for availing a shorter reach time. If the charging station is strategically built near an intersection of heavytraffic roads, then it will help the drivers be more at ease while driving since they will have a backup nearby.

•
Road potency (C 32 ): Hosseini and Sarder [85] studied the road potency for optimal site selection using a Bayesian network model. The higher the number of vehicles in the region, the greater the footfall will be, which will eventually increase the success rate of the charging station.

•
Parking areas (C 33 ): The increase in parking areas will lead to an increase in the use of electric vehicles, since one of the major thoughts which arises in a prospective buyer's mind is where one will park the vehicle. If the charging station is built near a parking area, then the vehicle owners can directly charge and park it.

Societal Factors (C 4 )
With the advancement in science and technology, the factors promoting wellness of human life have taken a back seat. The societal factors point out the major problems which may be looked at to improve the quality of the society at large.

•
Adverse impact of noise and electromagnetic fields (AI)(C 41 ): (Due to the construction of the electric vehicle charging station). An electric vehicle charging station has a constant aura of noise and an electromagnetic field surrounding it during the construction phase which might cause a certain category of people to develop problems. If proper measures can be taken in the initial stage, then this impact may be minimized since public health is of utmost importance.

•
Population density (C 42 ): This stands for the number of people living in each unit of area. When the population density in a locality is high, it shows that the area is overcrowded. There will be more consumption in such an area and the construction cost will also be high, but the quality of life will usually be low. The need of transportation in such an area is usually very high and an electric vehicle charging station constructed in such an area may just be what the people need.
Tables 1 and 2 describes the factors, sub-factors and alternatives considered for this research.    Figure 3 represents the hierarchical framework of the problem in the study. The rural, urban and total population for the census year starting from 2001 and 2011 has been represented in a graph (Figure 4). In this figure, the projected rural, urban and total population also been forecasted for the years 2021, 2031, 2041 and 2051 in a chronological way. The graph clearly indicates a consistent increase in the population from 4.23 million in 2001 to 6.60 million in 2051 (projected) with a sharp increase in the urban population, and the urban population growth is much higher than the rural, which indicates higher increasing infrastructural demand. The projected population made us think about the future transport services in the area.   The rural, urban and total population for the census year starting from 2001 and 2011 has been represented in a graph (Figure 4). In this figure, the projected rural, urban and total population also been forecasted for the years 2021, 2031, 2041 and 2051 in a chronological way. The graph clearly indicates a consistent increase in the population from 4.23 million in 2001 to 6.60 million in 2051 (projected) with a sharp increase in the  Urbanization is considered one of the most noteworthy anthropogenic inputs of the environmental framework, and thus the present study considered the spatio-temporal characteristics of urban growth and its inference in the transport of Howrah. The built-up land ( Figure 5) has been generated using the NDBI (normalized difference built-up index) with the following equation NDBI = MIR − (NIR/MIR) + NIR (Zha, 2003) [86]. Here, Urbanization is considered one of the most noteworthy anthropogenic inputs of the environmental framework, and thus the present study considered the spatio-temporal characteristics of urban growth and its inference in the transport of Howrah. The built-up land ( Figure 5) has been generated using the NDBI (normalized difference built-up index) with the following equation NDBI = MIR − (NIR/MIR) + NIR (Zha, 2003) [86]. Here, NIR is anear-infrared band such as ETM+ and TM, and LISSIII is a band no.4; MIR is a middle infrared band such as ETM+ and TM and LISSIII is a band no.5.The index is based on the unique spectral response of built-up lands that have a higher reflectance in the MIR wavelength range as compared to that in the NIR wave length range. Thereafter, the NDBI mapping has been prepared to understand the level of urbanization from 2000 ( Figure 5a) to 2010 (Figure 5b) in the study area. It helps to correlate the changes in land use patterns and its consequences to the water storage of the study area. The NDBI values range from The transport network of the Howrah district is mapped ( Figure 6) and the closest census towns (Figure 7) to the Howrah Municipal Corporation (HMC) were also plotted to understand the importance of daily communication with the Howrah station or surrounding areas. The railway and the road transport connectivity with the Howrah station, which is located in the HMC, is very good. A large number of daily commuters are traveling to the area, mostly for economic and educational purposes. It indicates the concentration of traffic in the area, and in turn increases the public transport connectivity as well as the requirement for improving the local transport system. The transport network of the Howrah district is mapped ( Figure 6) and the closest census towns (Figure 7) to the Howrah Municipal Corporation (HMC) were also plotted to understand the importance of daily communication with the Howrah station or surrounding areas. The railway and the road transport connectivity with the Howrah station, which is located in the HMC, is very good. A large number of daily commuters are traveling to the area, mostly for economic and educational purposes. It indicates the concentration of traffic in the area, and in turn increases the public transport connectivity as well as the requirement for improving the local transport system.  The ten selected points/locations ( Figure 8) were mapped in the HMC to understand the spatial coverage and important transport nodes in the area of the current study. The   The ten selected points/locations ( Figure 8) were mapped in the HMC to understand the spatial coverage and important transport nodes in the area of the current study. The The ten selected points/locations ( Figure 8) were mapped in the HMC to understand the spatial coverage and important transport nodes in the area of the current study. The population distribution of each of the 66 wards was also mapped (Figure 9) to observe the population pressure, which is also able to justify the present selection of the ten locations for the study. Most of the selected locations are in densely populated areas, where the public and local transport services are becoming very essential. Furthermore, the Howrah Maidan (S 4 ) and Salkia (S 8 ) are highly densely populated ( Figure 9) and thus very important in terms of transport services, whereas, Liluah (S 5 ), Belur (S 10 ) and Shibpur (S 7 ) belong to highly densely populated and Kadamtala (S 6 ), Bakultala (S 9 ), Belgachia (S 3 ) belong to moderately populated and Dasnagar (S 1 ) and Santragachi (S 2 ) belong to comparatively less populated among the ten points. But all these points are almost equally important and essential in terms of either population or transport or both. Sub criteria scores of each location has been represented in Table 3.
population distribution of each of the 66 wards was also mapped (Figure 9) to observe the population pressure, which is also able to justify the present selection of the ten locations for the study. Most of the selected locations are in densely populated areas, where the public and local transport services are becoming very essential. Furthermore, the Howrah Maidan ( 4 ) and Salkia ( 8 ) are highly densely populated ( Figure 9) and thus very important in terms of transport services, whereas, Liluah ( 5 ), Belur ( 10 ) and Shibpur ( 7 ) belong to highly densely populated and Kadamtala ( 6 ), Bakultala ( 9 ), Belgachia ( 3 ) belong to moderately populated and Dasnagar ( 1 ) and Santragachi ( 2 ) belong to comparatively less populated among the ten points. But all these points are almost equally important and essential in terms of either population or transport or both. Sub criteria scores of each location has been represented in Table 3.     Parking Areas In a scale of 1,3,5,7,9 Adverse impact of noise and electromagnetic field due to construction of electric vehicle charging station In a scale of 1, 3,5,7,9 Linguistic variables in HFN required for the comparison of factors and sub-factors are shown in Table 4. Comparison between factors in linguistic variables given by two DMs is presented in Table 5.
The above Table 6 represents preference of factors in defuzzified form using Equation (18).Normalized matrix is obtained using Equation (20); priority weight of factors are calculated using Equation (21). Societal factors obtain the maximum weight of 0.430, followed by environmental factor 0.37, followed by traffic factor 0.139 and economic factor 0.065.  The weight of factors obtained is represented in Table 7. Using Equation (22), C.I is calculated to be 0.08. For R.I, as n = 4, the value is 0.09. Thus using Equation (23), C.R = 0.08 0.9 = 0.09 < 0.1, hence the matrix is consistent. In the similar way, the comparison analysis of sub-factors with the help of two DMs has been calculated. The fuzzy weight of factors, sub-factors and global fuzzy weight are represented in Table 8.  The linguistic terms used in this study for the rating of alternatives with respect to sub-factors. Preference of alternatives with respect to sub-factors are expressed in linguistic terms, scores and crisp values are depicted in Table 9. Note 2. The sub-factors (C 22 ) and (C 23 ) are assigned score using Table 3; later, it is converted to HFN. A score of 9 implies "very high" and so on.

Ranking of Alternatives Using Fuzzy AHP-TOPSIS Method
Following the steps represented in Section 2.7, distance measure, relative closeness and the ranking of sites has been computed as depicted in Table 10. Table 10. Representation of distance measure, relative closeness and the ranking of sites.

Ranking of Alternatives UsingFuzzy AHP-COPRAS
Following the steps represented in Section 2.8, R +g , R −g , R g , R and ranking are computed as depicted in Table 11.

Comparison Analysis and Sensitivity Analysis
Two different MCDM techniques, fuzzy AHP-TOPSIS and fuzzy AHP-COPRAS were employed for the selection of the optimal site for aelectric vehicle charging station in and around the city of Howrah, West Bengal, India. Figure 10 denotes the comparative ranking obtained under the two methodologies, fuzzy AHP-TOPSIS and fuzzy AHP-COPRAS. We also tried to compare the alternatives ranking using the two different said proposed methods.  A sensitivity analysis was conducted to see the ranking obtained under different changing conditions. Figures 11 and 12 represents the clustered column chart to compare the ranking with the interchange of the sub-factors' weight. Two different cases are taken. In the first case, the sub-factors parking facilities 'C 15 ' and population density 'C 42 ' weights were interchanged. In the second case, land cost 'C 11 ' and generation of noise and air pollution 'C 21 ' weights were interchanged. For both these cases, two different methodologies were used in this study, i.e., fuzzy AHP-TOPSIS and fuzzy AHP-COPRAS. In the first case, the sub-factors parking facilities ′ 15 ′ and population density ′ 42 ′ weights were interchanged. In the second case, land cost ′ 11 ′ and generation of noise and air pollution′ 21 ′ weights were interchanged. For both these cases, two different methodologies were used in this study, i.e., fuzzy AHP-TOPSIS and fuzzy AHP-COPRAS. Figure 11, i.e., ranking obtained by fuzzy AHP-TOPSIS under sensitivity analysis shows that the alternative ( 4 ), ( 5 ), ( 6 ), ( 10 ), ( 3 ) and ( 2 ) are consistent with first, second, third, fourth, fifth and sixth position, respectively, under the considered two cases, whereas following Figure 12, i.e., ranking yield by fuzzy AHP-COPRAS under sensitivity analysis shows that the sites ( 4 ), ( 5 ), ( 2 ), ( 10 ), ( 3 ) and ( 9 ) scores the rank of first, second, fifth, seventh, ninth and tenth position, respectively. The other sites' variation of rank is noticed under the sensitivity analysis and depicted in the mentioned figures. Figure 11. Sensitivity analysis ranking obtained under fuzzy AHP-TOPSIS. Figure 11. Sensitivity analysis ranking obtained under fuzzy AHP-TOPSIS. Figure 11, i.e., ranking obtained by fuzzy AHP-TOPSIS under sensitivity analysis shows that the alternative (S 4 ), (S 5 ), (S 6 ), (S 10 ), (S 3 ) and (S 2 ) are consistent with first, second, third, fourth, fifth and sixth position, respectively, under the considered two cases, whereas following Figure 12, i.e., ranking yield by fuzzy AHP-COPRAS under sensitivity analysis shows that the sites (S 4 ), (S 5 ), (S 2 ), (S 10 ), (S 3 ) and (S 9 ) scores the rank of first, second, fifth, seventh, ninth and tenth position, respectively. The other sites' variation of rank is noticed under the sensitivity analysis and depicted in the mentioned figures.

Results and Discussion
This section discusses the results obtained by the methodology FAHP-TOPSIS, FAHP-COPRAS and the sensitivity analysis. The ranking obtained under the two MCDM techniques yield the site "Howrah Maidan" (S 4 ) as the best alternative for e-vehicle site selection followed by "Liluah" (S 5 ) and "Belur" (S 10 ). The FAHP-TOPSIS ranked the alternative "Belur" (S 10 ) and "Salkia" (S 8 ) equally at the third position. Ranking obtained for all the sites are presented in Tables 10 and 11. In the sensitivity analysis, where the sub-factors' weight are interchanged as discussed in Section 5, it is seen that the site selected "Howrah Maidan" (S 4 ) consistently remains in the first position. The rankings obtained under the sensitivity analysis by the two methods are depicted in Figures 11 and 12.

Conclusionsand Future Scope
Ease of commutation, a pollution-free mode of transport, as well as employment generation are direct and indirect benefits of the e-vehicle. For developing countries where the pollution level is quite high and proportion of the roads is lower, e-vehicles can be a game changer. According to our study, considering ten locations across the city of Howrah Maidan due to its proximity to India's largest railway station Howrah and various other attributes ranked it as number one, followed by Liluah and Belur. These three locations are consistent with rank one, two and three, respectively, irrespective of the two MCDM methodologies applied in this research. High population density, enhanced level of consumerism, and the presence of various public facilities with higher footfalls led to this higher ranking.
This paper used the GIS and MCDM tools FAHP, F-TOPSIS and F-COPRAS to obtain the optimal site selection for the e-vehicle charging station. HFN has been used to give a preferential rating of factors, sub-factors and alternatives. The ranking obtained using MCDM tools are logical and scientific. The present findings provide important references for future potential work and problem solving.
For the site selection, factors such as environmental, economic, traffic, societal are incorporated with their respective sub-factors. Through this research, it is observed that the societal factor is the most significant and out of the sub-factors, population density is the most important.
HFN is used here as it captures the hesitancy and vagueness in an efficient way. To practice the qualitative criteria evaluation for imprecise information, FAHP, F-TOPSIS and F-COPRAS are used. Comparative analysis which uses F-TOPSIS and F-COPRAS in our example has depicted consistent results. The reliability, robustness and efficiency of this methodology has also been tested through sensitivity analysis. HFN captures a wider range of linguistic terms but usage of HFN makes computation a bit longer.
In Howrah, India, a very old and unplanned city, creating e-vehicle charging infrastructure across the city is important. A site like Santragachi (S 2 ) in this research which ranks lower at present could acquire a significant position due to huge infrastructural investment by the government in that area. A big railway terminus is coming up, which has the potential to change the demography of the area. Thus futuristic sub-factors can be incorporated into future research. The absence of futuristic sub-factors is a limitation at present. A larger number of decision makers based on the administrative point of view can be explored in future research. The proposed methodology used in this research can be applied in different fields such as new vendor selection, and treatment selection for new diseases where more ambiguity and uncertainty is prevalent. The other MCDM tools such as PROMETHEE, VIKOR, and WASPAS can be used in the future with intuitionistic, neutrosophic, and hesitant fuzzy numbers to yield improved, robust and practical solutions.