ORIGINAL RESEARCH article

Front. Energy Res., 24 November 2025

Sec. Fuel Cells, Electrolyzers and Membrane Reactors

Volume 13 - 2025 | https://doi.org/10.3389/fenrg.2025.1659232

Energy-efficient parameter estimation of solid oxide fuel cells under varying pressure conditions using the black widow optimization algorithm

  • 1. Department of Chemical Engineering, Thapar Institute of Engineering and Technology, Patiala, India

  • 2. Chitkara University Institute of Engineering and Technology, Chitkara University, Rajpura, Punjab, India

  • 3. Department of Mechanical Engineering, Indian Institute of Technology (ISM), Dhanbad, Jharkhand, India

  • 4. Department of Multidisciplinary Engineering, The NorthCap University, Gurugram, Haryana, India

  • 5. Department of Mechanical Engineering, Faculty of Engineering, Universidad Tecnológica Metropolitana, Santiago, Chile

  • 6. State Key Laboratory of Intelligent Manufacturing Equipment and Technology, School of Mechanical Science and Engineering, Huazhong University of Science and Technology, Wuhan, China

Abstract

Solid oxide fuel cells (SOFCs) are highly efficient and fuel-flexible energy conversion devices, but accurately estimating their governing parameters remains a challenge due to the nonlinear behavior of electrochemical processes. This study presents the first application of the black widow optimization (BWO) algorithm for estimating six critical SOFC parameters—open-circuit potential (E0), Tafel slope (A), exchange current density (I0), concentration loss coefficient (B), limiting current density (Il), and ohmic resistance (Rohm)—under varying pressure conditions (1–5 atm). The objective was to minimize the mean squared error (MSE) between experimental and predicted polarization curves while ensuring computational efficiency. The proposed BWO framework achieved superior accuracy, with an MSE of 0.52 at 5 atm and convergence within 3.74 s, significantly outperforming benchmark metaheuristic algorithms such as particle swarm optimization (PSO), gray wolf optimization (GWO), and the whale optimization algorithm (WOA). Robustness was confirmed through cross-validation, where polarization curves predicted at unseen conditions deviated by less than 5% from experimental results. This demonstrates that the estimated parameters effectively capture intrinsic SOFC electrochemical behavior rather than overfitting specific datasets. Beyond numerical accuracy, the optimized parameters enhanced the predictive stability of voltage–current (V–I) and power–current (P–I) characteristics across all studied pressures, directly supporting improved operational reliability and long-term stack durability. The combination of higher precision, faster convergence, and strong generalizability positions BWO as a promising tool for real-time SOFC optimization. The findings establish a robust framework for parameter identification that not only reduces uncertainty in SOFC modeling but also contributes to practical advances in performance optimization and system longevity. Future extensions of this research will include real-time implementation under dynamic operating environments and integration with hybrid renewable energy systems to improve scalability, efficiency, and sustainability.

1 Introduction

