MPPT Based on Sinusoidal Extremum-Seeking Control in PV Generation

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.


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 International Journal of Photoenergy 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][11][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 powervoltage 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.

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 x, the output signal y oscillates around its average value.Both sinusoidal signals will be in phase if the input signal x is smaller than the maximum of the nonlinear map and in counterphase if the input signal x is larger than the maximum, as it is shown in Figure 1.It can be observed that when the input signal x reaches maximum, then output signal y doubles its frequency; also it can be appreciated that the amplitude of the ripple of the output signal y depends on the slope of the curve y = f (x).Also it can be noted that, when the signal y is multiplied by a sinusoidal of the same frequency and phase, the multiplier output g is positive before the maximum and negative at the right side of the maximum.
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 y.
In the following sections, we describe the detection block function and analyze the stability condition of the schema shown in Figure 2.

Detection Block.
Detection block output u is proportional to the gradient of the nonlinear map as it is shown in the following paragraphs.
Given a signal x at the output of the integrator block, the output of the nonlinear mapping y will correspond to y = f (x + x 0 sin(ω 0 t)). (1) Considering that the sinusoidal perturbation is small, namely, given the condition x 0 x, then expression (1) can be approximated by its Taylor development as Thus, using the trigonometric identity 2sin 2 (ω 0 t) = 1 − cos(2ω 0 t), the signal at the output of the multiplier block g can be approximated by ( Now, considering that the low-pass filter attenuates completely the first and second harmonics, the expression of the filter output u can be written as where * is the convolution operator, L −1 stands for the inverse Laplace transform, and L −1 {a/(s + a)} 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 u is proportional to the derivative of the nonlinear map output with respect to its input, namely, the gradient df (x)/dx.Therefore, denoting the constants terms as k 1 , that is, k 1 = kx 2 0 /2, the averaged dynamical behavior can be described as

Dynamical Stability
where the first row indicates the integrator differential equation and the second row the filter differential equation.We assume that f has an only maximum and is a concave function; thus the second derivative or Hessian is negative, that is, d 2 f /dx 2 < 0.
It can be observed that the time derivative of the gradient corresponds to Now, renaming the state variables as z 1 = df /dx and z 2 = u, we can express the averaged dynamical system as 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 V (z), we note that the existence of a continuous function m(z 1 ) such that m(z 1 ) • z 1 > 0, for all z 1 / = 0, implies that First-order filter  It can be corroborated given that it is true for positive and negative z 1 , that is, Now taking the function m(z 1 ) as m(z 1 ) )z 1 , it can be seen that, since a > 0, k 1 > 0, and h(z 1 ) < 0, then m(z 1 ) • z 1 > 0, for all z 1 / = 0. Therefore, z1 0 (−ak 1 (1/h(α 1 ))α 1 )dα 1 > 0, and, thus, V (z 1 , z 2 ) is positive definite.
In addition the time derivative of V (z 1 , z 2 ) 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 V (z 1 , z 2 ) is positive definite and V (z 1 , z 2 ) 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 V = 0 is the origin.
We should remark that the only assumption is that the nonlinear map has an only máximum, that is, d 2 f /dx 2 < 0. 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].

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 P SA = f 1 (v SA ) that relates the output power P SA of a solar array with respect its port voltage v SA .The function is a concave with uniform irradiation, namely, d 2 f 1 (v SA )/dv 2 SA < 0. It means that the function P SA = f 1 (v SA ) has a unique maximum at MPP.
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 v SA = f 2 (D) between the duty-cycle D and the input voltage for D ∈ [0, 1] in a time scale much slower than its switching frequency and its time constants.Then, it can be stated that v SA = f 2 (D) is a monotonically increasing function inside the duty-cycle interval, and the composite function [16], namely, d 2 f (D)/dD 2 < 0. 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.

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 opencircuit voltage v OC 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 V BAT 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.
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.

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 the figure, it can be derived that the MPPT efficiency is η MPPT = Mean Power/ Max Power = 37.66 W/38.38 W = 98%; 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 V SA , the power signal P SA passes through the maximum; therefore, the system is effectively around the optimal point.Battery voltage V BAT is 23.9 V, and mean values of voltage and current are 14.95 V and 2.52 A, respectively.
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 V SA , battery voltage V BAT , panel output current I SA , integrator output x, and power P SA .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.
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 V SA which recovers a value of 15 V after the transient.The second signal depicts the battery voltage V BAT when a pulse variation of 6 V appears.third signal corresponds to the power signal P SA , and it recovers the MPP in less than 20 ms.The fourth signal is the panel output current I SA , whereas the lower signal shows the integrator output x.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.

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 powervoltage 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.
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 x 0 , and we have assumed that the sinusoidal perturbation x 0 is small.Averaged signals that take part in the sinusoidal extremum-seeking circuit are x = (1/T) T 0 (x +x 0 sin(ω 0 t))dt, y = (1/T) T 0 y(t)dt, and u = (1/T) T 0 u(t)dt, being T = 2π/ω 0 .The averaged output of the nonlinear map is y = f (x), and, from expression (2), we can state that f (x) f (x), and, thus, dynamical relations of averaged signals can be represented by Figure3.

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

Figure 7 :Figure 8 :
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.