14èmesJournées Internationales de Thermique JITH2009

27-29 Mars, 2009, Djerba, Tunisie

L'effet de température et de pression sur la performance

de la pile à combustible PEMFC

[Link] a,*, Hocine .BENMOUSSA a , Noureddine BOURMADAb , [Link] b

a
Laboratoire d’étude des systèmes énergétiques industriels (LESEI)
b
Laboratoire de chimie et de chimie de l'environnement. Université de Batna, Algérie

djamel_hd2@[Link]

Résumé:
Un modèle bidimensionnel transitoire des transferts de chaleur et de masse de
l’oxygène humide est présenté dans la partie cathodique d’une pile à combustible de type
membrane échangeuse de proton PEMFC. Le modèle se comporte des différentes équations de
conservation tel que : l’équation de masse, l’équation de mouvement, l’équation d'énergie. Un
programme basé sur la méthode des volumes finis a été réalisé pour simuler le système
d’équations. L’objectif est de combiner ce modèle avec un autre modèle électrochimique, pour
déterminer le potentiel réel de la pile à partir de la courbe de polarisation, qui est fonction de la
fraction de l’oxygène, la température et la pression. Les résultas montrent que le potentiel réel
d’une pile à combustible décroît par rapport au potentiel de Nernst, à cause des polarisations
qui sont: Polarisation d’activation, concentration et ohmique. On a conclue finalement que
l’augmentation de la performance (la puissance) de la pile est influencée par les paramètres
opératoires, quand on augmenter la pression et la température du gaz.
Mots clés :
Transfert de masse, diffusion, convection, migration, volume finis, PEMFC

1- INTRODUCTION

Les piles à combustible à membrane échangeuse de protons (PEMFC) concentrent

actuellement les plus gros efforts de recherche dans le monde. Elles fonctionnent à des
températures relativement basses entre 80°C et 100°C et à des pressions de 1 à 5 bars. Les
électrodes sont exposées à un flux de gaz, soit l’hydrogène (ou autre combustible) et l’oxygène
(ou air ou autre comburant), elles sont généralement poreuses permettant la diffusion des gaz
pour approvisionner en réactifs les zones actives où le catalyseur est en contact avec le
conducteur ionique. Les couches de diffusion entourant les électrodes diffusent les gaz,
permettent le transfert des électrons et assurent la gestion de l’eau pour une humidification
optimale de la membrane et éventuellement une évacuation. Les plaques bipolaires permettent
la distribution des gaz, l’évacuation de l’eau, la collecte du courant et assurent la connexion
entre les piles élémentaires qui composent la pile à combustible. Elles sont constituées de
matériaux possédant un niveau de conductivité électrique suffisant et une imperméabilité au
gaz supérieure à la membrane et doivent contenir des canaux de distribution de gaz
généralement obtenus par usinage ou directement par moulage du polymère chargé sur une
empreinte sans phase ultérieure d’usinage.
Cependant selon les applications visées, de nombreux problèmes de PEMFC restent à
résoudre avant leur commercialisation. Ayant une compréhension des éléments qui affectent la
performance de ce type de pile et de la façon dont les changements de conception affectent ces
éléments, il est important de concevoir de meilleures piles à combustible. Un modèle de
simulation est nécessaire qui peut être employé pour examiner des conceptions optimales.
14èmesJournées Internationales de Thermique JITH2009
27-29 Mars, 2009, Djerba, Tunisie

