Hybrid Moth-Flame Optimization Algorithm and Incremental Conductance for Tracking Maximum Power of Solar PV / Thermoelectric System under Di ﬀ erent Conditions

: For an e ﬃ cient energy harvesting by the PV / thermoelectric system, the maximum power point tracking (MPPT) principle is targeted, aiming to operate the system close to peak power point. Under a uniform distribution of the solar irradiance, there is only one maximum power point (MPP), which easily can be e ﬃ ciently determined by any traditional MPPT method, such as the incremental conductance (INC). A di ﬀ erent situation will occur for the non-uniform distribution of solar irradiance, where more than one MPP will exist on the power versus voltage plot of the PV / thermoelectric system. The determination of the global MPP cannot be achieved by conventional methods. To deal with this issue the application of soft computing techniques based on optimization algorithms is used. However, MPPT based on optimization algorithms is very tedious and time consuming, especially under normal conditions. To solve this dilemma, this research examines a hybrid MPPT method, consisting of an incremental conductance (INC) approach and a moth-ﬂame optimizer (MFO), referred to as (INC-MFO) procedure, to reach high adaptability at di ﬀ erent environmental conditions. In this way, the combination of the two di ﬀ erent algorithms facilitates the utilization of the advantages of the two methods, thereby resulting in a faster speed tracking with uniform radiation distribution and a high accuracy in the case of a non-uniform distribution. It is very important to mention that the INC method is used to track the maximum power point under normal conditions, whereas the MFO optimizer is most relevant for the global search under partial shading. The obtained results revealed that the proposed strategy performed best in both of the dynamic and the steady-state conditions at uniform and non-uniform radiation. their low e ﬃ ciency and inconsistent power output due to the temperature ﬂuctuations. The TEG is a metal interconnected thermocouple consisting of a series p and n -types pellets.


Introduction
A rapid increase of utilized fossil fuels and the limitation on producing this fuel in a large number of countries around the world has led to the depletion of fossil fuel, increasing the price and causing environment problems such as acid rain, air pollution and global warming [1,2]. Photovoltaic (PV) system is considered an effective way to capture solar irradiance and directly converts it to electricity [3,4]. At the present time, PV systems are employed in different applications such as water pumping and charging vehicle batteries [5][6][7]. PV solar cell conversion efficiency is considered relatively low, within the range of 10-20% for commercial silicon cell and up to 40% for most multi-junction cells. For better overall PV energy efficiency generation, the application of cogeneration technique is widely attracting the researcher's attention, and this is achieved by capturing the waste heat as well and using it as an auxiliary energy production source. This can be done by incorporating a PV panel and thermoelectric generator (TEG) forming the so-called Hybrid photovoltaic/ thermoelectric generator (HPVTEG) system, which is increasingly attracting considerable interest among the issues of energy conversion efficiency improvement. HPVTEG system that integrates PV panel and TEG comprises double layer elements; the upper layer representing the PV panel, and a TEG at the bottom layer. TEG is fixed at the back of the PV module, which in turn utilizes the thermal waste to generate an additional electricity. The hot side of the TEG is connected to the backside of the PV panel and the cold side will be attached by heat sink as illustrated in Figure 1a. Figure 1b shows a typical schematic diagram of the HPVTEG. During the next few years, the TEGs are expected to attract profound applications as a prospective source of electricity [8]. They convert the embedded thermal energy into electricity, in such a way that, depending on the temperature difference across them, all possess a typical low conversion efficiency [9]. There are several TEGs arrangements in the electrical circuits to obtain the specified voltage and current. The major drawback to apply the TEG technology is their low efficiency and inconsistent power output due to the temperature fluctuations. The TEG is a metal interconnected thermocouple consisting of a series p and n-types pellets.

