Advanced
Developed MPPT Algorithm for Photovoltaic Systems without a Voltage Sensor
Developed MPPT Algorithm for Photovoltaic Systems without a Voltage Sensor
Journal of Power Electronics. 2013. Nov, 13(6): 1042-1050
Copyright © 2013, The Korean Institute Of Power Electronics
  • Received : April 30, 2013
  • Accepted : September 07, 2013
  • Published : November 30, 2013
Download
PDF
e-PUB
PubReader
PPT
Export by style
Share
Article
Author
Metrics
Cited by
TagCloud
About the Authors
Milad Momayyezan
School of Electrical and Computer Engineering, College of Engineering, University of Tehran, Tehran, Iran
Hossein Iman-Eini
School of Electrical and Computer Engineering, College of Engineering, University of Tehran, Tehran, Iran
imaneini@ut.ac.ir

Abstract
This paper presents a study of maximum power point tracking (MPPT) for photovoltaic arrays with only one current sensor. Initially, a review of MPPT methods with only a current sensor is performed with extension for a variety of dc/dc converters. Furthermore, the same topology is developed to achieve better performance in the presence of sensor offset and environmental noise. The proposed method is robust, cost effective, and behaves well dynamically and in the steady state. After a theoretical analysis of presented approach, its validity and effectiveness are verified by simulation and experimental results.
Keywords
I. INTRODUCTION
In recent years, renewable energy sources, especially photovoltaic (PV) panels and wind turbines have experienced increased attention in terms of power generation. They are pollution-free, reliable and considered to be economical sources of electricity for remote areas. However, the installation costs of PV systems are high [1] , and they suffer from low efficiency. Additionally, the nonlinear characteristics of solar cells and their dependency on environmental conditions have a major impact on the complexity of their control techniques [2] .
To overcome the above mentioned issues, a variety of MPPT algorithms has been proposed to date. The inherent characteristics of each method make them suitable for particular applications. For instance, [3] , [4] present new approaches to solve tracking confusion under rapidly changing climate conditions. Reference [5] makes a comparative study of different MPPT methods in terms of operating conditions, number of sensors, convergence speed, and system cost. Considering above parameters, each method may be considered only for specific applications. According to [6] , the Perturb and Observe (P&O) algorithm has received a great deal of attention due to its simplicity and ease of implementation [7] . A comparative study of low-cost topologies demonstrates the effectiveness of P&O [8] . Accordingly, a novel algorithm is presented in [9] . The proposed method is implemented by a simple microcontroller, utilizes only one current sensor and has a high overall efficiency [10] . On the other hand, it has been developed only for step-down buck converters and suffers from inherent oscillations around the MPP. According to [11] , the boost topology is superior for PV applications with its lower implementation cost and better dynamic response when compared with buck converters. Furthermore, in order to avoid the series connection of PV arrays, it is necessary to boost the input voltage. Therefore, in the first part of this paper, the concepts of the available topologies with only a current sensor are extended for different kinds of dc/dc converters. In the second part, the impact of sensor offset and environmental noise on the same methodologies are investigated. Referring to [12] , [13] , the measurement noise and sensor offset may deteriorate the performance of tracking algorithms. Hence, an analysis of environmental effects and proposing a robust algorithm against noise and offset errors would be a promising solution in practice and for critical applications.
II. MPPT WITH ONLY ONE CURRENT SENSOR
The major parameters that affect I-V characteristic of a PV module are temperature and irradiation. Under all conditions, there is only one operating point for a module with the highest power delivery. Using a power converter between the solar module and the load, it is possible to control cell voltage and track the maximum power point. Fig. 1 shows a simplified block diagram of such configuration with the application of a boost converter. The output voltage is almost constant due to the connection of battery packs in parallel with the load. Hence, the controller determines the converter duty cycle in a way that maximum power is delivered to the load at different temperature and irradiation levels.
PPT Slide
Lager Image
Application of boost converter as MPPT.
According to (1), the duty cycle represents the switching on-time within each cycle. Eqns. (2)-(4) demonstrate the relation between the cell input voltage, Vpv and the load voltage, Vbat , in buck, boost, and buck-boost converters, respectively.
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
Considering (2)-(4), Vpv is written as:
PPT Slide
Lager Image
where, D * is the equivalent of ( 1/D ) in the buck converter, (1-D) / D in the buck-boost converter, and (1-D) in the boost converter while the converter is operating in the continuous conduction mode (or CCM). Neglecting the system losses, the input power to the converter may be written as:
PPT Slide
Lager Image
Since Vbat is almost constant, the following definition is introduced for the equivalent power, P * .
PPT Slide
Lager Image
The principal operation of the P&O algorithm is defined in a way that the cell voltage Vpv is a little perturbed and the impact on the delivered power, Pin is continuously observed. This procedure is repeated until condition (8) is achieved.
PPT Slide
Lager Image
Substitution of (7) into (6) and derivation with respect to Vpv leads to:
PPT Slide
Lager Image
Due to the fact that the battery charging/discharging speed is much lower than the MPPT response, the value of
PPT Slide
Lager Image
is set to zero which yields:
PPT Slide
Lager Image
In different kinds of dc/dc converters, the parameter D * can be determined according to the duty cycle value, D . As a result, in general form, the parameter D is perturbed and the equivalent power, P * is observed. Using this approach, only one current sensor is needed to estimate P * while conventional methods need at least two sensors to estimate Pin .
III. IMPROVEMENT OF THE STEADY STATE AND DYNAMIC RESPONSE
In classic MPPT algorithms, a tradeoff should be made between the speed of the dynamic response and the steady state behavior by appropriate selection of the duty cycle steps. Therefore, in order to achieve a high performance under all conditions, an upadting procedure must be realized. In [14] , a modified adaptive method has been proposed for the Hill Climbing approach where the tracker automatically modifies the value of ΔD to improve the MPPT performance. However, in this paper a simpler algorithm is developed that uses a dynamical variation of ΔD within the tracking process. Fig. 2 illustrates a flowchart of the proposed method.
According to Fig.2 , the duty cycle value D and the step size ΔD are initialized with D0 and ΔD0 , respectively. Then, Ipv is measured by the current sensor and the paramtere P * is calculated according to (7). In the next step, ΔP *= P *( t + Δt )- P *( t ) is computed and compared with the extreme value “ e ” that defines the acceptable oscillation around the MPP. If the conditional case | ΔP *|> e is true, then the duty cycle is updated as follows.
PPT Slide
Lager Image
Algorithm flowchart for presented tracker.
When ΔP * is positive, the duty cycle is updated with D1 = D0 + K × ΔD0 , where K is a constant integer, e.g. K =1 . The negative sign on the other hand, means that the operating point has already passed the MPP. In this case, ΔD is reduced by a factor lower than unity for precise tracking of the MPP. The sign of K is then reversed and the value of the duty cycle is tuned for the next cycle.
To achieve fast tracking, the value of ΔD is chosen to be large enough at initial steps. Whenever ΔP * becomes negative, the sign of ΔD is reversed and its value is reduced.This procedure will be repeated until the value of | ΔP *|(or the oscillation magnitude around the MPP) becomes smaller than the extreme limit “ e ”. When the variation of the step sizeis stopped (or the converter operates at the steady state), it should be able to exit the current situation by the irradiation level or temperature changes. To achieve this goal, the controller updates ΔD with a larger value at the steady state and waits for environmental changes. Using this strategy, the tracking time, oscillations around the MPP, and the power loss are reduced.
IV. SIMULATION RESULTS
To evaluate the behaviour of the proposed algorithm, the PV system shown in Fig. 1 is modelled in the MATLAB/SIMULINK environment [15] , based on the configuration shown in Fig. 3 . In this case, two parallel modules (“solar module 36/45”) supply the load. The characteristic parameters of these modules are given in Table I .
PPT Slide
Lager Image
Simplified solar cell model.
UTILIZED PARAMETERS FOR PV MODULES
PPT Slide
Lager Image
UTILIZED PARAMETERS FOR PV MODULES
The simulation results for the controlled system with the proposed MPPT algorithm are shown in Fig. 4 . The represented waveforms are the input power Pin , the equivalent power P *, the PV current Ipv , the PV voltage Vpv , and the corresponding duty cycle as a function of time. The estimated maximum delivered power to the converter was approximately 41W, while a value of 90W, according to the module datasheet, was expected. This is due to the fact that the selected parameters, i.e. Isc =3.8A and Voc = 18V are chosen based on the real test conditions of Fig. 9 (a). In this experiment, the short circuit current and open circuit voltage were far from the ideal condition because of a low irradiation level and cell temperature.
As shown in Fig. 4 , the converter duty cycle is increased by the controller from the initial value ( D =0.5) to the point in which the maximum power is delivered to the load, i.e. D =0.59. It is also obvious that the PV output power P and the defined equivalent power P * demonstrate similar behaviour.
PPT Slide
Lager Image
Simulation results for boost converter with improved MPPT algorithm, (a) Output power of PV panels; (b) Defined equivalent power P*; (c) Output current of PV panels; (d) Output voltage of PV panels; and (e) Converter duty cycle versus time
V. INVESTIGATION OF THE SENSOR OFFSET AND NOISE EFFECT ON MPPT WITH ONLY ONE CURRENT SENSOR
In this section, the effect of erroneous current measurements on the deviation from the ideal MPP is investigated. As stated in [1], MPPT algorithms typicallyutilize the cell short circuit current or open circuit voltage to tune the duty cycle and track the MPP. Any erroneous measurements in such approaches result in incorrect MPPT, while the stated procedure with one current sensor overcomes this issue.
First, the error is considered to be a constant deviation (or offset) from the ideal value, which is common in current sensors. Electronic amplifiers and analog-to-digital converters may be sources of dc offset. Fig. 5 shows the evaluation of the tracker in the presence of dc offset.
PPT Slide
Lager Image
Evaluation of MPPT algorithm in presence of sensor offset
The module output power has been demonstrated for two different cases: considering the sensor offset and without any offset (or the ideal condition). According to Fig. 5 , the sensor offset causes incorrect estimation of the output power. The magnitude of the error depends on the cell operating point and increases by getting closer to the MPP. However, the mentioned issue will not affect the tracker operation since the corresponding duty cycle for MPP (in presence of the dc offset) is equal to the case where the offset equals zero. This advantage simplifies the implementation regarding the dc offset in measurements.
In the second part, the measurement error is considered to be a random type or (white noise). The tracker suffers from this type of error due to the surrounding environment, dc power supplies or the switching action of power converters. Like many conventional methods such as P&O and Hill Climbing, the presented method uses a differential term, i.e. ΔP *, to track the MPP. However, the application of differential terms will result in a low signal to noise ratio (or SNR) and a deterioration of tracker.
The following simulation evaluates the performance of P&O and the suggested method in the presence of white noise. In this simulation, White Gaussian Noise is added deliberately to the sensed current. The corresponding result for a noise with a power of -35db is demonstrated in Fig. 6 .
PPT Slide
Lager Image
Evaluation of noise effect on MPPT algorithm, (a) P&O algorithm; (b) suggested algorithm with one current sensor.
As illustrated in Fig. 6 , both algorithms, which are implemented based on differential terms, are not accurate enough and may fail. Dark circles 1 and 2 highlight the deviation of the tracker during transient and steady state conditions, respectively. Such deteriorations lead to higher losses and weak performance in real conditions. As an alternative, a partial derivative expression is used instead of differential terms. Although this modification adds some complexity, it may be a solution to reduce the noise effect and to improve the steady state performance. In addition, because of accurate MPP detection, there is no oscillation around the MPP which implies higher efficiency.
VI. PROPOSING A ROBUST ALGORITHM FOR NOISY ENVIRONMENTS
As stated in section 2, the following equation may be employed for tracking the MPP in different kinds of dc/dc converters:
PPT Slide
Lager Image
Inserting the value of P * according to (7) into (11), yields:
PPT Slide
Lager Image
Substituting the value of D * , e.g. for the boost converter in CCM, leads to the following condition:
PPT Slide
Lager Image
To solve the differential equation of (13), the partial term
PPT Slide
Lager Image
is rephrased. In this case, first Ipv is calculated by substitution of the following relations (derived from Fig. 3 and equation (4)) into (13):
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
Then Ipv is written as:
PPT Slide
Lager Image
As a result, the following relation for
PPT Slide
Lager Image
is derived:
PPT Slide
Lager Image
Or:
PPT Slide
Lager Image
Substituting (19) into (13) leads to (20) which is a mathematical equation without any differential terms.
PPT Slide
Lager Image
The parameter REF is used to represent the value of
PPT Slide
Lager Image
According to (20), the tracker should modify the duty cycle D and measures the current Ipv until the conditional case “ REF = 0” is achieved. In practice, reaching this condition is not applicable due to the discontinuity of the duty cycle. Therefore, the best operating point is where the magnitude of REF becomes minimum (as close as possible to zero). It should be noted that (20) is valid for the boost converter while it can be recalculated for typical dc/dc converters in a similar way.
VII. INVESTIGATION OF A ROBUST ALGORITHM BY SIMULATION
To evaluate the performance of the proposed method, the last simulation is repeated in the presence of white noise. In this case, the module output power and the value of REF (or
PPT Slide
Lager Image
) are demonstrated as a function of the duty cycle. Fig. 7 shows the simulation results while the system is asymptotically stable and any variations in the duty cycle at the steady state will be returned to the MPP by the tracker.
PPT Slide
Lager Image
Investigation of proposed MPPT algorithm based on eq. (20), (a) Module output power and the value of REF as a function of duty cycle; (b) Tracking points around MPP with and without considering white noise.
According to Fig. 7 (a), the maximum power point corresponds to a duty cycle in which the value of REF is almost zero. Fig. 7 (b) demonstrates the robustness of the proposed algorithm in terms of noisy current where tracking points are shown for two cases: considering white noise (blue points) and without noise (red points). Comparing the obtained results with those in Fig. 6 (a)-(b), the robustness of the tracker is confirmed. In the presented method with the concept of eq. (20), the parameters dependent on module temperature may affect the efficiency of the PV system, [16 - 17] . Hence, for successful MPPT under all conditions, the effects of temperature variations on VT and I0 need to be considered. Equations (21) and (22) describe the dependency of I0 and VT on temperature, respectively.
PPT Slide
Lager Image
PPT Slide
Lager Image
where, I0N and VTN are the saturation current and thermal voltage at standard temperature TN =300°k. Eg is the silicon band gap energy, k is the Boltzmann constant, and q is the electric charge of the electron. Referring to the literature, valuable contributions towards the model based methods for estimating the temperature of PV modules are available in [17 - 19] . In this paper, the NOCT model described in [17] is employed to estimate the cell temperature. The MPPT algorithm with the improved characteristics of being insensitive to noise, dc offset, and converter topology is shown in Fig. 8 . In the beginning, the cell temperature and the PV current Ipv are measured. Then, the mathematical expression of (20) is calculated as a reference (termed as REF). This value is compared to the extreme parameter “e”, which defines an acceptable deviation from the MPP. Then, the value of the duty cycle is updated within each period until the conditional case is satisfied. Like the explained algorithm in Fig. 2 , the duty cycle steps are reduced whenever the tracker passes the MPP. Moreover, to exit from the steady state in the case of environmental changes, the tracker looks for variations and modifies the parameter ΔD with a larger step size. The estimation of the cell temperature will be performed with a lower frequency due to the slower dynamics of temperature variations.
PPT Slide
Lager Image
Flowchart of new MPPT algorithm.
In order to analyse the performance of the explained algorithm, a setup similar to the structure in Fig. 1 has been implemented. The hardware setup includes two parallel modules (“solar module 36/45”), a low cost microcontroller, a boost converter, and three paralleled 60 Ah-12V lead-acid batteries with a resistive load. Fig. 11 shows an image of the hardware setup.
VIII. EXPERIMENTAL RESULTS
The first experiment was carried out at around 1:30 PM. Fig. 9 (b-d) show the tracking process of a controller based on the presented MPPT algorithm.
PPT Slide
Lager Image
Evaluation of proposed MPPT algorithm at 1:30 PM, (a) I-V characteristic of two parallel modules; (b) P-V characteristic along with tracking points; (c) three final points around MPP; and (d) PV voltage and PV current waveforms during MPP tracking.
As shown, the first three steps are larger to achieve a fast dynamic response. Then, the duty cycle step reduces to minimize the oscillations around the MPP and to improve the steady state performance. Fig. 9 (c) shows the last three points that the controller oscillates between. It should be noticed that the oscillation margin may be reduced more by decreasing the extreme value “e”. The variations of the PV voltage and current during the tracking process are shown in Fig 9 (d). According to this figure, the tracker found the MPP after 4 iterations (or 2 seconds) from system start up. As time passes, the voltage and current variations around the MPP reduce and consequently the power oscillations become lower.
The last experiment investigates the performance of the proposed algorithm under different conditions. This test was performed at around 6:30 PM where the irradiation level was lower. The results are shown in Fig.10 . It should be noted that the tracking is processed in existence of the dc offset from the current sensor.
Because of air pollution and air dust, the I-V characteristics of the solar panels were different from the ideal condition. According to the cell datasheet, the short circuit current for each module should be around 2.98A. However, in practice the maximum short circuit current was 2A. The measured I-V characteristics for the two parallel modules just before testing and for two different hours of a summer day are illustrated in Fig. 9 (a) and 10 (a).
PPT Slide
Lager Image
Evaluation of proposed MPPT algorithm at 6:30 PM, (a) I-V characteristic of two parallel modules; (b) P-V characteristic along with tracking points; and (c) three final points around MPP.
PPT Slide
Lager Image
Experimental setup, (a) boost converter with controller; (b) “36/45 solar modules” employed in the experiment.
IX. CONCLUSIONS
In this paper, after studying conventional algorithms that utilize only one current sensor, a novel MPPT method is proposed. The major characteristics of the presented algorithm are robustness against environmental noise with acceptable performance in the steady state and under transient conditions. Furthermore, it may be implemented for different kinds of dc/dc converters with the application of a simple and low-cost microcontroller. To improve the dynamic response and minimize the oscillations around the MPP, variable step size method is employed. The validity and performance of the new method are verified by simulation and experimental results.
BIO
Milad Momayyezan was born in Isfahan, Iran, in 1988. He received his B.S. in Electrical Engineering from the University of Tehran, Iran, in 2011. In the same year, he started his M.S. education at the Royal Institute of Technology (KTH), Stockholm, Sweden. He is currently with the Institute for Power Electronics and Electrical Drives (ISEA), Aachen, Germany. His current research interests include renewable energy systems, power convertors, and electric vehicles.
Hossein Iman-Eini received his B.S. and M.S. from the University of Tehran, Tehran, Iran, in 2001 and 2003, respectively, and his Ph.D. from both the University of Tehran and the Grenoble Institute of Technology, Grenoble, France, in 2009, all in Electrical Engineering. He is currently an Assistant Professor in the School of Electrical and Computer Engineering, University of Tehran. His current research interests include the modeling and control of power converters, multilevel converters and renewable energy systems.
References
Enslin J. H. R. , Wolf M. S. , Snyman D. B. , Swiegers W. 1997 “Integrated photovoltaic maximum power point tracking converter,” IEEE Trans. Ind. Electron 44 (6) 769 - 773    DOI : 10.1109/41.649937
Patel H. , Agarwal V. 2008 “Maximum power point tracking scheme for PV systems operating under partially shaded conditions,” IEEE Trans. Ind. Electron 55 (4) 1689 - 1698    DOI : 10.1109/TIE.2008.917118
Sera D. , Teodorescu R. , Hantschel J. , K.noll M. 2008 “Optimized maximum power point tracker for fast-changing environmental conditions,” IEEE Trans. Ind. Electron. 55 (7) 2629 - 2637    DOI : 10.1109/TIE.2008.924036
Fermia N. , Granozio D. , Petrone G. , Vitelli M. 2007 “Predictive & adaptive MPPT perturb and observe method,” IEEE Trans. Aerosp. Electron. Syst. 43 (3) 934 - 950    DOI : 10.1109/TAES.2007.4383584
Esram T. , Chapman P. L. 2007 “Comparison of photovoltaic array maximum power point tracking techniques,” IEEE Trans. Energy Convers. 22 (2) 439 - 449    DOI : 10.1109/TEC.2006.874230
Salas V. , Olías E. , Barrado A. , Lázaro A. 2006 “Review of the maximum power point tracking algorithms for stand-alone photovoltaic systems,” Solar Energy Materials & Solar Cells 90 (11) 1555 - 1578    DOI : 10.1016/j.solmat.2005.10.023
Femia N. , Petrone G. , Spagnuolo G. , Vitelli M. 2005 “Optimization of perturb and observe maximum power point tracking method,” IEEE Trans. Power Electron 20 (4) 963 - 973    DOI : 10.1109/TPEL.2005.850975
Hohm D. P. , Ropp M. E. 2003 “Comparative study of maximum power point tracking algorithms,” Journal of Progress in Photovoltaics: Research and Applications, Wiley Interscience 11 (1) 47 - 62    DOI : 10.1002/pip.459
Salas V. , Olías E. , Lázaro A. , Barrado A. 2005 “New algorithm using only one variable measurement applied to a maximum power point tracker,” Solar Energy Materials and Solar Cells 87 (1-4) 675 - 684    DOI : 10.1016/j.solmat.2004.09.019
Salas V. , Olias E. , Lazaro A. , Barrado A. 2005 “Evaluation of a new maximum power point tracker (MPPT) applied to the photovoltaic stand-alone systems,” Solar Energy Materials & Solar Cells 87 (1-4) 807 - 815    DOI : 10.1016/j.solmat.2004.07.053
Xiao W. , Ozog N. , Dunford W. G. 2007 “Topology study of photovoltaic interface for maximum power point tracking,” IEEE Trans. Ind. Electron 54 (3) 1696 - 1704    DOI : 10.1109/TIE.2007.894732
Al-Atrash H. , Batarseh I. , Rustom K. 2010 “Effect of measurement noise and bias on Hill-Climbing MPPT algorithms,” I IEEE Trans. Aerosp. Electron. Syst. 46 (2) 745 - 760    DOI : 10.1109/TAES.2010.5461654
Latham A. , Sullivan C. 2010 “Optimization of a continuous-time maximum power point tracking algorithm in the presence of noise,” in IEEE Workshop on Control and Modeling for Power Electronics
Xiao W. , Dunford W. G. 2004 “A modified adaptive hill climbing MPPT method for photovoltaic power systems,” in Proc. 35th Annu. IEEE Power Electron. Spec. Conf. 1957 - 1963
Walker G. 2001 “Evaluating MPPT converter topologies using a MATLAB PV model,” J. Electr. Electron. Eng. 21 (1) 49 - 56
Nishioka K. , Hatayama T. , Uraoka Y. , Fuyuki T. , Hagihara R. , Watanabe M. 2003 “Field-test analysis of PV-system-output characteristics focusing on module temperature,” Journal of Solar Energy Materials and solar cells 75 (3) 665 - 71    DOI : 10.1016/S0927-0248(02)00148-4
Mattei M. , Notton G. , Cristofari C. , Muselli M. , Poggi P. 2006 “Calculation of the polycrystalline PV module temperature using a simple method of energy balance,” Renewable Energy. 31 (4) 553 - 567    DOI : 10.1016/j.renene.2005.03.010
King D. L. , Boyson W. E. , Kratochvil J. A. 2004 “Photovoltaic Array Performance Model,” Sandia National Lab. SAND2004-3535
Krishna H. A. , Misra N. K. , Suresh M. S. 2009 “Solar cell as a temperature sensor for measuring temperature of solar panels in Satellites,” Journal of Instrum. Soc. of India 39 (1)