A Comparison between Particle Swarm and Grey Wolf Optimization Algorithms for Improving the Battery Autonomy in a Photovoltaic System

: This research focuses on a photovoltaic system that powers an Electric Vehicle when moving in realistic scenarios with partial shading conditions. The main goal is to ﬁnd an efﬁcient control scheme to allow the solar generator producing the maximum amount of power achievable. The ﬁrst contribution of this paper is the mathematical modelling of the photovoltaic system, its function and its features, considering the synthesis of the step-up converter and the maximum power point tracking analysis. This research looks at two intelligent control strategies to get the most power out, even with shading areas. Speciﬁcally, we show how to apply two evolutionary algorithms for this control. They are the “particle swarm optimization method” and the “grey wolf optimization method”. These algorithms were tested and evaluated when a battery storage system in an Electric Vehicle is fed through a photovoltaic system. The Simulink/Matlab tool is used to execute the simulation phases and to quantify the performances of each of these control systems. Based on our simulation tests, the best method is identiﬁed.


Introduction
Due to its availability and the good conversion factor, solar energy technology has advanced at an exponential rate in the last few years [1]. As a renewable energy source, generation from solar energy eliminates pollution caused by traditional energy industries by lowering air nocive emissions [2]. Moreover, electricity generation from this resource is quite viable for a variety of uses. In particular, the rapid development of solar energy instruments gives a complete kit of tools that can be directly applied into the field of Electric Vehicles (EV) [3,4]. Several studies have suggested that photovoltaic cells can be used to cover EVs' surfaces to store a significant amount of electricity in the storage system [5,6]. This would increase EVs' autonomy, which will in turn increase the use of EVs. Some additional benefits are also associated to solar-powered Electric Vehicles [7,8]. First, the load peaks may be reduced so that the grid management is easier. Second, a decrement of the costs of charging the EVs would also be perceived by the drivers.
In [9], researchers provided some statistics that prove that there is a huge free space to place photovoltaic (PV) cells in the car. As pointed out, these components can be used to help feed the car with electricity. It is possible to have 6 kW of electrical power in some buses or trucks with this kind of installation [10].
In order to benefit from this technology, the solar cells must be installed in a particular location on the vehicle to enhance the average autonomy of the vehicle. In the previously cited references, it is stated that the Audi car model can get around 600 W from solar equipment. If the EV model is pure electric, the PV cells can help to feed the vehicle with 1% of the total electric power when it is moving. In addition, the energy can be obtained while the EV is parked. It is proved that the PV solution is interesting for heavy trucks because there is more vacant area that can be used to generate renewable electricity. The PV system can provide 2% of their consumption when moving along a long road. However, the statistics are different in the city, as there are many shading areas where the efficiency of the PV system may be affected. The way to estimate the energy return is challenging since it depends on the vehicle's relative situation compared to the sun's light and the presence of road obstructions. Because the car can quickly change from one type of shading scenario to another, obtaining the maximum amount of energy in all these conditions is quite tough. This is even more complex if we take into account that the efficiency of PV energy production mainly depends on the dynamic conditions associated to solar energy. The most relevant features are radiation, temperature, and the state of the PV field surface (dirt, deterioration). These factors directly influence its photon absorption and therefore affect the productivity of PV panels. In addition, the phenomenon of partial shading is one of the problems that impair the proper functioning of a PV plant. Partial shading is a non-uniform distribution of illumination on photovoltaic modules, which is due to several reasons. Indeed, there are two types of shading [11,12] The distant shade corresponds to the disappearance of the sun behind the horizon line. Alternatively, close shading is often due to unavoidable obstacles such as power lines, trees, neighboring buildings or dirt.
This undesirable phenomenon affects the conversion efficiency and the ability to extract the maximum available power from the PV field by generating multiple local maxima in the PV curves. In addition, shading also disrupts the operation of PV cells, causing two problems. The first problem is the mismatch, which is due to the fact that the total current in a PV field is limited by the current of the shaded module (low power) if the current Isc (the current flowing through the photovoltaic cell in short-circuit) of the shaded module is greater than the Isc of the uniformly lit modules. The second problem is the onset of hot spot. This problem occurs when the Isc current of the shaded module is less than the Isc of the uniformly illuminated modules, so that the shaded module behaves like an energy receiver extracting energy from the other PV modules. This effect can be noticed in the PV curve of a shaded PV panel, depicted in Figure 1. In partially shaded conditions, when the PV system does not receive uniform irradiation, the P-V characteristics become more complex, with multiple local peak displayed and a single global peak. These points are referred to as Maximum Power Points (MPP). In addition to the deformation of the I-V curve, shading may lead to the heating of this module by dissipation of this energy. If the operating point of the shaded module reaches the breakdown voltage, the module will be destroyed by the avalanche effect [13], [14].
Thus, it is crucial to find an optimal operating voltage to extract the maximum amount of PV power when designing the PV system while avoiding the aforementioned problems. A Maximum Power Point Tracking (MPPT) management tool is a required In addition to the deformation of the I-V curve, shading may lead to the heating of this module by dissipation of this energy. If the operating point of the shaded module reaches the breakdown voltage, the module will be destroyed by the avalanche effect [13,14].
Thus, it is crucial to find an optimal operating voltage to extract the maximum amount of PV power when designing the PV system while avoiding the aforementioned problems. A Maximum Power Point Tracking (MPPT) management tool is a required control loop that helps to get the most power out of a PV system. The basic performance of this algorithm is to adjust the duty cycle of the power converter connecting the PV modules to the DC-link. Concerning this tracking technique, various methods were carried out and help improve the energetic performances of the system. The traditional process relies on the incremental method's principle (IC) [15,16]. It is based on the principle of the zero derivative of the output power P of the PV panel, with respect to its output voltage V at the point maximum MPP power. In the maximum point, it is positive on the left and negative on the right, so the algorithm tries to find the voltage where this condition holds [17].
The Perturb and Observe (P&O) approach was the focus of further studies [18,19]. Because of its simplicity and convenience of use, it is frequently implemented. The primary benefit of this method is its direct torque design and lesser rate of monitored parameters. However, it has a significant flaw in terms of chattering on the supplied power form.
Classic MPPT methods such as P&O and IC are based on moving the next operating point (OP) in the direction of increasing PV power. However, when partial shading is present, the P-V curve is no longer monastically increasing, as shown in Figure 1. Thus, these conventional methods can only achieve a local MPP and may not reach the global maximum [20,21]. Therefore, it is necessary to develop an appropriate MPPT algorithm that can get to the global maximum power regardless of the state of illumination on the modules.
To overcome this limitation, other strategies based on intelligent optimization, such as the fuzzy logic technique or the neural network method, have been designed in the same sector [22,23]. The two techniques lead to higher profitability, but their main issue is the database required to adapt these algorithms to PV systems. Optimization algorithms were also used to help resolve this issue of partial shading of the photovoltaic system [24,25]. In particular, the evolutionary algorithms aim to have an adaptable MPPT tracking method, based on the animal behavior to find food [26,27].
There is a variety of swarm algorithms, which have been applied in multiple systems such as in [28][29][30][31]. Among them, Particle Swarm Optimization (PSO) [25] and Grey Wolf Optimization (GWO) [32] have shown their reliability to solve real optimization problems where the objective function is not linear. In particular, the works in [17,33] only considered these two algorithms to configure a DC/DC power converter. The review presented in [34] show that PSO algorithms are still investigated to tune the power converters of microgrids. Moreover, the study elaborated by Mirjalili in [32] presents a comparison between multiple swarm algorithms. As a conclusion, they state that the better results were found for the PSO and the GWO algorithms. Indeed, the two algorithms are inspired by natural competence to reach high speed and precision. Based on these previous works, in this paper we evaluate their relative feasibility and performance of employing the swarm algorithms to configure the power converter of the PV panels in order to cope with different shading conditions. These two algorithms have already been applied and evaluated separately in PV systems. Particle Swarm Optimization algorithm can help to calculate the duty cycle of the power converter in the PV connection dynamically. Several works tested this solution for this application, as in [35] where the authors proved that this solution could be efficient if it is running offline. On the other side, the Grey Wolf Optimization (GWO) algorithm appeared as a useful solution for extracting energy from the PV system with maximum efficiency [36]. However, the two algorithms have several parameters and constants, which must be fixed initially to start the algorithms correctly.
The contribution of the paper is to apply and evaluate these two optimization algorithms for the same PV system, considering different partial shading conditions. The PV system model considered in this work is based on a commercial PV panel. So we can compare both performances to study their suitability in realistic implementations. The evaluation of each of these algorithms is based on the precision and the speed for tracking the global MPP with different partial shading conditions. Specifically, we have studied the two swarm-intelligence based algorithms for four shading conditions in a four-module PV system. This research paper is organized to initially present an introduction section, which describes the objective of the paper and explains the state-of-the-art in this application field. In the second section, the PV system is modeled with the necessary equations that regroup all the parameters and constants that define this physical system. In the third section, the two optimization algorithms are explained. Their flowcharts are exposed and the principle running of each one is described for the MPPT algorithm. In the next section, the simulation conditions and the obtained results are shown for each irradiance case. In the end, a conclusion section is formatted for resuming the paper and giving some perspectives of this work.

