Optimal Design of Standalone Photovoltaic System Based on Multi-Objective Particle Swarm Optimization: A Case Study of Malaysia

: This paper presents a multi-objective particle swarm optimization (MOPSO) method for optimal sizing of the standalone photovoltaic (SAPV) systems. Loss of load probability (LLP) analysis is considered to determine the technical evaluation of the system. Life cycle cost (LCC) and levelized cost of energy (LCE) are treated as the economic criteria. The two variants of the proposed PSO method, referred to as adaptive weights PSO ( AWPSO cf ) and sigmoid function PSO ( SFPSO cf ) , are implemented using MATLAB software to the optimize the number of PV modules in (series and parallel) and number of the storage battery. The case study of the proposed SAPV system is executed using the hourly meteorological data and typical load demand for one year in a rural area in Malaysia. The performance outcomes of the proposed AW / SFPSO cf methods give various conﬁgurations at desired levels of LLP values and the corresponding minimum cost. The performance results showed the superiority of SFPSO cf in terms of accuracy is selecting an optimal conﬁguration at ﬁtness function value 0.031268, LLP value 0.002431, LCC 53167 USD, and LCE 1.6413 USD. The accuracy of AW / SFPSO cf methods is veriﬁed by using the iterative method.

of the system was employed as an objective function without considering the availability of the system. However, the proposed GA compared with LPSP method and worst month method based on cost criteria. The results showed the superiority of the GA method. Chokri et al. [23] presented FL method to supply energy for a SAPV system in Safax, in the coast of Tunisia. Monthly average of the daily solar radiation was utilized as input of the proposed sizing algorithm. However, the cost and level of reliability of the system have not been mentioned. Moreover, the authors of [22] proposed GRNN method to predict optimal size of the a SAPV system in five Malaysian sites. Hourly meteorological data and load demand were used in this study. However, the proposed GRNN has been evaluated based on a single objective function which is the loss of load probability (LLP). Finally, Nur et al. [24] proposed firefly algorithm sizing algorithm (FLSA) for sizing of a SAPV system. LPSP was used as a technical parameter whereas the economic parameter was not mentioned. The FLSA was found faster than PSO, evolutionary programming (EP) and GA. Some of software tools have been utilized for sizing of a SAPV system based on economic evaluation, and planning and analysis tools such as PV.sys, Transient System Simulation (TRNSYS), Hybrid Optimization Model for Electrical Renewable (Homer), improved Hybrid Optimization by Genetic Algorithm (iHOGA), RetScreen, and System Advisor Model (SAM) [25]. Several studies have employed hybrid methods by combining two methods or more to find an optimal design of the SAPV system [26][27][28][29]. The hybrid methods overcome the drawbacks of AI methods by high ability of convergence to optimal solutions and reducing the execution time. In [30], a hybrid method was proposed to optimal size of a SAPV system by using differential evolution multi-objective optimization (DEMO) and multi-criteria decision making (MCDM). Set of configurations were conducted by using DEMO algorithm. Then, AHP and TOPSIS were obtained to select optimal configuration. Techno-economic criteria for minimizing two objectives was utilized to describe the performance of the proposed method. However, the complexity is one of the drawbacks of the proposed method. Furthermore, the numbers of series PV and battery are assumed to be constant during the optimization. Fabian et al. [31] proposed a computational method for designing off-grid SAPV system by analysing the suppressed demand (SD) effect in south-west highlands, Bolivia. Three types of electrical load profiles were used based on simulated geographical data. The LPSP was employed as technical evaluation using two different capacity of the lithium-ion battery. Whist, the cost of the components of the system was set at 2.5 USD/Wp. The major drawback of the proposed method in Ref. [31] is by supplying the required load demand only by energy conducted from storage battery. This will lead not only increasing the replacements of the storage battery but also raising the cost of the system.
Based on the reviewed literature, multi-objective PSO (MOPSO) method has been not applied based on three conflicts objectives for SAPV systems. The two variants of MOPSO methods are proposed to overcome the obstacles of accurate results of intuitive and analytical methods as well as the long execution time of the numerical methods. In addition, MOPSO method can deal with three objectives synchronously based on techno-economic criteria. Finally, MOPSO method offers a set of optimal configurations by considering the given priority weights to the individual objectives. Therefore, this methodology intends to suggest optimal number of PV array which comprises of parallel and series PV modules and optimal number of storage battery. The technical criteria is calculated according to the method of the loss of load probability (LLP) and the economic criteria is computed according to the method of life cycle cost (LCC) and levelized cost of energy (LCE). The three conflicting objectives: LLP, LCC, and LCE are weighted, normalized, and then aggregated by a mono-objective function which is required to be minimized to select an optimal solution as the best solution. The performance of the proposed AW/SFPSO c f are analysed using hourly meteorological data and typical load demand. The validation of the two variation of PSO algorithms are verified by using numerical method. The rest of the paper is organized as follow. Section 2 presents the definitions, concepts, steps of modeling the SAPV systems. The techno-economic evaluations criteria of the SAPV system are given in Section 3. The particle swarm optimization and the details of the proposed multi-objective optimization approach

