Abstract

In recent years, the solar energy has become one of the most important alternative sources of electric energy, so it is important to operate photovoltaic (PV) panel at the optimal point to obtain the possible maximum efficiency. This paper presents a new optimization approach to maximize the electrical power of a PV panel. The technique which is based on objective function represents the output power of the PV panel and constraints, equality and inequality. First the dummy variables that have effect on the output power are classified into two categories: dependent and independent. The proposed approach is a multistage one as the genetic algorithm, GA, is used to obtain the best initial population at optimal solution and this initial population is fed to Lagrange multiplier algorithm (LM), then a comparison between the two algorithms, GA and LM, is performed. The proposed technique is applied to solar radiation measured at Helwan city at latitude 29.87°, Egypt. The results showed that the proposed technique is applicable.

1. Introduction

In the last years, global warming and energy policies have become a hot topic on the international agenda as the developed countries are trying to reduce their greenhouse gas emissions. Renewable energy sources are considered as a technological option for generating clean energy. Among them, photovoltaic (PV) system has received a great attention as it appears to be one of the most promising renewable energy sources. PV power generation has an important role to play as it is a green source. The only emissions associated with PV power generation are those from the production of its components. However, the development for improving the efficiency of the PV system is still a challenging field of research. Maximum power point tracking (MPPT) algorithms are necessary in PV applications because the MPP of a solar module varies with the irradiation and temperature as shown in Figures 1 and 2. So, the use of MPPT is required to obtain the maximum output power from a module [1]. The optimization process can be defined as the finding values of variables that minimize or maximize the objective function and satisfy the constraints. The optimization problems are centered on three main factors which are as follows:(1)objective function which is to be minimized or maximized;(2)a set of unknown variables that have effect on the objective function;(3)A set of constraints that allow the unknown to take on certain values but exclude others.

The most popular literature optimization techniques of the PV power are based on the usage of the MPPT. The MPPT control method uses one estimate processes between every two Perturb processes in search for the maximum PV output (EPP) which is proposed in [2].

An intelligent approach for MPPT DC-DC Boost converter focused on Perturb and Observe (P&O) algorithm and compared to a designed fuzzy logic controller which is presented in [3]. A comparative study of two types of MPPT which are P&O and incremental conductance method is introduced in [4]. The implementation of fuzzy logic controller based on the change of the PV module power, , and its change with respect to the module voltage, is studied in [5], the proposed fuzzy logic determines the size of the voltage then obtains the corresponding Power. The performance of fuzzy logic with various membership functions (MFs) is tested to optimize the MPPT. A novel intelligent fuzzy logic controller for MPPT based on boost converter and single-phase grid-connected inverter is introduced in [6]. An intelligent control method for MPPT of a photovoltaic system under variable temperature and insolation conditions which uses a fuzzy logic controller applied to a DC-DC converter device is proposed in [7]. A fuzzy logic control to control MPPT for a PV system is proposed in [8]. The analysis of the optimal operation of PV panels as a function of the weather conditions (solar irradiation, temperatures), and the design of a PV system provided with MPPT command ensuring instantaneously optimal operation of photovoltaic panels is presented in [9]. MPPT simulation based on using particle swarm optimization (PSO) algorithm is presented in [10]. The algorithm is employed on a buck-boost converter and tested experimentally using a PV array simulator. The design and implementation of an optimized stand-alone solar pumping system are presented in [11].

This paper presents a new optimization approach to maximize the electrical power of the PV panel. The technique which is based on objective function represents the output power of the PV module and constraints, equality and inequality. First, the dummy variables that have effect on the output power are classified into two categories: dependent and independent. The proposed approach is a multistage one as GA is used to obtain the best initial population at optimal solution, and then the best initial population is fed to Lagrange multiplier algorithm (LM), then a comparison between two algorithms, GA and LM, is performed. The proposed technique is applied to solar radiation measured at Helwan city at latitude 29.87°, Egypt. The results showed that the proposed technique is applicable.

2. Mathematical Model of PV Module

The solar flux striking a PV panel can be resolved to three components; direct beam radiation that passes in a straight line through the atmosphere to the receiver, diffuse radiation that has been scattered by molecules in the atmosphere, and reflected radiation that has bounced off the ground or other surface in front of the collector as shown in Figure 3.

The total rate of radiation strikes a collector on a clear day as follows [12]: where is the air mass, is the altitude angle, is the solar azimuth angle, is the PV module azimuth angle, positive for eastsouth and negative for westsouth, is the PV module tilt angle, is reflection factor, and is sky diffuse factor and is given by the following equation [12]: where and are dependent on the day number and can be obtained by the following equations [1] and is the day number:

