Next Article in Journal
Research on Optimal Control of Excavator Negative Control Swing System
Next Article in Special Issue
Optimization of an Inter-Plant Hydrogen Network: A Simultaneous Approach to Solving Multi-Period Optimization Problems
Previous Article in Journal
Effect of Mixed Commercial Cold Flow Improvers on Flow Properties of Biodiesel from Waste Cooking Oil
Previous Article in Special Issue
Short-Term Wind Power Prediction Using GA-BP Neural Network Based on DBSCAN Algorithm Outlier Identification
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Development of an Oxygen Pressure Estimator Using the Immersion and Invariance Method for a Particular PEMFC System

by
Ángel Hernández-Gómez
1,†,
Victor Ramirez
1,2,*,† and
Belem Saldivar
2,3
1
Department of Renewable Energy, Centro de Investigación Científica de Yucatán (CICY), Yucatán P.C. 97205, Mexico
2
Cátedras CONACYT, Ciudad de México P.C. 03940, Mexico
3
Facultad de Ingeniería, Universidad Autónoma del Estado de México (UAEM), Toluca P.C. 50000, Mexico
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Processes 2020, 8(9), 1095; https://doi.org/10.3390/pr8091095
Submission received: 28 July 2020 / Revised: 25 August 2020 / Accepted: 27 August 2020 / Published: 3 September 2020
(This article belongs to the Special Issue Multi-Period Optimization of Sustainable Energy Systems)

Abstract

:
The fault detection method has been used usually to give a diagnosis of the performance and efficiency in the proton exchange membrane fuel cell (PEMFC) systems. To be able to use this method a lot of sensors are implemented in the PEMFC to measure different parameters like pressure, temperature, voltage, and electrical current. However, despite the high reliability of the sensors, they can fail or give erroneous measurements. To address this problem, an efficient solution to replace the sensors must be found. For this reason, in this work, the immersion and invariance method is proposed to develop an oxygen pressure estimator based on the voltage, electrical current density, and temperature measurements. The estimator stability region is calculated by applying Lyapunov’s Theorem and constraints to achieve stability are established for the oxygen pressure, electrical current density, and temperature. Under these estimator requirements, oxygen pressure measurements of high reliability are obtained to fault diagnosis without the need to use an oxygen sensor.

1. Introduction

Fuel cell (FC) system is an advanced power system necessary for a clean, sustainable, and environmentally friendly future, because FCs are promising candidates as an alternative to conventional fossil fuels, due to their higher energy density, energy efficiency, and very low emissions [1,2,3]. The main operation of the FCs is to transform gaseous fuel chemical energy into electricity. Besides, the FCs can be used as alternative stationary and mobile power source [4,5]. The main types of FCs are proton exchange membrane, direct methanol, solid oxide, molten carbonate, phosphoric acid, alkaline, and microbial [6].
In particular, the proton exchange membrane fuel cell (PEMFC) has attracted the attention of researchers in the last few decades due to its characteristics as low operating temperature, low noise, quick start-up capability, light mass, and high-power density [2,4,6,7]. The PEMFCs have recently passed the test phase and have slightly reached the commercialization stage due to the impressive research effort [8]. However, the two biggest limitations preventing the PEMFC system from further commercialization are its reliability and durability [7].
A lot of studies on PEMFC performance have been carried out, since three-dimensional simulation models to more detailed measurement techniques, such as electrochemical impedance spectroscopy [9,10]. To have a PEMFC diagnosis, the fault detection method has been used commonly to guarantee correct and safe operation in the PEMFC system [7,11,12]. However, to achieve such a diagnosis, several sensors have been used to measure different parameters like the mass flow, oxygen pressure, hydrogen pressure, compressor velocity, electrical current, water pressure, voltage, and temperature of the stack [11,13,14].
A lot of researchers have worked on the development of sensors with high reliability [15,16,17,18]. These devices must present characteristics, such as high sensitivity and selectivity, robustness, fast response time, operation at high temperature and low power consumption [19,20,21]. However, in real applications, the reliability of sensors during the system operation is variable. Thus, inaccurate sensor measurements can provide misleading results in PEMFC fault diagnosis, which can end in failures and damages of the PEMFC system [7,11]. To solve this problem, novel methods have been proposed to reduce errors in PEMFC fault diagnosis [9,10]. For this reason, an efficient method to replace the oxygen sensor is proposed in this work, since the oxygen management system is an important subsystem, which is used for supplying proper oxygen pressure in the PEMFC stack cathode. Besides, the complexity and nonlinearity of the oxygen pressure are difficult to model [22]. So, using the voltage, electrical current density, and temperature measurements and applying the immersion and invariance (gradient estimator) method it is possible to develop an oxygen pressure estimator for getting high-reliability oxygen measurements avoiding the use of oxygen sensor for PEMFC system fault diagnosis.
The paper is organized as follows, the formulation of a gradient estimator to develop the oxygen pressure estimator is described in Section 2. The PEMFC potential-current behavior is discussed in Section 3. The oxygen pressure estimator applied to a PEMFC system is presented in Section 4. The simulation and results are introduced in Section 5. Finally, some concluding remarks are presented in Section 6.

