Parameter Identiﬁcation of Optimized Fractional Maximum Power Point Tracking for Thermoelectric Generation Systems Using Manta Ray Foraging Optimization

: Thermoelectric generation systems (TEGSs) are used to convert temperature difference and heat ﬂow into DC power based on the Seebeck theorem. The basic unit of TEGS is the thermoelectric module (TEM). TEGSs have gained increasing interest in the research ﬁelds of sustainable energy. The output power from TEM is mostly reliant on differential temperature between the hot and cold sides of the TEM added to the value of the load. As such, a robust MPPT strategy (MPPTS) is required to ensure that the TEGS is operating near to the MPP while varying the operating conditions. Two main drawbacks may occur in the conventional MPPTSs: low dynamic response, such as in the incremental resistance (INR) method, and oscillations around MPP at steady state, such as in the hill climbing (HC) method. In the current research work, an optimized fractional MPPTS is developed to improve the tracking performance of the TEGS, and remove the two drawbacks of the conventional MPPTSs. The proposed strategy is based on fractional order control (FOC). The main advantage of FOC is that it offers extra ﬂexible time and frequency responses of the control system consent for better and robust performance. The optimal parameters of the optimized fractional MPPTS are identiﬁed by a manta ray foraging optimization (MRFO). To verify the robustness of the MRFO, the obtained results are compared with ten other algorithms: particle swarm optimization; whale optimization algorithm; Harris hawks optimization; heap-based optimizer; gradient-based optimizer; grey wolf optimizer; slime mould algorithm; genetic algorithm; seagull optimization algorithm (SOA); and tunicate swarm algorithm. The maximum average cost function of 4.92934 kWh has been achieved by MRFO, followed by SOA (4.5721 kWh). The lowest STD of 0.04867 was also accomplished by MRFO. The maximum efﬁciency of 99.46% has been obtained by MRFO, whereas the lowest efﬁciency of 74.01% was obtained by GA. Finally, the main ﬁndings proved the superiority of optimized fractional MPPTS compared with conventional methods for both steady-state and dynamic responses.


Introduction
A massive part of heat energy consumption is mostly wasted by the industrial sectors and transportation. Thus, the recovery system of wasted heat gives a significant portion of the daily needs from energy, specifically for the transportation sector [1,2]. The recovery of the wasted heat can be accomplished by utilizing thermoelectric generation (TEG) technology [3]. Thermoelectric generators convert heat directly into electricity. They are 1.
A new OFMPPTS to maximize the harvested the energy of the TEG system is developed. The developed approach utilizes the fractional control that offers better and robust performance.

2.
A new application of the MRFO algorithm is proposed for determining the optimal parameters of OFMPPTS. The suggested design by using integrated features of the MRFO, and fractional order control introduces a promising solution of MPPT in TEG systems. 3.
Improved efficiency of the TEG system for different applications is suggested in this research paper through simple and effective FOMPPTS by using only electrical measured signals from the TEG source. Compared to the other existing methods, reduced output power fluctuations in output power, voltage, in addition to current waveforms from the TEG system, are taken by the suggested FOMPPTS.
The rest of the paper is arranged as following. Section 2 presents the modelling of the thermoelectric generator. The manta ray foraging optimization and the suggested optimized fraction MPPT are detailed in Sections 3 and 4, respectively. The results and discussions for different tested operating points and case studies are introduced in Section 5. Lastly, Section 6 presents the main findings of the research paper.

Modelling of Thermoelectric Generator
Thermoelectric generator semiconductor devices use the Seebeck effect to produce a DC voltage. The generated voltage drives electrical current, and produces useful DC power at a load. The general view of the TEG is shown in Figure 1a. The fundamental building block of the TEG is a thermocouple that includes p-type and n-type semiconductors, as demonstrated in Figure 1b. To form a TEG module, some of the p-type and n-type groups are combined in series to raise the DC voltage and, accordingly, the output power. Simultaneously, these groups are thermally in parallel to reduce the thermal resistance [23,24]. The couples are planted between two parallel ceramic plates. These plates create a hot side and cold side, as described in Figure 1b, and give structural rigidity, a flat surface for mounting, and a dielectric layer to prevent electrical short circuits.
The main contributions in this paper are outlined as follows: 1. A new OFMPPTS to maximize the harvested the energy of the TEG system is developed. The developed approach utilizes the fractional control that offers better and robust performance. 2. A new application of the MRFO algorithm is proposed for determining the optimal parameters of OFMPPTS. The suggested design by using integrated features of the MRFO, and fractional order control introduces a promising solution of MPPT in TEG systems. 3. Improved efficiency of the TEG system for different applications is suggested in this research paper through simple and effective FOMPPTS by using only electrical measured signals from the TEG source. Compared to the other existing methods, reduced output power fluctuations in output power, voltage, in addition to current waveforms from the TEG system, are taken by the suggested FOMPPTS.
The rest of the paper is arranged as following. Section 2 presents the modelling of the thermoelectric generator. The manta ray foraging optimization and the suggested optimized fraction MPPT are detailed in Section 3 and Section 4, respectively. The results and discussions for different tested operating points and case studies are introduced in Section 5. Lastly, Section 6 presents the main findings of the research paper.