Actuellement, beaucoup de modèles existent et qui s’intéressent à cet effet. Mais, pour créer un
outil qui est utile aux concepteurs, le modèle de simulation doit tenir compte de beaucoup de
phénomènes physiques comme les transferts de chaleur et de masse, aussi bien les effets
électriques et électrochimiques qui se produisent et doit tenir compte des composants
importants de la pile à combustible, comme les plaques collectrices et la membrane. Le champ
d’écoulement de réactifs dans les électrodes, les couches de catalyseur, et la membrane doivent
être modelés ensemble de sorte que les effets des changements de conception soient reflétés à
travers le modèle. Un modèle complet pourrait encore identifier les distributions uniformes de
densité de
courant et pourrait aider à identifier les propriétés matérielles qui sont importantes pour la
performance de la pile à combustible.
Pour modéliser les piles à combustible de type PEMFC, il faut satisfaire les critères suivants du
modèle qui doit être:
- basé sur une description mathématique de la physique de la pile, puis résolu en une,
deux, et trois dimensions,
- utilisé comme outil d’étude parce qu’il facilite la compréhension des phénomènes
physiques qui régissent le comportement de la pile,
- utile comme outil de conception parce qu’il tient compte de l’évaluation de la
conception de la pile à combustible.
Ceci exige que le modèle soit flexible en ce qui concerne la géométrie et le matériel
approprié. Les résultats de la modélisation mathématique de la pile à combustible à membrane
échangeuse de protons (PEMFC) dépendent d’une série d’hypothèses. Ces dernières sont faites
de telles façons à simplifier la résolution numérique et non pas pour l’exactitude de l’étude du
modèle. Cependant, dans les modèles complets où on cherche à modéliser la pile entière, les
hypothèses qui négligent le transport de l’eau liquide et de traiter les couches de catalyseur
comme interfaces plutôt que comme volumes, diminuent considérablement l’exactitude et
l’utilité du modèle.
Le modèle complet devrait inclure le transport le long de la direction d'écoulement afin
d’étudier des effets importants comme l’épuisement de réactifs et la non uniformité de la
densité de courant et la température de pile. L’utilisation des propriétés physiques correctes est
aussi importante pour une nécessité mathématique, pour un modèle donné, plusieurs propriétés
physiques sont exigées qui peuvent être déterminées expérimentalement. Les propriétés qui ne
peuvent pas être déterminées expérimentalement sont typiquement traitées comme propriétés
estimées, certaines sont ajustées pour avoir les résultats modélisés en accord avec les données
expérimentales. Quand, le nombre et l’incertitude de propriétés estimées augmentent,
l’exactitude des modèles diminue.

2- MODELE MATHEMATIQUE

Il n’existe pas a priori de bons ou de mauvais modèle, mais qu’il soit adapté ou non à la
fonction que nous voulons lui faire remplir. Le modèle doit permettre de tester différentes
configurations ou choix des paramètres.
2.1- Modèle dynamique
Dans ce modèle, les équations mathématiques de transport, les termes sources et les
relations de fermeture sont présentés. La forme des équations peut changer d’un composant à
un autre de la pile à combustible, c'est-à-dire chaque élément possède sa propre équation. La
couche de diffusion, la couche de catalyseur et la membrane sont des milieux poreux.
• Equation de continuité:
14èmesJournées Internationales de Thermique JITH2009
27-29 Mars, 2009, Djerba, Tunisie

(
div ρε U = 0)


(1)
• Equation de mouvement :
Où l’équation qui expliquent la résistance de l’écoulement (c’est-à-dire le terme de Darcy).
K
U =− grad ( P )


(2)
ε .µ
• Equation de pression :
A partir des équations (1) et (2) on conclure l’équation (3)
K 
div  grad ( P )  = 0 (3)
µ 
• Equation conservation de espèce
La concentration de chaque espèce (oxygène, hydrogène et eau) peut être calculer a partir de
l’équation de transport qui s’écrire:
∂ε Ck
( )
+ div(ε UCk ) = div Dkeff grad ( Ck ) + SCK


(4)
∂t
Où k : O2 , H 2O, H 2
• Equation d’énergie :
L’énergie est transportée par la conduction et la convection dans les trois phases de la
pile combustible (polymère, liquide, gaz). Les effets des pertes ohmiques dans la membrane et
les effet des partes activation dans les couche de catalyseur sont pris en considération par un
terme additionnel de source dans l'équation d'énergie.
(
∂ ερ c pT ) + div(ερ c UT ) = div ( λ )
grad (T ) + ST

 eff
(5)
∂t
p

