International Journal of PhotoenergyVolume 2012 (2012), Article ID 672765, 7 pagesdoi:10.1155/2012/672765

Research Article

MPPT Based on Sinusoidal Extremum-Seeking Control in PV Generation

R. Leyva,¹ C. Olalla,¹ H. Zazo,¹ C. Cabal,² A. Cid-Pastor,¹ I. Queinnec,2,3 and C. Alonso2,3

¹DEEEA, Universitat Rovira i Virgili, Avinguda Paisos Catalans 26, 43007 Tarragona, Spain²LAAS, CNRS, 7 Avenue du Colonel Roche, 31077 Toulouse, France³Université de Toulouse, UPS, INSA, INP, ISAE, LAAS, 31077 Toulouse, France

Received 25 October 2011; Accepted 8 December 2011

Academic Editor: Jean-Louis Scartezzini

Copyright © 2012 R. Leyva et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The paper analyses extremum-seeking control technique for maximum power point tracking circuits in PV systems. Specifically, the paper describes and analyses the sinusoidal extremum-seeking control considering stability issues by means a Lyapunov function. Based on this technique, a new architecture of MPPT for PV generation is proposed. In order to assess the proposed solution, the paper provides some experimental measurements in a 100 W prototype which corroborate the effectiveness of the approach.

1. Introduction

Photovoltaic panels require adapting the output-port voltage to extract the maximum deliverable power when weather conditions change. Such energy adaptation should have a good performance; that is, the operation point should be ensured to be near to the maximum and the adapting electronic system should behave reliably and efficiently. A reliable behavior means that the adaptation mechanism is dynamically stable and does not induce great stress to its components. Efficiency of the adapting electronic system depends on the losses in switching converters which adapt the voltage levels between the PV panels and the loads [1].

There exist several approaches to implement maximum power point tracking (MPPT) in PV systems. Remarkable surveys on this subject can be found in [2, 3]. In their survey paper, Onat [2] classifies MPPT algorithms in two categories: indirect and direct methods. The indirect methods use prior information or a mathematical characterization of the PV panel and take measurements far from the desired operating point in order to estimate the voltage or current at maximum power point (MPP). The main drawback of this approach is that the PV system does not operate near to the MPP during the measuring time interval. Some instances of indirect methods are curve-fitting, look-up table, PV panel short-circuit, and PV panel open-circuit method. Another drawback of indirect methods is that they may cause large stress to components since measuring process may involve large changes in working conditions. The last two disadvantages have prompted some authors to use direct methods.

Direct methods take measurements of PV panel current and voltage and their corresponding time derivatives at the operating point in order to drive PV systems toward the maximum. Given that direct methods do not carry out abrupt changes in the operating point to make measurements, measuring frequency can be higher, and therefore these MPPT methods can faster track the optimum power point.

The most common direct MPPT methods are the P&O algorithm [4] and the Incremental Conductance method [5]. Other techniques as fuzzy control [6, 7] and neural networks are also being used to develop MPPT controllers. This is because of the nonlinear relationships involved in photovoltaic systems. Tuning MPPT circuits are a difficult task, and a bad tuned MPPT circuit becomes instable when large disturbances occur. It is important to point out that instabilities may occur in any MPPT circuit when inherent delays in digital implementation or filters in analog implementations are not taken into account in the MPPT circuit design.

In parallel, with the development of MPPT methods in renewable energy field, some researchers have studied the technique named extremum seeking control (ESC) in the field of automatic control. This technique is concerned with algorithms that seek the maximum or the minimum of a nonlinear map. A system governed by ESC autooscillates around the optimum or the oscillation is forced by a sinusoidal signal. Autooscillating ESC algorithm is reported by Morosanov [8] and has been adapted to PV systems in [9]. Sinusoidal ESC has also been successfully applied to track the MPP of PV systems. This technique uses a sinusoidal perturbation to estimate the gradient of the voltage-power curve. Using this gradient value, ESC drives the PV system to the MPP.

Nevertheless, sinusoidal ESC is still in an incipient stage in PV systems, despite the research of Brunton et al. [10] where the ripple of the switching converter is used to estimate the gradient. It is also worth to mention the complementary approaches cited in [11]. Moreover, a supervisory strategy that uses sinusoidal ESC for cascaded PV architectures is reported in [12]. It is worth to note that MPPT implementations in [10–12] use digital devices.

Stability issues about sinusoidal ESC can be found in [13, 14]. Nevertheless, the nonlinear map of previous references is approximated by a parabola. Given that a parabolic approximation is not a good approximation of the power-voltage curve of a PV system, the stability analysis is only valid for small disturbances.