Modelling of Thermoelectric Generator
Thermoelectric generator semiconductor devices use the Seebeck effect to produce a DC voltage. The generated voltage drives electrical current, and produces useful DC power at a load. The general view of the TEG is shown in Figure 1a. The fundamental building block of the TEG is a thermocouple that includes p-type and n-type semiconductors, as demonstrated in Figure 1b. To form a TEG module, some of the p-type and n-type groups are combined in series to raise the DC voltage and, accordingly, the output power. Simultaneously, these groups are thermally in parallel to reduce the thermal resistance [23,24]. The couples are planted between two parallel ceramic plates. These plates create a hot side and cold side, as described in Figure 1b, and give structural rigidity, a flat surface for mounting, and a dielectric layer to prevent electrical short circuits. When a temperature difference is occurred across the hot and cold sides, a voltage potential will appear based on the Seebeck effect. The Seebeck effect occurs due to the movement of charge carriers within the semiconductors. Figure 2 presents the equivalent-circuit of the TEG unit. The open circuit voltage as a source (V O.C ) is series connected with an internal resistance (R int ) [23,24]. When a temperature difference is occurred across the hot and cold sides, a voltage potential will appear based on the Seebeck effect. The Seebeck effect occurs due to the movement of charge carriers within the semiconductors. Figure 2 presents the equivalentcircuit of the TEG unit. The open circuit voltage as a source (VO.C) is series connected with an internal resistance (Rint) [23,24]. The value of open-circuit voltage is mainly dependent on the temperature difference (ΔT) across of the two sides of the TEG. It can be estimated using the following relation: where denotes the Seebeck coefficient, and denote the hot and cold side temperatures, respectively.
Based on the heat-energy-equilibrium theory [25], the heat flow from the high-temperature-side ( ) to the low-temperature-side ( ) can be estimated as follows: where K denotes the thermal-conductivity coefficient. The TEG output power, Pteg, is mainly dependent on the difference of the heat energy on both sides of the TEG. It can be calculated by the following equation: Considering Figure 2, and applying KVL, the output voltage of the TEG, , can be expressed by the following relation: In sum, the TEG output power, Pteg, can be expressed as follows:

MRFO Algorithm
Zhao et al. [26] proposed manta ray foraging optimization (MRFO) as an innovative technique. MRFO mimicked the manta ray's foraging tactic when catching their food. The creatures follow three stages until they reach their preys, including chain, cyclone, and somersault foraging. The chain foraging is the first step that may be followed, where the The value of open-circuit voltage is mainly dependent on the temperature difference (∆T) across of the two sides of the TEG. It can be estimated using the following relation: where α sb denotes the Seebeck coefficient, T hot and T cold denote the hot and cold side temperatures, respectively. Based on the heat-energy-equilibrium theory [25], the heat flow from the hightemperature-side (Q h ) to the low-temperature-side (Q c ) can be estimated as follows: where K denotes the thermal-conductivity coefficient. The TEG output power, P teg , is mainly dependent on the difference of the heat energy on both sides of the TEG. It can be calculated by the following equation: Considering Figure 2, and applying KVL, the output voltage of the TEG, V teg , can be expressed by the following relation: In sum, the TEG output power, P teg , can be expressed as follows:

MRFO Algorithm
Zhao et al. [26] proposed manta ray foraging optimization (MRFO) as an innovative technique. MRFO mimicked the manta ray's foraging tactic when catching their food. The creatures follow three stages until they reach their preys, including chain, cyclone, and somersault foraging. The chain foraging is the first step that may be followed, where the creatures target to find highly concentrated plankton positions. Therefore, each agent of the creatures moves toward the food location. Zhao et al. [26] modelled this stage mathematically as described below: where M xi (t) represents the i th creature location at iteration t, r n is a random number, M best (t) is the best position obtained so far, N is total creature number, and the symbol of α is a weighting factor, that is given by Equation (8) as below: Considering Equation (7), the position for i th creature, except that of the first one, is reliant on the position from (i − 1) th agent M t x i and the best one (M t best ). While the manta ray creatures recognize the position of the plankton spot, they amalgamate, forming long chains, then swim in the direction of the prey in a spiral shape, and towards the manta ray in front of them. That stage is called cyclone foraging. Zhao et al. [26] mathematically described the spiral shaped as below: where ω refers to a random number, the mathematical formula of the cyclone stage is modelled as below: where β is a weighting operator that is given below: where t and T are the current and total number of iterations, respectively. The symbol of r 1 refers to a random number. Zhao et al. [26] added an auxiliary equation to enhance the exploration stage of the MRFO optimizer. During the auxiliary equation, the agents were forced to discover a new random location in the search landscape. The auxiliary formula is described as below: where Ub and Lb denote the upper and lower limits, M xrand denotes a random position. For the third and last stage of somersault foraging, every creature swims around the food back and forth, and somersault (tumble) to new position, as expressed below: where S M is a factor of somersault utilized to assign the manta rays somersault range, r 2 and r 3 are random values. In MRFO's final stage, every individual can move to any location in the search space between its current position and an asymmetrical one around the food. Thus, the distance Mathematics 2021, 9, 2971 6 of 18 between the manta ray and the food location is minimized throughout this stage, which refers to a convergence process to the optimal solution. Consequently, the factor S M range is adaptively dismissed through the iterations. The main steps of MRFO applied in the maximization procedure are demonstrated in Algorithm 1.

Optimized Fractional MPPT Strategy
To improve the dynamic tracking behavior of the incremental conductance (INC) method, a fraction control is implemented. The mathematical model of the INC method is extracted from the fact that the derivative of the TEG power related to the TEG voltage (dP/dV) is zero at the MPP. Therefore, the slope of the TEG power characteristics is zero at the MPP, positive for values located on the left side of the MPP, and negative for amounts located on the right side. Consequently, the error signal can be defined by the following relation: where V and I denote the instant voltage and current, respectively, and n is the increment time.
The INC method can be implemented using a discrete integrator with the error signal as the input.
The optimized fractional MPPTS executes a fractional PI as an alternative of an integer integrator, as presented in Figure 3. Fractional order control includes a non-integer order that has some merits, including flexibility in design and high robustness. The fractional PI transfer function can be defined as follows [27]: where T c represents the controller transfer function, the integral and proportional are, respectively, k I and k p , and˘is the fractional order.
where represents the controller transfer function, the integral and proportional are, respectively, and , and is the fractional order.

Results and Discussion
To investigate the TEG power of altering the temperature differences, a computer MATLAB program software was developed. The TEG model of TEG1-12611-6.0 is used in the simulations. For a single TEG, the nominal current, voltage, and power are 3.4 A, 4.2 V, and 14.6 W, respectively, at a hot side temperature of 300 °C, and a cold side temperature of 30 °C. Figure 4a displays the power versus current of the TEG through changing the hot side temperature by maintaining the cold side at a temperature of 30 °C. Referring to this figure, the estimated datasets and manufacturing datasets are extremely satisfying. Furthermore, the location of MPP mostly varies with the difference temperature between both cold and hot sides of the TEG. Power versus current of the TEG with a changing cold side temperature at a fixed hot side temperature of 300 °C is given in Figure  4b. Figure 5 shows the matched load parameters (power, current, voltage, and resistance) when varying the hot side and cold side temperature values.

