# Accurate Simulation of MPPT Methods Performance When Applied to Commercial Photovoltaic Panels

Authors: Javier Cubas; Santiago Pindado; Ángel Sanz-Andrés
Journal: The Scientific World Journal
 (2015)

---

## Abstract

A new, simple, and quick-calculation methodology to obtain a solar panel model, based on the manufacturers’ datasheet, to perform MPPT simulations, is described. The method takes into account variations on the ambient conditions (sun irradiation and solar cells temperature) and allows fast MPPT methods comparison or their performance prediction when applied to a particular solar panel. The feasibility of the described methodology is checked with four different MPPT methods applied to a commercial solar panel, within a day, and under realistic ambient conditions.


---

## Body

## 1. Introduction

Oil prices soared in the second half of the 20th century, the main causes being sociopolitical instabilities such as revolutions or wars, shortage of supply, stock market speculation, or the increasing demand from emerging nations [1]. The economic problems derived from the unstable situation of the energy market forced both the development of renewable energy sources (photovoltaic energy, solar thermal energy, and wind energy) [2–4] and the increase in efficiency in terms of energy generation, transportation, and consumption. In addition, the demand for more efficient energy supply systems has been increased as societies all over the world are increasingly aware of the consequences of the global climate change produced by greenhouse gas emissions.Among the different renewable energy sources that have been developed and spread in the past half century, the photovoltaic energy is a good example of the above two tendencies. On the one hand, the installed power has increased enormously in some parts of the world, as it has been mainly supported by government investment [5, 6]. On the other hand, technology improvements have produced more efficient solar cells in terms of energy (from 11% efficiency of silicon solar cells in the 1950s [7], so that now it is possible to get gallium-arsenide solar cells of above 30% efficiency [8–11]), cheaper cells based on organic technologies [12–14], and better control processes to maximize the energy supply in every irradiation and temperature conditions.It is well known that photovoltaic systems are affected by factors that reduce their efficiency such as(i)
changes on irradiation,(ii)
changes on cells temperature,(iii)
impedance variations at the system output,(iv)
partial shading on the photovoltaic panel.All these factors change the behavior of the panel, which is normally defined by the output current-output voltage curve (hereinafter, theI-V curve); see Figure 2 further in the text, and, consequently, the maximum power point (MPP) of the panel is also modified. Therefore, if a photovoltaic system needs to be optimized in terms of power production, implementation on the system of a methodology to “follow” the changes of the MPP is required.To the authors’ knowledge, the first works to analyze and track the maximum power point of photovoltaic systems were carried out in the late 1960s [15]. Since then, the number of papers produced has broadly increased in line with the oil prices evolution. See in Figure 1 the number of papers related to photovoltaic maximum power point tracking (MPPT) methods that, according to Esram and Chapman, were published between 1968 and 2005, together with the oil prices per barrel (Brent and West Texas Intermediate). In addition, it can be noted that the number of papers per year regarding MPPT in photovoltaic systems has been continually rising as, taking Google Scholar as reference, it has increased from 252 in 2006 to 3220 in 2013. These figures point out the present relevance of improving the MPPT methods.Figure 1
Left axis: oil prices per barrel (Brent and West Texas Intermediate). Right axis: number of papers published between 1968 and August 2005 according to Esram and Chapman.Figure 2
I-V curve of a solar panel. The three characteristic points (short-circuit, maximum power, and open-circuit points) are indicated on the curve.Papers regarding MPPT methods are not only numerous but diverse, an interesting classification of them being shown by Reza Reisi et al. [16]. These authors divide the different MPPT methods into three categories as follows:(i)
offline methods, which are those that require information regarding the panelI-V curve, and also temperature and irradiation levels; among these methods it is possible to find
(a)
open-circuit voltage and short-circuit current methods, which define the MPP, respectively, as a fraction of the panel open-circuit voltage,Voc, or the short-circuit current, Isc;(b)
artificial intelligence methods, such as those based on neural networks or fuzzy logic;(ii)
online methods, which are those that require instantaneous measurements of the photovoltaic panel output current and output voltage; the main advantage of these methods is that no information on the panelI-V curve or regarding irradiation or temperature levels is required; some examples of these methods could be
(a)
perturbation and observation methods;(b)
incremental conductance methods;(c)
power peak seeking methods;(d)
incremental conductance method;(e)
ripple correlation control methods;(iii)
hybrid methods, which are those that combine two different methods, each one from one of the above categories; the offline method is used to get a quick approximation to the MPP whereas the online method is used to improve the result.Nevertheless, it should also be pointed out that this classification does not include all possible methods. As an example, the power curve fitting method described by Salas et al. and other authors [18, 19] requires an accurate knowledge of both the panel behavior (i.e., theI-V curve) and the current and voltage measurements at several working points. A more extensive description and classification of MPPT methods can be found at [15, 19–22].As expected, the high scientific production regarding MPPT methods makes periodical benchmarking necessary, in order to be able to select the best method for every photovoltaic system. Since the 1990s several works aiming at comparing the different MPPT methods can be found at the available literature [15, 16, 18–28]. However, the conclusions from those works should not be widely accepted, as the possible combinations of external factors like the irradiance and temperature variation speed with internal factors such as the system design are an unreliable figure. On the other hand, the specific use/purpose of the generated photovoltaic energy should be taken into account when choosing among the different possible MPPT methods for the considered power system. In line with this statement, one can distinguish between the massive energy supply from solar parks (grid-connected systems) and the supply for small autonomous (stand-alone) systems such as satellites and other spacecraft (although some of them can hardly be considered as small, that is, the International Space Station). In this sense, the differences of MPPT methods applied to mass energy generation or grid-connected systems and to power generation on spacecraft can be underlined. The grid-connected photovoltaic systems need to maximize the production at the lower cost impact, taking into account the fact that a small increase in efficiency can be translated into a huge growth of revenues after some years of operation. This is achieved through MPPT methods with high dynamic and tracking capacities. On the other hand, it is also true that spacecraft require the best possible energy efficiency, but only after considering the reliability of the system [29] and some other associated problems, such as the thermal control [30]. Therefore, in stand-alone systems, especially in spacecraft, simple but more reliable MPPT methods are normally chosen, thereby excluding more efficient, but also more complicated methods.Validation or comparison of MPPT methods is usually carried out either based on experimental testing [31, 32] results or based on computer simulations [25, 33, 34] although sometimes both kinds of procedures are combined [23, 24]. Experimental comparisons of MPPT methods are logically more realistic, but it should be pointed out that simulations have some important advantages in terms of fast results, cost, and versatility, which make them the best option for algorithms development or efficiency estimations.One of the most significant difficulties of computational analysis lies in the solar panel behavior simulation. The most common way to do it is through equivalent circuit models, as this has proven to be very accurate [35]. Nevertheless, it must be underlined that different sun irradiation and temperature levels on the panel, as is required by any MPPT validation, involve the recalculation of the different parameters which define the aforementioned equivalent circuit. In many cases, some simplifications of the equivalent circuit are considered, taking into account the fact that multiple recalculations are a quite complicated task. An example of these simplifications could be a circuit without series resistor [34], without parallel resistor [23], or without any resistor [25]. However, series and shunt resistances are directly related to the slope of theI-V curve at open-circuit and short-circuit points, respectively (dI/dVsc≈-Rs-1 and dI/dVoc≈-Rsh-1), and not taking them into account modifies the shape of theI-V curve and, therefore, the position of the MPP. These simplifications make the behavior of the MPPT numerical simulation slightly diverge from the real performance of the photovoltaic system.In the present paper a simple, direct, and analytical way to define the behavior of a photovoltaic system under any sun irradiation and cells temperature conditions is presented. Based on the proposed method, solar panel models that fit perfectly to manufacturers’ datasheet experimental data can be defined. The equivalent circuit model proposed consists of one diode and two resistors (one in series and the other one in parallel), with no further simplification considered.Through the methodology presented in this work, accurate simulations of MPPT strategies/algorithms can be achieved, bearing in mind different irradiance and solar cell temperature level and also different DC-DC convertors to be connected to the panel. The aforementioned methodology is based on the 5-parameter equivalent circuit model of a photovoltaic device (i.e., solar panel), the parameter-extraction being analytical (described in Section2). This analytical extraction to fit the equivalent circuit behavior of theI-V curves corresponding to solar cell has proven to be accurate in previous works, in which different photovoltaic technologies (monocrystalline/polycrystalline) [36] or materials (Si/GaAs) [35] were analyzed. Also, if the equations are adjusted to reflect the number of cells connected in series and in parallel, this model can be used to analyze multijunction cells, solar panels, or even groups of solar panels [37]. Nevertheless, it should also be mentioned that the low number of theI-V curve points used for extrapolating the curve itself is the greater limitation of the analytical methods, as these points need to be very accurately measured (any deviation at those points will be extrapolated along theI-V curve). On the other hand, some optimization of the analytical results can be introduced if some extra data is available [37].The text is organized as follows. In Section2 solar panel analytical modeling is described. In Section 3 the way to take into account ambient condition variations in the modeling is shown. In Section 4 four different MPPT methods are presented, whereas calculations and results regarding the MPPT methods comparison carried out with the proposed methodology are included in Section 5. Finally, a case study is included in Section 6, whereas the conclusions of the present work are summarized in Section 7.


## 2. Solar Panel Modeling

In Figure2, a typical solar cell/panelI-V curve is shown. The characteristic points, short-circuit [Isc,0], open-circuit [0,Voc], and maximum power [Imp,Vmp] points are indicated in the graph. This behavior is similar to the one from a circuit formed by a source of current, Ipv, in parallel to a diode (see Figure 3). Two resistors, one in parallel (shunt resistor), Rsh, and the other one in series (series resistor), Rs, are added to take into account losses, such as the ones produced in cell solder bonds, interconnection, and junction box, together with current leakage through the high conductivity shunts across the p-n junction [35].Figure 3
Typical 1-diode equivalent circuit of a solar panel.The equation that defines the behavior of the 1-diode/2-resistors equivalent circuit of Figure3 is(1)I=Ipv-I0exp⁡V+IRsaVT-1-V+IRsRsh,where, together with the above explained variables (Ipv, Rsh, and Rs), I0 and a are, respectively, the reverse saturation current and the ideality factor that takes into account the deviation of the diodes from the Shockley diffusion theory. In accordance with Jain and Kapoor (2005), several theories have been developed to explain the effects that are taken into account by the ideality factor. It is generally accepted that a value from a=1 and a=2 takes into account recombination effects inside the cell, whereas larger values take into account shunt resistance effects and nonuniformities on the aforementioned recombination [38]. VT is not an unknown parameter; it is the thermal voltage of the diode and depends on the charge of the electron, q; the Boltzmann constant, k; the number of cells in series, n; and the temperature, T:(2)VT=nkTq.The above expression (1) leaves 5 parameters to be adjusted. This parameter identification is normally performed by fitting the equivalent circuit performance to anI-V curve obtained for certain working conditions (sun irradiance and cells temperature). However, on many occasions the aformentionedI-V curve is not available, only information from the manufacturer’s datasheet being accesible. In this case, analytical models are probably the best way to approach the equivalent circuit parameter-extraction problem [39–41]. Making use of the characteristic points of theI-V curve (see Figure 2), the following expressions can be derived from (1) [35]:(i)
short-circuit equation:(3)Isc=Ipv-I0exp⁡IscRsaVT-1-IscRsRsh,
(ii)
open-circuit equation:(4)0=Ipv-I0exp⁡VocaVT-1-VocRsh,
(iii)
maximum power point equation:(5)Imp=Ipv-I0exp⁡Vmp+ImpRsaVT-1-Vmp+ImpRsRsh,
(iv)
zero derivative for the power at maximum power point equation:(6)-ImpVmp=-I0aVT1-ImpVmpRsexp⁡Vmp+ImpRsaVT-1Rsh1-ImpVmpRs.As four equations are defined for five parameters to be identified, some initial estimation with regard to any of those parameters or to any particular characteristic of theI-V curve is required. In the first case, calculations are generally started with an estimation of the ideality factor, a (which normally lies in the bracket [1,1.5] [39, 42]), as the error of using this starting estimation is reduced (which only affects the curvature of the curve around the maximum power point). Besides, some works in the open literature [43] give information to make such estimations depending on the solar cell technology (Si, Ga-AS, etc.), and, more importantly, this ideality factor can be used as an iteration parameter to obtain an improved approximation to the equivalent circuit [37]. In the second case, some authors use the slope of theI-V curve at short-circuit or open-circuit points (measured or estimated) to obtain additional boundary conditions equations [35, 41]. Therefore, once the ideality factor a has been estimated, the following expressions for the equivalent circuit parameters can be obtained, from (3), (4), (5), and (6), by using the Lambert W-function, W(z), where z=WzeW(z) [36]:(7)Rs=AW-1BeC-D+C,A=aVTImp,B=-Vmp2Imp-IscVmpIsc+VocImp-Isc,C=-2Vmp-VocaVT+VmpIsc-VocImpVmpIsc+VocImp-Isc,D=Vmp-VocaVT,Rsh=Vmp-ImpRsVmp-RsIsc-Imp-aVTVmp-ImpRsIsc-Imp-aVTImp,I0=Rsh+RsIsc-VocRshexp⁡Voc/aVT,Ipv=Rsh+RsRshIsc.The equivalent circuit obtained with the above equations can reproduce the consideredI-V curve. However, this curve changes as a function of the sun irradiation and the solar cells temperature, as previously stated. An increase of the temperature produces a higher short-circuit current, Isc, lower values of the open-circuit voltage, Voc, and the maximum available power, Pmp=ImpVmp. These variations of the characteristic points are approximately linear in relation to the temperature, T:(8)Voc,T=Voc,Tr1+βVocT-Tr100,Vmp,T=Vmp,Tr1+βVmpT-Tr100,llIsc,T=Isc,Tr1+αIscT-Tr100,lImp,T=Imp,Tr100+γT-Tr100+βVmpT-Tr,where Tr is the reference temperature; βVoc and βVmp are, respectively, the percentage variation of the open-circuit and maximum power point voltages when the temperature increases one degree; αIsc and αImp are the percentage variation of the short-circuit and maximum power point currents when the temperature increases one degree; finally, γ is the percentage variation of the maximum power with temperature.As aforementioned, sun irradiation variations also modify theI-V curve. However, manufacturers normally do not include any information regarding these variations in the solar panel datasheets. Commonly, the shape of theI-V curve is considered essentially invariant with irradiation levels within ranges around one solar constant, so this leads to the following equation, considering, respectively, linear and exponential variations of the short-circuit current, Isc, and the open-circuit voltage, Voc, with temperature, whereas Rs remains unaffected for temperature variations [44]. Those conditions lead to the following equation [39]:(9)Ipv,G=Ipv,GrGGr.In the above equation G is the irradiance on the cell/solar panel, Ipv,G is the photocurrent delivered by the current source of the equivalent circuit, and Gr and Ipv,Gr are the reference values.Taking all the above statements into account, a simple but accurate (i.e., strictly respecting the data from the manufacturer’s datasheet) way to model a solar panel for any irradiation and temperature levels could be summarized as follows.(i)
Estimate the ideality factora.(ii)
Establish the temperature range in which the solar panel behavior should be modeled and calculate theI-V curve characteristic points in that range by making use of (8).(iii)
Obtain from expressions (7) the equivalent circuit parameters within the selected temperature range.(iv)
Fit a polynomial expression to the variations of the equivalent circuit parameters as a function of temperature.(v)
Introduce the effect of the irradiance in the parameterIpv by using expression (9).The above explained process has been applied to the YL280C-30b solar panel manufactured by Yingli Solar (Baoding, China), which will be used in the MPPT simulations included in Section 5. Data regarding the characteristic points of theI-V curve from the manufacturer’s datasheet are included in Table 1. The temperature range considered for the simulations is from 10°C to 65°C.Table 1
Characteristic points of YL280C-30b solar panel (Yingli Solar), included in the manufacturer’s datasheets[ref] at Standard Test Conditions (STC): 1000 W/m2 irradiance, 25°C cell temperature, and AM1.5 g spectrum according to EN 60904-3 [17].
YL280C-30b (monocrystalline)
n
60
T
r (°C)
25
P
mp (W)
280
γ (%/°C)
−0.42
I
mp (A)
8.96
α
I
mp (%/°C)
—
V
mp (V)
31.3
β
V
mp (%/°C)
−0.41
I
sc (A)
9.50
α
I
sc (%/°C)
0.04
V
oc (V)
39.1
β
V
oc (mV/°C)
−0.31With regard to the ideality factor some additional considerations have been made. Although taking into account the silicon monocrystalline technology a valuea=1.2 should be selected [43], a lower value has finally been chosen, a=1.05, as this parameter decreases with temperature and the solar panel is supposed to operate at higher temperatures than the reference one, T=25°C, from the Standard Test Conditions (STC). After choosing the value for the ideality factor, the other parameters have been calculated using (7); see Figure 4. The polynomial fittings to the data have also been included in the graphs of the figure, the irradiance level, G, being considered in the case of the photocurrent, Ipv; see the following:(10)IpvT,G=9.50+3.89·10-3ΔTlllllllllllllllllll-1.09·10-6ΔT2-2.96·10-8ΔT3GGr,llllRsT=3.44·10-1+4.23·10-4ΔTllllllllllllllllll+5.90·10-6ΔT2+1.93·10-8ΔT3,lliRshT=8.83·10-4+2.57·10-5ΔTlllllllllllllllllllll-4.02·10-7ΔT2-7.56·10-9ΔT3-1,lllllI0T=exp⁡-2.19·101+1.56·10-1ΔTlllllllllllllllllllllllll-5.20·10-4ΔT2+1.52·10-6ΔT3.When evaluating an MPPT method the maximum possible power that could be extracted from the panel pmax⁡t=ImptVmpt has to be calculated in every instant, t. Then, the efficiency of the method can be estimated with the following expression [45]:(11)ηMPPT=∫0TTpMPPTtdτ∫0TTpmax⁡tdτ,where pMPPTt is the instantaneous power obtained from the panel using the selected MPPT method and TT is the total period of time in which the aforementioned MPPT method is evaluated.Figure 4
Calculated values of the equivalent circuit parametersRs, Rsh, Ipv, and I0, for YL280C-30b monocrystalline panel as a function of the temperature, T. Calculated points are indicated with symbols whereas the polynomial approximations fitted to those data (10) have been included in each case as solid lines.
(a)
(b)
(c)
(d)In case of evaluation through experimental procedures, the maximum possible power,pmax⁡t, cannot be directly obtained from the solar panel, being instead estimated from the sun irradiance and temperature levels combined with the panel efficiency at STC. On the other hand, in case of analytical modeling, voltage sampling can be applied to (1) in order to find the highest value of the productV·I [39], or the maximum power point variables, Imp and Vmp, can be obtained from expressions (5) and (6). As these expressions are implicit and coupled, their resolution needs to be done by using numerical calculation or through simplifications in order to reduce the complexity of the calculations [46]. To the authors’ knowledge, no method to uncouple variables Imp and Vmp from (5) and (6) has been proposed in the available literature. Nevertheless, this possibility does indeed exist by making the following changes of variable:(12)φ=-ImpVmp,θ=Vmp+RsImp.Then, from (5) and (6), it is possible to obtain uncoupled equations as a function of the above new variables:(13)1+φRs+Rsh1+φRsaVTRsh+I0exp⁡Rsφ-1Rs+Rshφ-1lllllllllllllllllll×1+φRs+Rsh1+φRs+Ipv+I0RshaVT=0,(14)θ=aVTRsφ-1Rs+Rshφ-1×1+φRs+Rsh1+φRs+Ipv+I0RshaVT.The variable φ can be calculated solving the implicit (13), and then θ can be obtained directly from (14). Then, unmaking the change of variables, Imp and Vmp, can be derived after only one implicit equation resolution:(15)Imp=φθRsφ-1,(16)Vmp=θ1-Rsφ.This procedure gives the values of current and voltage at the MPP in a very easy way, avoiding voltage sampling or simplifications that could reduce the accuracy of the results. In order to speed up the numerical resolution of (13) it might be useful to start it with the initial value of φ0=-Isc/Voc.


## 3. Ambient Conditions

Two different ambient conditions (sun irradiance and temperature levels) have been considered for the simulations carried out. The first one is a cloudy day, with low but very unstable irradiance level, and the second one is a sunny day; see graphs in Figure5. These data are, respectively, the irradiance level at the Goddard Space Flight Center (GSFC) on May 13 (cloudy day) and May 14, 1971 (sunny day), and have been extracted from the work by Thekaekara [47]. Temperature and wind speed conditions corresponding to those days have been obtained from the measurements done at Washington, DC, this city being quite close to GSFC [48].Figure 5
Ambient conditions on May 13, 1971 (a), (b), and May 14, 1971 (c), (d). Solar irradiance is represented in left side, and temperature and wind velocity are in right side.
(a)
(b)
(c)
(d)With the data of Figure5 regarding sun irradiance and ambient temperature it is possible to calculate temperature of the solar panel with the equation(17)Tp=Ta+G800NOCT-20°C,where Ta is the ambient temperature, G is the irradiance level, and NOCT is the open-circuit module operation temperature at 800 W/m2 irradiance, 20°C ambient temperature, and 1 m s−1 wind speed. This last parameter is normally included in the manufacturers’ datasheet. Although some authors take wind speed variations into account [49], in the present work they are not considered, as a constant thermal loss coefficient is assumed for wind speeds higher than 1 m s−1. Making use of (17), combined with the data from Figure 5, temperature of solar panels located at the Goddard Space Flight Center (GSFC) on March 13 and March 14, 1971, can be estimated; see Figure 6.Figure 6
Solar panel temperatures on May 13, 1971 (a), and May 14, 1971 (b), at the Goddard Space Flight Center (GSFC). These graphs were calculated according to ambient conditions for those days; see Figure5 and (17).
(a)
(b)Following the procedure described in the previous section and taking into account the aforementioned sun irradiance (Figure5) and panel temperature (Figure 6) data, it is possible to calculate theI-V curve of the selected solar panel (YL280C-30b) in each instant of the specified cloudy and sunny days. In Figures 7 and 8, the evolution of theI-V curve along those days is shown (for one value of time, t, the shape of theI-V curve (see Figure 2) can be appreciated). On the surface plotted in the figures, a thick line has been drawn to indicate the maximum power point (MPP) which, as can be observed, indicates variations regarding the current and voltage at the aforementioned MPP during the chosen days. These current and voltage levels at MPP have been separately plotted in Figures 7(b), 7(c), 8(b), and 8(c). These values of the maximum power available from the solar panel are used in the following section as a reference to compare the studied MPPT methods.Figure 7
Evolution of theI-V curve of YL280C-30b monocrystalline solar panel according to ambient conditions, that is, sun irradiance and panel temperature levels, on May 13, 1971 (a), at the Goddard Space Flight Center (GSFC). Evolution of current (b) and voltage (c) at maximum power point (MPP).
(a)
(b)
(c)Figure 8
Evolution of theI-V curve of YL280C-30b monocrystalline solar panel according to ambient conditions, that is, sun irradiance and panel temperature levels, on May 14, 1971 (a), at the Goddard Space Flight Center (GSFC). Evolution of current (b) and voltage (c) at maximum power point (MPP).
(a)
(b)
(c)


## 4. MPPT Methods Studied

As mentioned in Section1, the aim of the present work is to describe a simple methodology to model a solar panel operating within a wide range of ambient conditions and to analyze different MPPT options for photovoltaic systems (see in Figure 9 the block-diagram corresponding to a photovoltaic system equipped with a MPPT). In the following subsections four MPPT methods are studied as an example of the proposed method:(i)
constant voltage method,(ii)
open voltage method,(iii)
short-current pulse method,(iv)
perturb and observe method.In this comparison example the ideal Boost-converter has been selected due to its simplicity (see Figure 10). However, it should be underlined that other more sophisticated DC-DC converter possibilities could also be included in the simulation with the proposed methodology, this possibility being out of the scope of the present work, which is focused on the solar panel accurate modeling for MPPT analysis.Figure 9
General scheme of a photovoltaic system formed by the solar panels, DC-DC converter controlled by an MPPT method, and load.Figure 10
Sketch of the ideal Boost-converter considered.The principle of the Boost-converter consists of two operational modes, depending on the position of the switch,Tsc. When the switch is closed (on-state) the inductor L stores the energy of the source while the capacitor feeds the load. Otherwise, when the switch is opened (off-state) the only path available for the current is through the diode D and the source feeds the load and charges the capacitor. If the converter works with a commutation period, T0, the duty cycle, α, is the fraction of that period in which Tsc is connected; therefore α ranges between 0 (Tsc is never on) and 1 (Tsc is always on). For an ideal Boost-converter, the relationship between the input and output variables (current and voltage) is defined by(18)Vout=Vin1-α,(19)Iout=1-αIin.Therefore, bearing in mind (17) and (18), the effect of the Boost-converter can be considered as an equivalent resistor, whose value depends on the duty cycle:(20)Req=VinIin=1-α2Rload.Taking the above expression into account, Req increases if α decreases which, considering theI-V curve of a solar panel, entails an increase of the voltage output from the aforementioned panel (see Figure 11). Therefore, the task of an MPPT method would be to find the particular value of the duty cycle which maximizes the output power from the solar panel. As previously mentioned, in the following subsections four MPPT methods are analyzed using the proposed solar panel method, the corresponding algorithms being adapted from the work by Dolara [31].Figure 11
Effect of Boost-converter duty cycle (i.e., the equivalent resistor) variations on the operational point of a solar panel.


### 4.1. Constant Voltage Method

This is the simplest MPPT method among the four considered methods. It is based on maintaining the operating voltage of the panel,V, as close as possible to a reference value, Vref. This reference value can be based on the data supplied by the manufacturer. Although this method represents an improvement in terms of energy efficiency when compared to the case of no MPPT method, it is also true that it can be only optimized for one ambient condition. Therefore, the operational point will not be coincident with the maximum power point of theI-V curve for other sun irradiance and temperature conditions.As can be observed in the corresponding block-diagram of Figure12, if solar panel output voltage, V, is higher than the reference value, Vref, then the duty cycle, α, should be reduced, as indicated in the Boost-converter analysis. Otherwise, if V<Vref it will have to be increased to obtain the opposite effect.Figure 12
Block-diagrams of the different MPPT used in the present work.


### 4.2. Open Voltage Method

This method is based on the principle that the MPP voltage is always a constant fraction of the open-circuit voltage (i.e.,Vmp/Voc=koc), no matter under which sun irradiance and temperature levels the solar panel is operating. This consideration closely matches reality, although the value of the relationship depends on the solar panel being studied and, as previously suggested, is not entirely constant in the face of sun irradiance and temperature variations. A major drawback of this method is the need to know the open-circuit voltage for given conditions. This voltage level is normally measured by opening the circuit. Logically, this operation stops the flow of energy. Therefore, the number of measurements per unit of time should be adjusted to reduce the loss of energy, bearing in mind that too many measurements will capture variations of the ambient conditions but reduce the solar panel output power. The control principle of the Boost-converter is the same as the previously explained method. The duty cycle, α, should be increased so as to lower the panel output voltage and raised to increase it (see the corresponding flow chart in Figure 12).