Single Diode (SD) Model
The single-diode PV model offers high ability in reflecting the reality behavior when it is dealing with non-linearity and stochastic nature among other mathematical models [33]. The electrical equivalent circuit of a single-diode (SD) PV module model is shown in Figure 2. The circuit comprises an ideal dependent current source, which is the photocurrent ( ) of solar cell and has high sensitivity to randomly nature of the meteorological data. Furthermore, the output voltage of the solar cell is given by the connection of the parallel reverse mode diode with the current source. Parallel resistance and represent the losses of the PV solar cell. Therefore, the output current of the PV solar cell is represented by the following equation: where refers to the photocurrent, denotes the diode reverse saturation currents of the diode (in A), and are the series and shunt resistances (in Ω), respectively, is output voltage (in V) and is output current (in A) of the solar cell, respectively, and represents the diode thermal voltage (in V) which is given by, = where refers to the diode ideality factor that represents the diffusion current's components, is the Boltzmann's constant (1.3806503×10 −19 J/K), represents the cell temperature ( ) and refers to the electron charge (1.60217646×10 −19 Coulombs) where , , , and of Equation (1) are the five parameters that require great computations efforts to be extracted. These five parameters can also be sensitive to the operational conditions.

Single Diode (SD) Model
The single-diode PV model offers high ability in reflecting the reality behavior when it is dealing with non-linearity and stochastic nature among other mathematical models [33]. The electrical equivalent circuit of a single-diode (SD) PV module model is shown in Figure 2. The circuit comprises an ideal dependent current source, which is the photocurrent (I Ph ) of solar cell and has high sensitivity to randomly nature of the meteorological data. Furthermore, the output voltage of the solar cell is given by the connection of the parallel reverse mode diode with the current source. Parallel resistance and R s represent the losses of the PV solar cell. Therefore, the output current of the PV solar cell is represented by the following equation: where I Ph refers to the photocurrent, I o denotes the diode reverse saturation currents of the diode (in A), R s and R p are the series and shunt resistances (in Ω), respectively, V is output voltage (in V) and I is output current (in A) of the solar cell, respectively, and V t represents the diode thermal voltage (in V) which is given by, where d refers to the diode ideality factor that represents the diffusion current's components, KB is the Boltzmann's constant (1.3806503×10 −19 J/K), T c represents the cell temperature (K) and q refers to the electron charge (1.60217646×10 −19 Coulombs) where I Ph , I o , R s , R p and d of Equation (1) are the five parameters that require great computations efforts to be extracted. These five parameters can also be sensitive to the operational conditions.  The PV array comprises of series and parallel modules, respectively. Therefore, the PV array's output current can be expressed as follows: where; , , , and of Equation (1) are the five parameters that extracted by the proposed Improved electromagnetism-like algorithm (IEM) algorithm the output I-V and P-V characteristics curves are shown in Figure 3a,b, , are the output voltage (A) and current (V) of the PV array, respectively.  The PV array's efficiency can be calculated as follows: where is the amount of the solar irradiation fall on the PV array surface W/m 2 and represents the PV array's area. Kyocera (KC120-1) multi-crystalline silicon PV module are utilized in this study. The PV module's specification under standard test condition (STC) are given in Table 1. It's important to notice that the parameters of SD PV model using IEM algorithm were extracted at maximum power point (MPP) which cannot translate the actual performance of the PV system. Thus, it needs to shift to applicable irradiance and temperature levels such as: 840 (W/m 2 ) and 32.170 °C, respectively. Chenlo et al. [34] used five equations for translation curve from MPP to any new point. The root mean The PV array comprises of series N p and parallel N s modules, respectively. Therefore, the PV array's output current can be expressed as follows: where; I Ph , I 0 , R s , R p and d of Equation (1) are the five parameters that extracted by the proposed Improved electromagnetism-like algorithm (IEM) algorithm the output I-V and P-V characteristics curves are shown in Figure 3a  The PV array comprises of series and parallel modules, respectively. Therefore, the PV array's output current can be expressed as follows: where; , , , and of Equation (1) are the five parameters that extracted by the proposed Improved electromagnetism-like algorithm (IEM) algorithm the output I-V and P-V characteristics curves are shown in Figure 3a  The PV array's efficiency can be calculated as follows: where is the amount of the solar irradiation fall on the PV array surface W/m 2 and represents the PV array's area. Kyocera (KC120-1) multi-crystalline silicon PV module are utilized in this study. The PV module's specification under standard test condition (STC) are given in Table 1. It's important to notice that the parameters of SD PV model using IEM algorithm were extracted at maximum power point (MPP) which cannot translate the actual performance of the PV system. Thus, it needs to shift to applicable irradiance and temperature levels such as: 840 (W/m 2 ) and 32.170 °C, respectively. Chenlo et al. [34] used five equations for translation curve from MPP to any new point. The root mean The PV array's efficiency can be calculated as follows: where G T is the amount of the solar irradiation fall on the PV array surface W/m 2 and A represents the PV array's area. Kyocera (KC120-1) multi-crystalline silicon PV module are utilized in this study. The PV module's specification under standard test condition (STC) are given in Table 1. It's important to notice that the parameters of SD PV model using IEM algorithm were extracted at maximum power point (MPP) which cannot translate the actual performance of the PV system. Thus, it needs to shift to applicable irradiance and temperature levels such as: 840 (W/m 2 ) and 32.170 • C, respectively.
Chenlo et al. [34] used five equations for translation curve from MPP to any new point. The root mean square error (RMSE) difference value between the measurement and simulation data was found 1% by Hadj Arab et al. [35]. However, the used equations in [34,35] distorted the shape of simulated new point of the I-V curve and the simulation results demonstrated that solely two equations are required to establish the I-V curve at under any weather conditions: Table 1. Standalone photovoltaic (SAPV) system's key specification [36].

