Next Article in Journal
Potential for Biogas Production from Water Hyacinth and Banana Peels: A Case Study of Substrates Harvested from Lomé, Togo
Previous Article in Journal
Steam Reforming of Tar Impurities from Biomass Gasification with Ni-Co/Mg(Al)O Catalysts—Operating Parameter Effects
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Real Case: How to Combine Polarization Curve and EIS Techniques to Identify Problematic Cells in a Commercial PEM Stack

by
Guillermo Gómez
*,
Pilar Argumosa
and
Jesús Maellas
Energy and Environment Department, National Institute of Aerospace Technology (INTA), E-28850 Madrid, Spain
*
Author to whom correspondence should be addressed.
Fuels 2024, 5(3), 476-493; https://doi.org/10.3390/fuels5030026
Submission received: 17 July 2024 / Revised: 11 August 2024 / Accepted: 27 August 2024 / Published: 2 September 2024

Abstract

:
Nowadays, the mobility sector is assessing different technologies to substitute the internal combustion engines in order to reduce its CO2 emissions; one of these possible alternatives is the Polymer Electrolyte Membrane (PEM) fuel cell. So, the development of non-destructive diagnostic tools that could identify defective cells and/or any malfunctioning behavior and can be easily embarked on in any vehicle will expand the durability of PEM fuel cells, improve their performance, and enable them to carry out predictive maintenance. In this research, we use an in-house developed methodology that combines the polarization curve and electrochemical impedance spectroscopy (EIS) techniques to characterize different cells of a commercial PEM stack, identifying malfunctioning ones.

1. Introduction

Currently, the global mobility sector is sifting to electrification in order to mitigate the environmental effects partially caused by internal combustion engine vehicles, such as global warming and/or air pollution [1,2,3]. In this new context, different electric vehicle technologies have appeared, which can basically be classified into five categories: battery electric vehicles (BEVs), plug-in hybrid electric vehicles (PHEVs), hybrid vehicles (HVs), fuel cell electric vehicles (FCEVs), and extended-range electric vehicles (ER-EVs) [4], but only the BEVs and FCEVs are fully electric, as the other technologies combine electric and conventional combustible engines.
So, the fuel cells together with batteries are called to play a leading role to address the decarbonization of this sector [5], especially if hydrogen is produced from renewable energies [6,7,8,9]. And among the fuel cell technologies, the PEM fuel cells are the most suitable for mobility applications [10,11,12] because of their characteristics, operating temperature, easy scale-up, all solid components, fast dynamic response, high energy conversion efficiency, and power density.
Therefore, [13,14,15] development of non-intrusive and easily embarked diagnostic tools that could assess the performance of stack cells during real operation, identifying anomalous behavior, will ease the transition from combustion engines to this more friendly environmental technology, as these tools will enable to carry out predictive maintenance, expanding the durability of PEM fuel cells, and improving their performance.
According to [16], fuel cell failure diagnosis tools can be classified into two categories: off-line and on-line techniques.
The off-line category refers to those methodologies in which the analyzed fuel cell is operated under some conditions that do not necessarily match with the normal operating fuel cell conditions. This category includes electrochemical characterization techniques such as cyclic voltage (CV), galvanostatic charge (GSC), electrochemical impedance spectroscopy (EIS), etc.
On the other hand, the on-line diagnosis methodologies are able to infer the fuel cell state of health when it is operated in their intended application. In Refs. [13,16,17] these techniques use some mathematical algorithms to process in real time different sensor readings into a diagnosis. According to the used algorithms, those techniques can be classified into: model-based, data-driven, signal-based, and hybrid methods. Where the model-based methods consist in a differential analysis between real-time measurements and healthy fuel cell model outputs.
In this research, we use an in-house developed methodology, which is described in [18], to detect malfunctioning cells of a commercial stack. This method is based on a five-parameter equation, which is able to model the voltage of a cell in function of current. And, using the results of EIS and polarization curve tests to determine these five parameters for this cell, we transpose the test conditions (temperature, pressure, relative humidity, etc.) to the model and obtain for this cell its open circuit voltage (OCV), internal resistance, and an expression for its activation losses. And therefore, comparing these results of a cell with the results for a healthy one, we can try to determine if the assessed cell shows flaws.
According to the previous paragraphs, our methodology can be classified as model-based, but on the other hand, it cannot be considered on-line because the used EIS technique is sensitive to electric noise, so it should not be performed during transient events, and the normal driving activity is plenty of them. Even though the EIS technique used in this research is an off-line one, currently there are some efforts to develop on-line diagnosis tools based on the EIS technique [19] that may be used to implement the fuel cell analysis methodology described in this article.
The choice of the selected techniques to characterize a cell was based on the fact that the necessary equipment to perform tests should be easily embarked upon and not complex, and both techniques fulfill this requirement, and additionally, they are common fuel cell electrochemical characterization techniques and used to diagnostic (PEM) fuel cells and stacks [20,21,22,23,24,25,26].
The polarization curve test [27,28,29] probably is the most basic PEM fuel cell characterization test. In this test, the voltage of the fuel cell is measured after stationary behavior is reached, when the fuel cell is forced to work at specific current and test conditions (temperature, pressure, gas relative humidity, gas stoichiometry, etc.). Synthesizing the general performance of the fuel cell into a voltage signal can be used to determine the fuel cell state of health, but on the other hand, this integration process can blur the causes of the fuel cell behavior.
Electrochemistry Impedance Spectroscopy [30,31,32] is a technique that tries to cause resonance with subprocesses of an electrochemical system and study its response, obtaining extra-information that is not accessible for non-resonance techniques such as polarization curve tests. To this end, an electrical stimulus (a known tiny periodic voltage or current signal) is applied to electrodes at different frequencies, obtaining a response (a resulting current or voltage). And this extra information makes the EIS technique suitable to detect some internal flaws as hydrogen crossover flow, study the effects of operating conditions on the fuel cell internal resistance, and determine the hydric conditions of fuel cells [33,34,35,36].
The aim of this research is to test the sensitivity of the methodology described in [18] to identify flaws in a stack. To this end, we used a commercial PEM stack, in which we previously detected a gas leakage in the seal of one cell when the stack temperature dropped below 15 °C. Then, fixing the ambient temperature at 20 °C and applying our method to the cell with flaws and on the other two normal cells, and comparing results, we observed differences between the internal resistance and the expressions for activation losses of the unhealthy and healthy cells.

