Performance Enhancement of Photovoltaic Water Pumping System Based on BLDC Motor under Partial Shading Condition †

: The use of photovoltaic (PV) energy for water pumping is considered one of the most promising areas of photovoltaic applications. This work aims to improve the power extraction of a battery-less photovoltaic pumping system. The DC–DC converter is used to accomplish the maximum power point tracking (MPPT). In addition, the BLDC motor is used for maximum exploitation of the delivered power by the PV array, and to enhance the reliability of the pumping system. Furthermore, an MPPT strategy, based on the cuckoo swarm optimization technique, has been developed to improve the control performance under partial shading conditions. The effectiveness of the proposed system has been examined through simulation using MATLAB and Simulink software.


Introduction
Water resources are essential for satisfying human needs, ensuring food production and protecting health, as well as for social and economic development. Worldwide, water pumping is generally dependent on electricity or diesel-generated electricity. Diesel is expensive and scarce in the countryside of many developing countries, and even when it is available, transporting it to remote areas is difficult. This is because there are no roads or supporting infrastructure in most of the remote villages. In rural areas, the installation of a new transmission line and transformers is extremely expensive for irrigation or drinking water [1]. A solar water pumping system (SWPS) is made up of various components that can be classified as mechanical, electrical, or electronic. These components have different construction, operation, and performance characteristics [2].
PV pumping systems based on AC motors, especially induction motors (IM), are often favored because they are more dependable, cost less, and do not need regular maintenance. DC motors are commonly used in low-power solar PV water pumps [3]. However, DC motors with brushes possess low efficiency, and they need frequent maintenance, owing to the sliding brush contacts and the commutator. Three-phase brushless direct-current (BLDC) motors are suitable options, due to their high-efficiency capabilities, compactness, and easy-to-drive qualities [4]. Therefore, this motor has gained considerable attention for water pumping in the recent decade, due to its many advantages, which constitute its desirable characteristics, especially for this application [3].
The use of maximum power point tracking (MPPT) methods is critical to increase the efficiency of photovoltaic (PV) systems. Because of its good balance between complexity, precision, and dependability, the perturb and observe (P&O) approach is the most prevalent algorithm [5]. It does, however, have several disadvantages, such as a delayed response time, steady-state power oscillation around the MPP, and poor tracking under abrupt irradiance changes or shading effects [6]. Particle swarm optimization (PSO) and cuckoo search (CS) are two optimization strategies that have been proposed to increase the performance of MPPT [7]. Cuckoo search (CS) is an evolutionary optimization technique that was inspired by a type of bird, named the cuckoo [8], by Yang and Deb (2009). Cuckoos put their eggs in other birds' or species' nests [9]. They fly from one nest to the next, looking for the best one, the one that offers the greatest possibility of the eggs hatching successfully. Several papers [10,11] have discussed the applicability of CS to MPPT. It is shown that it is more stable, has greater convergence, and is more efficient.
In this paper, we present a method for tracking the maximum power point of solar water pumping systems under partial shading conditions. The proposed system is based on a BLDC motor to improve the reliability and efficiency of the system. Moreover, a cuckoo search algorithm is proposed to overcome the disadvantages of the conventional MPPT algorithm. The proposed system has been investigated using MATLAB/Simulink software. Figure 1 depicts the designed solar water pumping system's configuration. The photovoltaic generator (PVG) is a group of solar panels with connections that are selected to meet the power needs of the association's power electronic converters, motor, and centrifugal pump. The DC-DC converter is connected between a PVG and a pump to extract the maximum energy possible throughout the day. The peak power operating point of the PVG is tracked using an MPPT algorithm, under various conditions [12].

