[go: up one dir, main page]

Next Article in Journal
Deep Learning-Based Vision Systems for Robot Semantic Navigation: An Experimental Study
Next Article in Special Issue
Granular Weighted Fuzzy Approach Applied to Short-Term Load Demand Forecasting
Previous Article in Journal
PARKTag: An AI–Blockchain Integrated Solution for an Efficient, Trusted, and Scalable Parking Management System
Previous Article in Special Issue
Ensemble Learning for Nuclear Power Generation Forecasting Based on Deep Neural Networks and Support Vector Regression
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Red-Billed Blue Magpie Optimizer for Electrical Characterization of Fuel Cells with Prioritizing Estimated Parameters

by
Attia A. El-Fergany
1 and
Ahmed M. Agwa
2,*
1
Electrical Power and Machines Department, Faculty of Engineering, Zagazig University, Zagazig 44519, Egypt
2
Department of Electrical Engineering, College of Engineering, Northern Border University, Arar 73222, Saudi Arabia
*
Author to whom correspondence should be addressed.
Technologies 2024, 12(9), 156; https://doi.org/10.3390/technologies12090156
Submission received: 8 August 2024 / Revised: 31 August 2024 / Accepted: 3 September 2024 / Published: 8 September 2024
(This article belongs to the Collection Electrical Technologies)

Abstract

:
The red-billed blue magpie optimizer (RBMO) is employed in this research study to address parameter extraction in polymer exchange membrane fuel cells (PEMFCs), along with three recently implemented optimizers. The sum of squared deviations (SSD) between the simulated and measured stack voltages defines the fitness function of the optimization problem under investigation subject to a set of working constraints. Three distinct PEMFCs stacks models—the Ballard Mark, Temasek 1 kW, and Horizon H-12 units—are used to illustrate the applied RBMO’s feasibility in solving this challenge in comparison to other recent algorithms. The highest percentages of biased voltage per reading for the Ballard Mark V, Temasek 1 kW, and Horizon H-12 are, respectively, +0.65%, +0.20%, and −0.14%, which are negligible errors. The primary characteristics of PEMFC stacks under changing reactant pressures and cell temperatures are used to evaluate the precision of the cropped optimized parameters. In the final phase of this endeavor, the sensitivity of the cropped parameters to the PEMFCs model’s performance is investigated using two machine learning techniques, namely, artificial neural network and Gaussian process regression models. The simulation results demonstrate that the RBMO approach extracts the PEMFCs’ appropriate parameters with high precision.

1. Introduction

Among clean source alternatives, fuel cells (FCs) are thought to be a relatively new and quickly developing application of renewable energy-based conversion technology [1,2]. Owing to their interesting features, such as their small size, high power density, low operating pressures and temperatures, and the absence of dynamic elements, FCs have promising applications in smart grids, distributed generations, and transportation [3,4]. This electrochemical device, called FC for short, is a type of equipment that utilizes hydrogen gas as an input fuel to transform chemical energy into electric energy and heat, applying oxygen/air as an oxidant [5,6].
FCs are classified into various varieties, each with its own set of properties and applications. Here are a few common types: (i) proton exchange membrane FCs (PEMFCs) [7,8,9]; (ii) solid oxide FCs, which operate at high temperatures (600–900 °C) and are suited for stationary power generation and combined heat and power systems [10,11]; (iii) molten carbonate FCs, which run at high temperatures (600–700 °C), which are employed for large-scale stationary applications and can use natural gas as fuel [12]; (iv) phosphoric acid FCs, which operate at moderate temperatures (150–200 °C) and are frequently employed in stationary power applications, such as in commercial buildings [13]; (v) alkaline FCs, which function at relatively low temperatures (60–90 °C) and were formerly utilized in space applications and some ground vehicles [14]; (vi) direct methanol FCs, which are commonly used in portable applications and small electronics [15]; (vii) and reversible FCs, which can operate both as an FC and as an electrolyzer and are used for energy-storage applications [16].
One of the most common types of FCs is PEMFCs. Because the temperature and supply pressure can affect the output voltage, which can range from 0.9 to 1.23 V/cell, a series of PEMFCs is linked in order to increase the output voltage to the desired level [9,17,18,19]. A stack is a collection of PEMFCs connected in series/parallel for specific applications. Furthermore, because of polarization losses, their output voltage shows a non-linear relationship with the drawn load current. Stated otherwise, the activation voltage drop causes the PEMFCs’ output voltage to drop quickly at first; then, the ohmic voltage drop causes it to decrease linearly, and eventually, concentration losses cause it to decline exceedingly [20,21,22].
The model of PEMFCs has a set of unknown parameters in its mathematical form that are not specified in the fabricators’ datasheets. To accurately simulate the real behavior of PEMFCs, certain parameters must be optimally estimated. As a result, numerous attempts have been made to fully define the model’s ungiven parameters for the PEMFCs units. These efforts can essentially be divided into two categories: traditional and soft computing-based optimization frameworks [23,24]. Recently, machine learning has been used to achieve the same goal, as publicized in [25,26,27,28,29].
There are many optimization frameworks used to define the unknown parameters of PEMFCs stacks [4,30,31]. Among these recent optimizers are the artificial rabbits optimizer (ARO) [3], bonobo optimizer [32], chaotic Harris hawks optimizer (CHHO) [33], converged moth search algorithm [34], circle search algorithm (CSA) [35], transient search optimizer (TSO) [36], grey wolf optimizer [37,38,39], grasshopper optimizer (GHO) [40], whale optimization algorithm (WOA) [41], artificial bee colony differential evolution (DE) shuffled complex optimizer (ABDEO) [42], manta rays foraging optimizer (MRFO) [43], sine–cosine crow search algorithm [17], hybrid artificial bee colony DE optimizer [44], improved Archimedes optimizer [45,46], and Kepler optimizer [47]. In addition, recently, a number of cutting-edge, well-known optimization strategies have been made available to address the issue of extracting PEMFC’s parameters, including the atom search optimizer [48], equilibrium optimizer [49], gradient-based optimizer [50], shark smell optimizer [51], teamwork optimizer [52], modified farmland fertility optimizer (MFFO) [53], honey badger optimizer [54], two novel approaches reported in [55], human memory optimizer [56], and a reliable exponential distribution optimizer [57].
The “no free lunch” (NFL) theory [58] states that there is not a single algorithm that can handle every engineering optimization problem because every optimization technique has benefits and drawbacks for different tasks. There is no definitive solution yet, and choosing between the problems of optimization methods X and Y can be challenging depending on a number of factors, such as the degree of non-linearity, non-convexity, multi-modality, separability of the control variables, high dimensionality, etc. Until such a response is obtained in these endeavors, the attempts will persist. As shown above, there has been much success in determining these parameters; nonetheless, there is always room for improvement in order to more precisely address the ideal PEMFC stack model values.
The red-billed blue magpie optimizer (RBMO), developed in 2024 by Fu et al., can be used to solve problems involving continuous engineering optimization [59]. The mutually beneficial and effective predation habits of red-billed blue magpies were a model for the RBMO. The red-billed blue magpie’s hunting, locating, pursuing, attacking, and food-storage activities are all modeled mathematically in the RBMO’s procedures. It should be mentioned that the RBMO is an algorithm that uses swarm intelligence and is motivated by the red-billed blue magpie’s hunting strategy [59].
In this study, the RBMO is used to extract the parameters of PEMFCs in three test cases under various scenarios, along with implementing three other recent optimizers, namely, the dandelion optimization algorithm (DOA) [60], sinh-cosh optimizer (SCHO) [61], and growth optimizer (GO) [62]. The necessary verifications, including comparisons, are made using some specific measures. In addition to that, the principal performance of PEMFCs units under varied conditions is investigated and discussed.
The DOA is used to solve problems involving a continuous optimization. The DOA mimics the wind-powered long-distance flight of dandelion seeds, which can be divided into three stages: rising, descending, and landing. Seeds rise spirally, descend slowly, and land randomly, with Brownian motion and Levy random walk describing their trajectory [60]. The DOA has been recently applied successfully to solve electric power system problems [63,64,65,66,67]. The GO is a metaheuristic algorithm that solves both continuous and discrete global optimization problems. This technique mimics natural growth processes, relying on biologically inspired mechanisms to effectively explore and exploit the solution space [62]. The GO stands out as a powerful tool in the field of optimization, combining natural inspiration with advanced algorithmic techniques to solve complex problems effectively. GO has recently been used to address a few engineering issues, according to published reports [62,68,69]. The SCHO is a metaheuristic algorithm inspired by the mathematical properties of the hyperbolic sine and cosine functions [61,70]. It is designed to solve optimization problems by effectively balancing exploration and exploitation of the search space. The SCHO is a flexible and powerful tool that can be used to solve a wide range of optimization issues. It makes use of the mathematical characteristics of hyperbolic functions to effectively traverse challenging environments and is applied to solve a few engineering and medical problems [70,71,72,73,74].
Now, let us highlight the key contributions of this article: (i) assessing RBMO’s performance to optimally give the values of unknown parameters in PEMFCs’ model by implementing three recent optimizers, namely, DOA, GO, and SCHO; (ii) carefully examining three real-world study cases such as Ballard Mark V, Temasek 1 kW, and Horizon H-12 under various operating scenarios; and (iii) conducting numerous comparisons and uncertainty assessments for the obtained results.

