Simulation and Comparison of Mathematical Models of PV Cells with Growing Levels of Complexity

: The amount of energy generated from a photovoltaic installation depends mainly on two factors—the temperature and solar irradiance. Numerous maximum power point tracking (MPPT) techniques have been developed for photovoltaic systems. The challenge is what method to employ in order to obtain optimum operating points (voltage and current) automatically at the maximum photovoltaic output power in most conditions. This paper is focused on the structural analysis of mathematical models of PV cells with growing levels of complexity. The main objective is to simulate and compare the characteristic current-voltage (I-V) and power-voltage (P-V) curves of equivalent circuits of the ideal PV cell model and, with one and with two diodes, that is, equivalent circuits with ﬁve and seven parameters. The contribution of each parameter is analyzed in the particular context of a given model and then generalized through comparison to a more complex model. In this study the numerical simulation of the models is used intensively and extensively. The approach utilized to model the equivalent circuits permits an adequate simulation of the photovoltaic array systems by considering the compromise between the complexity and accuracy. By utilizing the Newton–Raphson method the studied models are then employed through the use of Matlab/Simulink. Finally, this study concludes with an analysis and comparison of the evolution of maximum power observed in the models.


Introduction
The Energy Union Framework Strategy is aiming to a serious transition from an economy dependent on fossil fuels to one more reliant on renewables [1] and among the available sources of renewable energy, solar energy is on the most abundant [2][3][4], which could be assertively harnessed, especially in the southern countries of Europe. According to Club of Rome study embracing the circular economy concept could signify up to 70% decrease in carbon emissions by 2030, of five European economies [5]. By targeting a renewable energy based economy and a circular economy at the same time could be the way achieve the Energy Union Framework Strategy targets [1]. Although free and available on a planetary scale, the role in the global energy mix is unobtrusive, competing not only with other forms of non-renewable energy, such as gas and coal [6], but also with its more direct rival-the wind renewable energy source [7]. Except for a very limited number of countries where proactive and generous income policies were implemented at the beginning of the last decade, there has been a more recent mobilization of European governments in this sector, legislating on specific instruments to stimulate production decentralized and small scale power [8]. In the last few years, the solar energy has been gaining importance in the worldwide energy evolution tendency due to a constantly increasing efficiency and lifespan, the decrease of the price of PV modules and by being environmentally friendly [9]. Solar photovoltaic (PV) systems are steadily becoming one of the main three electricity sources in Europe [10]. The entire installed PV capacity in 2016 reached 303 GWp and in that year Spain and Italy were responsible for 5.4 and 19.3 GWp, respectively [11].
In a typical PV system several photovoltaic modules are linked in series in order to create a PV string. The aim is to reach a certain voltage and power output. With the intention of accomplishing a greater power, such PV strings can be linked in parallel in order to make a PV array. For the duration of a constant irradiance condition, the power-voltage (P-V curve) characteristics of a PV string show a typical P-V curve peak. Such type of a peak embodies the maximum power of the PV string [4]. The P-V characteristics of a PV system are nonlinear and are affected by both the ambient temperature and solar irradiance, which in turn reveal distinct MPPs. With the purpose of optimizing the use of PV systems, conventional MPPT algorithms are often used [12].
In the planning of a photovoltaic power plant the electric power produced is strongly linked to the meteorological conditions (solar radiation and temperature) [13,14]. Due to the intermittent nature of solar energy, power forecasting is crucial for a correct interpretation of business profitability and payback time [15]. In the current market there is a great offer of manufacturers that, of course, have quite different technological production processes. All this leads to two modules with an identical technical sheet, under nominal test conditions, to differ in performance and produce very different results [16]. The actual operating conditions in both solar radiation and temperature will very rarely coincide with the combination of nominal meteorological variables. Thus, the broader characterization is of utmost importance for studying the differences [17]. In the end the main goal is to realistically quantify the performance, giving credibility to the estimation process in function of meteorological specificities of each season of the year. Ultimately it is desired that the process has enough resolution to reach the count up to the daily cycle.
The characterization requires the compilation of a large amount of data required for the application of appropriate mathematical model. In the concrete case of the photovoltaic cell the analytical model opens the doors for the detailed description in function of the external variables, which for all effects determine the general forms of the characteristic curves. However, to make modeling effective, it will be the model's intrinsic parameters which will more or less shape the link to the experimental data [18]. Molding is particularly critical at three points of operation. First, the predicted forecast of the peak electric power, then the open circuit operating points, that is, the maximum potential difference to be supported by the power electronics in the DC-AC conversation in the cut-off state, and finally the short-circuit, that is the maximum current to be supported by the electric cables in the event of a fault.
The models share in common the same electrical base model. The cell being a photoelectric device is modeled with a DC current source and a junction diode in parallel. From here all models are effectively variations with the introduction of more electrical elements. The elements may be of a series element of resistive nature by recreating the internal losses by Joule, or a parallel resistance simulating the internal leakage current, or a supplementary diode, which is normally associated with the losses by recombination of the carriers in the zone of the depletion layer [19].
Researchers have been increasingly focusing on MPPT techniques [20][21][22][23][24]. Authors in [25] have proposed a glowworm swarm optimization-based MPPT for PVs exposed to uneven temperature distribution and solar irradiation. A technique based on Radial Movement Optimization (RMO) for detecting the MPPT under partial shading conditions and then compared with the results of the particle swarm optimization (PSO) method is studied in [19]. Authors in [26] focus on the analysis of dynamic characteristic for solar arrays in series and MPPT based on optimal initial value incremental conductance strategy under partially shaded conditions. In [27] the authors optimize the MPPT with a model of a photovoltaic panel with two diodes in which the solution is implemented by Pattern Search Techniques. A PV source that was made by utilizing un-illuminated solar panels and a DC power supply that functions in current source mode is proposed in [18]. The authors in [28] address a simple genetic algorithm (GA)-based MPPT method and then compare the experimental and theoretical results with conventional methods. A direct and fully explicit method of extracting the solar cell parameters from the manufacturer datasheet is tested and presented in [29] and the authors base their method on analytical formulation which includes the use of the Lambert W-function with the aim of turning the series resistor equation explicit. The authors in propose a three-point weight method shared with fuzzy logic for increasing the speed of MPPT [30] and in this study the simulation was performed in Matlab and was experimentally validated.
The followed methodology was made for the comparison of the models in meteorological conditions as wide as possible. Extreme scenarios of incident solar radiation were simulated. The simulated temperature was considered suitable. The main goal of this study is to simulate and compare the characteristic curves of equivalent circuits of the ideal PV cell and, with one and with two diodes, respectively, namely equivalent circuits with five and seven parameters. The role of every parameter is assessed and compared. The ideal model of the PV cell is given in detail in [31]. The aim was to find areas of model intervention in which the modeling could lead to identical results. In this study the numerical simulation of the models is used intensively and extensively. The method used to model the equivalent circuits allows an adequate simulation of the photovoltaic array systems by taking into consideration the compromise between accuracy and complexity. By using the Newton-Raphson method the studied models are simulated through the use of Matlab/Simulink. All the simulations were carried out on the basis of a solar cell whose electrical specifications are given in [32].
The remainder of this paper is organized as follows. In Section 2 the equivalent circuit with five parameters is presented while in Section 3 the equivalent circuit with seven parameters is presented. The comparison between the one-diode model and the two-diode model is presented in Section 4. Finally, the conclusions are addressed in Section 5.

Representative Equations
The five-parameter circuit completes the frame of internal resistive losses. The fifth parameter corresponds to one more parasite resistance, referred to in this paper as the parallel resistance R p . Unlike the series resistance (R s ) it does not interfere directly with the power delivered to the load. However, it penalizes the operation of the cell by providing an alternative path for a portion of the photoelectric current. It is called a leakage current because it reduces the amount of current flowing at the PN junction [33], thereby affecting the voltage to the terminals of the photovoltaic cell. The five-parameter electrical circuit is the most widely used model in the analytical study of the photovoltaic cell. This model offers a good compromise in terms of complexity and performance [34], thus being the choice of several authors in this area of research [25,26,28,[35][36][37][38]. The model of the equivalent circuit with five parameters of the photovoltaic (PV) module can be observed in Figure 1. incremental conductance strategy under partially shaded conditions. In [27] the authors optimize the MPPT with a model of a photovoltaic panel with two diodes in which the solution is implemented by Pattern Search Techniques. A PV source that was made by utilizing un-illuminated solar panels and a DC power supply that functions in current source mode is proposed in [18]. The authors in [28] address a simple genetic algorithm (GA)-based MPPT method and then compare the experimental and theoretical results with conventional methods. A direct and fully explicit method of extracting the solar cell parameters from the manufacturer datasheet is tested and presented in [29] and the authors base their method on analytical formulation which includes the use of the Lambert W-function with the aim of turning the series resistor equation explicit. The authors in propose a three-point weight method shared with fuzzy logic for increasing the speed of MPPT [30] and in this study the simulation was performed in Matlab and was experimentally validated. The followed methodology was made for the comparison of the models in meteorological conditions as wide as possible. Extreme scenarios of incident solar radiation were simulated. The simulated temperature was considered suitable. The main goal of this study is to simulate and compare the characteristic curves of equivalent circuits of the ideal PV cell and, with one and with two diodes, respectively, namely equivalent circuits with five and seven parameters. The role of every parameter is assessed and compared. The ideal model of the PV cell is given in detail in [31]. The aim was to find areas of model intervention in which the modeling could lead to identical results. In this study the numerical simulation of the models is used intensively and extensively. The method used to model the equivalent circuits allows an adequate simulation of the photovoltaic array systems by taking into consideration the compromise between accuracy and complexity. By using the Newton-Raphson method the studied models are simulated through the use of Matlab/Simulink. All the simulations were carried out on the basis of a solar cell whose electrical specifications are given in [32].
The remainder of this paper is organized as follows. In Section 2 the equivalent circuit with five parameters is presented while in Section 3 the equivalent circuit with seven parameters is presented. The comparison between the one-diode model and the two-diode model is presented in Section 4. Finally, the conclusions are addressed in Section 5.

Representative Equations
The five-parameter circuit completes the frame of internal resistive losses. The fifth parameter corresponds to one more parasite resistance, referred to in this paper as the parallel resistance Rp. Unlike the series resistance (Rs) it does not interfere directly with the power delivered to the load. However, it penalizes the operation of the cell by providing an alternative path for a portion of the photoelectric current. It is called a leakage current because it reduces the amount of current flowing at the PN junction [33], thereby affecting the voltage to the terminals of the photovoltaic cell. The five-parameter electrical circuit is the most widely used model in the analytical study of the photovoltaic cell. This model offers a good compromise in terms of complexity and performance [34], thus being the choice of several authors in this area of research [25,26,28,[35][36][37][38]. The model of the equivalent circuit with five parameters of the photovoltaic (PV) module can be observed in Figure 1.  According to the junction or nodal rule the sum of currents is governed by the following condition: (1) and the voltage in the diode is equivalent to: Solving in order of I and replacing I D with the diode expression and I p with V d /R p the following equation is obtained: where I s is the current created by photoelectric effect, I is is the reverse saturation current, q is the charge of the electron, K is the Boltzmann constant (1.38 × 10 −23 J/ • K), T is the temperature of the junction, m is the reality parameter, R s is the parasite resistance in series and R p is the parallel parasite resistance. The value of R p is usually quite high in the manufactured photovoltaic cells. However, several authors with regard to this finding, consider useless the inclusion of this resistance [39][40][41][42][43][44]. On the other hand there are authors who consider R s negligible when the value is very low [45][46][47].
After obtaining the characteristic equation I-V the electrical power is calculated by: Deriving at the peak of power, one can find the voltage coordinate, as follows: The solution in order of V is only resolvable if applying an iterative numerical method.

Analytical Extraction of Parameters
Five equations are required. By consulting the manufacturer's information under nominal reference conditions the following equations are obtained: In practice the system is reduced to four algebraic equations. By observing Equation (3) it can be stated that: Energies 2018, 11, 2902

of 21
This means that it is possible to eliminate the term −1 without degrading the approximation given by the model to the I-V curve. This measure simplifies the analytical resolution of the four variables [48]. Thus, the system is limited to: The fourth equation is the expression of the power derivative in order of the voltage. The derivative can be decomposed as a function of V and I: which leads to: Since the Equation (3) is the type of I = f (I,V), the implied derivative as a function of I and V is: and dividing by dV it results in: By replacing Equation (18) in Equation (15) it is obtained: Solving the partial derivatives it is reached the explicit expression of the Equation (15): where the final presentation is: The system equations do not allow the separation of individual parameters I is , R s , R p , and m through the analytical solution. For this reason, appropriate numerical methods must be used.

Assessing Equations
The inverse saturation current is obtained by Equation (3) and it is referred to the open circuit operating point: In previous models only the equation of V in order to I required the Newton-Raphson method. If we try to derive the expression of V ca with Equation (3) set under open circuit conditions: The final result becomes: The assignment of one more parameter to the circuit structure renders impracticable the analytical resolution of the Equation (3). This means that it is not possible to separate and isolate the variables I and V in each member through elementary functions. Being the expression of the type I = f (I,V) the equation is commonly referred to as transcendental equation. In general, a transcendent equation does not have an exact solution [49]. The only way to find an approximate solution lies in the use of numerical calculation. In this context the Newton-Raphson algorithm was chosen. The Newton-Raphson method is a fairly fast (quadratic) convergence computational technique for calculating the roots of a function [50][51][52]. Due to its simplicity it lends itself perfectly to such problems. Then the Newton-Raphson method is used through its generic expression as shown in Equation (25): Being x n+1 the estimated value in the present iteration, x n the value obtained in the previous iteration, f (x) the function initialized with x n and the f'(x) the derivative initialized with x n . Accordingly, Equation (24) takes the form of a transcendental equation. Thus, by using the Newton-Raphson method through its generic expression (25) the voltage V ca can be assessed by: where I ca becomes a null value.
And the corresponding procedure for current I is: Knowing that the value of V is an input variable in the algorithm, then for each V there will be the corresponding I, computed iteratively by Equation (27). The convergence process ends when the following error criterion ε is satisfied: The nominal characteristic curves are obtained with Equations (4), (22), and (27). The remaining scenarios are supported by Equations (4), (26), (27), (29), and (30), where Equation (29) is a cubic relation between the inverse current and the temperature as proposed in [53,54]: where I is n is the inverse saturation current and T n is the temperature, both under Standard Test Conditions (STC) reference conditions. Additionally, in this study, the following simplification was taken into account, where G is the incident radiation in W/m 2 :

Comparison between Constant R s and Variable R p
Two comparative scenarios were designed for the characteristic curves at nominal reference conditions. In the first one the load is interconnected to a photovoltaic circuit dominated by resistive losses R p (R s = 0). The R p resistance was adjusted with 10 Ω, 200 Ω, and 1000 Ω, respectively. In the second scenario a fixed value of 10 mΩ was established for the R s resistor. The simulations can be observed in Figure 2. By observing the two graphs it is apparent that the resistance Rp does not interfere in the region of influence of the junction diode. In the region where the influence of the photoelectric current source predominates, the lowest value tested does not show a significant disturbance: the plot is very similar to the set of points estimated with Rp = ∞.
As the figures do not have sufficient detail, the curves were enlarged by a range of values close to the peak power. Figure 3 shows that the leakage current is virtually zero from 200 Ω. While for 10 Ω, the effect being visible is not at all significant. By observing the P-V curves in circuits with internal losses it can be verified that the lines are very similar, as can be in Figure 4. In other words, the leakage of current modelled by the resistance Rp in this range of values does not compromise the By observing the two graphs it is apparent that the resistance R p does not interfere in the region of influence of the junction diode. In the region where the influence of the photoelectric current source predominates, the lowest value tested does not show a significant disturbance: the plot is very similar to the set of points estimated with R p = ∞. As the figures do not have sufficient detail, the curves were enlarged by a range of values close to the peak power. Figure 3 shows that the leakage current is virtually zero from 200 Ω. While for 10 Ω, the effect being visible is not at all significant. By observing the P-V curves in circuits with internal losses it can be verified that the lines are very similar, as can be in Figure 4. In other words, the leakage of current modelled by the resistance R p in this range of values does not compromise the estimated maximum power. In this context of temperature and solar radiation this conclusion becomes valid. (STC).
By observing the two graphs it is apparent that the resistance Rp does not interfere in the region of influence of the junction diode. In the region where the influence of the photoelectric current source predominates, the lowest value tested does not show a significant disturbance: the plot is very similar to the set of points estimated with Rp = ∞.
As the figures do not have sufficient detail, the curves were enlarged by a range of values close to the peak power. Figure 3 shows that the leakage current is virtually zero from 200 Ω. While for 10 Ω, the effect being visible is not at all significant. By observing the P-V curves in circuits with internal losses it can be verified that the lines are very similar, as can be in Figure 4. In other words, the leakage of current modelled by the resistance Rp in this range of values does not compromise the estimated maximum power. In this context of temperature and solar radiation this conclusion becomes valid.

Characteristic Curves in Function of Temperature and Radiation
Using the same set of Rp resistors, three data set scenarios are established, as can be observed in Figures 5-7. Each scenario is simulated with a specific solar power, 100 W/m 2 , 500 W/m 2 , and 1000 W/m 2 , respectively, and having in common the same interval of test temperatures (10 °C, 25 °C, 50 °C, and 75 °C).

Characteristic Curves in Function of Temperature and Radiation
Using the same set of R p resistors, three data set scenarios are established, as can be observed in Figures 5-7. Each scenario is simulated with a specific solar power, 100 W/m 2 , 500 W/m 2 , and 1000 W/m 2 , respectively, and having in common the same interval of test temperatures (10 • C, 25 • C, 50 • C, and 75 • C). Using the same set of Rp resistors, three data set scenarios are established, as can be observed in Figures 5-7. Each scenario is simulated with a specific solar power, 100 W/m 2 , 500 W/m 2 , and 1000 W/m 2 , respectively, and having in common the same interval of test temperatures (10 °C, 25 °C, 50 °C, and 75 °C).      At the higher power range (1000 W/m 2 ) the I-V curves are literally identical. Consequently, the P-V curves do not experience significant changes between 10 • C and 75 • C. The peak power is then identical, regardless of whether the cell is manufactured with a high parallel resistance (1000 Ω) or with considerably low resistance (10 Ω).
At a medium power range (500 W/m 2 ) the performance is matched to that observed in the power ceiling of 1000 W/m 2 . With the solar radiation reduced to a tenth (100 W/m 2 ) of the highest power range, finally, there is some deviation in the I-V curve, characterized by R p = 10 Ω. In thermal terms, there is no correlation with R p : the difference with the versions with higher R p losses is apparently constant for the analyzed temperature scale.
Summarizing the Figures 4-6, it can be concluded that the impact of the parallel resistance on the performance of the cell, is barely expressive in the generality of the tested meteorological conditions, except for a slight disruption of the peak of power with 10 Ω in R p and under weak incident solar radiation.

Representative Equations
The seven-parameter electrical circuit is the next step in the electrical modeling of the photovoltaic cell. Equivalently to the mono-diode model (five parameters), the full version of two diodes brings together the complete set of losses. The seventh parameter is the leakage current modeled by the parallel resistance R p . The seven-parameter equivalent electrical circuit with two diodes can be observed in Figure 8. The seven-parameter electrical circuit is the next step in the electrical modeling of the photovoltaic cell. Equivalently to the mono-diode model (five parameters), the full version of two diodes brings together the complete set of losses. The seventh parameter is the leakage current modeled by the parallel resistance Rp. The seven-parameter equivalent electrical circuit with two diodes can be observed in Figure 8.
The voltage across the two diodes is equivalent to: Making the appropriate substitutions the final expression of I as a function of V is: where Is is the photoelectric equivalent current, Iis1 and Iis2 are the saturation currents of diode 1 and diode 2 respectively, m1 and m2 the ideality parameters of diode 1 and diode 2, respectively. As in the previous model, q represents the charge of the electron, K is the Boltzmann constant (1.38 × 10 −23 J/°K), T is the temperature of the junction, m is the reality parameter, Rs is the parasite resistance in series and Rp is the parallel parasite resistance. Despite being computationally more demanding, several authors argue that the approximation is more accurate than that achieved one with less complex models [45,[55][56][57][58]. For instance, for a low The sum of currents at the top node is: The voltage across the two diodes is equivalent to: Making the appropriate substitutions the final expression of I as a function of V is: where I s is the photoelectric equivalent current, I is1 and I is2 are the saturation currents of diode 1 and diode 2 respectively, m 1 and m 2 the ideality parameters of diode 1 and diode 2, respectively. As in the previous model, q represents the charge of the electron, K is the Boltzmann constant (1.38 × 10 −23 J/ • K), T is the temperature of the junction, m is the reality parameter, R s is the parasite resistance in series and R p is the parallel parasite resistance. Despite being computationally more demanding, several authors argue that the approximation is more accurate than that achieved one with less complex models [45,[55][56][57][58]. For instance, for a low radiation level, the two-diode model estimates with better approximation than the one-diode model [55,59].
Several authors simplify the identification of the parameters, reducing the number of effectively calculated variables. The most common is the reduction of seven to five variables by specifying fixed values. Usually this practice is related to the parameters m 1 and m 2 [55,[60][61][62]. Other authors opt for complete identification through elaborated methodologies such as particle examination optimization [63], the estimation based on neural networks [64], on genetic algorithms [65] or through algebraic relations as a function of temperature [62].
If the expression of I-V is identified the output electrical power P obeys to: From where by the derivative of the power peak it is possible to reach the value of V: Since the expression is transcendental the solution can only be found with a numerical algorithm that is able to extract the root.

Analytical Extraction of Parameters
Only six equations are required (the variable I s is excluded from the system since the linear dependence with temperature is known) the system is: As the role of parasite resistance R s is more pronounced in the vicinity of V ca , an orderly relation to this variable is determined through the derivative of the characteristic expression of the seven-parameter model: By rearranging this equation around R it becomes as follows: By replacing V with V ca and I with 0, the R s is as assessed as follows: where dV dI V ca = −R s is an initial estimation of the series resistance for the purposes of iterative numerical calculation.
By using Equation (40) de derived equation of R p is as follows: Meaning that in order of R p it becomes as follows: The R p in the vicinity of the short-circuit operating point is represented as follows: where dV dI I s = −R P is an approximate value of the parallel resistance for the purposes of iterative numerical calculation. It is usually estimated with the slope around the short-circuit operating point.
The still missing equation is the derivative of the power P as a function of V at the maximum electric power point. The developed equation takes the form of:

Assessing the Simulation Equations
The inverse saturation current is determined by Equation (33) at the open circuit operating point with the following analytical expression: The open circuit voltage V ca is given by: where I ca is the open circuit current. The current I is given by the following equation: The current and power strokes in nominal regime are estimated with Equations (34), (47), and (49). In a more comprehensive meteorological frame the calculation is carried out with the Equations (29), (30), (34), (48), and (49).

Characteristic Curves in Function of the Solar Radiation and the Parallel Resistance R p
The incorporation of the parallel resistance R p completes the number of variables that characterize the equivalent circuit of two diodes. Similarly to what was done with the equivalent representation of a diode, the importance of this resistive loss in the formation of the typical curves was examined, giving natural attention to the maximum power point. The structure was simulated with five different R p values, exposed to progressively higher levels of solar radiation, between the 100 W/m 2 and 1000 W/m 2 . Figures 9 and 10 show the generated curves.
The incorporation of the parallel resistance Rp completes the number of variables that characterize the equivalent circuit of two diodes. Similarly to what was done with the equivalent representation of a diode, the importance of this resistive loss in the formation of the typical curves was examined, giving natural attention to the maximum power point. The structure was simulated with five different Rp values, exposed to progressively higher levels of solar radiation, between the 100 W/m 2 and 1000 W/m 2 . Figures 9 and 10 show the generated curves.  In both a diode structure and a two-diode structure, the leakage current through Rp is very small if the resistance is simulated with 5000 Ω-from this value the parallel branch approaches an infinite resistance. Then, since the resistances of 10 Ω and 200 Ω lead to characteristic curves identical to those found with 5000 Ω, it can be stated that Rp is negligible if its value is equal to or greater than 10 Ω. The same is no longer true with Rp reduced to 1 Ω. The current I starts to decrease to values close to the short-circuit voltage instead of remaining constant until the measurements of the maximum power electrical coordinates.

Characteristic Curves as a Function of Temperature and Parallel Resistance Rp
In this section, the evolution of the curves as a function of Rp (10 Ω and 200 Ω) is analyzed. The series resistance is equal and constant in both models. The parameters m, m1, and m2 are initialized with 1.5, 1, and 2, respectively. The simulation was carried out with three scenarios of solar radiation with 100 W/m 2 , 500 W/m 2 and 1000 W/m 2 , respectively, and each with four levels of temperature, (10 °C, 25 °C, 50 °C, and 75 °C). The results can be witnessed in  With solar power at 1000 W/m 2 the open circuit voltage does not show any deviation between the two models. The same does not happen with 100 W/m 2 and 500 W/m 2 aggravating with the decrease of incident solar exposure. As for maximum power, the two-diode model is always at an advantage whatever the scenario. The difference is visible between 100 W/m 2 and 500 W/m 2 , with a tendency to increase in the downward direction of the sun exposure. Over the same range of solar radiation, the temperature tends to maintain the constant difference. In both a diode structure and a two-diode structure, the leakage current through R p is very small if the resistance is simulated with 5000 Ω-from this value the parallel branch approaches an infinite resistance. Then, since the resistances of 10 Ω and 200 Ω lead to characteristic curves identical to those found with 5000 Ω, it can be stated that R p is negligible if its value is equal to or greater than 10 Ω. The same is no longer true with R p reduced to 1 Ω. The current I starts to decrease to values close to the short-circuit voltage instead of remaining constant until the measurements of the maximum power electrical coordinates.

Characteristic Curves as a Function of Temperature and Parallel Resistance R p
In this section, the evolution of the curves as a function of R p (10 Ω and 200 Ω) is analyzed. The series resistance is equal and constant in both models. The parameters m, m 1 , and m 2 are initialized with 1.5, 1, and 2, respectively. The simulation was carried out with three scenarios of solar radiation with 100 W/m 2 , 500 W/m 2 and 1000 W/m 2 , respectively, and each with four levels of temperature, (10 • C, 25 • C, 50 • C, and 75 • C). The results can be witnessed in Figures 11-13.

Comparative Table of Peak Power in the Set of Models
In order to support the conclusions concerning the role played by the parallel resistance in the one diode and two diode models, the data referring to the maximum power were agglutinated. The results are organized according to three solar power levels and by test temperature families and can be observed in Figures 14-16.  With solar power at 1000 W/m 2 the open circuit voltage does not show any deviation between the two models. The same does not happen with 100 W/m 2 and 500 W/m 2 aggravating with the decrease of incident solar exposure. As for maximum power, the two-diode model is always at an advantage whatever the scenario. The difference is visible between 100 W/m 2 and 500 W/m 2 , with a tendency to increase in the downward direction of the sun exposure. Over the same range of solar radiation, the temperature tends to maintain the constant difference.

Comparative Table of Peak Power in the Set of Models
In order to support the conclusions concerning the role played by the parallel resistance in the one diode and two diode models, the data referring to the maximum power were agglutinated. The results are organized according to three solar power levels and by test temperature families and can be observed in Figures 14-16.

Comparative Table of Peak Power in the Set of Models
In order to support the conclusions concerning the role played by the parallel resistance in the one diode and two diode models, the data referring to the maximum power were agglutinated. The results are organized according to three solar power levels and by test temperature families and can be observed in Figures 14-16.    From the analysis of Figures 14-16 several conclusions can be assessed. At maximum or half solar exposure, both models lead to maximum values very close to the estimated power at infinite Rp. This trend remains with the temperature varying between 10 °C and 75 °C. Decreasing the light exposure to one tenth shows a significant deviation in any of the models with the Rp reduced to 10 Ω. The two-diode equivalent circuit is in any case more generous at peak power. The power deficit between the one-diode model versus the two-diode model is constant over the entire temperature range for the same level of radiation. However, it tends to worsen with the weakening of the incidence variable. The two-diode model tends to approximate the ideal one-diode model with the progressive reduction of incident radiation. With the limited incidence at 100 W/m 2 the equivalent circuit performance (Rs = 10 mΩ and Rp = 200 Ω) is comparable, as it can be observed in Figure 16.

Conclusions
In this paper the equivalent electrical circuits used in the modeling of non-organic photovoltaic cells was presented, paying particular attention to the modeling of silicon-made cells. Two equivalent circuits of models were analyzed and then compared with the ideal model of the PV cell. The first equivalent circuit consists of the model of one diode, the second equivalent circuit of two diodes. The results show that the two-diode equivalent circuit is more advanced than the diode circuit in modeling an internal leakage current. The second diode fulfills this role by describing the additional losses, associated with the recombination of carriers in the depletion layer. The results firmly reveal that the five-parameter model is more penalized with the decrease in radiation than the seven-parameter counterpart. The same trend was observed with the rise in temperature. With a meteorological frame characterized by 1000 W/m 2 of radiation and 10 °C of temperature, the  From the analysis of Figures 14-16 several conclusions can be assessed. At maximum or half solar exposure, both models lead to maximum values very close to the estimated power at infinite Rp. This trend remains with the temperature varying between 10 °C and 75 °C. Decreasing the light exposure to one tenth shows a significant deviation in any of the models with the Rp reduced to 10 Ω. The two-diode equivalent circuit is in any case more generous at peak power. The power deficit between the one-diode model versus the two-diode model is constant over the entire temperature range for the same level of radiation. However, it tends to worsen with the weakening of the incidence variable. The two-diode model tends to approximate the ideal one-diode model with the progressive reduction of incident radiation. With the limited incidence at 100 W/m 2 the equivalent circuit performance (Rs = 10 mΩ and Rp = 200 Ω) is comparable, as it can be observed in Figure 16.

Conclusions
In this paper the equivalent electrical circuits used in the modeling of non-organic photovoltaic cells was presented, paying particular attention to the modeling of silicon-made cells. Two equivalent circuits of models were analyzed and then compared with the ideal model of the PV cell. The first equivalent circuit consists of the model of one diode, the second equivalent circuit of two diodes. The results show that the two-diode equivalent circuit is more advanced than the diode circuit in modeling an internal leakage current. The second diode fulfills this role by describing the additional losses, associated with the recombination of carriers in the depletion layer. The results firmly reveal that the five-parameter model is more penalized with the decrease in radiation than the seven-parameter counterpart. The same trend was observed with the rise in temperature. With a meteorological frame characterized by 1000 W/m 2 of radiation and 10 °C of temperature, the From the analysis of Figures 14-16 several conclusions can be assessed. At maximum or half solar exposure, both models lead to maximum values very close to the estimated power at infinite R p . This trend remains with the temperature varying between 10 • C and 75 • C. Decreasing the light exposure to one tenth shows a significant deviation in any of the models with the R p reduced to 10 Ω. The two-diode equivalent circuit is in any case more generous at peak power. The power deficit between the one-diode model versus the two-diode model is constant over the entire temperature range for the same level of radiation. However, it tends to worsen with the weakening of the incidence variable. The two-diode model tends to approximate the ideal one-diode model with the progressive reduction of incident radiation. With the limited incidence at 100 W/m 2 the equivalent circuit performance (R s = 10 mΩ and R p = 200 Ω) is comparable, as it can be observed in Figure 16.

Conclusions
In this paper the equivalent electrical circuits used in the modeling of non-organic photovoltaic cells was presented, paying particular attention to the modeling of silicon-made cells. Two equivalent circuits of models were analyzed and then compared with the ideal model of the PV cell. The first equivalent circuit consists of the model of one diode, the second equivalent circuit of two diodes. The results show that the two-diode equivalent circuit is more advanced than the diode circuit in modeling an internal leakage current. The second diode fulfills this role by describing the additional losses, associated with the recombination of carriers in the depletion layer. The results firmly reveal that the five-parameter model is more penalized with the decrease in radiation than the seven-parameter counterpart. The same trend was observed with the rise in temperature. With a meteorological frame characterized by 1000 W/m 2 of radiation and 10 • C of temperature, the deviation approaches 6.1%. By reducing the exposure to one-tenth the deviation reaches 12%. At the other end of the temperature range, the deviation reaches 10.4% at full sun exposure, and worsens up to 20.3% with exposure limited to the maximum. The most important contribution deduced form this study is that the two-diode model tends to approximate to the ideal PV cell model (one-diode model) with the progressive reduction of incident radiation. With the incidence limited to 100 W/m 2 the equivalent circuit performance (R s = 10 mΩ and R p = 200 Ω) is almost identical to the ideal one-diode model. This means that, for regions were the solar incident radiation is lower, the ideal one-diode model behaves similarly to the more complex seven parameter equivalent circuit, thus allowing the user to opt for this circuit in detriment to the other more complex one which allows using a less complex software tool.