Optimal Parameter Identification of a PEM Fuel Cell Using Recent Optimization Algorithms

: The parameter identiﬁcation of a PEMFC is the process of using optimization algorithms to determine the ideal unknown variables suitable for the development of an accurate fuel-cell-performance prediction model. These parameters are not always available from the manufacturer’s datasheet, so they need to be determined to accurately model and predict the fuel cell’s performance. Five optimization methods—bald eagle search (BES) algorithm, equilibrium optimizer (EO), coot (COOT) algorithm, antlion optimizer (ALO), and heap-based optimizer (HBO)—are used to compute seven unknown parameters of a PEMFC. During optimization, these seven parameters are used as decision variables, and the ﬁtness function to be minimized is the sum square error (SSE) between the estimated cell voltage and the actual measured cell voltage. The SSE obtained for the BES algorithm was noted to be 0.035102. The COOT algorithm recorded an SSE of 0.04155, followed by ALO with an SSE of 0.04022 and HBO with an SSE of 0.056021. BES predicted the performance of the fuel cell accurately; hence, it is suitable for the development of a digital twin for fuel-cell applications and control systems for the automotive industry. Furthermore, it was deduced that the convergence speed for BES was faster compared to the other algorithms investigated. This study aims to use metaheuristic algorithms to predict fuel-cell performance for the development and commercialization of digital twins in the automotive industry.