Model of the PV System
In this section, we first describe the model for a solar PV cell. Then, we integrate it into the model of a PV system.

A. Solar PV Cell Model
A PV module consists of several solar cells connected in series and in parallel to achieve the desired voltage and current levels. A solar panel cell is essentially a p-n semiconductor junction. When exposed to light, a direct current is generated. For simplicity, the single diode model of Figure 2 is used in this document [37]. This model offers a good compromise between simplicity and precision with the basic structure.  In Figure 2, we apply Kirchhoff's law. The current will be obtained by the following Equation: where Iph the current generated by light or photocurrent and Ip the current flowing in the parallel resistor, which can be computed as: ID is the current of the diode, which is proportional to the saturation current. The following equation expresses the value of this magnitude . .

− 1
being Isd the reverse saturation current in amperes (A), q the electron charge (1.6 × 10 -19 C), K the Boltzmann constant (1.38 × 10 -23 J/K), T the cell temperature in Kelvin (K) and n the ideal factor. We replace the voltage-current characteristic equation of a solar cell in Equation (1). So we derive that: The equivalent circuit of the general model consists of a photo current (Iph), a diode, a parallel resistance (Rp) expressing a leakage current, and a series resistance (Rs) due to the contacts between the semiconductors and the metal parts. This equivalent circuit is depicted in Figure 2.
In Figure 2, we apply Kirchhoff's law. The current will be obtained by the following Equation: where I ph the current generated by light or photocurrent and I p the current flowing in the parallel resistor, which can be computed as: I D is the current of the diode, which is proportional to the saturation current. The following equation expresses the value of this magnitude being I sd the reverse saturation current in amperes (A), q the electron charge (1.6 × 10 -19 C), K the Boltzmann constant (1.38 × 10 -23 J/K), T the cell temperature in Kelvin (K) and n the ideal factor. We replace the voltage-current characteristic equation of a solar cell in Equation (1). So we derive that: The photocurrent depends mainly on the solar radiation and the operating temperature of the cell, which is described by the following Equation: where: I sc : is the short-circuit current of the cell at 25 • C and 1000 W/m 2 K i : the temperature coefficient of the cell short-circuit current, On the other hand, the cell saturation current varies with the temperature of the cell, which is described as follows: being: Irs: the reverse saturation current of the PV cell Eg: is the gap energy of the semiconductor used in the PV cell in electron-volt (eV) n: is the ideal factor, which depends on the PV technology listed in Table 1. The reverse saturation current is given by the following Equation (7): B. Solar PV Module Model A typical PV cell produces less than 2 W at about 0.5 V so that cells must be connected in series-parallel on a module to produce enough power [38]. A PV array, like presented in Figure 3, is a group of several PV modules electrically connected in series (N s cells) and in parallel (N p columns) to generate the required current and voltage.

B. Solar PV Module Model
A typical PV cell produces less than 2 W at about 0.5 V so that cells must be connected in series-parallel on a module to produce enough power [38]. A PV array, like presented in Figure 3, is a group of several PV modules electrically connected in series ( cells) and in parallel ( columns) to generate the required current and voltage. The voltage-current characteristic equation of a PV module becomes as follows: The following Table 2 gives the parameters of the PV module used in this work.  The voltage-current characteristic equation of a PV module becomes as follows: The following Table 2 gives the parameters of the PV module used in this work.  Figure 4 shows the proposed block diagram of the PV system. In this model. We include the intelligent MPPT algorithm (based on GWO and PSO), which will be described in Sections 3 and 4. The proposed algorithms aim to compensate for the limitations of conventional maximization algorithms, which cannot find the overall maximum, especially during a partial shading phenomenon. Regardless of the partial shading profile, it can find the maximum power point on the power-voltage characteristic (Ppv) of a photovoltaic system. This characteristic is usually referred to as Ppv = f (Vpv). The algorithm is applied to the semiconductor switch of the boost converter installed to connect the PV panels to the DC-link. Specifically, the duty cycle of this device is adjusted according to the output of the Intelligent MPPT algorithm. In our implementation, V out is the voltage that will be used for feeding the battery with the required power for its charge.
algorithm is applied to the semiconductor switch of the boost converter installed to connect the PV panels to the DC-link. Specifically, the duty cycle of this device is adjusted according to the output of the Intelligent MPPT algorithm. In our implementation, is the voltage that will be used for feeding the battery with the required power for its charge.

Optimization Algorithms
This paper focuses on two swarm-intelligence based algorithms. Specifically, they are the Particle Swarm Optimization and the Grey-Wolf Optimization algorithms. They have been selected for this PV system due to its demonstrated capability to adjust power converters.

A. The PSO algorithm
This algorithm was one of the solutions that was used for resolving a huge mathematical optimization problem. It is classified as an evolutionary metaheuristic that belongs to the larger class of evolutionary algorithms. It was exposed by Russel Eberhart (electrical engineer) and James Kennedy (sociopsychologist) in 1995 [39]. It was originally inspired by the social behavior of animals evolving in swarms, such as schools of fish and flight groups of birds. We can see highly complex movement dynamics in these organisms, despite the fact that each individual has minimal intelligence and only local knowledge of his place in the swarm.

Optimization Algorithms
This paper focuses on two swarm-intelligence based algorithms. Specifically, they are the Particle Swarm Optimization and the Grey-Wolf Optimization algorithms. They have been selected for this PV system due to its demonstrated capability to adjust power converters.

A. The PSO algorithm
This algorithm was one of the solutions that was used for resolving a huge mathematical optimization problem. It is classified as an evolutionary metaheuristic that belongs to the larger class of evolutionary algorithms. It was exposed by Russel Eberhart (electrical engineer) and James Kennedy (sociopsychologist) in 1995 [39]. It was originally inspired by the social behavior of animals evolving in swarms, such as schools of fish and flight groups of birds. We can see highly complex movement dynamics in these organisms, despite the fact that each individual has minimal intelligence and only local knowledge of his place in the swarm. [40]. Therefore, local information and the memory of each individual are used to decide their displacement. Simple rules, such as "stay close to another candidate", "go in the same direction, "or" to go at the same speed, "are sufficient to maintain the cohesion of the swarm and allow the implementation of complex adaptative collective behaviors [26].
The particle swarm is a population of simple agents called particles. Each particle is considered as a solution to the problem, where it has a position and a speed. In addition, each particle has a memory allowing it to remember its best performance and the best performance achieved by "neighboring" particles (informants). In addition, each particle has a group of informants, historically called its neighborhood [41,42]. A swarm of particles, which are potential solutions to the optimization problem, will look for the global optimum with their movements.
The following three components that influence the movement of a particle are: 1. An inertia component: the particle tends to follow its current direction of movement; 2.
A cognitive component: the particle tends to move towards the best site through which it has already passed; 3. A social component: particle tends to rely on the congeners' experience and, thus, to go to the best site already reached by its neighbors.
The strategy for moving a particle is shown in Figure 5, where the three previous trends are illustrated.