Introduction
A rapid increase of utilized fossil fuels and the limitation on producing this fuel in a large number of countries around the world has led to the depletion of fossil fuel, increasing the price and causing environment problems such as acid rain, air pollution and global warming [1,2]. Photovoltaic (PV) system is considered an effective way to capture solar irradiance and directly converts it to electricity [3,4]. At the present time, PV systems are employed in different applications such as water pumping and charging vehicle batteries [5,6,7]. PV solar cell conversion efficiency is considered relatively low, within the range of 10-20% for commercial silicon cell and up to 40% for most multijunction cells. For better overall PV energy efficiency generation, the application of cogeneration technique is widely attracting the researcher's attention, and this is achieved by capturing the waste heat as well and using it as an auxiliary energy production source. This can be done by incorporating a PV panel and thermoelectric generator (TEG) forming the so-called Hybrid photovoltaic/ thermoelectric generator (HPVTEG) system, which is increasingly attracting considerable interest among the issues of energy conversion efficiency improvement. HPVTEG system that integrates PV panel and TEG comprises double layer elements; the upper layer representing the PV panel, and a TEG at the bottom layer. TEG is fixed at the back of the PV module, which in turn utilizes the thermal waste to generate an additional electricity. The hot side of the TEG is connected to the backside of the PV panel and the cold side will be attached by heat sink as illustrated in Figure 1a. Figure 1b shows a typical schematic diagram of the HPVTEG. During the next few years, the TEGs are expected to attract profound applications as a prospective source of electricity [8]. They convert the embedded thermal energy into electricity, in such a way that, depending on the temperature difference across them, all possess a typical low conversion efficiency [9]. There are several TEGs arrangements in the electrical circuits to obtain the specified voltage and current. The major drawback to apply the TEG technology is their low efficiency and inconsistent power output due to the temperature fluctuations. The TEG is a metal interconnected thermocouple consisting of a series p and n-types pellets. The hybrid photovoltaic thermoelectric system is classified according to the way of their connection to a non-concentric and concentric, and their detailed classification is given in Figure 2. The hybrid photovoltaic thermoelectric system is classified according to the way of their connection to a non-concentric and concentric, and their detailed classification is given in Figure 2. Based on the temperature coefficient, and the concentration of the heat, the efficiency of the system will vary accordingly [10], and it was found that the hybrid model efficiency is better than that of a single PV solar cell only. The performance of HPVTEG system was investigated in several works of research [11][12][13][14][15][16][17][18][19]. Zhang et al. [11] presents the feasibility and features of the integrated system consisting of the PV solar cell and TEG. The obtained results confirm that the HPVTEG system is viable and the overall efficiency is enhanced. The theoretical simulation proved that the efficiency of hybrid system is improved by 30% by adding the TEG to the backside of the solar PV panel [18]. The enhancement of the conversion efficiency predicted by Sark was about 23% [17]. Experimentally, a 100 W prototype of a hybrid model was validated, tested, and implemented to an automobile [19]. In this work, the exhaust gas waste heat energy recovery for automobiles through TEG. New hybrid models were designed, constructed, into variances, with different solar cell and collectors, which were fixed on the hot and cold layers, respectively and then were tested. This PV-TEG-STC (solar collector) model used for producing electrical and thermal energy actually is facing two setbacks; the relatively lower TEG generated output power and the higher module cost [13]. These obstacles were overcome by increasing the solar collector exposure area and choosing a cheaper TEG. In addition, the hybrid system HPVTEG was implemented experimentally to hybrid electric vehicles proved that it was an efficient energy source [14]. With a sun-tracking arrangement, this model proved very efficient compared with the corresponding traditional units [15]. It is very important to notice that the hybrid model output is dependent on the temperature difference between the two sides of the TEG; once it is increased, the output will increase also. As a solution, Reference [16] suggested to optimize the coolant flow rate. The heat flux distributions and the heat transfer was numerically studied by using finite element method. The results show that the integrated design characteristic and performance was enhanced for both of TEG and solar cell. Table 1 summarizes the overall power and performance of different PV/TEG systems.  Based on the temperature coefficient, and the concentration of the heat, the efficiency of the system will vary accordingly [10], and it was found that the hybrid model efficiency is better than that of a single PV solar cell only. The performance of HPVTEG system was investigated in several works of research [11][12][13][14][15][16][17][18][19]. Zhang et al. [11] presents the feasibility and features of the integrated system consisting of the PV solar cell and TEG. The obtained results confirm that the HPVTEG system is viable and the overall efficiency is enhanced. The theoretical simulation proved that the efficiency of hybrid system is improved by 30% by adding the TEG to the backside of the solar PV panel [18]. The enhancement of the conversion efficiency predicted by Sark was about 23% [17]. Experimentally, a 100 W prototype of a hybrid model was validated, tested, and implemented to an automobile [19]. In this work, the exhaust gas waste heat energy recovery for automobiles through TEG. New hybrid models were designed, constructed, into variances, with different solar cell and collectors, which were fixed on the hot and cold layers, respectively and then were tested. This PV-TEG-STC (solar collector) model used for producing electrical and thermal energy actually is facing two setbacks; the relatively lower TEG generated output power and the higher module cost [13]. These obstacles were overcome by increasing the solar collector exposure area and choosing a cheaper TEG. In addition, the hybrid system HPVTEG was implemented experimentally to hybrid electric vehicles proved that it was an efficient energy source [14]. With a sun-tracking arrangement, this model proved very efficient compared with the corresponding traditional units [15]. It is very important to notice that the hybrid model output is dependent on the temperature difference between the two sides of the TEG; once it is increased, the output will increase also. As a solution, Reference [16] suggested to optimize the coolant flow rate. The heat flux distributions and the heat transfer was numerically studied by using finite element method. The results show that the integrated design characteristic and performance was enhanced for both of TEG and solar cell. Table 1 summarizes the overall power and performance of different PV/TEG systems. Terminating the HPVTEG with a load equivalent to the internal resistance of the system will facilitate to harvest a maximum power, and this is in compliance with the maximum power transfer principle, therefore, and accordingly the system operating point will be shifted to its peak value. The maximum power point tracker usually comprises a DC-to-DC power converter, controlled in such way that always maximizes the system output power, under the all-operating conditions. By this way a matching of the HPVTEG virtual load to its actual internal resistance by varying the duty cycle of the converter. If the HPVTEG has a direct connection with the load, then it's working point would be set by the load impedance, which indicates that the system will not produce its maximum power [32,33]. In general, this type of power tracking will enable an efficient interfacing of the HPVTEG system with the converter, to transfer a highest power at a fixed voltage. Various papers were dedicated to track the PV solar systems peak power point under variable environmental conditions. These techniques include; fractional open-circuit voltage, short-circuit current, the Artificial Neural Network (ANN) technique [33], Incremental Conductance Method, Hill Climbing Method, and the Fuzzy Logic control [32]. Conventional MPPT methods can easy extract the MPP under uniform distribution of solar radiation. Nevertheless, under partial shading condition (PSC), they cannot extract the global MPP since multiple local MPPs are exhibited on the power against voltage curve of HPVTEG system. Several MPPTs, based on modern optimization are proposed to solve this problem. These methods include; particle swarm optimization [34], cuckoo search [34], mine blast optimization [35], teaching learning based optimization [36], flower pollination and differential evolution [37]. Although the advantages of the several optimization algorithms have been expressed in the literature, it has been demonstrated by the No-Free-Lunch Theorem [38] that none of these algorithms can solve all optimization problems. Such theorem confirms the significance of recent optimizers in various applications since the efficacy of an optimizer to solve a set of problems does not guarantee its success in other application. The main drawback of soft computing based global MPPT techniques is that they are time consuming. Therefore, the main objective of this research is focused on building up a hybrid MPPT technique that combines incremental conductance (INC) and moth-flame optimizer by Mirjalili [39] to achieve better adaptability in various environments. The incorporation of the two different algorithms leads to the utilization of the advantages of both methods, thereby providing faster tracking speed with a uniform radiation distribution and high accuracy with a non-uniform radiation distribution.

HPVTEG System Components
The considered HPVTEG system mainly consists of a PV panel, TEG bank, DC/DC booster converter, control system, and the load as illustrated in Figure 3. As shown in Figure 3, to save cost, the TEG bank is connected in series with PV panel in order to use only one controller and one DC-DC converter instead of using two controllers with two DC-DC converters when TEG bank and PV panel are connected in parallel. The PV solar panel has three bypass diodes for minimizing and preventing the hotspot problems during partial shading. The MPP tracking methods are employed to increase the harvested energy from HPVTEG system under normal and abnormal conditions. They are used to adjust the DC-DC converter duty-cycle for controlling the operating voltage and current to extract maximum output from the hybrid system. TEG bank consisting of 72 units, divided into two sub-groups is adopted. Each group includes six strings, having six serially connected units each. Three switches, S1, S2 and S3, are used to switch between the serial and parallel connections of the two sub-groups inside the TEG bank. Figure 4 illustrates the flow chart of the proposed strategy. The proposed strategy uses the measured voltages across each bypass diode in order to control the configuration of TEG bank and identify the tracking mode of operation.
to extract maximum output from the hybrid system. TEG bank consisting of 72 units, divided into two sub-groups is adopted. Each group includes six strings, having six serially connected units each. Three switches, S1, S2 and S3, are used to switch between the serial and parallel connections of the two sub-groups inside the TEG bank. Figure 4 illustrates the flow chart of the proposed strategy. The proposed strategy uses the measured voltages across each bypass diode in order to control the configuration of TEG bank and identify the tracking mode of operation.