Introduction
For both small power uses and bigger industrial applications, clean energy sources are becoming increasingly necessary due to the rapid decline in fossil fuel sources and the growing demand for electricity [1].Harnessing energy from renewable sources such as solar and wind is often relied on, but these sources are affected by environmental conditions, which has led to the development of fuel cells to help supplement existing green energy sources.Fuel cells have traditionally been divided into stationary, portable, or transportation types [2].Fuel-cell technology has advanced quickly in the automotive sector due to the growing usage of fuel cells for heavy-duty land vehicles such as public buses.Stationary fuel cells for homes and businesses have also become more popular [3].
Energies 2023, 16, 5246 2 of 20 Stationary fuel cells have a variety of applications.Several businesses and researchers have developed a keen interest in fuel cells in recent years.The chemical energy produced by the reaction between oxygen and hydrogen, or natural air, can be quickly converted into electrical energy by fuel cells [4].Many fuel-cell types exist, such as alkaline, solid oxide, proton exchange membrane, and phosphoric acid fuel cells [5].Each of these types of fuel cells has its own use, but the most common type of fuel cell used in the automotive industry is the PEMFC [6,7].The advantages of PEMFCs include their high efficiency, low emissions, high power density, and the ability to be quickly refueled.PEMFCs can be used to power many types of applications, including transportation, standby power, residential, and industrial applications.Additionally, PEMFC technology is relatively simple and cost-effective when compared to other fuel-cell technologies.The No Free Lunch Theorem states that any optimization algorithm will have no performance advantage over any other algorithm when averaged over all possible problems.This follows from the fact that any algorithm must make assumptions and trade-offs to optimize a problem.Any such assumptions can be beneficial for some problems while detrimental for others.The accuracy of the mathematical model for a PEMFC is highly dependent on the amount of manufacturer information available.Therefore, the unknown parameters for PEMFCs must be clearly defined to establish a perfect agreement between the experimental and mathematical models.The chemical characteristic within the cell is a key determinant of the output voltage of the cell [8].There are currently some concepts established in the literature capable of estimating the exact PEMFC parameters.Most of these optimization-based approaches are considered as being simple, consistent, and robust [9].In terms of the optimal parameters in relation to PEMFCs, several investigations have been carried out in the literature.For instance, a multiverse optimizer was explored for the determination of the parameters for a PEMFC-equivalent circuit [10].The total number of parameters considered in the study was seven.Other authors equally reported various optimization parameters capable of ensuring convergence [11].There was a clear classification of the various methods, specifically involving evolutionary-based, swarm physics as well as nature.Hegazy et al. [12] also evaluated various metaheuristics in the optimization of microgrids.To reduce the sum squared error (SSE), adaptive sparrow search algorithms were also explored to evaluate the variation between the calculated and measured voltage [13].Yousri et al. [14] concluded that the fractional-order-modified Harris hawk optimizer was the best algorithm for the mathematical modeling of PEMFCs.To further reduce the disparity between the estimated and empirical results, other authors used a developed coyote optimization algorithm [15].Yuan et al. [15] carried out the investigation using two types of fuel cells.The application of an improved monarch butterfly optimizer has also been explored in determining unknown parameters to ensure a reduction in the integral time absolute error [16].With varying conditions around the cell, two types of cells were also studied (2 kW Nexa and 6 kW Nedstack PS6) using an improved chimp optimizer [17].The improved chimp optimizer was also evaluated for three commercial fuel cells.A 15-nature algorithm was used to validate the outcome of the study.Tabbi et al. [18] concluded that the artificial ecosystem-based algorithm presented a better result compared to the grey wolf optimizer, particle swarm optimizer, slime mold algorithm, and Harris hawk optimizer [18].A hybrid method comprising a combination between a vortex search algorithm and differential evolution for the identification of the optimal parameters for PEMFCs has also been presented in the literature [19].The SSE was considered as the ideal fitness function between the experimental voltages and those generated for the stack deduced mathematically.A monarch butterfly optimizer was also experimented on a 250 W PEMFC stack under varying conditions [20].A further investigation also evaluated the ideal parameters for PEMFCs by considering the SSE between experimental and numerical data using metaheuristics [21].A slime mold optimizer was also used by Gupta et al. for the determination of the unknown parameters for PEMFCs [22].An improved search optimizer for evaluating the parameters for PEMFCs was also investigated by Qin et al. [23].The SSE between the experimental and estimated data was considered as the target.A Bayesian-regularized neural network was adopted in extracting the ideal parameters for PEMFCs [24].Using a sunflower optimizer, an estimation of the unknown parameters for PEMFCs was carried out by Yuan et al. [25].A converged moth search algorithm has also been reported in the literature to be suitable in reducing the SSE between the measured and experimental voltages [26].BCS 500 stack and Nedstack PS6 data were used for the study.A bi-subgroup algorithm for extracting the unknown parameters from PEMFCs has also been reported in the literature [27].Syah et al. [28] used a balanced strider method in reducing the total squared variation between experimental and numerical voltages.The cost of the stack was also utilized as the target using a modified grass fibrous root algorithm [29].Using three types of fuel cells (BCS 500/250 and Nedstack PS6), a particle swarm optimizer was utilized in the identification of the best parameters for the fuel cell [30].Mossa et al. [31] used two optimizers (Harris hawk and atom search algorithm) mainly for estimating the unknown values for PEMFCs.The SSE value recorded for the measured and mathematical data was taken as the fitness function.Rezk et al. [32] utilized a gradient-based algorithm for the determination of the ideal parameters between three types of fuel cells and compared it with other algorithms.A chaos owl search optimizer has also been adopted to reduce the sum squared deviation for measured as well as mathematically determined voltages [33].A chaotic binary shark smell optimizer was also adopted for the estimation of the unknown data for a PEMFC [34].A heterogeneous comprehensive learning Archimedes optimizer was also adopted for determining the unknown parameters for fuel cells [35].A deer hunting algorithm has also been adopted in estimating the parameters for a PEMFC.The study also involved the application of artificial neural networks [36].A water cycle algorithm was equally harnessed in estimating the unknown values for PEMFCs [37].An enhanced bald eagle search algorithm has also been reported in the literature as being ideal for exploring the unknown values [38] using the SSE as the target for reducing the values between the measured and estimated data.A coyote optimizer was employed by Abaza et al. [39] for estimating the unknown parameters.The utilization of a semi-empirical model for PEMFCs has equally been reported using the SSE as the target [40].By reducing the integral absolute error between the various types of voltages, Lu et al. [41] explored the application of a crow search approach in the determination of unknown parameters.An equilibrium optimizer has equally been adopted for estimating the unknown values for fuel cells [42].A satin bowerbird algorithm has been captured in the literature as being suitable for the estimation of the ideal parameters for developing polarization curves [43].An L-SHADE-EpSin approach was used in the literature in the development of a model for PEMFC estimation [44].Isa et al. [45] equally explored the application of an antlion as well as a dragonfly optimizer for estimating the unknown parameters for PEMFCs.The identification of fuel PEMFC parameters has also been explored using Harris hawk algorithms by Song et al. [46].A barnacles mating optimizer was also explored in modeling the fuel-cell behavior accurately [47].An Archimedes optimizer was also used in the literature mainly for reducing the deviation between the data gathered from the experimental and mathematical models [48].A transient search optimizer was also reported as being suitable in accurately estimating the unknown parameters for PEMFCs [49].Calasan et al. [50] used the Lambert W function in reducing the deviation between measured as well as simulated voltages.This study, however, presents the accurate prediction of the seven unknown parameters of PEMFCs often not available in the manufacturer's data sheet.The study further corroborates the prospects of metaheuristic algorithms in the development of a digital twin for PEMFCs, paying critical attention to the sum squared error, which is the objective function, and the computational time for the entire predictive process.Section 2 of the investigation captures the detailed information regarding the electrochemical reaction of a PEMFC, while Section 3 highlights the information on the PEMFC-parameter-estimation approach and Section 4 delves more into the bald eagle search (BES) algorithm.A summary of the experimental procedure is found in Section 5, and the results are summarized and discussed in Section 6.The main findings are outlined in Section 7 to reflect on the overall outcome of the study.