Components Characteristics Value
Kyocera

Mathematical Model of the Battery
The storage battery capacity (kW) is simulated considering the required load demand and the days of the autonomy which can be given by [37], where E Load is the load demand, AD is the autonomy days (typically 3-5 days) [38], DOD is the depth of discharge (80%), ηb * ηinv are the battery (85%) and inverter efficiencies (95%). On the other hand, the state of chare of the storage battery (SOC) is modeled as the following [16]: where SOC(t) and SOC(t − 1) represent the SOC for the storage battery at final and initial points of charging and discharging, respectively, and E bat refers to hourly capacity of the battery. The minimum energy in the storage battery can be written as: The annual amount of energy considered to be stored in storage battery can be given by [38], where ηch denotes to the storage battery's charging efficiency.

Evaluations Criteria for the Sizing of Standalone PV System
In SAPV system, the chosen optimal configuration can significantly affect the performance of the system to meet the required load demand. In this case study, the optimal configuration is based on minimizing three objectives, namely, loss of load probability (LLP), life cycle cot (LCC), and levelized cost of energy (LCE).

Loss of Load Probability (LLP)
The LLP indicates the availability of standalone PV/battery system. It can be defined as the ratio of annual energy deficits to annual required load demands throughout specific time period [39], the LLP is expressed as below:

Life Cycle Cost (LCC)
The life cycle cost (LCC) of the SAPV system is considered in this research work to be the best criteria of economic profitability of the system cost analysis, following the outcomes of previous research. The LCC can be defined as the total cost of the all components of the standalone PV system. There are four main parts proposed in this system: PV array, storage battery, bidirectional converter, and charge regulator. The LCC takes into calculation the initial capital cost IC cap , the present value of replacement cost C rep , and the present value of operation and maintenance cost (C O&M ). The lifetime and unit cost of the system's components are tabulated in Table 2. Therefore, the mathematical equation can be written as [21]: The initial capital cost of each part of the system components is taken into account which comprises of the component price, installation, and the cost of civil works. The initial cost of the off-grid SAPV system IC cap is given by [40]: where C PV, C Unit,PV represent the total capacity (W) and unit cost ($/W) of the PV array, respectively; C Batt , C Unit,Batt represent the total capacity (W) and unit cost ($/W) of the battery, respectively; C Bidi , C Unit,Bidi represent the total capacity (W) and unit cost ($/W) of bidirectional inverter, respectively; C CH , C Unit,CH represent the total capacity (W) and unit cost ($/W) of the charge regulator, respectively; and (C O ) represents the total constant cost including both civil work and installation cost. The present value of the replacement cost of each part of the system component is defining as the present value of all the replacement costs during the system lifetime. Since the lifetime of the storage battery, charge controller, and the bidirectional inverter are shorter than PV arrays, the replacement cost of the components of the SAPV system with the inflation rate (FR) and the real interest rate of components replacements (IR) are included in this calculation. Therefore, the present value of replacement cost (C rep ) can be described as [40]; where C Unit is unit component cost of storage battery in ($/W), charge regulator in ($/W) and bidirectional inverter in ($/W); (C nom ) is the nominal capacity of the replacement system component such as storage battery in (Wh), charge controller in (W) and the bidirectional inverter in (W), and (N rep ) is the number of replacement of each component over the system life period (LP). The calculation of the present value of operation and maintenance cost (C O&M ) of the standalone PV system can be expressed as [41]: where C O&M0 represents the operation and maintenance cost in the first year in ($). The financial data of the proposed SAPV system are given in Table 2 [26,42]. Table 2. Shows the unit cost for standalone PV system's components [26,43].