The depletion of fossil fuels, environmental pollution caused by their use, and global dependence on these finite resources make the transition to clean, renewable alternatives one of the major challenges for future energy security. Consequently, much research has explored novel and eco-friendly energy technologies. Among this, fuel cells have emerged as highly significant in complementing renewable sources such as solar, wind, geothermal, and biomass (). Fuel cells and batteries share some similarities, but differ in one crucial aspect: while batteries store and subsequently release chemical energy with a lifetime limited by stored reactants, fuel cells can operate continuously as long as an external supply of fuel and oxidant is maintained (). Furthermore, various types of fuel cells can be powered by ethanol, hydrogen, methanol, natural gas, or even biogas. Solid oxide fuel cells (SOFCs) typically operate at high temperatures, whereas proton exchange membrane fuel cells (PEMFCs) function at a lower range of 50–100 °C (; ). Through electrochemical processes, these devices convert the chemical energy of fuels into electricity. However, the performance of SOFCs is strongly influenced by activation, concentration, and Ohmic polarization losses (). In turn, SOFC material properties such as electrode absorbency, tortuosity, and pore radius critically affect these polarizations (; ). The accurate experimental determination of these parameters remains challenging, thereby necessitating mathematical modeling to represent the system’s complex dynamics. Contemporary SOFC models often incorporate values drawn from widely varying sources (; ). As emphasized by , , and , a precise mathematical model is indispensable for predicting SOFC behavior and determining unknown parameters under diverse operating conditions. developed a 3-D dynamic model of a planar anode-supported SOFC to investigate transient responses under load variations, showing that excessive current density increases can lead to fuel starvation in the porous anode layer, accelerating degradation. This study further demonstrated that increasing anode support thickness mitigates starvation and improves efficiency, providing guidelines for SOFC design and control. Nevertheless, parameter estimation remains particularly complex because SOFC models are inherently nonlinear. Conventional numerical or linear programming methods have frequently failed to provide robust optimization, largely due to local optima entrapment and sensitivity to gradient initialization (; ; ). To address these challenges, metaheuristic optimization strategies have gained substantial popularity in recent years. For example, introduced a converged grass fibrous root optimizer, while applied marine-predator-based optimization for SOFC modeling. proposed an interior random search algorithm, and applied radial movement optimization. Although these approaches improved accuracy, they often suffered from slower convergence or limited robustness. , through their work on the optimal coupling of fuel cells and supercapacitors, further emphasized the importance of parameter optimization in enhancing both system cost-effectiveness and durability. In parallel, research on protonic ceramic fuel cells (PCFCs) has provided additional theoretical insights into electrochemical behavior under varying operating conditions. developed a theoretical framework demonstrating that hydrogen humidification exerts a strong, non-monotonic influence on current density, whereas air humidification plays only a minor role. They further suggested that optimizing hydrogen inlet humidity can reduce polarization losses, thereby improving device efficiency. Extending this research, investigated the combined effects of cell voltage and external operating conditions, showing that proton conductivity may either increase or decrease, depending on humidity, while hole conductivity consistently rises with voltage, reducing Faradaic efficiency. Importantly, their findings revealed that fundamental parameters such as conductivities, transport numbers, and EMF differ markedly between open-circuit and operating conditions. These contributions underscore the necessity of advanced modeling strategies and highlight the pivotal role of environmental and voltage-dependent factors in fuel cell performance optimization.

1.1 Novelty and contributions

While recent studies such as have successfully applied optimization methods for SOFC parameter estimation under different temperature conditions, the present research extends this concept to pressure-dependent operating environments.

The main novelty and contributions are

  • first application of the black widow optimization (BWO) algorithm for SOFC parameter estimation under varying pressure conditions;

  • comparison against leading metaheuristics (PSO, GWO, and WOA) demonstrating superior accuracy, convergence speed, and robustness;

  • improved predictive reliability of polarization curves across multiple operating pressures with < 5% deviation from experimental results;

  • contribution toward long-term SOFC stability and longevity by reducing parameter uncertainty and mismatch.

The present study aims to establish a robust and generalizable framework for accurate parameter estimation in SOFCs by employing the BWO algorithm. The specific objectives are to

  • develop and implement the BWO algorithm to extract six critical SOFC parameters (E0, A, I0, B, IL, and Rohm) that govern polarization behavior;

  • validate the proposed BWO framework across multiple operating pressures (1–5 atm) using experimental polarization data from a commercial 5-kW tubular SOFC stack;

  • perform cross-validation on unseen experimental conditions (e.g., testing the model at 5 atm while training on 1–2–4 atm) in order to verify the robustness and generalizability of the estimated parameters;

  • conduct a comparative analysis with established metaheuristic algorithms such as particle swarm optimization (PSO), gray wolf optimization (GWO), and the whale optimization algorithm (WOA), evaluating convergence speed, accuracy, and stability.

  • demonstrate the practical impact of precise parameter estimation by linking reduced model–experiment mismatch to improved predictive stability, better operational reliability, and the potential extension of SOFC longevity.

2 Methodology

2.1 SOFC concept

The fuel cell is a device for electrochemical energy conversion, continuously producing heat and electricity from chemical energy. In SOFCs, hydrogen acts as the fuel while oxygen (from air) is used as the oxidant. Unlike combustion-based engines, SOFCs directly convert chemical energy into electrical energy without intermediate combustion, thereby avoiding harmful flue gas emissions (). The anodic hydrogen oxidation, cathodic oxygen reduction, and the overall electrochemical reaction in an SOFC are represented by Equations 13.

Figure 1 displays the SOFC schematic. Reducing oxygen molecules at the cathode produces negative oxygen ions. The ionic conduction electrolyte conducts negative ions, while the thin electrolyte sheet blocks electrons ().

FIGURE 1

).

2.2 SOFC voltage description

The estimation of the output potential (voltage) of a single SOFC can be achieved by utilizing the thermodynamic potential, ENernst, obtained chemically under no load conditions, along with the potential losses occurring during the reaction. This estimation is denoted as Vcell. The cell voltage accounting for thermodynamic and loss components is defined in Equation 4. The losses encompassed in this context consist of the Ohmic potential drop (Vohm), the concentration potential drop (Vconc), and the activation potential drop (Vact) ().