2. Mathematical Formulation of PEMFCs’ Modeling

This section discusses the basic formulation of the PEMFCs’ model. The modeling approach for PEMFCs involves developing computer and mathematical models that represent the physical and electrochemical processes that take place inside the FCs. The mathematical model developed to simulate the performance of PEMFCs is referred to as Mann’s model [75]. This model is widely used in the PEMFCs’ research field. The PEMFCs stacks contain a variety of voltage drops, including concentration voltage ( V c ), ohmic or resistive voltage ( V Ω ), and activation voltage ( V a ).
Three different regions are identified on the I-V polarization curve of a single PEMFC: concentration, ohmic, and activation losses. Activation losses, a reflection of the slow electrochemical reactions at first, are the cause of the fast output voltage decay of the PEMFC when it is first starting up under light load. The entire resistance that the protons and electrons encounter, which is represented by ohmic losses, causes the output voltage to subsequently decrease linearly. Because of the increasing water content at higher load conditions, the output voltage drops quickly and lowers the reactant concentration in both electrodes. These voltage drops, which contribute to the overall voltage loss in the FC, may have a significant impact on the system’s performance and efficacy. A thorough understanding of these voltage losses, as well as specific attempts to reduce them, are required to improve the efficacy of PEMFCs stacks. Scientists and engineers use a range of strategies to achieve this goal, including catalyst development, improved flow field designs, and improved reactant gas management. As a result, the PEMFC’s terminal voltage is given as follows (1):
V f c = E N V Ω + V a + V c
For working temperatures below 100 °C, E N , which represents the reversible open-circuit voltage, may be computed using the formula revealed in (2).
E N = 1.2290 8.5 × 1 0 4 T c 298.15 + 4.385 × 1 0 5 T c ln P H 2 . P O 2
Equation (3) gives an estimate of the activation voltage loss ( V a ).
V a = [ ξ 1 + T c ξ 2 + ξ 3 ln C 02 + ξ 4 ln I f c ]
The concentration of oxygen (mol/cm3), which is defined as
C O 2 = P O 2 5.08 × 10 6 e 498 T
The equivalent resistance of the FC is used to calculate V Ω , which has the following definition:
V Ω = I f c ( R m + R c )
Equations (6) and (7) can be used for R m calculations.
R m = ρ m l A m
where
ρ m = 181.6 1 + 0.03 J + 0.062 T c 303 2 J 2.5 λ 0.634 3 J e 4.18 . 1 303 T c
The formula shown in (8) can be used to predict V c .
V c = β ln 1 J J m a x
Typically, the PEMFC stack is made up of a number of series n-cells, and the voltage across the stack is expressed as follows, assuming all cells behaved with the same performance.
V s t a c k = n . V f c
The aforementioned calculation is applied under the presumption that each cell behaves in the same way and that the resistors connecting the cells are ignored. In a closer look at the above mathematical PEMFCs’ model beginning from (1) to (9), it is clear from the above-mentioned formulas that for obtaining a fully defined electrochemical-based model, seven unknown parameters ( ξ 1 , ξ 2 , ξ 3 , ξ 4 , λ , R C , and β ) should be estimated. Model validation, parameter optimization, and refinement cycles are all parts of the iterative process of parameter estimation for Mann’s model. To obtain precise and trustworthy parameters values that accurately and reliably describe the behavior of the PEMFCs units under investigation, it is necessary to combine experimental data, computational modeling, and optimization techniques. Certainly parameter optimization approaches, whether applied in Mann’s model or any other mathematical model, pursue the recognized parameters values that minimize the discrepancy between model predictions and experimental results of voltage dataset points.

3. Problem Formulation and Constraints

The assessment criteria, in addition to the identification procedure, play an important role in precisely determining the best values of the unknown parameters and fitting the calculated I-V curve to the experimental one. A well-chosen fitness function (FF) streamlines the process of determining parameters and distinguishes between different model identification techniques in terms of the permissible ranges of results both quantitatively and qualitatively. The most common and widely used FF available in the literature are given in (10) for a fair comparisons with others. The sum of squared deviations (SSD), which is estimated using (11), is adapted to mathematically describe the FF used in this study [1,3,17,23,26,34,36,37,40,41,42,43,47,54].
F F = M i n i m i z e j = 1 N m V m , j V e , j 2
The above-mentioned FF is subjected to a set of operating and practical min/max constraints; they are mentioned in (11).
s . t .   ξ i m i n   ξ i ξ i m a x   ,   i 1 , 2,3 , 4 λ m i n   λ   λ m a x   R C m i n   R C   R C m a x β m i n     β   β m a x
As the reader can see, seven unknown parameters must be precisely defined in order to verify the simulation’s accuracy and provide a good representation of the behaviors of PEMFCs stack under study. It could be important to note that the RBMO is self-constrained all inequality constraints, confirming that there is no chance of breaks outside of the working min/max range, which means no need to penalize the FF.

4. Procedures of the RBMO

Swarm intelligence is used by the red-billed blue magpie (RBMO) algorithm to simulate its activities of hunting, pursuing, attacking, and storing food. Fruits, insects, and small vertebrates are the main food sources for this adaptable predator. When roaming in small groups, they locate food sources more easily and collaborate to become more productive. In addition, they store food for later use and obtain nutrition from the earth and trees. The phases of exploration and exploitation are included in the RBMO mathematical model. Like other meta-heuristic competing optimizers, the RBMO initializes a random candidate solution in the manner described below, as depicted in (12).
X i , j = L B , j + r a n d 1 H B , j L B , j ,   j D   and   j   P o p
Red-billed blue magpies forage for food in small groups (2–5 agents), as described in (13), or clusters (more than 10 agents), as revealed in (14), employing walking, leaping, and tree-searching techniques. Because of their adaptability, food is always available despite shifting environmental conditions and resource availability.
X i t + 1 = X i t + r 2 1 p m = 1 p X m t X r s t
X i t + 1 = X i t + r 3 1 C m = 1 C X m t X r s t
The red-billed blue magpie is an expert predator that uses a variety of hunting techniques, such as flying, leaping, and rapid pecking. Its main objectives are larger prey in clusters and small animals or plants in small groups, which is described mathematically in (15) and (16), respectively. Thus, it can secure food in a variety of situations thanks to its adaptable nature, as specified in (15) and (16).
X i t + 1 = X F t + w t . r 3 1 C m = 1 C X m t X r s t
X i t + 1 = X F t + w t . r 2 1 p m = 1 p X m t X r s t
The most crucial parameter in the RBMO, w t , is obtained and updated repeatedly, as explained in (17), and it is in charge of striking a balance between exploration and exploitation (this factor is modified from the original paper [59], which is 1 t T m a x . e x p 2 t T m a x ).
w t = e x p b t T m a x b
As illustrated in (18), red-billed blue magpies conceal food and save data for globally ideal values.
X i t + 1 = X i t     F F i o l d >   F F i n e w X i t + 1   O t h e r w i s e
The time complexity in the worst case of BIG O(…) of the RBMO comprises the following: (i) the initialization complexity and (ii) the main loop complexity. The reader may refer to Figure 1 for further illustrations. The initialization complexity of the RBMO is O P o p . D . The main loop complexity of the RBMO is O 2 . P o p . D . T m a x for exploration and exploitation stages. Thus, the overall time complexity, on the other hand, is equal to O P o p . D + 2 . P o p . D . T m a x O P o p . D . T m a x .
Further details about the RBMO and its distinguishing characteristics compared to other meta-heuristics, along with its pseudo-code, are available in [59]. Figure 1 depicts the general procedures of the RBMO, showing the algorithm’s framework.

5. Numerical Simulations, Demonstrations, and Validations