Components
Cost

Levelized Cost of Energy (LCE)
The third parameter to represent the economic aspects is levelized cost of energy (LCE) which refers to the ratio of the total cost of the system's component to the total energy conducted by a SAPV system during specific period of time, and it is expressed as following [44,45], where; TAC represents the total annual cost of the system components and can be expressed as: E total represents the total energy conducted by a SAPV system during specific life period.

Multi-Objective Particle Swarm Optimization (MOPSO)
The particle swarm optimization (PSO) algorithm is a stochastic optimization technique population-based search to choose the optimal solutions. It was developed in 1995 by Eberhart and Kennedy [46] who have been inspired by the social behavior of the bird flocking and fish schooling. PSO algorithm has several advantages such as; it can easily be programmed, shorter calculation time with high ability of convergence to the optimal values, and robustness to control parameters [47,48]. PSO algorithm presents high convergence rate for single-objective function and can be powerful in finding the Pareto front in multi-objective optimization [49].
In PSO algorithm, the group of random parameters are initiated (search space) to search the optimal parameter values by updating in each generation. At each iteration, a particle is developed by two values which are individual best and global best. The individual best is the best solution achieved by the particle. The global best is the best value obtained among the population in all previous iterations. Each particle has a position which represents the value of the variables and a velocity that moves the direction towards the individual and global bests. The fitness function ( f ) is used to select the best solution from all feasible solutions. In this study, the fitness function is given by obtaining the minimum values of LLP, LCC and LCE. The PSO algorithm comprises of the following steps (see Figure 4). Step 1: initialization The variables of the PSO algorithm are valued and the particles are initialized in the predetermined range of the search space. After that, the fitness function for each particle is calculated. Step 1: initialization The variables of the PSO algorithm are valued and the m particles are initialized in the predetermined range of the search space. After that, the fitness function for each particle is calculated.
Step 2: updating the individual and global best The particle who has best fitness function is chosen as optimal value and stored in pbest and gbest is chosen among the population and updating during iterations.
Step 3: updating the velocity and position The particle's velocity and position in the population is updated using the following equations: where x denotes to particle position and v represents the velocity of the particle in the iteration i. The velocity v is given by the following: where p G represents the best individual particle position and p g is defined as the best global particle position, c 1 and c 2 represent the learning factors which control the individual and global bests. r 1 and r 2 are uniform numbers between 0 and 1. v G i is called inertia and it moves the particle in the same direction with its velocity, c 1 r 1 p G i − x G i is called cognitive part and this part returns the particle to a previous position which has higher individual fitness, and c 2 r 2 p g i − x G i is called social part and this part moves the particle to the best region in the population and follows the best direction of the neighbor.
Step 4: checking the feasibility If the particle exceeds the predetermined range, it will replace by the previous position. The PSO algorithm will stop once the gbest is chosen as the best optimal solution. Otherwise, the steps 2-4 are repeated.
In this study, to avoid the premature convergence of the basic particle swarm optimization algorithm, the PSO c f was proposed to provide a better convergence to the optimal solutions. A contraction factor is used to ensure convergence to a stable point, avoiding velocity clamping.
The PSO c f can be expressed by following issue:

Adaptive Inertia Weight AWPSO c f
The value of w determines the overall convergence and exploitation performance of the proposed PSO algorithm. A big search step can speed up the convergence which may skip optimal solutions. In meanwhile, a small search step can ensure a better convergence accuracy. However, the tradeoff between the exploration and exploitation are necessary in designing nonlinear equations. The w value can be represented by follows: Max iter (22) where i denotes to current number of iterations, Max iter refers to maximum number of generations. From the above equation, the value of w is large at early stage. Then, it decreases as the iteration go on. By this way, the algorithm can hit more accurate solutions at end of iterations.