2. Formulation of Gradient Estimator

The immersion and invariance (gradient estimator) method has been proposed to solve problems of stabilization and adaptive control of nonlinear systems, which are present in any real practical problem [23,24,25,26]. The key step for the estimator development using this method is the construction of a monotone mapping, which explicitly depends on some of the estimator tuning parameters [27,28]. For these reasons, in this work, this method has been used to develop the oxygen pressure estimator.
The estimator design is formulated by proposing a function where the system behavior representation distinguishes between measurable and not measurable signals. As shown in [28,29], there is a general kind of function dependent on two variables θ and ξ expressed by
F ( θ , ξ ) = G ( θ ) + H ( ξ ) + K ( θ , ξ )
with θ > 0 and ξ > 0 , where ξ and θ are known and time-dependent variables, such that measurable signals F ( θ , ξ ) and H ( ξ ) are represented by
y ( t ) = F ( θ , ξ ) H ( ξ ) .
Indeed, the representation in the non-linear regression form will be
y ( t ) = ϕ ( θ , ξ ) ,
where
ϕ ( θ , ξ ) : = G ( θ ) + K ( θ , ξ ) .
Given this formulation, the following proposition can be stated.
Proposition 1.
Consider the function ϕ ( θ , ξ ) , where F ( θ , ξ ) and H ( ξ ) are known and the variable corresponding to the non-linear regression model satisfies that the partial derivative of ϕ ( θ , ξ ) with respect to θ is greater than zero. Then, the gradient estimator is given by
θ ^ ˙ = γ ( y ( t ) ϕ ( θ ^ , ξ ) )
with γ > 0 ensuring that
lim t θ ^ = θ ,
for all initial condition θ ^ 0 such as ϕ ( θ ^ 0 , ξ ) θ > 0 .
Proof. 
To show that the immersion and invariance estimator converges to the desired value, it is necessary to use the monotonicity property of the function ϕ ( θ , ξ ) concerning θ . Then, as:
ϕ ( θ , ξ ) θ > 0 ,
the function is strictly monotonically increasing and also fulfills
( θ ^ θ ) ϕ ( θ ^ , ξ ) ϕ ( θ , ξ ) > 0 θ ^ θ ,
taking the Lyapunov’s function candidate
V ( θ ^ ) = 1 2 γ ( θ ^ θ ) 2 ,
its time-derivative along the trajectories of (2)–(5) is given by
V ˙ = ( θ ^ θ ) [ ϕ ( θ ^ , ξ ) ϕ ( θ , ξ ) ] < 0 θ ^ θ .
Note that the negative definiteness of V ˙ immediately follows from (8). Then, the proof is completed by using Lyapunov’s Second Stability Theorem. □

3. PEMFC Potential-Current Behavior

