Abstract

This paper develops a maximum power point tracking (MPPT) algorithm to optimize photovoltaic (PV) array performance and make it more compatible to rapidly varying weather conditions. In particular, a novel extremum seeking controller (ESC), which uses a Bode ideal cutoff (BICO) filter in its structure, is designed and tested on a simulated PV array. The new algorithm is compared against the commonly used ESC MPPT algorithm with first-order filters. The BICO extremum seeking controller achieves transient rise to the MPP faster than the common extremum seeking MPPT, which is the faster and more robust method among all other methods. This claim has been discussed and proved mathematically in this paper, in addition to simulation illustrations. This faster extremum seeking algorithm enables PV systems to detect rapid variations in the environmental conditions like irradiation and temperature changes.

1. Introduction

The PV cells exhibit a nonlinear - characteristics as shown in Figure 1 and their power output mainly depends on the nature of the connected load. Since connecting the load directly to the PV system results in a poor overall efficiency, to minimize the life cycle cost of any PV system, increasing the efficiency by changing the operating point of the system using an intermediate maximum power point tracker (MPPT) is desirable.

MPPT controls the output current and voltage and consequently output power of the PV panel adaptively to maintain maximum efficiency and better performance in the presence of environmental variations. Typically MPPT algorithms are implemented on a solar array using a switching power converter; for instance, in a grid-tied inverter, the solar array charges a capacitor and then the current is switched out of the capacitor at an optimal varying duty cycle in order to extract maximum power from the PV array.

A number of solar power converter architectures with MPPT are discussed in the literature [1, 2]. As discussed in these works, convergence speed is one of the most important features among all different MPPT algorithms. Brunton et al. have pointed out in their paper that, “as irradiance decreases rapidly, the curve shrinks and the MPV and MPI decrease. If the MPPT algorithm does not track fast enough, the control current or voltage will fall off the curve” [3]. Consequently, any improvement in the rise time of MPPT improves the reliability of the system and increases the power extraction and efficiency of the whole system.

1.1. Maximum Power Point Tracking Algorithms

There are many different maximum power point tracking techniques for photovoltaic systems which are well established in the literature [411]. These techniques vary in many aspects as simplicity, convergence speed, digital or analogical implementation, sensors required, cost, range of effectiveness, and other aspects. In analog world, short current (SC), open voltage (OV), and CV are good options for MPPT; otherwise with digital circuits that require the use of microcontroller, Perturb and Observe (P&O), Incremental Conductance (IC), and temperature methods are easy to implement [12]. Table 1 and Figure 2 present the comparison among different MPPT methods considering the costs of sensors, microcontroller, and additional power components and also their efficiencies. In this table, A means absence, L means low, M means medium, and H means high.

Currently, the most popular and workhorse MPPT algorithm is P&O, because of its balance between performance and simplicity. However, it suffers from the lack of speed and adaptability which is necessary for tracking the fast transients under varying environmental conditions [3]. Recently, another adaptive algorithm, called extremum seeking control, has been developed [13] to overcome these weaknesses.

1.2. Extremum Seeking Control

A promising new robust MPPT algorithm is the method of extremum seeking (ES) control, which carries all P&O’s benefits like simplicity and performance, and in addition, improves its weaknesses [3]. Figures 3 and 4 present the comparison between ESC and P&O, where the inverter controls the current and voltage according to the set-points which are provided by these two MPPT algorithms.

The power, current, and voltage are plotted versus time for the ESC and P&O algorithms as well as the true maximum power of the system.

P&O and ESC methods oscillate closely around the real maximum power voltage, as seen in the power versus time plot. Obviously, the ESC method rises to the MPP orders of magnitude more rapidly than the P&O.

ESC MPPT has some advantages from hardware implementation point of view. Brunton et al. have mentioned in their paper that “the ripple-based ES algorithm has good MPPT performance over a range of inverter capacitor sizes. Typically, the choice of capacitor is expensive because it must be well characterized and large enough to maintain a small ripple. However, because the ES control signal exploits the natural inverter ripple, a smaller capacitor allows the tracking of rapid irradiance changes.”

ESC for peak power point tracking method has been successfully applied to biochemical reactors [14, 15], ABS control in automotive brakes [16], variable cam timing engine operation [17], electromechanical valves [18], axial compressors [19], mobile robots [20], mobile sensor networks [21, 22], optical fibre amplifiers [23], and so on. A good survey of the literature on this topic prior to 1980 can be found in [24] and a more recent overview can be found in [25]. Astrom and Wittenmark rated extremum seeking as one of the most promising adaptive control methods [26].

Since extremum seeking control has better features and performance compare to P&O which is the best known MPPT algorithm, in this paper, the improvement of better than the best MPPT algorithm has been investigated.

2. Basic Regular Extremum Seeking Algorithm

As shown in Figure 5, ESC algorithm employs a slow periodic perturbation, , which is added to the estimated signal, . If the perturbation is slow enough, the plant appears as a static map and its dynamics do not interfere with the peak seeking scheme [13]. If is on either side of , which is the optimal point, the perturbation signal , will create a periodic response of which is either in phase or out of phase with . The high-pass filter eliminates the “DC component" of . Thus and high-pass filter will be approximately two sinusoids which are in phase if or out of phase if .