This section evaluates the performance of the RBMO in conjunction with implementing three recent optimizers, named the DOA [60], the GO [62], and the SCHO [61], with other recently published results using three test cases under various scenarios. In these test scenarios, Ballard Mark V, Temasek 1 kW, and Horizon H-12 are the PEMFCs units that were used to attain the above-mentioned goal. Numerous autonomous attempts were conducted, and the most optimally cropped outcomes are documented in the upcoming subsection. The numerical simulations were performed using MATLAB Software version R2023b on a laptop running with Windows 10 Enterprise and outfitted with an Intel(R) Core(TM) i7-7700HQ CPU @ 2.80GHz, 2.80 GHz, and 16 GB of RAM. The implemented optimizers’ adjusted settings for every test case are 50 populations and a maximum of 100 iterations for the DOA, the SCHO, and the GO for fair comparisons. However, the adapted control settings of the RBMO are b = 2.2 ,   T m a x = 100 , and P o p = 50 . It might be important to note that these values were arrived at through a process of trial-and-error procedures. Furthermore, it should be noted that the simulation time is not regarded as a crucial component in the suggested strategy because the design parameters are decided offline. Only the RBMO’s findings (data tip) are displayed in plots in order to provide a clear vision and because the other optimizers’ results are extremely comparable to them except for convergence patterns. The following sections of the study go into greater depth on the numerical results and comparisons among the optimizers that were used and implemented.

5.1. Test Cases and Their Associated Results

The technical specifications of several commercial PEMFCs stacks are revealed in Table 1 for the test cases under study, and their principle data are obtained from [17,23,37,40,41,42,47,54].
Table 2 recapitulates the lower and higher bounds of the PEMFCs’ stack unknown parameters, which are frequently used in cutting-edge technology and in the state of the art [1,3,17,23,26,34,36,37,40,41,42,43,47,54] to ensure fair comparisons with other competitors.
The patterns of the RBMO, the DOA, the GO, and the SCHO’s convergence signatures for the test cases under study, which show the FF’s final cropped best values as the data tip by the RBMO only, are shown in Figure 2a–c. The reader may observe that the RBMO approach was able to obtain extremely competitive results with fine convergence and steady progress at iteration 50. When the parallel pool is enabled with four working processors, the CPU average processing times over 100 implementations of the RBMO approach for finalizing this optimization task are 2.1 s, 3.1 s, and 2.5 s for Ballard Mark V, Temasek 1 kW, and Horizon H-12, respectively.

5.1.1. Test Case 1: Ballard Mark V

The Ballard Mark-V 5 kW (its specs and datasheet are given in Table 1, the third row from the top of the table) PEMFC’s stack modeling is adapted to the RBMO-based technique plus DOA, GO, and SCHO; Table 3 provides the appraised model’s parameters and the solution’s fitness (the last row of the table). The outcomes of the competitors from the literature are placed alongside them. It was clear that the DOA method outperformed the closest best FF, with a value of 0.853608 V 2 . Figure 3a,b of the Ballard Mark-V shows the calculated model versus the actual readings of the V-I plots and internal voltage drops inside the stack, respectively. In this particular experiment, there are thirteen observations for voltage dataset points. It is evident that the computed voltage and the measured voltages fit each other accurately, as confirmed in Figure 3a.
With the purpose of investigating how temperature and pressure affect PEMFCs’ performance, the voltages are again calculated at 45, 60, and 75 degrees Celsius, with a constant pressure of ( P H 2 = / P O 2 = = 1.5/1 atm). Figure 3c shows the findings of the temperature's impact on the performance of the PEMFCs graphically. The PEMFC’s output voltage rises in conjunction with temperature. However, Figure 3d illustrates how altering partial pressures affects stack performance while maintaining a constant temperature of 343 Kelvin. The voltage rises in concert with the pressure.

5.1.2. Test Case 2: Temasek 1 kW

The Temasek 1 kW PEMFCs stack (its specs and datasheet are given in Table 1, the second row from the bottom of the table), modeling is adapted to the RBMO-based technique; Table 4 provides the estimated model parameters and the solution’s fitness (the last row of the table). The outcomes of the competitor techniques from the literature, plus the two other methods implemented, are positioned alongside them. It was clear that the RBMO-based framework outperformed the closest best FF with a value of 0.597504 V 2 , much less than the other competitors' results. Similar to the above-mentioned test case 1, Figure 4a,b show the calculated model versus the actual readings of the V-I plots and internal voltage drops inside the stack, respectively. In this particular experiment, there are twenty observations. It is evident that the computed voltage and the measured voltages fit each other accurately. With the purpose of investigating how temperature and pressure affect PEMFCs’ performance, the voltages are again calculated at 45, 60, and 75 degrees Celsius, with a constant pressure of ( P H 2 = / P O 2 = 0.5/0.5 atm). Figure 4c shows the findings of the temperature impact on the performance of the PEMFCs graphically. The PEMFC’s output voltage rises in conjunction with temperature. However, Figure 4d illustrates how altering partial pressures affects stack performance while maintaining a constant temperature of 323 Kelvin. The voltage rises in concert with the pressure.

5.1.3. Test Case 3: Horizon H-12

In this scenario, 13 FCs connected in series are investigated. The H-12 PEMFC has a thickness of 25 m and a surface area of 8.1 cm2. Its specifications and datasheet are given in Table 1, the last row of this table. This device has a nominal power of 12 W (7.8 V and 1.5 A), a hydrogen input pressure of 0.45–0.55 bar, and a maximum operating temperature of 55 °C Overall, 246.9 mA/cm2 is the greatest current density. Figure 3c depicts the RBMO’s convergence pattern, which indicates the final cropped value of the SSD over 100 distinct runs. Table 5 contains further information and comparisons between the RBMO, implementing algorithms, and other standing optimization techniques. In this situation, the RBMO-based method produced significantly superior results than the other analyzed algorithms. The RBMO method reduced the FF by roughly 37.5% less than the other analyzed techniques.
Meanwhile, Figure 5a depicts the comparison of the calculated and measured power. In this situation, there are fifteen experimental readings. These two figures ensure the accuracy of the theoretical representation and the usefulness of the RBMO approach in producing results that correspond to the actual operating performance of the PEMFC stack. The internal voltage drops (VDs) losses in the stack are depicted graphically in Figure 5b, comprising the total voltage drop (TVD). The losses are seen to grow as the current increases.
Several simulations are executed at different temperatures of 30 °C, 40 °C, and 55 °C, with pressure constants of P O 2 = 1 bar and P H 2 = 0.4935 bar. Figure 5c,d illustrate the values of the output voltage after altering the temperature. When the temperature rises, so do the output voltage and power. Furthermore, the simulations are performed at various hydrogen pressures ( P H 2 = 0.4 and 0.55 bar) while maintaining a constant temperature of 29 °C and an oxygen pressure of 1 bar. Figure 5c,d demonstrates the voltage and power changes as a function of pressure change. When the pressure increases, the voltage and power increase, although with a lower effect than when the temperature changes.

5.2. Validations by Some Measures

Table 6 shows the actual and estimated dataset voltage points using the RBMO. In this table, the average absolute voltage errors between actual and model calculations are 0.19815, 0.10947, and 0.05259 V for the Ballard Mark V, Temasek 1 kW, and Horizon H-12 test cases, respectively, as indicated. It is also worth noting that the number of readings for the Ballard Mark V, Temasek 1 kW, and Horizon H-12 test cases are 13, 20, and 15, respectively. It is likely to notice that the reported numbers only have four digits, which can lead to a little difference between the precise value of SSDs, as generated by the MATLAB calculations, and subsequent manual assessments. The biased voltage error V B E , j and mean absolute error M A E are computed as follows:
V B E , j = V m ,   j V e ,   j ,     j N m
M A E = j = 1 N m V B E , j N m
Figure 6a–c depicts the percentage of biased voltage against the drawn current at measured records for the three studied test cases. Upon closely examining Figure 6a–c, it becomes evident that the highest percentages of biased voltage per reading for Ballard Mark V, Temasek 1 kW, and Horizon H-12 are equal to +0.65%, +0.20%, and −0.14%, respectively.
In addition to the above-mentioned results, the mean absolute percentage errors (MAPE) for voltage dataset are estimated using (21).
M A P E = k = 1 N m V m , k V e , k V m , k N m × 100
For the purpose of quantification, the Ballard Mark V, Temasek 1 kW, and Horizon H-12 test cases’ respective MAPE values are 0.7598%, 0.8259%, and 0.6003%. It appears that the audiences are aware of how tiny and unimportant the reported particular values of MAEs and MAPEs are. The latter makes it abundantly evident that the calculated and real voltage datasets match well, demonstrating the applicability of the applied RBMO in producing the optimal values for the seven unknown parameters.
The Pearson correlation coefficient ρ V m , V e , between measured and simulated datasets, is calculated using the formula given in (22). The values of ρ V m , V e for the Ballard Mark V, Temasek 1 kW, and Horizon H-12 test cases are 0.9978, 0.9930, and 0.9915, respectively. The aforementioned discussion demonstrates the RBMO-based optimization framework’s dominance after adaptation. This has been performed for each of the three test scenarios that have been studied.
ρ V m   ,   V e   = 1 N m 1 . i = 1 N m V m , i V ¯ m σ m . V e , i V ¯ e σ e
Furthermore, a p-value with a 5% significance level is generated for the model’s derived voltage points versus the observed voltage datasets, which have values of 6.75 × 10 12 , 0.00, and 2.80 × 10 14 , for the studied PEMFCs units. These negligible numbers, which are very close to zero, indicate a significant joining and a low likelihood of detecting the null hypothesis.

5.3. Sensitivity Analysis of Extracted Parameters

Sensitivity analysis (SA) quantifies the variations in results due to parameter variations [80]. It aims to rank, screen, and map regions producing extreme output values [81]. Global SA (GSA) evaluates model sensitivity for input/parameter variations. The GSA can be used before model calibration to identify significant parameters and extreme values. The GSA methods are based on probabilistic or random sampling that can be applied to any mathematical model with multiple inputs and outputs. The main concept is to optimize the information extracted from an optimal number of model runs by defining the evaluated locations in the parameter space. These methods are tailored to models with moderate to large inputs and are ideal for identifying maximum or minimum outputs. In this effort, two models, namely, G P R and a n n models, are implemented to perform the GSA. First-order ( S i ) and total order ( S T i ) indices are collected for both G P R and A N N models as well. They are called g p r S i , g p r S T i , a n n S i , and a n n S T i . In most circumstances, knowing S i and S T i is sufficient to ascertain how sensitive the examined function is to each of the individual input parameters ( ξ 1 , ξ 2 , ξ 3 , ξ 4 , λ , R C , and β ). The sensitivity range is ±10% of the original estimated best parameters values by the RBMO, and the number of samples is equal to 100 when using the Sobol method and Jansen estimator.
With a closer look at Table 7, the ranking of the seven parameters based on their significant impact on the performance of the model, it can be concluded that, in accordance with g p r and a n n   i n d i c e s , the model is highly sensitive to variations in ε 2 ,   ε 1 , and ε 3 ; moderately sensitive to ε 4 ; and less sensitive to λ , β , and R c for Ballard Mark V and Temasek 1 kW PEMFC generating units. On the other hand, for Horizon H-12, almost the same conclusion can be summed up for ε 1 ,   ε 2 , and ε 3 , with moderate sensitivity to β and less sensitivity to λ , ε 4 , and R c . Bear in mind the H-12 unit has limited power and a limited operating range. One more note is that both g p r and a n n sensitivity-based models produce the same conclusion.

6. Conclusions and Suggestions

This study employs the RBMO to extract the unknown parameters, i.e., ( ξ 1 , ξ 2 , ξ 3 , ξ 4 , λ , R C , and β ) of the PEMFCs stacks under varying operating conditions. Furthermore, the fitness function is derived from the sum of square deviations between the computed and actual PEMFC stack voltages. The RBMO’s results on Ballard Mark V, Temasek 1 kW, and Horizon H-12 are three popular PEMFC units that were compared to three implemented algorithms, namely, DOA, GO, and SCHO, plus other recently published competitors. Some specific measures, including biased voltage errors, mean absolute voltage errors, mean absolute percentage errors, Pearson correlation coefficients, and p-values, are used to assess the outcomes of the RBMO’s procedure. It is evident that the highest percentages of biased voltage per reading for Ballard Mark V, Temasek 1 kW, and Horizon H-12 are, respectively, +0.65%, +0.20%, and −0.14%, which are negligible errors. Based on the three test cases under investigation, the comparison results show that the used RBMO is the best approach among the competing methods. In the last stage of this work, sensitivity analysis based on Gaussian process regression and artificial neural network models is used to prioritize the ranks of cropped defined parameters. This indicates that the PEMFCs’ model is highly sensitive to variations of ε 2 ,   ε 1 and ε 3 ; moderately sensitive to ε 4 ; and less sensitive to λ , β , and R c . Based on the encouraging outcomes of its use, the RBMO can be applied to other domains like the parameter estimation of solar cells, tuning load frequency controllers of microgrids comprising 100% renewables, wind speed forecasts, and other renewable energy applications.

Author Contributions

Conceptualization, A.M.A. and A.A.E.-F.; methodology, A.M.A.; software, A.A.E.-F.; validation, A.M.A. and A.A.E.-F.; formal analysis, A.M.A.; investigation, A.A.E.-F.; resources, A.M.A.; data curation, A.M.A.; writing—original draft preparation, A.M.A.; writing—review and editing, A.A.E.-F.; visualization, A.A.E.-F.; supervision, A.A.E.-F.; project administration, A.M.A.; funding acquisition, A.M.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Deanship of Scientific Research at Northern Border University, Arar, KSA, through the project number “NBU-FFR-2024-2968-03”.

Institutional Review Board Statement

Not applicable. This article does not contain any studies with human participants or animals.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

The authors extend their appreciation to the Deanship of Scientific Research at Northern Border University, Arar, KSA, for funding this research work through project number “NBU-FFR-2024-2968-03”.

Conflicts of Interest

The author confirms that there are no conflicts of interest.

Abbreviations

FCsFuel cells.
PEMFCsProton exchange membrane FCs.
SSDSum of squared deviations.
DOADandelion optimization algorithm.
RBMORed-billed blue magpie optimizer.
SSOShark smell optimizer.
CHHOChaotic Harris hawks optimizer.
QOBOQuasi oppositional bonobo optimizer.
MFFOModified farmland fertility optimizer.
MAEOModified artificial ecosystem optimizer.
CSACircle search algorithm.
NNONeural network optimizer.
MRFOManta rays foraging optimizer.
WOAWhale optimization algorithm.
GhOGrasshopper optimizer.
TSOTransient search optimizer.
AROArtificial rabbits optimizer.
HBOHoney badger optimizer.
DEDifferential evolution.
ABDEOArtificial bee colony DE optimizer.
AEOArtificial ecosystem-based optimizer.
ELBOEnhanced Lévy flight bat optimizer.
SCEShuffled complex evolution.
PFAPathfinder algorithm.
ASSAAdaptive sparrow search algorithm.
PSOParticle swarm optimizer.
GAGenetic algorithm.
IAHAImproved artificial hummingbird algorithm.
FFFitness function.
SCHOSinh cosh optimizer.
GOGrowth optimizer.
annArtificial neural network.
GPRGaussian process regression.

Nomenclature

V c Concentration voltage.
V Ω Ohmic or resistive voltage.
V a Activation voltage.
P O 2 and P H 2 Partial pressures of O 2 and H 2 , respectively.
T c Cell temperature (K).
I f c Drawn load current from the FCs.
ξ 1   t o   ξ 4 Parametric Coefficients.
C O 2 Concentration of oxygen (mol/ c m 3 ).
R C and R m Contact and the membrane resistances, respectively.
A m Active area of the cell ( c m 2 ).
J Actual current density (mA/ c m 2 ).
λ Membrane water content.
ρ m Resistivity of the membrane (Ω.cm).
β Constant number.
J m a x Ultimate current density (A/cm2).
n Number of series-n fuel cells.
V m , j Measured voltage value of j t h reading.
V e , j Estimated voltage value of j t h reading.
I m , j Measured voltage value of j t h reading.
N m Number of readings taken from the experimental setup.
L B Lower limits of the uncertain PEMFCs parameters.
H B Higher limits of the uncertain PEMFCs parameters.
ξ i m i n , ξ i m a x Min/max range of parametric coefficients ξ i .
λ m i n , λ m a x Min/max range of water content, λ parameter.
R C m i n , R C m a x Min/max range of contact resistance.
β m i n , β m a x Min/max range of β constant.
V ¯ m Average of experimental voltage points.
V ¯ e Average of estimated voltage points.
σ m Standard deviation of experimental voltage points.
σ e Standard deviation of estimated voltage points.
r 1 , r 2 , r 3 Random variables in the range [0, 1].
P o p Population size.
t Iteration counter.
T m a x Maximum iterations.
D Number of design variables.
X i t Position of the candidate in t t h iteration.
X i t + 1 New position of the agent.
p Number of agents in small groups.
C Number of agents in clusters.
X m t Agent chosen randomly.
X r s t Randomly selected agents.
X F t Position of the food.
F F i o l d Fitness value before the position updated for ith agent.
F F i n e w Fitness value after the position updated for ith agent.