Electrochemical Reaction Inside a PEMFC
The electrochemical reaction inside a PEMFC involves the oxidation of fuel molecules at the anode and the reduction in oxygen molecules at the cathode.Protons (H+) are generated at the anode and flow across an electrolyte membrane to the cathode [51].This electrochemical reaction produces water, heat, and electricity.Figure 1 presents the various compositions of a fuel cell.As can be observed, the various electrodes (anode and cathode) are separated by a membrane.The membrane is designed to be permeable to only protons.Furthermore, there is the presence of a catalyst to speed up the electrochemical reaction [52].The cell is designed to ensure that when the hydrogen fuel reaches the membrane, it dissociates into two ions.Due to the morphological composition of the membrane, the protons go through the membrane to the cathodic electrode, while the electrons go through an external circuit to produce electricity and water as a byproduct of the reaction.The overall reaction within the cell is summarized in Equations ( 1)-( 3) [51].

Electrochemical Reaction Inside a PEMFC
The electrochemical reaction inside a PEMFC involves the oxidation of fuel molec at the anode and the reduction in oxygen molecules at the cathode.Protons (H+) are erated at the anode and flow across an electrolyte membrane to the cathode [51].electrochemical reaction produces water, heat, and electricity.Figure 1 presents the ous compositions of a fuel cell.As can be observed, the various electrodes (anode cathode) are separated by a membrane.The membrane is designed to be permeab only protons.Furthermore, there is the presence of a catalyst to speed up the electroch ical reaction [52].The cell is designed to ensure that when the hydrogen fuel reache membrane, it dissociates into two ions.Due to the morphological composition of the m brane, the protons go through the membrane to the cathodic electrode, while the elect go through an external circuit to produce electricity and water as a byproduct of the tion.The overall reaction within the cell is summarized in Equations ( 1)-( 3) [51].Anode Complete chemical reaction

Fuel Cell Modeling Mathematically
For both small power uses and bigger industrial applications, clean energy sou are becoming increasingly necessary due to the rapid decline in fossil fuel sources an growing demand for electricity [50].Harnessing energy from renewables is often r on, but these sources are affected by environmental conditions, which has led to th velopment of fuel cells to help supplement existing green energy sources.Fuel cells traditionally been divided into stationary, portable, or transportation types [51,52].Anode Complete chemical reaction