The voltage, also known as the thermodynamic potential, is typically described by the Nernst equation for hydrogen and oxygen. This equation is applicable when operating at a temperature denoted as T (Kelvin) and at specific partial pressures of hydrogen (), oxygen (), and water () ().—where represents potential (standard), R is equal to 8.314 J/mol/K, and F is 96486 C/mol (Faraday constant).

The activation potential drop (Vact) is denoted by Butler–Volmer expression (), represented by Equation 6:—where is the load current density (mA/cm2), is the current density (exchange) (mA/cm2), and A is the Tafel line.

Potential drop (concentration) (Vconc) is shown as

The coefficient B is an unidentified parameter that is dependent on the specific operating conditions, while IL represents the restrictive current density.

The potential drop (Ohmic) (Vohm) is shown as—where Rohm is the ionic resistance expressed correctly as kΩ·cm2.

Hence, the cell voltage is shown as

For SOFC stack system Ncell, the potential drop (overall) is given as (). The total stack voltage for Ncell series-connected cells is defined in Equation 10.

With respect to Equations 9, 10 above, it is possible to calculate Vcell (or Vstack) and obtain the I-V and P-V characteristics of the SOFC model by assuming that , A, , B, , and have predetermined outcomes. These values in different working scenarios are not adequately comprehended in practical applications. Hence, the objective is to precisely delineate outcomes.

2.3 Objective function for optimization

SOFC stack parameters are estimated using an optimization problem. The goal is to optimize the gap between an SOFC system’s electrochemical model and observations. To forecast SOFC stack voltage, we optimize the six parameters , A, , B, , and in Equations 58. Minimum average square error represented in Equation 11 is required between observed and modeled stack voltages. BWO is a global meta-heuristic optimization method for solving the objective function. This is presented as—where N represents measured data points and is the vector of six unidentified variables. The constraint limits for the optimization problem are defined in Equation 12. Thus, the bounded range of these constraints is as follows:

—where is the minimum limit and denotes the maximum limit.

3 Black widow optimization algorithm

A 2020 research group presented a new metaheuristic optimization approach based on the spider’s mating behavior due to its simplicity and adaptability (), which has solved technological and scientific problems. This mating was the inspiration for the BWO algorithm. This stage causes early convergence because poorly suited species are excluded from the circle. BWO performs well in exploitation and exploration, avoids local optima, and converges quickly. BWO can also balance exploitation and discovery. BWO can study over wide ranges to obtain the optimum global answer, making it a viable choice for optimization with multiple local optima. Following are the steps involved in the BWO algorithm.

3.1 Initial population

Problem variables must be structured to answer the existing situation before an optimization situation can be resolved. GA and PSO call this arrangement a “chromosome” and “particle position,” respectively; BWO calls it a “widow.” The BWO algorithm represents each issue solution as a black widow spider. Each black widow spider displays the problem variables. The fitness solution (widow) is represented by Equation 13 that consists of Nvar multidimensional optimization problem’s solution as an array of values.

3.2 Procreate

The pairings mate separately to produce the next generation. In nature, every couple mates, distinct from others. In reality, every pair gives 1,000 eggs, but some spider progeny survive and are stronger. For replication, this approach requires a matrix alpha and a widow array with random values. Offspring are produced using Equation 14, which has x1 and x2 as parents and y1 and y2 as offspring:

This strategy should prevent duplicate random numbers by repeating it times.

3.3 Cannibalism

Three distinct types of cannibalism are considered within the BWO framework. The first occurs during mating, where the female black widow cannibalizes her partner, effectively distinguishing individuals based on relative fitness. The second involves sibling cannibalism, in which stronger spiderlings eliminate their weaker counterparts; survivability in this process is determined using a defined cannibalism rating (CR). The third type, maternal cannibalism, arises when offspring consume the mother. In all cases, spiderling strength is evaluated based on the calculated fitness values.

3.4 Mutation

In the mutation phase, a subset of individuals (“Mute pop”) is randomly selected from the population (). The selected variables are subjected to random swapping of two array entries. The size of the Mute pop is determined according to the defined mutation rate.

3.5 Convergence

These three stop conditions are similar to evolutionary optimization: (1) initiate iterations; (2) noting that the ideal fitness value has not altered throughout iterations; and (3) precision necessary.

3.6 Parameter conditioning