Results and Discussion
To investigate the TEG power of altering the temperature differences, a computer MATLAB program software was developed. The TEG model of TEG1-12611-6.0 is used in the simulations. For a single TEG, the nominal current, voltage, and power are 3.4 A, 4.2 V, and 14.6 W, respectively, at a hot side temperature of 300 • C, and a cold side temperature of 30 • C. Figure 4a displays the power versus current of the TEG through changing the hot side temperature by maintaining the cold side at a temperature of 30 • C. Referring to this figure, the estimated datasets and manufacturing datasets are extremely satisfying. Furthermore, the location of MPP mostly varies with the difference temperature between both cold and hot sides of the TEG. Power versus current of the TEG with a changing cold side temperature at a fixed hot side temperature of 300 • C is given in Figure 4b. Figure 5 shows the matched load parameters (power, current, voltage, and resistance) when varying the hot side and cold side temperature values.
where represents the controller transfer function, the integral and proportional are, respectively, and , and is the fractional order.

Results and Discussion
To investigate the TEG power of altering the temperature differences, a computer MATLAB program software was developed. The TEG model of TEG1-12611-6.0 is used in the simulations. For a single TEG, the nominal current, voltage, and power are 3.4 A, 4.2 V, and 14.6 W, respectively, at a hot side temperature of 300 °C, and a cold side temperature of 30 °C. Figure 4a displays the power versus current of the TEG through changing the hot side temperature by maintaining the cold side at a temperature of 30 °C. Referring to this figure, the estimated datasets and manufacturing datasets are extremely satisfying. Furthermore, the location of MPP mostly varies with the difference temperature between both cold and hot sides of the TEG. Power versus current of the TEG with a changing cold side temperature at a fixed hot side temperature of 300 °C is given in Figure  4b. Figure 5 shows the matched load parameters (power, current, voltage, and resistance) when varying the hot side and cold side temperature values.

Evaluation of MRFO
Ten optimization algorithms: WOA; HHO; SMA; PSO; HBO; GBO; GA; GWO; SOA; and TSA were used in comparison with the proposed MRFO to identify the best parameters of OFMPPTS. For the duration of the optimization process, the parameters of OFMP-PTS, including integration gain, proportional gain, and fraction order, were selected to be the decision variables, but the cost function that is necessary to be maximum is the harvested energy (kWh) of the TEG. For fair comparison, the population size, and the max. no. of iterations for all optimization algorisms are set to 10 and 25, respectively. In addition, in this paper, the default parameter settings for the MRFO and counterparts are utilized, which is a fair and appropriate practice. The optimal parameters of OFMPPTS using different optimization algorithms are given in Table 1. To confirm the consistency of the considered algorithms, everyone executed 30 runs. The statistical analysis of these optimizers is also shown in Table 2.