Fuel Cell Modeling Mathematically
For both small power uses and bigger industrial applications, clean energy sources are becoming increasingly necessary due to the rapid decline in fossil fuel sources and the growing demand for electricity [50].Harnessing energy from renewables is often relied on, but these sources are affected by environmental conditions, which has led to the development of fuel cells to help supplement existing green energy sources.Fuel cells have traditionally been divided into stationary, portable, or transportation types [51,52].The Energies 2023, 16, 5246 5 of 20 polarization curve for a fuel cell being operated at 80 • C is captured in Figure 2. It can be noticed that there are three main areas on the polarization curve.These areas are broadly known as activation losses, ohmic losses, and concentration losses [53].There is no linearity within the activation region.The activation region presents holistic information regarding the electrochemical reaction within the cell.The ohmic losses are commonly found within the membrane.The final section is the mass concentration losses due to variations in the concentration gradient within the cell [54].The total cell voltage is depicted in Equation ( 4) as V f c [53].
known as activation losses, ohmic losses, and concentration losses [53].There ity within the activation region.The activation region presents holistic informa ing the electrochemical reaction within the cell.The ohmic losses are comm within the membrane.The final section is the mass concentration losses due t in the concentration gradient within the cell [54].The total cell voltage is depict tion (4) as   [53].  denotes the open circuit voltage, while   represents the activatio tion;  ℎ is the ohmic loss, and   is the concentration loss [54].It is a that within the ohmic section, the output voltage is dependent on the current d slope is also obtained based on the ionic resistance of the electrolyte, as explai The concentration loss is due to mass transfer limitations leading to a sharp de voltage to zero.Increasing the total output voltage(  ) of the cell is dependent o ber of cells (  ) connected in a series, as depicted in Equation ( 5) [50].
=   ×     is basically the open circuit voltage, which is determined using the N tion [55], but other parameters that consider the variation in temperature surro cell are also taken into account, as depicted in Equation ( 6) [50].The temperature for the cell is represented as , while the partial pressur is   2 and that of hydrogen is   2 .Equations ( 7) and ( 8) denote the various p E cell denotes the open circuit voltage, while V act represents the activation polarization; V ohmic is the ohmic loss, and V conc is the concentration loss [54].It is also evident that within the ohmic section, the output voltage is dependent on the current density.The slope is also obtained based on the ionic resistance of the electrolyte, as explained earlier.The concentration loss is due to mass transfer limitations leading to a sharp decline in the voltage to zero.Increasing the total output voltage (V t ) of the cell is dependent on the number of cells (X n ) connected in a series, as depicted in Equation ( 5) [50].
E cell is basically the open circuit voltage, which is determined using the Nernst equation [55], but other parameters that consider the variation in temperature surrounding the cell are also taken into account, as depicted in Equation ( 6) [50].
Energies 2023, 16, 5246 6 of 20 The temperature for the cell is represented as T, while the partial pressure of oxygen is P O 2 and that of hydrogen is P H 2 .Equations ( 7) and ( 8) denote the various partial pressure parameters represented mathematically [51,52].
The anodic relative humidity is represented as RH a , while that of the cathode is RH c .The anode pressure at the inlet is P a , while that at the cathode is P c .The area of the cell is captured as A, while the current is i cell .The water vapor saturation is (P Sat H 2 O ), and this parameter has a direct correlation to temperature T, as captured in Equation ( 9).On the other hand, Equation ( 10) is utilized in the determination of the activation losses.ξ 1, ξ 2, ξ 3, ξ 4, are semi-empirical parametric coefficients, while oxygen concentration is highlighted as C O 2 and computed using Equation (11).Calculating the ohmic losses is achieved using Equation ( 12) [52].
C O 2 = P O 2 5.08 × 10 6 exp 498 T (11) R m and R c Ωcm −2 denote the electronic and ionic resistance, respectively.The electronic resistance is attributed to the slightest perturbations in relation to the current and voltage, and Equation ( 13) [50] is utilized to calculate this mathematically, while the membrane parametric coefficient is determined using Equation (14).Equation ( 15) [51] is utilized in computing the concentration polarization mathematically.b is the parametric coefficient, sometimes referred to as the diffusion parameter, and the maximum current density is J max , while J is the actual current density.

PEMFC-Parameter-Estimation Approach
In developing a computational model for PEMFCs mathematically, the calculation of the specific seven model parameters ( ξ 1, ξ 2, ξ 3, ξ 4, R, B, λ) is very important.The major setback here is the fact that these parameters are usually unknown.There is an equally significant variation in the model parameters based on the operating conditions.These phenomena often cause an effect on the developed IV curve being accurate.Most manufacturers of fuel cells do not usually provide this information, and identifying these parameters is challenging.A solution to mitigate this challenge is considering the problem using various optimization techniques.The introduction of an artificial intelligence approach for determining the unknown parameters for fuel cells is gradually becoming a primary research focus.The estimation of model parameters can be deduced from the experimental data using the RMSE as an objective function for the experimental and numerically determined datasets.From Equation ( 16), the experimental value is V a , while the predicted voltage is V i .N is the number of data points.Figure 3 highlights the various steps adopted in the parametric estimation of a fuel cell's unknown parameters.
predicted voltage is   .N is the number of data points.Figure 3 highlights the various steps adopted in the parametric estimation of a fuel cell's unknown parameters.

Bald Eagle Search Algorithm
The BES algorithm is an efficient heuristic search algorithm that uses a combination of depth-first and best-first searches.It works by first expanding a node according to the best-first search, then exploring each branch in the deepest possible manner.This algorithm has been found to be effective for a range of applications, such as finding paths in graphs and solving mazes.The algorithm involves three stages: selecting the space with the most potential prey, searching within that space, and swooping from the best-found position to an optimal hunting spot.The first stage is modeled as follows: where α is a constant [1.5, 2], and r is a random value.The second stage can be modeled as follows: ) where X and Y are directional coordinates, calculated as follows:

Bald Eagle Search Algorithm
The BES algorithm is an efficient heuristic search algorithm that uses a combination of depth-first and best-first searches.It works by first expanding a node according to the best-first search, then exploring each branch in the deepest possible manner.This algorithm has been found to be effective for a range of applications, such as finding paths in graphs and solving mazes.The algorithm involves three stages: selecting the space with the most potential prey, searching within that space, and swooping from the best-found position to an optimal hunting spot.The first stage is modeled as follows: where α is a constant [1.5, 2], and r is a random value.The second stage can be modeled as follows: where X and Y are directional coordinates, calculated as follows: where β 1 is a constant [5,10], and R is a constant gain [0.5, 2].The last stage can be expressed as follows: Energies 2023, 16, 5246 The BES flowchart is illustrated in Figure 4.
The BES flowchart is illustrated in Figure 4.

Experimental Procedure
Avista SR-12 PEMFC is considered in this work.The current deduced from the cell varied between 0 and 34 A, and the variation in terms of voltage was between 23 and 43 VDC.The fuel for the cell was passed through an in-house-developed chamber to ensure the hydrogen gas was humidified before entering the cell.This step was critical in ensuring the membrane was well-humidified to allow an increase in protonic conductivity but a reduction in resistance to the electrolyte.Attached to the cell is also a boost converter and a battery, as well as a load cell.The hydrogen gas entering the cell was varied in terms of pressure and flow rate (Figure 5a) to ensure the determination of the effect of the operating conditions on the overall cell performance.The cell was made up of 48 cells, with the active area denoted as 62.5 cm 2 .The highest cell current was 42 A, while the stack temperature varied between 65 and 80 °C.The cell stacks were cooled using air (Figure 5a,b).

Experimental Procedure
Avista SR-12 PEMFC is considered in this work.The current deduced from the cell varied between 0 and 34 A, and the variation in terms of voltage was between 23 and 43 VDC.The fuel for the cell was passed through an in-house-developed chamber to ensure the hydrogen gas was humidified before entering the cell.This step was critical in ensuring the membrane was well-humidified to allow an increase in protonic conductivity but a reduction in resistance to the electrolyte.Attached to the cell is also a boost converter and a battery, as well as a load cell.The hydrogen gas entering the cell was varied in terms of pressure and flow rate (Figure 5a) to ensure the determination of the effect of the operating conditions on the overall cell performance.The cell was made up of 48 cells, with the active area denoted as 62.5 cm 2 .The highest cell current was 42 A, while the stack temperature varied between 65 and 80 • C. The cell stacks were cooled using air (Figure 5a,b).

Results and Discussion
The optimal parameters of an SR-12 PEM FC have been defined by using five recent optimization methods, namely, bald eagle search algorithm (BES), equilibrium optimizer (EO), coot algorithm (COOT), antlion optimizer (ALO), and heap-based optimizer (HBO).The laptop specifications were HP OMEN 17, Core i9, 32 GB, 1 TB SSD.The simulations were performed using MATLAB software version 2020a.The specifications of SR-12 PEM 500 are shown in Table 1.

Results and Discussion
The optimal parameters of an SR-12 PEM FC have been defined by using five recent optimization methods, namely, bald eagle search algorithm (BES), equilibrium optimizer (EO), coot algorithm (COOT), antlion optimizer (ALO), and heap-based optimizer (HBO).The laptop specifications were HP OMEN 17, Core i9, 32 GB, 1 TB SSD.The simulations were performed using MATLAB software version 2020a.The specifications of SR-12 PEM 500 are shown in Table 1.In order to guarantee an equal comparison, the number of populations was kept at 25, while the maximum number of iterations (n max ) was selected as 250 and 500 to show the effect of the number of iterations on the algorithms' performance.The maximum and minimum limits of the unknown parameters are outlined in Table 2. Table 3 shows the optimal values of the PEMFC parameters after 30 runs.The absolute error for the obtained parameters using different methods compared to BES is presented in Table A1.The results were evaluated statistically, as shown in Table 4.   Considering Table 4, with a maximum number of 500 iterations, the mean cost function values range from 0.035102 to 0.056021.The minimum mean cost function value of 0.035102 is obtained by BES, followed by COOT (0.04155), whereas the largest mean cost function value of 0.056021 is obtained by HBO.The standard deviation values range between 1.15 × 10 −5 and 0.019082.The minimum standard deviation value of 1.15 × 10 −5 is achieved by BES, followed by COOT (0.006312), whereas the largest standard deviation value of 0.019082 is obtained by HBO.In the case of a maximum number of 250 iterations, the mean cost function values range from 0.035794 to 0.089071.The minimum mean cost function value of 0.035794 is obtained by BES, followed by ALO (0.04758), whereas the largest mean cost function value of 0.089071 is obtained by HBO.The standard deviation values range between 0.001557 and 0.045321.The minimum standard deviation value of 0.001557 is achieved by BES, followed by ALO (0.010739), whereas the largest standard deviation value of 0.045321 is obtained by HBO.So, BES is better than ALO, EO, COOT, and HBO.The details of the objective function values for the different algorithms throughout 30 runs is demonstrated in Table 5. Figure 6 shows the current density, voltage, and power characteristics of the PEMFC using BES.There is an excellent agreement between the estimated and measured data.This proves the superiority of the BES in determining the unknown parameters of the PEMFC.The absolute error in the cell voltage with different optimization algorithms is presented in Table 6 and Figure 7.For the first run with 500 iterations, as presented in Tables 5 and 6, the SSE values range from 0.0351 to 0.06104.The minimum SSE value of 0.0351 is obtained using the BES, followed by COOT (0.0419), whereas the worst SSE of 0.06104 is obtained by ALO.The best MAE of 0.03251 is obtained by BES, and the worst MAE of 0.04247 is obtained by ALO.Referring to Figure 7a, the maximum absolute error of 0.13831 is obtained by ALO.In the case of the first run with 250 iterations, as demonstrated  5 and 6, the SSE values range from 0.03559 to 0.06954.The minimum SSE value of 0.03559 is obtained using the BES, followed by EO (0.03907), whereas the worst SSE of 0.06104 is obtained by ALO.The best MAE of 0.03251 is obtained by BES, and the worst MAE of 0.06954 is obtained by ALO.Referring to Figure 7b, the maximum absolute error of 0.24777 is obtained by HBO.
presented in Table 6 and Figure 7.For the first run with 500 iterations, as presented in Tables 5 and 6, the SSE values range from 0.0351 to 0.06104.The minimum SSE value of 0.0351 is obtained using the BES, followed by COOT (0.0419), whereas the worst SSE of 0.06104 is obtained by ALO.The best MAE of 0.03251 is obtained by BES, and the worst MAE of 0.04247 is obtained by ALO.Referring to Figure 7a, the maximum absolute error of 0.13831 is obtained by ALO.In the case of the first run with 250 iterations, as demonstrated in Tables 5 and 6, the SSE values range from 0.03559 to 0.06954.The minimum SSE value of 0.03559 is obtained using the BES, followed by EO (0.03907), whereas the worst SSE of 0.06104 is obtained by ALO.The best MAE of 0.03251 is obtained by BES, and the worst MAE of 0.06954 is obtained by ALO.Referring to Figure 7b, the maximum absolute error of 0.24777 is obtained by HBO.Figures 8 and 9 demonstrate the objective function variation with 500 and 250 iterations, respectively.As demonstrated in Figure 8a and Table 5, for the first run with 500 iterations, the SSE values converge to 0.061045, 0.035099, 0.041897, 0.043371, and 0.061045, respectively, for ALO, BES, COOT, EO, and HBO.BES catches the optimal solution rapidly, whereas HBO needs more time to reach its best solution.During Run 14, as demonstrated in Figure 8b and Table 5, the SSE values converge to 0.069544, 0.03559, 0.066005, 0.039074, and 0.069544, respectively, for ALO, BES, COOT, EO, and HBO.BES catches the optimal solution rapidly, whereas for the second time, HBO needs more time to reach its best solution.Figures 8 and 9 demonstrate the objective function variation with 500 and 250 iterations, respectively.As demonstrated in Figure 8a and Table 5, for the first run with 500 iterations, the SSE values converge to 0.061045, 0.035099, 0.041897, 0.043371, and 0.061045, respectively, for ALO, BES, COOT, EO, and HBO.BES catches the optimal solution rapidly, whereas HBO needs more time to reach its best solution.During Run 14, as demonstrated in Figure 8b and Table 5, the SSE values converge to 0.069544, 0.03559, 0.066005, 0.039074, and 0.069544, respectively, for ALO, BES, COOT, EO, and HBO.BES catches the optimal solution rapidly, whereas for the second time, HBO needs more time to reach its best solution.ANOVA and Tukey tests were conducted to support the investigation.An ANOVA test is a statistical method used to compare two or more groups.This test measures the mean differences between groups, providing an assessment of how likely each observed difference is to have arisen by random chance.An ANOVA helps answer the question of   ANOVA and Tukey tests were conducted to support the investigation.An ANOVA test is a statistical method used to compare two or more groups.This test measures the mean differences between groups, providing an assessment of how likely each observed difference is to have arisen by random chance.An ANOVA helps answer the question of ANOVA and Tukey tests were conducted to support the investigation.An ANOVA test is a statistical method used to compare two or more groups.This test measures the mean differences between groups, providing an assessment of how likely each observed difference is to have arisen by random chance.An ANOVA helps answer the question of Energies 2023, 16, 5246 15 of 20 whether the difference between the means of two or more samples is statistically significant.It works by comparing the variability between groups and the variability within groups to determine whether the differences between group means could have arisen by chance.Table 7 gives the ANOVA test results, and Figure 10 illustrates the corresponding ranking.The BES algorithm can deliver the best performance regarding mean fitness and variations, as demonstrated by Figure 10.The Tukey test is a statistical procedure used to compare the means of two or more independent samples or datasets.It is a type of post hoc test that takes into account the similarity of group means when evaluating the differences between them.The test is based on the assumption that all datasets have the same variance and that their distribution can be approximated by a normal distribution.It is a powerful tool for detecting significant differences between groups.
The Tukey HSD test, as shown in Figure 11, was performed to approve the two types of ANOVA test results.The ALO, EO, and HBO groups significantly differ from BES for the first case.The COOT mean results indicate that the COOT can provide the second-best results.For the second case, two groups (EO and HBO) have means significantly different from BES, whereas ALO and COOT can provide near performance.However, the BES algorithm offers the best performance for the two cases.The ultimate results of the BES confirm that it can effectively solve this issue.
whether the difference between the means of two or more samples is statistically significant.It works by comparing the variability between groups and the variability within groups to determine whether the differences between group means could have arisen by chance.Table 7 gives the ANOVA test results, and Figure 10 illustrates the corresponding ranking.The BES algorithm can deliver the best performance regarding mean fitness and variations, as demonstrated by Figure 10.The Tukey test is a statistical procedure used to compare the means of two or more independent samples or datasets.It is a type of post hoc test that takes into account the similarity of group means when evaluating the differences between them.The test is based on the assumption that all datasets have the same variance and that their distribution can be approximated by a normal distribution.It is a powerful tool for detecting significant differences between groups.
The Tukey HSD test, as shown in Figure 11, was performed to approve the two types of ANOVA test results.The ALO, EO, and HBO groups significantly differ from BES for the first case.The COOT mean results indicate that the COOT can provide the second-best results.For the second case, two groups (EO and HBO) have means significantly different from BES, whereas ALO and COOT can provide near performance.However, the BES algorithm offers the best performance for the two cases.The ultimate results of the BES confirm that it can effectively solve this issue.

Conclusions
The optimal parameter identification process of the SR-12 PEM fuel-cell model has been investigated in this research work using different recent optimization algorithms.Five optimization methods, namely, bald eagle search algorithm (BES), equilibrium optimizer (EO), coot algorithm (COOT), antlion optimizer (ALO), and heap-based optimizer (HBO), have been considered.Two different numbers of iterations were chosen: 250 and

Figure 2 .
Figure 2. Various types of losses in a fuel cell.

Figure 2 .
Figure 2. Various types of losses in a fuel cell.

Figure 3 .
Figure 3.A diagrammatic representation of the various approaches used in estimating a cell's unknown parameters.

Figure 3 .
Figure 3.A diagrammatic representation of the various approaches used in estimating a cell's unknown parameters.

Figure 6 .
Figure 6.Current density, voltage, and power characteristics of the PEMFC using BES.

Figure 6 .
Figure 6.Current density, voltage, and power characteristics of the PEMFC using BES.

Figure 7 .
Figure 7. Absolute error in the cell voltage with different optimization algorithms: (a) n max 500; (b) n max 250.

Table 2 .
Maximum and minimum limits of the parameters.

Table 3 .
Optimal values of the PEMFC model parameters using different optimization methods.
t is the average time for a single simulation run.

Table 4 .
Statistical assessments for the considered optimization algorithms.

Table 5 .
Objective function values for the different algorithms throughout 30 runs.

Table 6 .
Absolute error in the cell voltage with different optimization algorithms.

Table 6 .
Absolute error in the cell voltage with different optimization algorithms.