A PV cell can be simulated by a real diode in parallel with an ideal current source, , which depends on striking radiation. The generalized equivalent circuit of the PV cell is given in Figure 4 [13, 14].

One can derive the following equations for current and voltage:

The reverse saturation current, , is dependent on the temperature and is given by the following equation:

The short circuit current, , depends on the solar radiation and cell temperature as follows: where is the cell output current, is the cell output voltage, and is the cell series resistance ().

is the cell parallel resistance (), is the diode ideality factor, is the Boltzmann constant ( J/K), is the cell temperature (°C), is the reference temperature of the PV cell, is the cell reverse saturation current at reference temperature, is the band gap of semiconductor used in the cell, is the cell short circuit current at reference temperature and radiation, is the short circuit current temperature coefficient, and is the solar radiation strikes a tilted module in W/m2.

The cell temperature is calculated as follows: where is ambient temperature and is cell temperature in a module when ambient is 20°C.

A PV module consists of series cells, , and parallel branches, , which are given in Figure 5.

The PV module’s current under variable operating conditions can be described as follows [15]:

3. Proposed Technique

3.1. The Proposed Optimization Problem

The mathematical model of any continuous optimization problem can be described as follows: where is an objective function to be minimized or maximized subjected to equality and inequality constraints as follows:(1)equality constraint(2)parametric inequality constraint(3)functional inequality constraintwhere is a vector of state, dependant, variables and is a vector of control (independent) variables. First the dummy variables that have effect on the PV module power are classified into two categories: independent () and dependent variables ()— and .

The proposed objective function is expressed in the following form:

Subjected to the following proposed parametric constrained: the proposed equality constraint is given as

The limits of independent variables are selected according to the following aspects.(1)When , the PV module becomes horizontal and produces power, while for , the PV module becomes vertical and produces zero power; so, the selected limits are assumed between 0° and 80°.(2)The PV module azimuth angle is positive for east of south and becomes negative for west of south; so, the limits are selected as .

3.2. Genetic Algorithm

Genetic algorithm, GA, is gradient free, parallel optimization algorithms that use a performance criterion for evaluation and a population of possible solutions to the search for a global optimum. GA is capable of handling complex and irregular solution spaces and has been applied to various difficult optimization problems [16]. It is inspired by the biological process of Darwinian evolution where selection, mutation, and crossover play a major role. The manipulation is done by the genetic operators that work on the chromosomes in which the parameters of possible solutions are encoded. In each generation of the GA, the new solutions replace the solutions in the population that is selected for removal. The main elements of GA are populations of chromosomes, selection according to fitness, crossover to produce new offspring, and random mutation of new offspring. A simple GA flow chart is shown in Figure 6.

The simplest form of genetic algorithm involves three types of operators: selection, crossover (single point), and mutation.

Selection

It selects chromosomes in the population for reproduction.

Crossover

It randomly chooses a locus and exchanges the subsequences before and after that locus between two chromosomes to create two offsprings.

Mutation

It randomly flips some of the bits in a chromosome.

3.3. The Method of Lagrange Multipliers

Lagrange multipliers, LM, play a vital role in the study of constrained optimization. It can be interpreted as the rate of change in the objective function with respect to changes in the associated constraint function [16]. The main formula of this method can be described as follows: where and are the respective gradients. Constant in (16) are called the Lagrange multipliers of the constrained problem.

By solving (15) one can get the value of and then obtain the value of vector and vector at optimal solution and the optimal power extracted from the PV module. The proposed analysis uses a real data for solar radiation and ambient temperature measured by solar radiation and meteorological station located at National Research Institute of Astronomy and Geophysics Helwan, Cairo, Egypt, which is located at latitude 29.87°N and longitude 31.30°E. The station is over a hill top of about 114 m height above sea level.

The initial condition of any optimization problem plays an important role in finding the optimal solution so GA program is performed under varying the initial condition of the independent variables until the optimal, maximum, power is obtained then Lagrange multiplier algorithm program is used and compared with GA to investigate the validity and quality of our proposed method.

The main steps of the proposed technique can be summarized as follows.(1)Using the real data for solar radiation and ambient temperature as an input to genetic algorithm program that obtain the optimal output power from the PV module.(2)Studying the effect of changing the initial population on the value of fitness function in order to obtain the optimal (maximum) possible power from the PV module.(3)The initial population at the optimal PV output power is fed to the Lagrange multiplier algorithm and comparing the results with those obtained by GA. Figure 7 shows these proposed technique steps.

4. Numerical Analysis

The analysis uses bpsx150 PV module which has 72 cells connected in series and its electrical characteristic is shown in Table 1 [17].

The bpsx150 PV module is simulated by MATLAB/Simulink toolbox as shown in Figure 8. Referring to Figure 8, there are three main inputs which are the module voltage, the solar radiation, and the ambient temperature; these items are labeled in Figure 9 with numbers 1, 2, and 3. Simulink is a software package of MATLAB program for modeling, simulating, and analyzing dynamical systems. It supports linear and nonlinear systems, modeled in continuous time, sampled time, or a hybrid of the two. For modeling, Simulink provides a graphical user interface, GUI, for building models as block diagrams, using click and drag mouse operations. The detailed configuration of the PV module is given in Figure 9. The analysis of the proposed algorithm is performed on measured solar radiation and ambient temperature as inputs to GA program. These data are measured at Helwan city, for the sunny day of June 10, 2012, and start from 6:10 AM to 5:50 PM. The direct beam solar radiation and the corresponding ambient temperature, Helwan city, at latitude 29.87°, Egypt, are given in Figure 10.

The initial population of GA program which contains plays an important role in finding the optimal power, so they are changed as follows.(1) The initial tilt angle, , is changed at acceptable limits, fixing both and at acceptable values until the maximum power is obtained as shown in Table 2.

The optimal output power for this case is equal to at which is fixed at this value.(2) Changing the initial PV module azimuth angle, , at acceptable limits with fixed acceptable initial module voltage until the optimal power is obtained as shown in Table 3.

The optimal output power for this case is equal to at which is fixed at this value.(3) Changing the initial PV module voltage until the optimal power is obtained as shown in Table 4.

The final optimal output power is equal to at .

Finally, one can derive that the optimal obtained power by GA at 6:10 AM is equal to and is occurred at initial population which is equal to . The curve shows the change of PV module power with changing the initial population is shown in Figure 11. GA program response is shown in Figure 12 that contains the fitness function, the current best individuals, the average distance between individuals and the number of individuals for each iteration. The constant dependant values at optimal power are , , , , and .

The final population at optimal solution is equal to , , and .

The optimal power and the corresponding tilt angle for 24 hours of the day are shown in Figure 13 and Table 5.

Once the bpsx150 PV module has a fixed number of series cells, , and a fixed number of parallel branches, , then they are not be considered as variables.

The initial condition at optimal power obtained from the GA is fed to the Lagrange multiplier algorithm to find the optimal power. The MATLAB function fmincon is used to perform the Lagrange multiplier analysis.

The Lagrange multiplier algorithm response is shown in Figure 14. A comparison between the GA and Lagrange multiplier algorithm is performed and is shown in Figure 15 and Table 6.

A comparison between the GA and Lagrange multiplier algorithm is performed and is shown in Figure 14 and Table 6.

It is clear that the difference between two proposed solution methods is very small as the initial condition of them is the same; the operating conditions of the PV module do not change. So, one can derive that the error between two methods is acceptable and the proposed method is applicable.

The benefits of this proposed approach are as follows.(i)A new objective function of the PV module output power subjected to constraints is presented.(ii)A simple optimization methodology to obtain the optimal power extracted from the PV module and corresponding tilt angle based on real data is formulated.(iii)Two suggested solution techniques are manipulated and the error between them is very small.

5. Conclusion

In this work a new optimization approach to maximize the PV module output power is presented. A new proposed objective function of PV module power and constraints is also presented. Genetic algorithm is performed under varying the initial population of independent variables, then Lagrange multiplier algorithm is simulated. The optimal PV module power obtained from both algorithms is compared. The analysis is based on real measured data of solar radiation and ambient temperature at Helwan city at latitude 29.87°, Egypt. The analysis is done in complete sunny day and the results showed that the error between two proposed solution techniques is acceptable and the proposed constrained optimization method is applicable.

Acknowledgment

Cordial thanks and deep respect are offered to the Solar and Space Research Department, National Research Institute of Astronomy and Geophysics, Helwan, Egypt, for supporting this work and providing the authors with the required real data.

Supplementary Materials

This paper presented a new optimization technique to maximize the PV module output power based on two different solution methods; GA and Lagrange Multiplier. The analysis is performed for actual measured data for both solar radiation and ambient temperature.

  1. Supplementary Material