References

  1. Ismael, I.; El-Fergany, A.A.; Gouda, E.A.; Kotb, M.F. Cooperation search algorithm for optimal parameters identification of SOFCs feeding electric vehicle at steady and dynamic modes. Int. J. Hydrogen Energy 2023, 50, 1395–1407. [Google Scholar] [CrossRef]
  2. Martsinchyk, A.; Milewski, J.; Dybiński, O.; Szczęśniak, A.; Siekierski, M.; Świrski, K. Experimental investigation of novel molten borate fuel cell supported by an artificial neural network for electrolyte composition selection. Energy 2023, 279, 127921. [Google Scholar] [CrossRef]
  3. Riad, A.J.; Hasanien, H.M.; Turky, R.A.; Yakout, A.H. Identifying the PEM Fuel Cell Parameters Using Artificial Rabbits Optimization Algorithm. Sustainability 2023, 15, 4625. [Google Scholar] [CrossRef]
  4. Abdel-Basset, M.; Mohamed, R.; Abouhawwash, M. On the facile and accurate determination of the highly accurate recent methods to optimize the parameters of different fuel cells: Simulations and analysis. Energy 2023, 272, 127083. [Google Scholar] [CrossRef]
  5. Hernández-Gómez, Á.; Langarica-Cordoba, D.; Martinez-Rodriguez, P.R.; Guilbert, D.; Ramirez, V.; Saldivar, B. PEM Fuel Cell Emulators: A Review. Electronics 2023, 12, 2812. [Google Scholar] [CrossRef]
  6. Wang, Y.; Zhang, Y.; Zhang, C.; Zhou, J.; Hu, D.; Yi, F.; Fan, Z.; Zeng, T. Genetic algorithm-based fuzzy optimization of energy management strategy for fuel cell vehicles considering driving cycles recognition. Energy 2023, 263, 126112. [Google Scholar] [CrossRef]
  7. Wang, Y.; Seo, B.; Wang, B.; Zamel, N.; Jiao, K.; Adroher, X.C. Fundamentals, materials, and machine learning of polymer electrolyte membrane fuel cell technology. Energy AI 2020, 1, 100014. [Google Scholar] [CrossRef]
  8. Ding, R.; Zhang, S.; Chen, Y.; Rui, Z.; Hua, K.; Wu, Y.; Li, X.; Duan, X.; Wang, X.; Li, J.; et al. Application of machine learning in optimizing proton exchange membrane fuel cells: A review. Energy AI 2022, 9, 100170. [Google Scholar] [CrossRef]
  9. Tian, P.; Liu, X.; Luo, K.; Li, H.; Wang, Y. Deep learning from three-dimensional multiphysics simulation in operational optimization and control of polymer electrolyte membrane fuel cell for maximum power. Appl. Energy 2021, 288, 116632. [Google Scholar] [CrossRef]
  10. Abaza, A.; El Sehiemy, R.A.; El-Fergany, A.; Bayoumi, A.S.A. Optimal parameter estimation of solid oxide fuel cells model using bald eagle search optimizer. Int. J. Energy Res. 2022, 46, 13657–13669. [Google Scholar] [CrossRef]
  11. Rizk-Allah, R.M.; El-Hameed, M.A.; El-Fergany, A.A. Model parameters extraction of solid oxide fuel cells based on semi-empirical and memory-based chameleon swarm algorithm. Int. J. Energy Res. 2021, 45, 21435–21450. [Google Scholar] [CrossRef]
  12. Lu, Z.; Zhang, H.; Duan, L.; Zhang, R.; Lu, C.; Desideri, U. Experimental analysis of ammonia as input fuel for molten carbonate fuel cell. Fuel 2024, 357, 130032. [Google Scholar] [CrossRef]
  13. Verma, V.; Choudhury, S.R.; Rathour, V.; Choudhury, S.R.; Ganesan, V. Hollow core mesoporous carbon spheres as catalyst support for improved platinum utilization in phosphoric acid fuel cells. Microporous Mesoporous Mater. 2024, 367, 113005. [Google Scholar] [CrossRef]
  14. Hamada, A.T.; Orhan, M.F.; Kannan, A.M. Alkaline fuel cells: Status and prospects. Energy Rep. 2023, 9, 6396–6418. [Google Scholar] [CrossRef]
  15. Okech, G.; Emam, M.; Mori, S.; Ahmed, M. Experimental study on the effect of new anode flow field designs on the performance of direct methanol fuel cells. Energy Convers. Manag. 2024, 301, 117988. [Google Scholar] [CrossRef]
  16. Yang, Y.; Lei, J.; Huang, X.; Liao, Z.; Liu, Y.; Tu, Z. Recent Development in Reversible Solid Oxide Fuel Cells: Theory, Integration and Prospective. ChemElectroChem 2024, 11, e202300593. [Google Scholar] [CrossRef]
  17. Alizadeh, M.; Torabi, F. Precise PEM fuel cell parameter extraction based on a self-consistent model and SCCSA optimization algorithm. Energy Convers. Manag. 2021, 229, 113777. [Google Scholar] [CrossRef]
  18. Shaheen, M.A.; Ullah, Z.; Hasanien, H.M.; Tostado-Véliz, M.; Ji, H.; Qais, M.H.; Alghuwainem, S.; Jurado, F. Enhanced transient search optimization algorithm-based optimal reactive power dispatch including electric vehicles. Energy 2023, 277, 127711. [Google Scholar] [CrossRef]
  19. Alharbi, A.G.; Olabi, A.; Rezk, H.; Fathy, A.; Abdelkareem, M.A. Optimized energy management and control strategy of photovoltaic/PEM fuel cell/batteries/supercapacitors DC microgrid system. Energy 2024, 290, 130121. [Google Scholar] [CrossRef]
  20. Sultan, H.M.; Menesy, A.S.; Kamel, S.; Turky, R.A.; Hasanien, H.M.; Al-Durra, A. Optimal Values of Unknown Parameters of Polymer Electrolyte Membrane Fuel Cells Using Improved Chaotic Electromagnetic Field Optimization. IEEE Trans. Ind. Appl. 2021, 57, 6669–6687. [Google Scholar] [CrossRef]
  21. Pei, P.; Meng, Y.; Chen, D.; Ren, P.; Wang, M.; Wang, X. Lifetime prediction method of proton exchange membrane fuel cells based on current degradation law. Energy 2023, 265, 126341. [Google Scholar] [CrossRef]
  22. Mian, S.H.; Nazir, M.S.; Ahmad, I.; Khan, S.A. Optimized nonlinear controller for fuel cell, supercapacitor, battery, hybrid photoelectrochemical and photovoltaic cells based hybrid electric vehicles. Energy 2023, 283, 129121. [Google Scholar] [CrossRef]
  23. Gouda, E.A.; Kotb, M.F.; El-Fergany, A.A. Investigating dynamic performances of fuel cells using pathfinder algorithm. Energy Convers. Manag. 2021, 237, 114099. [Google Scholar] [CrossRef]
  24. Korkmaz, S.A.; Çetinkaya, S.A.; Yuksel, O.; Konur, O.; Erginer, K.E.; Colpan, C.O. Comparison of various metaheuristic algorithms to extract the optimal PEMFC modeling parameters. Int. J. Hydrogen Energy 2023, 51, 1402–1420. [Google Scholar] [CrossRef]
  25. Han, E.; Ghadimi, N. Model identification of proton-exchange membrane fuel cells based on a hybrid convolutional neural network and extreme learning machine optimized by improved honey badger algorithm. Sustain. Energy Technol. Assess. 2022, 52, 102005. [Google Scholar] [CrossRef]
  26. Yang, B.; Zeng, C.; Wang, L.; Guo, Y.; Chen, G.; Guo, Z.; Chen, Y.; Li, D.; Cao, P.; Shu, H.; et al. Parameter identification of proton exchange membrane fuel cell via Levenberg-Marquardt backpropagation algorithm. Int. J. Hydrogen Energy 2021, 46, 22998–23012. [Google Scholar] [CrossRef]
  27. Huo, W.; Li, W.; Zhang, Z.; Sun, C.; Zhou, F.; Gong, G. Performance prediction of proton-exchange membrane fuel cell based on convolutional neural network and random forest feature selection. Energy Convers. Manag. 2021, 243, 114367. [Google Scholar] [CrossRef]
  28. Legala, A.; Zhao, J.; Li, X. Machine learning modeling for proton exchange membrane fuel cell performance. Energy AI 2022, 10, 100183. [Google Scholar] [CrossRef]
  29. Su, D.; Zheng, J.; Ma, J.; Dong, Z.; Chen, Z.; Qin, Y. Application of Machine Learning in Fuel Cell Research. Energies 2023, 16, 4390. [Google Scholar] [CrossRef]
  30. Wilberforce, T.; Rezk, H.; Olabi, A.; Epelle, E.I.; Abdelkareem, M.A. Comparative analysis on parametric estimation of a PEM fuel cell using metaheuristics algorithms. Energy 2023, 262, 125530. [Google Scholar] [CrossRef]
  31. Abdel-Basset, M.; Mohamed, R.; El-Fergany, A.; Chakrabortty, R.K.; Ryan, M.J. Adaptive and efficient optimization model for optimal parameters of proton exchange membrane fuel cells: A comprehensive analysis. Energy 2021, 233, 121096. [Google Scholar] [CrossRef]
  32. Fathy, A.; Rezk, H.; Alharbi, A.G.; Yousri, D. Proton exchange membrane fuel cell model parameters identification using Chaotically based-bonobo optimizer. Energy 2023, 268, 126705. [Google Scholar] [CrossRef]
  33. Menesy, A.S.; Sultan, H.M.; Selim, A.; Ashmawy, M.G.; Kamel, S. Developing and applying chaotic harris hawks optimization technique for extracting parameters of several proton exchange membrane fuel cell stacks. IEEE Access 2019, 8, 1146–1159. [Google Scholar] [CrossRef]
  34. Sun, S.; Su, Y.; Yin, C.; Jermsittiparsert, K. Optimal parameters estimation of PEMFCs model using converged moth search algorithm. Energy Rep. 2020, 6, 1501–1509. [Google Scholar] [CrossRef]
  35. Qais, M.H.; Hasanien, H.M.; Turky, R.A.; Alghuwainem, S.; Loo, K.-H.; Elgendy, M. Optimal PEM Fuel Cell Model Using a Novel Circle Search Algorithm. Electronics 2022, 11, 1808. [Google Scholar] [CrossRef]
  36. Hasanien, H.M.; Shaheen, M.A.; Turky, R.A.; Qais, M.H.; Alghuwainem, S.; Kamel, S.; Tostado-Véliz, M.; Jurado, F. Precise modeling of PEM fuel cell using a novel Enhanced Transient Search Optimization algorithm. Energy 2022, 247, 123530. [Google Scholar] [CrossRef]
  37. Ali, M.; El-Hameed, M.; Farahat, M. Effective parameters’ identification for polymer electrolyte membrane fuel cell models using grey wolf optimizer. Renew. Energy 2017, 111, 455–462. [Google Scholar] [CrossRef]
  38. Miao, D.; Chen, W.; Zhao, W.; Demsas, T. Parameter estimation of PEM fuel cells employing the hybrid grey wolf optimization method. Energy 2020, 193, 116616. [Google Scholar] [CrossRef]
  39. Diab, A.A.Z.; Abdul-Ghaffar, H.I.; Ahmed, A.A.; Ramadan, H.A. An effective model parameter estimation of PEMFCs using GWO algorithm and its variants. IET Renew. Power Gener. 2022, 16, 1380–1400. [Google Scholar] [CrossRef]
  40. El-Fergany, A.A. Electrical characterisation of proton exchange membrane fuel cells stack using grasshopper optimiser. IET Renew. Power Gener. 2018, 12, 9–17. [Google Scholar] [CrossRef]
  41. El-Fergany, A.A.; Hasanien, H.M.; Agwa, A.M. Semi-empirical PEM fuel cells model using whale optimization algorithm. Energy Convers. Manag. 2019, 201, 112197. [Google Scholar] [CrossRef]
  42. Hachana, O. Accurate PEM fuel cells parameters estimation using hybrid artificial bee colony differential evolution shuffled complex optimizer. Int. J. Energy Res. 2022, 46, 6383–6405. [Google Scholar] [CrossRef]
  43. Selem, S.I.; Hasanien, H.M.; El-Fergany, A.A. Parameters extraction of PEMFC’s model using manta rays foraging optimizer. Int. J. Energy Res. 2020, 44, 4629–4640. [Google Scholar] [CrossRef]
  44. Hachana, O.; El-Fergany, A.A. Efficient PEM fuel cells parameters identification using hybrid artificial bee colony differential evolution optimizer. Energy 2022, 250, 123830. [Google Scholar] [CrossRef]
  45. Sun, X.; Wang, G.; Xu, L.; Yuan, H.; Yousefi, N. Optimal estimation of the PEM fuel cells applying deep belief network optimized by improved archimedes optimization algorithm. Energy 2021, 237, 121532. [Google Scholar] [CrossRef]
  46. Fathy, A.; Babu, T.S.; Abdelkareem, M.A.; Rezk, H.; Yousri, D. Recent approach based heterogeneous comprehensive learning Archimedes optimization algorithm for identifying the optimal parameters of different fuel cells. Energy 2022, 248, 123587. [Google Scholar] [CrossRef]
  47. Ashraf, H.; Elkholy, M.M.; Abdellatif, S.O.; El-Fergany, A.A. Accurate emulation of steady-state and dynamic performances of PEM fuel cells using simplified models. Sci. Rep. 2023, 13, 19532. [Google Scholar] [CrossRef]
  48. Agwa, A.; El-Fergany, A.; Sarhan, G. Steady-state modeling of fuel cells based on atom search optimizer. Energies 2019, 12, 1884. [Google Scholar] [CrossRef]
  49. Seleem, S.I.; Hasanien, H.M.; El-Fergany, A.A. Equilibrium optimizer for parameter extraction of a fuel cell dynamic model. Renew. Energy 2021, 169, 117–128. [Google Scholar] [CrossRef]
  50. Rezk, H.; Ferahtia, S.; Djeroui, A.; Chouder, A.; Houari, A.; Machmoum, M.; Abdelkareem, M.A. Optimal parameter estimation strategy of PEM fuel cell using gradient-based optimizer. Energy 2022, 239, 122096. [Google Scholar] [CrossRef]
  51. Rao, Y.; Shao, Z.; Ahangarnejad, A.H.; Gholamalizadeh, E.; Sobhani, B. Shark Smell Optimizer applied to identify the optimal parameters of the proton exchange membrane fuel cell model. Energy Convers. Manag. 2019, 182, 1–8. [Google Scholar] [CrossRef]
  52. Syah, R.; Guerrero, J.W.G.; Poltarykhin, A.L.; Suksatan, W.; Aravindhan, A.; Bokov, D.O.; Abdelbasset, W.K.; Al-Janabi, S.; Alkaim, A.F.; Tumanov, D.Y. Developed teamwork optimizer for model parameter estimation of the proton exchange membrane fuel cell. Energy Rep. 2022, 8, 10776–10785. [Google Scholar] [CrossRef]
  53. Menesy, A.S.; Sultan, H.M.; Korashy, A.; Kamel, S.; Jurado, F. A modified farmland fertility optimizer for parameters estimation of fuel cell models. Neural Comput. Appl. 2021, 33, 12169–12190. [Google Scholar] [CrossRef]
  54. Ashraf, H.; Abdellatif, S.O.; Elkholy, M.M.; El-Fergany, A.A. Honey badger optimizer for extracting the ungiven parameters of PEMFC model: Steady-state assessment. Energy Convers. Manag. 2022, 258, 115521. [Google Scholar] [CrossRef]
  55. Ćalasan, M.; Micev, M.; Hasanien, H.M.; Aleem, S.H.E.A. PEM fuel cells: Two novel approaches for mathematical modeling and parameter estimation. Energy 2024, 290, 130130. [Google Scholar] [CrossRef]
  56. Shaheen, A.M.; Alassaf, A.I.; Alsaleh, I.; El-Fergany, A.A. Enhancing model characterization of PEM Fuel cells with human memory optimizer including sensitivity and uncertainty analysis. Ain Shams Eng. J. 2024, 103026. [Google Scholar] [CrossRef]
  57. Ali, H.H.; Fathy, A. Reliable exponential distribution optimizer-based methodology for modeling proton exchange membrane fuel cells at different conditions. Energy 2024, 292, 130600. [Google Scholar] [CrossRef]
  58. Joyce, T.; Herrmann, J.M. A review of no free lunch theorems, and their implications for metaheuristic optimisation. Nat.-Inspired Algorithms Appl. Optim. 2018, 744, 27–51. [Google Scholar]
  59. Fu, S.; Li, K.; Huang, H.; Ma, C.; Fan, Q.; Zhu, Y. Red-billed blue magpie optimizer: A novel metaheuristic algorithm for 2D/3D UAV path planning and engineering design problems. Artif. Intell. Rev. 2024, 57, 1–89. [Google Scholar] [CrossRef]
  60. Zhao, S.; Zhang, T.; Ma, S.; Chen, M. Dandelion Optimizer: A nature-inspired metaheuristic algorithm for engineering applications. Eng. Appl. Artif. Intell. 2022, 114, 105075. [Google Scholar] [CrossRef]
  61. Bai, J.; Li, Y.; Zheng, M.; Khatir, S.; Benaisa, B.; Abualigah, L.; Wahab, M.A. A Sinh Cosh Optimizer. Knowl.-Based Syst. 2023, 82, 111081. [Google Scholar] [CrossRef]
  62. Zhang, Q.; Gao, H.; Zhan, Z.-H.; Li, J.; Zhang, H. Growth Optimizer: A powerful metaheuristic algorithm for solving continuous and discrete global optimization problems. Knowl.-Based Syst. 2023, 261, 110206. [Google Scholar] [CrossRef]
  63. Abbassi, R.; Saidi, S.; Abbassi, A.; Jerbi, H.; Kchaou, M.; Alhasnawi, B.N. Accurate Key Parameters Estimation of PEMFCs’ Models Based on Dandelion Optimization Algorithm. Mathematics 2023, 11, 1298. [Google Scholar] [CrossRef]
  64. Elhammoudy, A.; Elyaqouti, M.; Hmamou, D.B.; Lidaighbi, S.; Saadaoui, D.; Choulli, I.; Abazine, I. Dandelion Optimizer algorithm-based method for accurate photovoltaic model parameter identification. Energy Convers. Manag. X 2023, 19, 100405. [Google Scholar] [CrossRef]
  65. Aribowo, W.; Suprianto, B.; Prapanca, A. A novel modified dandelion optimizer with application in power system stabilizer. Int. J. Artif. Intell. 2023, 12, 2033–2041. [Google Scholar] [CrossRef]
  66. Sajid, I.; Gautam, A.; Sarwar, A.; Tariq, M.; Liu, H.-D.; Ahmad, S.; Lin, C.H.; Sayed, A.E. Optimizing Photovoltaic Power Production in Partial Shading Conditions Using Dandelion Optimizer (DO)-Based MPPT Method. Processes 2023, 11, 2493. [Google Scholar] [CrossRef]
  67. Nguyen, H.D.; Pham, L.H. Solutions of economic load dispatch problems for hybrid power plants using Dandelion optimizer. Bull. Electr. Eng. Inform. 2023, 12, 2569–2576. [Google Scholar] [CrossRef]
  68. Aribia, H.B.; El-Rifaie, A.M.; Tolba, M.A.; Shaheen, A.; Moustafa, G.; Elsayed, F.; Elshahed, M. Growth optimizer for parameter identification of solar photovoltaic cells and modules. Sustainability 2023, 15, 7896. [Google Scholar] [CrossRef]
  69. Chen, L.; Qiu, Z.; Wu, Y.; Tang, Z. Optimizing k-coverage in energy-saving wireless sensor networks based on the Elite Global Growth Optimizer. Expert Syst. Appl. 2024, 256, 124878. [Google Scholar] [CrossRef]
  70. Wang, X.; Wei, Y.; Guo, Z.; Wang, J.; Yu, H.; Hu, B. A Sinh–Cosh-Enhanced DBO Algorithm Applied to Global Optimization Problems. Biomimetics 2024, 9, 271. [Google Scholar] [CrossRef]
  71. Abualigah, L.; Ekinci, S.; Izci, D. Aircraft Pitch Control via Filtered Proportional-Integral-Derivative Controller Design Using Sinh Cosh Optimizer. Int. J. Robot. Control Syst. 2024, 4, 746–757. [Google Scholar]
  72. Emam, M.M.; Houssein, E.H.; Samee, N.A.; Alkhalifa, A.K.; Hosney, M.E. Optimizing cancer diagnosis: A hybrid approach of genetic operators and Sinh Cosh Optimizer for tumor identification and feature gene selection. Comput. Biol. Med. 2024, 180, 108984. [Google Scholar] [CrossRef] [PubMed]
  73. Izci, D.; Rizk-Allah, R.M.; Snášel, V.; Ekinci, S.; Migdady, H.; Daoud, M.S.; Altalhi, M.; Abualigah, L. Refined sinh cosh optimizer tuned controller design for enhanced stability of automatic voltage regulation. Electr. Eng. 2024, in press. [Google Scholar] [CrossRef]
  74. Ibrahim, R.A.; Aly, M.A.S.; Moemen, Y.S.; El Sayed, I.E.T.; Elaziz, M.A.; Khalil, H.A. Boosting Sinh Cosh Optimizer and arithmetic optimization algorithm for improved prediction of biological activities for indoloquinoline derivatives. Chemosphere 2024, 359, 142362. [Google Scholar] [CrossRef]
  75. Mann, R.F.; Amphlett, J.C.; Hooper, M.A.; Jensen, H.M.; Peppley, B.A.; Roberge, P.R. Development and application of a generalised steady-state electrochemical model for a PEM fuel cell. J. Power Sources 2000, 86, 173–180. [Google Scholar] [CrossRef]
  76. Fawzi, M.; El-Fergany, A.A.; Hasanien, H.M. Effective methodology based on neural network optimizer for extracting model parameters of PEM fuel cells. Int. J. Energy Res. 2019, 43, 8136–8147. [Google Scholar] [CrossRef]
  77. Sultan, H.M.; Menesy, A.S.; Hassan, M.; Jurado, F.; Kamel, S. Standard and Quasi Oppositional bonobo optimizers for parameter extraction of PEM fuel cell stacks. Fuel 2023, 340, 127586. [Google Scholar] [CrossRef]
  78. Menesy, A.S.; Sultan, H.M.; Korashy, A.; Banakhr, F.A.; Ashmawy, M.G.; Kamel, S. Effective parameter extraction of different polymer electrolyte membrane fuel cell stack models using a modified artificial ecosystem optimization algorithm. IEEE Access 2020, 8, 31892–31909. [Google Scholar] [CrossRef]
  79. Zhu, Y.; Yousefi, N. Optimal parameter identification of PEMFC stacks using adaptive sparrow search algorithm. Int. J. Hydrogen Energy 2021, 46, 9541–9552. [Google Scholar] [CrossRef]
  80. Herrera, P.A.; Marazuela, M.A.; Hofmann, T. Parameter estimation and uncertainty analysis in hydrological modeling. Wiley Interdiscip. Rev. Water 2022, 9, e1569. [Google Scholar] [CrossRef]
  81. Al, R.; Behera, C.R.; Zubov, A.; Gernaey, K.V.; Sin, G. Meta-modeling based efficient global sensitivity analysis for wastewater treatment plants–An application to the BSM2 model. Comput. Chem. Eng. 2019, 127, 233–246. [Google Scholar] [CrossRef]