In the paper, we review the sinusoidal ESC technique. Also, we present a new stability demonstration based on Lyapunov analysis where we only assume that the nonlinear map is concave. This approach ensures the stability and therefore the reliability for all the range of operation of the PV system. The paper also presents a 100 W prototype which allows us to evaluate the effectiveness of the approach. The prototype differs from that of [10], since we use low cost analog devices. We report in detail experimental measurements that corroborate the effectiveness of the approach.

The paper is organized as follows. In Section 2, we review extremum-seeking control and specifically we analyze sinusoidal ESC taking into account stability aspects by means of a Lyapunov function. The global stability property ensures a correct behavior in front of abrupt changes in operating conditions. Section 3 describes a novel architecture for PV generation in which the MPPT circuit is based on sinusoidal ESC. In Section 4, we describe an experimental verification which allows us to assess the proposed MPPT solution. And, finally, we summarize the main ideas in Section 5.

2. Sinusoidal Extremum-Seeking Control

The objective of ESC is to force the operating point to be as close as possible to the optimum for a system described by an unknown nonlinear map with an only extremum (i.e., a maximum or a minimum).

Sinusoidal ESC principle, shown in Figure 1, can be summarized as follows: given a nonlinear input-output map, if a sinusoidal signal of little amplitude is added to the input signal , the output signal oscillates around its average value. Both sinusoidal signals will be in phase if the input signal is smaller than the maximum of the nonlinear map and in counterphase if the input signal is larger than the maximum, as it is shown in Figure 1. It can be observed that when the input signal reaches maximum, then output signal doubles its frequency; also it can be appreciated that the amplitude of the ripple of the output signal depends on the slope of the curve . Also it can be noted that, when the signal is multiplied by a sinusoidal of the same frequency and phase, the multiplier output is positive before the maximum and negative at the right side of the maximum.

Figure 1: Sinusoidal ESC principle.

The method can be implemented by the schema of Figure 2. The schema consists of the nonlinear input-output map, an integrator and a detection block, and a small sinusoidal signal that is added to the map input. The detection block function demodulates the output signal .

Figure 2: Sinusoidal ESC schema.

In the following sections, we describe the detection block function and analyze the stability condition of the schema shown in Figure 2.

2.1. Detection Block

Detection block output is proportional to the gradient of the nonlinear map as it is shown in the following paragraphs.

Given a signal at the output of the integrator block, the output of the nonlinear mapping will correspond to

Considering that the sinusoidal perturbation is small, namely, given the condition , then expression (1) can be approximated by its Taylor development as

Thus, using the trigonometric identity , the signal at the output of the multiplier block can be approximated by

Now, considering that the low-pass filter attenuates completely the first and second harmonics, the expression of the filter output can be written as where is the convolution operator, stands for the inverse Laplace transform, and represents the filter impulsional response.

Therefore, under the assumptions of small amplitude of the sinusoidal perturbation and enough attenuation of the harmonics of multiplier output, it can be stated that the output signal of the detection block is proportional to the derivative of the nonlinear map output with respect to its input, namely, the gradient .

2.2. Dynamical Stability Analysis

The following analysis is carried out by means of averaged signals. Averaged signals and real signals differ only in a magnitude which depends on , and we have assumed that the sinusoidal perturbation is small.

Averaged signals that take part in the sinusoidal extremum-seeking circuit are , , and , being . The averaged output of the nonlinear map is , and, from expression (2), we can state that , and, thus, dynamical relations of averaged signals can be represented by Figure 3.

Figure 3: Averaged-signal extremum-seeking block diagram.

Therefore, denoting the constants terms as , that is, , the averaged dynamical behavior can be described as where the first row indicates the integrator differential equation and the second row the filter differential equation.

We assume that has an only maximum and is a concave function; thus the second derivative or Hessian is negative, that is, .

It can be observed that the time derivative of the gradient corresponds to

Now, renaming the state variables as and , we can express the averaged dynamical system as being .

Then, we consider the following candidate Lyapunov function to prove the stability of the nonlinear system (7),

First, in order to verify the positive definiteness of , we note that the existence of a continuous function such that , for all, implies that . It can be corroborated given that it is true for positive and negative , that is,

Now taking the function as , it can be seen that, since , , and , then , for all. Therefore, , and, thus, is positive definite.

In addition the time derivative of corresponds to

And substituting the state variables derivatives according to (7) in the time derivative of the Lyapunov function, (10) can be rewritten as

Hence, given that is positive definite and is seminegative definite, the stability of the system is proven. It is also straight to establish the asymptotical stability of the system (7), by using the Lasalle invariance principle [15] since the only invariant of system (7) for which is the origin.

We should remark that the only assumption is that the nonlinear map has an only máximum, that is, . The nonlinear map corresponding to the power-voltage curve of a PV panel has an only maximum. In the following section, we adapt the analyzed algorithm as MPPT circuit for a PV system.

