A Simple and Efficient MPPT Method for Low-Power PV Cells

Small-size PV cells have been used to power sensor nodes.These devices present limited computing resources and so low complexity methods have been used in order to extract themaximumpower from the PV cells. Among them, the fractional open circuit voltage (FOCV)method has been widely proposed, where the maximum power point of the PV cell is estimated from a fraction of its open circuit voltage. Here, we show a generalization of the FOCV method that keeps its inherent simplicity and improves the tracking efficiency. First, a single-diode model for PV cells was used to compute the tracking efficiency versus irradiance. Computations were carried out for different values of the parameters involved in the PV cell model. The proposed approach clearly outperformed the FOCVmethod, specially at low irradiance, which is significant for powering sensor nodes. Experimental tests performed with a 500mW PV panel agreed with these results.


Introduction
The advances in electronics and communication protocols have led to a widespread use of wireless sensor networks (WSN).In most applications, the sensor nodes of the WSN are required to be wireless both for communication and powering.As for the power supply, the use of small-size PV cells or modules has been proposed.The power-voltage (P-V) curve of a photovoltaic (PV) cell or panel presents a maximum power point (MPP) that changes with temperature and irradiance.To extract the maximum power under varying conditions an MPP tracker can be used.Several MPP tracking (MPPT) methods have been proposed in the literature [1][2][3][4].Since sensor nodes present limited computing resources, low complexity MPPT methods are preferred for this particular application.Because the location of the sensor nodes is mostly determined by the application, a wide range of irradiance can be expected at the sensor placement.Thus, a high efficiency is desirable over a wide range of irradiance, and specially for low irradiance as the power source is scarce.
One of the simplest and most popular MPPT methods is the fractional open circuit voltage (FOCV) technique, which estimates the MPP voltage ( MPP ) from a fraction of the open circuit voltage ( OC ); that is, where  MPP,est is the estimated value of the actual  MPP and  is an empirical constant whose value should be set following a thorough characterization of the PV panel under varying meteorological conditions (irradiance and cell temperature). OC is either measured periodically (by momentarily opening the output of the PV panel) or by using a pilot cell (i.e., an additional solar cell of the same type configured in open circuit voltage configuration).Typical reported values for  range from 0.73 to 0.8 depending on the PV panel type and characteristics [2,3].Because of its simplicity, the FOCV method has been recently applied to small-size PV cells in order to power autonomous sensors [5][6][7][8][9].
In this work, we propose to generalize (1) in order to estimate  MPP by using a linear fit; that is, where  and  are empirical coefficients.The use of (2) will be referred to as the linear open circuit voltage (LOCV) method.In fact, the FOCV method can be considered a particular case of the proposed LOCV method with  = 0.Both computed and experimental results of the proposed approach will be presented and compared with the FOCV method.As will be shown, the LOCV method significantly improves the performance of the FOCV method, especially at low irradiance, while maintaining its inherent simplicity.The work presented here builds upon [10], where we first presented (2) and some initial results.

Solar Cell Model
Different equivalent circuits have been used in the literature in order to model the current/voltage (I-V) characteristic of a silicon PV cell [11][12][13][14][15].Among them, the single-diode model, shown in Figure 1, offers a good compromise between simplicity and accuracy [13], whereby  PH is the photogenerated current,  is the cell current,  is the cell voltage, and  s and  p are, respectively, the series and shunt resistances.This model will be used here in order to generate computed data of the I-V curve of a PV cell.
The corresponding expression of the I-V characteristic is given by [16] where  PH has been approximated by  SC , the short circuit current of the cell;  O is the saturation current of the diode;  is the electron charge;  d is the ideality factor of the diode, which, for silicon, is usually between 1 and 2 [5,7];  is the Boltzmann constant; and  cell is the cell temperature in Kelvin.
By considering open circuit conditions ( = 0 and  =  OC ) in (3), we can write the parameter  O as The parameters  SC and  OC in (3) and (4) change with the irradiance and temperature as where  cell is the cell temperature;  is the incident irradiance (in W/m 2 );  and  are the current and voltage temperature coefficients, respectively;  SC r and  OC r are given at a reference irradiance ( r ) and reference cell temperature ( r ); and  OC and  OC are two empirical constants used to model the significant variation of  OC at low .Typically,  r = 1000 W/m 2 (=100 mW/cm 2 ) and  r = 25 ∘ C. Values of  OC = −0.04 and  OC = 1000 W/m 2 are adequate for many silicon PV cells [17].When directly illuminated, solar cells heat up above the ambient temperature ( a ), which is known as the self-heating effect, and  cell can be obtained from where  cell,NOCT , known as the nominal operating cell temperature (NOCT), is the temperature of the cell when exposed to 800 W/m 2 at  a = 20 ∘ C and wind speed of 1 m/s.It is empirically determined, and for silicon solar cells range between 42 ∘ C and 48 ∘ C.