The proposed algorithm’s settings are key to better results. These qualities include cannibalism, mutation, and procreation. This study’s factor values are in Table 1. Controlling the right number of parameters balances the exploitation and discovery stages. The BWO algorithm has three important regulatory parameters: PP, CR, and MR. The procreating rate (PP) determines the number of procreators. This parameter promotes differentiation and improves space search accuracy by regulating offspring birth. The cannibalism operative’s CR excludes the incorrect people. The correct search variables can ensure good exploitation stage performance between local and global minima. A proportion of mutation participants is MR. This parameter ensures exploitation and discovery balance. This option controls how search agents migrate from global to local and points them to the best answer. Figure 2 shows the proposed algorithm’s flowchart ().

TABLE 1

AlgorithmsOperational criterion
Particle swarm optimization (PSO)Swarm = 200; w = 0.8
Ideal factor, c1 = 1.0& c2 = 2.0
Overall repetitions = 2000
Grey wolf optimization (GWO)Size = 200
Constant (c) = 1.0
Overall repetitions = 2000
Whale optimization algorithm (WOA)Size = 200; a = −1 to −2 (default)
Overall repetitions = 2000
Black widow optimization (BWO)(PP) = 0.50; (CR) = 0.40; (MR) = 0.40
Overall repetitions = 2000

Operational parameters of studied algorithms.

FIGURE 2

).

The suggested algorithmic approach has been implemented using a mathematical model of an SOFC. Specifically, a 5-kW dynamic tubular SOFC system, as reported by , was considered for this study. The operational data of the investigated stack are summarized in Table 2, while Table 3 presents the defined investigation bounds for the six unidentified parameters to be estimated.

TABLE 2

Operational data and conditionsValues
No. of cells, Ncells96
Power (kW)5
Temperature1,027.15 K
Stack’s operating pressures1 atm; 2 atm; 4 atm; 5 atm
Partial pressures of hydrogen () and oxygen ()0.91 and 0.21 respectively

Data sheet test stack ().

TABLE 3

VariablesMinimum boundMaximum bound
(V)01.2
(V)01
(mA/cm2)0100
(V)01
(mA/cm2)010000
(kΩ.cm2)01

Working condition of the system.

4 Results and discussion

4.1 Standard benchmark functions

The proposed BWO algorithm was benchmarked against eight well-established test functions to validate global optimization capacity. Table 4 presents the standard benchmark functions employed for performance testing, detailing each function’s mathematical formulation, dimensionality (Dim), search range (R), and global minimum value fmin(x). In all cases, BWO consistently achieved zero or near-zero best cost values, while competing methods (PSO, GWO, and WOA) frequently converged to local optima with higher residual error (Table 5).

TABLE 4

No.Function nameFormulationDimRfmin(x)
f1Sphere30[-5.12, 5.12]0
f2Schwefel 2.2130[-100, 100]0
f3Schwefel 2.2230[-100, 100]0
f4Quartic30[-1.28, 1.28]0
f5Rastrigin30[-5.12, 5.12]0
f6Ackley30[-35,35]0
f7Griewank30[-600, 600]0
f8Camel three hump02[-5,5]0

Performance testing benchmark functions.

TABLE 5

FunctionMeasurePSOGWOWOABWO
f1Best Cost8.76E-152.56E-301.73E-090
Median1.65E-142.95E-301.97E-090
f2Best Cost2.38E-051.34E-062.66E-073.6E-200
Median0.0000531.33E-061.08E-080
f3Best Cost2.38E-054.62E-071.32E-060
Median0.0000533.6E-076.24E-090
f4Best Cost0.003240.001710.002650
Median0.002390.000678.98E-190
f5Best Cost2.75E-111.17E-131.73E-120
Median5.14E-116.95E-142.88E-140
f6Best Cost201.03E-138.62E-148.88E-16
Median0.001451.21E-142.61E-292.04E-31
f7Best Cost1.99E-110.004663.03E-160
Median2.64E-110.009693.15E-160
f8Best Cost1.5E-1037.2E-1929.9E-2000
Median4.3E-103000

Statistical outcomes of studied functions.

These findings confirm that BWO maintains stronger balance between exploration and exploitation. The Wilcoxon rank test (Table 6) further demonstrates the statistical superiority of BWO, with p values < 0.05 across all comparisons. This ensures that BWO is not only faster but also significantly more reliable than competing methods—an essential feature for highly nonlinear SOFC models.

TABLE 6

Related algorithmsBWO vs. PSOBWO vs. GWOBWO vs. WOA
p-values3.08E-334.68E-222.55E-26

Wilcoxon rank test (p ⩽ 0.05).

4.2 Model performance at varying pressures