Solar Panel
For this research a Trina solar panel type (TSM-205DA01A.05) was adopted, which consists of 72 multi-crystalline silicon solar cells, connected in series with a power rating of 205 W at 1000 W/m 2 . Each 24 solar cells are shunted by a bypass diode [40]. These diodes are reversed biased during the uniform solar radiation distribution while under shadowing effect; they are forward biased and carry the current instead of the PV panel. Trina solar panel electrical specifications are given in Table 2. The

Solar Panel
For this research a Trina solar panel type (TSM-205DA01A.05) was adopted, which consists of 72 multi-crystalline silicon solar cells, connected in series with a power rating of 205 W at 1000 W/m 2 . Each 24 solar cells are shunted by a bypass diode [40]. These diodes are reversed biased during the uniform solar radiation distribution while under shadowing effect; they are forward biased and carry the current instead of the PV panel. Trina solar panel electrical specifications are given in Table 2. The power vs. voltage curves under uniform and non-uniform solar radiation distributions are shown in Figure 5, while the test results of a PV panel is shown in Table 3. Figure 5 shows that the power-voltage characteristics have so many maximum points, equal to the number of irradiance levels incident on the PV module in the case of the partial shading condition. Whereas it has one single MPP at uniform solar radiation distribution.

Thermoelectric Generator Modeling
TEG consists of a thermocouple, made of a metal interconnected p and n-type pellets. Normally the thermocouples are arranged in serially connected arrays, to increase the operating voltage, and thermally in parallel to decrease the thermal resistance, and they are interlapped between two ceramic heat sinks at its two ends for a uniform thermal expansion. One of the two sinks has a higher temperature, and will be known as the hot side of the TEG; while the other has a lower temperature,

Thermoelectric Generator Modeling
TEG consists of a thermocouple, made of a metal interconnected p and n-type pellets. Normally the thermocouples are arranged in serially connected arrays, to increase the operating voltage, and thermally in parallel to decrease the thermal resistance, and they are interlapped between two ceramic heat sinks at its two ends for a uniform thermal expansion. One of the two sinks has a higher temperature, and will be known as the hot side of the TEG; while the other has a lower temperature, designated as the cold side [41].
The phenomenon whereby the presence of a temperature difference between any two semiconductors, a thermos-electromotive force will be induced, and this is known as a Seebeck Effect, which is considered to be the basic principle of the TEG [42]. The open-circuit thermo-electro motive force V OC is given by the following equation: where: T h and T c -the hot and cold side temperatures.
The Seebeck coefficient characterizes the induced e.m.f, when the temperature at the junction is increased by 1 • C. Usually metals possess smaller coefficient values. Another physical effect, known as the Peltier Effect, states that if a direct current flows through a junction between two different metal layers, one junction will be heated and the other will be cooled, the heat absorption and dissipation direction depends on the current polarity. The Peltier coefficient is then defined as; where: P P denotes heat-transfer rate; I teg denotes DC current flowing in the TEG.
The Peltier coefficient characterizes the amount of the junction heat either at heating or cooling, meanwhile, the reversible absorption or liberation of heat in a homogeneous material simultaneously exposed to a thermal gradient upon the flow of an electric current, defines the so-called Thomson Effect.
The conductor's heat dissipation when the current flows toward the higher temperature is given by: where: τ denotes Thomson coefficient.
Both the Seebeck and Peltier coefficients characterize the metal junctions, so they can be determined only for a pair of junctions, while the Thomson coefficient characterizes an individual conductor. Peltier coefficient can be formulated as follows: The current flow through the TEG will cause an additional heat generation in the thermocouples, and this effect of heat dissipation is observed in both sides at different temperatures, but with the same amount of energy as follows: where Q j is the Joule heating.
The steady-state analysis at the both sides of the TEG, based on the energy equilibrium concept, the absorbed heat generated by the thermal load and the liberated heat removed by the heat sink can be represented respectively by the following equations [42]: where: k tc is the thermal conductivity.
The electrical power is equal to the difference between heat flow at the hot and the cold sides [42]: where: R int and I teg -the TEG electrical resistance and current respectively.
With reference to Kirchhoff's voltage law, the TEG terminal voltage (V teg ) will be given by: The datasheet of the commercial TEGs includes the following parameters: T h and T c , hot and cold-sides temperatures; P m , the matched load power; R int , matched load resistance; V m the load voltage at the matched load. Obviously, it is very easy to calculate the equivalent circuit parameters from the datasheet. The internal resistance R int and the Seebeck coefficient α of a TEG can be given by: In this work, a 1 Watt thermoelectric module has been used. Table 4 shows the electrical specifications and dimensions of a TE-MOD-1W2V-40S:

Conventional MPPT based on Incremental Conductance
The INC technique is widely considered the most convenient algorithm for the maximum power tracking under uniform distribution of solar irradiance [44]. Under normal operation of PV system, the PV power versus PV voltage graph contains single the peak point. This point actually represents the MPP at which the PV module produces its maximum power. The core idea of the INC method is that PV power derivative w.r.t its voltage is zero at the MPP [44]. The derivative of power against voltage can be represented by Equation (13) and accordingly, the error signal may be calculated as in Equation (14) [45]. Figure 6 shows INC method flowchart.
and the error signal will be given as: Accordingly, tracking the MPP needs the following strategy: where, K is the integrator gain. It is preferable to start the tracking process with large step size for a quick allocation of the peak power point, then gradually reduce it, when the working point is near to the specified, to reduce oscillations around it. INC−MPPT is implemented with the help of an integrator, which gain is K and fed by the error signal.
where, K is the integrator gain. It is preferable to start the tracking process with large step size for a quick allocation of the peak power point, then gradually reduce it, when the working point is near to the specified, to reduce oscillations around it. INC−MPPT is implemented with the help of an integrator, which gain is K and fed by the error signal.

Moth-Flame Optimizer
For many applications such as the optimal power flow problems, MFO has proved to be the most relevant, compared with many other optimization techniques [46]. The navigating mechanism of the moths, known as the transverse orientation principle, is simply described by the MFO. MFO that proposed by Mirjalili [39] has been used effectively in the field of tracking MPP by authors in Reference [47]. Moths are night-flying insects, sensing the moonlight, due to their unique navigation instinct. They maintain a constant angle in reference to the moon location, during flight, keeping a straight-line trajectory. This is named transverse orientation, but in the case of non-natural light sources, the moths do not go straight, but perform a spiral pass, although this source now is deemed a new attracting point for the moths.

Moth-Flame Optimizer
For many applications such as the optimal power flow problems, MFO has proved to be the most relevant, compared with many other optimization techniques [46]. The navigating mechanism of the moths, known as the transverse orientation principle, is simply described by the MFO. MFO that proposed by Mirjalili [39] has been used effectively in the field of tracking MPP by authors in Reference [47]. Moths are night-flying insects, sensing the moonlight, due to their unique navigation instinct. They maintain a constant angle in reference to the moon location, during flight, keeping a straight-line trajectory. This is named transverse orientation, but in the case of non-natural light sources, the moths do not go straight, but perform a spiral pass, although this source now is deemed a new attracting point for the moths.
The MFO technique is considered to be a population-based algorithm, so it can be simulated by the matrix below: where n is the number of moths and d is the size of the search space in which the moths and flames positions vectors operate. The corresponding moths' objective functions OM are given by the following matrix: Similarly, the following flames matrix can be set as: The MFO procedure includes three-raw approximation functions indicated below: where: "I" is the initial function that yields un-specified moth population, while their corresponding objective functions are: M(i, j) = (ul(i) − ll(i)) · rand(n, d) + ll(i) (20) rand: random distribution function.
where: ul and ll denote the upper and lower thresholds of the variables, respectively. According to the optimization procedure, this function should be first run to the end, followed by the population function "P" until the termination criteria "T" is met. Therefore, the function "P" will guide the moths around the search location. The main moth position updating mechanism with respect to the flame is the logarithmic spiral function.
This position reference to the flame is updated as follows: The logarithmic spiral function can be represented by the following formula: where: M i is the ith moth order; F j is the order of the jth flame; and S denotes the spiral function, D i is the distance between the ith moth and the jth flame; b is a constant; and t is a random number in [r, 1]. Then, r is the adaptive convergence constant that linearly decreases from −1 to −2 to fast-track convergence around the flames over the path of iterations. D is given by: Exploration occurs at the moment when the consecutive position lies outside the enclosed space between the moth and the flame. It is said that the exploitation is converged, once the next position lies inside the space, shown by the arrow indicated by Equation (25) shows that the number of flames reasonably decreases over the iterations for assuring the balance between the exploration and exploitation. Therefore, the moth positions are updated only according to the best flame in the final iteration steps defined by: where: "l" is the current number of iteration, N f is the maximum number of flames, and T max represents the maximum number of iterations.
The converter duty cycle is the most important key to control the DC-DC converter for maximizing the harvested energy from the HPVTEG system. Therefore, it has been selected to be the decision variable during the optimization process to maximize the output power of HPVTEG system. The relationship between the input voltage and load voltage of DC-DC converter can be formulated as follows; (26) where V HPVTEG is the HPVTEG voltage, and V load is the load voltage, D is the duty cycle. The output power of HPVTEG system that represents the objective function can be estimated as follows; The procedure of MFO optimizer to extract the global maximum power of HPVTEG system under partial shading condition can be summarized as follows. Since MFO contains both moths and the flames populations that are moving in the search space. Equations (16) and (18) are used to show that moths and the flames are grouped in the matrixes M and F with same dimension. MFO optimizer starts the optimization process by n number of moths for every decision variable. For the case study, only one variable is considered. This variable is the duty cycle of DC-DC converter. The moths' positions initialization are based on Equation (20). For each value of moths (duty cycles), the HPVTEG system is operated and the corresponding system voltage and current are sensed. Then based on the measured current and voltage, the system power (OF) is estimated by Equation (27). Updating the moths' position with considering the flames positions will be carried out using Equation (23).

Solar PV Panel Thermal Modeling
In order to select the proper type of TEG, the temperature difference should be determined. Computational Fluid Dynamics (CFD) commercial package ANSYSFLUENT 15 is used to predict the temperature difference between the solar panel and the ambient. The solar panel geometry used in this study is based on the TALLMAX framed 72 cell module produce by Trina Solar Company. Only one solar panel is considered, with panel dimensions of 1.956 × 0.922 × 0.004 m 3 . The maximum solar radiation can be harvested by facing the PV panel toward the south with inclination angle of 20.5 • [3]. CFD predictions as a solution of a three-dimensional Navier-Stokes equations with Analysis of Systems (ANSYS) Fluent 15 is a very strong and effective tool. The computational model was specified to be three-dimensional, therefore a simple algorithm is used to couple fluid pressure and velocity. The residual converged solution of the continuous velocity component, turbulence kinetic energy is for energy below 10 −6 , and the turbulence dissipation rate is 10 −3 . The mesh is a uniform hexahedral element, meshing was specified along the boundary and swept later to cover the entire model volume. The computations in the present paper have mainly been carried out using k-ε model, (under the two-equation model category). The temperature boundary condition, based on the average highest temperature during the year in Wadi Addwaser city as a case study, Riyadh Governorate, Kingdom of Saudi Arabia is (20.504 • latitude, 45.2 • longitude) according to Reference [48] the average highest temperature is 30 • C (303 K). Therefore, the four sidewalls of the solar panel are set to constant temperature 30 • C (303 K), while the upper face is subjected to solar radiation of 800 W/m 2 . To ensure that the results are independent of the mesh size, mesh dependence study is conducted by monitoring the temperature at three different positions across the center of the solar panel at (x = 0.85, y = 0.461, z = 0.002), (x = 1.5, y = 0.461, z = 0.002) and (x = 0.4, y = 0.461, z = 0.002). Three different number of mesh cells are adopted: 12568, 21216 and 25376. The temperature values at the three selected positions for different number of mesh cells are plotted in Figure 7; it is clear that there is a complete matching between the results of the last two meshes indicating that the solution is independent of mesh size for either of these two mesh, namely 21216 and 25376. The results obtained in this study are based on number of mesh cells of 21216. panel at (x = 0.85, y = 0.461, z = 0.002), (x = 1.5, y = 0.461, z = 0.002) and (x = 0.4, y = 0.461, z = 0.002). Three different number of mesh cells are adopted: 12568, 21216 and 25376. The temperature values at the three selected positions for different number of mesh cells are plotted in Figure 7; it is clear that there is a complete matching between the results of the last two meshes indicating that the solution is independent of mesh size for either of these two mesh, namely 21216 and 25376. The results obtained in this study are based on number of mesh cells of 21216.