Configuration of the PV Pumping System
The use of maximum power point tracking (MPPT) methods is critical to increase the efficiency of photovoltaic (PV) systems. Because of its good balance between complexity, precision, and dependability, the perturb and observe (P&O) approach is the most prevalent algorithm [5]. It does, however, have several disadvantages, such as a delayed response time, steady-state power oscillation around the MPP, and poor tracking under abrupt irradiance changes or shading effects [6]. Particle swarm optimization (PSO) and cuckoo search (CS) are two optimization strategies that have been proposed to increase the performance of MPPT [7]. Cuckoo search (CS) is an evolutionary optimization technique that was inspired by a type of bird, named the cuckoo [8], by Yang and Deb (2009). Cuckoos put their eggs in other birdsʹ or speciesʹ nests [9]. They fly from one nest to the next, looking for the best one, the one that offers the greatest possibility of the eggs hatching successfully. Several papers [10,11] have discussed the applicability of CS to MPPT. It is shown that it is more stable, has greater convergence, and is more efficient.
In this paper, we present a method for tracking the maximum power point of solar water pumping systems under partial shading conditions. The proposed system is based on a BLDC motor to improve the reliability and efficiency of the system. Moreover, a cuckoo search algorithm is proposed to overcome the disadvantages of the conventional MPPT algorithm. The proposed system has been investigated using MATLAB/Simulink software. Figure 1 depicts the designed solar water pumping systemʹs configuration. The photovoltaic generator (PVG) is a group of solar panels with connections that are selected to meet the power needs of the associationʹs power electronic converters, motor, and centrifugal pump. The DC-DC converter is connected between a PVG and a pump to extract the maximum energy possible throughout the day. The peak power operating point of the PVG is tracked using an MPPT algorithm, under various conditions [12]. To control the BLDC motor, a voltage source inverter must be included, which provides three-phase voltages. The boost duty cycle is obtained by the MPPT algorithm, which mainly depends on irradiance.

Photovoltaic Module
For simplicity and accuracy, the single-diode model is the most studied and used in this work. The single-diode model is shown in Figure 2, and includes a photo current source, a parallel diode, and series and shunt resistors, Rs and Rsh, respectively [13]. To control the BLDC motor, a voltage source inverter must be included, which provides three-phase voltages. The boost duty cycle is obtained by the MPPT algorithm, which mainly depends on irradiance.

Photovoltaic Module
For simplicity and accuracy, the single-diode model is the most studied and used in this work. The single-diode model is shown in Figure 2, and includes a photo current source, a parallel diode, and series and shunt resistors, R s and R sh , respectively [13]. The model in Figure 2 can be mathematically described by the following equation: where Vt is the thermal voltage, given by the following: k is the Boltzmann constant (1.38 × 10 −23 J/K), T is the absolute temperature (k), and q is the electronic charge (1.6 × 10 −19 C).
The PV array output current and voltage can be computed according to the paralleland series-connected modules. The PV module considered in this work is the TSM-250PA05.08 PV module, Trinasolar, China.

DC-DC Boost Converter
The step-up chopper can be applied to MPPT systems in which the output voltage needs to be greater than the input voltage. The relation between the output and input voltages (Vout and Vin, respectively) is as follows: where D is the duty cycle; the duty cycle always ranges from zero to one.

Brushless DC Motor
The BLDC motorʹs mathematical model can be formed in the same way as a threephase synchronous machine. Similarly, with the BLDC motor, the armature winding model is stated as follows [14]:

Centrifugal Pump
A centrifugal pump is a mechanical device that uses centrifugal force to transfer mechanical energy into pressure in a fluid. Its efficiency is rather good, and it can pump a huge amount of water [1]. The pump's torque is proportional to the rotor's squared speed, as follows: where is the proportionality constant, and it is given by the following: The model in Figure 2 can be mathematically described by the following equation: where V t is the thermal voltage, given by the following: k is the Boltzmann constant (1.38 × 10 −23 J/K), T is the absolute temperature (k), and q is the electronic charge (1.6 × 10 −19 C).
The PV array output current and voltage can be computed according to the paralleland series-connected modules. The PV module considered in this work is the TSM-250PA05.08 PV module, Trinasolar, China.

DC-DC Boost Converter
The step-up chopper can be applied to MPPT systems in which the output voltage needs to be greater than the input voltage. The relation between the output and input voltages (V out and V in , respectively) is as follows: where D is the duty cycle; the duty cycle always ranges from zero to one.

Brushless DC Motor
The BLDC motor's mathematical model can be formed in the same way as a threephase synchronous machine. Similarly, with the BLDC motor, the armature winding model is stated as follows [14]:

Centrifugal Pump
A centrifugal pump is a mechanical device that uses centrifugal force to transfer mechanical energy into pressure in a fluid. Its efficiency is rather good, and it can pump a huge amount of water [1]. The pump's torque is proportional to the rotor's squared speed, as follows: where K P is the proportionality constant, and it is given by the following: The pump's head and available mechanical power at the rotating impeller determine the pump's water flow and pressure. Affinity laws, which only require the pump ratings and actual input parameters, rotor speed and torque [1], might simplify the estimation of the pump's output characteristics.
where H, Q and P are the rated parameters of the pump at the rated speed, and Ω rn , H , Q and P are the parameters of the pump at speed Ω r , which is different to the rated speed.

Maximum Power Point Tracking
The perturb and observe (P&O) method is a basic iterative approach. It perturbs the operational point of the system, causing the PV array terminal voltage to oscillate about the MPP voltage, even though the solar irradiation and the cell temperature are constant. Figure 3 depicts the P&O algorithm flowchart [15].
The pump's head and available mechanical power at the rotating impeller determine the pump's water flow and pressure. Affinity laws, which only require the pump ratings and actual input parameters, rotor speed and torque [1], might simplify the estimation of the pump's output characteristics.
where H, Q and P are the rated parameters of the pump at the rated speed, and Ω , H′, Q′ and P′ are the parameters of the pump at speed Ω , which is different to the rated speed.

Maximum Power Point Tracking
The perturb and observe (P&O) method is a basic iterative approach. It perturbs the operational point of the system, causing the PV array terminal voltage to oscillate about the MPP voltage, even though the solar irradiation and the cell temperature are constant. Figure 3 depicts the P&O algorithm flowchart [15].

Commutation and BLDC Drive
In the BLDC motor, Hall effect sensors are required. The signal is generated by the motor's Hall effect sensor, and must be decoded in order to determine the right phaseswitching sequence [4]. Figure 4 shows the control scheme of the BLDC motor; the reference current Iref is delivered by an outer PI speed controller. Then, the current is controlled using hysteresis controllers.

Commutation and BLDC Drive
In the BLDC motor, Hall effect sensors are required. The signal is generated by the motor's Hall effect sensor, and must be decoded in order to determine the right phaseswitching sequence [4]. Figure 4 shows the control scheme of the BLDC motor; the reference current I ref is delivered by an outer PI speed controller. Then, the current is controlled using hysteresis controllers.

Partial Shading Effect
Shading in a PV module reduces power, but also creates current mismatch inside a PV series, and string and voltage mismatch between parallel strings. In order to reduce the power damage to the panel, a bypass diode is connected in parallel with the PV

Partial Shading Effect
Shading in a PV module reduces power, but also creates current mismatch inside a PV series, and string and voltage mismatch between parallel strings. In order to reduce the power damage to the panel, a bypass diode is connected in parallel with the PV module to provide an alternate current channel. A module's cells, on the other hand, no longer carry the same current. As a result, there are numerous maxima on the power-voltage curve. Unfortunately, most standard MPPT algorithms cannot tell the difference between local and global maxima. To investigate the performance of the SWPS under partial shade conditions, the PV modules were adjusted to varied irradiance values (400, 800, and 1000). The P-V characteristic exhibits different maxima, due to the influence of bypass diodes ( Figure 5).

Partial Shading Effect
Shading in a PV module reduces power, but also creates current mismatch inside a PV series, and string and voltage mismatch between parallel strings. In order to reduce the power damage to the panel, a bypass diode is connected in parallel with the PV module to provide an alternate current channel. A module's cells, on the other hand, no longer carry the same current. As a result, there are numerous maxima on the powervoltage curve. Unfortunately, most standard MPPT algorithms cannot tell the difference between local and global maxima. To investigate the performance of the SWPS under partial shade conditions, the PV modules were adjusted to varied irradiance values (400, 800, and 1000). The P-V characteristic exhibits different maxima, due to the influence of bypass diodes ( Figure 5).