2. Materials and Methods

In this section, we will briefly describe the mathematical expression used to model the cell behavior, the polarization curve and EIS test procedure, and the test setup with a description of the main equipment.

2.1. Mathematical Model Description

To model the cell behavior, we use Equation (1)
V = V 0 C t e I I b + d r · I
where V is the cell voltage; I is the current; and V 0 , C t e , b , d and r are the fitting parameters.
To obtain b and d , we use the equivalent circuit of Figure 1 to fit the EIS results, where R Ω is the ohmic resistance, R t c is charge transfer resistance, and CPE = constant phase element. Now, we plot R t c as function of current and fit the resulting curve by Equation (2), obtaining b and d. V is a fitting parameter of Equation (2).
R t c = V I b + d
Then, substituting b and d in Equation (1) and using this expression to fit the polarization curve, we achieve the rest of parameters, obtaining the following correspondences V 0 = open circuit voltage (OCV), r = internal resistance, and C t e I b + d mathematical expression for activation losses.
It is convenient to highlight that whenever we have used this methodology, the fitting parameters of Equation (1) have always had physical meaning if previously we have obtained b and d through Equation (2). Conversely, if we use Equation (1) without restrictions, the fitting is slightly better, but in this case, the fitting parameters have no physical sense.
A more comprehensive explanation of the model and the development of the mathematical expressions are included in [18].

2.2. Test Setup and Main Equipment Description

In this research, we tested an used commercial fuel cell system of 2500 W, which has a stack composed of 80 cells connected in series. It is an autonomous system that, for operation, only requires an external hydrogen supply, and it is designed to power applications such as aerial drones, ground robots, portable power packs, etc.
We selected this used unit because a tiny H2 leakage was observed in the vicinity of Cell 79 when the stack was tempered below 15 °C before starting up, providing us the opportunity to check if the methodology described in [18] was able to detect any anomalous behavior in Cell 79. Then, to this end, we applied this methodology to Cell 2, Cell 40, and Cell 79 with the aim of performing a differential analysis.
In Figure 2, there is a scheme of the test setup used to perform the tests. It basically consisted of a H2 supply system, an electronic load to simulate the power demand, a power supply to feed the Fuel Cell System BOP, a control station to command the Fuel Cell System and the electronic load, and EIS equipment. This EIS equipment was composed of the following:
  • Potentiostat: Solartron SI1287 (Electrochemical Interface).
  • Frequency analyzer: Solartron SI1255 HF (Frequency Response Analyzer).
Figure 2. Test setup.
Figure 2. Test setup.
Fuels 05 00026 g002
To analyze the EIS results, we used the Zplot and Zview software version 3.5h by Scribner.
These tests basically consisted of forcing the fuel cell system to work at some currents, and after fuel cell voltage stabilization, we launched a frequency profile, obtaining the impedance response for each of the selected cells. To perform these tests, the stack was tempered at 20 °C.

3. Results

In this section, we analyze the results of applying the methodology described in [18] on Cell 2, Cell 40, and Cell 79.

3.1. Test Results for Cell 79

Figure 3 plots the polarization and EIS curves obtained for Cell 79.
To better appreciate the EIS curves, we represent them individually and their fitting curve in Figure 4. To fit the EIS curves, we use the equivalent circuit of Figure 1.
If we analyze Figure 4, we can observe the existence of anomalous points. To discriminate these points, we calculate the distance between the measured and fitted values for each point through Equation (3).
d i s = Z m e a s i Z f i t t i max Z m e a s min Z m e a s 2 + Z m e a s i Z f i t t i max Z m e a s min Z m e a s 2
where Z m e a s i is the measured impedance real part of i-point; Z m e a s i is the measured impedance imaginary part of i-point; Z f i t t i is the fitted impedance real part of i-point; Z f i t t i is the fitted impedance imaginary part of i-point; max ( Z m e a s ) is the maximum value of all measured impedances real parts; min ( Z m e a s ) is the minimum value of all measured impedances real parts; max ( Z m e a s ) is the maximum value of all measured impedances imaginary parts; min ( Z m e a s ) is the minimum value of all measured impedances imaginary parts.
Then, define the outlier as the point that is 1.25 from the mean distance in Equation (4), where N is the total number of points for an EIS curve.
d i s > 1.25 · 1 N i = 1 N d i s
And suppressing the points that accomplish this condition, we obtain Figure 5. This criterion is not applied to the EIS curve at I = 0 A.
Now, we use the equivalent circuit of Figure 1 to fit the EIS curves of Figure 5, obtaining the following parameterization, Table 1.
Now we represent Rtc against current and fit the resulting curve by Equation (2), obtaining the following:
R t c = 0.06034 I 0.9084 + 0.07066 e r r r e l = 27 %
Defining the mean relative error by Equation (6)
< e r r r e l > = 1 N i = 1 N x i _ m e a s s x i _ f i t t x i _ m e a s s
where x i _ m e a s s are measured values, x i _ f i t t are fitting values, and N is the number of measurements.
In Figure 6, we represent the R t c versus current and its fitting curve; the left picture corresponds to the whole current range, while the right picture only shows the curves from I = 0.5 A to I = 50 A to better appreciate the discrepancies between measurement and fitting curve, observing that Equation (5) only describes correctly the measurements for low currents below 10 A.
Then, substituting b = 0.9084 and d = 0.07066 in Equation (1) and fitting the polarization curve of Figure 3, we obtain Equation (7):
V = 0.9585 0.13133 I 0.9084 + 0.07066 · I 0.00253 · I   ;   e r r r e l = 0.009
Now, in Figure 7, we plot together the polarization curve of Figure 3 and its fitting curve and observe that they are practically coincident, with a mean relative error lower than 1%.
Therefore, we have the following characteristic parameterization for Cell 79:
  • Open circuit voltage ( V 0 ) : 0.9585 V.
  • Internal ohmic resistance ( r ) : 0.00253 Ω.
  • Mathematical expression of activation voltage loss: 0.13133 I 0.9084 + 0.07066 .