Où les terme sources SCK et ST peut être modélisées par le modèle électrochimique.
2.2- Modèles électrochimiques :
Le modèle de l’assemblage retenu est celui présenté par [[Link] et coll ,2000], ce
modèle électrochimique semi-empirique simple et complet est le modèle
généralisé[Link] qui regroupe les modèles (analytique et empirique (1995). Ce dernier
peut être utilisé pour décrire le comportement dynamique d’un assemblage de pile à
combustible PEMFC, il permet de prendre en compte différents paramètres dont leurs
définitions sont essentielles pour un bon résultat, la régulation des performances de
l’assemblage en fonction de la densité de courant, et par conséquent décrire un comportement
proche de la réalité. Cet auteur propose par ailleurs une écriture théorique de la caractérisation
courant-tension d’une PEMFC. Cette caractéristique théorique est alors calée sur des résultats
expérimentaux.
Son objectif est de combiner les approches mécanistique et empirique du modèle. La
première nécessite de nombreux paramètres qui ne sont pas disponibles, la deuxième approche
permet par la caractérisation des corrélations sur les cas particuliers de déterminer dans une
gamme de matériaux définie de présenter un fonctionnement énergétique.

Le potentiel réel d’une pile à combustible décroît par rapport au potentiel de Nernst, à
cause des polarisations qui sont : Polarisation d’activation, Polarisation ohmique et Polarisation
de concentration.
V Pile = E Nernest − ηact − ηohm − ηconc (-6)
V Pile : Tension réelle d’une pile unitaire (Volt)
E Nernest : Potentiel thermodynamique (idéal) de chaque pile et représente la tension
14èmesJournées Internationales de Thermique JITH2009
27-29 Mars, 2009, Djerba, Tunisie

réversible(Volt).
ηact : Polarisation d’activation de l’anode et de la cathode (Volt).
ηohm : Polarisation ohmique de l’anode et de la cathode (Volt).
ηconc : Polarisation de concentration de l’anode et de la cathode (Volt).

• Potentiel idéal :

l’expression de l’équation de Nernst d’après [[Link] et coll , 1995] dans le cas

de la production de l’eau liquide dans la PEMFC.

( ) 1
( )

E Nernest = 1, 229 − 0,85.10−3. (T − 298,15 ) + 4,31.10−5.T . ln PH*2 + ln PO*2 
 2 
(7)

Où T est la température opératoire absolue de la pile, PH*2 et PO*2 sont les pressions partielles à
l’interface respectivement de l’hydrogène et de l’oxygène, dans le cas de la production de l’eau
liquide pour la PEMFC .
Étant donné que nous utiliserons de l’oxygène pur (pas de l’air).
PO*2 = Pc ,cath . 1 − χ sat H 2O  (8)
Pcath Pression de l’oxygène dans l’interface de la surface catalytique de la cathode.
χ sat H O : Fraction molaire de saturation d’eau dans le gaz humidifié.
2

Étant donné que nous utiliserons de l’hydrogène pur qui ne contient pas du monoxyde de
carbone à l’anode.
 1 
PH*2 = Pa ,cath . 1 − χ sat H 2O  (9)
 2 
Pa ,cath Pression de l’hydrogène à l’anode dans de la surface catalytique de la anode.
E Nernest = 1, 229 − 0,85.10−3. (T − 298,15 ) +
  (10)
( ) 1 1
2 2
1
4,31.10−5.T . ln Pcath .(1 − χ sat H 2O ) + ln  Panode .(1 − χ sat H 2O )  

 2

• Les phénomènes de polarisations

1. Polarisation d’activation

Sachant que la polarisation d’activation augmente avec la densité de courant, son

expression peut être écrite sous la forme:
( )
Vact = − ξ1 + ξ 2 .T + ξ3 .T .ln CO2 + ξ 4 .T .ln ( I FC )  (11)
Où I FC est le courant de fonctionnement de la pile et ξ1 , ξ 2 , ξ3 , ξ 4 sont des coefficients
paramétriques pour chaque modèle de pile:
ξ1 = −0.948 (12)
ξ 2 = 0, 00286 + 0, [Link] ( A ) + ( 4, 3.10−5 ) ln ( CH 2
) (13)
ξ3 = 7, 6.10−5 (14)
ξ 4 = −1, 93.10 −4
(15)
Où CO2 est la concentration de l’oxygène dans l’interface de la surface catalytique cathodique
et CH 2 est la concentration de l’hydrogène dans l’interface de la surface catalytique anodique
14èmesJournées Internationales de Thermique JITH2009
27-29 Mars, 2009, Djerba, Tunisie

elle est exprimée en fonction de la pression par la loi de Henry pour la diffusion dans un film
d’eau [34] :
P*O2
CO2 = 498
(16)
−
6
5, 08.10 .e T
P* H 2
CH 2 =  77 
(17)
 
T 
1, 09.106.e

2. Polarisation ohmique

Les pertes Ohmiques dans l’électrolyte obéissent à la loi d’ohm, nous pouvons exprimer les
pertes ohmiques par l’équation suivante :
V ohm = i * r (18)
r est la résistance équivalente de la membrane à la conduction des protons, elle est calculée à
partir de la relation suivante :
r = ρ M .hm (19)
Où hm est l’épaisseur de la membrane, ρ M est la résistance spécifique de la membrane
( Ω .cm ) , elle est obtenue par la relation :
  T 
2
2,5 
181, 6. 1 + 0, 03. ( i ) + 0, 062.   .(i ) 
  303  
ρM = (20)
 T − 303  
λH 2OSO3 − 0, 634 − 3. ( i )  .exp  4,18  
  T 
Où le terme est la résistance spécifique de la membrane ( Ω .cm ) à i = 0 et à une température
de 30 °C (T=303K), le terme exponentiel dans le dénominateur est le facteur de correction de
la température si cette dernière n’est pas à 30 °C, λH 2OSO3 est la Contenance de la membrane en
eau qui admet une valeur minimale de 0 et maximale de 22. A est la surface active de la pile.
l’expression de la polarisation ohmique devient :
  T 
2
2,5 
i .181, 6. 1 + 0, 03. ( i ) + 0, 062.   . ( i )  .hm
  303  
ηohm = (21)
  T − 303  
λH 2OSO3 − 0, 634 − 3. ( i )  .exp  4,18  
  T 

3. Polarisation de concentration (transfert de matière)

Ces pertes sont données par la relation suivante:

R .T  i 
ηcon = − 1 −  (22)
2.F  i max 
i : est densité de courant de fonctionnement.
i max : est la densité de courant maximale.
3. RESULTA :
14èmesJournées Internationales de Thermique JITH2009
27-29 Mars, 2009, Djerba, Tunisie

L’effet des différents paramètres opérationnels tels que la température, la pression la

teneur en eau, l’épaisseur de la membrane, et la surface active de la cathode sur la courbe de
polarisation (courant tension) ont été étudiés en détail. La densité de puissance de la pile a été
déduite de cette courbe.
0.7
1,2 T=353 °K 1.2
0,6
T=333 °K 0.6
P=5 bar
1,0
1.0
T=313 °K 0.5

la puissance (W/cm )
P=3 bar

la puissance (W/cm2)
0,8
la tension (V)

0,4

la tension (V)
0.8 0.4

0,6 P=1 bar

P=5 bar
0.6 0.3

0,4
0,2 0.2
T=353 °K

2
0.4
T=313 °K
P=1 bar
0,2 0.1
T=333 °K P=3 bar
0.2

0,0 0,0 0.0

0,0 0,2 0,4 0,6 0,8 1,0 1,2 1,4 1,6 1,8 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8
2 2
densité de courant (A/cm ) densité de courant (A/cm )

(a) (b)
Figure 1: Courbe de polarisation: (a) effet de la température (b) effet de la pression

1,0
1,4 0,7
1,4 hm=0.125E-1

1,2 λ=14 0,6

1,2 0,8

la puissance (W/cm2)
1,0 0,5 1,0 hm=0.25E-1
la puissant (W/cm )

0,6
la tension (V)
la tension (V)

0,8 λ=10 0,4 0,8 hm=0.125E-1

0,6 λ=12 0,3 0,6 0,4

hm=0.25E-1
0,4 0,2 0,4
2

λ=14 hm=0.5E-1 0,2

0,2
λ=12 0,1 0,2 hm=0.5E-1

λ=10
0,0 0,0 0,0 0,0
0,0 0,2 0,4 0,6 0,8 1,0 1,2 1,4 1,6 1,8 0,0 0,2 0,4 0,6 0,8 1,0 1,2 1,4 1,6
2 2
densité de courant (A/cm ) densité de courant (A/cm )

(a) (b)
Figure 2: Courbe de polarisation: (a) effet de l’épaisseur de la membrane (b) effet de la teneur
en eau dans la membrane

a- Effet de température :
Notre objectif est d’analyser l’importance de l’effet de température sur la performance
de la pile. Une compréhension de base de l’influence directe de la température sur les
différents paramètres du modèle est nécessaire. Les paramètres physiques qui dépendent de la
température sont les suivants :
• La densité de courant d'échange i 0 de la réaction de réduction de l’oxygène augmente
rapidement avec la température à cause de l’accélération de la cinétique de réaction.
• La conductivité de la membrane augmente car une température élevée conduit aussi à
une grande diffusion des protons d‘hydrogène dans la membrane, par la réduction de sa
résistance.
• Bien que l'équation de Nernst montre une diminution du potentiel de référence E 0 avec
une température croissante, les résultats expérimentaux indiquent une augmentation,
14èmesJournées Internationales de Thermique JITH2009
27-29 Mars, 2009, Djerba, Tunisie

qui peut être expliquée avec une diffusivité plus élevée de l'hydrogène avec
l'augmentation de la température.
• Concernant le coefficient de diffusion D, une augmentation de la température mène à
une augmentation de la diffusion des gaz.
Les courbes de polarisation ainsi que la densité de puissance pour différentes températures
sont montrées sur la figure (1.a). La performance augmente avec l'augmentation de la
température de la pile.
b- Effet de la pression :
De manière semblable à la température, la pression de fonctionnement affecte de
nombreux paramètres qui sont important pour le fonctionnement de la pile à combustible parmi
lesquels on cite :
• La composition en gaz d'entrée : Puisque la pression de saturation de la vapeur d'eau
dépend seulement de la température en assumant que les gaz d'entrée sont entièrement
humidifiés. Le changement de la pression de fonctionnement mène à un changement de
la composition en gaz d'entrée.
• Selon l'équation de Nernst, une augmentation de la pression mène à une augmentation
du potentiel E 0 de référence.
• Le coefficient de diffusion D : l'effet de la pression sur la diffusion des gaz est l'inverse
de celui de la température.

L'effet de la pression sur la performance de la pile est représenté sur la figure (1.b). On
remarque que l'influence de ce paramètre est la même que la température. L'amélioration de
puissance avec la variation de la pression n'est pas proportionnelle.

c- Effet de la teneur en eau et l'épaisseur de la membrane :

L’électrolyte de la pile à combustible PEMFC est une membrane polymérique dense
conductrice de protons mince de type acide. Son importance réside dans ses propriétés de
conduction, qui sont conditionnées par :
• De bonnes propriétés de résistance mécanique et chimique dans les conditions de
fonctionnement.
• Imperméables aux gaz.
• Maintient d'humidification adéquate à haute température : cette dernière peut
se faire par l’humidification des gaz réactifs, parce que l’utilisation de gaz sec peut entraîner un
assèchement de la membrane ce qui augmente sa résistance ionique et peut même bloquer le
passage des protons. La fonction d’humidification peut être localisée à l’intérieur ou à
l’extérieur du module et elle peut être réalisée.
1. à l’intérieur : Par passage dans un humidificateur à membrane (auto
humidification),
2. à l’extérieur : Par injection de gouttelettes d’eau directement dans le
courant gazeux.
Les propriétés de la membrane utilisée sont celle de la membrane de type Nafion 117, qui a
pour structure de base poly-tetrafluoéthylène (PTFE), est inerte dans un environnement
oxydant ou réducteur, et peut supporter des pression de 1-5 bars. Les surfaces nécessaires dans
une pile sont fonction de sa puissance.
Les pertes Ohmiques dans l’électrolyte obéissent à la loi d’ohm, nous pouvons les exprimés par
l’équation suivante : ηohm = i * r
Elles sont dues à la résistance que rencontre le flux protonique en traversant l’électrolyte. Ces
pertes peuvent être réduites en diminuant l’épaisseur de l’électrolyte et en améliorant sa
14èmesJournées Internationales de Thermique JITH2009
27-29 Mars, 2009, Djerba, Tunisie

conductivité ionique par l'augmentation de la teneur en eau de la membrane. La conductivité de

celle-ci dépend de ce paramètre selon expression suivante :
Les figures (2.a) et (2.b) représentent l'effet des deux paramètres lies à la membrane; la teneur
en eau et l'épaisseur sur la performance de pile. Augmentation de la teneur en eau conduit à
une amélioration de la densité de puissance, contrairement celle-ci diminue avec des
membrane épaisse.
4. CONCLUSION
Les résultats obtenus par le modèle de la courbe de la polarisation et la densité de puissance
montre que la performance de la pile PEMFC dépend des paramètres suivants :
épaisseur et teneur en eau de la membrane, pression, température et autres ….
La densité de puissance de la pile est importante pour des températures et pressions élevées.
L’effet de la membrane sur la performance est lié à son épaisseur et sa teneur en eau ; les
faibles épaisseurs et grande teneur en eau donne des densités de puissance élevées.

5. REFERENCES BIBLIOGRAPHIQUES
[1] B. Eaton, M.R. Von Spakovsky, M.W. Ellis, D.J. Nelson, B. Olsommer, N. Siegel, “One-
Dimensional, Transient Model of Heat, Mass, and Charge Transfer in a Proton Exchange Membrane,”
Proceedings of IMECE, 2001.

[2] S. Um, C.Y. Wang, K.S. Chen, “Computational Fluid Dynamics Modeling of Proton Exchange
Membrane Fuel Cells,” J. Electrochem. Soc, Vol 147, No. 12, pp. 485- 4493, 2000.

[3] T.E. Springer, S. Gottesfeld, “Pseudohomogeneous Catalyst Layer Model for Polymer Electrolyte
Fuel Cell”, Proc. of the Symposium on Modeling of Batteries and Fuel Cells, R.E. White, M.W.
Verbrugge, and J.F. Stockel, Editors, PV 91-10, pp. 197-208, The Electrochemical Society Softbound
Proceedings Series, Pennington, NJ, (b).

[4] K. Broka, P. Ekdunge, “Modelling the PEM fuel cell cathode,” J. Applied Electrochemistry, Vol.
27, pp. 281-289, 1997.

[5] C. Marr, X. Li, “Composition and Performance modeling of catalyst layer in a proton exchange
membrane fuel cell,” J. Power Sources, Vol. 77, pp. 17-27, 1999.

[6] Z. H. Wang, C. Y. Wang, “Two-Phase Flow and Transport in the Interdigitated Air Cathode of
Proton Exchange Membrane Fuel Cells,” Proceedings ASME Heat Transfer Division, HTD-Vol.366-
1, pp. 27-33, 2000.

[7] D. Natarajan, T. Van Nguyen, “A Two-Dimensional, Two-Phase, Multicomponent, Transient Model

of the Cathode of a Proton Exchange Membrane Fuel Cell Using Conventional Gas Distributors,” J.
Electrochem. Soc., Vol. 148, No. 12, 2001, pp. A1324-A1335.

[8] F. Jaouen, G. Lindbergh, G. Sundholm, “Investigation of Mass-Transport Limitations in the Solid

Polymer Fuel Cell Cathode,” J. [Link]., Vol.149, No.4, pp. A437-A447, 2002.

[9] D. Genevey, M.R. von Spakovsky, M.W. Ellis, D.J. Nelson, B. Olsommer, [Link], N. Montel, N.P.
Siegel, “Transient Model of the Heat, Mass and Charge Transfer as well as Electrochemistry in the
Catalyst Layer of a PEMFC” International Mechanical Engineering Congress and Exposition–
IMECE’2002, ASME IMECE Paper N° 33322, N. Y., N. Y., November, 2002.

[10] D. M. Bernardi, M.W. Verbrugge, “A Mathematical Model of the Solid-Polymer-Electrolyte Fuel

Cell,” J. Electrochem. Soc., Vol. 139, No. 9, pp. 2477-2490.