The integrator approximates the gradient update law which tunes to [13]. System of Figure 5 can be summarized in mathematical equations as follows: where is the convolution operator and is the inverse Laplace transform operator. The transfer functions for HHPF and HLPF in the regular SISO ESC scheme are and respectively, where [25]. This model will be used for stability analysis in this paper.

In the following sections, after introducing BICO filter [27, 28], the advantages of using this filter in the ESC algorithm from the stability and robustness point of view will be discussed.

3. BICO Extremum Seeking Control MPPT

In this section we are going to introduce a filter which was strongly favored by Bode [29], called Bode’s ideal cutoff characteristic (BICO) filter [27, 28]. The general transfer function of low-pass BICO filter is,

Its corresponding time response in the special case (rad/s) was derived by Oberhettinger and Baddi [30], but the time response of (2) is presented as, where is the th order Bessel function.

By replacing by , high-pass BICO filter is obtained. Figure 6 compares the frequency response of high-pass BICO filter and regular high-pass filter () with the same cutoff frequency of  rad/s.

As seen in this figure, BICO filter has a sharp edge in its cutoff frequency. This great feature causes almost no attenuation for frequencies higher than and a large attenuation in the lower frequencies. Therefore the behavior of this filter is close to an “ideal” filter. This sharp edge presents in the low-pass BICO filter. By combining high-pass and low-pass BICO filters, the band-pass BICO filter with sharp edges in both sides can be obtained as well.

4. Stability Analysis of Extremum Seeking Control

4.1. Mathematical Modeling of ESC Scheme

The stability analysis of ESC algorithm has been investigated in [13, 3134]. In these literatures traditional extremum seeking control with regular first-order filters has been considered. In this paper, we have considered the similar stability analysis approach as [34] to compare the stability and robustness of BICO-ESC and regular SISO ESC. The difference between this work and [34] is that in this paper, the stability analysis has been done by considering both low-pass and high-pass filters in the ESC structure but in [34] authors have investigated the stability of simplified ESC with only low-pass filter in its structure.

According to the previous discussion, the nonlinear map in the ESC scheme is considered to be concave and is assumed to have only one extremum point. Since in the PV system applications, MPPT is employed to extract maximum amount of power from PV panels, therefore, extremum point in this case is maximum point (), where at this point , and also . In ESC algorithm, the output of the nonlinear map is, where and are the amplitude and angular frequency of perturbation signal. Since the perturbation signal assumes to be small, therefore the Taylor expansion of (4) is where H.O.T. stands for higher order terms and is the approximation of . The high-pass filtered signal will be,

High-pass filter acts as a derivative opereator in series with a low-pass filter (). By multiplying the modulation signal to the resulted signal from the high-pass filter and passing the modulated signal through the low-pass filter and since , the output signal of integrator will be

Under the assumptions that the amplitude of the sinusoidal perturbation is small and the harmonics of high-pass filter are attenuated by low-pass filter, output of low-pass filter is proportional to the gradient of the nonlinear map with respect to its input and time. Therefore in the neighborhood of the extremum point, the amplitude of the estimated signal is small, since the gradient is small. It can be seen that this amplitude depends on and .

4.2. Stability Analysis of Averaged ESC Scheme

The averaging method is typically used to analyze the periodic steady state solutions of weakly nonlinear systems. Since the amplitude of perturbation in the ESC scheme is small, this system can be evaluated by its averaged model. The averaged form of signal is, where . Therefore the averaged model of ESC scheme is

On the other hand, in the neighborhood of the extremum point, can be approximated as,

If the difference between the extremum point and averaged point is defined by and since , thus,

By defining , then, . Therefore, . Substituting this relationship in (9) gives,

Without loss of generality and by assuming that is the output signal from a high-pass filter, HHPF, which its cutoff frequency is higher than cutoff frequency of the low-pass filter, HLPF, then (12) can be rewritten as

This system can be considered as a feedback system as shown in Figure 7, which can be considered as a perturbation. The loop gain of this system depends on the demodulation gain, , the integral gain, , and the curvature of nonlinear map, . This result completely matches with the results obtained from other analysis method in the simplified case [13, 3133].

5. Analysis of ESC MPPT on PV Panels

5.1. PV Cell Model

The PV cell - curve is modeled using the light emitting diode equations [3]:

Values and definitions for further simulations are shown in Table 2.

5.2. BICO ESC vsersus Regular SISO ESC

In order to compare the qualitative behavior of BICO MPPT with the regular SISO ESC MPPT, the Root-Locus (RL) analysis has been chosen. Clearly, other stability analysis methods are applicable at this point.

Since Root-Locus (RL) method is serving for linear systems, it is important to point out that the curvature constant, , is an uncertain parameter which depends on different factors like the PV panel manufacturer and weather conditions. According to [33], .