Figure 1. Procedures of the RBMO framework.
Figure 1. Procedures of the RBMO framework.
Technologies 12 00156 g001
Figure 2. Convergence patterns of all studied test cases.
Figure 2. Convergence patterns of all studied test cases.
Technologies 12 00156 g002
Figure 3. The principal performance of Ballard Mark V.
Figure 3. The principal performance of Ballard Mark V.
Technologies 12 00156 g003
Figure 4. Principal performance of Temasek 1 kW.
Figure 4. Principal performance of Temasek 1 kW.
Technologies 12 00156 g004
Figure 5. Principal performance of Horizon H-12 unit.
Figure 5. Principal performance of Horizon H-12 unit.
Technologies 12 00156 g005
Figure 6. Percentage of biased voltage.
Figure 6. Percentage of biased voltage.
Technologies 12 00156 g006aTechnologies 12 00156 g006b
Table 1. Technical specifications of various commercial PEMFCs units.
Table 1. Technical specifications of various commercial PEMFCs units.
PEMFCs’ UnitSpecs
n ,   c e l l s A m ,   c m 2 l ,   μ m J m a x ,   m A / c m 2 T ,   K e l v i n P H 2 ,   a t m P O 2 ,   a t m
Ballard Mark V 3550.617815003431.01.0
Temasek 1 kW201505115003230.50.5
Horizon H-12138.125246.9302.150.49351.0
Table 2. Limits of the PEMFC’s 7 unknown parameters.
Table 2. Limits of the PEMFC’s 7 unknown parameters.
Limits ξ 1 ( V ) ξ 2   ×   10 3   ( V / K ) ξ 3   ×   10 5 (V/K) ξ 4   ×   10 5 (V/K) R C   ( m Ω )   λ β (V)
L B −1.19970.80003.6000−26.00000.100013.00000.0136
H B −0.85326.00009.8000−9.54000.800023.00000.5000
Table 3. RBMO’s performance in comparison to other recent optimizers for Ballard Mark V.
Table 3. RBMO’s performance in comparison to other recent optimizers for Ballard Mark V.
Parameter ε 1   ( V ) ε 2 . 10 3   ( V / K ) ε 3 . 10 5   ( V / K ) ε 4 . 10 5   ( V / K ) λ R c   ( m Ω ) β   ( V ) S S D   ( V 2 )
Optimizer
CSA [35]1.181303.5690963.9929−16.283023.000.10000.01360.853601
MRFO [41]1.089803.82497.7306−16.283023.000.10000.01360.853661
NNO [76]0.979973.69409.0870−16.280023.000.10000.01360.863697
WOA [41]1.1978004.41839.7214−16.273023.000.10020.01360.853766
GhO [40]0.853003.41709.8000−15.950022.840.1000 0.01360.853661
ETSO [36]−0.853402.55913.6100−16.287023.000.10000.01360.853600
ARO [3]−1.1588593.52084.0526−16.725123.99 *0.10000.0158840.813912
IAHA [4]−1.01304.00008.9800−16.300023.00000.10000.01360.853608
HBO [54]−1.19974.334539.20688−16.28323.00000.10000.01360.853608
TSO [54]−0.85522.722274.86143−16.283123.00000.10000.01360.853608
ESMO [54]−0.85322.540553.60422−16.282423.00000.10000.01360.853608
GO−1.192633.525573.61898−16.0029622.999670.124420.013680.865333
SCHO−1.193923.982117.07842−14.8057022.420480.121540.014630.957596
DOA−0.857452.932646.31831−16.2830423.000000.100000.013600.853608
RBMO−1.0006183.836439.79532−16.2829622.999990.100000.013600.853608
* Out of the range of λ (13 to 23), as stated in Table 2; thus, it is unfeasible answer for SSD’s value.
Table 4. RBMO’s performance in comparison to other recent optimizers for Temasek 1 kW.
Table 4. RBMO’s performance in comparison to other recent optimizers for Temasek 1 kW.
Parameter ε 1   ( V ) ε 2 . 10 3   ( V / K ) ε 3 . 10 5   ( V / K ) ε 4 . 10 5   ( V / K ) λ R c   ( m Ω ) β   ( V ) S S D   ( V 2 )
Optimizer
SSO [51]−1.02992.41054.00000−9.540010.00050.108700.12741.6481
CHHO [33]−1.09444.42828.76560−21.465018.63920.189100.10160.80234
MFFO [53]−0.90353.82678.47510−22.934713.32510.100100.07050.791000
QOBO [77]−1.19973.82203.60000−22.950013.00000.100000.06800.783040
MAEA [78]−0.85443.57667.88880−22.925813.00170.100000.06830.79096
MAEO [78]−1.117063.82904.56767−2.2630922.52320.101920.110190.79243
HHO [78]−0.853203.499107.21118−2.4404922.99630.182650.071240.80553
GO−0.914023.416439.73056−9.5400013.032310.100760.161170.597749
SCHO−0.873442.801626.26615−10.0982718.000000.100000.185410.731239
DOA−1.168733.612335.66745−9.5400022.296250.100000.199620.601922
RBMO−0.919213.331839.04072−9.5400013.000000.100000.161050.597504
Table 5. RBMO’s performance in comparison to other recent optimizers for H-12.
Table 5. RBMO’s performance in comparison to other recent optimizers for H-12.
Parameter ε 1   ( V ) ε 2 . 10 3   ( V / K ) ε 3 . 10 5   ( V / K ) ε 4 . 10 5   ( V / K ) λ R c   ( m Ω ) β ( V ) S S D   ( V 2 )
Optimizer
WOA [41]−1.18702.669703.6000−9.540013.82400.80000.15980.1160
MRFO [43]−1.06302.364104.3272−9.540019.81500.28530.18290.0966
ASSA [79]−1.13002.440003.5700−9.540018.79000.714018.17000.0970
PFA [23]−1.11132.057303.6000−9.540022.99990.10580.18680.0965
ETSO [36]−1.0322852.7296777.7200−9.540022.998980.1124170.1868880.09653
TSO [36]−10.85321.5718523.6100−9.540013.024370.3278740.1752740.09685
HHO [36]−0.8614881.932106.0300−9.540013.542380.19740.1740670.09657
PSO [36]−1.0347542.544906.3200−9.540023.00000.80000.1827040.09658
ABDEO [42]−0.8543520.096136.76763−9.5400023.000000.100000.186850.096536
AEO [42]−1.0027827.449078.53668−9.5400018.139000.100000.182300.096536
SCE [42]−0.8557515.775363.60145−9.5400023.000000.100000.186850.096536
ELBO [42]−0.9771223.737876.46218−9.5400023.000000.100000.186850.096536
GO−1.186683.381348.80335−9.5400022.577000.101860.035100.0616502
SCHO−1.048622.503215.75161−9.7594314.805010.466640.031530.0672907
DOA−0.864281.757984.76938−9.5400021.809410.131550.0347970.601255
RBMO−0.977381.982963.68744−9.5400022.999700.100000.035050.0616371
Table 6. Detailed dataset points of actual and estimated voltage points determined by the RBMO.
Table 6. Detailed dataset points of actual and estimated voltage points determined by the RBMO.
Ballard Mark VTemasek 1 kWHorizon H-12
I m , j   A V m , j   V V e , j   V V B E , j I m , j   A V m , j   V V e , j   V V B E , j I m , j   A V m , j   V V e , j   V V B E , j
5.060033.250032.93920.31080.832218.251718.5784−0.32670.10409.58009.7198−0.1398
10.626030.800031.0698−0.26983.744817.522217.5809−0.05870.20009.42009.4477−0.0277
16.192029.750029.8076−0.05766.380017.149317.1880−0.03870.30909.25009.2524−0.0024
20.240028.700029.0379−0.33799.570016.857516.8590−0.00150.40309.20009.12350.0765
27.830028.000027.75030.249714.147016.468416.5030−0.03450.51009.09008.99990.0901
34.408026.600026.7082−0.108219.001416.192816.1966−0.00380.61408.95008.89380.0562
37.444026.250026.23400.016024.271815.884815.9076−0.02280.70308.85008.80970.0403
43.010025.200025.3594−0.159428.987515.690315.67200.01830.80608.74008.71710.0229
48.070024.500024.5439−0.043933.148415.495715.47570.02000.90808.65008.62810.0219
56.166023.800023.15830.641738.002815.317415.25630.06111.07608.45008.4823−0.0323
61.226022.050022.2109−0.160943.689315.139115.00800.13111.12708.41008.4373−0.0273
67.298021.000020.93020.069848.821114.944514.78920.15541.28808.20008.2885−0.0885
71.852019.600019.7503−0.150353.675414.782414.58500.19741.39008.12008.1858−0.0658
S u m V B E , V0.198258.529814.587914.38220.20571.45008.11008.1206−0.0106
63.384214.377114.17980.19731.57808.05007.96350.0865
68.654714.101513.95960.1420 S u m V B E , V0.0526
73.231613.842213.76710.0750
78.086013.534113.5612−0.0271
81.692113.226113.4067−0.1806
S u m V B E , V0.1095
Table 7. Sensitivity indices based on GPR and ANN.
Table 7. Sensitivity indices based on GPR and ANN.
PEMFCs UnitIndex ε 1   ( V ) ε 2 . 10 3   ( V / K ) ε 3 . 10 5   ( V / K ) ε 4 . 10 5   ( V / K ) λ R c   ( m Ω ) β   ( V )
Ballard Mark 5 g p r S i −0.020110.14612−0.04515−0.02258−0.02196−0.02004−0.01996
g p r S T i 0.736940.894280.127010.053360.016400.000020.00023
Rank2134576
a n n S i 0.053740.199680.072010.096200.092910.086870.08374
a n n S T i 0.807680.912630.139910.053260.016620.000020.00028
Rank2134576
Temasek 1 kW g p r S i −0.011860.12336−0.04709−0.02796−0.02626−0.02619−0.02688
g p r S T i 0.724540.870720.206210.013150.000560.000020.00143
Rank2134675
a n n S i 0.046190.204950.076470.084170.087680.086730.07861
a n n S T i 0.791040.885830.226180.012960.000560.000020.00164
Rank2134675
Horizon H−12 g p r S i 0.254050.03101−0.04507−0.03886−0.03850−0.03864−0.03890
g p r S T i 0.979690.678890.051940.000340.000010.000000.00107
Rank1235674
a n n S i 0.291410.078410.107730.103320.104580.104520.09923
a n n S T i 1.001660.722350.057220.000320.000010.000000.00121
Rank1235674
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

El-Fergany, A.A.; Agwa, A.M. Red-Billed Blue Magpie Optimizer for Electrical Characterization of Fuel Cells with Prioritizing Estimated Parameters. Technologies 2024, 12, 156. https://doi.org/10.3390/technologies12090156

AMA Style

El-Fergany AA, Agwa AM. Red-Billed Blue Magpie Optimizer for Electrical Characterization of Fuel Cells with Prioritizing Estimated Parameters. Technologies. 2024; 12(9):156. https://doi.org/10.3390/technologies12090156

Chicago/Turabian Style

El-Fergany, Attia A., and Ahmed M. Agwa. 2024. "Red-Billed Blue Magpie Optimizer for Electrical Characterization of Fuel Cells with Prioritizing Estimated Parameters" Technologies 12, no. 9: 156. https://doi.org/10.3390/technologies12090156

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop