[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Next Article in Journal
Theoretical Identification of Coupling Effect and Performance Analysis of Single-Source Direct Sampling Method
Next Article in Special Issue
An Optimized Triggering Algorithm for Event-Triggered Control of Networked Control Systems
Previous Article in Journal
2DOF IMC and Smith-Predictor-Based Control for Stabilised Unstable First Order Time Delayed Plants
Previous Article in Special Issue
Dynamic Stability Performance of Autonomous Microgrid Involving High Penetration Level of Constant Power Loads
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Mathematics
Volume 9
Issue 9
10.3390/math9091066
mathematics-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Parameter Estimation Methodology of Solid Oxide Fuel Cell Using Modern Optimization

by
Hesham Alhumade
1,2,*,
Ahmed Fathy
3,4,
Abdulrahim Al-Zahrani
1,
Muhyaddin Jamal Rawa
2,5 and
Hegazy Rezk
6,7
1
Department of Chemical and Materials Engineering, Faculty of Engineering, King Abdulaziz University, Jeddah 21589, Saudi Arabia
2
Center of Research Excellence in Renewable Energy and Power systems, King Abdulaziz University, Jeddah 21589, Saudi Arabia
3
Electrical Engineering Department, Faculty of Engineering, Jouf University, Sakaka 72388, Saudi Arabia
4
Electrical Power and Machine Department, Faculty of Engineering, Zagazig University, Zagazig 44519, Egypt
5
Department of Electrical and Computer Engineering, Faculty of Engineering, King Abdulaziz University, Jeddah 21589, Saudi Arabia
6
College of Engineering at Wadi Addawaser, Prince Sattam Bin Abdulaziz University, Al-Kharj 11911, Saudi Arabia
7
Electrical Engineering Department, Faculty of Engineering, Minia University, Minia 61517, Egypt
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(9), 1066; https://doi.org/10.3390/math9091066
Submission received: 27 March 2021 / Revised: 21 April 2021 / Accepted: 4 May 2021 / Published: 10 May 2021
(This article belongs to the Special Issue Advanced Modeling, Control, and Optimization Methods in Power Hybrid Systems - 2021)
Browse Figures
Figure 1
<p>Optimization process of EO.</p> "> Figure 2
<p>The proposed steps incorporating EO.</p> "> Figure 3
<p>The measured and calculated polarization curves of SOFC operated at 1073 K obtained via EO at (<b>a</b>) current density-voltage, (<b>b</b>) current density-power.</p> "> Figure 4
<p>(<b>a</b>) Current-voltage curve, (<b>b</b>) Current-power curve of SOFC operated at 1073 K obtained via other approaches.</p> "> Figure 5
<p>The variation of fitness function during iterative process for all employed optimizers applied for steady-state SOFC model.</p> "> Figure 6
<p>The measured and calculated (<b>a</b>) current density-voltage, (<b>b</b>) current density-power of SOFC operated at 1173 K, 1213 K, and 1273 K obtained via the proposed EO.</p> "> Figure 7
<p>The measured and calculated polarization curves of SOFC dynamic-state model operated at 1273 K obtained via EO (<b>a</b>) current-voltage, (<b>b</b>) current-power.</p> "> Figure 8
<p>The measured and calculated polarization curves of SOFC dynamic-state model operated at 1273 K obtained via MPA, HBO, SOA, and MRFO.</p> "> Figure 9
<p>The variation of fitness function during iterative process for all employed optimizers applied for dynamic-state SOFC model.</p> "> Figure 10
<p>First load disturbance applied on SOFC stack dynamic model.</p> "> Figure 11
<p>Second load disturbance applied on SOFC stack dynamic model.</p> ">
Versions Notes

Abstract

:
An optimal parameter estimation methodology of solid oxide fuel cell (SOFC) using modern optimization is proposed in this paper. An equilibrium optimizer (EO) has been used to identify the unidentified parameters of the SOFC equivalent circuit with the assistance of experimental results. This is presented via formulating the modeling process as an optimization problem considering the sum mean squared error (SMSE) between the observed and computed voltages as the target. Two modes of the SOFC-based model are investigated under variable operating conditions, namely, the steady-state and the dynamic-state based models. The proposed EO results are compared to those obtained via the Archimedes optimization algorithm (AOA), Heap-based optimizer (HBO), Seagull Optimization Algorithm (SOA), Student Psychology Based Optimization Algorithm (SPBO), Marine predator algorithm (MPA), Manta ray foraging optimization (MRFO), and comprehensive learning dynamic multi-swarm marine predators algorithm. The minimum fitness function at the steady-state model is obtained via the proposed EO with value of 1.5527 × 10−6 at 1173 K. In the dynamic based model, the minimum SMSE is 1.0406. The obtained results confirmed the reliability and superiority of the proposed EO in constructing a reliable model of SOFC.

1. Introduction

There is a growing demand for energy to meet the requirements of continuous industrial development and modern civilization. In parallel, there is a growing concern about the depletion of traditional energy sources such as fossil fuel and drawbacks of continuous consumption of fossil fuel such as climate change [1]. Indeed, a recent study expected that future energy demands might exceed the limits of current energy systems [2]. Moreover, the increasing global energy demands and consumption of fossil fuel will escalate the emissions of greenhouse gases and other toxic air pollutants. Therefore, alternative sources of energy such as renewable energy have earned significant attention in the recent decades. In particular, fuel cell is among the power generation systems that can deliver environmentally friendly quality energy with great energy conversion efficiency. Furthermore, fuel cell has a great potential in power delivery for stationary and movable applications compared to other storage technologies [3,4,5,6,7]. Other remarkable features of fuel cell over other energy alternatives include lower fuel oxidation temperature and reduced emissions [8]. Solid oxide fuel cell (SOFC) and polymer composites-based electrolyte fuel cell represent the most attractive types of fuel cell for a wide range of applications. A growing effort is made to deliver a model that can predict the performance of fuel cell over steady state or dynamics operating environmental conditions [9,10,11,12]. In fact, appropriate model identification requires feeding accurate input parameters to the governing equations that encompass physical and chemical properties of the cell, where the modelling methodologies can be empirical, semi-empirical, or theoretical [13,14,15,16].
The parameter extraction of the fuel cell model plays an important role in the simulation, evaluation, control, and optimization of a fuel cell system. The voltage drops in SOFC are mainly reliant on the parameters associated with the chemical processes inside SOFC [17,18]. Several methods were used to identify the accurate parameters of SOFC. Among these methods, the metaheuristic optimization-based methodologies were superior in resolving the SOFC parameter estimation problem due to their reliability, robustness, and simplicity. Shi et al. [19] proposed a strategy, Converged Grass Fibrous Root, to determine the best parameters of the SOFC model. Both temperature and pressure variation are considered. During the optimization process, seven parameters are assigned to be decision variables: the standard potential, the current limitation density, the Tafel line slope, a constant depends on the operating state of SOFC, the area-specific resistance, the anode exchange current density, and the cathode exchange current density. El-Hay et al. [10] suggested a methodology based on an interior search algorithm to estimate the steady state and transient parameters of SOFC. A proportional-integral controller is integrated with the dynamic model to enhance its performance throughout transient disturbances. A similar study also carried by the same authors based on Satin Bowerbird Optimizer was conducted [9]. During the optimization process, the decision variables are represented by the unknown parameters of SOFC, whereas the cost function is represented by mean squared deviations between experimental data and estimated SOFC voltages. In the same direction, Yousri et al. [20] proposed a modified algorithm called comprehensive learning dynamic multi-swarm marine predators to determine both static and dynamic parameters of SOFC. During the optimization process, the mean squared error between the experimental data and estimated SOFC voltage is used as the objective function that is required to be minimum. Nassef et al. [21] used the radial movement optimization algorithm (RMOA) to determine the best parameters of the SOFC model. The model of SOFC was created using a neural network. During the optimization process, four parameters, including electrolyte thickness, cathode interlayer thickness, anode porosity, and anode support layer thickness, are used as decision variables; in contrast, the objective function is represented by the power density of SOFC. By using the RMOA, the power density was increased by 17.28% compared with the genetic algorithm. In the same direction, Fathy et al. [22] suggested a methodology based on the moth-flame optimization algorithm (MFOA). The power density of SOFC using the MFOA was improved. It was increased by 18.92% and 5.56% compared to the genetic algorithm and RMOA, respectively.
The contribution of the current research work can be summarized as follows:
  • A novel approach based on Equilibrium Optimizer (EO) is suggested to determine the optimal parameters of the SOFC-based model.
  • The suggested methodology is validated through both steady-state and dynamic-state models of SOFC with the changing of the operational conditions.
  • A comprehensive comparison with previous works and other programs of the Archimedes optimization algorithm (AOA), Heap-based optimizer (HBO), Seagull Optimization Algorithm (SOA), Student Psychology Based Optimization Algorithm (SPBO), Marine predator algorithm (MPA), and Manta ray foraging optimization (MRFO).
  • The superiority and reliability of the suggested EO-based strategy in solving the SOFC parameter determination problem is verified.
The rest of the paper is organized as follows: The mathematical model of SOFC is illustrated in Section 2. Section 3 presents an overview about main aspects of the equilibrium optimizer. Then, the suggested optimization problem and solution methodology are explained in Section 4. Section 5 presents a detailed discussion of the obtained results and a comparative study with other methods. Finally, the main findings and the future work are outlined in Section 6.

2. SOFC Mathematical Model

In this section, the authors illustrate the static- and dynamic-based models of SOFC.

2.1. Steady-State Model

The output voltage of SOFC can be estimated using the following relation considering the activation loss, ohmic loss, and concentration loss [23]:
V c e l l = E n − V a − V o − V c
where En denotes the reversible voltage of the cell, and Va, Vo, and Vc denote the activation, ohm, and concentration voltage drops, respectively.
The activation loss, concentration loss and ohmic loss can be estimated using the following relations:
V a = a ⋅ ln J 2 J 0
V o = ∑ k i × R k
V c = − b ⋅ ln 1 − J J max
where a and b are constants; J denotes the current density; J0 denotes the exchange current density; Jmax is the maximum current density; Rk is the sum of ionic (electrolyte) and electronic resistances; and i denotes the output SOFC current.
To increase the rating voltage of SOFC, the number of cells are connected in series. Therefore, the total stack output voltage can be estimated using the following relation.
V s = n c × V c e l l = n c × E n − V a − V o − V c
where vs. is the stack voltage and nc is the number of cells. The reversible voltage can be written as follows [24]:
E n = E 0 + R T 2 F ln P H 2 P O 2 P H 2 O
where E0 denotes the reference voltage at unit activity and atmospheric pressure; T denotes operating temperature (K); PH2, PO2, and PH2O denote the hydrogen, oxygen, and water partial pressures, respectively; R is the universal gas constant with a value of 8.314 kJ (kmol K)−1; and F is the Faraday constant.

2.2. SOFC Dynamic Model

The gas molar flow in SOFC is reliant on hydrogen and oxygen partial pressures as follows [24]:
q H 2 P H 2 = k a n M H 2 = K H 2
q O 2 P O 2 = k a n M O 2 = K O 2
where qH2 is the molar flow of hydrogen; qO2 is the molar flow of oxygen; KH2 is the hydrogen molar constant; KO2 is the oxygen molar constant; kan is the anode valve constant; MH2 is the molar masses of hydrogen; and MO2 is the molar masses of oxygen.
The partial pressure derivative can be estimated using the following relation:
d P H 2 d t = R T V a n q H 2 i n − q H 2 o u t − q H 2 r
where Van denotes the anode volume; qH2in is the hydrogen input flow rate; qH2out is the hydrogen output flow rate; and qH2r denotes the reacted hydrogen flow rate.
The reacted hydrogen flow rate can be estimated based on the following relation:
q H 2 r = n c i 2 F = 2 K r i
where Kr is constant.
Considering the above equations and the Laplace transform, hydrogen and oxygen partial pressures can be formulated as follows.
P H 2 ( s ) = 1 K H 2 1 + τ H 2 q H 2 i n − 2 K r i
P H 2 O ( s ) = 1 K H 2 O 1 + τ H 2 O 2 K r i
P O 2 ( s ) = 1 K O 2 1 + τ O 2 q O 2 i n − K r i
where τH2, τH2O, and τO2 denote the flow time constants of hydrogen, water, and oxygen, respectively.
Ultimately, the dynamic model of SOFC voltage is given by the following relation.
V s = n c E 0 + R T 2 F ln P H 2 P O 2 P H 2 − a ⋅ ln J 2 J 0 + r × i − b ⋅ ln 1 − J J max

3. Overview of Equilibrium Optimizer

Equilibrium optimizer (EO) is a recent algorithm that was proposed by Faramarzi et al. in 2020 [25]. The core idea of the EO is extracted from the control volume mass balance models. During the optimization process of the EO, the particles and positions are assigned to the solutions and concentrations, respectively. The details explanations about the inspiration, mathematical model, and algorithm of the EO can be found in [25]. The mass balance formula is expressed as follows.
V d C d t = Q C e q − Q C − G
where C denotes the concentration of the control volume; V d C d t denotes the changing rate of the mass; Q denotes the flow rate; Ceq denotes the concentration at the balance state; and G denotes the mass generation rate.
By integration over time and rearranged with the above formula, the following relation can be used to express the concentration of the control volume.
C = C e q + ( C 0 − C e q ) F + G λ V ( 1 − f )
where λ denotes the turnover rate ( λ = Q V ) and f = e − λ ( t − t 0 ) , and t0 and C0 denote the initial time and concentration, respectively.
The time t is decreasing with increasing the number of iterations as follows:
t = ( 1 − I z ) a 2 I z
t 0 → = 1 λ → ln ( − a 1 s i g n ( r → − 0.5 ) [ 1 − e λ t → ] ) + t
F → = a 1 s i g n ( r → − 0.5 ) [ 1 − e λ t → ]
where I and z are the current iteration and maximum number of iterations, respectively. a1 and a2 are constants. r → is a random vector in range [0, 1].
Considering Equation (16), there are three sections describing the updating process for particles. The first section represents the equilibrium concentration. It represents the optimum solutions arbitrarily chosen from a pool. The second section is related to the concentration variations between a particle and the equilibrium state. The last section is related to the generation rate. It is mainly performing the role of an exploiter. The equilibrium state is the final convergence state of the EO optimization process. The equilibrium pool can be represented as follows.
C → e q . p o o l = C → e q ( 1 ) , C → e q ( 2 ) , C → e q ( 3 ) , C → e q ( 4 ) , C → e q ( a v e )
During the first iteration, the particle modifies the concentration using C → e q ( 1 ) whereas it is uses C → e q ( a v e ) with other iterations. The generation rate is defined as follows.
G → = G 0 → e − k → ( t − t 0 )
where G0 is the initial value, and k denotes a decay constant (k = λ).
G → = G 0 → e − λ → ( t − t 0 ) = G 0 → F →
G 0 → = G C P → ( C → e q − λ → C → )
G C P → = 0.5 r 1 r 2 ≥ G P 0 r 2 < G P
where r1 and r2 are random variables in range [0, 1], and GCP is a parameter that controls the generation rate.
The addition of memory saving helps each particle to save its coordinates in the search space. Moreover, it informs its fitness function value. The fitness function related to a particular particle in the ongoing iteration is compared with the previous one; then, the updating process is placed if it reaches better fit. This action enhances the exploitation phase. The optimization process of EO is illustrated in Figure 1.

4. The Proposed Methodology

In this section, the proposed methodology is explained via formulating the problem of SOFC parameter estimation as an optimization problem, and an objective function and corresponding constraints are also introduced. Moreover, the proposed solution methodology is also presented.

4.1. The Proposed Objective Function

In this section, the formulation of the identification process of SOFC circuit’s parameters as an optimization problem is explained via specifying the fitness function, the corresponding constraints, and the proposed methodology of the solution incorporated by the EO. The main target of the work is to construct a reliable equivalent circuit of SOFC by identifying six parameters, E0, a, Jo, r, b, and Jmax; this is achieved with the aid of experimental data of the current and voltage of FC. The fitness function represented in this work is the sum mean squared error between the measured and calculated terminal voltages of SOFC, and this can be described as follows:
M i n i m i z e S M S E = ∑ k = 1 N 1 N V m e a s , k − V c a l , k 2
where Vmeas,k and Vcal,k are the kth observed and computed voltages, respectively, and N is the number of measured datasets. The constraints related to the variables to be designed can be described as follows:
E 0 min ≤ E 0 ≤ E 0 max a min ≤ a ≤ a max J o min ≤ J o ≤ J o max r min ≤ r ≤ r max b min ≤ b ≤ b max J max min ≤ x 6 ≤ J max max
where min denotes the minimum limit and max denotes the maximum limit. In this work the authors only considered the disturbance on the SOFC output voltage. However, in future work, they will consider the disturbance not only on the SOFC output but also on the SOFC input such as those given in [26,27,28,29].

4.2. The Proposed EO Based Methodology

The equilibrium optimizer is selected due to many advantages: It is simple in implementation; it achieves balance between the exploration and exploitation phases; and there is diversity between the population individuals. These features enable the algorithm to be applicable for many optimization problems. Six parameters are required to be identified such that the SMSE is minimized. The proposed methodology incorporating the EO begins by defining the specifications of SOFC and the recorded measured data of the terminal voltage. Then, an initial population with a dimension of npop × dim, npop is the population size, and dim is the problem dimension, which is constructed with the aid of the minimum and maximum limits defined by the user. The initial corresponding fitness function (SMSE) is calculated, and the iterative process is implemented by calculating the fitness function of each particle. The obtained fitness function is checked with those of the equilibrium pool to decide the updating action of each particle. After that, the condition of the last particle is investigated, and the average of the equilibrium pool is calculated, which helps in estimating the generation rate. The process is continued until the constraint of maximum iteration is achieved. At this moment, the optimal results are obtained and can be printed. The proposed methodology incorporating EO is shown in Figure 2.

5. Numerical Analysis

The analysis is performed on two modes of the SOFC operation which are steady-state and dynamic-state. Both of them are investigated under variable-operating conditions. The commercial SOFC, which is manufactured by Siemens [30], is employed in steady-state analysis. In such a case, four measured datasets are recorded at temperatures of 1073, 1173, 1213, and 1273 K, where the proposed EO size of population is assigned as 50, and the number of iterations is selected as 100. The population-based approach presented in this work has some difficulty, such as getting premature and local optima, and the authors take into consideration this problem by performing the approach with 50 independent runs, and the best one is selected as global optima. This action minimizes the problem of falling in local optima. Other metaheuristic approaches are implemented and compared to the proposed EO; these algorithms are the Archimedes optimization algorithm (AOA), Heap-based optimizer (HBO), Seagull Optimization Algorithm (SOA), Student Psychology Based Optimization Algorithm (SPBO), Marine predator algorithm (MPA), Manta ray foraging optimization (MRFO), and comprehensive learning dynamic multi-swarm marine predators algorithm CLDMMPA [20]. Table 1 shows the obtained optimal parameters of SOFC operated at 1073 K via the proposed EO and the others. The proposed approach succeeded in achieving a fitness function of 2.6906 × 10−6 which is the same obtained via CLDMMPA. However, the CLDMMPA is complex in construction; moreover, the proposed EO consumes only 272.198102 s, which is the best compared to the others. The measured and calculated polarization curves obtained via the proposed EO are shown in Figure 3. Both curves are closely converged. Moreover, Figure 4 shows the estimated polarization curves obtained via the other approaches and the measured ones. Furthermore, the performance of each optimizer during the iterative process is shown in Figure 5. It is clear that the EO performance is the best compared to the others.
The optimal parameters of the SOFC steady-state based model at 1173, 1213, and 1273 K obtained via the proposed EO and the others are tabulated in Table 2, Table 3 and Table 4. Regarding the obtained results at 1173 K, the best fitness function is 1.5527 × 10−6 obtained via the proposed EO while CLDMMPA comes in the second rank with a SMSE of 1.5529 × 10−6. On the other hand, the worst approach is SOA with a fitness function of 3.1657 × 10−6. Moreover, during the operation at 1213 K, the EO outperformed the others in terms of elapsed time and fitness function. The reader can see that during operation at 1273 K, the proposed EO achieved a SMSE of 2.2995 × 10−6, which is the best compared to the others.
The polarization curves of the measured data and calculated data obtained via the proposed EO for the steady-state SOFC based model at 1173, 1213, and 1273 K are shown in Figure 6. The curves confirm the matching between the experimental and calculated data.
It is important to investigate the performance of each optimizer via calculating the statistical parameters, which include the best, worst, mean, median, variance, and standard deviation after 50 independent runs. These data are calculated and tabulated in Table 5. The proposed EO gives acceptable statistical parameters compared to the others.
The obtained results confirmed the superiority and reliability of the proposed methodology incorporating EO in identifying the optimal parameters of the SOFC steady-state based model.
It is important to confirm the availability of the presented approach in a dynamic/transient-based model of SOFC. Therefore, a 100 kW stack with specifications given in Table 6 is modeled in a dynamic-state model subjected to variable load disturbances. At the beginning, the proposed EO is applied to identify the optimal parameters of a 100 kW SOFC stack operated at 1273 K; the obtained parameters are tabulated in Table 7 in comparison to those obtained by the others. Regarding the obtained results, the proposed EO outperformed the others, achieving the minimum SMSE with a value of 1.0406. MRFO comes in the second rank with a fitness function of 1.0775, and then CLDMMPA achieves an SMSE of 1.3204 and comes in the third rank. Figure 7 shows the measured and calculated polarization curves obtained via the EO, and both curves are closely converged. However, Figure 8 shows the polarization curves obtained via MPA, HBO, SOA, and MRFO. The statistical parameters of all optimizers in such cases are calculated and tabulated in Table 8, where the best parameters are obtained by the proposed EO. Figure 9 shows the performance of each optimizer during implementing the iterative process. The performance of the proposed EO is confirmed to be better than the others.
After identifying the parameters of the transient-state based model of SOFC, the model of fuel cell is implemented in Simulink/Matlab, and two load disturbances are applied to the model. The first disturbance is shown in Figure 10 (1st graph), the power is changed from 30 kW to 60 kW at a time of 300 sec., and given the identified parameters via the proposed EO and the stack output power plotted with the load disturbance, it is clear that they are closely matched, this means that the EO succeeded in extracting the correct parameters of the SOFC dynamic model. Moreover, the load current of the constructed model is closely converged to the disturbance current (3rd graph, Figure 10). Moreover, the terminal voltage (2nd graph) and the voltage drops occur inside the stack (4th graph), and they are shown in Figure 10. The constructed model via the proposed EO succeeded in tracking the changes in the load power. Moreover, a second disturbance is applied on the dynamic model in which the load power has two variations; as shown in Figure 11 (1st graph), the load power changes from 20 kW to 40 kW at 200 s and then changes again to 60 kW at 400 s on the same graph. The output power from the constructed model with identified parameters via the proposed EO is given, and both curves are converged. The terminal voltage, the corresponding current, and the voltage drops are shown in Figure 11 (2nd graph, 3rd graph, and 4th graph, respectively).
Finally, it can be concluded that the proposed methodology incorporating the EO is reliable, superior, and efficient over other employed approaches in constructing a reliable model of the SOFC-based model operated under either steady-state or dynamic-state modes.

6. Conclusions

This work introduces a new methodology based on a new metaheuristic approach named equilibrium optimizer (EO) to estimate the optimal parameters of a solid oxide fuel cell (SOFC) model. This is achieved with the aid of experimental datasets of the fuel cell polarization curves. The sum squared error difference between the cell experimental and computed voltages is selected as the fitness function to be minimized. The work investigates two operating modes of FC, which are steady- and dynamic-states models under altering operating conditions. In the first model, the parameters are estimated at four temperatures via the recorded measured polarization curves at them. In the dynamic model, two load power disturbances are investigated after identifying the parameters via the proposed EO. The obtained results via the proposed EO are compared to those obtained by the Archimedes optimization algorithm (AOA), Heap-based optimizer (HBO), Seagull Optimization Algorithm (SOA), Student Psychology Based Optimization Algorithm (SPBO), Marine predator algorithm (MPA), Manta ray foraging optimization (MRFO), and comprehensive learning dynamic multi-swarm marine predators algorithm. In the case of the SOFC steady-state model, the proposed EO succeeded in achieving the best (minimum) fitness function of 2.6906 × 10−6, 1.5527 × 10−6, 2.6809 × 10−6, and 2.2995 × 10−6 at operating temperature of 1073 K, 1173 K, 1213 K, and 1273 K, respectively. The corresponding standard deviations in the four studied cases obtained via the proposed EO are 9.15673 × 10−7, 1.63132 × 10−6, 3.36081 × 10−6, and 2.2119 × 10−6. Regarding the obtained results of the SOFC dynamic-state model, the proposed EO outperformed the others, achieving the minimum SMSE with a value of 1.0406; the MRFO comes in the second rank with a fitness function of 1.0775, and then the CLDMMPA achieves a SMSE of 1.3204 and comes in the third rank. The proposed EO succeeded in achieving a variance of 0.02264 and a standard deviation of 0.15048 in this studied case. The findings of this study demonstrate the superiority and reliability of the proposed approach in constructing a good-performance model that converges to the real one.

Author Contributions

All authors collaborated and contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

This research received funding from the Institutional Fund Projects under grant no. (IFPHI-173-135-2020) supported by the Ministry of Education and King Abdulaziz University, Deanship of Scientific Research (DSR), Jeddah, Saudi Arabia.

Acknowledgments

This research work was funded by the Institutional Fund Projects under grant no. (IFPHI-173-135-2020). Therefore, authors gratefully acknowledge technical and financial support from the Ministry of Education and King Abdulaziz University, DSR, Jeddah, Saudi Arabia.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Capodaglio, A.G.; Cecconet, D.; Molognoni, D. An integrated mathematical model of microbial fuel cell processes: Bioelectrochemical and microbiologic aspects. Processes 2017, 5, 73. [Google Scholar] [CrossRef] [Green Version]
  2. Bizon, N.; Thounthong, P. Energy efficiency and fuel economy of a fuel cell/renewable energy sources hybrid power system with the load-following control of the fueling regulators. Mathematics 2020, 8, 151. [Google Scholar] [CrossRef] [Green Version]
  3. Wilberforce, T.; El-Hassan, Z.; Khatib, F.; Al Makky, A.; Baroutaji, A.; Carton, J.G.; Olabi, A.G. Developments of electric cars and fuel cell hydrogen electric cars. Int. J. Hydrogen Energy 2017, 42, 25695–25734. [Google Scholar] [CrossRef] [Green Version]
  4. El-Hay, E.A.; El-Hameed, M.A.; El-Fergany, A.A. Performance enhancement of autonomous system comprising proton exchange membrane fuel cells and switched reluctance motor. Energy 2018, 163, 699–711. [Google Scholar] [CrossRef]
  5. Rokni, M. Addressing fuel recycling in solid oxide fuel cell systems fed by alternative fuels. Energy 2017, 137, 1013–1025. [Google Scholar] [CrossRef] [Green Version]
  6. Olabi, A. Renewable energy and energy storage systems. Energy 2017, 136, 1–6. [Google Scholar] [CrossRef]
  7. Ramadhani, F.; Hussain, M.; Mokhlis, H.; Hajimolana, S. Optimization strategies for Solid Oxide Fuel Cell (SOFC) application: A literature survey. Renew. Sustain. Energy Rev. 2017, 76, 460–484. [Google Scholar] [CrossRef]
  8. Chowdhury, S.; Chowdhury, S.P.; Crossley, P. Microgrids and Active Distribution Networks; The Institution of Engineering and Technology: London, UK, 2009. [Google Scholar]
  9. El-Hay, E.; El-Hameed, M.; El-Fergany, A. Steady-state and dynamic models of solid oxide fuel cells based on Satin Bowerbird Optimizer. Int. J. Hydrogen Energy 2018, 43, 14751–14761. [Google Scholar] [CrossRef]
  10. El-Hay, E.; El-Hameed, M.; El-Fergany, A. Optimized parameters of SOFC for steady state and transient simulations using interior search algorithm. Energy 2019, 166, 451–461. [Google Scholar] [CrossRef]
  11. Van Biert, L.; Godjevac, M.; Visser, K.; Aravind, P. Dynamic modelling of a direct internal reforming solid oxide fuel cell stack based on single cell experiments. Appl. Energy 2019, 250, 976–990. [Google Scholar] [CrossRef]
  12. Gong, W.; Yan, X.; Hu, C.; Wang, L.; Gao, L. Fast and accurate parameter extraction for different types of fuel cells with decomposition and nature-inspired optimization method. Energy Convers. Manag. 2018, 174, 913–921. [Google Scholar] [CrossRef]
  13. Ettihir, K.; Boulon, L.; Becherif, M.; Agbossou, K.; Ramadan, H. Online identification of semi-empirical model parameters for PEMFCs. Int. J. Hydrogen Energy 2014, 39, 21165–21176. [Google Scholar] [CrossRef]
  14. Ang, S.M.C.; Brett, D.J.; Fraga, E.S. A multi-objective optimisation model for a general polymer electrolyte membrane fuel cell system. J. Power Sources 2010, 195, 2754–2763. [Google Scholar] [CrossRef] [Green Version]
  15. Tahmasbi, A.A.; Hoseini, A.; Roshandel, R. A new approach to multi-objective optimisation method in PEM fuel cell. Int. J. Sustain. Energy 2013, 34, 283–297. [Google Scholar] [CrossRef]
  16. Petrescu, S.; Petre, C.; Costea, M.; Malancioiu, O.; Boriaru, N.; Dobrovicescu, A.; Feidt, M.; Harman, C.; Stanciu, C. A methodology of computation, design and optimization of solar Stirling power plant using hydrogen/oxygen fuel cells. Energy 2010, 35, 729–739. [Google Scholar] [CrossRef]
  17. Virkar, A.; Williams, M.C.; Singhal, S. Concepts for ultra-high power density solid oxide fuel cells. ECS Trans. 2007, 5, 401–421. [Google Scholar] [CrossRef]
  18. Zhu, L.; Zhang, L.; Virkar, A.V. A parametric model for solid oxide fuel cells based on measurements made on cell materials and components. J. Power Sources 2015, 291, 138–155. [Google Scholar] [CrossRef] [Green Version]
  19. Shi, H.; Li, J.; Zafetti, N. New optimized technique for unknown parameters selection of SOFC using Converged Grass Fibrous Root Optimization Algorithm. Energy Rep. 2020, 6, 1428–1437. [Google Scholar] [CrossRef]
  20. Yousri, D.; Hasanien, H.M.; Fathy, A. Parameters identification of solid oxide fuel cell for static and dynamic simulation using comprehensive learning dynamic multi-swarm marine predators algorithm. Energy Convers. Manag. 2021, 228, 113692. [Google Scholar] [CrossRef]
  21. Nassef, A.M.; Fathy, A.; Sayed, E.T.; Abdelkareem, M.A.; Rezk, H.; Tanveer, W.H.; Olabi, A. Maximizing SOFC performance through optimal parameters identification by modern optimization algorithms. Renew. Energy 2019, 138, 458–464. [Google Scholar] [CrossRef]
  22. Fathy, A.; Rezk, H.; Ramadan, H.S.M. Recent moth-flame optimizer for enhanced solid oxide fuel cell output power via optimal parameters extraction process. Energy 2020, 207, 118326. [Google Scholar] [CrossRef]
  23. Wang, X.; Huang, B.; Chen, T. Data-driven predictive control for solid oxide fuel cells. J. Process. Control 2007, 17, 103–114. [Google Scholar] [CrossRef]
  24. Larminie, J.; Dicks, A.; McDonald, M.S. Fuel Cell Systems Explained; John Wiley & Sons: Chichester, UK, 2003; Volume 2. [Google Scholar]
  25. Faramarzi, A.; Heidarinejad, M.; Stephens, B.; Mirjalili, S. Equilibrium optimizer: A novel optimization algorithm. Knowl. Based Syst. 2020, 191, 105190. [Google Scholar] [CrossRef]
  26. Wei, Z.; Zou, C.; Leng, F.; Soong, B.H.; Tseng, K.-J. Online model identification and state-of-charge estimate for lithium-ion battery with a recursive total least squares-based observer. IEEE Trans. Ind. Electron. 2018, 65, 1336–1346. [Google Scholar] [CrossRef]
  27. Wei, Z.; Zhao, J.; Xiong, R.; Dong, G.; Pou, J.; Tseng, K.J. Online estimation of power capacity with noise effect attenuation for lithium-ion battery. IEEE Trans. Ind. Electron. 2019, 66, 5724–5735. [Google Scholar] [CrossRef]
  28. Wei, Z.; Zhao, D.; He, H.; Cao, W.; Dong, G. A noise-tolerant model parameterization method for lithium-ion battery management system. Appl. Energy 2020, 268, 114932. [Google Scholar] [CrossRef]
  29. Wei, Z.; Meng, S.; Xiong, B.; Ji, D.; Tseng, K.J. Enhanced online model identification and state of charge estimation for lithium-ion battery with a FBCRLS based observer. Appl. Energy 2016, 181, 332–341. [Google Scholar] [CrossRef]
  30. Pierre, J. Siemens energy. In Proceedings of the 11th Annual SECA Workshop, Pittsburgh, PA, USA, 27–29 July 2010. [Google Scholar]
  31. Xu, D.; Jiang, B.; Liu, F. Improved data driven model free adaptive constrained control for a solid oxide fuel cell. IET Control. Theory Appl. 2016, 10, 1412–1419. [Google Scholar] [CrossRef]
Mathematics 09 01066 g001 550
Figure 1. Optimization process of EO.
Figure 1. Optimization process of EO.
Mathematics 09 01066 g001
Mathematics 09 01066 g002 550
Figure 2. The proposed steps incorporating EO.
Figure 2. The proposed steps incorporating EO.
Mathematics 09 01066 g002
Mathematics 09 01066 g003 550
Figure 3. The measured and calculated polarization curves of SOFC operated at 1073 K obtained via EO at (a) current density-voltage, (b) current density-power.
Figure 3. The measured and calculated polarization curves of SOFC operated at 1073 K obtained via EO at (a) current density-voltage, (b) current density-power.
Mathematics 09 01066 g003
Mathematics 09 01066 g004 550
Figure 4. (a) Current-voltage curve, (b) Current-power curve of SOFC operated at 1073 K obtained via other approaches.
Figure 4. (a) Current-voltage curve, (b) Current-power curve of SOFC operated at 1073 K obtained via other approaches.
Mathematics 09 01066 g004
Mathematics 09 01066 g005 550
Figure 5. The variation of fitness function during iterative process for all employed optimizers applied for steady-state SOFC model.
Figure 5. The variation of fitness function during iterative process for all employed optimizers applied for steady-state SOFC model.
Mathematics 09 01066 g005
Mathematics 09 01066 g006 550
Figure 6. The measured and calculated (a) current density-voltage, (b) current density-power of SOFC operated at 1173 K, 1213 K, and 1273 K obtained via the proposed EO.
Figure 6. The measured and calculated (a) current density-voltage, (b) current density-power of SOFC operated at 1173 K, 1213 K, and 1273 K obtained via the proposed EO.
Mathematics 09 01066 g006
Mathematics 09 01066 g007 550
Figure 7. The measured and calculated polarization curves of SOFC dynamic-state model operated at 1273 K obtained via EO (a) current-voltage, (b) current-power.
Figure 7. The measured and calculated polarization curves of SOFC dynamic-state model operated at 1273 K obtained via EO (a) current-voltage, (b) current-power.
Mathematics 09 01066 g007
Mathematics 09 01066 g008 550
Figure 8. The measured and calculated polarization curves of SOFC dynamic-state model operated at 1273 K obtained via MPA, HBO, SOA, and MRFO.
Figure 8. The measured and calculated polarization curves of SOFC dynamic-state model operated at 1273 K obtained via MPA, HBO, SOA, and MRFO.
Mathematics 09 01066 g008
Mathematics 09 01066 g009 550
Figure 9. The variation of fitness function during iterative process for all employed optimizers applied for dynamic-state SOFC model.
Figure 9. The variation of fitness function during iterative process for all employed optimizers applied for dynamic-state SOFC model.
Mathematics 09 01066 g009
Mathematics 09 01066 g010 550
Figure 10. First load disturbance applied on SOFC stack dynamic model.
Figure 10. First load disturbance applied on SOFC stack dynamic model.
Mathematics 09 01066 g010
Mathematics 09 01066 g011 550
Figure 11. Second load disturbance applied on SOFC stack dynamic model.
Figure 11. Second load disturbance applied on SOFC stack dynamic model.
Mathematics 09 01066 g011
Table
Table 1. The optimal parameters of SOFC steady-state based model operated at 1073 K.
Table 1. The optimal parameters of SOFC steady-state based model operated at 1073 K.
CLDMMPA [20]MPAHBOSOAMRFOThe Proposed EO
E0 (V)0.907540.91270.912140.911010.918270.91056
a (V)0.0107410.0110580.0109290.0205020.01160.010724
Jo (A/cm2)0.0986270.0599940.0633210.0481270.0359180.074522
r (kΩ∙cm2)1.01.01.01.00.999191.0
b (V)0.0441040.0427840.0435020.00.0368860.044165
Jmax (A/cm2)1.01.01.00.642970.91831.0
Elapsed time (sec.)NA468.526403.607278.254560.753272.198102
SMSE2.6906 × 10−62.692 × 10−62.7003 × 10−64.123 × 10−62.7213 × 10−62.6906 × 10−6
Table
Table 2. The optimal parameters of SOFC steady-state based model operated at 1173 K.
Table 2. The optimal parameters of SOFC steady-state based model operated at 1173 K.
CLDMMPA [20]MPAHBOSOAMRFOThe Proposed EO
E0 (V)0.891030.89129 0.890830.879560.890870.89108
a (V)3.671 × 10−136.1748 × 10−133.23334 × 10−130.00775679.8137 × 10−63.14568 × 10−8
Jo (A/cm2)0.0951270.0183580.0411940.0873630.0686840.09998
r (kΩ∙cm2)0.404730.415930.399520.211390.40270.40610
b (V)0.188410.177410.192120.284970.188850.18741
Jmax (A/cm2)1.00.988621.00.94650.99680.99999
Elapsed time (sec.)NA495.653992396.181972303.170716612.841366303.185040
SMSE1.5529 × 10−61.5594 × 10−61.5557 × 10−63.1657 × 10−61.5563 × 10−61.5527 × 10−6
Table
Table 3. The optimal parameters of SOFC steady-state based model operated at 1213 K.
Table 3. The optimal parameters of SOFC steady-state based model operated at 1213 K.
CLDMMPA [20]MPAHBOSOAMRFOThe Proposed EO
E0 (V)0.861690.861890.861340.856220.86176 0.86164
a (V)5.0588 × 10−134.2805 × 10−283.4456 × 10−303.1223 × 10−293.3451 × 10−57.12575 × 10−8
Jo (A/cm2)0.088710.0545680.0119750.0478850.0634680.07768
r (kΩ∙cm2)0.159820.166330.148580.0001240.166290.15873
b (V)0.285290.280320.293510.40120.279180.28603
Jmax (A/cm2)1.00.999461.01.00.997320.99999
Elapsed time (sec.)NA464.968854367.890819236.154633571.548052273.530
SMSE2.6811 × 10−62.6846 × 10−62.6887 × 10−64.3804 × 10−62.6896 × 10−62.6809 × 10−6
Table
Table 4. The optimal parameters of SOFC steady-state based model operated at 1273 K.
Table 4. The optimal parameters of SOFC steady-state based model operated at 1273 K.
CLDMMPA [20]MPAHBOSOAMRFOThe Proposed EO
E0 (V)0.84780.847820.84802 0.840140.84802 0.8478
a (V)2.887 × 10−143.128 × 10−215.1251 × 10−62.238 × 10−203.2243 × 10−56.7797 × 10−12
Jo (A/cm2)0.0618160.0145230.0862830.0172850.0244370.0160
r (kΩ∙cm2)0.215640.216340.223230.15630.220210.2169
b (V)0.205750.205240.200090.360460.202260.2047
Jmax (A/cm2)1.01.00.999841.000.998881.0
Elapsed time (sec.)NA465.757529352.634515199.242260555.636660317.980527
SMSE2.2997 × 10−62.2996 × 10−62.3031 × 10−65.4888 × 10−62.3058 × 10−62.2995 × 10−6
Table
Table 5. Statistical parameters (best, worst, mean, median, variance, and standard deviation) of all optimizers used for steady-state SOFC model.
Table 5. Statistical parameters (best, worst, mean, median, variance, and standard deviation) of all optimizers used for steady-state SOFC model.
At T = 1073 K
Cldmmpa [20]MPAHBOSOAMRFOThe Proposed EO
Best2.6906 × 10−62.69203 × 10−62.70035 × 10−64.12297 × 10−62.7213 × 10−62.6906 × 10−6
Worst2.696 × 10−61.23064 × 10−54.49932 × 10−60.000344.8504 × 10−66.8151 × 10−6
Mean2.6926 × 10−64.75170 × 10−63.11043 × 10−69.38355 × 10−53.1121 × 10−63.3645 × 10−6
Median2.6917 × 10−63.61749 × 10−62.96160 × 10−69.33996 × 10−52.9395 × 10−62.6906 × 10−6
Variance 4.1347 × 10−186.38059 × 10−121.75759 × 10−139.20085 × 10−92.2231 × 10−138.38802 × 10−13
Std. deviation2.0334 × 10−92.52598 × 10−64.19236 × 10−79.59210 × 10−54.7150 × 10−79.15673 × 10−7
At T = 1173 K
CLDMMPA [20]MPAHBOSOAMRFOThe proposed EO
Best1.5529 × 10−61.5594 × 10−61.5556 × 10−63.1656 × 10−61.5563 × 10−61.55279 × 10−6
Worst1.5752 × 10−65.9028 × 10−62.1420 × 10−60.011022.6603 × 10−66.00407 × 10−6
Mean1.5593 × 10−63.0125 × 10−61.6979 × 10−60.000261.7105 × 10−62.31364 × 10−6
Median1.5572 × 10−62.4861 × 10−61.6532 × 10−61.26093 × 10−51.6493 × 10−61.58630 × 10−6
Variance 3.8987 × 10−171.9116 × 10−121.8382 × 10−142.42265 × 10−63.7587 × 10−142.66121 × 10−12
Std. deviation6.244× 10−91.3826 × 10−61.3558 × 10−70.001551.9387 × 10−71.63132 × 10−6
At T = 1213 K
CLDMMPA [20]MPAHBOSOAMRFOThe proposed EO
Best2.6811 × 10−62.6846 × 10−62.6887 × 10−64.3804 × 10−62.68961 × 10−62.68099 × 10−6
Worst2.7269 × 10−61.2325 × 10−53.4295 × 10−60.000323.22001 × 10−61.30287 × 10−5
Mean2.6921 × 10−64.6204 × 10−62.8357 × 10−62.37574 × 10−52.85206 × 10−64.02438 × 10−6
Median2.6867 × 10−63.3119 × 10−62.7703 × 10−61.31251 × 10−52.81486 × 10−62.78147 × 10−6
Variance 1.9891 × 10−169.4614 × 10−122.7256 × 10−143.77968 × 10−91.72574 × 10−141.12951 × 10−11
Std. deviation1.4103 × 10−83.0759 × 10−61.6509 × 10−76.14791 × 10−51.3136 × 10−73.36081 × 10−6
At T = 1273 K
CLDMMPA [20]MPAHBOSOAMRFOThe proposed EO
Best2.2997 × 10−62.2995 × 10−62.30310 × 10−65.4888 × 10−62.30579 × 10−62.2995 × 10−6
Worst2.3392 × 10−67.4728 × 10−63.0608 × 10−60.0002472.84248 × 10−67.6115 × 10−6
Mean2.3156 × 10−63.6859 × 10−62.4184 × 10−61.758508 × 10−52.40213 × 10−63.4890 × 10−6
Median2.3141 × 10−63.1707 × 10−62.3647 × 10−67.71186 × 10−62.37664 × 10−62.3085 × 10−6
Variance 1.3448 × 10−162.1494 × 10−121.9816 × 10−142.25858 × 10−98.64055 × 10−154.8928 × 10−12
Std. deviation1.1597 × 10−81.466 × 10−61.4077 × 10−74.7524 × 10−59.29545 × 10−82.2119 × 10−6
Table
Table 6. SOFC stack specifications [30,31].
Table 6. SOFC stack specifications [30,31].
ParameterValue
Prated (W)100 kW
nc384
E0 (V)1.18
T (K)1273
KH2 (kmol/s/atm)8.43 × 10−4
KO2 (kmol/s/atm)2.52 × 10−3
KH2O (kmol/s/atm)2.81 × 10−4
rH-O1.145
τH2 (s)26.1
τO2 (s)2.91
τH2O (s)78.3
Te (s)0.8
Table
Table 7. The optimal parameters of SOFC dynamic-state based model operated at 1273 K.
Table 7. The optimal parameters of SOFC dynamic-state based model operated at 1273 K.
CLDMMPA [20]MPAHBOSOAMRFOThe Proposed EO
E0 (V)1.14051.1998471.1949390.894612711.1133371.144398
a (V)0.0374490.048823350.047561760.0000.027290010.02982692
Jo (A/cm2)0.0954420.099870290.09395040.0906555850.031437290.0202932
r (kΩ∙cm2)0.00018294.3684 × 10−54.873756 × 10−50.0000.00026720470.0002547476
b (V)0.103860.15514250.15150450.326728950.082006610.08350225
Jmax (A/cm2)0.8367865.6602859.63171000825.4696825.7578
Elapsed time (sec.)NA433.946266331.304354222.830499593.725886327.35802
SMSE1.32043.08873.348634.16921.07751.0406
Table
Table 8. Statistical parameters (best, worst, mean, median, variance, and standard deviation) of all optimizers used for dynamic-state SOFC model.
Table 8. Statistical parameters (best, worst, mean, median, variance, and standard deviation) of all optimizers used for dynamic-state SOFC model.
CLDMMPA [20]MPAHBOSOAMRFOThe Proposed EO
Best1.32043.088673.348634.16921.07751.04061
Worst3.983525.113510.4735401.87253.41642.00140
Mean3.24628.549645.62091025.80991.71151.1352
Median3.52416.075934.883534.34734.51411.0816
Variance 0.5102736.35472.75283,915,661.9570.233600.02264
Std. deviation0.714346.029481.65921978.8030.152840.15048
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Alhumade, H.; Fathy, A.; Al-Zahrani, A.; Rawa, M.J.; Rezk, H. Optimal Parameter Estimation Methodology of Solid Oxide Fuel Cell Using Modern Optimization. Mathematics 2021, 9, 1066. https://doi.org/10.3390/math9091066

AMA Style

Alhumade H, Fathy A, Al-Zahrani A, Rawa MJ, Rezk H. Optimal Parameter Estimation Methodology of Solid Oxide Fuel Cell Using Modern Optimization. Mathematics. 2021; 9(9):1066. https://doi.org/10.3390/math9091066

Chicago/Turabian Style

Alhumade, Hesham, Ahmed Fathy, Abdulrahim Al-Zahrani, Muhyaddin Jamal Rawa, and Hegazy Rezk. 2021. "Optimal Parameter Estimation Methodology of Solid Oxide Fuel Cell Using Modern Optimization" Mathematics 9, no. 9: 1066. https://doi.org/10.3390/math9091066

APA Style

Alhumade, H., Fathy, A., Al-Zahrani, A., Rawa, M. J., & Rezk, H. (2021). Optimal Parameter Estimation Methodology of Solid Oxide Fuel Cell Using Modern Optimization. Mathematics, 9(9), 1066. https://doi.org/10.3390/math9091066

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