Evaluation of MRFO
Ten optimization algorithms: WOA; HHO; SMA; PSO; HBO; GBO; GA; GWO; SOA; and TSA were used in comparison with the proposed MRFO to identify the best parameters of OFMPPTS. For the duration of the optimization process, the parameters of OFMPPTS, including integration gain, proportional gain, and fraction order, were selected to be the decision variables, but the cost function that is necessary to be maximum is the harvested energy (kWh) of the TEG. For fair comparison, the population size, and the max. no. of iterations for all optimization algorisms are set to 10 and 25, respectively. In addition, in this paper, the default parameter settings for the MRFO and counterparts are utilized, which is a fair and appropriate practice. The optimal parameters of OFMPPTS using different optimization algorithms are given in Table 1. To confirm the consistency of the considered algorithms, everyone executed 30 runs. The statistical analysis of these optimizers is also shown in Table 2.  Referring to Table 1, the performance of MRFO is the best among all considered algorithms. The average cost function values varied between 4.92934 kWh and 3.59264 kWh. The maximum average cost function of 4.92934 kWh has been achieved by MRFO, followed by SOA (4.5721 kWh). The STD values varied between 0.04867 and 1.18114. The lowest STD was also accomplished by MRFO. Regarding the efficiency, the maximum efficiency of 99.46% has been obtained by MRFO, whereas the lowest efficiency of 74.01% was obtained by GA. The details of 30 runs using different algorithms are presented in Table 2. Using these data, a Wilcoxon sign rank test and Friedman rank test were computed, as illustrated in the last lines of Table 2. Referring to the number of ranks where the MRFO outperforms the other peers (R+) and the reverse (R−), one can observe that the MRFO has a remarkable performance. The p-values based on the Wilcoxon sign rank are less than 0.05, therefore, the null hypothesis is rejected (h = No), meaning there exists a significant difference between optimizers in favor of MRFO (R+ > R−). The average rank based on the Friedman rank test is reported at the last time to order the algorithms-as it is obvious the MRFO has the least mean rank, therefore, it has the first rank as a final classification. The GWO occupies the second position as a final rank with a Friedman average rank of 4.6833, and subsequently, MRFO is considered the most reliable and consistent optimizer. For further evaluation of the MRFO performance, Figure 6 shows the cost function variation using different algorithms.
classification. The GWO occupies the second position as a final rank with a Friedman average rank of 4.6833, and subsequently, MRFO is considered the most reliable and consistent optimizer. For further evaluation of the MRFO performance, Figure 6 shows the cost function variation using different algorithms. On the other hand, an analysis of the variance test (ANOVA) has been done to verify the consistency of the MRFO in determining the parameters of OFMPPT and the design problematic. ANOVA test results are presented in Table 3, and distribution of the observations through the number of runs is presented in Figure 7. From these results, the pvalue is much higher than the F value, which confirms the difference between the provided performances. As illustrated in Figure 7, the MRFO can provide the best results in terms of mean fitness. Moreover, it has the lowest variation range, where the fitness function is varied from 4.75 to 4.95, on the other hand, the SOA has outliers over a wide range from 2.2 to 4.89. Accordingly, the results demonstrate the MRFO robustness and consistency. To confirm the results obtained by the ANOVA test, and to approve the high performance of the proposed algorithm, a Tukey's HSD (honestly significant difference) test, which is used to confirm significant differences between the results means for each optimizer, has been performed. The results are presented in Figure 8, and, similar to the findings of the ANOVA test, the MRFO provides the higher mean fitness, thus, it has the ability to solve this kind of problem more efficiently. On the other hand, an analysis of the variance test (ANOVA) has been done to verify the consistency of the MRFO in determining the parameters of OFMPPT and the design problematic. ANOVA test results are presented in Table 3, and distribution of the observations through the number of runs is presented in Figure 7. From these results, the p-value is much higher than the F value, which confirms the difference between the provided performances. As illustrated in Figure 7, the MRFO can provide the best results in terms of mean fitness. Moreover, it has the lowest variation range, where the fitness function is varied from 4.75 to 4.95, on the other hand, the SOA has outliers over a wide range from 2.2 to 4.89. Accordingly, the results demonstrate the MRFO robustness and consistency.   To confirm the results obtained by the ANOVA test, and to approve the high performance of the proposed algorithm, a Tukey's HSD (honestly significant difference) test, which is used to confirm significant differences between the results means for each optimizer, has been performed. The results are presented in Figure 8, and, similar to the findings of the ANOVA test, the MRFO provides the higher mean fitness, thus, it has the ability to solve this kind of problem more efficiently.  To validate the superiority of the MRFO-based OFMPPTS, the tracking performance is compared with incremental conductance and hill climbing approaches. For the incremental conductance, gain of the discrete integrator is assumed to be 0.7. But for OFMP-PTS, parameters of fractional order PI are 00.03526, 7.03875, and 0.98753, respectively, for proportional gain, integral gain, and fractional order as determined using MRFO, and provided in Table 1. A layout of the TEG system is presented in Figure 9. The TEG system consists of five series connected TEG units, a DC converter working in the continuous current mode with 30 kHz as the switching frequency, an input inductance with 1 mH, and an output capacitor with 47 μF. The load value is 15 Ω. The location of MPP is mainly dependent on the difference in temperature of the two sides of the TEG, in addition to the To validate the superiority of the MRFO-based OFMPPTS, the tracking performance is compared with incremental conductance and hill climbing approaches. For the incremental conductance, gain of the discrete integrator is assumed to be 0.7. But for OFMPPTS, parameters of fractional order PI are 00.03526, 7.03875, and 0.98753, respectively, for proportional gain, integral gain, and fractional order as determined using MRFO, and provided in Table 1. A layout of the TEG system is presented in Figure 9. The TEG system consists of five series connected TEG units, a DC converter working in the continuous current mode with 30 kHz as the switching frequency, an input inductance with 1 mH, and an output capacitor with 47 µF. The load value is 15 Ω. The location of MPP is mainly dependent on the difference in temperature of the two sides of the TEG, in addition to the load demanded, so five different operating scenarios are considered to test the performance of OFMPPTS. These scenarios include varying the temperature for the hot side, varying the temperature for the cold side, and varying the load value. The details parameters for the investigated scenarios are presented in Table 4.
Mathematics 2021, 9, x FOR PEER REVIEW 12 of 1 load demanded, so five different operating scenarios are considered to test the perfor mance of OFMPPTS. These scenarios include varying the temperature for the hot side varying the temperature for the cold side, and varying the load value. The details param eters for the investigated scenarios are presented in Table 4. Figure 9. A complete TEG system diagram.