Computed Results
We employed (3) to (7) with the following parameter values typical of monocrystalline solar cells [17]:  SC r = 35 mA/cm 2 ;  OC r = 0.6 V; / = 12.5 A/cm 2 / ∘ C;  = −2 mV/ ∘ C; and  cell,NOCT = 42 ∘ C. The parameters  SC r and  in (5) can be, respectively, obtained by multiplying  SC r and / by the area of the cell.For the computations, a single solar cell with an area of 1 cm 2 was used.Nevertheless, as will be justified at the end of this section, the results and the derived conclusions are equally valid to cells of any area and PV panels composed of an arbitrary number of solar cells disposed in parallel and series.Figure 2 shows the computed I-V and P-V curves at three different levels of irradiance for the particular case of  s = 0,  p = ∞,  d = 1.5, and  a = 25 ∘ C. To obtain the data we simulated the PV cell model in SPICE.The power values were obtained by multiplying  by  at each data point.As can be seen, both  OC and  MPP slightly decrease at the highest irradiance, which is due to the self-heating effect of the PV cell.
From the data of the I-V and P-V curves, several parameters can be obtained such as  OC ,  MPP , and  MPP (power at the MPP).Table 1 shows numerical values of those parameters for an irradiance range from .20 W/m 2 to 1000 W/m 2 .Fourteen points of irradiance, logarithmically equally spaced, were selected to provide a dynamic range around 100 in  MPP .Other cases that will be discussed throughout this section are also shown in Table 1. Figure 3 Table 1: Computed V OC , V MPP ; and P MPP data at fourteen points of irradiance, logarithmically equally spaced, and for different values of the parameters of the PV cell model.

G (W/m 2 )
d = 1.5 (Figure 1) 6)  p = 10 (Figure 8) Figure 2: Generic I-V and P-V plots at several values of  and at  a = 25 ∘ C for a single solar cell with an area of 1 cm 2 .A curve joining the MPPs is also plotted.
represents the fourteen computed points (diamonds) of  MPP versus  OC and two least-squares regression lines fitted to the computed data corresponding, respectively, to the FOCV method, that is, (1) with  = 0.809, and the LOCV method, that is, (2) with  = 0.894 and  = −0.041.As can be seen, the regression line corresponding to the LOCV method better fits the computed data.The inferred parameters of the regression lines (, , and ) were used to obtain  MPP,est at the fourteen irradiance points for each of the two methods, by using (1) and (2), respectively.
The corresponding power values at  MPP,est ,  MPP,est were inferred from the computed P-V curves in order to obtain the tracking efficiency, which is given by  This parameter is used in the literature to compare the performance among different algorithms.Obviously, a value of 1 (100%) is the ultimate goal.Figure 4 shows the computed values of  MPPT versus  at  a = 25 ∘ C for the fourteen irradiance points.We added the results at two more temperatures, 0 ∘ C and 50 ∘ C. For these temperatures, the P-V curves were recalculated but we still used the same regression lines of Figure 3.This makes sense, as a PV panel can be characterized at a single temperature, for example, 25 ∘ C, and the calculated regression lines used for International Journal of Photoenergy  the full working temperature range.As can be seen, at low irradiance the LOCV method clearly outperforms the FOCV method.At higher irradiance, a rather high value of  MPPT (>99%) is achieved by both methods, although the FOCV method presents the lowest efficiency from . 100 W/m 2 to 1000 W/m 2 at  a = 0 ∘ C. The value of  MPPT for the LOCV method was always higher than 99.8% at the three computed temperatures.More computations were carried out at  a = 25 ∘ C for  d = 1 and  d = 2 (see Table 1).Again, better linear fits were obtained with the LOCV method (not shown).Figure 5 shows the corresponding computed values of  MPPT .For each of the cases, the parameters of the corresponding regression lines are provided.Again, the LOCV method clearly outperforms the FOCV method at low irradiance and slightly at medium irradiance.
Finally, computations were performed for  d = 1.5,  a = 25 ∘ C, nonzero values of  s , and finite values of  p .In [13] normalized values for  s and  p were defined as This normalization allows for an immediate comparison among different PV modules (of different area and characteristics).Based on [13], in our work we considered the following values: from 0.01 to 0.1 for  s and from 100 to 10 for  p .The performance for both  s = 0.01 and  p = 100 was almost identical to that shown in Figure 2 (for  a = 25 ∘ C), with the LOCV method outperforming the FOCV method.So these V MPP, est = 0.926V OC − 0.037 V MPP, est = 0.864V OC − 0.043 results are not shown here.Table 1 shows the data for the limiting cases  s = 0.1 (with  p = ∞) and  p = 10 (with  s = 0).For the case of  p = 10, the data for the lowest two irradiance levels were not used as they provided negligible values of  MPP (well below of 1% of the resulting  MPP at 1000 W/m 2 ).As for  s = 0.1, Figure 6 shows the computed values of  MPP versus  OC and the fitted regression lines.Due to the high value of  s (highest limit), the data values present a folded form at the highest irradiance levels.So the regression lines of the FOCV and LOCV methods cannot fit the data corresponding to the high irradiance levels as well as that in Figure 3. Otherwise, both lines are very similar in this case.Consequently, the computed values of  MPPT , shown in Figure 7 ( s = 0.1), are quite similar (and indeed relatively high) for both methods.As for  p = 10, Figure 8 again shows the computed values of  MPP versus  OC and the fitted regression lines.Due to the low relative value of  p (lowest limit) the regression line of the LOCV method cannot fit the data corresponding to the low irradiance levels as well as that in Figure 3.Even so, the computed values of  MPPT , also shown in Figure 7, still present a high efficiency, outperforming the FOCV method at all the irradiance levels, but specially at the low ones.Finally, we computed  MPPT for  s = 0.1 and  p = 10 (not shown).In that case, the LOCV method also outperformed the FOCV method at low and medium irradiance levels.
Increasing the PV cell area or adding identical PV cells in parallel will scale up the values of currents and thus of powers but the values of  OC and  MPP will remain the same and so the derived tracking efficiencies.Tracking efficency will also   remain constant by adding PV cells in series: both  MPP and  OC will scale up by the number of serial cells but their ratio will remain constant and so the derived tracking efficiencies.

Experimental Results
The LOCV method was tested with a 500 mW ( SC = 160 mA,  OC = 4.6 V) PV panel (MSX-005, Solarex) and compared with the FOCV method.These low-power panels are used, for example, to power autonomous sensors [5][6][7][8][9][10].In order to achieve reproducible results, we implemented a PV panel simulator by connecting a current source (GS610, Yokogawa) in parallel with the PV panel, which was coated with an opaque cover (Figure 9).In this way, the short circuit current ( SC ) of the PV panel was adjusted by the current source, emulating different levels of irradiance.Since the panel was not illuminated,  celL,NOCT = 20 ∘ C (i.e., the panel is not overheated).The current source was configured to cover the full range of the PV panel, varying from 5 mA to 158 mA in 9 mA steps.The PV panel simulator was characterized by using the GS610's measurement unit to measure the panel's voltage, a 2001 multimeter (Keithley) to measure the current of the panel, and a programmable voltage source (Agilent E3631A) in parallel with a 10 Ω/1 W resistor acting as a load.Figure 9 shows the experimental setup.All the instruments were controlled via the GPIB bus with a dedicated program using the graphical development environment LabVIEW.For each current value ( SC ), the voltage of the E3631A was increased from 0 V to 5 V in 0.1 V steps.PV output voltages and currents were measured and the corresponding power values were calculated in order to obtain the I-V and P-V curves.From each P-V curve the values of  OC ,  MPP , and  MPP were obtained.The limit values for  MPP were, respectively, 8.2 mW ( SC = 5 mA) and 545.9 mW ( SC = 158 mA).
Figure 10 represents the experimental data of  MPP versus  OC and two fitted least-squares regression lines corresponding to the FOCV and LOCV methods.As can be seen, the regression line corresponding to the LOCV method better fits the experimental data.From the two regression lines, the values of  MPP,est corresponding to the FOCV and LOCV methods were derived.Then, from the P-V curves, the values of  MPP,est and  MPPT were obtained.Figure 11 V OC (V) FOCV method at low irradiance, with  MPPT always being higher than 99.96%.The efficiency increase at low irradiance is of significant importance in order to power sensor nodes.

Conclusion
PV cells have been proposed in the literature in order to power the sensor nodes of WSN.Because of the limited computing capabilities of the sensor nodes, simple MPPT methods have to be used.Among them, the FOCV method has been widely proposed and used.Tracking efficiencies, though, are lower than that achieved with more complex methods.In this work, we have proposed the LOCV method, which outperforms the FOCV method while maintaining its inherent simplicity.Computations show that the LOCV method achieves a high efficiency for all the irradiance range whereas the FOCV method fails in achieving a high efficiency at low irradiance levels for most of the cases.The presence of extremely low values of shunt resistance of the PV cell negatively impacts the achieved efficiency on both methods but specially that of the FOCV method.Finally, experimental data from a low-power 500 mW PV panel confirmed the good performance of the LOCV method for a wide range of irradiance, which is of significant value for powering sensor nodes.

Figure 1 :
Figure 1: Single-diode model of a solar cell with series ( s ) and shunt ( p ) resistances.

Figure 3 :
Figure 3: Computed  MPP versus  OC for a single PV cell.Two least-square regression lines are also represented: the grey line corresponds to the FOCV method with  = 0.809 whereas the black line corresponds to the LOCV method with  = 0.894 and  = −0.041.The parameters of the PV cell model are  a = 25 ∘ C,  d = 1.5,  s = 0, and  p = ∞.

Figure 4 :
Figure 4: Computed  MPPT versus  by using the FOCV and LOCV methods for a single PV cell at different values of  a and with  d = 1.5,  s = 0, and  p = ∞.

Figure 5 :
Figure 5: Computed  MPPT versus  by using the FOCV and LOCV methods for a single PV cell at different values of  d and with  a = 25 ∘ C,  s = 0, and  p = ∞.

Figure 6 :n d = 1 . 5
Figure 6: Computed  MPP versus  OC for a single PV cell with  a = 25 ∘ C,  d = 1.5,  p = ∞, and  s = 0.1.Two regression lines corresponding to the FOCV and LOCV methods are also represented.

Figure 7 :
Figure 7: Computed  MPPT versus  by using the FOCV and the LOCV methods for a single PV cell with  a = 25 ∘ C, n  = 1.5, and for both  s = 0.1 (and  p = ∞) and  p = 10 (and  s = 0).

6 rFigure 8 :
Figure 8: Computed  MPP versus  OC for a single PV cell with  a = 25 ∘ C,  d = 1.5,  s = 0, and  p = 10.Two regression lines corresponding to the FOCV and LOCV methods are also represented.

Figure 9 :
Figure 9: The PV array simulator and the setup used for its characterization.

Figure 10 :Figure 11 :
Figure 10:  MPP versus  OC for a 500 mW panel simulator at  a = 25 ∘ C and two fitted least-square regression lines corresponding to the FOCV and LOCV methods. )