This study presents a novel search algorithm of maximum power point tracking for photovoltaic power generation systems. The I-V characteristics and the P-V power output under specific irradiation and temperature conditions are simulated. The performance of the algorithm under fully shaded and sudden partially shaded conditions as well as variable insulations levels is investigated. The developed algorithm performs a wide-range search in order to detect rapidly changing weather conditions, and keeps the simulated stand-alone or grid-connected systems continuously operating close to the maximum power point. The performance of the developed algorithm, under extremely changing environmental conditions, is found to be superior compared to that of other conventional algorithms. The results of this study show that, under uniform radiations conditions, the developed algorithm takes only half of the time required by the Perturbation and Observe algorithms to reach maximum power point MMP. Furthermore, when PV is subjected to sudden partial shading conditions, the algorithm rapidly detects these changes and reaches the new MMP in less than a second.
1. Introduction
It is now widely accepted that the nonrenewable sources in the world are finite and it is only a matter of time before reserves will essentially be consumed [1, 2]. It has been proven that the use of nonrenewable energy sources has severe effect on the environment. Due to environmental awareness and technological advancement, high oil price, and government support, the renewable electricity generation capacity has reached an estimated 240 gigawatts (GW) worldwide in 2007 while it was 160 GW in 2004 [3]. The solar photovoltaic (PV) power system is attractive renewable energy source due to its availability and economic feasibility [4, 5]. Stand-alone PV systems are found suitable for powering remote areas [4].
The power produced by a PV module depends on the operating temperature, the amount of falling solar irradiance over the PV Cells array, and the load connected [5, 6]. The power output of PV cells depends on the nonlinear current voltage (I-V) characteristics relationship. Because of this nonlinear relationship between the current and the voltage of the PV cell, there is a unique maximum power point at particular weather conditions, and this maximum power point keeps changing with the irradiance levels and ambient temperature. Therefore, a maximum power point tracking (MPPT) algorithm is commonly used, to obtain the maximum possible power under varying weather conditions and loads. Because of the nonlinear I-V relation, the power versus voltage P-V relation has more complicated behavior especially when the weather conditions change. This complex behavior makes analytical solution very difficult and force researchers to develop numerous techniques for finding MPP. MPPT algorithms are implemented using digital signal processors (DSP) or microcontrollers. They control DC/DC converters to drive the PV panel to operate at their MPP.
The main methods used to achieve MPPT for PV cells are Estimation Methods, Heuristic Method, and Search Algorithms. System modeling method [7–9], curve-fitting method [10, 11], open circuit method [12, 13], and short circuit method [14, 15] are examples of Estimation Methods. They run PV cells to track an Estimated Maximum Power Point rather than the Actual Maximum Power Point. In general, Estimation Methods depend on an approximated mathematical model to calculate an estimated MPP. PV cell current-voltage data, irradiance, and temperature levels are the required inputs for these methods. The main advantages of Estimation Methods are their simple implementation and fast response. On the other hand, they are expensive and inaccurate, require the use of many sensors, demand large computational power, and fail under rapidly changing atmospheric conditions.
Heuristic methods are recently developed to overcome the problem associated with the inaccuracy of the PV cell mathematical model. The succesful development of these methods is attributed to the recent advances in nonlinear control method. Fuzzy Logic Control methods [16, 17], Neural Network methods [18–22] and Genetics Algorithm methods are examples of Heuristic Methods developed to achieve MPPT of PV cells. The Fuzzy Logic Control method effectively tracks the MPP under various weather conditions. However, the performance highly depends on the expertise of the rule-based system designer which might lead to the failure of the controller in tracking the MPP under partial shading condition. The outcomes of MPPT using Neural Network methods are highly related to the accuracy and efficiency of the designed algorithm, the size of the training database, and the network training quality. Furthermore, they require gathered data for various conditions and multiple locations to guarantee a better performance.
Unlike estimation methods, Search Algorithms track the actual MPP rather than an estimated MPP. However, they continuously search for the MPP by increasing or decreasing the PV cell output voltage. Many methods are developed implementing the search algorithms concept. Some of those methods are the differentiation method [23, 24], the Perturbation and Observe methods, or the modified Perturbation and Observe methods [25–30], the Incremental Conductance method [31], the Parasitic Capacitance method [32], the Fibonacci Search method [33, 34], the Slide Control method [35], and the Particle Swamp Optimization method [36, 37]. The main advantage of these methods can be summarized as follows. No previous knowledge about the PV cell characteristics is required, no or minimal use of sensors is required to show acceptable response to rapidly changing conditions. On the other hand, search algorithms waste energy as they continuously oscillate around the MPP. Unfortunately, solving the oscillation problem requires complex computational power. Furthermore, the performance of these methods degrades significantly under partial shading conditions and fails under sudden partial shading conditions. Unfortunately, none of the developed search algorithms can successfully deal with all extreme weather conditions, namely, rapidly changing, partial shading, and sudden partial shading conditions.
All the developed methods found in the literature still suffer from weak performance under extreme weather conditions. It is the objective of this study to develop an efficient novel search algorithm that deals with rapidly changing, partial shading and sudden partial shading conditions to track the MPP of PV cells. Minimal user interface for wider adaptation to different solar panels and different environmental conditions is considered in developing the proposed search algorithm.
2. PV Cells Model
The modern solar cells are fabricated from a p-n junction and prepared in a small thickness semiconductor layer. These cells create electric current when subjected to sunlight. The PV cells act as a form of diode in which the parameters of this diode define the circuit model. Walker [38] proposed a simple approximate circuit for the PV cells as shown in Figure 1. He modeled the PV cells as a constant current source connected in parallel with a diode. The model includes the necessary series resistance but neglects the shunt resistance. The characteristics equation of the current I produced by the PV cell is proportional to the intensity of the radiation fallen on the PV cells [38] as follows:
(1)I=Ipv-I0[eq(V+IRs)/akT-1],
where Ipv is the current (A) generated by the photovoltaic cell due the incident light, k=1.3806503×10-23 (J/K) is the Boltzmann’s constant, T is the temperature in Kelvins (K) at the p-n junction, q=1.60217646×10-23 (C) is the electron charge, I0 is the diode saturation or leakage current (A), V is the PV cell output voltage (Volt), and a is the diode ideality factor.
Simple PV equivalent circuit model.
Villalva et al. suggested adding several parameters in order to capture the behavior of actual panels which consist of a number of connected photovoltaic cells [39]. Figure 2 demonstrates the equivalent circuit of a photovoltaic cell with the additional parameters for the improved model. The improved I-V characteristic of a PV array is given as [39]
(2)I=Ipv-I0[eq(V+IRs)/aVt-1]-V+RsIRp,
where Vt=(NskT/q) is the thermal voltage of the array with Ns cells connected in series and Rs and Rp are the equivalent series and shunt resistances (Ω) of the array, respectively. The series resistance is the sum of several internal and structural resistances within the device. The parallel resistance models the current leakage through the p-n junction.
Single-diode model of the theoretical PV with additional parameters for the improved PV model.
Figure 3 shows a MATLAB/Simulink model for the improved PV panel mathematical model. Figure 4 shows a MATLAB/Simulink model of the PV circuit built to obtain the I-V characteristics according to (2). The photogenerated current is function of solar radiation and temperature of the p-n junction is given by [39] as follows:
(3)Ipv=(Ipv,n+kiΔT)GGn,
where Ipv,n is the photogenerated current at the nominal conditions (25°C, 1000 W/m2), ΔT=T-Tn, Tn is the nominal temperature in kelvin, G is the radiation fallen on the device, and Gn is the nominal radiation. Figure 5 shows the Simulink model for calculating Ipv. The saturation current equation derived by [39] is given as
(4)I0=Isc,n+kiΔT(e((Voc,n+KvΔT)/aVt)-1),
where Isc,n is the short circuit current at nominal conditions and ki and kv are the current and voltage coefficients, respectively. Voc,n is the open circuit current at nominal conditions. Figure 6 shows the Simulink model used to calculate saturation current I0. The values of Isc,n, ki, kv, and Voc,n are supplied by the PV panel manufacturer. Villalva et al. suggested a method to adjust the values of Rs and Rp based on the fact that there is only pair {Rs,Rp}, which makes the simulated maximum power point Pmax,m=Vmp×Imp, equal to the experimental maximum power point Pmax,e [39]. The method yields the following equations:
(5)Pmax,m=Vmp{Vmp+RsImpRpIpv-I0[exp(qKTVmp+RsImpaNs)]-Vmp+RsImpRp},Rp=Vmp(Vmp+ImpRs)×({VmpIpv-VmpI0[exp(qKTVmp+RsImpaNs)]+VmpI0-Pmax,eexp(qKTVmp+RsImpaNs)})-1.
MATLAB/Simulink model of improved PV panel.
MATLAB/Simulink model for calculating the I-V characteristics.
PV cell MATLAB/Simulink model for calculating Ipv.
MATLAB/Simulink model for calculating I0.
From the above equations, it is obvious that output power of a PV module depends on the solar irradiance values and ambient temperature.
3. Partial Shading Problem
The temperature, the irradiation levels, and the shading of the system affect the performance of photovoltaic cells. Partial shading problems arise due to the existence of clouds or building shadows. This problem makes the photovoltaic power characteristics more complicated with multiple peaks in power. This reduces the efficiency of most MPPT techniques. The effect of partial shading problem appears significantly for large arrangement of panels. Under partially shaded condition, it has been found that the I-V curves have multiple stairs while the P-V curves have multiple peaks [40].
4. Boost Converter Design
In this paper, we utilized a boost convertor to change the PV panel operating point to its MPP. The operation of the DC-DC converters is controlled by the MPPT algorithm making the power output of the panel operate at its the maximum level. The MPPT algorithms are usually implemented using either digital signal processors (DSP) or a microcontroller.
DC-DC boost converter shown in Figure 7 transforms an unstable VI voltage source into a higher-level stable output voltage VO. Over the past years, these converters have shifted from the conventional analog control to pulse width modulated PWM digital control. These converters use solid-state components including MOSFETS, transistors, and diodes to operate a digital switch (on-off switching) to control the resultant output of these DC-DC converters. These converters employ capacitive and inductance elements to store and transfer energy and eliminate the noise (lowpass filter). Unlike analog control methods, the digital control methods of DC-DC converters are quite resilient and flexible in changing the software. Also, they offer the advantage of implementing more complex control algorithms.
Circuit diagram of boost converter.
Boost converters have two modes of operation. The closed-switch mode starts the diode reverse bias mode. It causes the input power source to store energy in the inductor as well as the capacitor to discharge into the storage battery, while the open-switch mode starts the diode forward bias mode causing higher output energy supply from both the power source and the inductor to the capacitor and the load. The equation defining the ratio between the input voltage and the output voltage is given as
(6)VinVout=11-D,
where D=Ton/(Ton+Toff) is the duty cycle of the input PWM signal. The microcontroller sends the proper PWM signal to the boost converter switch to reach the MPP advised by the MPPT algorithm.
In order to design an appropriate boost converter, the following steps are carried out.
The boost converter specifications listed in Table 1 are selected to satisfy continuous operation under all conditions.
A 200 kHz switching frequency is chosen to reduce the size of the boost converter components and decrease the power loss.
The boost inductor value plays a key role in determining the operational mode of the system. In order to have continuous operation mode, the inductance should satisfy the following equation:
(7)L≥R×ts2×(M-1)M3,
where R is the load resistance, L is the inductance, M is the maximum voltage gain which is equal to (Vout/Vin), and ts is the switching period. We chose the closest available conductor value of 68 μH.
The maximum duty cycle is calculated according to (6) and is found to be 0.846.
The capacitance value is given by
(8)C=Ipk2×L2×ΔV×(V0-Vi),
where Ipk is given by
(9)Ipk=ViL×D×ts=866×10-6×0.846×12×10-6=0.5A.
The allowed ripple in the output voltage is required to be 50 mV. Hence, the desired capacitance value obtained from (8) is as
(10)C=0.52×662×0.05×(13-2)=153μF.
The equivalent series resistance needed to limit the output ripple to 5 mV is calculated by
(11)ESR=ΔV0ΔI0=0.050.5=10mΩ.
Specification
Min
Max
Units
Input voltage
0
30
V
Output voltage
50
51
V
Output power
0
150
W
Operation frequency
1
MHz
Voltage allowed output ripple
50
mV
To sooth the signal, almost twice of the calculated capacitor value is used, that is, 330 μF. Using all the design calculated and chosen specifications, the Simulink model of the PV panel system with the boost converter is shown in Figure 8.
Full simulator model.
5. Theory of the Proposed MPPT Algorithm
There are several factors to consider when developing and choosing the techniques for performing MPPT, such as the ability of an algorithm to detect multiple maxima, costs, and convergence speed. MPPT is naturally a maxima-finding process. The proposed Maximum Power Point tracking algorithm implements the search algorithms concept. The reasons behind this choice as mentioned previously are no previous knowledge about the PV cell characteristics is required, simple implementation, and guaranteed convergence. The main disadvantage of search algorithms are that; they waste energy as they continuously oscillate around the MPP and they show inadequate response under partial shading conditions and fail under sudden partial shading conditions.
Most search algorithms model the data as a 1-D function and go about a Brute-Force method of finding the maxima of the function. These kinds of algorithms require a large amount of processing time. Other algorithms like the Shubert algorithm [40–42] rely on Lipschitz continuity. They use weighing parameters that emphasize on the local search versus the global search for the optima. However, these constants may not exist or could not be easily computed especially for optimizing nonlinear control system which is the case in this study. Also, these constants are required to be large enough to exceed the rate of change of I-V curve. This might lead to a large number of iterations as the rate of convergence towards the optimal point slows down. Regardless of these draw backs, the Lipschitz method remains highly attractive due to the ability to bound the rate of change of the function; thus, searching algorithms can be easily implemented and one parameter is required to be specified, that is, Lipschitz constant [41].
Eliminating the need to specify the constant and making the algorithm consider both local and global search are the criteria for developing the new algorithm. In this study, we followed Jones et al. methodology [41] by utilizing the advantages of Shubert algorithm, mainly the bounding of the rate of change. Jones et al. presented an algorithm called DIRECT algorithm, that is, a modification of the standard Lipschitzian to overcome the problems of normal Shubert’s algorithm [40] (Figure 11). In order to clarify the logical development and the features of the proposed algorithm, we begin by reviewing Shubert’s method and discussing its main drawbacks; then we will present Jones’s algorithm that eliminates the need to specify a Lipschitz constant. Finally, we present the newly developed algorithm.
5.1. Lipschitzian Optimization
The goal is to find the maximum functional value of f(x). The normal Shubert method can be summarized as follows.
Lipschitz continuity states that a function f(x) defined on the closed interval [l,u] is called Lipschitz continuous on [l,u] if there exists a positive constant, the Lipschitz constant, such that
(12)|f(x)-f(x′)|<α|x-x′|,∀x,x′∈[l,u].
Let us take a hypothetical function f(x) defined on [a,b]. If we substitute a and b for x′ into the definition of Lipschitz-continuity, we get the following two inequalities for f(x), where x∈[a,b],
(13)f(x)>f(a)-α(x-a),f(x)>f(b)+α(x-b).
The inequalities (13) form V-shaped formed from the two lines with slopes -α and +α with the intersection occur below f(x) as shown in Figure 9.
An initial lower bound for f using the Lipschitz constant.
The point of intersection for the two lines x1 is easy to calculate and is given as
(14)x1=(a+b)2+(f(a)-f(b))2α.
The initial lower bound of f is denoted by f1 and given by
(15)f1=(f(a)-f(b))2α-α(b-a).
The Shubert algorithm uses this straightforward idea in finding the minimum or the maximum of f(x). Shubert’s algorithm is an iterative algorithm that continues to perform the same operation on the regions [a,x1] and [x1,b], eliminating the higher f1 value intervals to finally reach the global minimum value as demonstrated by Figure 10.
Iterations of the Shubert algorithm in dividing the intervals of minimum f1 [41].
Dividing strategy of DIRECT algorithm.
Evaluating the function at the center of any interval rather than the bounds of the interval is the main idea of DIRECT algorithm developed by Jones et al. [41]. Mathematically, this can be expressed as
(16)f(x)>f(c)+K(x-c),forx≤c,f(x)>f(c)-K(x-c),forx≥c,
where c=(a+b)/2. Thus, the lower bound equation has to take into account the function value at the center of the interval:
(17)Lowerbound=f(c)-α(b-a)2.
Figure 13 shows the interval-dividing strategy of the DIRECT algorithm when a sampling interval [a,b] has been specified. Assume that the algorithm has already taken the sample c at the center of [a,b] in the previous step. This interval is then divided into three intervals [a1,b1], [a2,b2], and [a3,b3] resulting in two new center points to be evaluated c1, c2. The sample c simply becomes the center of the new middle interval. The algorithm then evaluates the three samples to decide the next sampling interval. It is clear that only two new samples in each dividing iterations are required for evaluation. Further subdivision for the potential interval containing optima is carried out until the optimal point is found. Theoretical details are found in reference [41]. The fact that in convex hull functions local optima are global optima is used to select potentially optimal interval. Suppose that we have partitioned the interval [l,u] into intervals [ai,bi] with midpoints ci% for i=1,…,m. Let ε>0 be a positive constant and let fmax be the current best function value. Interval j is said to be potentially optimal if there exists some rate-of-change constant K~>0 such that [41]
(18)f(cj)+K~(bj-aj)2≥f(cj)+K~(bi-ai)2,∀i=1,m,(19)f(cj)+K~(bj-aj)2≥fmin+ε|fmin|.
The inequality (18) selects intervals that would improve the current function value. For intervals with the same length, the interval with the highest function value at its center point is chosen to be the potentially optimal interval (POI). The inequality (19) ensures that the POIs exceed the current best solution by a nontrivial amount ε|fmax|. Figure 12 demonstrates how convex hull sets help choose POIs that satisfy both (18) and (19) [41]. If we construct convex hull from the function values at the center points, the intervals that make up convex hull are considered POIs. Graham’s Scan is efficient algorithm used to create a convex hull out of the set of center points [43]. Graham’s scan is a phase algorithm that can be summarized as follows, given a set of points Sp.
Find the point in Sp with the maximum value. If two or more have the same value, use one with the lowest x coordinate. Call it P0.
Calculate the angles in radians that each of the points makes with P0; then sort them in increasing order, and push them onto a stack.
If P0 forms a left turn with the last two points in the stack, we push P0 onto the stack; else, we discard and make the next point in the stack P0 and repeat.
Repeat step number (3) until you encounter P0 again.
Set of potentially optimal intervals.
Change in power under partially shaded condition identification in (a) I-V and (b) P-V.
A simplified approach is to calculate the direction cross product of the two vectors formed from three points P0-P1 and P1-P2. If the value is positive, it is a left turn and thus we keep the point and the interval. If it is negative, then we discard the interval all together.
5.2. Golden Section Search (for Rapidly Changing Conditions)
The golden section algorithm is used to detect the environmental change by continuously oscillating around the maximum power point. The Golden Section Search method is used to find the maximum or the minimum of a unimodal function by calculating the function at three different points. In this study, Golden Section Search (GSS) MPPT algorithm uses the voltage as the search variable. The main advantage of GSS algorithm is its fast convergence compared to many other MPPT algorithms. The MPPT algorithm is developed with the limiting parameters for fast convergence. The main steps in GSS algorithm are as follows.
Initialization
Determine xl and xu which is known to contain the maximum of the function f(x).
Determine two intermediate points x1 and x2 such that
(20)x1=xl+5-12(xu-xl),x2=x2-5-12(xu-xl).
If f(x1)>f(x2), then new points xl, x1, x2, and xu are updated as
(21)xl=x2,x2=x1,xu=xu,x1=xl+5-12(xu-xl).
If f(x1)<f(x2), then the new points xl, x1, x2, and xu are updated as
(22)xl=xl,xu=x1,x1=x2,x2=xu+5-12(xu-xl).
If xu-xL<ɛ (a predefined condition), then the maximum occurs at (xu+xL)/2 stop iteration; else, go to step 2.
The intermediate points x1 and x2 are chosen such that the ratio of the distance from these points to the boundaries of the search region is equal to the golden ratio. The golden ratio which is equal to 1.61803398 makes the algorithm converge at a constant speed.
5.3. Implementation for PV System
Theoretically, VMMP falls between 0 V and Voc_max. In the present study, a DC/DC converter is used to vary the optimal point OP of the PV system. Hence, the duty cycle would be in the range of (0,1). However, this range can be controlled to be much smaller in practice. The smaller range will increase the convergence speed. The maximum duty cycle is calculated according to (6) and found to be 0.846.
As mentioned previously, the I-V curves have multiple stairs while the P-V curves have multiple peaks under partially shaded condition as shown in Figure 15. To explain the main idea of the new algorithm, assume that DIRECT has successfully found the maximum power point as shown in Figure 13(a). When weather condition changes happen, the OP will move to a different point due to the change of the I-V curve. Since the duty cycle is not changed, the power of the PV is decreased from POP to POP′ as explained in Figure 13(b). To detect environment changes on the PV arrays, GSS algorithm continuously oscillates around the current VMMP. The GSS is chosen due to its rapid local searching for optimal point. When executing the GSS algorithm within small interval around the current VMMP, the GSS has the ability within four iterations to decide whether environment changes have occurred or not. If OP is not found when GSS iterated four times, DIRECT algorithm is called to search for global OP. This will increase the response time of the proposed algorithm. Figure 14 shows the simplified flowchart of the DIRECT search algorithm incorporated with GSS.
Flowchart of the proposed search algorithm.
Masked Simulink model of PV cells.
6. Results
The simulation results are carried out using MATLAB/Simulink to validate the performance of the proposed MPPT algorithm.
6.1. Photovoltaic Model Simulation Results
The developed photovoltaic cell Simulink model built and shown in Figure 15 is used to simulate the performance of PV cells under different temperatures and irradiance levels. Figure 16 shows the I-V characteristic curve of a practical photovoltaic device under different weather conditions. It can be seen that the I-V characteristics are dependent on the levels of irradiance and the temperature of PV cell. Figure 17 shows the power-voltage variations under different weather condition. It is clear that the P-V curve has single peak that could be easily found by conventional searching method.
I-V characteristic curves of a practical photovoltaic device under different weather conditions.
P-V curves plotted for different weather conditions.
To investigate partial shading conditions, first, a masked model of a single PV cell is built as shown in Figure 18. Then, a model for a panel consisting of three PV cells connected in parallel is built as shown in Figure 19. Each cell is subjected to different solar irradiation to allow producing different values of photovoltaic current Ipv.
Masked Simulink model to calculate Ipv.
Simulink model for partial shading.
In order to clarify the complexity associated with partial shading, sample simulations are carried out and their results are shown in Figures 20 and 21.It can be seen in Figures 20 and 21 that there are two local maxima in the power-voltage curves. As mentioned previously, many MPPT algorithms are incapable of dealing with the effects of partial shading and might mistakenly drive the system to its local maximum instead of the desired global maximum.
The variation of the I-V under varying irradiance (600, 200, and 100 W/m2).
The variation of the P-V under partially shaded condition (600, 200, and 100 W/m2).
6.2. Performance under Uniform Weather Condition
The Performance and operation of the proposed search algorithm have been evaluated using MATLAB/Simulink. The sampling time is chosen to be 0.05 s. For the implemented proposed MPPT algorithms, the simulation results have been obtained during starting up of the system. The results have been obtained for a solar irradiance value of the proposed system that is tested under two uniform radiation levels: 1000 W/m2 and 2000 W/m2. As shown in Figures 22 and 23, the proposed MPPT algorithm found the global maximum in a relatively short time, that is, in less than 0.8 seconds, with small oscillation in steady state. IT is also observed that the power loss from oscillation is insignificant. To evaluate the effectiveness of the proposed algorithm, its performance is compared with that of the Perturbation and Observe algorithm [25–31]. The result in Figure 24 shows that Perturbation and Observe needed 1.9 seconds to reach the MPP.
The simulated power curves for the 1000 W/m2, 25°C fully shaded designed algorithm.
The simulated power curves for the 2000 W/m2, 25°C.
The simulated power curves for the 2000 W/m2, 25°C Perturbation and Observe [25–31].
6.3. Proposed MPPT Algorithm under Partially Shaded
The proposed scheme for MPPT algorithm is tested under partially shaded conditions. The simulations were conducted with two consecutive scenarios. In the first scenario, the PV panels are subjected to uniform insolation condition. This condition is maintained for 0.4 s before it is changed to partially shaded condition. The cells temperatures are kept constant at 25°C. Figure 25 shows that the MMPT algorithm maintains the MPP until the radiation level varies at 0.4 s. It can be seen from Figure 25 that the proposed algorithm immediately detects the weather changes and starts immediately searching for the new MPP. As it can be seen, in Figure 25, the proposed MPPT found the new MPP after 0.7 s.
MMP tracking under partial shading (600, 200, and 100 W/m2).
7. Conclusion
A novel algorithm of maximum power point tracking for photovoltaic power generation system is presented. A mathematical model of the PV panel is presented based on the theory of photovoltaic. The V-I characteristics and the P-V power output under several irradiation levels and temperature conditions are simulated. The proposed algorithm is implemented in a PV panel connected to DC-DC boost converter with resistive loading. A full Simulink MATLAB model is built to simulate the performance of the proposed algorithm. The proposed MPPT algorithm is evaluated under extreme weather conditions. The results show that various advantages are gained with the proposed scheme compared to perturbation and observe algorithms. The response time is smaller and the oscillations around the MPP were reduced to obtain steady state maximum power output. It is found that the proposed algorithm quickly identifies the MPP of the solar panels under extreme weather conditions. Furthermore, the MMPT algorithm sustains its performance when subjected to sudden changes in the insulation levels. In addition, it is capable of finding the global maximum points under both fully and partially shaded conditions.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
BatainehK.FayezN.Thermal performance of building attached sunspace in Jordan climateProceedings of the 1st International Nuclear and Renewable Energy Conference (INREC '10)March 2010Amman, Jordan2-s2.0-7795327352110.1109/INREC.2010.5462554BatainehK.Numerical simulations for average temperature differential stirling engine201323BatainehK.DalalahD.Assessment of wind energy potential for selected areas in Jordan201359758110.1016/j.renene.2013.03.034BatainehK. M.DalalahD.Optimal configuration for design of stand-alone PV system201232BatainehK.HamzehA.Efficient maximum power point tracking algorithm for photovoltaic cellsProceedings of the 1st WSEAS International Conference on Industrial and Manufacturing Technologies2013Athens, GreeceAdel HamdyM.A new model for the current-voltage output characteristics of photovoltaic modules1994501-211202-s2.0-0028423029TakashimaT.TanakaT.AmanoM.AndoY.Maximum output control of photovoltaic (PV) arrayProceedings of the 35th Intersociety Energy Conversion Engineering Conference and Exhibit (IECEC '00)July 2000Las Vegas, Nev, USA3803832-s2.0-0033711351TakeharaN.KurokamiS.Power control apparatus and method and power generating system using themPatent US 5,654,883, 1997NishiokaK.SakitaniN.KurobeK.-I.YamamotoY.IshikawaY.UraokaY.FuyukiT.Analysis of the temperature characteristics in polycrystalline Si solar cells using modified equivalent circuit model20034212717571792-s2.0-1242332881PhangJ. C. H.ChanD. S. H.PhillipsJ. R.Accurate analytical method for the extraction of solar cell model parameters198420104064082-s2.0-0021424925LaffertyD.Coupling network for improving conversion efficiency of photovoltaic power sourceUS 4,873,480, 1989ChettyP.Maximum power transfer system for a solar cell arrayUS 4,604,567, 1986MasoumM. A. S.DehboneiH.Optimal power point tracking of photovoltaic system under all operating conditionsProceedings of the 17th Congress of the World Energy Council1998Houston, Tex, USAAlghuwainemS. M.Matching of a dc motor to a photovoltaic generator using a step-up converter with a current-locked loop1994911921982-s2.0-002839888610.1109/60.282492NoguchiT.TogashiS.NakamotoR.Short-current pulse-based adaptive maximum-power-point tracking for a photovoltaic power generation system2002139165722-s2.0-003708982010.1002/eej.1147TakunP.KaitwanidvilaiS.JettanasenC.Maximum power point tracking using fuzzy logic control for photovoltaic systemsProceedings of the International MultiConference of Engineers and Computer Scientists (IMECS '11)March 2011Hong Kong9869902-s2.0-79960585041CheikhM. S. A.LarbesC.KebirG. F. T.ZerguerrasA.Maximum power point tracking using a fuzzy logic control scheme20071032387395HiyamaT.KouzumaS.ImakuboT.Identification of optimal operating point of PV modules using neural network for real time maximum power tracking control19951023603672-s2.0-002932751210.1109/60.391904RoK.RahmanS.Two-loop controller for maximizing performance of a grid-connected photovoltaic-fuel cell hybrid power plant19981332762812-s2.0-003216559510.1109/60.707608HusseinA.HirasawaK.HuJ.MurataJ.The dynamic performance of photovoltaic supplied DC motor fed from DC-DC converter and controlled by neural networksProceedings of the International Joint Conference on Neural Networks (IJCNN '02)May 20026076122-s2.0-0036079684SunX.WuW.LiX.ZhaoQ.A research on photovoltaic energy controlling system with maximum power point trackingProceedings of the Power Conversion Conference200282282610.1109/PCC.2002.997626ZhangL.BaiY.Al-AmoudiA.GA-RBF neural network based maximum power point tracking for grid-connected photovoltaic systemsProceedings of the International Conference on Power Electronics, Machines and DrivesApril 200218232-s2.0-0036039763BavaroL. T. W.Power regulation utilizing only battery current monitoringPatent US 4,794,272, 1988HuaC.LinJ. R.DSP-based controller application in battery storage of photovoltaic systemProceedings of the IEEE 22nd International Conference on Industrial Electronics, Control, and Instrumentation (IECON '96)August 1996170517102-s2.0-0030396192EnslinJ. H. R.WolfM. S.SnymanD. B.SwiegersW.Integrated photovoltaic maximum power point tracking converter19974467697732-s2.0-003133157110.1109/41.649937Al-AmoudiA.ZhangL.Optimal control of a grid-connected PV system for maximum power point tracking and unity power factorProceedings of the 7th International Conference on Power Electronics and Variable Speed DrivesSeptember 199880842-s2.0-0032311716KasaN.IidaT.IwamotoH.Maximum power point tracking with capacitor identificator for photovoltaic power systemProceedings of the 8th International Conference on Power Electronics and Variable Speed DrivesSeptember 20001301352-s2.0-0033645810ZhangL.Al-AmoudiA.BaiY.Real-time maximum power point tracking for grid-connected photovoltaic systemsProceedings of the 8th International Conference on Power Electronics and Variable Speed DrivesSeptember 20001241292-s2.0-0033652013XiaoW.DunfordW. G.PalmerP. R.CapelA.Application of centered differentiation and steepest descent to maximum power point tracking2007545253925492-s2.0-4704911660510.1109/TIE.2007.899922EnriqueJ. M.AndújarJ. M.BohórquezM. A.A reliable, fast and low cost maximum power point tracker for photovoltaic applications201084179892-s2.0-7204911983310.1016/j.solener.2009.10.011HusseinK. H.MutaI.HoshinoT.OsakadaM.Maximum photovoltaic power tracking: an algorithm for rapidly changing atmospheric conditions1995142159642-s2.0-002919547710.1049/ip-gtd:19951577BrambillaA.GambararaM.GaruttiA.RonchiF.New approach to photovoltaic arrays maximum power point trackingProceedings of the 30th Annual IEEE Power Electronics Specialists Conference (PESC '99)July 19996326372-s2.0-003334327110.1109/PESC.1999.785575MiyatakeM.KounoT.NakanoM.Maximum power point tracking control employing fibonacci search algorithm for photovoltaic power generation systemProceedings of the International Conference of Power Electronics (ICPE '01)October 2001Seoul, Republic of Korea622625AhmedN. A.MiyatakeM.A novel maximum power point tracking for photovoltaic applications under partially shaded insolation conditions20087857777842-s2.0-3974909721610.1016/j.epsr.2007.05.026ZhangM.WuJ.ZhaoH.The application of slide technology in PV maximum power point tracking systemProceedings of the 5th World Congress on Intelligent Control and Automation (WCICA '04)June 2004559155942-s2.0-4444260635MiyatakeM.ToriumiF.EndoT.FujiiN.A novel maximum power point tracker controlling several converters connected to photovoltaic arrays with particle swarm optimization techniqueProceedings of the European Conference on Power Electronics and Applications (EPE '07)September 20072-s2.0-5124909145810.1109/EPE.2007.4417640ChowdhuryS. R.SahaH.Maximum power point tracking of partially shaded solar photovoltaic arrays200934102093210010.1016/j.renene.2009.01.006WalkerG.Evaluating MPPT converter topologies using a matlab PV model200121149552-s2.0-0035173089VillalvaM. G.GazoliJ. R.Ruppert FilhoE.Modeling and circuit-based simulation of photovoltaic arraysProceedings of the Brazilian Power Electronics Conference (COBEP '09)October 2009Mato Grosso do Sul, Brazil124412542-s2.0-7795000501110.1109/COBEP.2009.5347680PatelH.AgarwalV.MATLAB-based modeling to study the effects of partial shading on PV array characteristics20082313023102-s2.0-4084914255010.1109/TEC.2007.914308ShubertB. O.A sequential method seeking the global maximum of a function1972933793882-s2.0-0015401431JonesD. R.PerttunenC. D.StuckmanB. E.Lipschitzian optimization without the Lipschitz constant19937911571812-s2.0-002767853410.1007/BF00941892GalperinE. A.The cubic algorithm198511226356402-s2.0-0022189230