3.2. Test Results for Cell 40

In Figure 8, we represent the polarization and EIS curves obtained for Cell 40.
As to Cell 2, we represent individually the EIS and their fitting curves in Figure 9. To this end, we use the equivalent circuit of Figure 1.
Now, using Equations (3) and (4) to suppress the anomalous points from the analysis and fitting the resulting curves by the equivalent circuit of Figure 1, we obtained the following parameters for the filtered EIS curves, Table 2.
Now we represent Rtc against current and fit the resulting curve by Equation (2), obtaining
R t c = 0.06365 I 0.83056 + 0.05473 e r r r e l = 18 %
And proceeding as Section 3.2, we fit the polarization curve of Figure 8, obtaining
V = 0.96483 0.14259 I 0.83056 + 0.05473 · I 0.00118 · I   ;   e r r r e l = 0.004
Then, if we plot together the polarization and the fitting curves in Figure 10, we can observe they can be considered overlapped:
Therefore, the characteristic parameters of Cell 40 are as follows:
  • Open circuit voltage ( V 0 ) : 0.96483 V.
  • Internal ohmic resistance ( r ) : 0.00118 Ω.
  • Mathematical expression of activation voltage loss: 0.14259 I 0.83056 + 0.05473 .

3.3. Test Result for Cell 2

Figure 11 shows the polarization and EIS curves for Cell 2.
Now, as usual, we represent individually the EIS curves and their fitting curves in Figure 12. To fit the EIS curves, we use the equivalent circuit of Figure 1.
As before, we remove the anomalous points from the analysis through Equations (3) and (4), and using the equivalent circuit of Figure 1, we fit the filtered EIS curves, obtaining the parameters of Table 3.
As to the other cells, we represent Rtc against current and fit this curve by Equation (2), resulting in
R t c = 0.06416 I 0.86177 + 0.05516 e r r r e l = 20 %
Now proceeding as usual, we fit the polarization curve of Figure 11, obtaining
V = 0.9631 0.14453 I 0.86177 + 0.05516 · I 0.00165 · I   ;   e r r r e l = 0.006
In Figure 13, we plot together the polarization curve of Figure 11 with its fitting curve, obtaining a mean relative error between curves of 0.006.
Therefore, according to the obtained parametrization, Cell 2 has the following characteristics:
  • Open circuit voltage ( V 0 ) : 0.9631 V
  • Internal ohmic resistance ( r ) : 0.00165 Ω
  • Mathematical expression of activation voltage loss: 0.14453 I 0.86177 + 0.05516

4. Discussion

In Figure 14, we plot together the polarization curves of Cell 2, 40, and 79.
If we analyze Figure 14, we can observe that they show similar behavior; their differences are slightly with a relative mean standard deviation <2%. So, we can conclude the polarization curve test is not able to detect any flaws in Cell 79; even Cell 70 shows the best performance.
On the other hand, if we compare the values of open circuit voltage (OCV), internal resistance, and activation loss expression for each cell, we can observe discrepancies between cell behavior.
Before continuing with discussion, it is convenient to notice that when we use Equation (1) to fit the polarization curves, the fitting parameters V 0 , C t e , and r are completely free, and we have not imposed any restrictions on them.

4.1. Open Circuit Voltage ( V 0 ) Discussion

To analyze the V0, we plot together the measured and fitted open circuit voltages for each cell, obtaining Figure 15:
Then, it is easy to observe that the measured and fitted V0 can be considered coincident, so the methodology described in [18] provides very good results for open circuit voltage.

4.2. Internal Ohmic Resistance ( r ) Discussion

In this case we do not have a direct measure of the cell internal ohmic resistance; the R (Ω) of the EIS test includes the effect of the contact resistance between electrodes and stack. Then, taking into consideration that the internal ohmic resistance is predominant in the linear part of the polarization curve, we use the next mathematical expression to calculate an estimation of it for each cell.
r e s t i m a t e d = V 10 A V 50 A 50   A 10   A
where V 10 A is the cell voltage at 10 A and V 50 A is the cell voltage at 50 A. Now in Figure 16, we plot the estimated and fitted internal ohmic resistance for each cell.
In this case we can observe the estimated and fitted internal ohmic resistances are not coincident, but they have the same magnitude order, and the estimated resistances are always higher than the fitted ones, as it would be expected because in the estimated resistance is included the effect of the activation voltage loss, while in the fitted ones this effect has been suppressed.
But the most interesting results come out of the fitting internal resistances in comparison. We can observe that the lower value corresponds to Cell 40, and it makes sense as Cell 40 is internal, and it would be expected that the contact between internal cells was more intimate than external ones. But the internal resistance of Cell 79 is excessive; it multiplies by 1.53 the internal resistance of Cell 2 and by 2.14 the internal resistance of Cell 40, and this fact is consistent with a H2 leakage in Cell 79; bigger resistances could imply the existence of not correct tightened cells, and this could cause H2 leakage.

4.3. Activation Voltage Loss ( A V L ) Discussion

According to dimensional analysis, the kinetic voltage loss is an ohmic resistance, and it should be related to the transfer charge resistance ( R t c ). To check this relation, and taking into consideration that the effect of the AVL on the cell voltage is more determinant for low currents, we compare the AVL and R t c values for the [0–5] A range, Table 4 and Figure 17.
Now, comparing the AVL and R t c values, we observe they are different, but they show similar. Therefore, the AVL mathematical expression reproduces the behavior of transfer charge resistance.
But the most interesting result is that the behavior of the AVL for Cell 79 is clearly different to Cell 40 and Cell 2, while the AVL for Cell 40 and Cell 2 can be considered coincident. Cell 79 shows the best performance as it shows the lowest resistance, or, in other words, Cell 79 is more activated than Cell 2 and Cell 40. Then, if we calculate the AVL value of all cells when I = 0 A, we obtain
A V L c e l l 2 _ I = 0 A = 2.62   Ω ;   A V L c e l l 40 _ I = 0 A = 2.61   Ω ;   A V L c e l l 79 _ I = 0 A = 1.86   Ω
Now, taking into consideration that the AVL value decreases with current, we can calculate at which current the AVL value of Cells 2 and 40 is equal to 1.86   Ω .
A V L c e l l 2 = 0.14453 I 0.86177 + 0.05516 = 1.86 I = 0.012   A
A V L c e l l 40 = 0.14259 I 0.83056 + 0.05473 = 1.86 I = 0.010   A
Then, according to these results and our interpretation, the standby state of Cell 79 corresponds to the state of Cell 2 and Cell 40 when they are producing a current of I = 0.012 A and I = 0.010 A, respectively. This may be interpreted as the Cell 79 is more activated than the rest of the cells to compensate for the voltage loss caused by the internal resistance, and this extra work probably accelerates the degradation processes of Cell 79. If this interpretation is correct, this methodology may be a useful tool for predictive maintenance, but due to the current stage of our research, we are not in a position to forecast the durability of any defective cell. To this end, we would need to study how the model parameters evolve according to how the stack is used for its intended purpose, but this task will be faced in future research.
In light of these results, the usefulness of this methodology could be extended to other topics of fuel cell systems, such as the implementation of intelligent control logics in this kind of system. According to [37], a fuel cell system is a compound of several subsystems, and implementing an intelligent control would be desirable to obtain an accurate and reliable fuel cell model that could be able to provide the static and dynamic response of a FC in real operation. Currently, to obtain this model, there would be two approaches: physical and empirical models. The physical models, which try to model the physical phenomena of a fuel cell, are complex and require large computational resources; on the other hand, the empirical models are simpler, although they can provide accurate behaviors, their parameters could not be interpreted from a physical point of view.
In this research, we have developed an empirical and simple mathematical fuel cell model whose parameters have physical meaning. But currently our methodology is not able to reproduce the dynamic behavior of a fuel cell, as the used EIS technique does not provide correct parameterization in non-stationary tests. So the development of EIS techniques that could be used during dynamic testing could make our methodology a promising candidate for fuel cell control models.
Additionally, if we analyze the fuel cell model implemented through Equation (1), we find that probably its most remarkable characteristic is that the fitting parameters have physical meaning, so they could be useful to obtain an insight into a fuel cell internal phenomena. Conforming to [38] and [39], one of the biggest concerns with PEM Fuel Cell is improving the kinetics of cathode catalysts because the reaction kinetics of ORR are much slower than the HOR. So currently there are some efforts to obtain high-performance electrocatalysts for ORR with a reduced quantity of Pt. In this context, our research could play a role, as the AVL term of Equation (1) may be useful to define a metric to quantify the performance for different ORR electrocatalyst solutions.

5. Conclusions

Then, summarizing the most remarkable results of this research, we conclude that:
According to this research, the polarization curve test is not able to detect this type of flaw, as Cell 79 shows the best performance (Figure 14). But on the other hand, the methodology described in [18] is sensitive, the internal ohmic resistance of Cell 79 is clearly bigger, and the behavior described by AVL expression of this cell is different from the others. Therefore, this methodology could be a useful off-line tool for finding incorrectly tightened cells and/or detecting over-activated ones.
Conforming to the results, this technique provides excellent values for OCV and at least acceptable values for r and the AVL expression reproduces correctly the behavior of R t c for low currents below 5 A.
Equation (1) consistently provides good approximations for the polarization curves, but if the parameters b and d take the values calculated using the EIS results, the fitting parameters of Equation (1) have always physical sense. According to the article, the following are true:
  • V 0 : Open circuit voltage.
  • r : Internal ohmic resistance.
  • C t e I b + d : Mathematical expression for activation voltage loss (AVL).
The physical meaning of AVL parameters is more controversial. The Cte parameter can be thought of as a product of a proportional constant (PC) by V , or, in other words
A V L = C t e I b + d = P C · R t c = P C · V I b + d
PC is a proportional constant that provides the effective loss of voltage caused by R t c . The rest of parameters V , b, and d are used to model the R t c dependence with current. According to our interpretation of the fuel cell model represented in Figure 1, the AVL is modeled by the CPE and R t c , and when the fuel cell is operated at OCV, only a tiny residual current would cross the cathode double layer (CPE element), but as the power grows, more current would cross the double layer in parallel pathways, as it is represented in Figure 18.
Then the total charge transfer resistance can be interpreted as the sum of parallel resistances of the ionic current paths that cross the double layer, whose number increases with current:
1 R T = 1 R 1 + 1 R 2 + + 1 R n 1 R T = R 2 · R 3 · · · R n + R 1 · R 3 · · · R n + + R 1 · R 2 · · · R n 1 R 1 · R 2 · · · R n
where RT is the total charge transfer resistance, and Ri is the resistance of i-th ionic current path. Now, considering that R1 = R2 = ⋯ = Rn without loss of generality, we obtain a rational function for total charge transfer resistance that decreases with current:
1 R T = n · R n 1 R n R T = R n R T = R n ( I )
In light of this, V , b, and d are the fitting parameters of a rational function to describe the behavior of R t c with current. Then, according to our interpretation, d would be related to the existence of some tiny residual current that is independent of the power demand, and the parameter b is used to fit the dependence of R t c with current. Because the resistance of all ionic paths has not to be equal, b may be a measure of the double layer homogeneity, representing b = 1 as a homogeneous double layer. It is convenient to notice that the R t c behavior modeling is not a closed issue, and other mathematical prototypes can be checked. We use Equation (2) for its simplicity and because it works.

Author Contributions

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

Funding

This research received no external funding.

Data Availability Statement

No new data were created.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

AVLActivation Voltage Loss
BEVsBattery Electric Vehicles
BOPBalance of Plant
CVCyclic Voltammetry
CPEConstant Phase Element
EISElectrochemical Impedance Spectroscopy
ER-EVsExtended-range Electric Vehicles
FCEVsFuel Cell Electric Vehicles
GSCGalvanostatic Charge
HORHydrogen Reduction Reaction
HVsHybrid Vehicles
OCVOpen Circuit Voltage
ORROxygen Reduction Reaction
PEMPolymer Electrolyte Membrane
RtcCharge transfer resistance (Ω)

References

  1. Mishra, A.; Talreja, P.; Shrivastava, A. Performance Evaluation of Electric Vehicle Stocks: Paving a Way Towards Green Economy. In Global Economic Revolutions: Big Data Governance and Business Analytics for Sustainability, Proceedings of the ICGER 2023. Communications in Computer and Information Science, Sharjah City, United Arab Emirates, 27–28 February 2023; Musleh, M.A., Al-Sartawi, A., Helmy Abd Wahab, M., Hussainey, K., Eds.; Springer: Cham, Switzerland, 2024. [Google Scholar] [CrossRef]
  2. Taljegard, M.; Göransson, L.; Odenberger, M.F.; Johnsson, F. Impacts of electric vehicles on the electricity generation portfolio—A Scandinavian-German case study. Appl. Energy 2019, 235, 1637–1650. [Google Scholar] [CrossRef]
  3. Shaukat, N.; Khan, B. E-Mobility: Transportation Sector in Transition. In Handbook of Climate Change Mitigation and Adaptation; Lackner, M., Sajjadi, B., Chen, W.Y., Eds.; Springer: Cham, Switzerland, 2022. [Google Scholar] [CrossRef]
  4. Sanguesa, J.A.; Torres-Sanz, V.; Garrido, P.; Martinez, F.J.; Marquez-Barja, J.M. A Review on Electric Vehicles: Technologies and Challenges. Smart Cities 2021, 4, 372–404. [Google Scholar] [CrossRef]
  5. Fragiacomo, P.; Piraino, F.; Genovese, M.; Corigliano, O.; De Lorenzo, G. Strategic Overview on Fuel Cell-Based Systems for Mobility and Electrolytic Cells for Hydrogen Production. Procedia Comput. Sci. 2022, 200, 1254–1263. [Google Scholar] [CrossRef]
  6. Burkhardt, J.; Patyk, A.; Tanguy, P.; Retzke, C. Hydrogen mobility from wind energy—A life cycle assessment focusing on the fuel supply. Appl. Energy 2016, 181, 54–64. [Google Scholar] [CrossRef]
  7. Bragatto, T.; Bucarelli, M.A.; Carere, F.; Cavadenti, A.; Santori, F. Optimization of an energy district for fuel cell electric vehicles: Cost scenarios of a real case study on a waste and recycling fleet. Int. J. Hydrogen Energy 2022, 47, 40156–40171. [Google Scholar] [CrossRef]
  8. Oyewole, O.L.; Nwulu, N.I.; Okampo, E.J. Optimal design of hydrogen-based storage with a hybrid renewable energy system considering economic and environmental uncertainties. Energy Convers. Manag. 2024, 300, 117991. [Google Scholar] [CrossRef]
  9. Syré, A.M.; Shyposha, P.; Freisem, L.; Pollak, A.; Göhlich, D. Comparative Life Cycle Assessment of Battery and Fuel Cell Electric Cars, Trucks, and Buses. World Electr. Veh. J. 2024, 15, 114. [Google Scholar] [CrossRef]
  10. Mohammed, H.; Al-Othman, A.; Nancarrow, P.; Tawalbeh, M.; El Haj Assad, M. Direct hydrocarbon fuel cells: A promising technology for improving energy efficiency. Energy 2019, 172, 207–219. [Google Scholar] [CrossRef]
  11. Jayakumar, A.; Chalmers, A.; Lie, T.T. A review of Prospects for adoption of fuel electric vehicles in New Zealand. IET Electr. Syst. Transp. 2017, 7, 259–266. [Google Scholar] [CrossRef]
  12. Liu, P.; Xu, S. A review of low-temperature proton exchange membrane fuel cell degradation caused by repeated freezing start. Int. J. Hydrogen Energy 2023, 48, 8216–8246. [Google Scholar] [CrossRef]
  13. Aubry, J.; Steiner, N.Y.; Morando, S.; Zerhouni, N.; Hissel, D. Fuel cell diagnosis methods for embedded automotive applications. Energy Rep. 2022, 8, 6687–6706. [Google Scholar] [CrossRef]
  14. Cadet, C.; Jemeï, S.; Druart, F.; Hissel, D. Diagnostic tools for PEMFCs: From conception to implementation. Int. J. Hydrogen Energy 2014, 39, 10613–10626. [Google Scholar] [CrossRef]
  15. Pahon, E.; Hissel, D.; Jemei, S.; Steiner, N.Y. Signal-based diagnostic approach to enhance fuel cell durability. J. Power Sources 2021, 23, 230223. [Google Scholar] [CrossRef]
  16. Zhou, S.; Jervis, R. A Review of Polymer Electrolyte Fuel Cells Fault Diagnosis: Progress and Perspectives. Chem. Methods 2024, 4, e202300030. [Google Scholar] [CrossRef]
  17. Zheng, Z.; Petrone, R.; Péra, M.C.; Hissel, D.; Becherif, M.; Pianese, C.; Steiner, N.Y.; Sorrentino, M. A review on non-model based diagnosis methodologies for PEM fuel cell stacks and systems. Int. J. Hydrogen Energy 2013, 38, 8914–8926. [Google Scholar] [CrossRef]
  18. Gómez, G.; Argumosa, P.; Correro, A.; Maellas, J. Proposal of a New Technique to Obtain Some Fuel Cell Internal Parameters Using Polarization Curve Tests and EIS Results. Energies 2021, 14, 7161. [Google Scholar] [CrossRef]
  19. Manganiello, P.; Petrone, G.; Giannattasio, M.; Monmasson, E.; Spagnuolo, G. FPGA implementation of the EIS technique for the on-line diagnosis of fuel-cell systems. In Proceedings of the 2017 IEEE 26th International Symposium on Industrial Electronics (ISIE), Edinburgh, UK, 19–21 June 2017; pp. 981–986. [Google Scholar] [CrossRef]
  20. Wu, J.; Yuan, X.Z.; Wang, H.; Blanco, M.; Martin, J.J.; Zhang, J. Diagnostic tools in PEM fuel cell research: Part I Electrochemical techniques. Int. J. Hydrogen Energy 2008, 33, 1735–1746. [Google Scholar] [CrossRef]
  21. Pivac, I.; Radica, G.; Barbir, F.; Benouioua, D.; Harel, F.; Candusso, D. Deliverable D1.4 Diagnostic Methods for Automotive Fuel Cell Systems, GIANTLEAP Project no 700101. Fuel Cells and Hydrogen 2 Join Undertaking, 29 December 2017. Available online: https://ec.europa.eu/research/participants/documents/downloadPublic?documentIds=080166e5b77c06fa&appId=PPGMS (accessed on 30 June 2024).
  22. Najafi, B.; Bonomi, P.; Casalegno, A.; Rinaldi, F.; Baricci, A. Rapid Fault Diagnosis of PEM Fuel Cells through Optimal Electrochemical Impedance Spectroscopy Tests. Energies 2020, 13, 3643. [Google Scholar] [CrossRef]
  23. Zhang, C.; Li, W.; Hu, M.; Cheng, X.; He, K.; Mao, L. A Comparative Study of Using Polarization Curve Models in Proton Exchange Membrane Fuel Cell Degradation Analysis. Energies 2020, 13, 3759. [Google Scholar] [CrossRef]
  24. Aabid, S.E.; Turpin, C.; Regnier, J.; D’Arbigny, J.; Horde, T.; Pessot, A. Monitoring of Ageing Campaigns of PEM Fuel Cell Stacks Using a Model-Based Method. Int. J. Renew. Energy Res. 2021, 11, 92–100. [Google Scholar]
  25. Hu, Z.; Xu, L.; Li, J.; Gan, Q.; Xu, X.; Song, Z.; Shao, Y.; Ouyang, M. A novel diagnostic methodology for fuel cell stack health: Performance, consistency and uniformity. Energy Convers. Manag. 2019, 185, 611–621. [Google Scholar] [CrossRef]
  26. Pessot, A.; Turpin, C.; Jaafar, A.; Soyez, E.; Rallières, O.; Gager, G.; d’Arbigny, J. Contribution to the modelling of low temperature PEM fuel cell in aeronautical conditions by design of experiments. Math. Comput. Simul. 2019, 158, 179–198. [Google Scholar] [CrossRef]
  27. Malkov, T.; Pilenga, A.; Tsotridis, G.; de Marco, G. EU Harmonised Polarization Curve Test Method for Low Temperature Water Electrolysis. European Commission, Joint Research Centre, Directorate C: Energy, Transport and Climate. © European Union 2018. Available online: https://publications.jrc.ec.europa.eu/repository/bitstream/JRC104045/kjna29182enn_online.pdf (accessed on 30 June 2024).
  28. Mitzel, J.; Kabza, A.; Hunger, J.; Malkow, T.; Piela, P.; Laforet, H. Test Module P-08: Polarization Curve. Development of PEM Fuel Cell Stack Reference Test Procedures for Industry. STACK-TEST Project no. 303445. Fuel Cells and Hydrogen Joint Undertaking (FCH). 2015. Available online: http://stacktest.zsw-bw.de/fileadmin/stacktest/docs/Information_Material/Performance/TMs/TM_P-08_Polarisation_Curve.pdf (accessed on 30 June 2024).
  29. Mohsin, M.; Raza, R.; Mohsin-ul-Mulk, M.; Yousaf, A.; Hacker, V. Electrochemical characterization of polymer electrolyte membrane fuel cells and polarization curve analysis. Int. J. Hydrogen Energy 2020, 45, 24093–24107. [Google Scholar] [CrossRef]
  30. Piela, P.; Malkow, T.; de Marco, G. Test Module P-10c: H2-PEMFC and DMFC Stack Electrochemical Impedance Spectroscopy. Development of PEM Fuel Cell Stack Reference Test Procedures for Industry. STACK-TEST Project no. 303445. Fuel Cells and Hydrogen Joint Undertaking (FCH). 2015. Available online: http://stacktest.zsw-bw.de/fileadmin/stacktest/docs/Information_Material/Performance/TMs/TM_P-10c_Electrochemical_Method_Impedance_Spectroscopy.pdf (accessed on 30 June 2024).
  31. Macdonald, J.R.; Johnson, W.B. Fundamentals of Impedance Spectroscopy. In Theory, Experiment, and Aplications, 3rd ed.; Barsoukov, E., Macdonald, J.R., Eds.; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2018; Chapter 1; pp. 1–20. ISBN 9781119333173. [Google Scholar]
  32. Yuan, X.-Z.; Song, C.; Wang, H.; Zhang, J. Electrochemical Impedance Spectroscopy in PEM Fuel Cells Fundamentals and Applications; Springer: London, UK, 2010. [Google Scholar]
  33. Mérida, W.; Harrington, D.A.; le Canut, J.M.; McLean, G. Characterization of proton exchange membrane fuel cell (PEMFC) failures via electrochemical impedance spectroscopy. J. Power Sources 2006, 161, 264–274. [Google Scholar] [CrossRef]
  34. Mousa, G.; Golnaraghi, F.; DeVaal, J.; Young, A. Detecting proton exchange membrane fuel cell hydrogen leak using electrochemical impedance spectroscopy method. J. Power Sources 2014, 246, 110–116. [Google Scholar] [CrossRef]
  35. Tang, Z.; Huang, Q.A.; Wang, Y.J.; Zhang, F.; Li, W.; Li, A.; Zhang, L.; Zhang, J. Recent progress in the use of electrochemical impedance spectroscopy for the measurement, monitoring, diagnosis and optimization of proton exchange membrane fuel cell performance. J. Power Sources 2020, 468, 228361. [Google Scholar] [CrossRef]
  36. le Canut, J.M.; Abouatallah, R.M.; Harrington, D.A. Detection of Membrane Drying, Fuel Cell Flooding, and Anode Catalyst Poisoning on PEMFC Stacks by Electrochemical Impedance Spectroscopy. J. Electrochem. Soc. 2006, 153, A857–A864. [Google Scholar] [CrossRef]
  37. Guo, X.; Dong, Z.; Shen, J.; Xu, Y.; He, Q.; Zhao, X.; Ding, Z. Towards intelligent and integrated architecture for hydrogen fuel cell system: Challenges and approaches. Natl. Sci. Open 2023, 2, 20220038. [Google Scholar] [CrossRef]
  38. Kong, Z.; Wu, J.; Liu, Z.; Yan, D.; Wu, Z.P.; Zhong, C.J. Advanced electrocatalysts for fuel cells: Evolution of active sites and synergistic properties of catalysts and carrier materials. Exploration 2024, 4, 20230052. [Google Scholar] [CrossRef]
  39. He, Y.; Tan, Q.; Lu, L.; Sokolowski, J.; Wu, G. Metal-Nitrogen-Carbon Catalysts for Oxygen Reduction in PEM Fuel Cells: Self-Template Synthesis Approach to Enhancing Catalytic Activity and Stability. Electrochem. Energy Rev. 2019, 2, 231–251. [Google Scholar] [CrossRef]
Figure 1. Equivalent circuit used to fit the EIS curves (CPE = constant phase element).
Figure 1. Equivalent circuit used to fit the EIS curves (CPE = constant phase element).
Fuels 05 00026 g001
Figure 3. Left picture, Cell 79 polarization curve. Right picture, EIS curves.
Figure 3. Left picture, Cell 79 polarization curve. Right picture, EIS curves.
Fuels 05 00026 g003
Figure 4. Individual representation of the EIS curves, Cell 79.
Figure 4. Individual representation of the EIS curves, Cell 79.
Fuels 05 00026 g004
Figure 5. Individual representation of the filtered EIS curves, Cell 79.
Figure 5. Individual representation of the filtered EIS curves, Cell 79.
Fuels 05 00026 g005
Figure 6. R t c versus current and its fitting curve for Cell 79.
Figure 6. R t c versus current and its fitting curve for Cell 79.
Fuels 05 00026 g006
Figure 7. Cell 79 polarization curve and its fitting curve.
Figure 7. Cell 79 polarization curve and its fitting curve.
Fuels 05 00026 g007
Figure 8. Left picture, Cell 40 polarization curve. Right picture, EIS curves.
Figure 8. Left picture, Cell 40 polarization curve. Right picture, EIS curves.
Fuels 05 00026 g008
Figure 9. Individual representation of the EIS curves, Cell 40.
Figure 9. Individual representation of the EIS curves, Cell 40.
Fuels 05 00026 g009
Figure 10. Cell 40 polarization curve and its fitting curve.
Figure 10. Cell 40 polarization curve and its fitting curve.
Fuels 05 00026 g010
Figure 11. Left picture, Cell 2 polarization curve. Right picture, EIS curves.
Figure 11. Left picture, Cell 2 polarization curve. Right picture, EIS curves.
Fuels 05 00026 g011
Figure 12. Individual representation of the EIS curves, Cell 2.
Figure 12. Individual representation of the EIS curves, Cell 2.
Fuels 05 00026 g012
Figure 13. Cell 2 polarization curve and its fitting curve.
Figure 13. Cell 2 polarization curve and its fitting curve.
Fuels 05 00026 g013
Figure 14. Polarization curve comparison. Cell 2, Cell 40, and Cell 79.
Figure 14. Polarization curve comparison. Cell 2, Cell 40, and Cell 79.
Fuels 05 00026 g014
Figure 15. Measured and fitted open circuit voltage comparison, Cell 2, Cell 40, and Cell 79. EIS test.
Figure 15. Measured and fitted open circuit voltage comparison, Cell 2, Cell 40, and Cell 79. EIS test.
Fuels 05 00026 g015
Figure 16. Estimated and fitted internal ohmic resistance comparison, Cell 2, Cell 40, and Cell 79. EIS test.
Figure 16. Estimated and fitted internal ohmic resistance comparison, Cell 2, Cell 40, and Cell 79. EIS test.
Fuels 05 00026 g016
Figure 17. Representation of the Rtc (Ω) and AVL (Ω) values of Cell 2, Cell 40, and Cell 79 for the [0–5] current range.
Figure 17. Representation of the Rtc (Ω) and AVL (Ω) values of Cell 2, Cell 40, and Cell 79 for the [0–5] current range.
Fuels 05 00026 g017
Figure 18. Representation of CPE and R t c behaviors at different currents, according to [18].
Figure 18. Representation of CPE and R t c behaviors at different currents, according to [18].
Fuels 05 00026 g018
Table 1. Fitting parameter values for the equivalent circuit of Figure 1, Cell 79. T (F) and P are typical parameters of a CPE impedance.
Table 1. Fitting parameter values for the equivalent circuit of Figure 1, Cell 79. T (F) and P are typical parameters of a CPE impedance.
I (A)R (Ω)Rtc (Ω)T (F)P
0.00 ± 0.010.106 ± 0.0010.854 ± 0.0103.134 ± 0.0320.854 ± 0.010
0.50 ± 0.010.107 ± 0.0000.102 ± 0.0022.611 ± 0.0630.816 ± 0.010
1.00 ± 0.010.107 ± 0.000(5.309 ± 0.070)·10−22.279 ± 0.0680.833 ± 0.010
2.00 ± 0.010.107 ± 0.000(2.907 ± 0.036)·10−21.971 ± 0.0750.854 ± 0.011
5.00 ± 0.010.107 ± 0.000(1.414 ± 0.002)·10−21.600 ± 0.0780.898 ± 0.013
10.00 ± 0.010.106 ± 0.000(8.679 ± 0.201)·10−31.396 ± 0.1350.945 ± 0.023
20.00 ± 0.010.106 ± 0.000(6.478 ± 0.212)·10−31.287 ± 0.1860.955 ± 0.032
30.00 ± 0.010.106 ± 0.000(5.915 ± 0.312)·10−31.961 ± 0.4380.911 ± 0.052
40.00 ± 0.010.106 ± 0.000(6.721 ± 0.352)·10−33.028 ± 0.6090.839 ± 0.050
50.00 ± 0.010.106 ± 0.000(7.401 ± 0.424)·10−33.117 ± 0.7000.863 ± 0.056
Table 2. Fitting parameter values for the equivalent circuit of Figure 1, Cell 40. T (F) and P are typical parameters of a CPE impedance.
Table 2. Fitting parameter values for the equivalent circuit of Figure 1, Cell 40. T (F) and P are typical parameters of a CPE impedance.
I (A)R (Ω)Rtc (Ω)T (F)P
0.00 ± 0.01(9.37 ± 0.06)·10−21.163 ± 0.0182.973 ± 0.0360.844 ± 0.007
0.50 ± 0.01(9.44 ± 0.03)·10−20.104 ± 0.0022.336 ± 0.0680.834 ± 0.012
1.00 ± 0.01(9.35 ± 0.02)·10−2(5.971 ± 0.119)·10−22.327 ± 0.0860.815 ± 0.012
2.00 ± 0.01(9.20 ± 0.02)·10−2(3.278 ± 0.044)·10−21.875 ± 0.0760.856 ± 0.012
5.00 ± 0.01(9.00 ± 0.01)·10−2(1.586 ± 0.020)·10−21.521 ± 0.0730.891 ± 0.012
10.00 ± 0.01(8.91 ± 0.02)·10−2(9.480 ± 0.218)·10−31.192 ± 0.1090.948 ± 0.022
20.00 ± 0.01(8.81 ± 0.01)·10−2(6.478 ± 0.212)·10−31.287 ± 0.1860.955 ± 0.032
30.00 ± 0.01(8.80 ± 0.01)·10−2(5.738 ± 0.233)·10−31.512 ± 0.2700.942 ± 0.040
40.00 ± 0.01(8.76 ± 0.02)·10−2(5.950 ± 0.321)·10−31.760 ± 0.4200.928 ± 0.054
50.00 ± 0.01(8.73 ± 0.03)·10−2(6.310 ± 0.363)·10−32.408 ± 0.5500.861 ± 0.054
Table 3. Fitting parameter values for the equivalent circuit of Figure 1, Cell 2. T (F) and P are typical parameters of a CPE impedance.
Table 3. Fitting parameter values for the equivalent circuit of Figure 1, Cell 2. T (F) and P are typical parameters of a CPE impedance.
I (A)R (Ω)Rtc (Ω)T (F)P
0.00 ± 0.01(9.37 ± 0.06)·10−21.163 ± 0.0182.973 ± 0.0360.844 ± 0.007
0.50 ± 0.01(9.90 ± 0.03)·10−20.107 ± 0.0022.337 ± 0.0740.829 ± 0.013
1.00 ± 0.01(9.80 ± 0.03)·10−2(5.969 ± 0.103)·10−22.231 ± 0.0800.826 ± 0.012
2.00 ± 0.01(9.65 ± 0.02)·10−2(3.236 ± 0.043)·10−21.895 ± 0.0780.850 ± 0.012
5.00 ± 0.01(9.43 ± 0.02)·10−2(1.514 ± 0.022)·10−21.460 ± 0.0810.904 ± 0.015
10.00 ± 0.01(9.36 ± 0.01)·10−2(9.246 ± 0.181)·10−31.225 ± 0.1010.950 ± 0.020
20.00 ± 0.01(9.29 ± 0.01)·10−2(6.447 ± 0.186)·10−31.261 ± 0.1600.960 ± 0.028
30.00 ± 0.01(9.27 ± 0.01)·10−2(5.700 ± 0.265)·10−31.559 ± 0.3250.944 ± 0.046
40.00 ± 0.01(9.21 ± 0.02)·10−2(5.431 ± 0.299)·10−31.561 ± 0.3890.963 ± 0.056
50.00 ± 0.01(9.18 ± 0.05)·10−2(7.266 ± 0.732)·10−34.209 ± 1.50310.779 ± 0.091
Table 4. Rtc (Ω) and AVL (Ω) values of Cell 2, Cell 40, and Cell 79 for the [0–5] current range.
Table 4. Rtc (Ω) and AVL (Ω) values of Cell 2, Cell 40, and Cell 79 for the [0–5] current range.
Rtc (Ω)AVL (Ω)
Cell 2Cell 40Cell 79 Cell 2Cell 40Cell 79
0.0 A1.16301.16300.85400.0 A2.62022.60531.8586
0.5 A0.10700.10400.10200.5 A0.23870.23110.2176
1.0 A0.05970.05970.05311.0 A0.13700.13520.1227
2.0 A0.03240.03280.02912.0 A0.07720.07780.0674
5.0 A0.01510.01590.01415.0 A0.03560.03690.0299
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

Gómez, G.; Argumosa, P.; Maellas, J. A Real Case: How to Combine Polarization Curve and EIS Techniques to Identify Problematic Cells in a Commercial PEM Stack. Fuels 2024, 5, 476-493. https://doi.org/10.3390/fuels5030026

AMA Style

Gómez G, Argumosa P, Maellas J. A Real Case: How to Combine Polarization Curve and EIS Techniques to Identify Problematic Cells in a Commercial PEM Stack. Fuels. 2024; 5(3):476-493. https://doi.org/10.3390/fuels5030026

Chicago/Turabian Style

Gómez, Guillermo, Pilar Argumosa, and Jesús Maellas. 2024. "A Real Case: How to Combine Polarization Curve and EIS Techniques to Identify Problematic Cells in a Commercial PEM Stack" Fuels 5, no. 3: 476-493. https://doi.org/10.3390/fuels5030026

APA Style

Gómez, G., Argumosa, P., & Maellas, J. (2024). A Real Case: How to Combine Polarization Curve and EIS Techniques to Identify Problematic Cells in a Commercial PEM Stack. Fuels, 5(3), 476-493. https://doi.org/10.3390/fuels5030026

Article Metrics

Back to TopTop