2.
A cognitive component: the particle tends to move towards the best site through which it has already passed; 3. A social component: particle tends to rely on the congeners' experience and, thus, to go to the best site already reached by its neighbors.
The strategy for moving a particle is shown in Figure 5, where the three previous trends are illustrated. At the start of the algorithm, the swarm particles are randomly initialized in the search space of the problem. Then, at each iteration, each particle moves, linearly combining the three components mentioned above. Indeed, at iteration t + 1, the speed vector and the position vector are calculated from Equations (9) and (10), respectively, as follows: , = , + , ∈ 1,2. . . , where ω is a constant, called the coefficient of inertia; c1 and c2 are two constants, called acceleration coefficients; r1 and r2 are two random numbers drawn uniformly in the interval [0 1], at each iteration t and for each dimension j. The three components mentioned above (i.e., inertia, cognitive and social) are represented in Equation (9) by the following terms: In an N-dimensional search space, particle i of the swarm is modeled by its position vector → x i = (x i1 , x i2 , . . . , x iN ) and by its speed vector The value of the objective function determines the quality of its position at this point. The particle keeps in memory the best position through which it has ever passed, which we denote by → P best_i (p best_i1 , p best_i2 , . . . , p best_iN ). The best position reached by the particles of the swarm is denoted by → G best = (g best_1 , g best_2 , . . . , g best_N ).
At the start of the algorithm, the swarm particles are randomly initialized in the search space of the problem. Then, at each iteration, each particle moves, linearly combining the three components mentioned above. Indeed, at iteration t + 1, the speed vector and the position vector are calculated from Equations (9) and (10), respectively, as follows: where ω is a constant, called the coefficient of inertia; c1 and c2 are two constants, called acceleration coefficients; r1 and r2 are two random numbers drawn uniformly in the interval [0 1], at each iteration t and for each dimension j.
The three components mentioned above (i.e., inertia, cognitive and social) are represented in Equation (9) by the following terms: ωv t i,j corresponds to the inertia component of the displacement, where the parameter ω controls the influence of the direction of displacement on the future movement. c 1 r t 1,t p best t i,j − x t i,j corresponds to the cognitive component of the displacement, where the parameter c1 controls the cognitive behavior of the particle. c 2 r t 2i,j g best t j − x t i,j corresponds to the social component of the displacement, where the parameter c2 controls the social aptitude of the particle.
Once the particles have moved, the new positions are evaluated. The two vectors P best_i and G best are updated, at iteration t + 1, according to the two Equations (11) and (12). This procedure is presented in Algorithm 1, where M is the number of particles in the swarm.
The corresponding flowchart that describes this organization can be found in Figure 6 and explained as follows: The hunting process of the wolf pack involves three main steps: chasing, surroundi the prey, and attacking the prey. The algorithm starts with a given number of grey wolv with their positions generated arbitrarily.
The The prey position Xp (t + 1) update is calculated by averaging the positions of gr wolves α, β and ∆ (three temporarily optimal solutions). The following average functi is used for this purpose (1

B. The GWO Algorithm
Similarly to the particle swarm Optimization algorithm, one of the new algorithms, which are based on the metaheuristic principle is the Grey Wolf Optimization GWO algorithm. The researcher, Mirjalili, was one of the first researchers who developed this algorithm and exposed its running principle in 2014 [32,43]. To obtain the optimum solution of the problem to be optimized, the algorithm principle uses social authority, which is represented by the behaviour of the wolves when surrounding a victim. During the operation of hunting for the victim, this algorithm simulates the hierarchical supremacy of grey wolves until their movements end. It works in a similar way to population-based algorithms in which it simulates the natural behavior of grey wolves foraging for food in their social lives. Four types of grey wolf groups can be used to compose hierarchical commands. Figure 6 shows this hierarchy, with the following three levels: a.
The first level represents the group's command. A wolf in this level is called (alpha) (α). The alpha is responsible for deciding to hunt and orders the other wolves in the pack. Therefore, it can be considered as the correct solution. b.
The next level in the chain is called (β). Wolves in this level help the alpha to make decisions and monitor the actions of other groups. They can replace the alphas when they die or get older. c.
The lower level contains the delta and omega, which are the lowest ranks and who eat last after the wolves of the upper levels have finished.
The hunting process of the wolf pack involves three main steps: chasing, surrounding the prey, and attacking the prey. The algorithm starts with a given number of grey wolves with their positions generated arbitrarily.
The following equations determine the encircling behavior of each group of the pack: where → X(t) is the vector position of the grey wolf, → X p (t) is the vector prey position and → A and → C are the vectors gives by the following Equations: being t the current iteration, T max the total number of iterations and r 1 and r 2 are random vectors chosen in the interval [0,1]. The prey position Xp (t + 1) update is calculated by averaging the positions of grey wolves α, β and ∆ (three temporarily optimal solutions). The following average function is used for this purpose where: Equation (13) represents the distance from the current position, which should be minimized as much as possible so that the next position represented by Equation (14) gets closer and closer to the position of the prey. This will imply that the algorithm will get to the correct solution of the problem X P (t).
The parameter "a" used in this algorithm decreases linearly in the interval [2,0] for the successive iterations using Equation (15). Thus, it will model the behavior of the wolves when approaching the victim (exploration phase). For this phase, if the condition |A|< 1 is verified, the wolves attack the victim.
The alpha group are said to have the best possible knowledge of the location of prey. Once the position of the prey is determined, the hunt will be guided by the alpha group followed by the beta and delta wolves. The latter two groups participate in the hunt occasionally. The rest of the group is limited to take care of the injured wolves of the pack. When the prey stops moving, the wolves attack and finish the hunt [43].
The flowchart of this algorithm is illustrated in Figure 7. The choice of these parameters was fixed after several online simulation tests, where the goal is to find the best combination in terms of the algorithm running speed and the best performances. Specifically, we have varied the number of particles/wolves and the maximum number of iterations. The first test had the configuration of 150 iterations and 10 particles. The corresponding simulation time was evaluated to 30 min when using an I5 laptop with 8-GB as RAM memory. The resulting performances were found perfect. The second configuration was fixed to 50 iteration and 8 particles, and then the corresponding simulation time was evaluated to 24 min, but the resulting energetic performances were not so good. There were some problems with the stability of the output power.
Many other tests (more than 5 combinations) were also applied, and the best combination was found as it is indicated in Table 3. For the selected configuration, the simulation time was 26 min, and we get a good performance in terms of extracted power and stability.

C. The PSO and GWO MPPT controllers
To apply the previous evolutionary algorithms, we need to define the functions and parameters for the proposed problem. The objective is to maximize the energy extraction from the PV cells installed on EVs, even when there are shadowing areas involved in the scenario to consider.
For the PSO algorithm, the function to optimize is expressed in Equation (17).
In the GWO, each wolf position corresponds to the duty cycle factor applied on the MPPT. Thus, the Gmax represents the global best position of all of the wolves and Pmax is the corresponding best position of the corresponding wolf. If no partial shading exists, the global best wolf best position and the local wolf best position will be the same. However, for the partial shading case, there is several (Pmax) wolf best positions, and the global best position (Gmax) will be related to the max of the best wolf position Gmax = max (Pmax).
The proposed GWO algorithm aims to compensate for the limitations of conventional maximization algorithms, which cannot find the overall maximum, especially during a partial shading phenomenon. Regardless of the partial shading profile, this algorithm can find the maximum power point on the power-voltage characteristic of a photovoltaic system. Towards this goal, the objective function for this algorithm is defined as it is in Equation (18).
The proposed MPPT algorithm is based on the application of GWO to control the duty cycle. The expected consequence is that this will turn into decreasing the steady-state oscillations presented by the conventional MPP tracking algorithms. Consequently, the power loss due to oscillation is reduced, which increases the photovoltaic system efficiency.
For the implementation of GWO MPPT algorithm, duty cycle d is defined as a grey wolf. Therefore, the Equation (14) can be rewritten as follows: The fitness function of the GWO algorithm is formulated to have the objective of function 20. It is for the two optimization algorithms.
where P represents the PV power, d is the duty cycle of the boost converter, j is the current grey wolves number, and k is the iteration number. The flowchart and the parameters of the PSO and GWO MPPT algorithms are detailed in Figure 7 and Table 3, respectively. The choice of these parameters was fixed after several online simulation tests, where the goal is to find the best combination in terms of the algorithm running speed and the best performances. Specifically, we have varied the number of particles/wolves and the maximum number of iterations. The first test had the configuration of 150 iterations and 10 particles. The corresponding simulation time was evaluated to 30 min when using an I5 laptop with 8-GB as RAM memory. The resulting performances were found perfect. The second configuration was fixed to 50 iteration and 8 particles, and then the corresponding simulation time was evaluated to 24 min, but the resulting energetic performances were not so good. There were some problems with the stability of the output power.
Many other tests (more than 5 combinations) were also applied, and the best combination was found as it is indicated in Table 3. For the selected configuration, the simulation time was 26 min, and we get a good performance in terms of extracted power and stability.

Implementation and Simulation Results
The algorithms previously studied are applied to a PV system composed of 4 photovoltaic modules. They are connected in a combination of serial and parallel cells and the characteristics are given in Table 2. Considering the size of the PV modules, this system could be implemented in electric vehicle application or even an isolated grid for feeding some isolating area using a solar energy. Some recharge stations can also benefit from this solution and the proposed method can improve energy yield.
The simulation bloc can be illustrated as it is Figure 8, all the PV cells are connected to the DC bus, which can give information about the global voltage and current, in relation to the solar radiation factor. Then the MPPT algorithm can extract the maximum of power using the PSO or the GWO algorithm. A DC converter is used for having information about the outputted power delivered to the load.

Implementation and Simulation Results
The algorithms previously studied are applied to a PV system composed of 4 photovoltaic modules. They are connected in a combination of serial and parallel cells and the characteristics are given in Table 2. Considering the size of the PV modules, this system could be implemented in electric vehicle application or even an isolated grid for feeding some isolating area using a solar energy. Some recharge stations can also benefit from this solution and the proposed method can improve energy yield.
The simulation bloc can be illustrated as it is Figure 8, all the PV cells are connected to the DC bus, which can give information about the global voltage and current, in relation to the solar radiation factor. Then the MPPT algorithm can extract the maximum of power using the PSO or the GWO algorithm. A DC converter is used for having information about the outputted power delivered to the load. Indeed, to compare the efficiency of MPPT PSO and GWO we considered four simulation tests so that a wide variety of conditions are modelled, Table 4, gives the maximum power, voltage and current parameters for each shading condition. During the first test, Figure 9a, the irradiation was kept constant and uniform for the four modules (1000 W/m2). The second test is characterized by partial shading on two of the four modules of around 40%, represented in Figure 9b. It corresponds to the following lighting distribution on the four modules (600 W/m 2 , 600 W/m 2 , 1000 W/m 2 , 1000 W/m 2 ). The third test (in Figure 9c) resembles the second one, but the partial shading on the two modules is varied in the order of 50% on the first and 20% for the second (500 W/m 2 , 800 W/m 2 , 1000 W/m 2 , 1000 W/m 2 ). Three modules operate under partial shade during the fourth Indeed, to compare the efficiency of MPPT PSO and GWO we considered four simulation tests so that a wide variety of conditions are modelled, Table 4, gives the maximum power, voltage and current parameters for each shading condition. During the first test, Figure 9a, the irradiation was kept constant and uniform for the four modules (1000 W/m 2 ). The second test is characterized by partial shading on two of the four modules of around 40%, represented in Figure 9b. It corresponds to the following lighting distribution on the four modules (600 W/m 2 , 600 W/m 2 , 1000 W/m 2 , 1000 W/m 2 ). The third test (in Figure 9c) resembles the second one, but the partial shading on the two modules is varied in the order of 50% on the first and 20% for the second (500 W/m 2 , 800 W/m 2 , 1000 W/m 2 , 1000 W/m 2 ). Three modules operate under partial shade during the fourth Test, in  Test, in Figure 9d. This corresponds to the following lighting distribution (200 W/m 2 , 300 W/m 2 , 700 W/m 2 , 1000 W/m 2 )        C. Second case of partial shading.  C. Second case of partial shading.   Table 5 presents a comparative study between the two MPPT algorithms for di types of shading distribution on the panels. The essential criteria for a jud comparison between the two algorithms are: -Track of the point of maximum power under the different levels of partial sha -Rapid convergence towards the point of the global maximum. As can be observed, the simulation results from Figure 10 to Figure 13 sho evolution of the power, current, voltage and duty cycle of the PV system for the four of shading distribution we have tested on the panels. From these results, it c observed that the two solutions ensure a good MPP tracking. The advantage of th MPPT over the GWO MPPT is related to two issues: (i) the amplitude of the oscill at the transient state and (ii) the accuracy to track the point of maximum power. A oscillation exists for the case of GWO, which can be one the weaknesses of this algo There is also a small oscillation when executing the PSO at the beginning, wh radiation form changes. However, this will not cause a problem as in the real situ the modification of the radiation comes very slow. So, we a look a better performa a real situation.   Table 5 presents a comparative study between the two MPPT algorithms for different types of shading distribution on the panels. The essential criteria for a judicious comparison between the two algorithms are: As can be observed, the simulation results from Figure 10 to Figure 13 show the evolution of the power, current, voltage and duty cycle of the PV system for the four types of shading distribution we have tested on the panels. From these results, it can be observed that the two solutions ensure a good MPP tracking. The advantage of the PSO MPPT over the GWO MPPT is related to two issues: (i) the amplitude of the oscillations at the transient state and (ii) the accuracy to track the point of maximum power. A high oscillation exists for the case of GWO, which can be one the weaknesses of this algorithm. There is also a small oscillation when executing the PSO at the beginning, when the radiation form changes. However, this will not cause a problem as in the real situation, the modification of the radiation comes very slow. So, we a look a better performance in a real situation. The application of these algorithms in real time requires the use of the high-speed processor given the large number of operations to be carried out in one second (processing and control measurement). Therefore, the time needed to converge towards the best response depends on the speed of the algorithm used and the material available (essentially the speed of the processor). In addition, the presence of high amplitude oscillations during the transient phase is a harmful phenomenon for electrical systems and can cause a variety of problems. According to the simulation results and the criteria indicated above, the PSO MPPT algorithm shows itself well for the real-time application.
We have also studied how these algorithms impact on the battery State-of-Charge (SoC). In particular, we have studied the effects on a representative lithium-ion battery characterized by the parameters in Table 6. Figure 14 shows the efficiency of PSO MPPT versus GWO MPPT in terms of storage charge in the battery under uniform irradiation. A gain of about 0.0014% SoC for 1 s corresponds to almost 5.04% of battery charge for one hour. We have also studied how these algorithms impact on the battery State-of-Charge (SoC). In particular, we have studied the effects on a representative lithium-ion battery characterized by the parameters in Table 6. Figure 14 shows the efficiency of PSO MPPT versus GWO MPPT in terms of storage charge in the battery under uniform irradiation. A gain of about 0.0014% SoC for 1 s corresponds to almost 5.04% of battery charge for one hour.

Conclusions
This work attempts to study and discuss two MPPT techniques based on two metaheuristic optimization algorithms, i.e., PSO and GWO. These new techniques of MPPT overcome the problems of classic MPPT strategies (e.g., perturb and observe and incremental) when tracking the maximum power point, even in the presence of sudden changes of irradiation and shadows on the photovoltaic modules. The MPPT techniques studied show good behavior and better performance. A comparative study of simulation results for a different type of shading shows PSO-MPPT effectiveness compared to GWO-MPPT from the point of view of speed and oscillation during the transient state. In addition, a simulation test shows the efficiency of PSO MPPT versus GWO MPPT in terms of storage charge in the battery under uniform irradiation. As a future work, we would like to analyze how to adapt the algorithm's parameters (which are now constant) to the dynamic lighting conditions. Moreover, one of the future endeavors of this work is to compare more algorithms performances and search the best combination that can be used for such an optimization problem. Therefore, incremental algorithm, perturb and observe algorithm, the fuzzy solution and other population-based metaheuristic algorithms as bio-inspired algorithms, evolutionary algorithm and physics-based algorithm will be studied, tested and evaluated.

Conclusions
This work attempts to study and discuss two MPPT techniques based on two metaheuristic optimization algorithms, i.e., PSO and GWO. These new techniques of MPPT overcome the problems of classic MPPT strategies (e.g., perturb and observe and incremental) when tracking the maximum power point, even in the presence of sudden changes of irradiation and shadows on the photovoltaic modules. The MPPT techniques studied show good behavior and better performance. A comparative study of simulation results for a different type of shading shows PSO-MPPT effectiveness compared to GWO-MPPT from the point of view of speed and oscillation during the transient state. In addition, a simulation test shows the efficiency of PSO MPPT versus GWO MPPT in terms of storage charge in the battery under uniform irradiation. As a future work, we would like to analyze how to adapt the algorithm's parameters (which are now constant) to the dynamic lighting conditions. Moreover, one of the future endeavors of this work is to compare more algorithms performances and search the best combination that can be used for such an optimization problem. Therefore, incremental algorithm, perturb and observe algorithm, the fuzzy solution and other population-based metaheuristic algorithms as bio-inspired algorithms, evolutionary algorithm and physics-based algorithm will be studied, tested and evaluated.