As seen in Figure 7, the characteristic polynomial of averaged ESC system is , and to analyze the behavior of this system, Root-Locus (RL) analysis method can be employed. Since there is no command in MATLAB to plot Root-Locus (RL) for BICO transfer function, this equation has been solved for different and the roots have been plotted. To compare the behavior of regular SISO ESC with BICO ESC, the roots for this system has been plotted with the same method instead of using “rlocus” command.

Figure 8 shows the comparison between the root locus of averaged ESC using BICO and regular first-order filters. In this root locus plot, the constant gain is assumed to be and cutoff frequencies in both filters are assumed to be  rad/s. Clearly, by using first-order filters in the averaged ESC scheme, for some values of , system has complex poles which cause an oscillatory behavior in the response of the system. On the other side, by using BICO filters, roots of the characteristic polynomial (system poles) are always real numbers for any value of .

Besides the location of the poles, as can be seen in this figure, for the same value of , BICO filter has farther poles with respect to the origin compared to the first-order filter case. Therefore, the bandwidth of BICO filter is higher than the first-order filter with similar cutoff frequency. Higher bandwidth means faster response for BICO ESC. This feature is justifiable by looking at the frequency response of both filters. Since the bandwidth of BICO is close to the bandwidth of an ideal filter, BICO ESC becomes the fastest achievable MPPT algorithm.

6. Simulation Illustrations

Figures 9 and 10 show the maximum power point of a PV panel with the defined parameters in Table 2. As can be seen in these figures, environmental conditions and especially variations in the sun irradiation will change the nonlinear behavior of PV panels. Shadows, cloudy or dusty weather, and temperature variations cause moving of optimal operating point in PV systems. When temperature increases, the maximum output power of PV panels decreases and vice versa. In Figure 9 the variations of optimal operating point by changing the environmental temperature from 20°C to 30°C are illustrated. As Figure 10 presents, variation of peak output power happens in the wider range when sun irradiation changes. Generally, in PV panels, when sun irradiation increases, the output power increases, but when temperature increases, the maximum extractable power reduces.

To compare the proposed MPPT method which is called BICO MPPT with the ESC MPPT, the working conditions of both algorithms have been defined to be same. For all simulations, the ambient temperature is 25°C and the irradiation is assumed to be . Also the cutoff frequency of high-pass filters is  rad/s, and for low-pass filter, this frequency is  rad/s. Under these conditions, from Figure 9, the maximum amount of power which can be extracted from the simulated PV panel is 48 Watts and this maximum happens around 17 Volts. To implement the BICO MPPT, the discrete approximation of this filter in MathWorks, Inc.’s website has been used.

Figure 11 shows the outputs of extremum seeking algorithms resulted from two different MPPT methods. Figure 12 represents the maximum voltage tracking in BICO MPPT algorithm and ESC with first-order filters. As expected and proved before, BICO MPPT converges to the peak power point two times faster than the regular ESC MPPT algorithm. Faster convergence speed is due to the higher bandwidth of the averaged BICO system.

Figures 13 and 14 show the performance of the proposed MPPT approach compare to the ESC in the presence of white noise. As can be seen in these figures, in the presences of a white noise (noise power ), which is considered as the variations in the nonlinear map behavior, , BICO MPPT preforms better than ESC MPPT from attenuation and tracking point of view.

Another parameter which is considered in these simulations is the robustness of these two methods to the variations of the system gain. Figures 15(a) and 15(b) illustrate the performance of BICO MPPT and regular ESC MPPT to the gain variations, respectively. From these figures, it can be concluded that, BICO MPPT tolerates higher gain variations, but by increasing the gain of the integrator in ESC MPPT, it starts oscillating.

7. Conclusions

The ubiquity of using MPPT in the renewable energy systems has stimulated the persistent development of various MPPT algorithms. According to the comparison in this paper, one of the most popular methods of maximum power tracking method is the P&O method. However, this method fails in those situations when there is a need to track the MPPs rapidly.

ESC is a new robust MPPT algorithm which carries all P&O’s benefits and improves its weaknesses. In this paper, a novel ES algorithm has been presented which is called BICO MPPT. This algorithm employs BICO filter in its structure to improve the ESC rise time response even further. BICO filter is a forgotten part of Bode’s research that for the first time in this paper was used for PV applications.

BICO filter is the closest filter to the ideal filter and this feature of BICO improves the performance of ESC algorithm. Also, the advantages of BICO has been discussed in this paper. As discussed in this paper, using BICO in the ESC structure increases the bandwidth of the system, which improves the system response.

In addition, applying BICO filter to the ESC algorithm moves all the roots of characteristic polynomial of averaged BICO ESC system to the real axis. Consequently, this system has no pole with imaginary part, and therefore theoretically the system has no oscillatory response by increasing its integrator gain. These poles, however, can have imaginary part by increasing the gain of the system in the regular ESC. This feature increases the robustness of the BICO MPPT compared to regular ESC.

As can be seen in the results, BICO MPPT not only can follow the maximum power point faster than regular ESC, but also it shows more robustness in the presence of disturbance in the system or gain variations.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.