Impact of Di ﬀ erent Photovoltaic Models on the Design of a Combined Solar Array and Pumped Hydro Storage System

: The impact of di ﬀ erent photovoltaic models for a combined solar array and pumped hydro storage system was investigated. Al-Wehda dam located in Harta city in the northern of Jordan was used to validate the approach. The two-diode (TD), single-diode (SD), and ideal single-diode (ISD) solar models were evaluated in terms of the solar array size, reliability, and ecological e ﬀ ects. The impoundment of Al-Wehda dam was taken as the upper reservoir of the pumped hydro facility of the proposed renewable energy system. It was found that the PV power is more accurately modelled by considering the recombination loss in the TD solar model. This leads to a more realistic sizing and precise system evaluation. Results were obtained using the particle swarm optimization (PSO) algorithm and the whale optimization algorithm (WOA) for validation purposes. For instance, the PSO results showed that the realistic TD model is reliable, with an index of reliability of 98.558%. Further, it is the most ecological solution with an annual emissions reduction of 21.5198 Gg. The optimized values are 44,840 solar panels and 65.052 M.m 3 of the lower reservoir volume for the TD model. The number of PV panels are reduced by 16.67% and 7.93%, respectively, with the ISD and SD relative to the TD model.


Introduction
Renewable power plants are the long-term solution for oil importing countries. This is especially true for countries struggling to satisfy their energy needs, such as Jordan. It is well known that photovoltaic (PV) renewable energy systems (RESs) are clean, modular, and quiet resources of energy. Further, the World Climate Conference that was held in Paris in 2015 set the limit for the global average temperature rise to be 0.2 • C per decade, between the years 2000 and 2100. Furthermore, hybrid systems that couple PV with other RES topologies have been considered worldwide [1][2][3][4][5]. Moreover, PV arrays, combined with an energy storage technology, such as pumped hydro storage (PHS), provide a reliable source of electricity [6][7][8][9][10]. Therefore, an in-depth investigation should be system (ESS). Different load demand and weather conditions were studied with the aim of minimizing the system's operation cost. In [24], the authors presented an optimal energy management model of a standalone system, including a solar PV model using an ISD, diesel generator, and battery storage system. Their goal was to meet the load demand completely whilst satisfying the system constraints. In [25], the authors presented a method to find the optimal size for a hybrid solar PV and wind farm combined with battery storage. The PV output power was determined using the ISD-equivalent circuit. However, their goal was to minimize the total cost of the system, while maximizing power supply reliability. Furthermore, different constraints on the standalone system operation were considered while minimizing the fluctuations of power supplied to the utility grid. In [2], the researchers optimized a small-scale HRES using simulated annealing as a heuristics global approach. The system consists of solar and wind plants with batteries as ESS. The objective was to minimize the total cost of the HRES. Other reliability indicators, such as the loss of load probability (LOLP) and index of reliability (IR), were found at the optimum point in order to confirm the reliability of the system.
In summary, much of the research of HRES studies use the ISD or SD models for PV system design, for simplicity. However, the PV system designed in this paper was done using three different PV models for comparison purposes on a power system level. In other words, the PV plant, in this investigation, was modelled using the ISD, SD, and TD models. Then, the plant size, the index of reliability (IR), and detailed reliability results were computed for a grid-connected system, for each PV model for comparison purposes. It was shown that the TD solar model provides realistic and accurate sizing of the PV array, and hence the size and reliability of the entire system can be precisely found.
In this paper, an on-grid system including a PV array combined with PHS, as shown in Figure 1, was optimized using the PSO algorithm and whale optimization algorithm (WOA). The main components of the grid-connected system include the PV array, PHS sub-system showing the upper and lower reservoir, AC and DC bus with the appropriate power electronics converters, and the load demand, which is connected to the AC bus.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 3 of 21 generator, and battery bank as an energy storage system (ESS). Different load demand and weather conditions were studied with the aim of minimizing the system's operation cost. In [24], the authors presented an optimal energy management model of a standalone system, including a solar PV model using an ISD, diesel generator, and battery storage system. Their goal was to meet the load demand completely whilst satisfying the system constraints. In [25], the authors presented a method to find the optimal size for a hybrid solar PV and wind farm combined with battery storage. The PV output power was determined using the ISD-equivalent circuit. However, their goal was to minimize the total cost of the system, while maximizing power supply reliability. Furthermore, different constraints on the standalone system operation were considered while minimizing the fluctuations of power supplied to the utility grid. In [2], the researchers optimized a small-scale HRES using simulated annealing as a heuristics global approach. The system consists of solar and wind plants with batteries as ESS. The objective was to minimize the total cost of the HRES. Other reliability indicators, such as the loss of load probability (LOLP) and index of reliability (IR), were found at the optimum point in order to confirm the reliability of the system. In summary, much of the research of HRES studies use the ISD or SD models for PV system design, for simplicity. However, the PV system designed in this paper was done using three different PV models for comparison purposes on a power system level. In other words, the PV plant, in this investigation, was modelled using the ISD, SD, and TD models. Then, the plant size, the index of reliability (IR), and detailed reliability results were computed for a grid-connected system, for each PV model for comparison purposes. It was shown that the TD solar model provides realistic and accurate sizing of the PV array, and hence the size and reliability of the entire system can be precisely found.
In this paper, an on-grid system including a PV array combined with PHS, as shown in Figure 1, was optimized using the PSO algorithm and whale optimization algorithm (WOA). The main components of the grid-connected system include the PV array, PHS sub-system showing the upper and lower reservoir, AC and DC bus with the appropriate power electronics converters, and the load demand, which is connected to the AC bus.

System's Measured Data
Hourly data is required to run the optimization algorithms for the RES shown in Figure 1. The data was taken in 2018 for Harta city in the north of Irbid, Jordan, as shown in Figures 2 and 3, respectively. The sample time of the data acquisition is one hour. Note that the optimization algorithms are hourly testing the system in order to find the optimal system's configuration. Further, the duration of campaign acquisition for the measured data was one year, i.e., 8760 h.

System's Measured Data
Hourly data is required to run the optimization algorithms for the RES shown in Figure 1. The data was taken in 2018 for Harta city in the north of Irbid, Jordan, as shown in Figures 2 and 3, respectively. The sample time of the data acquisition is one hour. Note that the optimization algorithms are hourly testing the system in order to find the optimal system's configuration. Further, the duration of campaign acquisition for the measured data was one year, i.e., 8760 h.

System's Measured Data
Hourly data is required to run the optimization algorithms for the RES shown in Figure 1. The data was taken in 2018 for Harta city in the north of Irbid, Jordan, as shown in Figures 2 and 3, respectively. The sample time of the data acquisition is one hour. Note that the optimization algorithms are hourly testing the system in order to find the optimal system's configuration. Further, the duration of campaign acquisition for the measured data was one year, i.e., 8760 h.   The measured load demand (in MW) was obtained from the Jordanian National Electric Power Company (NEPCO), which is responsible for the availability and safety of power transmission. Further, the volumetric flow rate (in m 3 /s) input values into the PHS unit for Al-Wehda Dam was obtained from the Jordan Valley Authority (JVA). JVA is responsible for the public water supply, wastewater services, and development of water resources and dams in Jordan. Moreover, the solar incident irradiance, G(t) in W/m 2 , and the solar ambient temperature, T amb (t) in • C, were obtained from the royal scientific society (RSS), which is responsible for testing services locally and internationally. Note that G(t) and T amb (t) are used to compute the photo current (I Ph ) of the solar module as discussed in [6]. It is worth noting that the solar ambient temperature is the temperature of the solar panel surroundings. They were measured in front of the solar panel. Furthermore, RSS performs research studies related to RES in Jordan, such as the optimal tilt angle of a fixed PV array. Note that the designed PV array is to be positioned toward south in a fixed position. A tilt angle of 28 • is optimal for PV installation in Jordan [26]. Note that the three institutions provided the hourly data, shown in Figures 2 and 3, after requesting official letters from the employer of the authors.

Mathematical Modeling and Design Methodology
In order to truly size HRES, the proper design of the PV system is important to determine the correct number of PV modules [15]. Therefore, a basic introduction of the principles of PV cells and their mathematical relations is essential. Therefore, modelling and simulation of the PV modules and arrays are presented here.
As was mentioned, the objective of this work was a comparative analysis of the ISD, SD, and TD models of the PV module in a PV array, combined with PHS. This is to enhance the conversion efficiency of power from the PV system. For this, the optimal system size, ecological impact, and reliability will be studied in the three PV scenarios.

Modeling of a PV Array Output Power Generation
The PV cell is a semiconductor device whose p-n junction is exposed to light, which produces an electric current [27]. However, the design of the PV system depends on the load power. The mathematical models and the corresponding circuit models for ISD, SD, and TD, are investigated here. Furthermore, the physical behavior of each model and their characteristic solar curves are studied.

ISD Model
This model is the most simplified form of an ideal PV module as shown in Figure 4. The output current as a function of the voltage is expressed by (1) [28]: Appl. Sci. 2020, 10, x FOR PEER REVIEW 5 of 21 The measured load demand (in MW) was obtained from the Jordanian National Electric Power Company (NEPCO), which is responsible for the availability and safety of power transmission. Further, the volumetric flow rate (in m 3 /s) input values into the PHS unit for Al-Wehda Dam was obtained from the Jordan Valley Authority (JVA). JVA is responsible for the public water supply, wastewater services, and development of water resources and dams in Jordan. Moreover, the solar incident irradiance, G(t) in W/m 2 , and the solar ambient temperature, Tamb(t) in °C, were obtained from the royal scientific society (RSS), which is responsible for testing services locally and internationally. Note that G(t) and Tamb(t) are used to compute the photo current (IPh) of the solar module as discussed in [6]. It is worth noting that the solar ambient temperature is the temperature of the solar panel surroundings. They were measured in front of the solar panel. Furthermore, RSS performs research studies related to RES in Jordan, such as the optimal tilt angle of a fixed PV array. Note that the designed PV array is to be positioned toward south in a fixed position. A tilt angle of 28˚ is optimal for PV installation in Jordan [26]. Note that the three institutions provided the hourly data, shown in Figure 2 and Figure 3, after requesting official letters from the employer of the authors.

Mathematical Modeling and Design Methodology
In order to truly size HRES, the proper design of the PV system is important to determine the correct number of PV modules [15]. Therefore, a basic introduction of the principles of PV cells and their mathematical relations is essential. Therefore, modelling and simulation of the PV modules and arrays are presented here.
As was mentioned, the objective of this work was a comparative analysis of the ISD, SD, and TD models of the PV module in a PV array, combined with PHS. This is to enhance the conversion efficiency of power from the PV system. For this, the optimal system size, ecological impact, and reliability will be studied in the three PV scenarios.

Modeling of a PV Array Output Power Generation
The PV cell is a semiconductor device whose p-n junction is exposed to light, which produces an electric current [27]. However, the design of the PV system depends on the load power. The mathematical models and the corresponding circuit models for ISD, SD, and TD, are investigated here. Furthermore, the physical behavior of each model and their characteristic solar curves are studied.