A Logistic Sigmoid Function SFPSO c f
The second proposed inertia weight uses a logistic sigmoid function which is expressed by the following: where; l represents the maximum value of the curve (l = 1), R defines the steepness of the curve, and s 0 describes the sigmoid midpoint on the x-axis (s 0 = zero). The differences of the fitness function between the current and previous generations is defined by s which is randomly weighted as expressed by the following: where rand is randomly selected between the range of [0,1], X G best and X G−1 best are the best individual vectors for the current and previous generations, respectively.
In this work, the LLP, LCC, and LCE objectives are converted by aggregation function to a mono-objective which treats the multi-objective optimization problems as a mono-objective problem as written by the following: where x represents the decision variables vector, which relates to the search space, k that belongs to the individual objective function number and it is an aggregated function. The range of weights coefficients are between 0 > w i < 1 which represents the relative importance of the k objective function of the problem and it is assumed as: However, the three objective functions are non-scalable. Thus, it's important to normalize and implement the objective function as giving as below [50,51].
where f max i (x) and f min i (x) are the upper and lower boundaries of the ith individual objective function, respectively.

Energy Flow Modeling for SAPV System
In this research, the operation process of energy management of a SAPV 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 PV module, efficiency of charging and discharging of the storage battery, the efficiency of the inverter, converter and wires. Processes 2019, 7, x FOR PEER REVIEW 12 of 24 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: • Case 1: The energy of PV array equals to the energy of the required load demand. • Case 2: The energy of PV array is larger than the energy of the required load demand. • Case 3: The energy of PV array is less than the energy of the required load demand.
In the first case, the battery energy is not used and there is no damped or deficit energy which is equal to zero. In the second case, the energy produced by PV array is larger than the load energy. There is an excess energy and depending on the state of charge of the battery (SOC). 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 to charge the battery and the new value of SOC will be calculated by Equation (9) and the deficit energy is zero. In case 3, if the SOC of the storage battery is more than the minimum value _ as in equation (10) then the storage battery is able to fulfill the load demand. The load will be supplied with/without the energy of PV array and storage battery. Else if SOC of the battery is less 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: • Case 1: The energy of PV array equals to the energy of the required load demand. • Case 2: The energy of PV array is larger than the energy of the required load demand. • Case 3: The energy of PV array is less than the energy of the required load demand.
In the first case, the battery energy is not used and there is no damped or deficit energy which is equal to zero. In the second case, the energy produced by PV array is larger than the load energy. There is an excess energy and depending on the state of charge of the battery (SOC). 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 to charge the battery and the new value of SOC will be calculated by Equation (9) 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 as in equation (10) then the storage battery is able to fulfill the load demand. The load will be supplied with/without the energy of PV array and storage battery. Else if 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 LLP, LCC, and LCE are calculated in each iteration.

Results and Discussion
In this paper, AWPSO c f and SFPSO c f are proposed for finding optimal configuration of the SAPV system. The hourly meteorological data for one year in Klang Valley recorded by Subang Meteorological Station is utilized in the optimization process with latitude (3.12 • ) north and longitude (101.6 • ) east. The meteorological data were obtained from a monitoring system that used solar radiation transmitter of high-stability silicon PV detector model WE300 with an accuracy of (± 1%), while the (± 0.25 • C) is the accuracy of the temperature sensor for the surface of the PV arrays of the same model, air temperature sensor model WE700 with an accuracy of (± 0.1 • C) and range (−50 • C-+50 • C), and current transducer model: CTH-050 with input and output range of (0-50 A (DC), and 4-20 mA) [52]. The proposed SAPV system is used to supply the load demand for house in a remote area. A typical household is illustrated in Figure 6. Annual hourly meteorological data of Kuala Lumpur city, Malaysia are obtained for this purpose. The monthly average of daily solar irradiation and ambient temperature are illustrated in Figure 7.
Processes 2019, 7, x FOR PEER REVIEW 13 of 24 than the minimum value, the system cannot meet the load demand and there is deficit energy. The three-objective function LLP, LCC, and LCE are calculated in each iteration.

Results and Discussion
In this paper, and are proposed for finding optimal configuration of the SAPV system. The hourly meteorological data for one year in Klang Valley recorded by Subang Meteorological Station is utilized in the optimization process with latitude ( 3.12°) north and longitude (101.6°) east. The meteorological data were obtained from a monitoring system that used solar radiation transmitter of high-stability silicon PV detector model WE300 with an accuracy of (± 1%), while the (± 0.25 ℃) is the accuracy of the temperature sensor for the surface of the PV arrays of the same model, air temperature sensor model WE700 with an accuracy of (± 0.1 ℃) and range (−50 °C-+50 °C), and current transducer model: CTH-050 with input and output range of (0-50 A (DC), and 4-20 mA) [52]. The proposed SAPV system is used to supply the load demand for house in a remote area. A typical household is illustrated in Figure 6. Annual hourly meteorological data of Kuala Lumpur city, Malaysia are obtained for this purpose. The monthly average of daily solar irradiation and ambient temperature are illustrated in Figure 7.   Figure 8 illustrates the correlation between the computed and experimental values of the current using Improved electromagnetism (IEM), slap swarm optimization (SSA), EM, repaired adaptive differential evolution (Rcr-IJADE) [53], adaptive differential evolution (IADE) [54], and penalty-based differential evolution (PDE) [55] algorithms. The results demonstrate that the IEM algorithm is very close the to the experimental values. Moreover, the RMSE value is small when the number of data set points is reduced. Finally, the IEM algorithm can predict the system's output energy accurately.   Figure 8 illustrates the correlation between the computed and experimental values of the current using Improved electromagnetism (IEM), slap swarm optimization (SSA), EM, repaired adaptive differential evolution (Rcr-IJADE) [53], adaptive differential evolution (IADE) [54], and penaltybased differential evolution (PDE) [55] algorithms. The results demonstrate that the IEM algorithm is very close the to the experimental values. Moreover, the RMSE value is small when the number of data set points is reduced. Finally, the IEM algorithm can predict the system's output energy accurately. Table 3 Figure 8 illustrates the correlation between the computed and experimental values of the current using Improved electromagnetism (IEM), slap swarm optimization (SSA), EM, repaired adaptive differential evolution (Rcr-IJADE) [53], adaptive differential evolution (IADE) [54], and penaltybased differential evolution (PDE) [55] algorithms. The results demonstrate that the IEM algorithm is very close the to the experimental values. Moreover, the RMSE value is small when the number of data set points is reduced. Finally, the IEM algorithm can predict the system's output energy accurately.   In this study, three decision variables are obtained which are the parallel (N p ) and series N s numbers of PV modules and the number of storage battery (Bat). For fair comparison, the search space of N s and N p decision variables is (4,80) and (4,60) for Bat are considered same in both AWPSO c f and SFPSO c f algorithms. The size of the population is set to be 50. On the other hand, the maximum number of iterations is set to 50 in order to clarify stability of optimal solutions. The value of ∅ is selected to be 2. Figure 9 demonstrates the development of evaluation of aggregation function with maximum iteration of the proposed methods to optimize the size of SAPV system. The SFPSO c f has minimum f value which is 0.031249 at 12 generation. It can be observed that, the SFPSO c f and AWPSO c f algorithms have trivial difference between the fitness function values, which means that the two algorithms are very comparative.  algorithms. The size of the population is set to be 50. On the other hand, the maximum number of iterations is set to 50 in order to clarify stability of optimal solutions. The value of ∅ is selected to be 2. Figure 9 demonstrates the development of evaluation of aggregation function with maximum iteration of the proposed methods to optimize the size of SAPV system. The has minimum value which is 0.031249 at 12 generation. It can be observed that, the and algorithms have trivial difference between the fitness function values, which means that the two algorithms are very comparative.  Table 4 indicates the initializing weights ( , , ), total number of PV module ( ), number of parallel module ( ), number of series module ( ), and number of storage battery ( ). The range of weights sets are [1,9] with step size is 0.1. The acceptable level of LLP is less than 1%. The tradeoff between the defined level of reliability and minimum cost are required in choosing an optimal configuration of a SAPV system. Thus, the optimal values of , , and are 0.4, 0.5, and 0.1 for and algorithms. The number of , and are (47 and 60), (5 and 4), and  Table 4 indicates the initializing weights (W 1 , W 2 , W 3 ), total number of PV module (N), number of parallel module (N p ), number of series module (N s ), and number of storage battery (Bat). The range of weights sets are [1,9] with step size is 0.1. The acceptable level of LLP is less than 1%. The trade-off between the defined level of reliability and minimum cost are required in choosing an optimal configuration of a SAPV system. Thus, the optimal values of W 1 , W 2 , and W 3 are 0.4, 0.5, and 0.1 for AWPSO c f and SFPSO c f algorithms. The number of N p , N s and Bat are (47 and 60), (5 and 4), and (45,44) for the SFPSO c f and AWPSO c f algorithms, respectively. From Table 4, the SFPSO c f and AWPSO c f are very close to each other in terms of variables with bold face. The number of N PV module is increased from 192 to 357 and 348 at end levels of optimization with various sets of weights, which means that the value of f strongly depends on the weight's value. The performance results of the standalone PV system for the three objectives and fitness values are explained in Table 5. According to Table 5 for SFPSO c f and AWPSO c f algorithms. It's worth to note that increasing the value of weight is directly proportional to the relative individual objective and inversely proportional to the decreasing value of weight. The LLP value is dramatically decreased from 0.038 at W 1 = 1 to zero and 0.00007 at W 1 = 8 for the two algorithms corresponding with increasing LLC value as the number of N and Bat increased. The LCE is directly affected by the given weight value which is chosen to be high because of the balance between LLP and LCC objectives.  Figure 10 shows the development of N and Bat with generation for SFPSO c f and AWPSO c f algorithms based on optimal configuration which is denoted by bold face in Tables 4 and 5. Figure 11 demonstrates the optimal configuration of the two algorithms at all sets of weights for N and Bat with favourable level of LLP value. Figure 10 shows the development of and with generation for and algorithms based on optimal configuration which is denoted by bold face in Tables 4 and 5. Figure  11 demonstrates the optimal configuration of the two algorithms at all sets of weights for and with favourable level of LLP value.  algorithms based on optimal configuration which is denoted by bold face in Tables 4 and 5. Figure  11 demonstrates the optimal configuration of the two algorithms at all sets of weights for and with favourable level of LLP value.   Figure 12 is graphed of swarm's motion in 50 particles generations which can be viewed that 50 particles fly from different sets weights. During the optimization process, the particle best and global best will be adapted until achieving high reliability and minimum cost at N = 235 and 240 and Bat = 45 and 44, respectively. The chosen point is called global best that combine the three objectives domain which is a complex task in designing such system. It can be noticed from Figure 12 that particles become close to each other and the convergence is reached to the optimal value for the two SFPSO c f and AWPSO c f algorithms. Figure 11. Evaluation of and at desired level of loss of load probability (LLP) value with different weights sets. Figure 12 is graphed of swarm's motion in 50 particles generations which can be viewed that 50 particles fly from different sets weights. During the optimization process, the particle best and global best will be adapted until achieving high reliability and minimum cost at = 235 and 240 and = 45 and 44, respectively. The chosen point is called global best that combine the three objectives domain which is a complex task in designing such system. It can be noticed from Figure 12 that particles become close to each other and the convergence is reached to the optimal value for the two and algorithms. Figure 12. Evaluation of , , and value with different weights sets.
From Figure 13, the optimal configurations based on and algorithms are close to each other with small difference in the three objectives because of the variation in values of inertia weights during the optimization. It is worth to mentioned that the results for the independent runs can guarantee to convergence to the same optimal fitness value which shows the reliability of proposed two algorithms. From Figure 13, the optimal configurations based on SFPSO c f and AWPSO c f algorithms are close to each other with small difference in the three objectives because of the variation in values of inertia weights during the optimization. It is worth to mentioned that the results for the independent runs can guarantee to convergence to the same optimal fitness value which shows the reliability of proposed two algorithms. However, the drawbacks of aggregation method are that identify only one solution based on given weights for each individual objective. Furthermore, the weights initialization, which are aggregated by the aggregation function, is the main obstacle of aggregation method. It is worth mentioning that the results of the optimization problem that solved by aggregation method strongly depends on the weights. Therefore, other method can be proposed to overcome the limitations of the However, the drawbacks of aggregation method are that identify only one solution based on given weights for each individual objective. Furthermore, the weights initialization, which are aggregated by the aggregation function, is the main obstacle of aggregation method. It is worth mentioning that the results of the optimization problem that solved by aggregation method strongly depends on the weights. Therefore, other method can be proposed to overcome the limitations of the proposed method, and other nonlinear individual objective functions can be suggested and added for aggregation function as future work.
The results of the proposed methods for sizing of a SAPV system are validated by results of the numerical method to demonstrate the accuracy of the presented two methods. The range of the search space are defined by using intuitive method for every possible N s , N p , and Bat within range (1,40), (1,40) and (1,50), respectively. The desired value for LLP is larger than 0.001. some of optimal configurations based on numerical method that corresponding to the results of the SFPSO c f and AWPSO c f algorithms are tabulated in Table 6. It's clear that the numerical method is almost able to meet the optimal configurations of the two mentioned methods except that it's over the range of search space of the numerical method as well as if the optimal configuration exceeds the LLP value. The execution time by using numerical method is 10,295.7187 s and the number of configurations is 59,503. Therefore, a high consistency has been achieved for the two proposed methods when their compared with iterative method. The operation performance of the standalone PV system is analysed based on N s = 5, N p = 47 and Bat = 45 of the SFPSO c f algorithm as it has lower f than the AWPSO c f algorithm. Two methods are selected based on maximum and minimum solar radiation to show the performance of the system. The first month is based on the maximum solar radiation which was accrued in March. The monthly average daily solar radiation of about 5131.5 W/m 2 . The energy generated by the PV modules is about 3000.2154 KWh and the deficit energy is about 0 KWh. Therefore, the LLP value is 0, LCC is 53,167 USD and LCE is 16.5348 USD. The hourly performance of the proposed SAPV system of maximum solar radiation is illustrated in Figure 14. and = 45 of the algorithm as it has lower than the algorithm. Two methods are selected based on maximum and minimum solar radiation to show the performance of the system. The first month is based on the maximum solar radiation which was accrued in March. The monthly average daily solar radiation of about 5131.5 W/m 2 . The energy generated by the PV modules is about 3000.2154 KWh and the deficit energy is about 0 KWh. Therefore, the LLP value is 0, LCC is 53,167 USD and LCE is 16.5348 USD. The hourly performance of the proposed SAPV system of maximum solar radiation is illustrated in Figure 14. Whilst, the minimum solar radiation was registered in December. The monthly average daily solar radiation of about 3779.8 W/m 2 . The energy conducted by the PV modules is 4447.0 KWh and the deficit energy in this month is 20.18 KWh. Thus, the LLP value is 0.0191, LCC is 53,167 USD and LCE is 22.1186 USD. The hourly performance of the proposed SAPV system of minimum solar radiation is illustrated in Figure 15. According to Figures 14 and 15, we can observe that the energy conducted by the PV array are strongly dependent on the availability of the meteorological data. Thus, in March the SOC of the battery is almost at maximum values. Whilst, the SOC in December exceeded the minimum values of the storage battery in only three days. The monthly average daily performance of the proposed SAPV system throughout one year based on optimal configuration can be given in Table 7. The Whilst, the minimum solar radiation was registered in December. The monthly average daily solar radiation of about 3779.8 W/m 2 . The energy conducted by the PV modules is 4447.0 KWh and the deficit energy in this month is 20.18 KWh. Thus, the LLP value is 0.0191, LCC is 53,167 USD and LCE is 22.1186 USD. The hourly performance of the proposed SAPV system of minimum solar radiation is illustrated in Figure 15.
methods are selected based on maximum and minimum solar radiation to show the performance of the system. The first month is based on the maximum solar radiation which was accrued in March. The monthly average daily solar radiation of about 5131.5 W/m 2 . The energy generated by the PV modules is about 3000.2154 KWh and the deficit energy is about 0 KWh. Therefore, the LLP value is 0, LCC is 53,167 USD and LCE is 16.5348 USD. The hourly performance of the proposed SAPV system of maximum solar radiation is illustrated in Figure 14. Whilst, the minimum solar radiation was registered in December. The monthly average daily solar radiation of about 3779.8 W/m 2 . The energy conducted by the PV modules is 4447.0 KWh and the deficit energy in this month is 20.18 KWh. Thus, the LLP value is 0.0191, LCC is 53,167 USD and LCE is 22.1186 USD. The hourly performance of the proposed SAPV system of minimum solar radiation is illustrated in Figure 15. According to Figures 14 and 15, we can observe that the energy conducted by the PV array are strongly dependent on the availability of the meteorological data. Thus, in March the SOC of the battery is almost at maximum values. Whilst, the SOC in December exceeded the minimum values of the storage battery in only three days. The monthly average daily performance of the proposed SAPV system throughout one year based on optimal configuration can be given in Table 7. The According to Figures 14 and 15, we can observe that the energy conducted by the PV array are strongly dependent on the availability of the meteorological data. Thus, in March the SOC of the battery is almost at maximum values. Whilst, the SOC in December exceeded the minimum values of the storage battery in only three days. The monthly average daily performance of the proposed SAPV system throughout one year based on optimal configuration can be given in Table 7. The proposed algorithm AWPSO c f exhibits excellent performance results in term of reliability, cost-effective, and CPU-execution time. From Table 7, the difference between the conducted energy by PV modules, deficit energy, and dump energy are reasonable and promising for real investment of the SAPV system especially in rural area.

Conclusions
In this research, SFPSO c f and AWPSO c f methods for optimally sizing SAPV system were proposed. Three constraint objectives are utilized as techno-economic criteria which are normalized, weighted, and then aggregated by mono-objective function. The optimum design was chosen by minimizing objective function based on LLP, LCC, and LCE criteria. The results showed the superiority of SFPSO c f in term of accuracy which has minimum objective function with value 0.031249. The optimal configuration of SAPV system for SFPSO c f method was 5 connected in series and 47 in parallel, respectively. The number of storage battery was 45. The performance was analysed for one year using hourly meteorological data in Malaysia. The LLP value is 0.002431, LCC is 53,166.9 USD and LCE is 1.6431 USD. March and December, the months of the maximum and the minimum solar radiation periods, respectively, were selected to show the operational performance of SAPV system. The results of operational performance show high reliability of the standalone PV system; hence, the acceptable level of LLP value is less than 1%.
Author Contributions: C.G. and H.H. were involved in the paper editing and reviewing, conceptualization, investigation, and supervision. H.M.R. collected the data, methodology, software implementation, and writing original draft. M.A. and H.M.R. conducted formal analysis, conceptualization, and investigation. All authors have read and agreed to the published version of the manuscript.