Cuckoo Search Methodology
Firstly, each cuckoo lays one egg at a time and dumps its egg into a randomly chosen nest. Then, the best nest, with high-quality eggs, will carry over to the next generation. The number of available hosts' nests is fixed, and the egg laid by a cuckoo is discovered by the host bird with a probability Pa, such that 0 ≤ Pa ≤ 1. In this case, the host bird can either throw the egg away or abandon the nest and build a completely new nest. For simplicity, this last assumption can be approximated by the fraction, Pa, of nests, n, that are replaced by new nests [9].

Lévy Flight
Lévy flight is a random walk, where step sizes are extracted from Lévy distribution according to the power law y = l λ [10], where l is the flight length and λ is the variance.

Cuckoo Search Methodology
Firstly, each cuckoo lays one egg at a time and dumps its egg into a randomly chosen nest. Then, the best nest, with high-quality eggs, will carry over to the next generation. The number of available hosts' nests is fixed, and the egg laid by a cuckoo is discovered by the host bird with a probability Pa, such that 0 ≤ Pa ≤ 1. In this case, the host bird can either throw the egg away or abandon the nest and build a completely new nest. For simplicity, this last assumption can be approximated by the fraction, Pa, of nests, n, that are replaced by new nests [9].

Lévy Flight
Lévy flight is a random walk, where step sizes are extracted from Lévy distribution according to the power law y = l −λ [10], where l is the flight length and λ is the variance. Since 1 < λ < 3, y has infinite variance. In CS, the nest searching steps of the cuckoo are characterized by Lévy flight. Figure 6 depicts an example of Lévy flight in a 2D plane. Since 1 < λ < 3, y has infinite variance. In CS, the nest searching steps of the cuckoo are characterized by Lévy flight. Figure 6 depicts an example of Lévy flight in a 2D plane. To generate a new solution, xt + 1, for a cuckoo, Lévy flight is performed, as dictated by the following expression [10]: where: To generate a new solution, xt + 1, for a cuckoo, Lévy flight is performed, as dictated by the following expression [10]: where: ⊕: Entry wise multiplication λ is the Levy exponent, x i is the current location (samples/eggs).
Based on the constraints imposed by the optimization problem, it is important to tune the value of α to obtain the desired step size. In most cases, α is used as follows: where α 0 is the initial step change.

MPPT for PV Using CS
To use CS for designing MPPT, appropriate variables have to be selected for the search. Firstly, the samples are used, which are the PV voltages in our case, i.e., V i (i = 1, 2, . . . , n). The total number of samples is defined as n. Then, the step size is denoted by the fitness function J, which is the value of PV power at the MPP. Since J is dependent on the PV voltage, J = f (V). Initially, the generated samples are applied to the PV modules and the power is set as the initial fitness value. The maximum power provided by its corresponding voltage is considered as the current best sample. Thereafter, Lévy flight is performed; consequently, new voltage samples are generated, as follows: where α = α 0 (V best − V i ). A simplified scheme of the Lévy distribution is presented as follows: where k is the Lévy multiplying coefficient. In the following equation, u and v are determined from the normal distribution curves: If Γ denotes the integral gamma function, then the variables σ u and σ v are defined as follows: The greatest power provided by the voltage is chosen as the new best sample by comparing the power values. Aside from this best sample, the rest are destroyed at random, with a probability of Pa. This procedure mimics the behavior of the host bird finding and then destroying the cuckoo's eggs. Then, to replace the destroyed samples, new random samples are created. As a result, all the sample powers are measured again, and the current best is chosen by assessing J. The process is repeated until all of the samples have arrived at the MPP.

Simulation Results
The proposed SWPS is simulated using MATLAB/Simulink 2018a. A PV array with four solar modules, connected 2 × 2 in series and in parallel, is considered. In normal conditions, the PV array would be expected to produce 2500 W, with an MPPT voltage and current of 310 V and 8.5 A, respectively, at 1000 W/m 2 solar irradiance and 25 • C.
The simulation results of the PV side in SWPS, under shaded conditions, using P&O and CS MPPT algorithms, give the resulting waveforms shown in Figures 7 and 8, where "a-b" indicate P&O and "c-d" indicate CS.
comparing the power values. Aside from this best sample, the rest are destroyed at random, with a probability of Pa. This procedure mimics the behavior of the host bird finding and then destroying the cuckoo's eggs. Then, to replace the destroyed samples, new random samples are created. As a result, all the sample powers are measured again, and the current best is chosen by assessing J. The process is repeated until all of the samples have arrived at the MPP.

Simulation Results
The proposed SWPS is simulated using MATLAB/Simulink 2018a. A PV array with four solar modules, connected 2 × 2 in series and in parallel, is considered. In normal conditions, the PV array would be expected to produce 2500 W, with an MPPT voltage and current of 310 V and 8.5 A, respectively, at 1000 W/m 2 solar irradiance and 25 °C.
The simulation results of the PV side in SWPS, under shaded conditions, using P&O and CS MPPT algorithms, give the resulting waveforms shown in Figures 7 and 8, where "a-b" indicate P&O and "c-d" indicate CS. As can be observed from the results, the P&O algorithm failed to track the global MPP (Figure 7a,b); it converged to the third peak in Figure 5. The output power from the PV array is reduced. Under partial shading conditions, the output power from the PV array is reduced. On the other hand, the CS method (Figure 7c,d) captured the global peak effectively, avoiding the local maximums, under partial shading conditions, which indicates maximum extraction of the PV power. Compared to P&O, which tracked a local peak, the CS exhibits a faster convergence speed, less oscillation around the MPP under steady-state conditions, and no divergence from the MPP during varying weather conditions. Figure 8 shows the dynamic performance of the BLDC motor under shaded conditions, using P&O and cuckoo search. The BLDC motor shows better performance (i.e., torque and speed) when using the CS algorithm, since it runs on the maximum possible power, as the CS tracks the MPP of the PV array. Table 1 presents the comparison between the two MPPT algorithms, in regards to the power delivered to the pump and the rotor speed. As can be observed from the results, the P&O algorithm failed to track the global MPP (Figure 7a,b); it converged to the third peak in Figure 5. The output power from the PV array is reduced. Under partial shading conditions, the output power from the PV array is reduced. On the other hand, the CS method (Figure 7c,d) captured the global peak effectively, avoiding the local maximums, under partial shading conditions, which indicates maximum extraction of the PV power. Compared to P&O, which tracked a local peak, the CS exhibits a faster convergence speed, less oscillation around the MPP under steady-state conditions, and no divergence from the MPP during varying weather conditions. Figure 8 shows the dynamic performance of the BLDC motor under shaded conditions, using P&O and cuckoo search. The BLDC motor shows better performance (i.e., torque and speed) when using the CS algorithm, since it runs on the maximum possible power, as the CS tracks the MPP of the PV array. Table 1 presents the comparison between the two MPPT algorithms, in regards to the power delivered to the pump and the rotor speed.

Conclusions
In this work, a BLDC motor-driven water pumping system, based on PV generation, was designed. The power generated from such systems does not only depend on the climatic conditions, but also highly depends on the applied control technique. MPPT control techniques were used in order to improve the total system efficiency and extract the maximum PV energy under different operating conditions.

Conclusions
In this work, a BLDC motor-driven water pumping system, based on PV generation, was designed. The power generated from such systems does not only depend on the climatic conditions, but also highly depends on the applied control technique. MPPT control techniques were used in order to improve the total system efficiency and extract the maximum PV energy under different operating conditions.
The proposed system was simulated and tested using P&O and cuckoo MPPT techniques in MATLAB/Simulink. The performance of the two algorithms was compared under partial shading conditions, where P&O failed to track the global MPP, while the cuckoo search did successfully. Moreover, the results show the improvement in performance of the BLDC motor under non-uniform irradiance conditions.
The parameters of the BLDC motor are as follows: Rated torque 4.2 N·m, power 1.5 kW, speed 2889 rpm, number of poles 10, phase resistance 0.7 Ω, phase inductance 4 mH, torque constant K t = 0.7 N·m/A, and moment of inertia 0.0028 kg/m 3 .
The parameters of the DC-DC converter are as follows: C = 2.34 × 10 −5 F and L = 0.023 H.