ISD Model
This model is the most simplified form of an ideal PV module as shown in Figure 4. The output current as a function of the voltage is expressed by (1) [28]: However, this model does not provide accurate and realistic I-V and P-V characteristics. This is because Equation (1) does not represent the actual characteristics due to neglecting the two diode However, this model does not provide accurate and realistic I-V and P-V characteristics. This is because Equation (1) does not represent the actual characteristics due to neglecting the two diode losses. Actual PV modules are composed of several connected PV cells, in series with losses associated with leakage currents to the ground. These losses are modelled by a series resistance (R s ) and a shunt resistance (R sh ), which will be presented in Section 3.1.2.

SD Model
The SD circuit model of a PV panel is shown in Figure 5. In the literature [2,4,21,27,29,30], this is also called a five-parameter model. The model takes into account the series and leakage current to ground losses, which were not considered in the ISD model. They are represented by the series and shunt resistances R s and R sh .
Appl. Sci. 2020, 10, x FOR PEER REVIEW 6 of 21 losses. Actual PV modules are composed of several connected PV cells, in series with losses associated with leakage currents to the ground. These losses are modelled by a series resistance ( ) and a shunt resistance ( ), which will be presented in Section 3.1.2.

SD Model
The SD circuit model of a PV panel is shown in Figure 5. In the literature [2,4,21,27,29,30], this is also called a five-parameter model. The model takes into account the series and leakage current to ground losses, which were not considered in the ISD model. They are represented by the series and shunt resistances and . The shunt resistance, , is inversely proportional to the leakage current to the ground, and usually has a large value. Resistor represents the series current losses and the internal resistance of the PV module and is, by design, very small. The PV cells are added together to form a PV module and then interconnected in parallel series to form a PV array [17,31].
This model has been used in grid-connected and standalone PV systems [17,31]. This model is based on the well-known Shockley diode equation, which is also known as the five-parameter model. The relationship between the output current (I) and voltage (V) of the PV panel for a single unit, according to the SD model is expressed by (2) [6,32]: The five key parameters, , , , , and , are given under standard test conditions (STCs) which are 1 kW/m 2 air mass 1.5 (AM1.5) and 25 °C. These parameters can be computed using analytical or numerical solutions, based on the specifications provided by the manufacturer. Further, , , , and are given in [6,17,33]. However, the SD model neglects the recombination losses of the diode, which is a phenomenon that should be considered to accurately model a PV panel. The single diode model assumes a constant value for the ideality factor. The ideality factor is a measure of how closely the diode follows the ideal diode equation of the diffusion process [27,34]. It is usually considered as a function of the voltage across the device. Furthermore, at high voltage recombination in the device, electron-hole pairs are dominant at the surface of the bulk regions. This means that the ideality factor is close to one. However, at lower voltages, recombination in the junction dominates and the ideality factor approaches two [34]. The two exponential diodes model (i.e., TD model) takes into account the recombination phenomena and will be discussed in Section 3.1.3.

TD Model
The TD model, shown in Figure 6, of a PV panel takes into account both the diffusion and the recombination processes by including diodes and [35]. The first diode ( ) simulates the diffusion process of the minority carriers into the depletion layer, while diode ( ) represents the recombination in the space charge region of the P-N junction [36]. The shunt resistance, R sh , is inversely proportional to the leakage current to the ground, and usually has a large value. Resistor R s represents the series current losses and the internal resistance of the PV module and is, by design, very small. The PV cells are added together to form a PV module and then interconnected in parallel series to form a PV array [17,31].
This model has been used in grid-connected and standalone PV systems [17,31]. This model is based on the well-known Shockley diode equation, which is also known as the five-parameter model. The relationship between the output current (I) and voltage (V) of the PV panel for a single unit, according to the SD model is expressed by (2) [6,32]: The five key parameters, I ph , I o , R s , R sh , and V T , are given under standard test conditions (STCs) which are 1 kW/m 2 air mass 1.5 (AM1.5) and 25 • C. These parameters can be computed using analytical or numerical solutions, based on the specifications provided by the manufacturer. Further, I Ph , I o , V T , and I o stc are given in [6,17,33].
However, the SD model neglects the recombination losses of the diode, which is a phenomenon that should be considered to accurately model a PV panel. The single diode model assumes a constant value for the ideality factor. The ideality factor is a measure of how closely the diode follows the ideal diode Equation of the diffusion process [27,34]. It is usually considered as a function of the voltage across the device. Furthermore, at high voltage recombination in the device, electron-hole pairs are dominant at the surface of the bulk regions. This means that the ideality factor is close to one. However, at lower voltages, recombination in the junction dominates and the ideality factor approaches two [34].
The two exponential diodes model (i.e., TD model) takes into account the recombination phenomena and will be discussed in Section 3.1.3.

TD Model
The TD model, shown in Figure 6, of a PV panel takes into account both the diffusion and the recombination processes by including diodes D 1 and D 2 [35]. The first diode (D 1 ) simulates the diffusion process of the minority carriers into the depletion layer, while diode (D 2 ) represents the recombination in the space charge region of the P-N junction [36].
The TD model is a modified and improved version of the SD circuit. The SD circuit did not take into account the effect of recombination as shown in Figure 5. Many RES studies use the SD model for a PV system design [2,4,21,29,30]. However, the TD model is more realistic and accurate, because it takes into account both the diffusion effects and recombination losses. This results in a more precise and realistic mathematical description of a solar panel and leads to a correct size for the PV array and also the size, ecological impact, and reliability of the entire system [14][15][16].
After adding a second diode ( ) in parallel with the first one ( ), Equation (2) becomes as shown in Equation (3) [37,38]. The values of , , and are computed as in [12,16,37]. Thus, the output PV current for the TD model is updated to be as in Equation (3): Figure 6. Equivalent circuit of a PV panel using the TD model.
The optically generated current of the PV module is linearly dependent on the solar irradiance and is also influenced by temperature, according to the equation found in [6,14,38]. The hourly measured values of the solar irradiance and ambient temperature are depicted in Figure 2.
The diode ( or ) saturation currents depend on the temperature and can be expressed as shown in [16,37].
Some research studies provide another approach to get the value of the saturation current of or in the TD model. In fact, the diode saturation current is very small and to simplify the model, in these studies, both of the reverse saturation currents and were set equal as in [16,37,38]. The diode ideality factors and represent the diffusion effects and recombination losses. The recombination has an impact on both the short circuit current and open circuit voltage. The recombination region should be minimized in order to ensure that the p-n junction is convenient to collect the light generated current, and to minimize the losses in the current due to recombination. So, the recombination is a negative effect, and thus, it is a loss. Regarding the diffusion effect, the carrier must be generated within a diffusion length of the junction, so that it will be able to diffuse to the junction before recombining. When a carrier exists with a concentration in the semiconductor, it will move from the higher to the lower concentration area in the diffusion process. The impact of diffusion on the solar cell, will increase the carrier concentrations in a device with time. Additionally, a low diffusion length means that minority carriers disappear from the junction edge quickly due to recombination, thus allowing more carriers to cross and increasing the forward bias current. Consequently, to minimize recombination and achieve a high voltage, a high diffusion length is required. So, the diffusion and recombination effects are interrelated. However, the aforementioned models (i.e., SD and ISD), do not take into account the recombination losses as the TD model does, which is a physical phenomenon.
In accordance with Shockley's diffusion theory, is unity [14,37,39]. However, the value of is flexible. It was found that if (1.2 ≤ ≤ 2), the best match between their proposed model and the practical I-V curve is obtained. These studies state that = 1 and =1. It follows that the variable P can be chosen to be greater than or equal to 2.2 and less than 3. Under the aforementioned considerations, can be simplified as in [14,37,39].

ID1 Rs ID2
Rsh V I Figure 6. Equivalent circuit of a PV panel using the TD model.
The TD model is a modified and improved version of the SD circuit. The SD circuit did not take into account the effect of recombination as shown in Figure 5. Many RES studies use the SD model for a PV system design [2,4,21,29,30]. However, the TD model is more realistic and accurate, because it takes into account both the diffusion effects and recombination losses. This results in a more precise and realistic mathematical description of a solar panel and leads to a correct size for the PV array and also the size, ecological impact, and reliability of the entire system [14][15][16].
After adding a second diode (D 2 ) in parallel with the first one (D 1 ), Equation (2) becomes as shown in Equation (3) [37,38]. The values of I D1 , I D2 , and I Sh are computed as in [12,16,37]. Thus, the output PV current for the TD model is updated to be as in Equation (3): The optically generated current of the PV module is linearly dependent on the solar irradiance and is also influenced by temperature, according to the I Ph Equation found in [6,14,38]. The hourly measured values of the solar irradiance and ambient temperature are depicted in Figure 2.
The diode (D 1 or D 2 ) saturation currents I o1 or 2 depend on the temperature and can be expressed as shown in [16,37].
Some research studies provide another approach to get the value of the saturation current of D 1 or D 2 in the TD model. In fact, the diode saturation current is very small and to simplify the model, in these studies, both of the reverse saturation currents I o 1 and I o 2 were set equal as in [16,37,38].
The diode ideality factors α 1 and α 2 represent the diffusion effects and recombination losses. The recombination has an impact on both the short circuit current and open circuit voltage. The recombination region should be minimized in order to ensure that the p-n junction is convenient to collect the light generated current, and to minimize the losses in the current due to recombination. So, the recombination is a negative effect, and thus, it is a loss. Regarding the diffusion effect, the carrier must be generated within a diffusion length of the junction, so that it will be able to diffuse to the junction before recombining. When a carrier exists with a concentration in the semiconductor, it will move from the higher to the lower concentration area in the diffusion process. The impact of diffusion on the solar cell, will increase the carrier concentrations in a device with time. Additionally, a low diffusion length means that minority carriers disappear from the junction edge quickly due to recombination, thus allowing more carriers to cross and increasing the forward bias current. Consequently, to minimize recombination and achieve a high voltage, a high diffusion length is required. So, the diffusion and recombination effects are interrelated. However, the aforementioned models (i.e., SD and ISD), do not take into account the recombination losses as the TD model does, which is a physical phenomenon.
In accordance with Shockley's diffusion theory, α 1 is unity [14,37,39]. However, the value of α 2 is flexible. It was found that if (1.2 ≤ α 2 ≤ 2), the best match between their proposed model and the practical I-V curve is obtained. These studies state that α 1 +α 2 P = 1 and α 1 = 1. It follows that the variable P can be chosen to be greater than or equal to 2.2 and less than 3. Under the aforementioned considerations, I o can be simplified as in [14,37,39].

Solar Array Model
A single PV panel does not produce enough power to satisfy common applications. A configuration of cells is called a module. In other words, a panel consists of one or more PV cells assembled and pre-wired for field installation [37]. A similar series-parallel configuration of PV panels composes the array or a solar PV plant. The wattage rating of a PV array is usually the maximum DC power output under STC. Figure 7 represents the circuit model of TD for the solar PV array.

Solar Array Model
A single PV panel does not produce enough power to satisfy common applications. A configuration of cells is called a module. In other words, a panel consists of one or more PV cells assembled and pre-wired for field installation [37]. A similar series-parallel configuration of PV panels composes the array or a solar PV plant. The wattage rating of a PV array is usually the maximum DC power output under STC. Figure 7 represents the circuit model of TD for the solar PV array.
A PV array is composed of a number of strings connected in parallel (NPar). The number of PV modules connected in series per string is (Nser). Note that (NPar×Nser) is the number of PV modules (NPV) of the PV array, which is one of the design variables to be optimized by using the proposed methodologies and optimization algorithms in Sections 3.4.3-3.4.5, respectively. The optical and saturation currents for a PV array can be expressed as in [16]. The formula that describes the relationship between the output current (I) and voltage (V) of a PV array is derived as in [6,37,39]. The output power from a PV array can be obtained as shown in Equation (4). Moreover, the values of the series and parallel resistances of the PV array are given by the equations described in [16]:

Characteristics of a PV Panel
In order to accurately size hybrid renewable energy systems (HRESs), including a solar PV array combined with pumped hydro storage, the proper design of the PV system is important to determine the realistic number of PV modules [17]. A "KC200GT" polycrystalline PV panel was considered in this paper. Table 1 presents the characteristics of this PV panel as provided by the manufacturer. Further, the I-V and P-V characteristic curves of the mentioned PV panel for the ISD, SD, and TD models were simulated as shown in Figure 8.
Here, 25 °C is considered as a special case in Figure 8 in order to depict the I-V and P-V characteristics under standard test conditions (STCs), and to show the differences between the ISD, SD, and TD photovoltaic models. Moreover, Figure 8 shows that for a single PV panel, the power that can be extracted is 189.149 W in the case of the TD model. This power increases by 8.43% with the ISD model. So, the TD model gives more realistic and accurate results because it takes into account the realistic phenomenon of recombination losses. This will lead to a true size of the HRES, especially with the large-scale power system, whose load is shown in Figure 3.
In summary, many RES studies use the ISD model for PV system design [32,33]. However, the TD model (shown in Figure 6) is more accurate, because it takes into account both the diffusion and the recombination effects. This results in a more precise mathematical description of the solar panel, leading to a true size of the PV array and true reliability and environmental estimates of the entire system [17]. Figure 7. Equivalent circuit of a PV array using the TD model.
A PV array is composed of a number of strings connected in parallel (N Par ). The number of PV modules connected in series per string is (N ser ). Note that (N Par × N ser ) is the number of PV modules (N PV ) of the PV array, which is one of the design variables to be optimized by using the proposed methodologies and optimization algorithms in Sections 3.4.3-3.4.5, respectively. The optical and saturation currents for a PV array can be expressed as in [16].
The formula that describes the relationship between the output current (I) and voltage (V) of a PV array is derived as in [6,37,39]. The output power from a PV array can be obtained as shown in Equation (4). Moreover, the values of the series and parallel resistances of the PV array are given by the equations described in [16]:

Characteristics of a PV Panel
In order to accurately size hybrid renewable energy systems (HRESs), including a solar PV array combined with pumped hydro storage, the proper design of the PV system is important to determine the realistic number of PV modules [17]. A "KC200GT" polycrystalline PV panel was considered in this paper. Table 1 presents the characteristics of this PV panel as provided by the manufacturer. Further, the I-V and P-V characteristic curves of the mentioned PV panel for the ISD, SD, and TD models were simulated as shown in Figure 8.

Modeling of Hydropower Facilities
There can be different topologies for a PHS system. However, the volume of the reservoir and Here, 25 • C is considered as a special case in Figure 8 in order to depict the I-V and P-V characteristics under standard test conditions (STCs), and to show the differences between the ISD, SD, and TD photovoltaic models. Moreover, Figure 8 shows that for a single PV panel, the power that can be extracted is 189.149 W in the case of the TD model. This power increases by 8.43% with the ISD model. So, the TD model gives more realistic and accurate results because it takes into account the realistic phenomenon of recombination losses. This will lead to a true size of the HRES, especially with the large-scale power system, whose load is shown in Figure 3.
In summary, many RES studies use the ISD model for PV system design [32,33]. However, the TD model (shown in Figure 6) is more accurate, because it takes into account both the diffusion and the recombination effects. This results in a more precise mathematical description of the solar panel, leading to a true size of the PV array and true reliability and environmental estimates of the entire system [17].

Modeling of Hydropower Facilities
There can be different topologies for a PHS system. However, the volume of the reservoir and the head elevation are the essential design parameters [40]. Other parameters of the PHS, such as the penstock and the selection of pumps and turbines, are more important in the mechanical design [41]. In this paper, the Al-Wehda impoundment was assumed to be the UR of the PHS energy system in this study. This dam is near Harta city, in northern Jordan. The actual head of this dam is 110 m. In this design, the volume of the lower reservoir (LR) of the PHS is one of the design variables to be optimally determined in order to truly size the hydro power plant facility. The two main operating modes to be considered in the energy management strategy and the mathematical modeling of the PHS system are the generating (discharging) and the pumping (charging) modes.

Generating (Discharging) Mode
In this paper, a PHS plant was operated based on the hourly measured volumetric flow rate . q T (t), as seen in Figure 2, so that the peak periods of load demand can be satisfied. This is due to the flexibility and the fast response of the PHS plant in serving the load and pumping. This PHS system is combined with a PV array to get a reliable resource of energy suppling the load demand of Harta shown in Figure 3. When the PV power is insufficiently satisfying the load demand, the discharging mode starts, and water will be released down from the upper reservoir to the turbine. The output from the hydro turbine generating unit can be expressed using Equations (5) and (6) [40]:

Pumping (Charging) Mode
This mode is in effect when there is surplus power from the PV array. The water will be pumped to the upper reservoir of Al-Wehda dam. Therefore, the flow rate of water extracted from LR in the charging mode can be expressed as shown in Equations (7) and (8) [6]:

Modeling of a PHS Plant
In order to model the PHS system, the quantity of water stored in LR must be known. The water quantity at time t is determined as shown in Equation (9) [6,8]. The water level in the LR can be modelled by the percentage of water in the LR (POW), which is between 0 and 1 and it is calculated using Equation (10): The volume of water of the PHS plant expressed in Equation (9) is subject to the constraints in Equation (11). The optimal value of V LR is also subject to the constraint POW as shown in Equation (12). Thus, the energy balance in Equation (20) will be optimally met in every time step [8]:

System Design Methodology
Numerous PV renewable energy studies and existing methodologies use only one solar model, which is the ISD or SD model in their studies. However, in our approach, designing the PV plant was done three times using three different PV models, ISD, SD, and TD. In this case, a comprehensive methodology was performed in this paper in comparison with other existing methodologies that use only one solar PV model. Furthermore, a detailed comparison in Section 4.1 for both PSO and WOA optimization algorithms for the three aforementioned PV models will be shown.

System Reliability
Reliability evaluation is considered as a good measure for the planning and operation stages of the power system, by precisely computing many reliability indices [42]. These performance indices evaluate the process of meeting the load demand for an HRES. Reliability is a major concern to meet the load demand over a given time period. This is due to the irregular nature of the power supply of RES, which significantly influences the HRES power generation and supply. This will also impact the load demand in standalone and on-grid HRES applications [5].
In many studies, the reliability of the system is taken as either a constraint or a goal to be achieved in the optimization process. In this paper, the value of the index of reliability (IR), shown in Equation (13) [5,29], will be optimized using the PSO algorithm and WOA to find the optimal size of the HRES configuration. Other reliability indices, such as loss of power supply probability (LPSP), loss of load probability (LOLP), and energy expected not served (EENS), will be computed corresponding to the optimal value of IR.
Here, LPSP is the probability, over a given time period, in which the power demand is not fully satisfied by the generated power. It is widely used in the literature to assess the system reliability due to the RES' intermittent nature. It is mathematically represented as shown in Equation (14) for a one-hour time step [29,30]. Furthermore, the power system reliability can be evaluated based on the LOLP index. The LOLP index is defined as the total failure hours of the energy deficit divided by the number of sample hours for the entire period, which is 8760 h for the year of 2018 in this paper [43]. This index is widely used in order to evaluate whether the system is able to serve the load demand, by indicating the percentage of time that the system is insufficient at supplying the load. This index can be computed as shown in Equation (15) [6,25]. Moreover, EENS is a probabilistic index for the reliability analysis of a power system. It represents the expected energy not supplied to the load due to the energy deficit as in Equation (16) [5,44]:

Ecological Impacts
The indirect and direct anthropogonic greenhouse gases (GHG) have been estimated to be 28,717 gigagrams in Jordan. A national inventory was performed by the Jordanian Ministry of Environment to estimate the GHG emissions in Jordan in 2006 [45].
Equation (17) is used to compute the total amount of GHG emissions. This is done for each kWh purchased from the grid and contains conventional brown energy sources [46]. Table 2 shows the constants used to compute emissions in Jordan: Another indicator for the GHG emission is the emissions reduction, which represents the difference between the baseline and the project emissions, as shown in Equations (18) and (19): The baseline emission is the amount of emissions generated by the grid-connected power plants when generating the same amount of power as the RES. The project emission is defined as the amount of emission generated by the RES, such as solar and hydro. This indicator is considered zero for the applicable RES projects.

System Operation and Optimization
The objective of the system optimization in this paper was to determine the optimal reliable system configuration. The objective function is the index of reliability (IR). Further, the decision variables are the number of PV modules (N PV ) and the volume of the lower reservoir (V LR ), with respect to technical constraints. Moreover, PSO and WOA are hourly testing of all possible configurations in order to meet the load consumption and determines the one that has the optimal value of IR. The energy balance Equation of the entire system is hourly tested and is expressed by Equation (20): P PV (t) N PV f PV η inv + P T (t) + P grid (Brown) (t) = P PVP (t) + P L (t) + P grid (sales) (t) (20) The system operation flow chart is shown in Figure 9. It describes the energy management strategy of the two modes of operation. In other words, it traces hydro operation with respect to PV generation as modelled by the energy balance Equation of Equation (20). Furthermore, Figure 9 helps researchers to simply follow the design methodology and the operation strategy modeling of the system described in Figure 1.
The system operation flow chart is shown in Figure 9. It describes the energy management strategy of the two modes of operation. In other words, it traces hydro operation with respect to PV generation as modelled by the energy balance equation of Equation (20). Furthermore, Figure 9 helps researchers to simply follow the design methodology and the operation strategy modeling of the system described in Figure 1.

Particle Swarm Optimization (PSO) Algorithm
The PSO algorithm is one of the popular nature-inspired metaheuristic optimization algorithms. It was developed by James Kennedy and Russell Eberhart in 1995 [47]. Moreover, it has been revitalized several times in order to solve practical issues related to engineering and power system optimization [48][49][50]. Recently, PSO has emerged as one of the promising algorithms in solving various optimization problems in the field of science and engineering as well [51][52][53][54][55].
As such, PSO is inspired by the social and cooperative behavior displayed by various species to look for their needs in the search space. The algorithm is guided by personal experience (P best ), overall experience (G best ), and the present movement of the particles to decide their next positions in the search space. Further, the experiences are accelerated by two factors c 1 and c 2 , and two random numbers generated between [0,1]. Whereas, the present movement is multiplied by an inertia factor, w, varying between [w min : w max ].
The initial population (swarm) of size N and dimension D is denoted by X = [X 1 , X 2 , . . . , X N ] T . Each individual (particle) is given by X i = [X i,1 , X i,2 , . . . , X i,D ]. Additionally, the initial velocity of the population is denoted by V = [V 1 , V 2 , . . . , V N ] T . The index i varies from 1 to N whereas the index j varies from 1 to D. The velocity and position of the particles are computed using Equations (21) and (22), respectively [56]:

Whale Optimization Algorithm (WOA)
The WOA is one of the nature-inspired meta-heuristic optimization algorithms. It simulates the social behavior of humpback whales. It imitates the bubble-net hunting mechanism of humpback whales. It was proposed by Mirjalili and Lewis in 2016 [57]. WOA is designed based on the spiral bubble-net feeding strategy, which is mathematically modelled in order to perform optimization.
Humpback whales recognize the prey's location and surround them. The WOA algorithm presume that the current supreme candidate solution is the targeted victim, or it is close to the optimal solution. After the supreme search is determined, the whales update their positions towards the supreme search. This behavior is given by Equation (23) and the position vector is computed as in Equation (24) [57]: where n is the iteration's number and X * is the position vector of the obtained best solution until now.
Note that this position is updated at each iteration if there is a superior solution. Furthermore, A and C are coefficients of the vectors and can be computed as in [57].

Results and Discussion
The mathematical models, design, simulation, methodology, and optimization algorithms were implemented to supply the load demand of Harta city, north of Irbid in Jordan. The average load was 4.067 MW/day, with a total energy of 38.3475 GWh in 2018. The hourly load values are shown in Figure 3.
Note that the results were obtained and computed using Matlab software. This was done after converting the energy management and operation strategy, shown in Figure 9, of our PV array combined with PHS grid-connected system into a Matlab code. Then, both the PSO algorithm and WOA were run with the inputs of the hourly measured data in 2018 shown earlier in Figures 2 and 3. Afterwards, the optimal value of the index of reliability (IR) was found for each PV model. Further, the corresponding detailed reliability results were computed at the optimal value of IR.

Impacts of PV Modeling
The variations of the objective function of IR as a percentage %, using PSO and WOA, are shown in Figures 10 and 11, respectively. Figure 10 shows that the PSO optimal values of IR are 98.553%, 98.527%, and 98.558% for the ISD, SD, and TD, respectively. Furthermore, the WOA optimal values of IR are 98.557%, 98.524%, and 98.565%, respectively, for the ISD, SD, and TD solar models. It can be noticed that the optimization results are approximately the same for the two algorithms. Further, the percentage difference between the PSO and WOA solutions is negligible. This means that the PSO and WOA are very close to the global optimum solution. Table 3 shows the corresponding decision variables of the system configuration at the optimal value of IR. The number of generations of PSO or WOA is 100 for all cases.
in Figures 10 and 11, respectively. Figure 10 shows that the PSO optimal values of IR are 98.553%, 98.527%, and 98.558% for the ISD, SD, and TD, respectively. Furthermore, the WOA optimal values of IR are 98.557%, 98.524%, and 98.565%, respectively, for the ISD, SD, and TD solar models. It can be noticed that the optimization results are approximately the same for the two algorithms. Further, the percentage difference between the PSO and WOA solutions is negligible. This means that the PSO and WOA are very close to the global optimum solution. Table 3 shows the corresponding decision variables of the system configuration at the optimal value of IR. The number of generations of PSO or WOA is 100 for all cases.   98.527%, and 98.558% for the ISD, SD, and TD, respectively. Furthermore, the WOA optimal values of IR are 98.557%, 98.524%, and 98.565%, respectively, for the ISD, SD, and TD solar models. It can be noticed that the optimization results are approximately the same for the two algorithms. Further, the percentage difference between the PSO and WOA solutions is negligible. This means that the PSO and WOA are very close to the global optimum solution. Table 3 shows the corresponding decision variables of the system configuration at the optimal value of IR. The number of generations of PSO or WOA is 100 for all cases. Figure 10. Progress charts of index of reliability (IR) in %, using particle swarm optimization (PSO) for the ideal single diode (ISD), single diode (SD), and two diode (TD) models. Figure 11. Progress charts of IR in %, using whale optimization algorithm (WOA) for the ISD, SD, and TD models. Figure 11. Progress charts of IR in %, using whale optimization algorithm (WOA) for the ISD, SD, and TD models. It can be seen that there is a slight difference in the value of IR obtained by PSO and WOA in the ISD, SD, and TD models. However, the TD model has the highest value of IR. This is because the TD model takes into account the recombination losses, which implies the need for more PV to satisfy the load demand. So, neglecting the recombination losses in both ISD and SD models will actually increase the number of solar panels, i.e., the size of the PV array. The percentage decrease of the number of PV panels is 16.67% and 7.93% for the ISD and SD models for the PSO algorithm, respectively, as shown in Table 3. Moreover, the corresponding percentage decrease in the case of WOA is 17.74% and 11.76% for the ISD and SD models, respectively. Tables 4 and 5 show the reliability indicators and energy results at the optimal value of IR. Table 4 shows that the most reliable model is the TD PV model in terms of EENS, LOLP, and LPSP. Table 5 shows that the PV energy, obtained using PSO, of the ISD and SD models decreased by 9.3% and 3.23%, respectively, from the TD realistic model. This means that the SD model is more accurate than the ISD model and close to the realistic TD model. Furthermore, the PV energy, obtained using WOA, of the ISD and SD models decreased by 10.47% and 7.26%, respectively, from the TD model. Thereby, the TD realistic model should be considered in large-scale RES solar applications. This is because it provides true and accurate size, energy, and reliability results. This will help the decision makers in their managerial environmental feasibility and reliability RES studies.
Here, the amount of GHG emissions and the corresponding emissions reduction were computed as shown in Table 6. The solutions obtained using PSO and WOA show that the TD solar model has the lowest GHG emissions and the highest ER compared with the ISD and SD solar models. For instance, it can be noticed regarding the PSO algorithm's solution that using the TD realistic solar model will mitigate the GHG emission by 1.72% and 7.2%, compared with the SD and ISD models, respectively. Further, the corresponding ER of the TD model increased by 2.07% and 6.34%, respectively, compared with both the SD and ISD solar models. This reinforces the recommendation that the TD model used in all PV-related studies will yield results in a more ecological RES.

Uncertainty Analysis
Uncertainty analysis measures the robustness of the system when one or more of the input data are varying. This includes different parameters, such as temperature, diode ideality factors, series, and leakage loss resistances. In this subsection, the uncertainty of the aforementioned parameters was tested for the TD realistic model as shown in Tables 7 and 8. This is to present the impact of increasing or decreasing these parameters on RES. In Section 4.1, it was concluded that the solutions of both PSO and WOA are close to each other. Therefore, uncertainty assessment will be performed using PSO as a case validation.  Table 7 shows that the diode ideality factor has a major impact on the design of the RES. When α for diode 2 is equal to 1.8, the PV panel percentage difference will be 4.63%. Moreover, the series loss influences the RES more than the R sh as shown in Table 8. This is because the second diode current, affected by the R s as shown in Equation (4), will make the loss much larger if this resistance becomes larger. Furthermore, when the temperature increases, the PV energy will slightly decrease and vice versa.
Moreover, Table 9 shows the percentage change of the TD model at α 2 = 1.5 and α 2 = 1.8 from the other two models. As can be seen, this slight change in α, due to recombination losses, has a large impact on the system sizing. The N PV and PV energy difference between the TD model at α 2 = 1.8 and the ISD model are 20.36% and 7.65%. Additionally, the GHG emissions will increase by 9.19%, when using the ISD solar PV model. Therefore, neglecting the recombination loss will achieve inaccurate sizing and under the estimation ecological estimates for the RES.

Conclusions
In this paper, the impact of different PV models in a hybrid power system was investigated. The models considered were the two-diode (TD), single-diode (SD), and the ideal single-diode (ISD) models of the solar PV renewable energy system (RES). This study was also performed in terms of the reliability and emissions performance indicators. The RES considered included a PV array combined with a pumped hydro storage (PHS), which was investigated for the Al-Wehda dam, Irbid, Jordan. This impoundment was treated as the upper reservoir (UR) of the PHS of the proposed RES. Precise mathematical modelling was done on the output PV power. This was done by taking into account the recombination losses in the TD PV model. This led to a more realistic sizing and more accurate system evaluation. Furthermore, the hourly measured values of the load demand, solar irradiance, ambient temperature, and volumetric water flow rate were obtained from related institutions in Jordan.
The results were obtained using the PSO and WOA optimization algorithms. The results showed that the TD realistic model is more reliable than the other models. The PSO optimal value of IR was 98.558%. The corresponding number of panels were 44,840 with a 65.052 M·m 3 volume of the lower reservoir. It was shown that the number of PV panels decreased by 16.67% and 7.93%, respectively, for the ISD and SD compared with the TD realistic model. These high percentage differences could lead to improper sizing and underestimation of the RES. Moreover, the WOA optimal value of IR was 98.565%. The percentage changes in N PV for WOA compared with PSO were 1.78%, 4.64%, and −0.5%, for the ISD, SD, and TD models, respectively.
Further, the ecological impacts of the proposed models were assessed. It showed that the TD model is more ecological, using the two optimization algorithms, compared with the other models, resulting in an emissions reduction of 21.519 and 21.529 Gg/year using PSO and WOA, respectively. These emissions reductions increased by 2.07% and 6.34%, respectively, for the SD and ISD solar PV models referring to the PSO solution. Furthermore, the GHG emissions were found to be 6.6879 and 6.64125 Gg/year in the case of the TD realistic model for PSO and WOA, respectively. These emissions increased for the SD and ISD models, respectively.
Finally, uncertainty analysis on the recombination losses, series and leakage losses, and temperature were performed in order to assess the robustness of the system using the PSO algorithm. This uncertainty analysis showed that the percentage difference of the number of PV panels for the realistic TD model decreased by 12% and 20.36% from the SD and ISD models, respectively, when the α 2 values were 1.5 and 1.8. Further, the PV energy of the ISD and SD models decreased by 10.47% and 4.48%, respectively, from the TD realistic model, when the α 2 was 1.5 and 1.8. This means that the SD model is more accurate than the ISD model and closer to the realistic TD model.

STC
Standard test conditions. f PV PV derating factor in %.

G(t)
Hourly solar irradiance in kW/m 2 . T amb (t) Hourly ambient temperature in • C.

I Ph
Photo current in A. I D1 /I D2 Shockley diode currents due to diffusion and recombination mechanisms in A.

I sh
Current flowing in the shunt resistance in A. I o1 /I o 2 Diodes reverse saturation currents in A.

I o STC
Diodes reverse saturation currents at STC in A. α 1 /α 2 Diodes ideality factors due to diffusion and recombination. V T Diode thermal voltages in V. R s /R sh Series and shunt resistances in Ω.

N ser
Number of PV modules connected in series per string.

N Par
Number of PV strings connected in parallel. P PV Rated PV power in W. P PV P Surplus PV power transferred to the hydro pumps. η inv Efficiency of the inverter . q T (t) Hourly hydro turbine flow rate in m 3 /s. h net Effective head altitude in m.

C T
Turbine generating coefficient in kg/m·s 2 . C P Turbine pumping coefficient in m 3 ·s/kg.

V LR
Volume of water stored in lower reservoir in m 3 .

POW
Percentage of water in the lower reservoir (LR) in %. α Evaporation and leakage loss in %. η T Overall efficiency of the turbine-generator unit in %. η P Overall pumping efficiency in %. ρ Density of water (1000 kg/m 3 ). g Acceleration due to gravity (9.8 m/s 2 ). q charge Quantity of water released in the charging mode in m 3 . q discharg Quantity of water stored in the discharging mode in m 3 .

E PVg
Annual PV energy generated in MWh.

E Tg
Annual hydro energy generated in MWh.

E grid sales
Annual energy sales to the grid in MWh.

E PVp
Annual pumped PV energy in MWh.

E brown
Annual energy purchased from the utility grid in MWh.