1st Scenario
In this scenario, as presented in Table 4, two different values of the difference tem perature (ΔT) between the hot and cold sides of the TEG are used: 270 °C and 220 °C. Th load resistance and cold side temperature are 15 Ω and 30 °C, respectively, for both cases The hot side temperature values are 300 °C and 250 °C, respectively. Under these cond tions, the power, voltage, and current from the TEG at MPP are 73 W, 21 V, and 3.4 A respectively, for a difference temperature of 270 °C, whereas the power, voltage, and cur rent are: 57.86 W; 18.36 V; and 3.14 A, respectively, for a difference temperature of 220 °C The tracking system must preserve the maximum power regardless of the load value. Th tracking performance using OFMPPTS, HC, and INC methods is presented in Figure 10 for a difference temperature of 270 °C, whereas in Figure 10b, for a difference temperatur of 220 °C. At the beginning, for both cases, the suggested OFMPPTS achieves a fast re sponse to achieve the maximum power. It needs around 0.02 s to catch the MPP, wherea the INC requires more than 0.2 s. As such, thanks to the proposed OFMPPTS, the dynami tracking performance has been improved. Another superior criterion for the suggeste OFMPPTS represents the low fluctuations under the steady-state condition compare Figure 9. A complete TEG system diagram. Table 4. The details parameters for the investigated scenarios.

Scenario Time
Step

1st Scenario
In this scenario, as presented in Table 4, two different values of the difference temperature (∆T) between the hot and cold sides of the TEG are used: 270 • C and 220 • C. The load resistance and cold side temperature are 15 Ω and 30 • C, respectively, for both cases. The hot side temperature values are 300 • C and 250 • C, respectively. Under these conditions, the power, voltage, and current from the TEG at MPP are 73 W, 21 V, and 3.4 A, respectively, for a difference temperature of 270 • C, whereas the power, voltage, and current are: 57.86 W; 18.36 V; and 3.14 A, respectively, for a difference temperature of 220 • C. The tracking system must preserve the maximum power regardless of the load value. The tracking performance using OFMPPTS, HC, and INC methods is presented in Figure 10a for a difference temperature of 270 • C, whereas in Figure 10b, for a difference temperature of 220 • C. At the beginning, for both cases, the suggested OFMPPTS achieves a fast response to achieve the maximum power. It needs around 0.02 s to catch the MPP, whereas the INC requires more than 0.2 s. As such, thanks to the proposed OFMPPTS, the dynamic tracking performance has been improved. Another superior criterion for the suggested OFMPPTS represents the low fluctuations under the steady-state condition compared with the HC method. Therefore, the dynamic tracking performance of OFMPPTS is superior when compared with INC and HC approaches for the fast tracking, as well as the diminished fluctuation at steady-state conditions. with the HC method. Therefore, the dynamic tracking performance of OFMPPTS is superior when compared with INC and HC approaches for the fast tracking, as well as the diminished fluctuation at steady-state conditions.

2nd Scenario
In the second scenario, the cold side temperature and the value of the load are maintained at 30 • C and 15 Ω, respectively. On the other side, the hot side temperature is reduced from 300 • C to 250 • C at a time of 0.6 s. Under such a condition, TEG power is decreased from 73 W to 57.86 W, due to the decreasing of the temperance difference between both the two TEG sides from 270 • C to 220 • C. The dynamic tracking performance using OFMPPTS, HC, and INC methods for the second scenario is demonstrated in Figure 11. After the sudden change in the hot side temperature, the suggested OFMPPTS reaches the maximum power of 57.86 W rapidly, and with the lowest fluctuations compared with the INC and HC. The INC needs a large time to come back to the optimal operating point, and still the HC suffers from the high oscillations around the MPP under steady-state.
In the second scenario, the cold side temperature and the value of the load are maintained at 30 °C and 15 Ω, respectively. On the other side, the hot side temperature is reduced from 300 °C to 250 °C at a time of 0.6 s. Under such a condition, TEG power is decreased from 73 W to 57.86 W, due to the decreasing of the temperance difference between both the two TEG sides from 270 °C to 220 °C. The dynamic tracking performance using OFMPPTS, HC, and INC methods for the second scenario is demonstrated in Figure  11. After the sudden change in the hot side temperature, the suggested OFMPPTS reaches the maximum power of 57.86 W rapidly, and with the lowest fluctuations compared with the INC and HC. The INC needs a large time to come back to the optimal operating point, and still the HC suffers from the high oscillations around the MPP under steady-state. Figure 11. Dynamic tracking performance using OFMPPTS, HC, and INC methods for 2nd scenario.