Results and Discussions
As mentioned before, in this research a TEG bank consisting of 72 units, divided into two subgroups is adopted. Each group includes six strings, having six serially connected units each. A MATLAB code is developed in order to simulate the bank and test it at different configurations. These two groups can be connected in series or in parallel based on system tracking mode of operation. Based on the temperature distribution of a solar panel shown in Figure 7, the hot side temperature of the TEG is assumed to be equal the highest temperature of 73 °C occupying the central part of the solar panel. Whereas the cold side of TEG is equal to the atmosphere temperature of 30 °C. Therefore, the difference temperature between the two sides is 43 °C. For this difference in temperature, when the two groups of TEGs are connected in series, the voltage, current and power at MPP will be 13.80 V, 2.11 A and 29.12 W respectively. Table 5 and Figure 9 show the power, current and voltage of TEG bank under different conditions.

Results and Discussions
As mentioned before, in this research a TEG bank consisting of 72 units, divided into two sub-groups is adopted. Each group includes six strings, having six serially connected units each. A MATLAB code is developed in order to simulate the bank and test it at different configurations. These two groups can be connected in series or in parallel based on system tracking mode of operation. Based on the temperature distribution of a solar panel shown in Figure 7, the hot side temperature of the TEG is assumed to be equal the highest temperature of 73 • C occupying the central part of the solar panel. Whereas the cold side of TEG is equal to the atmosphere temperature of 30 • C. Therefore, the difference temperature between the two sides is 43 • C. For this difference in temperature, when the two groups of TEGs are connected in series, the voltage, current and power at MPP will be 13.80 V, 2.11 A and 29.12 W respectively. Table 5 and Figure 9 show the power, current and voltage of TEG bank under different conditions.  Four different cases are used to validate the proposed strategy. In the first case, the solar radiation distribution is uniform and equal to 800 W/m 2 . In this case, and based on Table 6 and Figure  10, the PV panel current is 4.2 A. Therefore, the parallel connection of the two TEG sub-groups is better than the series connection. The situation is reversed in case of 400 W/m 2 uniform distribution, the series connection of the two TEG sub-groups is better than the parallel connection. In this case the power increased by 7.4 %. During partial shading condition, the PV panel current is reduced. Accordingly, the series connection of the two TEG groups is better than the series connection.  Four different cases are used to validate the proposed strategy. In the first case, the solar radiation distribution is uniform and equal to 800 W/m 2 . In this case, and based on Table 6 and Figure 10, the PV panel current is 4.2 A. Therefore, the parallel connection of the two TEG sub-groups is better than the series connection. The situation is reversed in case of 400 W/m 2 uniform distribution, the series connection of the two TEG sub-groups is better than the parallel connection. In this case the power increased by 7.4%. During partial shading condition, the PV panel current is reduced. Accordingly, the series connection of the two TEG groups is better than the series connection.   For the first case, the solar radiation is uniform so there is only single MPP of 192.44 W in the power against voltage curve. Based on the proposed strategy, the voltage across each bypass diode is approximately same. Therefore, the INC method is called and started the tracking process. Figure 11  For the first case, the solar radiation is uniform so there is only single MPP of 192.44 W in the power against voltage curve. Based on the proposed strategy, the voltage across each bypass diode is approximately same. Therefore, the INC method is called and started the tracking process. Figure  11 illustrated the time variation of a HPVTEG system power, voltage, current and boost converter duty cycle, from which it is clear that the INC method is effectively, allocates the maximum power point with a tracking time of 0.1 s. For the fourth configuration, three different solar irradiances levels of (800, 500 and 300 W/m 2 ) are subjected to the PV panel. Under this situation, three MPPs located on the power/voltage curves since the number of maximum power point equal to the number of levels of solar irradiance. One global maximum power point of 98.11 W is allocated on the center power/voltage curve. Figure 12 shows the dynamic performance of INC method and MFO method under partial shading condition. From this figure, it can be noted that the MFO bypasses the first local peak power point (83.5 W) and fixed the optimal one of 98.11 W. The conventional INC-based tracker caught the first local point of 83.5 W, due to its disability to differentiate between the local and optimal maximum power point. At this situation the tracking efficiency increases by 17.5 % using MFO in comparison with the INC method. The boost convertor duty cycle updating for each iteration in the case of MFO method is shown in Figure 13. Considering this figure, one can see that the optimizer starts the process by three values of duty cycles (0.3, 0.6 and 0.8). These values are updated for each iteration until the optimal duty cycle corresponding to global maximum power is determined. More case studies can be found For the fourth configuration, three different solar irradiances levels of (800, 500 and 300 W/m 2 ) are subjected to the PV panel. Under this situation, three MPPs located on the power/voltage curves since the number of maximum power point equal to the number of levels of solar irradiance. One global maximum power point of 98.11 W is allocated on the center power/voltage curve. Figure 12 shows the dynamic performance of INC method and MFO method under partial shading condition. From this figure, it can be noted that the MFO bypasses the first local peak power point (83.5 W) and fixed the optimal one of 98.11 W. The conventional INC-based tracker caught the first local point of 83.5 W, due to its disability to differentiate between the local and optimal maximum power point. At this situation the tracking efficiency increases by 17.5% using MFO in comparison with the INC method. The boost convertor duty cycle updating for each iteration in the case of MFO method is shown in Figure 13. Considering this figure, one can see that the optimizer starts the process by three values of duty cycles (0.3, 0.6 and 0.8). These values are updated for each iteration until the optimal duty cycle corresponding to global maximum power is determined. More case studies can be found in Appendix A (Table A1, Figures A1 and A2).

Conclusion
Improving energy conversion efficiency of PV system through thermoelectric generator (TEG) integration is presented in this paper. This is achieved by capturing the waste heat and using it as an auxiliary energy production source. The temperature difference between cold and hot side of the TEG is predicted by using CFD commercial package ANSYS FLUENT 15. The modeling and performance evaluation of system are carried out using the MATLAB software (version 2018, mathworks: MA, US). The system maximum power is harvested through a proposed a hybrid MPPT algorithm (INC-MFO) combining incremental conductance (INC) and moth-flame optimizer (MFO) to achieve better adaptability in various environment. The incorporation of two different algorithms combines the advantages of both, thereby providing faster tracking speed with uniform radiation distribution and high accuracy with non-uniform radiation distribution. The proposed strategy calls the INC method to track MPP under normal condition and calls MFO optimizer to catch the optimal MPP under nonuniform distribution. The obtained results reveal that the proposed strategy performed best in both

Conclusion
Improving energy conversion efficiency of PV system through thermoelectric generator (TEG) integration is presented in this paper. This is achieved by capturing the waste heat and using it as an auxiliary energy production source. The temperature difference between cold and hot side of the TEG is predicted by using CFD commercial package ANSYS FLUENT 15. The modeling and performance evaluation of system are carried out using the MATLAB software (version 2018, mathworks: MA, US). The system maximum power is harvested through a proposed a hybrid MPPT algorithm (INC-MFO) combining incremental conductance (INC) and moth-flame optimizer (MFO) to achieve better adaptability in various environment. The incorporation of two different algorithms combines the advantages of both, thereby providing faster tracking speed with uniform radiation distribution and high accuracy with non-uniform radiation distribution. The proposed strategy calls the INC method to track MPP under normal condition and calls MFO optimizer to catch the optimal MPP under nonuniform distribution. The obtained results reveal that the proposed strategy performed best in both

Conclusion
Improving energy conversion efficiency of PV system through thermoelectric generator (TEG) integration is presented in this paper. This is achieved by capturing the waste heat and using it as an auxiliary energy production source. The temperature difference between cold and hot side of the TEG is predicted by using CFD commercial package ANSYS FLUENT 15. The modeling and performance evaluation of system are carried out using the MATLAB software (version 2018, mathworks: MA, US). The system maximum power is harvested through a proposed a hybrid MPPT algorithm (INC-MFO) combining incremental conductance (INC) and moth-flame optimizer (MFO) to achieve better adaptability in various environment. The incorporation of two different algorithms combines the advantages of both, thereby providing faster tracking speed with uniform radiation distribution and high accuracy with non-uniform radiation distribution. The proposed strategy calls the INC method to track MPP under normal condition and calls MFO optimizer to catch the optimal MPP under non-uniform distribution. The obtained results reveal that the proposed strategy performed best in both dynamic response and steady-state in cases of uniform and non-uniform irradiance distribution.

Conflicts of Interest:
The authors declare no conflicts of interest.        Figure A2. The performance of HPVTEG under different scenarios.