### 4.3. Short-Current Pulse Method

This method is based on a similar principle to the previous one. In this case the current at MPP is considered to be a fraction of the short-circuit current (i.e.,Imp/Isc=ksc). This fraction is supposed to remain constant in every ambient condition. As with the previous case, this closely matches reality, although there are slight variations depending on sun irradiance and temperature changes. In this case part of the power is lost during the short-circuit current measurements. In this case an increase of the duty cycle, α, raises the solar panel output current, the decrease of the duty cycle produces the opposite effect (see the corresponding flow chart in Figure 12).


### 4.4. Perturb and Observe Method

This algorithm is based on continuous modifications of the solar panel operational voltage, checking after each perturbation if the generated power has increased or decreased. If the power increases the next voltage perturbation will go in the same direction and is changed if the power decreases. This process is indicated in the corresponding block-diagram in Figure12.


## 4.1. Constant Voltage Method

This is the simplest MPPT method among the four considered methods. It is based on maintaining the operating voltage of the panel,V, as close as possible to a reference value, Vref. This reference value can be based on the data supplied by the manufacturer. Although this method represents an improvement in terms of energy efficiency when compared to the case of no MPPT method, it is also true that it can be only optimized for one ambient condition. Therefore, the operational point will not be coincident with the maximum power point of theI-V curve for other sun irradiance and temperature conditions.As can be observed in the corresponding block-diagram of Figure12, if solar panel output voltage, V, is higher than the reference value, Vref, then the duty cycle, α, should be reduced, as indicated in the Boost-converter analysis. Otherwise, if V<Vref it will have to be increased to obtain the opposite effect.Figure 12
Block-diagrams of the different MPPT used in the present work.


## 4.2. Open Voltage Method

This method is based on the principle that the MPP voltage is always a constant fraction of the open-circuit voltage (i.e.,Vmp/Voc=koc), no matter under which sun irradiance and temperature levels the solar panel is operating. This consideration closely matches reality, although the value of the relationship depends on the solar panel being studied and, as previously suggested, is not entirely constant in the face of sun irradiance and temperature variations. A major drawback of this method is the need to know the open-circuit voltage for given conditions. This voltage level is normally measured by opening the circuit. Logically, this operation stops the flow of energy. Therefore, the number of measurements per unit of time should be adjusted to reduce the loss of energy, bearing in mind that too many measurements will capture variations of the ambient conditions but reduce the solar panel output power. The control principle of the Boost-converter is the same as the previously explained method. The duty cycle, α, should be increased so as to lower the panel output voltage and raised to increase it (see the corresponding flow chart in Figure 12).


## 4.3. Short-Current Pulse Method

This method is based on a similar principle to the previous one. In this case the current at MPP is considered to be a fraction of the short-circuit current (i.e.,Imp/Isc=ksc). This fraction is supposed to remain constant in every ambient condition. As with the previous case, this closely matches reality, although there are slight variations depending on sun irradiance and temperature changes. In this case part of the power is lost during the short-circuit current measurements. In this case an increase of the duty cycle, α, raises the solar panel output current, the decrease of the duty cycle produces the opposite effect (see the corresponding flow chart in Figure 12).


## 4.4. Perturb and Observe Method

This algorithm is based on continuous modifications of the solar panel operational voltage, checking after each perturbation if the generated power has increased or decreased. If the power increases the next voltage perturbation will go in the same direction and is changed if the power decreases. This process is indicated in the corresponding block-diagram in Figure12.


## 5. Results

Simulations were carried out with MATLAB. The solar panels were modeled following indications from Section2. The ambient conditions considered are those described in Section 3, the temperature of the solar panels being calculated with (16). In all cases a 10 ms time step has been considered for data acquisition and control, for every MPPT method analyzed. Regarding the short-current pulse method and the open voltage method a measurement of Isc and Voc was taken every 3 seconds (300 time steps). Duty cycle was varied in steps of Δα=0.005 for all methods, the initial value in every calculation being α0=0.12. The control parameters corresponding to the studied methods were optimized, for the solar panel considered in the present work, following these criteria.(i)
Constant voltage method: for the reference voltage value,Vref, the voltage level at MPP, Vmp, indicated by the manufacturer at 35°C, was selected:(21)Vref=31.07-0.41·10=27.2V.(ii)
Open voltage method: as control constant,koc, the ratio of voltage level at MPP, Vmp, to the open-circuit voltage, Voc, at STC, was chosen:(22)koc=31.339.1=0.8.(iii)
Short-current pulse method: as control constant,ksc, the ratio of current level at MPP, Imp, to the short-circuit current, Isc, at STC, was chosen:(23)ksc=8.969.5=0.94.Results are included in Table 2 and Figure 13, together with the maximum extractable power from the solar panel obtained with expressions (13) and (14). The high efficiency obtained for all MPPT methods can be explained as in every case an ideal DC-DC converter (without energy losses) has been considered. The best performance of MPPT method (for the studied conditions) seems to be the perturb and observe method, the reason being the noninterrupted power extraction from the solar panel. The worst performance of MPPT method is the constant voltage method, as it shows the larger influence from the ambient conditions. The open voltage and short-current pulse methods show a quite good performance, the grey zone below the plots being produced by the switching that disconnect the solar panel to perform the measurements of Voc and Isc. The better efficiency showed by the short-current pulse method when compared to the open voltage method could be explained as the ratio Imp/Isc is less dependent on the temperature than the ratio Vmp/Voc.Table 2
Energy obtained from YL280C-30b monocrystalline solar panel calculated with ambient conditions measured at the Goddard Space Flight Center (GSFC) on May 13, 1971 (cloudy day), and May 14, 1971 (sunny day).
MPPT method
Sunny day
Cloudy day
Energy[W h]
η
MPPT
Rank
Energy[W h]
η
MPPT
Rank
Ideal
1880
1
—
526
1
—
Constant voltage
1826
0.971
4
496
0.943
4
Open voltage
1860
0.989
3
513
0.976
3
Short-current pulse
1867
0.993
2
518
0.985
2
Perturb and observe
1871
0.995
1
522
0.992
1Figure 13
Power produced by YL280C-30b monocrystalline solar panel calculated with the studied MPPT methods, for ambient conditions measured at the Goddard Space Flight Center (GSFC) on May 13, 1971 (cloudy day), and May 14, 1971 (sunny day). Maximum extractable power (ideal) has been also included.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)Logically, the efficiency of the studied MPPT methods is worse in case of cloudy days than in case of sunny days. This is clear due to the faster variations on the ambient conditions.


## 6. Case Study

In order to check the proposed methodology, a real photovoltaic facility, whose behavior was measured, together with the ambient conditions, by Houssamo et al. [50], has been studied. This photovoltaic facility is formed by eight Q125PI solar panels manufactured by Conergy, organized in four two-panel series connected in parallel. The characteristics of the panels are included in Table 3, whereas the sun radiation and the panel temperature during the measurements are included in Figure 14.Table 3
Characteristic points of Q125PI solar panel (Conergy), included in the manufacturer’s datasheets at Standard Test Conditions (STC): 1000 W/m2 irradiance, 25°C cell temperature, and AM1.5 g spectrum according to EN 60904-3 [17].
Conergy Q125Pl (polycrystalline)
n
36
T
r (°C)
25
P
mp (W)
125
γ (%/°C)
−0.426
I
mp (A)
7.36
α
I
mp (%/°C)
—
V
mp (V)
17
β
V
mp (%/°C)
−0.352
I
sc (A)
7.94
α
I
sc (%/°C)
0.035
V
oc (V)
21
β
V
oc (V/°C)
−0.074Figure 14
Ambient conditions (irradiance and solar panels temperature) during operation of the photovoltaic facility studied from [50].The procedure described in Section4 was followed in order to model the studied photovoltaic facility. The extracted power was calculated using the perturb and observe (P&O) MPPT algorithm. This power was compared to the measured one from the facility [50], which was programmed following two different MPPT algorithms, perturb and observe (P&O) and incremental conductance (INC). The results from the simulation show a higher extracted power when compared to the behavior measured directly on the facility, this discrepancy between results being explained as no power losses (wiring, connections, dirt over the panels, degradation, losses at the Boost-converter, etc.) were taken into account in the simulation carried out. Bearing in mind that no information regarding these losses was included in [50], the results were scaled down by multiplying by a constant (0.64) and, as a result, a much better correlation was obtained (see Figure 15). Also, an equivalent resistor can be then considered in order to take into account the power losses. The results corresponding to the simulation carried out with a 1.375 Ω resistor connected in series with each pair of solar panels are included in Figure 15. A quite good correlation between the results obtained with the present methodology and the ones measured by Houssamo et al. [50] can be observed in the figure. To obtain a better approximation to the measured results, a combination of the scaled results (losses proportional to the extracted power like the ones produced by dirt in the panels) and the losses resulting from electrical current (taken into account with equivalent resistors) should be considered. This result has been included in Figure 15 combining both kinds of power losses at 50% each one. An excellent correlation with the measured extracted power can be observed, the higher deviation being located where the effect of the MPPT algorithm selection is relevant. The highest differences in extracted power between both MPPT algorithms selected by the authors of [50] were measured in period from 20 s to 25 s, which is precisely the period where the higher deviation of the results obtained by present methodology from the measurement results is observed.Figure 15
Power from an 8-panel photovoltaic facility calculated with the proposed methodology and P&O MPPT algorithm (scaled results and results with power losses considered). The results from the testing measurements [50] (considering P&O and INC MPPT algorithms) are also included in the figure (P&O: blue line; INC: red line).


## 7. Conclusions

In the present work a solar panel modelling to analyze MPPT methods is described. The most significant conclusions derived from the work are the following.(i)
It represents a simple and quick-calculation methodology to obtain a solar panel model to perform MPPT simulations, from the manufacturers’ datasheet.(ii)
The developed method takes into account for the calculations possible variations on the ambient conditions (sun irradiation and solar cells temperature).(iii)
The proposed analytical methodology allows fast MPPT methods comparison.(iv)
Based on the developed methodology the performance of a specific MPPT method can be foreseen when applied to a particular solar panel from which only manufacturer’s datasheet information is available.(v)
Simple, accurate, and low calculation resources demanding way to obtain the maximum extractable power from solar panels are presented.(vi)
The feasibility of the described methodology has been checked with four different MPPT methods applied to a commercial solar panel, within a day, and under realistic ambient conditions. Also a case study was analyzed, the results being compared to testing measured with good correlation.(vii)
Finally, it should also be mentioned that although the parameter-extraction analytical model which is the core of the proposed methodology has proven to be as accurate as several numerical methods [35], the most recently developed algorithms can be an improvement in terms of accuracy (at least if the photovoltaic device is partially shaded).



---

*Source: 914212-2015-03-22.xml*