3rd Scenario
In the third scenario, the temperature for the hot side and the value of the load are maintained at 300 °C and 15 Ω, respectively. The temperature for the cold side is raised from 30 °C to 50 °C at a time of 0.6 s. Under such condition, TEG power is reduced from 73 W to 60.57 W, due to the decreasing of the temperance difference between both two TEG sides from 270 °C to 250 °C. The dynamic tracking performance using OFMPPTS, HC, and INC methods for the third scenario is presented in Figure 12. The suggested OF-MPPTS is capable of reaching the maximum power (60.57 W) rapidly after 0.02 s, but the INC still needs more time for updating the duty cycle, and returns to the optimal operating point.

3rd Scenario
In the third scenario, the temperature for the hot side and the value of the load are maintained at 300 • C and 15 Ω, respectively. The temperature for the cold side is raised from 30 • C to 50 • C at a time of 0.6 s. Under such condition, TEG power is reduced from 73 W to 60.57 W, due to the decreasing of the temperance difference between both two TEG sides from 270 • C to 250 • C. The dynamic tracking performance using OFMPPTS, HC, and INC methods for the third scenario is presented in Figure 12. The suggested OFMPPTS is capable of reaching the maximum power (60.57 W) rapidly after 0.02 s, but the INC still needs more time for updating the duty cycle, and returns to the optimal operating point.

4th Scenario
In this scenario, there is a sudden change simultaneously for the hot side temperature and the cold side temperature, respectively, when keeping the load constant at 15 Ω. At 0.6 s, the cold side temperature is increased from 50 °C to 80 °C, and, at the same time, the

4th Scenario
In this scenario, there is a sudden change simultaneously for the hot side temperature and the cold side temperature, respectively, when keeping the load constant at 15 Ω. At 0.6 s, the cold side temperature is increased from 50 • C to 80 • C, and, at the same time, the hot side temperature is increased from 250 • C to 300 • C. With this scenario, TEG power is decreased from 47 W to 44.85 W. The dynamic tracking performance using OFMPPTS, HC, and INC methods for the fourth scenario is given in Figure 13. After considering the sudden change in the temperature, the suggested OFMPPTS rapidly succeeded to catch the maximum power of 44.85 W after 0.02 s, but the INC still needs around 0.2 s to reach to the same point. The high oscillations around the MPP are very clear when using the HC method. Such oscillations are removed thanks to using the OFMPPTS.

4th Scenario
In this scenario, there is a sudden change simultaneously for the hot side temperature and the cold side temperature, respectively, when keeping the load constant at 15 Ω. At 0.6 s, the cold side temperature is increased from 50 °C to 80 °C, and, at the same time, the hot side temperature is increased from 250 °C to 300 °C. With this scenario, TEG power is decreased from 47 W to 44.85 W. The dynamic tracking performance using OFMPPTS, HC, and INC methods for the fourth scenario is given in Figure 13. After considering the sudden change in the temperature, the suggested OFMPPTS rapidly succeeded to catch the maximum power of 44.85 W after 0.02 s, but the INC still needs around 0.2 s to reach to the same point. The high oscillations around the MPP are very clear when using the HC method. Such oscillations are removed thanks to using the OFMPPTS.

5th Scenario, Change in the Load
Evaluating the tracking performance is important when a sudden change occurs in the load demand. In the fifth scenario, the temperature values for both hot and cold sides are fixed to 300 • C and 80 • C, respectively. In this case, the nominal values for voltage current and power are 15.2 V, 2.93 A, and 44.85 W, respectively. At a time of 0.4 s, a sudden change in the load demand is implemented. An additional resistance of 25 Ω is connected in parallel with the refence load of 25 Ω. The dynamic tracking performance using OFMPPTS, HC, and INC methods for the fifth scenario is given in Figure 14. Following the sudden change in the load, the suggested OFMPPTS rapidly succeeded to catch the maximum power of 44.85 W after 0.01 s, but the INC still needs around 0.2 s to reach to the same point of 44.85 W. The high oscillations around the MPP are very clear when using the HC method. Such oscillations are removed thanks to using the OFMPPTS. change in the load demand is implemented. An additional resistance of 25 Ω is connected in parallel with the refence load of 25 Ω. The dynamic tracking performance using OFMP-PTS, HC, and INC methods for the fifth scenario is given in Figure 14. Following the sudden change in the load, the suggested OFMPPTS rapidly succeeded to catch the maximum power of 44.85 W after 0.01 s, but the INC still needs around 0.2 s to reach to the same point of 44.85 W. The high oscillations around the MPP are very clear when using the HC method. Such oscillations are removed thanks to using the OFMPPTS.

Conclusions
To improve the transient response of incremental conductance-based MPPT, and remove the steady-state fluctuations of the hill climbing approach, an optimized fractional MPPT strategy (OFMPPTS) has been developed to maximize the harvested energy from thermoelectric generator systems. The suggested OFMPPTS utilizes the fractional order control so that it provides more adjustable time and frequency responses of the control system, allowing better and more robust performance. The best parameters of the OFMP-PTS have been determined firstly using MRFO. Ten optimization algorithms: WOA; HHO; SMA; PSO; HBO; GBO; GA; GWO; SOA; and TSA were used in comparison with the proposed MRFO to identify the best parameters of OFMPPTS. The main finding proved that the performance of MRFO is the best among all considered algorithms. The average cost function values varied between 4.92934 kWh and 3.59264 kWh. The maximum average cost function of 4.92934 kWh has been achieved by MRFO, followed by SOA (4.5721 kWh). The STD values varied between 0.04867 and 1.18114. The lowest STD was also accomplished by MRFO. Regarding the efficiency, the maximum efficiency of 99.46% has been obtained by MRFO, whereas the lowest efficiency of 74.01% has been obtained by GA. Also, two tests were performed to demonstrate the robustness of MRFO: analysis of the variance test (ANOVA) and Tukey's (HSD). Finally, to validate the superiority of the MRFO-based OFMPPTS, the tracking performance was compared with incremental

Conclusions
To improve the transient response of incremental conductance-based MPPT, and remove the steady-state fluctuations of the hill climbing approach, an optimized fractional MPPT strategy (OFMPPTS) has been developed to maximize the harvested energy from thermoelectric generator systems. The suggested OFMPPTS utilizes the fractional order control so that it provides more adjustable time and frequency responses of the control system, allowing better and more robust performance. The best parameters of the OFMPPTS have been determined firstly using MRFO. Ten optimization algorithms: WOA; HHO; SMA; PSO; HBO; GBO; GA; GWO; SOA; and TSA were used in comparison with the proposed MRFO to identify the best parameters of OFMPPTS. The main finding proved that the performance of MRFO is the best among all considered algorithms. The average cost function values varied between 4.92934 kWh and 3.59264 kWh. The maximum average cost function of 4.92934 kWh has been achieved by MRFO, followed by SOA (4.5721 kWh). The STD values varied between 0.04867 and 1.18114. The lowest STD was also accomplished by MRFO. Regarding the efficiency, the maximum efficiency of 99.46% has been obtained by MRFO, whereas the lowest efficiency of 74.01% has been obtained by GA. Also, two tests were performed to demonstrate the robustness of MRFO: analysis of the variance test (ANOVA) and Tukey's (HSD). Finally, to validate the superiority of the MRFO-based OFMPPTS, the tracking performance was compared with incremental conductance and hill climbing approaches. Five different operating scenarios were considered to test the performance of OFMPPTS. These scenarios included varying the hot side temperature, varying the cold side temperature, and varying the load value. The results emphasized the superiority of the optimized fractional MPPTS compared with conventional methods for both steady-state and dynamic responses. Future studies should employ the MRFO algorithm to address the hardware validation of the proposed MPP algorithm. Furthermore, the authors should consider the newest meta-heuristics for red fox optimization algorithm, polar bear algorithm, and black widow optimization algorithm to solve the problem of partial shading of the PV array surface.