1. Introduction
With the continuous increase in the need for electrical energy and the continuous shortage of fossil fuels and the impact of geopolitical problems on energy supplies and the environmental impact of excessive use, the need for renewable energies, especially the energy generated from PV cells, has increased. Most of the world’s nations realize this problem and started ambitus programs to completely rely on renewable energy sources by 2050 [
1]. Statistics indicate a significant rise in the use of PV in the production of electric energy, as the worldwide capacity of PV cells increased to 1300 megawatts, exceeding the capacity generated from wind energy by 400 megawatts [
2]. With the rapid progress in modern energy storage systems (ESS) and smart grid systems [
3,
4], the problem of intermittency in the generated power as a result of climate change has been overcome. The ESS can save the extra energy greater than the load needs and serve this stored energy when there is a deficiency in the extracted power from renewable energy sources (RES) compared to the power of the load. Moreover, the smart grid system can control the loads by different smart grid concepts to a level near the available generation from RES.
The PV systems are used to directly generate electricity from sunlight. The extracted power from the PV array is directly proportional to the light intensity, operating temperature, and the output voltage of the PV array. Connecting many modules in series and parallel are required to increase the voltage and current of the PV array. The relation between the generated power and the terminal voltage is nonlinear and it has only one peak at about 0.8 of the open circuit voltage (
Voc) of the PV array in case of uniform irradiance. In case of non-uniform irradiance falling on the PV modules, different generated power will be generated from these PV modules. For extracting the maximum power from these modules, each module should work with its optimal voltage and current which is not the case in real PV systems because modules are connected in series and parallel. This means that the current in each series branch is the same in all series modules, meanwhile, the terminal voltage of each module is different. A negative voltage may be generated at the terminal of some modules in some severe partial shading conditions (PSC). The negative voltage of shaded modules occurs when these modules act as a load on other modules due to PSC. The occurrence of the negative voltage on some of the shaded modules generates heat inside the module which may destroy it. This phenomenon is called the hot-spot phenomenon [
5]. For this reason, a bypass diode should be added in parallel with each module for hotspot protection. When the shunt diode is activated, the generated power from these modules is wasted, and the PV system loses this quantity of energy. Because of the PSC and shunt diodes, the P-V characteristics of the PV array will have fewer than or equal to the number of series modules in the PV array that have varying irradiances. The global peak (GP) is having the highest power among these peaks, whereas the other local peaks (LPs) have lower power than the GP. Several ways have been developed to track the maximum power point (MPP) during real-time operation with various PSCs. Several strategies have been introduced to track the MPP during their real-time operation with different PSCs. As a result, a dc-dc converter was utilized to track the MPP of PV systems by manipulating the power electronics switches with logic created by maximum power point tracker (MPPT) approaches. As illustrated in
Figure 1, MPPT techniques are utilized to extract the maximum power provided by a PV system by managing the on/off times of power electronic switch/switches. Some typical procedures employed the incremental change in voltage to track the MPP, such as hill climbing (HC), perturb and observe (P&O), and incremental conductance (In.Con.) [
6]. Other smart techniques, such as using fuzzy logic controllers [
7] and artificial neural networks [
8], have been used as an MPPT of PV systems, but all of these strategies fall into the conventional strategies category because they cannot track the GP and they may stick at one of the LPs in the event of PSC. As a result, typical MPPT techniques are not suggested for usage with PSC-equipped PV systems. The metaheuristic optimization algorithms (MOAs) can follow the GP and prevent the PV system from becoming caught in one of the LPs. Several MOA techniques, including particle swarm optimization (PSO) [
9], bat algorithm (BA) [
10], grey wolf optimization (GWO) [
11], and musical chairs algorithm (MCA) [
12], among others, have been employed as MPPT of PV systems. All of these MOAs have several problems, including extended convergence times, premature convergence, and particle stagnation at one of the LPs. The majority of recent studies on this subject have been proposed to overcome these challenges [
13,
14,
15,
16,
17,
18]. Still, additional efforts are needed in this sector to lower convergence time while maintaining GP tracking accuracy.
Various strategies have been used in the literature to overcome the long convergence time problem. The majority of these studies are centered on making changes to current MOAs to capture the GP quicker [
13,
14,
15,
16,
17,
18]. To overcome the random aspect of the PSO in tracking the MPP of PV systems, a deterministic approach was used to modify it [
13,
14]. The fundamental concept behind this approach is to replace the random values that should be multiplied by the acceleration factors to estimate particle velocity. The accelerated parameters are replaced by 1.0 in this investigation, and the random numbers are deleted. As a result, just the inertia weight parameter has to be adjusted. This strategy has been compared with conventional PSO and shows better performance [
13]. The main shortcoming of this optimization algorithm is the random initialization which may cause premature convergence to one of the LPs and a long convergence time that can be avoided with better initialization algorithms [
19,
20]. The strategy used in [
14] improved the random initialization of particles by initializing these particles at the predicted position of peaks. Moreover, it reduced the swarm size to reduce the convergence time. This strategy is reduced the convergence time but it should be trained for different operating voltages due to the particles using the terminal voltage, not the duty ratio [
14]. The predicted positions of peaks used in this strategy are based on the anticipated peaks placed at 0.8
Voc which is not accurate as has been discussed in [
20]. Another technique employed a linear drop in inertia weight value from 0.9 to 0.4 to increase global search at the start of optimization and improve local search at the end [
21]. This method lowered the convergence time and steady-state oscillations, but it still has to be improved. Another strategy suggested the variation of the inertia weight from 0.8 to 0.1 [
15] for the same purpose. Some other studies introduced a dynamic inertia weight in which the value of inertial weight will change based on the convergence performance [
16,
17,
18]. Another approach for linearly adjusting the acceleration parameters and inertia weight is provided [
22,
23]. All these modifications are implemented based on try-and-error bases without an optimal determination of the MOAs’ control parameters. To circumvent the use of trial and error procedures in obtaining the control parameters of MOAs, an intriguing strategy for calculating these optimal control parameters for PSO [
9] and BA [
24] is presented. In this technique, two nested optimization loops are used: the inner one to track the MPP of the PV system, and the outer one to optimize the control parameters for the internal one for the shortest time of convergence and zero failure rate. These MOAs have been used with photovoltaic systems with varying number of peaks to identify the ideal swarm size, inertia weight, and acceleration parameters. These strategies significantly enhance performance while maintaining a quick time of convergence and great accuracy.
The success of catching the GP and the convergence time will rise as the number of particles increases, and vice versa. This suggests that the time of convergence and failure rate are related to the swarm size. As a result, it is critical to choose the number of particles that provides the quickest time of convergence and zero failure rate. Some solutions employed three search agents [
25,
26], five search agents [
27], and six search agents [
28], among others. Other algorithms calculated the appropriate number of particles based on the number of peaks for the shortest time of convergence and zero failure rate [
28].
Another strategy is used to reduce the time of convergence while maintaining a zero failure rate using hybrid MPPT techniques (HMTs) [
29,
30,
31,
32,
33,
34,
35]. The idea behind the use of HMTs is the use of an effective GP searching strategy to determine its position at the beginning of the optimization, then use the fast local search and low ripple technique to accurately capture the GP. Some hybrid strategies used MOA at the beginning of optimization and conventional MPPT after that [
29,
30,
31,
32,
33]. There are other HMTs used two MOAs such as [
34,
35]. A detailed discussion of the HMT techniques is introduced in [36].
In terms of the time of convergence and failure rate, the MPPT's success depends on the initial placements of search agents in all MOAs. For the search agents, the majority of the MOAs employed random position initialization [
13]. Random initialization raises the failure rate and increases the time of convergence and should be avoided in MPPT applications. Several strategies are used to replace the random initialization by dividing the search area (voltage or duty ratio) to equal distances and initializing the search agents at these distances [
19]. This strategy is better than random initialization but still, the convergence time can be further reduced using initialization at predicated positions of peaks [
20]. This strategy has the fastest time of convergence and the lowest rate of failure than random initialization, but the swarm size should be equal to the number of peaks which may limit the flexibility of the MPPT algorithms. This point can be avoided by selecting a swarm size equal to the peaks and the rest of the particles can be randomly distributed.
Another strategy using the skipping model algorithm to reduce the time of convergence while maintaining a zero failure rate is introduced [
37,
38,
39,
40,
41]. The idea of this strategy is to avoid the search within certain values and concentrate on other areas that probably contain the GP. This strategy reduced the convergence time but it increased the calculation time which may limit the operating frequency and sampling time which consequently increases the convergence time. A detailed discussion of these algorithms is shown in [
42].
Another issue that all MOAs have when utilized as an MPPT of a PV system is termed search agent stagnation in one of the local peaks. This issue was resolved by initializing the search agents whenever the change in extracted power exceeded the present tolerance, as stated in Eqn. (1). The high value of the predefined tolerance may cause the system to be insensitive to critical changes in shading patterns and leave the search agents at one of the LPs and lose the GP, especially in gradual changes in shading patterns. Meanwhile, a low value of the specified tolerance may lead the system to reinitialize without necessity, increasing the oscillations of the PV system waveforms. The predefined tolerance is used between 5% [
43] to 10% [
44]. Some strategies avoid the dependency of re-initialization based on Eqn. (1) by re-initialization of the search agents every certain time [
45] or by using scanning search agents re-initialization at certain periods [
46,
47,
48].
where,
Pi and
Pi-1 are the extracted power from the photovoltaic system at iteration
i and
i-1, respectively,
ε is a predetermined tolerance.
1.1. Motivation
Because of the long convergence time associated with the usage of MOAs in MPPT of photovoltaic energy systems applications, researchers sought to employ novel MOAs or improve current ones. Nonetheless, the long time of convergence and high rate of failure necessitates greater work due to their relevance in PV system functioning. As a result, it is critical to assess and compare some of the most current MOAs in MPPT PV system applications with previous ones. Due to this, the dandelion optimization algorithm (DOA), a recently developed and promising optimization algorithm [
49] is introduced in this paper to evaluate its performance compared to superior MOAs used before for this purpose such as PSO [
9], GWO [
11], and MCA [
12]. Moreover, optimum initialization, optimal design of the dc-dc converter, optimal swarm size, and avoidance of search agent stagnation in LPs are tactics used to optimize the performance of MOAs when employed as an MPPT of PV systems.
1.2. Innovation and Contribution
Several MOAs have been employed in PV system MPPT applications. Several of these MOAs have shown greater performance, but additional efforts should be made to test novel MOAs to further reduce the time of convergence and rate of failure, which may be translated into an improvement in extracted power and efficiency of photovoltaic systems. As a result, the recently developed dandelion optimization algorithm (DOA) [
49] has been employed for the first time in the MPPT of PV systems. This research also provides a unique strategy for significantly reducing convergence time and avoiding search agent stagnation in LPs. The innovation and contribution involved in this paper are listed below:
Evaluation of the application of the DOA in a photovoltaic MPPT as a function of conversion time and failure rate.
Calculate the best swarm size to achieve the shortest time of convergence maintaining zero failure rate.
Evaluating the performance of the MPPT with different initialization strategies.
Using a novel strategy for avoiding the stagnation of search agents in LPs.
1.3. Paper Outlines
The remainder of this study provides a full discussion of the PV array modeling in
Section 2.
Section 3 has a full overview of the DOA and how it may be employed in the MPPT of photovoltaic energy applications.
Section 4 introduces the simulation experiments that were performed to compare the proposed DOA MPPT algorithm to alternative MOAs techniques.
Section 5 introduces the experimental work performed to validate the simulated results.
Section 6 introduces the findings of this investigation.
2. PV Array Modelling
The photovoltaic cell, which is composed of two semiconductor layers (P-N layers), is the smallest component of the PV array. The sunlight falls on the N-layer which has free electrons in its atom’s outer layer that can be easily moved from its atom if it has enough energy to move. The photon energy has adequate energy that can give this free electron the energy to move from the N-layer to the P-layer which has a free hole. The N-layer atoms turn into positive ions as the electron goes from the N-layer to the P-layer, while the P-layer atoms turn into negative ions, which might result in a voltage difference. The produced energy from the PV cell may be transmitted from the PV cell to the load after the load is linked between the P and N-layers. The PV cells should be arranged in parallel and series to get the required current and voltage of the PV modules. For the same objective, the modules should also be linked in parallel and in series. The simplest photovoltaic cell model is called the single diode model (SDM), which is the simplest way to represent the PV cell performance, is depicted in
Figure 2 [
50] which is used to represent the PV cell used in this study. Another model with higher accuracy when more than one diode shunt to the first diode to well represent the charge diffusion and recombination components charge of the PV cell [
51]. Some other studies recommend using three diodes in the PV cell model to get more accurate results [
52]. The main problem of increasing the number of diodes will increase the calculation burden of the model without a substantial improvement in the accuracy compared to the SDM [
53]. The SDM is providing adequate accuracy with a reasonable calculation burden and for this reason it is used in the modeling of this study.
From the above discussion, the PV cell can be modeled as a current generator in a shunt with a diode. The PV cell output current can be obtained from Eqn. (2) [
53].
Where,
Ig is the current source value,
K is the Boltzmann constant,
a is the diode ideality constant (
a=0.95194),
T is the temperature of PV cells (
oK).
I and
V are the terminal current and voltage of the PV modules, respectively,
Rsh and
Rs are the shunt and series resistances of the photovoltaic cell model, respectively.
The current is used to represent one PV cell shown in Eqn. (2) should be modified to model the current in the PV array as expressed in Eqn. (3).
Where,
Np and
Ns are the number of PV cells in each branch and the number of series PV cells in each branch, respectively.
The current of the source current is directly proportional to the solar irradiation and also functions in the operating temperature of the PV cell, as shown in Eqn. (4).
Where,
Ign is the light-generated current,
Tn and
Gn are the standard test temperature (25
oC) and standard solar irradiance (1000 W/m
2), respectively, and
KI is the current temperature coefficient (0.12499 %/
oC).
The diode saturation current Io can be obtained from Eqn. (5).
Where,
Eg is the semiconductor’s band-gap energy, and
I0n is the rated saturation current at standard test condition which can be obtained from Eqn. (6)
From Eqn. (5) and (6), the diode saturation current can be obtained from Eqn. (7).
Where,
KV is the voltage temperature coefficient (-0.349 %/
oC).
3. Dandelion Optimization Algorithm
Modern optimization methods must be utilized in conjunction with PV system MPPT to precisely predict the GP in a short time. The dandelion optimization algorithm (DOA) has been used in several applications, including Extreme Learning Machine (ELM) for biomedical classification problems [
49,
54], traffic flow prediction [
55], parameter estimation of PEMFCs' models [
56], the speed reducer problem of a mechanical device [
57], AVR-LFC architecture for a multi-area power system employing hybrid fractional-order PI and PIDD controllers [
58], and reactive power dispatch optimization with DG unit uncertainty [
59], and credit card fraud detection [
60]. Because the DOA performs well in these applications, it has been employed in the MPPT of photovoltaic energy systems in this research. The DOA, which was launched in 2017, was inspired by the life cycles of dandelion plants [
49]. The dandelion seeds can be spread for a long distance by wind. The structure of the seed enables it to travel with the wind that can carry the seeds due to the vortexes above it which can lift the dandelion seeds (DSs) in the rising stage. Once the rain occurs or the humidity increases, the DSs gain more weight and land in different locations. The landed seeds may be able to plant again and some others cannot plant again. The plants can plant again and will be used to generate a new generation. The same concept may be used to track the best solution to many optimization challenges. The DOA is divided into three stages: ascending, mutation, and selection. The objective is to model these three steps and apply them to find optimum solutions to optimization issues, as detailed in the following subsections. As indicated in Eqn. (8), the optimization technique is utilized to maximize the power supplied by the photovoltaic system by regulating the dc-dc converter’s duty ratio.
Where
dopt is the duty ratio corresponding to maximum power,
d is the duty ratio,
P is the extracted power from the photovoltaic energy system.
Dandelions are classified into two categories: core (CDs) and assistant dandelions (ADs). The CD has the greatest amount of power (Pmax), meanwhile, the ADs are the rest of the dandelions.
The mathematical modeling for the breeding cycle of the DSs is shown in the following subsections:
3.1. Rising Stage
Due to the vortices above the DSs, the lift force is created and it can carry the seeds for a distance depending on the wind speed and the humidity. The radius of sowing of the CD is representing the radius of the dandelions and it can be obtained from Eqn. (9).
Where
U and
L are the upper and lower duty ratio values, respectively, and
e and
g are the fade and growth factors, respectively, and
a is a factor termed the cross trend that may be calculated from Eqn. (10) [
60].
Where,
and
are the maximum power at previous and current iterations, respectively. Meanwhile,
ε is a specified tolerance to prevent a denominator value of zero.
The sowing radius of the DAs is given in Eqn. (11).
Where,
and
are the position of CD and AD of search agent i at iteration t, respectively.
ω is the weight factor used to enhance the stability of the search agents and it can be obtained from Eqn. (12) [
60].
Where, PE is the ratio of the number of calls to the goal function to the total number of calls. The total number of calling the objective function is not known since the optimization continuously works in real time. For this reason, similar values are used in [
49]. The value of the inertia factor is shown in Eqn. (12) starts with 1.0 and gradually reduced to zero when PE=PEmax and stays at zero till the end of the simulation. The re-initialization of search agents of the optimization algorithm is setting the inertia factor with 1.0 again and reducing it again with the progress of the optimization. The inertia factor enhances the effect of the previous radius of the ADs on the current radius and gradually reduces this effect and makes it depending on the difference between the positions of the CD and AD as shown in Eqn. (11).
3.2. Mutation Sowing
The ADs search particles will move toward the CD search agent and it will search for GP during their journey. A mutation approach should be employed with the CD to prevent early convergence or the ability of the search agents to become caught in one of the local peaks. This mutation strategy is done based on the Levy flight as shown in Eqn. (13).
Where, Where Levy() is a random duty ratio value derived from the Levy flight distribution with 𝛽 = 1.5 [
60].
3.3. Selection Stage
The search agents should be evaluated in terms of their fitness value in comparison to the other search agent's fitness values. Based on this assessment, a selection strategy is used to select the seeds (search agents) that will be used in the next iteration and the seeds will be removed from the search agents' swarm size. The probability of the fitness value of a certain search agent compared to the other search agents is shown in Eqn. (14), or it can be calculated from the difference between the fitness value and the average value as shown in Eqn. (15).
Reference [
49] proposes selecting search agents with low and high probabilities and removing search agents with medium probabilities to improve the DOA's exploration performance and avoid becoming caught in one of the local peaks. This technique is extremely effective at the beginning of the optimization to improve exploration, but after capturing the position of the GP, it should eliminate the search agent with a low probability to improve the exploitation of the DOA utilized in this study.
3.4. Improved DOA for MPPT of PV Systems
The suggested approach in this study is designed to improve DOA exploration and exploitation. Several solutions have been proposed in the literature to increase the exploitation performance of the MOAs, including:
1- Reducing the swarm size gradually [
61,
62,
63], where the MOA is started with a high number of search agents to increase the exploration and gradually reduces the search agents to enhance exploitation.
2- Enhancing local search pressure in which an adaptive scale factor for local search is introduced to enhance the differential evaluation’s local search [
64,
65].
3- Hybrid optimization methods utilize MOA with high exploration at the start of the optimization and MOA with strong exploitation at the end to improve exploitation performance. This method has been used with differential evolution [
55,
66,
67].
4- Dynamic variation of the control parameter, where the control parameters change during the optimization iterations [
9,
21,
24,
68,
69,
70].
The above improvement strategies have been used with the modified strategy called a guided probability-based DOA (GDOA) [
55]. In this strategy, a learning factor is introduced to learn from the CD based on the fitness value in which the highest fitness value will get a higher enhanced learning factor to enhance the exploitation performance of the DOA. Moreover, the middle search agents will be removed at the start of the optimization to improve exploration; however, after each iteration, the worst AD search agent (the one with the lowest generated power) will be removed from the swarm size in each iteration to improve the proposed algorithm's exploitation performance. The swarm size that started the simulation is called
SSmax and the minimum value of swarm size is called
SSmin. The logic used in the proposed algorism is shown in
Figure 3. The position of each search agent should be selected and the fitness values of these search agents will be determined. Moreover, the best power generated from the PV system should be compared with the previous one based on Eqn. (1). In case the condition is shown in Eqn. (1) is validated, the DOA should be reinitialized and the optimization started again due to the substantial change (
ε >0.1) in the shading patterns. Meanwhile, in case the condition is shown in Eqn. (1) is not verified, the search agents' positions should be adjusted depending on the fitness values given by the previous iteration.
The swarm size changes throughout optimization, and it can be calculated using Eqn. (16).