An accurate mathematical model to represent the PEMFC potential V c has been reported in [30], where V c is a depending function of stack current, cathode pressure, reactant partial pressures, PEMFC temperature, and membrane humidity using a combination of physical and empirical relationships, and can be expressed in terms of the Nernst’s potential E t h and the three main types of potential drops; activation V a c t , ohmic V o h m , and concentration V c o n .
V c ( θ , ξ ) = E t h ( θ ) V o h m ( ξ ) V a c t ( θ , ξ ) V c o n ( θ , ξ ) ,
where θ denotes the oxygen pressure ( a t m ), and ξ the electrical current density in the cell ( A · cm 2 ). Nernst’s potential E t h . The Nernst’s potential or open-circuit potential is the maximum power obtained by one cell corresponding to exchange Gibbs free energy as a result of the difference between reactant products and Gibbs’s free energy. It can be described by the following equation [30,31,32].
E t h ( θ ) = E 0 + B 1 ( T 0 T ) + B 2 T ln P H 2 θ 1 / 2 P H 2 O ,
where T 0 and T are the initial temperature and the cell temperature, respectively (K), P H 2 is a positive constant that represents the hydrogen pressure ( a t m ), and E 0 is the reference potential (V). B 1 and B 2 are positive constants that depend on stack temperature and potential ( V / K ) [30]. Water pressure is represented by P H 2 O ( a t m ).
Ohmic potential drop V o h m . The ohmic potential drop arises from the resistance of the polymer membrane to the transfer of protons and from the resistance of the electrode and the collector plate to the transfer of electrons [30,31,32].
V o h m ( ξ ) = ξ A f c R o h m ,
where R o h m > 0 is the internal electrical resistance ( Ω ) and A f c is the cell active area. Besides, the ohmic resistance can be expressed as a function of the membrane conductivity (cm 1 · Ω 1 ), σ m .
R o h m = t m σ m ,
where t m is the thickness of the membrane ( c m ), and σ m is a function of membrane water content λ m and the cell temperature T.
σ m = b 1 exp b 2 1 303 1 T ,
where b 1 is a function of membrane water content and b 2 is a constant [30].
b 1 = b 11 λ m b 12 ,
where b 2 , b 11 , and b 12 are usually determined empirically. In this work, the values for b 2 , b 11 , and b 12 are taken from [33].
Activation potential drop V a c t . The activation potential drop comes when the movement of electrons needs to break and form chemical bonds in the anode and cathode (i.e., part of the available energy is lost in driving the chemical reaction that transfers the electrons to and from the electrodes). Although the activation overvoltage occurs at both PEMFC electrodes, the reaction of hydrogen oxidation at the anode is faster than the reaction of oxygen [30,31,32].
V a c t ( θ ) = V o + V a ( θ ) ( 1 exp [ c 1 ξ ] ) ,
where c 1 is a constant. The functions V o and V a are both dependent on oxygen pressure and temperature. They have been calculated empirically by
V o = V o 0 + B 1 ( T 0 T ) 1.07551 B 2 T + 3 B 2 T 2 ln P c a P s a t P a t m ,
where V o 0 is the initial potential drop (V) at zero current density. P c a and P a t m are the pressures of the cathode and atmospheric, respectively ( a t m ). The water saturation pressure P s a t ( m P a ) is expressed as
log 10 P s a t = 1.69 10 10 T 4 + 3.85 10 7 T 3 3.39 10 4 T 2 + 0.143 T 20.92 .
The function V a is given as:
V a ( θ ) = B 3 θ 0.1173 + P s a t 2 + B 4 θ 0.1173 + P s a t + B 5 ,
where the constants B 3 , B 4 , and B 5 are dependent on the stack temperature and the voltage ( V / K ) and usually are determined empirically [30].
Concentration of potential drop V c o n . The concentration of potential drop corresponds to the concentration gradients formed due to mass diffusions from the flow channels to the reaction sites (catalyst area). The factors underlying this potential drop are high current densities, slow transportation of reactants and products, and water film covering the catalyst surfaces to the anode and cathode [30,31,32].
V c o n ( θ , ξ ) = ξ c 2 ξ I m a x c 3 ,
where c 3 R + is a constant, I m a x is the maximum electrical current density in the cell and c 2 is an oxygen pressure function [30].
c 2 = B 6 θ 0.1173 + P s a t + B 7 i f   h ( θ ) < 0 , B 8 θ 0.1173 + P s a t + B 9 i f   h ( θ ) 0 ,
where
h ( θ ) = θ + 0.1173 P s a t 0.2346   a t m ,
and B 6 , B 7 , B 8 , and B 9 are constant values that depend on the stack temperature and are usually determined empirically.
Lemma 1.
The discontinuous function c 2 defined in (22) can be approximated by the continuous function C 2 given below.
C 2 = D 1 θ 0.1173 + P s a t + D 2 ,
where
D 1 = 1 2 B 6 + 1 2 B 8 + 1 2 ( B 6 B 8 ) tanh [ h ( θ ) ] ,
and
D 2 = 1 2 B 7 + 1 2 B 9 + 1 2 ( B 7 B 9 ) tanh [ h ( θ ) ] .
The parameters values of the PEMFC voltage model are taken from [30] (see Table 1).

4. Application of Oxygen Pressure Estimator to a PEMFC System

The oxygen pressure estimator presented in this section is derived from the results presented in Section 2 and Section 3. The measurable signal is defined by applying the Equations (2)–(11).
y ( t ) = V c ( θ , ξ ) + V o h m ( ξ ) ,
where
ϕ ( θ , ξ ) : = E t h ( θ ) V a c t ( θ , ξ ) V c o n ( θ , ξ ) .
Now a proposition related to the PEMFC system is presented.
Proposition 2.
Consider the function ϕ ( θ , ξ ) , with C 2 and θ are greater than zero, such that, inequality (26) is satisfied
B 2 T ( P c a P s a t 3 θ ) 2 V a θ θ ( P c a P s a t ) 0 .
Then, ξ can be expressed in terms of θ and T as follows:
ξ = D B 2 T 2 P c a P s a t 3 θ θ ( P c a P s a t ) V a θ 1 c 3 + 1 ,
where
D = I m a x c 3 c 3 C 2 c 3 1 C 2 θ 1
and ϕ ( θ , ξ ) θ > 0 .
Proof. 
The proof starts with the partial derivative of ϕ with respect to θ , which is given by
ϕ ( θ , ξ ) θ = B 2 T 2 P c a P s a t 3 θ θ ( P c a P s a t ) V a θ 1 exp [ c 1 ξ ] C 2 θ c 3 C 2 c 3 1 ξ c 3 + 1 I m a x c 3 .
Now, taking the set of values ( θ , ξ ) , that satisfy inequality (28),
0 B 2 T 2 P c a P s a t 3 θ θ ( P c a P s a t ) V a θ 1 exp [ c 1 ξ ] C 2 θ c 3 C 2 c 3 1 ξ c 3 + 1 I m a x c 3 .
Since 0 < 1 exp [ c 1 ξ ] < 1 for ξ and c 1 > 0 , then
0 B 2 T 2 P c a P s a t 3 θ θ ( P c a P s a t ) V a θ C 2 θ c 3 C 2 c 3 1 ξ c 3 + 1 I m a x c 3 < ϕ ( θ , ξ ) θ .
So, the admissible limit values ( θ , ξ ) that satisfy inequality (28) can be found when this is equal to zero.
0 = B 2 T 2 P c a P s a t 3 θ θ ( P c a P s a t ) V a θ C 2 θ c 3 C 2 c 3 1 ξ c 3 + 1 I m a x c 3 ,
setting
D = I m a x c 3 c 3 C 2 c 3 1 C 2 θ 1 ,
thus,
ξ = D B 2 T 2 P c a P s a t 3 θ θ ( P c a P s a t ) V a θ 1 c 3 + 1 .
As C 2 θ > 0 , then,
B 2 T 2 P c a P s a t 3 θ θ ( P c a P s a t ) V a θ 0 ,
since θ > 0 and P c a P s a t > 0 , then,
B 2 T ( P c a P s a t 3 θ ) 2 V a θ θ ( P c a P s a t ) 0 .
 □
Now the following proposition is introduced as a result of the combination of Proposition 1 and Proposition 2. This result shows the estimator and its stability using Lyapunov’s functions.
Proposition 3.
Consider the function ϕ ( θ , ξ ) , C 2 and θ are greater than zero, such that, inequality (26) is satisfied and with ξ expressed as:
ξ = D B 2 T 2 P c a P s a t 3 θ θ ( P c a P s a t ) V a θ 1 c 3 + 1 .
Then, the gradient estimator of oxygen pressure is given by
θ ^ ˙ = γ ( y ( t ) ϕ ( θ ^ , ξ ) ) ,
with γ > 0 , ensuring that
lim t θ ^ = θ .
Proof. 
For the values of θ and ξ stated in the hypothesis of Proposition 2, it is obtained that the partial derivative of ϕ ( θ , ξ ) with respect to θ is greater than zero. Then, by Proposition 1, the gradient estimator of oxygen pressure is given by
θ ^ ˙ = γ ( y ( t ) ϕ ( θ ^ , ξ ) ) ,
with γ > 0 , ensuring that
lim t θ ^ = θ .
 □