SOFC stack performance was studied at four pressures (1, 2, 4, and 5 atm) while maintaining a constant temperature of 1,027.15 K. A dataset of 100 experimental points was compared with model predictions. Figure 3 presents the optimizer-specific MSE minima across pressures. BWO achieved lowest MSE in all cases (0.30–0.52), with deviations within 5% of experimental values. In contrast, PSO and GWO exhibited higher variability and slower convergence.

FIGURE 3

The resultant convergence contours for the analyzed methods at different pressures are presented in Figure 4. It can be observed that the proposed strategy successfully avoids local optima and consistently outperforms the comparative algorithms under the conditions examined. The approach demonstrates faster convergence throughout the evolutionary process. Comparative analysis further indicates that by maintaining an appropriate balance between exploration and exploitation, the BWO algorithm is capable of significantly accelerating convergence rates.

FIGURE 4

The polarization plots (Figures 5, 6) show that as the pressure increases, the electrode current density grows while the cell voltage decreases. This trend is consistent with physical expectations since higher pressure improves reactant availability (enhancing current density) but simultaneously increases overpotential losses, lowering output voltage.

FIGURE 5

FIGURE 6

To further interpret the physical meaning of the optimized parameters, we emphasize that E0 represents the thermodynamic potential defining maximum electrochemical efficiency, A governs activation polarization linked to electrode reaction kinetics, I0 reflects the intrinsic electro-catalytic activity of the electrodes, B characterizes diffusion-controlled concentration losses, Il defines the limiting reactant-transport capacity, and Rohm quantifies ionic/electronic resistances across the electrolyte and interconnects. The consistency of these parameters across pressure conditions indicates that the proposed BWO framework captures realistic electrochemical behavior rather than over-fitting numerical data.

4.3 Parity analysis and robustness

Figure 7 (parity plot) illustrates that BWO-based predictions remain within a ±5% error band for all measured pressure conditions. This demonstrates the robustness of the algorithm across a wide operating envelope. To further validate generalizability, a cross-validation experiment was conducted. Model parameters were estimated using experimental data at 1 , 2 , and 4 atm, after which the fitted parameters were applied to predict the polarization behavior at the unseen 5 atm condition. The results confirmed excellent predictive ability, with deviations ≤ 6% in both V–I and P–I polarization curves. This provides strong evidence that the six estimated parameters (E0, A, I0, B, IL, and Rohm) do not merely overfit the dataset but instead capture the intrinsic physical behavior of the SOFC system.

FIGURE 7

Comparison with international literature:

  • reported ∼8% error using grass fibrous root optimization without cross-validation;

  • achieved ∼6–7% error with marine predator algorithm.

The BWO framework, validated through unseen test prediction at 5 atm, was found to maintain an error of ≤6%, thereby confirming its superior accuracy and generalizability. The proposed framework is therefore not only accurate but also transferable to operating scenarios beyond the training dataset, effectively addressing a key concern in parameter estimation research.

4.4 Longevity and stability justification

From a practical standpoint, the accurate optimization of Rohm and Il mitigates local overheating and fuel-starvation zones in the anode, while improved estimation of A and I0 ensures balanced current distribution during transient loading. These effects collectively reduce electrode delamination risks and enhance stack reliability under elevated-pressure operation. Accurate parameter estimation minimizes the deviation between predicted and actual operating points. This allows more precise load management and thermal control, preventing hotspots and fuel starvation in the anode. By ensuring correct estimation of ohmic resistance (Ω·cm2) and limiting current density, degradation processes (electrolyte stress and anode re-oxidation) are slowed. The literature shows that improved control based on accurate modeling can extend SOFC lifetime by 15%–20% (; ). These findings substantiate the conclusion that reduced model–experiment mismatch ensures smoother operation over extended durations. The contribution to system longevity is therefore not hypothetical but directly associated with minimized deviation and enhanced predictive stability.

4.5 Computational efficiency

Figure 8 presents average computational times across 20 runs. BWO consistently achieved convergence within ∼3.7 s, compared to ∼6–8 s for PSO, ∼7 s for GWO, and ∼9 s for WOA. This demonstrates not only accuracy but also real-time feasibility—essential for online SOFC monitoring.

FIGURE 8

All computations were executed on an Intel Core i7 (12th Gen, 3.2 GHz) workstation with 16 GB RAM using MATLAB R2023a. Each full BWO run (2000 iterations, 200 agents) required ∼3.7 s CPU time. For large-scale SOFC stacks, the computational demand grows nearly linearly with the number of cells; therefore, high-performance or parallel implementations are recommended for real-time or multi-stack applications.

4.6 Comparative analysis versus published literature

Table 7 provides a comparative summary of reported optimization approaches for SOFC parameter estimation. Previous studies, such as and , achieved prediction errors of 6%–8% with moderate convergence times, while demonstrated ∼6% error under temperature-dependent conditions. In contrast, the present BWO-based framework achieves ≤ 5% error with significantly faster convergence (∼3.7 s), thereby establishing a new benchmark for pressure-dependent SOFC modeling. Its enhanced accuracy and robustness highlight the framework’s potential to improve stack stability and extend operational lifetime.

TABLE 7

Study (Year)Method/AlgorithmOperating conditionReported error/DeviationComputational efficiencyKey remarks
Genetic algorithms (GA)Transient equivalent model∼10–12%ModerateEarly application; limited to simplified transient circuit models.
Radial movement optimizationSteady state, small-scale SOFC∼8–10%∼6–7 s convergenceReasonable accuracy but slow convergence.
Grass fibrous root optimizationPlanar SOFC, transient models∼7–8%Moderate (∼6–7 s)Convergence improved but accuracy limited.
Marine predator algorithmStatic and dynamic SOFC∼6–7%∼6–8 sGood robustness, but deviations >5%.
Henry gas solubility optimization (HGSO)Temperature-dependent SOFC∼6%∼6 sNovel temperature modeling; no pressure-dependent studies.
Present Study (2025)Black widow optimization (BWO)Pressure-dependent SOFC (1–5 atm, 1027 K)≤5% error (MSE = 0.52 at 5 atm)3.7 s convergence (40%–60% faster)First to model SOFC parameters under varying pressures; superior accuracy, convergence, and robustness demonstrated.

Comparative analysis of optimization algorithms applied for SOFC parameter estimation.

5 Conclusion

This study provides a comprehensive evaluation of the black widow optimization (BWO) algorithm for SOFC parameter estimation under pressure-dependent operating conditions. The major conclusions follow.

5.1 Key findings

  • First application of BWO to SOFC modeling under varying pressure conditions, extending beyond temperature-based optimization methods ().

  • Accurate parameter estimation achieved with MSE as low as 0.52 at 5 atm, consistently outperforming PSO, GWO, and WOA in both accuracy and convergence speed.

  • Robust predictive capability, with model–experiment deviations within ±5% across all operating conditions, demonstrating strong reliability.

  • Computational efficiency enhanced: BWO converged in ∼3.7 s on average, reducing processing time by ∼40–60% compared to competing algorithms.

Future research may extend this work by integrating the BWO framework with AI-based predictive models to enable adaptive and real-time control of SOFC systems. Validation under dynamic operating environments, including transient load changes and fuel composition variations, will further strengthen its applicability. In addition, coupling the framework with hybrid renewable energy systems could enhance efficiency, stability, and scalability in practical deployments. These directions will build upon the present study, providing a pathway toward advancing SOFC technology in clean energy applications and contributing to global efforts for sustainable and efficient energy solutions.

Statements

Data availability statement

No external dataset was used or generated in this study; all data supporting the findings are presented within the manuscript.

Author contributions

PS: Conceptualization, Data curation, Formal Analysis, Investigation, Software, Visualization, Writing – original draft. AS: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Project administration, Software, Writing – original draft. YK: Conceptualization, Formal Analysis, Investigation, Methodology, Project administration, Validation, Writing – original draft. RR: Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Visualization, Writing – review and editing. PB: Data curation, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Writing – review and editing. AG: Data curation, Formal Analysis, Investigation, Methodology, Project administration, Resources, Supervision, Writing – review and editing.

Funding

The author(s) declare that no financial support was received for the research and/or publication of this article.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Glossary

  • ENernst

    Thermodynamic potential [V]

  • Vohm

    Ohmic voltage drop [V]

  • Vconc

    Concentration voltage drop [V]

  • Vact

    Activation voltage drop [V]

  • Partial pressures of hydrogen [bar]

  • Partial pressures of oxygen [bar]

  • Partial pressures of water [bar]

  • T

    Temperature [K]

  • Load current density [mA/cm2]

  • Exchange current density [mA/cm2]

  • Limiting current density [mA/cm2]

Abbreviation

  • SOFC

    Solid oxide fuel cell

  • BWO

    Black widow optimization

  • MSE

    Mean squared error

  • PSO

    Particle swarm optimization

  • GWO

    Gray wolf optimization

  • WOA

    Whale optimization algorithm

  • V-I

    Voltage-current

  • P-I

    Power-current

  • EIS

    Electrochemical impedance spectroscopy

  • PEMFC

    Proton exchange membrane fuel cell

  • CR

    Cannibalism rate

  • MR

    Mutation rate

  • PP

    Procreation probability

  • ECSA

    Electrochemical surface area

  • IL

    Limiting current density

  • Ncell

    Number of cells

  • F

    Faraday constant

  • R

    Universal gas constant

  • atm

    Atmosphere (pressure unit)

  • Kilo-ohm

  • mA/cm2

    Milliampere per square centimeter

  • K

    Kelvin (temperature unit)

  • AI

    Artificial intelligence

  • Eq.

    Equation

References

  • 1

    AmiriA.VijayP.TadéM. O.AhmedK.IngramG. D.PareekV.et al (2016). Planar SOFC system modelling and simulation including a 3D stack module. Int. J. Hydrogen Energy41 (4), 29192930. 10.1016/j.ijhydene.2015.12.076

  • 2

    AnyenyaG. A.SullivanN. P.BraunR. (2017). Modeling and simulation of a novel 4.5 kW multi-stack solid-oxide fuel cell prototype assembly for combined heat and power. Energy Convers. Manag.140, 247259. 10.1016/j.enconman.2017.02.071

  • 3

    El-HayE.El-HameedM.El-FerganyA. (2019). Optimized parameters of SOFC for steady state and transient simulations using interior search algorithm. Energy166, 451461. 10.1016/j.energy.2018.10.038

  • 4

    FallahM.MahmoudiS. M. S.YariM. (2018). A comparative advanced exergy analysis for a solid oxide fuel cell using the engineering and modified hybrid methods. Energy Convers. Manag.168, 576587. 10.1016/j.enconman.2018.04.114

  • 5

    FanX.SunH.YuanZ.LiZ.ShiR.RazmjooyN. (2020). Multi-objective optimization for the proper selection of the best heat pump technology in a fuel cell-heat pump micro-CHP system. Energy Rep.6, 325335. 10.1016/j.egyr.2020.01.009

  • 6

    FangQ.BlumL.PetersR.PeksenM.BatfalskyP.StoltenD. (2015). SOFC stack performance under high fuel utilization. Int. J. Hydrogen Energy40 (2), 11281136. 10.1016/j.ijhydene.2014.11.094

  • 7

    HayyolalamV.KazemA. A. P. (2020). Black widow optimization algorithm: a novel meta-heuristic approach for solving engineering optimization problems. Eng. Appl. Artif. Intell.87, 103249. 10.1016/j.engappai.2019.103249

  • 8

    KangS.AhnK.-Y. (2017). Dynamic modeling of solid oxide fuel cell and engine hybrid system for distributed power generation. Appl. Energy195, 10861099. 10.1016/j.apenergy.2017.03.077

  • 9

    KhattarN.SidhuJ.SinghJ. (2019). Toward energy-efficient cloud computing: a survey of dynamic power management and heuristics-based optimization techniques. J. Supercomput.75, 47504810. 10.1007/s11227-019-02764-2

  • 10

    LeeK.KangS.AhnK.-Y. (2017). Development of a highly efficient solid oxide fuel cell system. Appl. Energy205, 822833. 10.1016/j.apenergy.2017.08.070

  • 11

    MirM.ShafieezadehM.HeidariM. A.GhadimiN. (2020). Application of hybrid forecast engine based intelligent algorithm and feature selection for wind signal prediction. Evol. Intell.11 (4), 559573. 10.1007/s12530-019-09271-y

  • 12

    MirjaliliS.LewisA. (2016). The whale optimization algorithm. Adv. Eng. Softw.95, 5167. 10.1016/j.advengsoft.2016.01.008

  • 13

    NassefA. M.FathyA.SayedE. T.AbdelkareemM. A.RezkH.TanveerW. H.et al (2019). Maximizing SOFC performance through optimal parameters identification by modern optimization algorithms. Renew. Energy138, 458464. 10.1016/j.renene.2019.01.072

  • 14

    NeratM. (2017). Modeling and analysis of short-period transient response of a single, planar, anode supported, solid oxide fuel cell during load variations. Energy138, 728738. 10.1016/j.energy.2017.07.133

  • 15

    PutilovL. P.DeminA. K.TsidilkovskiV. I.TsiakarasP. (2019). Theoretical modeling of the gas humidification effect on the characteristics of proton ceramic fuel cells. Appl. Energy242, 14481459. 10.1016/j.apenergy.2019.03.096

  • 16

    PutilovL. P.TsidilkovskiV. I.DeminA. K. (2020). Revealing the effect of the cell voltage and external conditions on the characteristics of protonic ceramic fuel cells. J. Mater. Chem. A8 (25), 1264112656. 10.1039/d0ta03935g

  • 17

    ShiH.LiJ.ZafettiN. (2020). RETRACTED: new optimized technique for unknown parameters selection of SOFC using converged Grass Fibrous Root Optimization Algorithm. Energy Rep.6, 14281437. 10.1016/j.egyr.2020.05.024

  • 18

    SinghP.SandhuA. (2022). Optimal parameter extraction of proton exchange membrane fuel cell using Henry gas solubility optimization. Int. J. Energy Res.46 (13), 1821218224.

  • 19

    SinghP.SandhuA. (2023). Parameter identification of a solid oxide fuel cell model under different temperature conditions. Chem. Eng. & Technol.46, 23752385. 10.1002/ceat.202300009

  • 20

    SinglaM. K.GuptaJ.AlsharifM. H.JahidA. (2023a). Optimizing integration of fuel cell technology in renewable energy-based microgrids for sustainable and cost-effective energy. Energies16 (11), 4482. 10.3390/en16114482

  • 21

    SinglaM. K.GuptaJ.NijhawanP.AlsharifM. H.KimM.-K. (2023b). Sustainable development of fuel cell using enhanced weighted mean of vectors algorithm. Heliyon9 (3), e14578. 10.1016/j.heliyon.2023.e14578

  • 22

    SinglaM. K.GuptaJ.SinghB.NijhawanP.AbdelazizA. Y.El-ShahatA. (2023c). Parameter estimation of fuel cells using a hybrid optimization algorithm. Sustainability15 (8), 6676. 10.3390/su15086676

  • 23

    VirkarA.WilliamsM. C.SinghalS. (2007). Concepts for ultra-high power density solid oxide fuel cells. ECS Trans.5 (1), 401421. 10.1149/1.2729021

  • 24

    WangC.NehrirM. H. (2007). A physically based dynamic model for solid oxide fuel cells. IEEE Trans. Energy Convers.22 (4), 887897. 10.1109/tec.2007.895468

  • 25

    WangK.HisselD.PéraM.-C.SteinerN.MarraD.SorrentinoM.et al (2011). A review on solid oxide fuel cell models. Int. J. Hydrogen Energy36 (12), 72127228. 10.1016/j.ijhydene.2011.03.051

  • 26

    WuY.GaoH. (2006). Optimization of fuel cell and supercapacitor for fuel-cell electric vehicles. IEEE Trans. Veh. Technol.55 (6), 17481755. 10.1109/tvt.2006.883764

  • 27

    XiongG.ZhangJ.ShiD.YuanX. (2020). A simplified competitive swarm optimizer for parameter identification of solid oxide fuel cells. Energy Convers. Manag.203, 112204. 10.1016/j.enconman.2019.112204

  • 28

    YousriD.HasanienH. M.FathyA. (2021). Parameters identification of solid oxide fuel cell for static and dynamic simulation using comprehensive learning dynamic multi-swarm marine predators algorithm. Energy Convers. Manag.228, 113692. 10.1016/j.enconman.2020.113692

  • 29

    YuD.ZhangT.HeG.NojavanS.JermsittiparsertK.GhadimiN. (2020). RETRACTED: energy management of wind-PV-storage-grid based large electricity consumer using robust optimization technique. J. Energy Storage27, 101054. 10.1016/j.est.2019.101054

  • 30

    ZhouJ.LiuQ.ZhangL.PanZ.ChanS. H. (2016). Influence of pore former on electrochemical performance of fuel-electrode supported SOFCs manufactured by aqueous-based tape-casting. Energy115, 149154. 10.1016/j.energy.2016.08.093

  • 31

    ZhuL.ZhangL.VirkarA. V. (2015). A parametric model for solid oxide fuel cells based on measurements made on cell materials and components. J. Power Sources291, 138155. 10.1016/j.jpowsour.2015.04.155

Summary

Keywords

solid oxide fuel cell, hydrogen, optimization algorithms, mathematical modeling, black widow optimization

Citation

Singh P, Sandhu A, Khan Y, Raman R, Barmavatu P and Garg A (2025) Energy-efficient parameter estimation of solid oxide fuel cells under varying pressure conditions using the black widow optimization algorithm. Front. Energy Res. 13:1659232. doi: 10.3389/fenrg.2025.1659232

Received

03 July 2025

Revised

05 October 2025

Accepted

06 October 2025

Published

24 November 2025

Volume

13 - 2025

Edited by

A. Brouzgou, University of Thessaly, Greece

Reviewed by

Dimitris Ipsakis, Technical University of Crete, Greece

Zhiqiang Niu, Loughborough University, United Kingdom

Updates

Copyright

*Correspondence: Roshan Raman,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics