Optimum Design of a Standalone Solar Photovoltaic System Based on Novel Integration of Iterative-PESA-II and AHP-VIKOR Methods

Solar energy is considered one of the most important renewable energy resources, and can be used to power a stand-alone photovoltaic (SAPV) system for supplying electricity in a remote area. However, inconstancy and unpredictable amounts of solar radiation are considered major obstacles in designing SAPV systems. Therefore, an accurate sizing method is necessary to apply in order to find an optimal configuration and fulfil the required load demand. In this study, a novel hybrid sizing approach was developed on the basis of techno-economic objectives to optimally size the SAPV system. The proposed hybrid method consisted of an intuitive method to estimate initial numbers of PV modules and storage battery, an iterative approach to accurately generate a set of wide ranges of optimal configurations, and a Pareto envelope-based selection algorithm (PESA-II) to reduce large configuration by efficacy obtaining a set of Pareto front (PF) solutions. Subsequently, the optimal configurations were ranked by using an integrated analytic hierarchy process (AHP) and vlsekeriterijumskaoptimizacija i kompromisonoresenje (VIKOR). The techno-economic objectives were loss of load probability, life cycle cost, and levelized cost of energy. The performance analysis results demonstrated that the lead?acid battery was reliable and more cost-effective than the other types of storage battery. Record Type: Published Article Submitted To: LAPSE (Living Archive for Process Systems Engineering) Citation (overall record, always the latest version): LAPSE:2020.0510 Citation (this specific file, latest version): LAPSE:2020.0510-1 Citation (this specific file, this version): LAPSE:2020.0510-1v1 DOI of Published Version: https://doi.org/10.3390/pr8030367 License: Creative Commons Attribution 4.0 International (CC BY 4.0) Powered by TCPDF (www.tcpdf.org)


Introduction
The electricity demand is rapidly increasing due to growth in population and the risk of increasing electricity bills and tariffs. These issues lead to encourage energy system designers to a transformation of technology in terms of "leaving the grid" or "living in off-grid" [1,2]. The standalone photovoltaic (SAPV) system is considered one of the most important renewable, sustainable, clean, and environmentally friendly energy sources. However, SAPV systems need to be optimally designed in order to maximize their reliability and minimize the total cost of the system [3][4][5][6]. Therefore, an efficient sizing of a SAPV system is essential in order to supply the demanded electricity. The optimal design of the SAPV system strongly depends on the availability of the meteorological and hourly load demand.

Optimization Methods
Indeed, the intuitive method utilizes simplified calculations without considering the fluctuation of meteorological data and the nonlinearity relationship between subsystems [10]. In addition to this, the output-generated power of the PV arrays is calculated during the designing period that is related to safety factor, which the safety factor is chosen by the user [27,28]. Ghafoor et al. [29] proposed a sizing method to optimal the size of an off-grid PV system in Pakistan. A simple mathematical method was used to estimate the output power of PV arrays. Moreover, the average daily meteorological data were used to minimize the annualized life cycle cost (ALCC), and the technical criteria was not revealed.
To overcome the limitations of the intuitive methods, analytical methods are used for this purpose. In an analytical method, a relationship is proposed as being compatible with the reliability and size of the SAPV system. The performance of the SAPV system can be evaluated for a various set of the feasible size of system components. Then, the optimal configuration is selected by comparing single or multiple performance indices of various sets of configurations. The sizing of the PV system can be characterized by a simple calculation. Meanwhile, the major drawback is in finding the coefficients of the equation of the standalone PV system [30][31][32][33]. Riza in 2015 [34] presented a sizing methodology of PV panel and battery capacity in the SAPV system with a lighting load in Malaysia. The design space approach was obtained for sizing the curve of the PV system and storage battery configuration using three objectives: LPSP, LCC, and excess energy. According to the authors, the proposed method in [34] offered an advantage because the user can select the realistic available PV array and storage battery size configuration. However, simple component models for the system model were used, which affected the result of the output because these models cannot give the variation of the realistic metrological data. Thus, the system may be over/under-sizing, leading to increase in the capital costs or decrease in the reliability of the entire system.
In general, the most used method for sizing of the SAPV system is the numerical method due to it is accuracy and ability to discover all the design space [10]. This method takes a long period of time in implementing the system simulation. The numerical method can be classified into two parts, which are stochastic (deterministic) and probabilistic approaches [10,35]. The probabilistic approach is preferred because the uncertainty of the meteorological data is considered, in which the performance of the system can be more reliable and accurate than the stochastic approach. Moreover, the energy availability of the system can be represented in a quantitative manner [15,[36][37][38]. In the numerical method, the various number of the PV arrays and storage battery configurations are conducted in the design space. Then, each pair of the PV arrays and storage batteries are simulated on the basis of hourly meteorological data and the objective function to estimate the desired configuration at a defined level of reliability of the system. Then, the lower cost from the set of configurations is selected on the basis of the objective function as an optimal solution [39]. Daily meteorological data and simple PV modules were used in these research works [13][14][15][16][17][18][19]40] to find the current-voltage (I-V) data curve at different operation conditions. The daily meteorological data were used in [13][14][15][16][17][18]40] lead to give unsatisfactory performance of the system.
Meanwhile, other studies used heuristic techniques in order to find the optimal size of a SAPV system using techno-economic objective functions such as artificial bee colony (ABC) [41], genetic algorithm (GA) [42,43], generalized regression neural network (GRNN) [44], firefly (FL) [45], and particle swarm optimization (PSO) [46]. The main advantage of these searching algorithms is the ability to converge the optimal solution in a short time. Thus, heuristic techniques can overcome the drawback of numerical algorithms. In addition to that, some of the studies combined two or more in order to enhance the converge to the optimal values and reduce the execution time, such as the research works presented in [43,47,48]. Hussein et al. [46] proposed an optimal design of the SAPV system using multi-objective particle swarm optimization (MOPSO) in Malaysia. Two variants of the PSO algorithm were presented with reference to a sigmoid function PSO (SFPSO c f ) and adaptive weights PSO (AWPSO c f ) using techno-economic criteria. The performance results demonstrated that the SFPSO c f algorithm has a trivial superiority in selecting an optimal configuration in terms of accuracy. However, the main demerit of the method proposed in reference [46] the use of a non-scale (NS) approach that can conduct only one solution on the basis of given weights. In addition, the proposed algorithm cannot cover all the Pareto front (PF) solution when all possible status of weights within [0,1] are considered.
On the basis of the literature review, the proposed numerical methods in [13,15,16,18,[36][37][38][39][40]49] do not simultaneously solve multi-objective optimization for the sizing of the SAPV system, which is considered as a theoretical gap in the presented study. Therefore, the real challenge is to obtain a set of optimal solutions called the Pareto optimal solutions front. However, the challenge is in providing solutions close to the true PF (convergence), finding well-distributed solutions on the PF (diversity), and covering maximum PF (coverage). In [43], a set of optimal configurations was generated using the NSGA-II algorithm for the SAPV system based on two types of storage batteries: lead-acid and lithium-ion. However, NSGA-II may fail in establishing an accurate PF solution due to the obtained aggregation method [24,50]. Moreover, the authors in [50] proposed three scenarios in order to optimally determine the size of the SAPV system on the basis of a salp swarm algorithm considering techno-economic criteria. The results demonstrated that the multi-objective salp swarm non-dominated roulette wheel (MOSS-NDRW) scenario outperformed others in terms of converge and convergence.
As a matter of fact, the aim of establishing a set of PF solutions is to choose the final "best solution" alternative. Therefore, MCDM techniques have been applied in order to provide a preferred structure according to the relative importance of criteria [51]. To rank alternatives, the technique for order preference by similarity to ideal solution (TOPSIS) and vlsekeriterijumskaoptimizacija i kompromisonoresenje (VIKOR) are suitable in order to determine the best feasible solution for the decision-maker [52]. Both methods are efficient in order to use when the quantity data and objective are given. In TOPSIS method, the alternatives can be ordered based on the shortest and farthest distances from the ideal to the negative ideal solutions, respectively. From this view, in [53], the authors proposed a population-based differential evolution (DE) method in order to construct a set of optimal solutions based on NSGA-II, then AHP integrated with TOPSIS was used to select the optimal solution from a set of PF solutions. However, the relative importance is not considered by TOPSIS [54]. On the other hand, VIKOR is functionally associated with discrete alternative problems, making it a more practical technique for tackling real-world problems [55]. Moreover, VIKOR can rapidly select the best alternative from many alternatives and attributes. However, the main obstacle of the VIKOR approach is in the selection of appropriate weight and checking judgment consistency. Thus, AHP is a convenient technique for acquiring the relative importance of various criteria.
In this study, a novel integration of the Iterative-PESA-II and AHP-VIKOR method was developed in order to select the optimum size of the SAPV system. Firstly, the intuitive method was obtained for rough calculation of numbers of PV modules and storage batteries. Then, the iterative approach computed the number of PV modules that were connected in series and parallel. Moreover, the size of the SAPV system was optimally considered on the basis of the techno-economic criteria of LLP, LCC, and levelized cost of energy (LCE). In this proposed sizing approach, the favorable level of reliability was defined in order to find a set of SAPV system configurations using hourly meteorological data for one year. After that, the Pareto envelope-based selection algorithm (PESA-II) method was utilized in order to nominate a set of optimal configurations to construct optimal PF solutions. Subsequently, the integrated AHP-VIKOR technique was proposed in order to rank and choose the optimum configuration of the SAPV system. The proposed hybrid approach overcame all the drawbacks of intuitive, analytical, numerical, and stochastic methods that were faced in the sizing of the SAPV system. Finally, an accurate single-diode (SD) PV model was used to establish the I-V characteristic curve, which can give a realistic evaluation for the extracted energy of the PV array. The proposed method showed a high level of reliability and low cost in designing a SAPV system.

Standalone PV System
A typical SAPV system mainly comprises the PV arrays that convert the sun's rays to a DC current; the storage battery can be used for storing the excess energy and power conditioner units containing a DC-DC converter and a DC-AC inverter. Figure 1 depicts a typical standalone PV system. characteristic curve, which can give a realistic evaluation for the extracted energy of the PV array. The proposed method showed a high level of reliability and low cost in designing a SAPV system.

Standalone PV System
A typical SAPV system mainly comprises the PV arrays that convert the sun's rays to a DC current; the storage battery can be used for storing the excess energy and power conditioner units containing a DC-DC converter and a DC-AC inverter. Figure 1 depicts a typical standalone PV system. Figure 1. The block diagram of a generalized stand-alone photovoltaic (SAPV) system.

Single-Diode Model
The most commonly used representation of the performance of the output current of the PV module is the single-diode model. This is due to it is simplicity, accuracy, and ability to reflect the actual behavior of the PV cells, as illustrated in Figure 2. For instance, the single-diode PV model offers a high ability in reflecting reality behavior when it is dealing with non-linearity and stochastic nature among other mathematical models [56]. The single-diode model consists of current and voltage as the outputs current and voltage of the PV module; photocurrent ℎ , which has a high sensitivity to environmental conditions; the diode reverse saturation (in A); and the parallel and series resistances and , which represent the losses of the PV solar cell (Ω). Thus, the following equation represents the output current of the single-diode PV module: where refers to the diode thermal voltage (in V) and can be represented as where is the ideality factor of the diode, refers to the Boltzmann constant (1.3806503 × 10 −19 J/K), is the cell temperature ( ), and represents the electron charge (1.60217646 × 10 −19 Coulombs). The five parameters ( ℎ , , , , and ) need to be extracted accurately.

Single-Diode Model
The most commonly used representation of the performance of the output current of the PV module is the single-diode model. This is due to it is simplicity, accuracy, and ability to reflect the actual behavior of the PV cells, as illustrated in Figure 2. For instance, the single-diode PV model offers a high ability in reflecting reality behavior when it is dealing with non-linearity and stochastic nature among other mathematical models [56]. The single-diode model consists of current I and voltage V as the outputs current and voltage of the PV module; photocurrent I Ph , which has a high sensitivity to environmental conditions; the diode reverse saturation I o (in A); and the parallel and series resistances R p and R s , which represent the losses of the PV solar cell (Ω). Thus, the following equation represents the output current of the single-diode PV module: where V t refers to the diode thermal voltage (in V) and can be represented as where d is the ideality factor of the diode, KB refers to the Boltzmann constant (1.3806503 × 10 −19 J/K), T c is the cell temperature (K), and q represents the electron charge (1.60217646 × 10 −19 Coulombs). The five parameters (I Ph , I o , R s , R p , and d) need to be extracted accurately.  Furthermore, the PV array consists of parallel and series modules. Therefore, the output current of the PV array is expressed as follows: where and are the output current (A) and voltage (V) of the PV array, respectively. The five parameters ( ℎ , , , , and ) are extracted by the improved electromagnetism-like (IEM) algorithm [57]. The accuracy of the IEM algorithm can be represented by establishing current-voltage (I-V) and power-voltage (P-V) characteristic curves, which are shown in Figure 3a,b. The efficiency of the PV array is computed by the following equation: where denotes to the PV array's area (m 2 ), and represents the amount of the solar irradiation fall on the surface of the PV array (W/m 2 ). The type of the used PV module is a Kyocera KC120-1 multicrystalline silicon, which was used in this research study. To simulate the output power of the PV array, a maximum power point tracker (MPPT) was employed in this study as shown in Table 1.  Furthermore, the PV array consists of parallel N p and series N s modules. Therefore, the output current of the PV array is expressed as follows: where I and V are the output current (A) and voltage (V) of the PV array, respectively. The five parameters (I Ph , I o , R s , R p , and d) are extracted by the improved electromagnetism-like (IEM) algorithm [57]. The accuracy of the IEM algorithm can be represented by establishing current-voltage (I-V) and power-voltage (P-V) characteristic curves, which are shown in Figure 3a,b. The efficiency of the PV array is computed by the following equation: where A denotes to the PV array's area m 2 , and G represents the amount of the solar irradiation fall on the surface of the PV array W/m 2 . The type of the used PV module is a Kyocera KC120-1 multi-crystalline silicon, which was used in this research study. To simulate the output power of the PV array, a maximum power point tracker (MPPT) was employed in this study as shown in Table 1.

Battery Mathematical Model
Various types of storage battery capacity (kW) were proposed in this study in order to verify those most effective and suitable for the SAPV system. Therefore, the days of autonomy are given by the following [58]: where represents the capacity of the storage battery; refers to number of the autonomy days (typically between 3-5 days) [13]; represents the required load demand; and * *

Battery Mathematical Model
Various types of storage battery capacity (kW) were proposed in this study in order to verify those most effective and suitable for the SAPV system. Therefore, the days of autonomy are given by the following [58]: where C bat represents the capacity of the storage battery; AD refers to number of the autonomy days (typically between 3-5 days) [13]; E Load represents the required load demand; and DOD * ηb * ηinv are the depth of discharge (80%), efficiency of the storage battery (85%), and efficiency of the inverter (95%). The state of chare of the battery (SOC) can be expressed by [36] where P pv (t) and P load (t) refer to the hourly output energy generated by PV modules and energy required by the load, SOC(t − 1) and SOC(t) represent the status of the state of charge (SOC) for the battery at initial and final points of charging and discharging, respectively; and E bat denotes the battery hourly capacity. For a more reliable and longer lifetime for the storage battery, the minimum energy stored in the battery must be pointed as given by the following: The required energy supposed to be stored in the battery can be written by the following equation [13]: where ηch represents the battery's charging efficiency.

Load Profile
In this work, the hourly load demand of a typical house in a rural area in Malaysia was considered, which can be found in details in the load demand in [46]. The proposed hybrid method was analyzed using hourly meteorological data for one year. Furthermore, the data used was taken from Subung meteorological station in Subung Jaya, Klang Valley, Malaysia (latitude 3.12 • north and longitude 101.6 • east) [46].

Proposed Criteria for the Sizing of the SAPV System
In order to obtain the optimal configuration of the SAPV system, the number of PV modules (in series (N s ) and parallel (N s )) and the number of the storage battery (Bat) should be chosen accurately by considering techno-economic criteria. Therefore, three important criteria, namely, loss of load probability (LLP), life cycle cost (LCC), and levelized cost of energy (LCE), were chosen to increase the reliability and minimize the total cost of the system. In this study, the objective function (F(X)) of the problem could be obtained by simultaneously minimizing LLP, LCC, and LCE. Thus, the multi-objective optimization (MOOP) to optimally size of the SAPV system can be described by the following equations [59]: where X is the vector of the decision variables that should be optimized by the proposed method.

Loss of Load Probability
The availability of the system can be measured by using LLP, which is defined as the ratio of the energy deficits to the required load demand during a specific period of time [60]. The range of the LLP's values are between 0 and 1. Thus, high reliability could be chosen for optimal sizing of the SAPV system when it was close to 0. LLP can be represented by the following equation:

Life Cycle Cost
The best indicator to evaluate the capital cost of the system with it is components was the life cycle cost (LCC). Thus, LCC was used to compute four parts of the SAPV system which were PV panels, converter, storage battery, and inverter. Moreover, LLC considered the initial total cost values (IC tot ), replacement cost values C rep , and operation, as well as maintenance cost values (C O&M ) during a period of time. The unit cost and the lifetime of the system's components are given in Table 2. Therefore, the mathematical equations can be written by [41]: The initial total cost of all components of the SAPV system was considered, which included the price of each part of the system, the installation, and the amount of the civil works. The initial total cost can be given by (IC tot ), which is given by [61]: where C PV, C Unit,PV are total capacity and unit cost of the PV array, respectively; C Batt , C Unit,Batt are the total capacity and unit cost of the battery, respectively; C Bidi , C Unit,Bidi are the total capacity and unit cost of the inverter, respectively; C CH , C Unit,CH are the total capacity and unit cost of the converter, respectively; and (C O ) represents the total constant cost with including both civil work and installation cost. In this study, three components required replacement during the system lifetime, these being the storage battery, converter, and inverter. The values of the inflation rate (FR) and interest rate (IR) of the components of the system were considered in this study. Thus, the replacement cost (C rep ) can be given as [61] C rep = C Unit × C nom where C Unit represents the unit component cost of the storage battery, converter, and inverter; (C nom ) denotes to the nominal capacity of the replacement system component; and (N rep ) refers to the number of replacement of each component over the system life period (LP). The operation and maintenance cost (C O&M ) calculation can be described as [62] where C O&M0 is the operation and maintenance cost in the first year. The detailed units' costs and services of the proposed SAPV system can be found in [61,63], whereas the units' cost of the lead-acid, absorbent glass mat (AGM), and lithium-ion batteries are tabulated in Table 2.

Levelized Cost of Energy (LCE)
Another parameter used to represent the economic aspects of the SAPV system was LCE. LCE is defined as the ratio of the total cost of the all components of the SAPV system to the total of the conducted energy by PV arrays throughout one year, which can be given as follows [64,65]: where TAC is the total cost of the system components and E tot is the total annual energy generated by the system during a specific period of time and can be represented by the following equation:

The Proposed Optimization Sizing Method
In this research work, the proposed optimization of the SAPV system was based on four combined methods-the intuitive and numerical, PESA-II, and integrated AHP-VIKOR methods. The descriptions of these methods are shown in Figure 4. Another parameter used to represent the economic aspects of the SAPV system was LCE. LCE is defined as the ratio of the total cost of the all components of the SAPV system to the total of the conducted energy by PV arrays throughout one year, which can be given as follows [64,65]: where is the total cost of the system components and is the total annual energy generated by the system during a specific period of time and can be represented by the following equation:

The Proposed Optimization Sizing Method
In this research work, the proposed optimization of the SAPV system was based on four combined methods-the intuitive and numerical, PESA-II, and integrated AHP-VIKOR methods. The descriptions of these methods are shown in Figure 4.

Intuitive Method
The intuitive method is used to predict the number of PV modules and storage batteries. However, it was observed from the literature that the intuitive method may result in over/under-sizing of the SAPV system. Therefore, the intuitive method was applied in order to determine the ranges of search space. The number of PV modules can be calculated on the basis of the intuitive method, as follows: where E load represents the daily energy load demand, PSH is the peak sunshine hours, η inv and η bat are the efficiencies of the inverter and system, A is the area of the PV array, and S F is the design safety factor. The initial capacity value of the storage battery is calculated on the basis of Equation (5). In addition to this, the number of the storage batteries in series and parallel can be calculated by the following equations: where V bus and V battery are the voltage bus of the system and battery, respectively, and C b,Wh represents the days of autonomy and W-hour capacity of the battery. Thus, the total number of the storage battery is given by multiplying the series N Sbat and parallel N Pbat numbers of batteries, and is expressed as follows: The N Tbat is used to determine the design space of the storage battery. Afterward, the numerical method is used for finding a set of optimal configurations of the system.

Iterative Method
After determining the bounders of the PV modules and storage batteries numbers, the iterative method is employed by increasing the numbers of PV modules (in parallel and series) and storage battery from minimum to maximum based on predefined ranges. The operation process of energy management of SAPV on the basis of the iterative method is illustrated in Figure 5. The first step of the algorithm can be performed by obtaining the specification of the system such as the efficiency of the PV module, efficiency of charging, and discharging of the storage battery, and the efficiency of the inverter, converter, and wires. After obtaining the required specifications, which are tabulated in Table 3, the simulation starts to calculate the E_net by subtracting the predicted output energy of the PV model from the energy of the load. Then, the net energy is classified in three cases, which are In the first case, the battery energy is not used and there is no damped or deficit energy that is equal to zero. In the second case, the energy produced by the PV array is larger than the load energy. There is an excess energy and depending on the state of charge (SOC) of the battery. If the SOC of the battery is full, the amount of excess energy will be damped. In contrast, the amount of excess energy will be used in order to charge the battery and the new value of SOC will be calculated by Equation (7) and the deficit energy is zero. In case 3, if the SOC of the storage battery is more than the minimum value E bat_min (Equation (8)) then the storage battery is able to meet the load demand. The load will be supplied with/without the energy of the PV array and storage battery. Or else, if the SOC of the battery is less than the minimum value, the system cannot meet the load demand and there is deficit energy. The three-objective function of LLP, LCC, and LCE are calculated in each iteration. battery is less than the minimum value, the system cannot meet the load demand and there is deficit energy. The three-objective function of LLP, LCC, and LCE are calculated in each iteration.  Figure 5. Flow chart of the proposed sizing numerical method of the SAPV system. Figure 5. Flow chart of the proposed sizing numerical method of the SAPV system.

PESA-II Method
Indeed, the numerical methods are most accurate due to their ability to discover all possible configurations in the design space. However, the iterative methods suffer not only from prolonged execution time but also in selecting the desired optimal configuration (global solution) due to the generated large number of configurations. The real challenge is to select a set of optimal solutions called Pareto front solutions. Therefore, the PESA-II method is abstained to handle the three objectives LLP, LCC, and LCE simultaneously. The Pareto envelope-based selection algorithm (PESA) was proposed by Corne et al. [66]. In PESA, there are two types of population: internal (IP) and external population (EP), the EP is archived to store the non-dominated solution in which a set of PF solutions can be established. In our proposed method, PESA-II concept can be explained by the following steps: Step 1: Reading the saved numbers of PV modules and storage batteries, then store them in the internal population (IP).
Step 2: Reading the LLP, LCC, LCE objectives, then saving them in the objective function (OF) array.
Step 3: Comprise the non-dominated members of IP into EP.
Step 4: The next member of the EP (achieve) being chosen on the basis of region-based selection by assigning the three objectives' fitness in objective space to all the hyperboxes.
Step 5: The fewer individuals in the hyperbox being chosen as a winner for the next iteration. In this way, PESA-II can provide better convergences and diversity in a solution space than SPEA, PAES, and PESA algorithms.
In step 3, the EP always updates the non-dominated solutions and its size is set as 200. The new solution enters the EP if it is non-dominated within IP and should be not dominated by current solutions in the EP because the dominated solution will be removed from EP. According to steps 4 and 5, the selective fitness method is applied to all hyperboxes in order to choose fewer individuals on the hyperbox using binary tournament selection. By using binary tournament selection, the chance of selecting most isolated individuals (fewer individuals) will be 1 − P−l P 2 , assuming there are two hyperboxes in the current approximation to the PF, as illustrated in Figure 6. This means that increasing the number of individuals leads to having a high possibility of selection. Meanwhile, the chance of selecting the best individual is m P 2 . Therefore, the chance of selecting an individual from the least crowded box is more than crowded individual which can be calculated regarding to the ratio of this probability 2Pl − l 2 /m 2 , where l and m are the least and most crowded boxes, respectively. In cases where the new non-dominated solutions try to enter EP when EP is over-full, the new non-dominated solutions can occupy less crowded members than the other solutions. The main advantage of PESA-II can be pointed out by using the binary tournament selection to adaptive range equalization and normalization of the OF values. Thus, this selection technique is oriented towards obtaining PF in a well-distributed manner [67]. In this method, there is no mutation or crossover phases; this is because all possible configurations have been investigated previously and there is no need to search for new regions in the search space. Therefore, the execution is performed only one time.

Integrated AHP-VIKOR Method
The main goal of proposing an integrated AHP-VIKOR method is to select one optimal solution (optimum solution) from the overall PF. Therefore, the AHP method was used in this study to calculate the weights for each criterion in a decision matrix (DM). Afterward, the VIKOR method was employed to a rank optimal configuration on the basis of preference criteria considered for practical judgment. Figure 7 demonstrates the structure of the integrated AHP-VIKOR method. The next subsection will discuss the integrated AHP-VIKOR method in detail.

Integrated AHP-VIKOR Method
The main goal of proposing an integrated AHP-VIKOR method is to select one optimal solution (optimum solution) from the overall PF. Therefore, the AHP method was used in this study to calculate the weights for each criterion in a decision matrix (DM). Afterward, the VIKOR method was employed to a rank optimal configuration on the basis of preference criteria considered for practical judgment. Figure 7 demonstrates the structure of the integrated AHP-VIKOR method. The next subsection will discuss the integrated AHP-VIKOR method in detail.

Integrated AHP-VIKOR Method
The main goal of proposing an integrated AHP-VIKOR method is to select one optimal solution (optimum solution) from the overall PF. Therefore, the AHP method was used in this study to calculate the weights for each criterion in a decision matrix (DM). Afterward, the VIKOR method was employed to a rank optimal configuration on the basis of preference criteria considered for practical judgment. Figure 7 demonstrates the structure of the integrated AHP-VIKOR method. The next subsection will discuss the integrated AHP-VIKOR method in detail.

AHP Method
The AHP method was developed by Saaty et al. [68], and is most widely used to demonstrate the priority importance of the criteria. The proper weights of criteria can be determined by using a pairwise comparison matrix on the basis of the opinion of three experts. The steps of the integrated AHP-VIKOR method are given by the following: Step 1 The AHP method can be started after defining a hierarchy. The hierarchy includes the criteria and decision goal. The criteria used to select the optimal configuration of SAPV system are LLP, LCC, and LCE. The number of pairwise comparisons can be computed by n × (n − 1)/2, where n is the number of criteria. The pairwise comparison matrix is tabulated in Table 4 to show the importance of criteria on the basis of three evaluators. In this study, a pairwise comparison matrix is constructed on the basis of three experts. To initialize the weights, a pairwise comparison matrix can be built in the next step.
Step 2 The preferences on the three criteria obtained in AHP are considered. Normalization of each element is required in order to construct the normalized matrix [69], where each element in the (prefer order) matrix is the corresponding each element of the same matrix divided by the sum of all columns. Afterward, the aggregation vector's elements are filled by the sum of elements of the corresponding row of the normalized matrix.
Step 3 The weights can be computed by dividing each element in the aggregation vector by taking the sum of all elements in the aggregation vector. The summation of weights for three criteria of each evaluator must not exceed 1.
Step 4 In this step, the consistency ratio (CR) is used to check if the elements of the preferred order matrix are consistent or not. The equation of CR can be written as follows: where RI is a random index, which was chosen as 0.58 for this study, and CI is the consistency index, which is computed by the following equation: where λmax is the maximum value of the preferred order matrix and n is the order of the matrix.
It is important to notice that, if the value of the CR is larger than 0.1, the pairwise comparison matrix is not consistent, and it will be ignored. The pairwise comparison matrix, normalized matrix, aggregation vector, weights vector, and consistency ratio are illustrated in Table 5.

VIKOR Method
After assigning the preferences of the three distinct objectives by three evaluators using the AHP technique, the need for ranking the optimal configuration of the SAPV system can be performed by using the VIKOR method. This is because the VIKOR method has the ability to identify the best alternatives in situations with multiple criteria [55]. The mathematical steps of the VIKOR methods are given in Figure 8.

VIKOR Method
After assigning the preferences of the three distinct objectives by three evaluators using the AHP technique, the need for ranking the optimal configuration of the SAPV system can be performed by using the VIKOR method. This is because the VIKOR method has the ability to identify the best alternatives in situations with multiple criteria [55]. The mathematical steps of the VIKOR methods are given in Figure 8.

Results and Discussion
In this study, the proposed methodology was obtained in order to find a set of optimal configurations of the SAPV system to supply the load demand for a house in a remote area. Firstly, the intuitive method was employed to determine the search space of PV modules (series and parallel) and the three types of storage batteries (lead-acid, AGM, and lithium-ion). The ranges of the search space, which consisted of the number of PV modules and number of storage batteries used in the iterative method, were chosen by utilizing the intuitive method from 1 to 40 for series and parallel PV modules, whereas the range of the number of the storage battery was chosen from 3 to 50. The

Results and Discussion
In this study, the proposed methodology was obtained in order to find a set of optimal configurations of the SAPV system to supply the load demand for a house in a remote area. Firstly, the intuitive method was employed to determine the search space of PV modules (series and parallel) and the three types of storage batteries (lead-acid, AGM, and lithium-ion). The ranges of the search space, which consisted of the number of PV modules and number of storage batteries used in the iterative method, were chosen by utilizing the intuitive method from 1 to 40 for series and parallel PV modules, whereas the range of the number of the storage battery was chosen from 3 to 50. The safety factor was set as 2 in the intuitive method. Afterward, the iterative method was used to construct sets of optimum configurations using a wide range of reliability and costs for the three types of storage batteries. The main advantage of the numerical approach is that the global optimal solutions can be reached by investigating all possible configurations in the design space. The numerical approach started by increasing N s and N p one by one in the first loop and the second loop, respectively. Meanwhile, the number of the storage battery increased in the third loop. Therefore, the numbers of configurations based on the numerical method were 59,503, 45,806, and 40,098 and the central processing unit (CPU) executions times were 10,295.71 s, 8151.07 s, and 8656.25 s for the lead-acid, AGM, and lithium-ion batteries, respectively. Figure 9 shows the correlation between the LLP, LCC, and LCE starting from the minimum value of each objective and ending at the maximum value by using the numerical approach for lead-acid, AGM, and lithium-ion storage batteries. According to [13], the acceptable level of LLP is equal to 1%. Thus, the favorable level of LLP was set at ≥0.1%. For convenience, the maximum and minimum values of LLP, LCC, and LCE objectives with their numbers of PV modules and storage batteries are tabulated in Table 6. safety factor was set as 2 in the intuitive method. Afterward, the iterative method was used to construct sets of optimum configurations using a wide range of reliability and costs for the three types of storage batteries. The main advantage of the numerical approach is that the global optimal solutions can be reached by investigating all possible configurations in the design space. The numerical approach started by increasing and one by one in the first loop and the second loop, respectively. Meanwhile, the number of the storage battery increased in the third loop. Therefore, the numbers of configurations based on the numerical method were 59,503, 45,806, and 40,098 and the central processing unit (CPU) executions times were 10,295.71 s, 8151.07 s, and 8656.25 s for the leadacid, AGM, and lithium-ion batteries, respectively. Figure 9 shows the correlation between the LLP, LCC, and LCE starting from the minimum value of each objective and ending at the maximum value by using the numerical approach for lead-acid, AGM, and lithium-ion storage batteries. According to [13], the acceptable level of LLP is equal to 1%. Thus, the favorable level of LLP was set at ≥0.1%. For convenience, the maximum and minimum values of LLP, LCC, and LCE objectives with their numbers of PV modules and storage batteries are tabulated in Table 6.   It is worth mentioning that selecting the optimal configuration from large numbers of configurations using the numerical method is a complex task. It is necessary to nominate a set of optimal solutions to be represented as optimal PF of multi-objectives by simultaneously minimizing the LLP, LCC, and LCE. On this basis, the PESA-II method was employed to perform this task. As mentioned in the literature, the difficulty of several multi-objective algorithms is in comparing two conflicting solutions and finding the fitness of the most probable solution. Moreover, the obtained PF must achieve three objectives, which are convergence, diversity, and coverage. In this research study, PSEA-II is used to obtain a set of optimal PF solutions based on three conflicting objectives, which are LLP, LCC, and LCE. The correlation between LLP, LCC, and LCE objectives using three types of batteries is illustrated in Figure 10. The aim of obtaining sets of PF solutions is to minimize the three objectives, and the trade-off is shown in Figure 10. From Figure 11, the PESA-II method can efficiently obtain sets of PF solutions for lead-acid, AGM, and lithium-ion batteries in terms of good spread, convergence, and coverage. The CPU execution time of the three sets PF solutions for the lead-acid, AGM, and lithium-ion batteries were 19,810.81 s, 4783.03 s, and 1347.06 s, respectively. For more convenience, the minimum and maximum values of the three objectives for the three types of batteries using the PESA-II method are tabulated in Table 7.
Processes 2020, 8, x; doi: FOR PEER REVIEW www.mdpi.com/journal/processes It is worth mentioning that selecting the optimal configuration from large numbers of configurations using the numerical method is a complex task. It is necessary to nominate a set of optimal solutions to be represented as optimal PF of multi-objectives by simultaneously minimizing the LLP, LCC, and LCE. On this basis, the PESA-II method was employed to perform this task. As mentioned in the literature, the difficulty of several multi-objective algorithms is in comparing two conflicting solutions and finding the fitness of the most probable solution. Moreover, the obtained PF must achieve three objectives, which are convergence, diversity, and coverage. In this research study, PSEA-II is used to obtain a set of optimal PF solutions based on three conflicting objectives, which are LLP, LCC, and LCE. The correlation between LLP, LCC, and LCE objectives using three types of batteries is illustrated in Figure 10. The aim of obtaining sets of PF solutions is to minimize the three objectives, and the trade-off is shown in Figure 10. From Figure 11, the PESA-II method can efficiently obtain sets of PF solutions for lead-acid, AGM, and lithium-ion batteries in terms of good spread, convergence, and coverage. The CPU execution time of the three sets PF solutions for the lead-acid, AGM, and lithium-ion batteries were 19,810.81 s, 4783.03 s, and 1347.06 s, respectively. For more convenience, the minimum and maximum values of the three objectives for the three types of batteries using the PESA-II method are tabulated in Table 7.   It is evident that the optimal PF solutions were efficiently constructed on the basis of the PESA-II method for three types of storage batteries. Therefore, the integrated AHP-VIKOR method was employed in order to arrange the preference of the optimal solutions. According to Tables 4 and 5, the first evaluator gave extreme importance to LLP over LLC, and very strong importance to LCC over LCE. Meanwhile, equality was given between LCC and LCE. The second evaluator assumed equality between LLP and LCC, and very strong importance of LLP over LCE. Contrastingly, it gave extreme importance to LCC over LCE. The third evaluator proposed moderate and strong importance of LLP over LCC and LCE, respectively, whereas it gave moderate importance to LCC over LCE. Thus, values of weights for LLP, LCC, and LCE objectives were (0.79, 0.1, and 0.1) for the first evaluator, (0.48, 0.44, and 0.07) for the second evaluator, and (0.63, 0.26, and 0.1) for the third evaluator as illustrated in Table 5.
The consistency ratios (CR) for the three evaluators were 0.006, 0.01, and 0.033. The minimum CR value was recommended and meant the values of the given weights to the three conflicting objectives were stable and perfectly consistent. The results of the AHP-VIKOR method based on three evaluators for the lead-acid battery are given in Table 8A-C. The optimal configurations were ranked by the Qi value in descending order and Friedman rank for the first 10 to 200 configurations. The first optimal configuration of the SAPV system based on the lead-acid battery for the first evaluator was composed of 250 numbers of PV modules (25 in series and 10 in parallel) and 40 numbers of storage batteries at an LLP equal to 0.00324, as well as the values of LCC and LCE being USD 54,032.19 and USD 1.5679, respectively. The Qi value decreased from 1 to 0.0022 at the first rank from 200 configurations on the basis of the Friedman rank. The CPU execution for the first evaluator was 0.656 s. Meanwhile, the first optimal configuration for the second evaluator based on the lead-acid battery consisted of 190 PV modules (10 in series and 19 in parallel) and 34 storage batteries at 0.03 LLP, and the values of LCC and LCE were $43,276.91 and $1.6524, respectively. The value of Qi at first rank was 0.0103 from 200 solutions using the Friedman rank. The CPU execution time was 0.4062 s. The results of optimal configurations of the third evaluator using the lead-acid battery were given by 216 PV modules (24 in series and 9 in parallel) and 38 storage batteries at LLP equal to 0.0109, and the values of LCC and LCE were $48,330.29 and $1.6232, respectively. The Qi value was 0.0059 at first rank of 200 configurations on the basis of the Friedman rank method. The CPU execution of the third evaluator for the lead-acid battery was 0.5156 s. The results of the AHP-VIKOR method based on the laid battery for three evaluators can be found in Table 8A-C. It can be observed that the value of Qi was largely dependent on the given weight values by the evaluators for the three objectives. Moreover, the solution gained with min Si is a maximum group utility and the solution gained Ri is with a minimum individual of regret of the "opponent". Thus, it had a proportional relationship with Si and Ri. The minimum value of Qi was registered in the first evaluator at 0.0022, then followed by third and second evaluators. This confirms that each configuration was ranked according to the weight values. The PESA-II method showed its ability to achieve a trade-off between LLP, LCC, and LCE objectives. This fact can be seen in the values of Qi within the first 10 configurations of rank, as shown in Table 8A-C. In Table 8B, the values of LLP were reduced from 0.03 at 1st rank to 0.003 at 10th rank, and the LCC values were increased from $43,276.91 to $52,713.91. The LCE values were changed on the basis of the energy conducted by the PV modules and the capacity of the storage batteries during the period of time.     Figure 11. Scattering of N and Bat at predefined levels of LLP using the PESA-II method for three types of batteries.  In a similar manner to Tables 9 and 10, the Qi values had minimum values with the first evaluator at a high level of accuracy according to the given weight value to LLP. The Qi values at 1 st rank for the first evaluator for the AGM and lithium-ion storage batteries were 0.0019 and 0.0038, respectively. Meanwhile, the Qi values at the first rank for the second and third evaluators based on AGM and lithium-ion batteries were 0.0094, 0.0086, 0.0181, and 0.0087, respectively. From Table 10 (part A), it can be noted that there were four configurations with the rank of 4.5; this happened due to interdependence in total numbers of PV modules and storage batteries. The configurations in Table 10, some configurations of Table 9, and some configurations in Table 8 (part B) had undesirable levels of LLP. This was because the second evaluator gave moderate importance between LLP and LCC objectives, as observed from Qi values that were larger than others. Therefore, a higher priority must be given to the LLP objective in designing of SAPV system. The CPU_execution times for the three evaluators based on AGM and lithium-ion storage batteries were 0.5468, 0.4531, 0.6093, 0.5781, 0.3437, and 0.4218 s. It was clear that the computational times of the three evaluators based on the lithium-ion battery were less than AGM and lead-acid. This was because of the large capacity of the lithium-ion battery, which is about 6.6 kWh compared with other batteries. However, the lithium-ion battery is very expensive compared with AGM and lead-acid batteries at the same level of reliability, as illustrated in Tables 8-10 (part A). Moreover, the levelized cost of energy values presented by the lead-acid battery were less than other batteries, which makes the lead-acid battery more suitable in the designing of a SAPV system. However, more replacement times in the lead-acid battery are required. In this regard, it is prudent to choose the lead-acid battery for the three evaluators at first rank to analyze the performance of the SAPV system for one year.   The daily performance of the SAPV system under the optimal configuration, which was composed of 250 PV modules (25 in series and 10 in parallel) and 40 lead-acid storage batteries, is illustrated in Figure 12. The first subplot shows solar irradiation (G) using hourly meteorological data for one year. Meanwhile, the output power of PV modules (P_PV) and load demand (E_L) are illustrated in the second subplot; the average G and P_PV over one year were around 4363.8 Wh/m 2 /day and 2938.4 KWh, respectively. The state of charge (SOC), deficit energy (E_deficit), and dumped energy (E_dump) are shown in the third, fourth, and fifth subplots, respectively. According to Figure 12, the SAPV system was usually able to fulfill the required load demand through the year, except for 12 days, as demonstrated in Table 11. The deficit energy for one year was about 72 kWh/year, which is equal to 28 h under the chosen optimal configuration at 0.0032 LLP. In the third subplot, the SOC was almost at high levels, except for the days listed in Table 11. The reason for this is that level of solar irradiance is low, especially in December, as shown in the fourth subplot. The last subplot indicated excess energy for one year. From Table 11, the lowest average daily solar irradiance was registered on 24 December at about 1591.7 Wh/m 2 /day. Therefore, the output power of PV modules, SOC, and deficit energy were 36 KWh, 471.8 KWh, and 4.5 KWh, respectively. In fact, the shortages almost happened in the early morning before the appearance of sunshine due to the low solar irradiance before one day, and the amount of energy stored in batteries was very little. On the other hand, the maximum solar irradiance occurred on 1 March and was about 6638.8 Wh/m 2 /day. Thus, the output power of PV modules, SOC, and excess energy were 140 KWh, 975 KWh, and 79.6 KWh, respectively.     Figure 13 shows the monthly state of charge, the output power of PV modules, load demand, excess energy, deficit energy, and ambient temperature for one year. As stated previously, the maximum solar irradiance occurred in March, whereas the minimum solar irradiance was registered in December, as illustrated in Figure 13. Last but not least, several observations were noted, as follows: • Employing the intuitive method can significantly reduce the computational time of the iterative approach, which was found to be very accurate compared with other methods.

•
The PESA-II method demonstrated its capability not only in obtaining PF solutions considering three conflicting objectives, but also by reducing the big solutions in the design space using binary tournament selection.

•
An effective and integrated AHP-VIKOR method allowed for the reaching of the final optimum configuration and selected the most appropriate storage battery for the SAPV system. • The performance of the proposed SAPV system was analyzed on the basis of hourly metrological data for one year. Last but not least, several observations were noted, as follows: • Employing the intuitive method can significantly reduce the computational time of the iterative approach, which was found to be very accurate compared with other methods.

•
The PESA-II method demonstrated its capability not only in obtaining PF solutions considering three conflicting objectives, but also by reducing the big solutions in the design space using binary tournament selection.

•
An effective and integrated AHP-VIKOR method allowed for the reaching of the final optimum configuration and selected the most appropriate storage battery for the SAPV system.

•
The performance of the proposed SAPV system was analyzed on the basis of hourly metrological data for one year. • Overall, our proposed novel iterative-hybrid method tended to be very accurate compared with metaheuristic multi-objective optimization methods due to its random nature, which led to the omission of several important solutions.
It is worth mentioning that there were some limitations that should be considered in future work, such as selecting optimal configurations on the basis of PESA-II to establish set PF solutions, which can be improved by using another method that has a better capacity to deal with large data conducted by a numerical method. Moreover, the given weights' values by using AHP method can be improved because the AHP method requires manual calculations.

Conclusions
This paper presented a novel approach by integrating iterative-PESA-II and AHP-VIKOR methods to find an optimum configuration of the SAPV system on the basis of techno-economic criteria and in considering various types of storage battery. The results demonstrated that the desirable and most suitable configuration for the SAPV system was chosen according to the first expert based on lead-acid battery, which consisted of 250 PV modules (25 in series and 10 in parallel) and 40 storage batteries. The values of Qi, LLP, LCC, and LCE of the optimum configuration were 0.0022, 0.0032, 54032.19, and 1.56 USD, respectively. Therefore, the performance analysis study based on the most favorable configuration was also presented for one year. The optimal configuration of the SAPV system consisted of 872 PV modules and 28 storage batteries. The proposed new hybrid method showed high accuracy in finding the optimal configuration compared with meta-heuristic and numerical methods. The proposed hybrid method showed that the SAPV system could not meet the required load demand only 28 h of one year, which is a promising method to tackle optimization problems.