5. Simulations and Results

The Runge–Kutta fourth-order algorithm, described in [34], and the values of the parameters given in the Table 1 were used to perform the simulations. The first step was to determine the stability region for the estimator under the established constrains of the Propositions 2 and 3. The estimator stability region is given within the interval (0 atm, 0.45 atm) and the simulation results of such constraints are shown in Figure 1 and Figure 2. The behavior of the partial derivative of ϕ with respect to θ as a function of θ and ξ for different temperatures is shown in Figure 1.
The behavior of θ and ξ considering the established constraints for different temperatures is shown in Figure 2.
Within stability region, the oxygen pressure estimator and the PEMFC potential-current simulations were performed using oxygen pressure equal to 0.3 a t m and different values for θ 0 and γ . The oxygen pressure estimator shows an asymptotic convergence to the proposed value for oxygen pressure. The estimator behavior can be appreciated for different values θ 0 in Figure 3, and different values of γ in Figure 4.
The electrical current density calculated based on the estimator proved an asymptotic convergence to the electrical current density calculated for oxygen pressure equal to 0.3 a t m , the simulation is shown for different values of θ 0 in Figure 5, and for different values of γ in Figure 6.
The cell potential calculated based on the estimator evidenced an asymptotic convergence to the potential calculated for oxygen pressure equal to 0.3 a t m , the simulation results for different values of θ 0 are shown in Figure 7 and for different values of γ in Figure 8.
Finally, the power or potential-current performance curve based on the estimator demonstrated an asymptotic convergence to the power for oxygen pressure equal to 0.3 atm, the simulation is shown for different values of θ 0 in Figure 9, and for different values of γ in Figure 10. This curve has proved to be of vital importance for the PEMFC system fault diagnosis [35].

6. Conclusions

To avoid oxygen sensors for PEMFCs, an oxygen pressure estimator has been developed based on the immersion and invariance (gradient estimator) method, and its stability conditions are established using Lyapunov’s Theorem. Additionally, in this work, the PEMFC electrical current density has been characterized in terms of oxygen pressure and temperature under certain constraints.
The oxygen pressure estimator presents an absolute convergence within the stability region to the measurable value of oxygen pressure. However, the corresponding working condition can be different because it is directly related to laboratory environmental conditions. So, the next step is to evaluate the performance of the proposed estimator under different PEMFC conditions to improve the oxygen pressure estimator.

Author Contributions

Conceptualization, Á.H.-G. and V.R.; methodology, V.R.; software, Á.H.-G.; validation, Á.H.-G. and V.R.; formal analysis, Á.H.-G.; investigation, Á.H.-G. and V.R.; resources, V.R.; writing—original draft preparation, Á.H.-G. and V.R.; writing—review and editing, Á.H.-G., V.R., and B.S.; visualization, Á.H.-G.; supervision, V.R.; project administration, V.R.; funding acquisition, V.R. All authors have read and agreed to the published version of the manuscript.

Funding

This work has been supported by CONACyT-Mexico under the project 2015-01-786 (National Problems).

Acknowledgments

Thanks to the CONACYT-Mexico program Becas Nacional (Tradicional) 2018–2 and the scholarship 2018–000068–02NACF.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sopian, K.; Daud, W.R.W. Challenges and future developments in proton exchange membrane fuel cells. Renew. Energy 2006, 31, 719–727. [Google Scholar] [CrossRef]
  2. Shen, M.; Scott, K. Power loss and its effect on fuel cell performance. J. Power Sources 2005, 148, 24–31. [Google Scholar] [CrossRef]
  3. Leonardi, S.G.; Bonavita, A.; Donato, N.; Neri, G. Development of a hydrogen dual sensor for fuel cell applications. Int. J. Hydrogen Energy 2018, 43, 11896–11902. [Google Scholar] [CrossRef]
  4. Pukrushpan, J.T.; Stefanopoulou, A.G.; Peng, H. Chapter one-background and introduction. In Control Fuel Cell Power System; Advances in Industrial Control; Grimble, M.J., Johnson, M.A., Eds.; Springer: London, UK, 2004; pp. 1–13. [Google Scholar] [CrossRef]
  5. Larminie, J.; Dicks, A. Fuel Cell Systems Explained, 2nd ed.; John Wiley & Sons, Ltd.: Chichester, UK, 2013; Chapter 1; pp. 1–24. [Google Scholar] [CrossRef]
  6. Daud, W.; Rosli, R.; Majlan, E.; Hamid, S.; Mohamed, R.; Husaini, T. PEM fuel cell system control: A review. Renew. Energy 2017, 113, 620–638. [Google Scholar] [CrossRef]
  7. Mao, L.; Jackson, L.; Huang, W.; Li, Z.; Davies, B. Polymer electrolyte membrane fuel cell fault diagnosis and sensor abnormality identification using sensor selection method. J. Power Sources 2020, 447, 227394. [Google Scholar] [CrossRef]
  8. Wee, J.-H. Applications of proton exchange membrane fuel cell systems. Renew. Sustain. Energy Rev. 2007, 11, 1720–1738. [Google Scholar] [CrossRef]
  9. Zhang, X.; Zhou, J.; Chen, W. Data-driven fault diagnosis for PEMFC systems of hybrid tram based on deep learning. Int. J. Hydrogen Energy 2020, 45, 13483–13495. [Google Scholar] [CrossRef]
  10. Abbaspour, A.; Yen, K.K.; Forouzannezhad, P.; Sargolzaei, A. An Adaptive Resilient Control Approach for Pressure Control in Proton Exchange Membrane Fuel Cells. IEEE Trans. Ind. Appl. 2019, 55, 6344–6354. [Google Scholar] [CrossRef]
  11. Li, S.; Aitouche, A.; Wang, H.; Christov, N. Sensor fault estimation of PEM fuel cells using Takagi Sugeno fuzzy model. Int. J. Hydrogen Energy 2020, 45, 11267–11275. [Google Scholar] [CrossRef]
  12. 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]
  13. Higgins, S.R.; Ewan, J.; St-Pierre, J.; Severa, G.; Davies, K.; Bethune, K.; Goodarzi, A.; Rocheleau, R. Environmental sensor system for expanded capability of PEM fuel cell use in high air contaminant conditions. Int. J. Hydrogen Energy 2018, 43, 22584–22594. [Google Scholar] [CrossRef]
  14. Arama, F.Z.; Mammar, K.; Laribi, S.; Necaibia, A.; Ghaitaoui, T. Implementation of sensor based on neural networks technique to predict the PEM fuel cell hydration state. J. Energy Storage 2020, 27, 101051. [Google Scholar] [CrossRef]
  15. Jung, S.-W.; Lee, E.K.; Kim, J.H.; Lee, S.-Y. High-concentration nafion-based hydrogen sensor for fuel-cell electric vehicles. Solid State Ion. 2020, 344, 115134. [Google Scholar] [CrossRef]
  16. Xiao, N.; Wu, R.; Huang, J.J.; Selvaganapathy, P.R. Development of a xurographically fabricated miniaturized low-cost, high-performance microbial fuel cell and its application for sensing biological oxygen demand. Sens. Actuators Chem. 2019, 127432. [Google Scholar] [CrossRef]
  17. Lee, C.-Y.; Lin, J.-T.; Chen, C.-H.; Lee, S.-J.; Wang, Y.-S. Development of a four-in-one sensor for low temperature fuel cell. Renew. Energy 2019, 135, 1452–1465. [Google Scholar] [CrossRef]
  18. He, L.; Liu, Q.; Zhang, S.; Zhang, X.; Gong, C.; Shu, H.; Wang, G.; Liu, H.; Wen, S.; Zhang, B. High sensitivity of TiO2 nanorod array electrode for photoelectrochemical glucose sensor and its photo fuel cell application. Electrochem. Commun. 2018, 94, 18–22. [Google Scholar] [CrossRef]
  19. Montpart, N.; Baeza, M.; Baeza, J.A.; Guisasola, A. Low-cost fuel-cell based sensor of hydrogen production in lab scale microbial electrolysis cells. Int. J. Hydrogen Energy 2016, 41, 20465–20472. [Google Scholar] [CrossRef]
  20. Lavanya, N.; Sekar, C.; Fazio, E.; Neri, F.; Leonardi, S.; Neri, G. Development of a selective hydrogen leak sensor based on chemically doped SnO2 for automotive applications. Int. J. Hydrogen Energy 2017, 42, 10645–10655. [Google Scholar] [CrossRef]
  21. Hayakawa, I.; Iwamoto, Y.; Kikuta, K.; Hirano, S. Gas sensing properties of platinum dispersed-TiO2 thin film derived from precursor. Sens. Actuators Chem. 2000, 62, 55–60. [Google Scholar] [CrossRef]
  22. Yang, D.; Wang, Y.; Chen, Z. Robust fault diagnosis and fault tolerant control for PEMFC system based on an augmented LPV observer. Int. J. Hydrog. Energy 2020, 45, 13508–13522. [Google Scholar] [CrossRef]
  23. Astolfi, A.; Ortega, R. Immersion and invariance: A new tool for stabilization and adaptive control of nonlinear systems. IEEE Trans. Autom. Control. 2003, 48, 590–606. [Google Scholar] [CrossRef] [Green Version]
  24. Hu, J.; Zhang, H. Immersion and invariance based command-filtered adaptive backstepping control of VTOL vehicles. Automatica 2013, 49, 2160–2167. [Google Scholar] [CrossRef]
  25. Zhu, R.; Wang, H.; Yin, G.; Ding, Z. High performance nonlinear adaptive control of temperature in cryogenic wind tunnel. Int. J. Robust Nonlinear Control 2019, 25, 5118–5136. [Google Scholar] [CrossRef]
  26. Ortega, R.; Nikiforov, V.; Gerasimov, D. On modified parameter estimators for identification and adaptive control, A unified framework and some new schemes. Annu. Rev. Control. 2020. [Google Scholar] [CrossRef]
  27. Liu, X.; Ortega, R.; Su, H.; Chu, J. Immersion and invariance adaptive control of nonlinearly parameterized nonlinear systems. IEEE Trans. Autom. Control 2010, 55, 2209–2214. [Google Scholar] [CrossRef]
  28. Ortega, R.; Liu, X.; Su, H.; Chu, J. Immersion and invariance adaptive control of nonlinearly parameterized nonlinear systems *. IFAC Proc. Vol. 2010, 43, 641–646. [Google Scholar] [CrossRef]
  29. Liu, X.; Ortega, R.; Su, H.; Chu, J. On adaptive control of nonlinearly parameterized nonlinear systems: Towards a constructive procedure. Syst. Control Lett. 2011, 60, 36–43. [Google Scholar] [CrossRef]
  30. Pukrushpan, J.T.; Stefanopoulou, A.G.; Peng, H. Chapter three-fuel cell system model: Fuel cell stack. In Control Fuel Cell Power System; Advances in Industrial Control; Grimble, M.J., Johnson, M.A., Eds.; Springer: London, UK, 2004; pp. 31–56. [Google Scholar] [CrossRef]
  31. Musio, F.; Tacchi, F.; Omati, L.; Stampino, P.G.; Dotelli, G.; Limonta, S.; Brivio, D.; Grassini, P. PEMFC system simulation in matlab-simulink® environment. Int. J. Hydrogen Energy 2011, 36, 8045–8052. [Google Scholar] [CrossRef]
  32. Sankar, K.; Aguan, K.; Jana, A.K. A proton exchange membrane fuel cell with an airflow cooling system: Dynamics, validation and nonlinear control. Energy Convers. Manag. 2019, 183, 230–240. [Google Scholar] [CrossRef]
  33. Springer, T.E.; Zawodzinski, T.A.; Gottesfeld, S. Polymer electrolyte fuel cell model. J. Electrochem. Soc. 1991, 138, 2334–2342. [Google Scholar] [CrossRef]
  34. Burden, R.L.; Faires, J.D.; Burden, A.M. Chapter Five-Problemas de Valor Inicial Para Ecuaciones de Diferenciales Ordinarias. Available online: https://latinoamerica.cengage.com/ls/analisis-numerico-2/ (accessed on 27 August 2020).
  35. Shen, M.; Meuleman, W.; Scott, K. The characteristics of power generation of static state fuel cells. J. Power Sources 2003, 115, 203–209. [Google Scholar] [CrossRef]
Figure 1. Behavior of the partial derivative of ϕ concerning θ .
Figure 1. Behavior of the partial derivative of ϕ concerning θ .
Processes 08 01095 g001
Figure 2. Behavior of the partial derivative of ϕ concerning θ .
Figure 2. Behavior of the partial derivative of ϕ concerning θ .
Processes 08 01095 g002
Figure 3. Estimator behavior with different values θ 0 .
Figure 3. Estimator behavior with different values θ 0 .
Processes 08 01095 g003
Figure 4. Estimator behavior with different values γ .
Figure 4. Estimator behavior with different values γ .
Processes 08 01095 g004
Figure 5. Simulation of electrical current density stability with different values θ 0 .
Figure 5. Simulation of electrical current density stability with different values θ 0 .
Processes 08 01095 g005
Figure 6. Simulation of electrical current density stability with different values γ .
Figure 6. Simulation of electrical current density stability with different values γ .
Processes 08 01095 g006
Figure 7. Simulation of potential stability with different values θ 0 .
Figure 7. Simulation of potential stability with different values θ 0 .
Processes 08 01095 g007
Figure 8. Simulation of potential stability with different values γ .
Figure 8. Simulation of potential stability with different values γ .
Processes 08 01095 g008
Figure 9. Simulation of power stability with different values θ 0 .
Figure 9. Simulation of power stability with different values θ 0 .
Processes 08 01095 g009
Figure 10. Simulation of power stability with different values γ .
Figure 10. Simulation of power stability with different values γ .
Processes 08 01095 g010
Table 1. Parameters for the PEMFC voltage model [30].
Table 1. Parameters for the PEMFC voltage model [30].
SymbolParameterValue
A f c Cell Active Area100 cm 2
E 0 Reference Potential1229 × 10 3 V
I m a x Maximum Current Density2.2 A·cm 2
P a t m Atmospheric Pressure1 atm
P c a Cathode Pressure2 atm
P H 2 Hydrogen Pressure1 atm
P H 2 O Water Pressure 396 × 10 3 atm
t m Membrane Thickness 125 × 10 4 cm
T 0 Initial Temperature 298.15 K
V o 0 Initial Potential Drop 279 × 10 3 V
b 1 Membrane Humidity Function 686 × 10 4 cm 1 · Ω 1
b 2 Constant350 K
b 11 Constant 513 × 10 5 SO 3 · H 2 O 1 · cm 1 · Ω 1
b 12 Constant 326 × 10 5 cm 1 · Ω 1
c 1 Constant10 A 1
c 3 Constant2 (-)
B 1 Temperature and Potential Function 85 × 10 5 V·K 1
B 2 Temperature and Potential Function 43085 × 10 6 V· K 1
B 3 Temperature and Potential Function 1618 × 10 8 T V· K 1 + 1618 × 10 5 V
B 4 Temperature and Potential Function 18 × 10 4 T V·K 1 166 × 10 3 V
B 5 Temperature and Potential Function 58 × 10 5 T V· K 1 + 5736 × 10 4 V
B 6 Temperature Function 358 × 10 5 T K 1 622 × 10 3
B 7 Temperature Function 725 × 10 6 T K 1 + 1.68
B 8 Temperature Function 433 × 10 6 T K 1 68 × 10 3
B 9 Temperature Function 8 × 10 4 T K 1 + 54 × 10 2
λ m Membrane Humidity14 H 2 O· (SO 3 ) 1

Share and Cite

MDPI and ACS Style

Hernández-Gómez, Á.; Ramirez, V.; Saldivar, B. Development of an Oxygen Pressure Estimator Using the Immersion and Invariance Method for a Particular PEMFC System. Processes 2020, 8, 1095. https://doi.org/10.3390/pr8091095

AMA Style

Hernández-Gómez Á, Ramirez V, Saldivar B. Development of an Oxygen Pressure Estimator Using the Immersion and Invariance Method for a Particular PEMFC System. Processes. 2020; 8(9):1095. https://doi.org/10.3390/pr8091095

Chicago/Turabian Style

Hernández-Gómez, Ángel, Victor Ramirez, and Belem Saldivar. 2020. "Development of an Oxygen Pressure Estimator Using the Immersion and Invariance Method for a Particular PEMFC System" Processes 8, no. 9: 1095. https://doi.org/10.3390/pr8091095

APA Style

Hernández-Gómez, Á., Ramirez, V., & Saldivar, B. (2020). Development of an Oxygen Pressure Estimator Using the Immersion and Invariance Method for a Particular PEMFC System. Processes, 8(9), 1095. https://doi.org/10.3390/pr8091095

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