Moreover, we remark this prove ensures the stability and therefore the reliability for all the range of operation of the PV system. An analysis that only considers the stability in front of small signal disturbances can be found in [14].

3. Electronic Architecture of the Sinusoidal ESC as MPPT Circuit

The solar panel model, as that depicted in Figure 4, can be represented by the function that relates the output power of a solar array with respect its port voltage . The function is a concave with uniform irradiation, namely, . It means that the function has a unique maximum at MPP.

Figure 4: Proposed MPPT PV generator schema.

We adapt the PV panel output voltage to the load (in our case a battery) by means of a boost DC-DC switching converter, as shown in Figure 4. In addition, the relation between the duty-cycle and the input voltage in the boost converter corresponds to for in a time scale much slower than its switching frequency and its time constants. Then, it can be stated that is a monotonically increasing function inside the duty-cycle interval, and the composite function is a concave function [16], namely, . Therefore, signals in the schema of Figure 4 tend to the equilibrium, and the equilibrium corresponds to the maximum output power point, as stated in previous section.

4. Experimental Verification

In this section, we describe an electronic implementation of the schema depicted in Figure 4. We have built a 100 W prototype to verify the proposed approach. The prototype features are first described; then, we depict the obtained waveforms.

4.1. Prototype Description

Figure 5 illustrates the experimental prototype based on sinusoidal ESC. The DC-DC boost-like converter processes the power from an 85 Wp monocrystalline photovoltaic array, whose nominal open-circuit voltage is 22.1 V and the nominal voltage value at the maximum power point is 18 V. The power is stored in an acid-lead battery whose nominal voltage is 24 V. The capacitor values are 2 μF, and inductor value is 30 μH in the dc-dc converter, as Figure 5 depicts. The converter switching frequency is 300 KHz. Power measures use an IC AD633 analog multiplier.

Figure 5: Prototype of PV generator with MPPT based on sinusoidal ESC.

Besides the PV module, the dc-dc converter and its PWM circuit, the battery, and the power calculating multiplier, there exists also an adder block which sums the sinusoidal dithering before the PWM block. The key parameters of the sinusoidal dithering generator are its amplitude and frequency. The generator amplitude must be much larger than the switching ripple but small enough since MPPT efficiency depends on it. In the prototype, the generator amplitude and frequency are, respectively, 0.5 Vpp and 2 kHz, whereas the PWM sawtooth amplitude is 7 Vpp. The demodulating or detecting multiplication is also carried out by the IC AD633. A first-order filter is located after the detecting multiplier whose cut-off frequency is 200 Hz aprox. A tradeoff appears in the election of the crossover frequency since a higher value allows a faster behavior; nevertheless a lower value filters undesired harmonics better. The integrator block values must ensure that no saturation appears, and the controller gain must ensure stability and a suitable degree of damping. Their values are shown in Figure 5.

4.2. Measured Signals

In this subsection, we explain the behavior of the prototype in steady state, also when a 3 V pulse is put in series with photovoltaic panel, and when a 6 V pulse is put in series with the battery. Figure 6 shows the steady-state waveforms. Given maximum power and the mean power displayed in the figure, it can be derived that the MPPT efficiency is ; a better efficiency can still be achieved, but we show cleaner waveforms with this efficiency level. It can also be noted that, during a half period of the voltage signal , the power signal passes through the maximum; therefore, the system is effectively around the optimal point. Battery voltage is 23.9 V, and mean values of voltage and current are 14.95 V and 2.52 A, respectively.

Figure 6: Steady-state waveforms.

Figure 7(a) depicts the behavior when a voltage pulse appears in series with the PV panel. Signals from top to down are the following: panel output voltage , battery voltage , panel output current , integrator output , and power . It can be appreciated that when the system loses the MPP due to a pulsed variation, the system recovers its steady state very fast. Figure 7(b) shows a detail of the transition when a voltage pulse is added to the panel voltage.

(a)

(b)

Figure 7: (a) Transient waveforms for voltage-pulsed signal in series with the PV panel. (b) Detail of transient for voltage pulse added to panel voltage.

Figure 8(a) shows the transient in front of a pulse variation of the battery voltage. The upper signal in the oscilloscope picture is the panel output voltage which recovers a value of 15 V after the transient. The second signal depicts the battery voltage when a pulse variation of 6 V appears. The third signal corresponds to the power signal , and it recovers the MPP in less than 20 ms. The fourth signal is the panel output current , whereas the lower signal shows the integrator output . As can be viewed, in Figure 8(b), after a step decrease of the battery voltage, the sinusoidal oscillation of the power seems to fold near the MPP, changing in 180° the phase of the main harmonic of the signal.

(a)

(b)

Figure 8: (a) Transient waveform when voltage pulse signal is added to battery voltage. (b) Detail of the transient waveforms during a decrease step of the battery voltage.

5. Conclusions

The paper reviews the control technique named sinusoidal ESC. A novel architecture for MPPT photovoltaic generation is proposed based on this technique. Such technique provides a high efficient method to track the maximum power point, since it results in a small oscillation around the MPP. This oscillation depends on the amplitude of the sinusoidal modulator signal and is proportional to the slope of the power-voltage curve of the PV panel. Therefore, given that the slope is small near optimum, then the oscillation is small near the maximum and the MPPT has a very good performance. The paper also analyzes the stability of sinusoidal ESC for PV systems by means a Lyapunov function. In addition, in order to verify the effectiveness of the approach, we describe the implementation a prototype for a 100 Wp solar generator, which provides accurate measures which are in good agreement with previous derivation.

References

S. Sladić, S. Skok, and D. Nedeljković, “Efficiency considerations and application limits of single-phase active power filter with converters for photoenergy applications,” International Journal of Photoenergy, vol. 2011, Article ID 643912, 8 pages, 2011.

N. Onat, “Recent developments in maximum power point tracking technologies for photovoltaic systems,” International Journal of Photoenergy, vol. 2010, Article ID 245316, 11 pages, 2010.

C. Cabal, Optimisation énergétique de l'étage d'adaptation électronique dédié à la conversion photovoltaïque, Thèse de Doctorat, Université de Toulouse, 2008, (In French).

N. Femia, G. Petrone, G. Spagnuolo, and M. Vitelli, “Optimization of perturb and observe maximum power point tracking method,” IEEE Transactions on Power Electronics, vol. 20, no. 4, pp. 963–973, 2005.

K. H. Hussein, I. Muta, T. Hoshino, and M. Osakada, “Maximum photovoltaic power tracking: an algorithm for rapidly changing atmospheric conditions,” IEE Proceedings: Generation, Transmission and Distribution, vol. 142, no. 1, pp. 59–64, 1995.

Subiyanto, A. Mohamed, and H. Shareef, “Hopfield neural network optimized fuzzy logic controller for maximum power point tracking in a photovoltaic system,” International Journal of Photoenergy, vol. 2012, Article ID 798361, 13 pages, 2012.

M. Taherbaneh, A. H. Rezaie, H. Ghafoorifard, K. Rahimi, and M. B. Menhaj, “Maximizing output power of a solar panel via combination of sun tracking and maximum power point tracking by fuzzy controllers,” International Journal of Photoenergy, vol. 2010, Article ID 312580, 13 pages, 2010.

I. S. Morosanov, “Methods of extremum control,” Automation and Remote Control, vol. 18, pp. 1077–1092, 1957.

R. Leyva, C. Alonso, I. Queinnec, A. Cid-Pastor, D. Lagrange, and L. Martínez-Salamero, “MPPT of photovoltaic systems using extremum seeking control,” IEEE Transactions on Aerospace and Electronic Systems, vol. 42, no. 1, pp. 249–258, 2006.

S. L. Brunton, C. W. Rowley, S. R. Kulkarni, and C. Clarkson, “Maximum power point tracking for photovoltaic optimization using ripple-based extremum seeking control,” IEEE Transactions on Power Electronics, vol. 25, no. 10, Article ID 5460938, pp. 2531–2540, 2010.

A. M. Bazzi and P. T. Krein, “Concerning maximum power point tracking for photovoltaic optimization using ripple-based extremum seeking control,” IEEE Transactions on Power Electronics, vol. 26, no. 6, Article ID 5959287, pp. 1611–1612, 2011.

A. I. Bratcu, I. Munteanu, S. Bacha, D. Picault, and B. Raison, “Cascaded DC—DC converter photovoltaic systems: power optimization issues,” IEEE Transactions on Industrial Electronics, vol. 58, no. 2, Article ID 5416284, pp. 403–411, 2011.

K. B. Ariyur and M. Krstić, Real-Time Optimization by Extremum-Seeking Control Hoboken, Wiley-Interscience, Hoboken, NJ, USA, 2003.

R. Leyva, P. Artillan, C. Cabal, B. Estibals, and C. Alonso, “Dynamic performance of maximum power point tracking circuits using sinusoidal extremum seeking control for photovoltaic generation,” International Journal of Electronics, vol. 98, no. 4, pp. 529–542, 2011.

H. K. Khalil, Nonlinear Systems, Prentice-Hall, Englewood Cliffs, NJ, USA, 2nd edition, 1996.

S. Boyd and L. Vandenberghe, Convex Optimization, Cambridge University Press, Cambridge, UK, 2004.