See discussions, stats, and author profiles for this publication at: https://www.researchgate.

net/publication/357326937

Tutoriel sur l'utilisation du code ANSYS/Fluent pour modéliser une cellule de

pile à membrane échangeuse de protons (PMEP)

Article · December 2021

CITATIONS READS

0 1,249

2 authors:

Yassine Amadane H. Mounir

Mohammed V University of Rabat Mohammadia School of Engineers, Mohammed V university in Rabat
19 PUBLICATIONS   69 CITATIONS    78 PUBLICATIONS   289 CITATIONS   

SEE PROFILE SEE PROFILE

Some of the authors of this publication are also working on these related projects:

Thermo-fluid performance of a (PEFC) with serpentine flow field using Al2O3, SiO2, and ZnO nanofluids View project

modelling and experimenting the ageing of CSP materials View project

All content following this page was uploaded by Yassine Amadane on 25 December 2021.

The user has requested enhancement of the downloaded file.

 Tutoriel sur l’utilisation du code ANSYS/Fluent pour modéliser
une cellule de pile à membrane échangeuse de protons (PMEP)

Y. Amadane,a* H. Mounir,a

a
Research Team EMIsys, Research Center Engineering 3S, Ecole Mohammadia des

ingénieurs ,Université Mohammed V, Rabat , Maroc.

*E-mail : [email protected] Tel : +212663186697

[email protected]

[email protected]
Résumé :

ANSYS est un code de simulation numérique qui permet la réalisation de géométrie et du maillage
jusqu'à l’obtention d’une solution. Le présent tutoriel s’intéressera à la modélisation d’une pile à
membrane échangeuse de protons en utilisant le code ANSYS-Fluent. On se limite à une résolution
sur une cellule, et le système choisi est modélisé en régime stationnaire. Ce travail est destiné aux
chercheurs débutants pour bien assimiler les notions de base de modélisation ainsi la compréhension
des phénomènes physicochimiques qui s’y produisent au sein d’une pile à membrane échangeuse de
protons. Finalement, il faut noter que c’est le premier tutoriel en français qui traite la pile (PMEP) par
ANSYS.

Mots-Clés : ANSYS-Fluent, pile à membrane, régime stationnaire.

 1. Introduction

La pile à membrane échangeuse de protons (PMEP) ( en anglais Proton Exchange Membrane Fuel
Cell- PEMFC) est un système électrochimique qui permet la conversion de l’énergie chimique en
énergie purement électrique avec la production de l’eau et de la chaleur. Les protons libérés dans la
réaction chimique vont traverser la membrane pendant que les électrons sont passés entre les deux
électrodes[1].

Une cellule de pile à membrane est alors composée de plaques d’alimentation et de couches poreuses
et d’une membrane. La membrane joue un rôle très fondamental dans ce type de piles ; elle conduit
les protons et ne laisse jamais passer les électrons et les gaz. Les électrodes sont responsables sur les
réactions chimiques et sur le transport des protons. Les couches poreuses sont responsables sur le
transport des électrons et elles doivent avoir un bon comportement thermique et mécanique. Les
plaques bipolaires lorsqu’il s’agit d’un stack , permettent d’assurer l’écoulement des fluides et de
conduire l’électricité [2] .

Récemment, Plusieurs chercheurs ont utilisé des logiciels numériques pour faire des calculs
bidimensionnels et tridimensionnels de l’écoulement des fluides tout en tenant compte de l’effet de la
température et la viscosité etc. En effet, CFD( Computational fluid dynamic) a été considéré comme
un outil important pour le développement des modèles plus avancés dans différents domaines. Parmi
les différents codes connus actuellement on peut trouver Fluent, CFX, Numeca…etc [3.4].

Le code Fluent permet le calcul des écoulements des fluides et il utilise une interface graphique. Ce
code de simulation contient de nombreux modèles y compris le modèle de la pile à membrane
échangeuse de protons (PMEP) (Fuel Cells and Electrolysis-PEMFC) [5] . Nous présentons les
étapes de calcul à suivre :

 Préparer une géométrie en 3D

 Discrétiser le domaine en volumes élémentaires(Maillage)
 Intégration des équations sur volumes élémentaires, un algorithme itératif est fait jusqu’à la
convergence
 Visualiser les résultats
L’objectif de ce tutoriel est de modéliser une cellule de pile à membrane échangeuse de protons selon
trois dimensions en utilisant le code ANSYS-Fluent. Ce document a pour but de permettre au
chercheurs qui travaillent sur les piles à combustibles de se familiariser avec le code ANSYS-Fluent
c'est-à-dire préparer la géométrie, discrétiser le domaine, choisir le modèle physique et finalement
analyser les résultats.
 2. La création de la géométrie sur le logiciel Solidworks

Dans cette étape nous allons voir comment on peut dessiner la géométrie d’une pile à membrane
échangeuse de protons.

Suivre les instructions suivantes pour préparer la géométrie [7] .

Les étapes à suivre pour préparer le dessin sur SolidWorks

1) Cliquer sur Sketch et sélectionner Top Plane

On sélectionne Rectangle et on fait entrer les paramètres choisis .
Par la suite, on choisit Features , et puis on clique sur Extruded Boss/Base

2) On fait entrer l’épaisseur dans Thickness, et puis on clique sur green checkmark.
Pour créer un canal cliquer sur Sketch Tab et utiliser de nouveau corner Rectangle
Après enter les paramètres du canal et cliquer sur green checkmark

3) Après ces deux étapes, on clique sur insert components puis browse pour ajouter la couche
poreuse (Gas diffusion Layer- GDL)

4) Cliquer sur le bouton mate puis mate selections et finalement sur the vertex of the current
collector on répète la même démarche pour catalyst layer et membrane .

5) Cliquer sur Move component pour faire la rotation , utiliser Delta XYZ puis citer 180
degrees autour de l’axe z et finalement cliquer sur Apply

6) Vérifier que vous avez formé sept parties et enregistrer la géométrie en utilisant l’extension
IGS(.igs)

7) Importer la géométrie sur DesignModeler

Fig.1
Vous trouverez la géométrie montré dans la figure suivante

Fig. 2

Tableau 1 : Les différents noms utilisés pour les parties de la géométrie

Parties de la cellule Appellation

Partie 1 Cathode current collector
Partie 2 Cathode gas diffusion layer
Partie 3 Cathode catalyst layer
Partie 4 Membrane
Partie 5 Anode catalyst layer
Partie 6 Anode gas diffusion layer
Partie 7 Anode current collector
Partie 8 Cathode channel
Partie 9 Anode channel

3. Maillage sur Workbench

Pour ouvrir le logiciel de maillage, il est nécessaire que l’étape précédente soit validée.

Suivre les instructions suivantes pour discrétiser le domaine [8] .

Les étapes à suivre pour faire le maillage

1) Lancer le logiciel de maillage depuis Workbench, en double cliquant sur Maillage il faut
toujours vérifier que l’ étape précédente est validée, voir la figure 5.

2) Identifier et nommer différentes parties de la géométrie pour définir les conditions aux
limites et les conditions sur le maillage
Inlets : inlet_a , inlet_c.
Outlets : outlet_a, outlet_c
Terminals : terminal_a, terminal_c

3) On sélectionne toutes les Sélections Nommées dans l’arborescence , comme vous

remarquez sur la figure 6.

4) On sélectionne Dimensionnement dans Maillage , on choisit le nombre de division puis on

sélectionne les arrêts désirés :

Fig. 3

On répète la même démarche pour tous les arrêts.

Finalement on clique sur Générer un maillage

5) Cliquer sur mettre à jour

Fig. 4
 6) Vérifier dans Workbench que la case maillage est bien validée

Fig. 5 et Fig. 6

Vous trouverez le maillage montré dans la figure suivante

Fig. 7

4. La solution

Premièrement il faut ouvrir le logiciel ANSYS-Fluent : Démarrer/Tous les

programmes/ANSYS/Workbench, puis en double cliquant sur Mécanique des fluides(Fluent)
comme vous voyez dans la figure suivante [9,10]

Fig. 8

Quand vous cliquez sur Configuration , la figure suivante apparaîtra :

Fig. 9

Pour Dimension , Options et Processing Options on coche 3D, Double Precision et Serial
respectivement puis OK .

4.1. Le chargement du modèle

Pour charger le modèle de la pile à membrane échangeuse de protons (PEMFC) on utilise la

commande suivante : /define/models/addon-module et on choisit le chiffre 3 [11].

Fig. 10

Le modèle a été bien chargé comme vous remarquez sur le message suivant

En double cliquant sur Models puis on choisit Fuel Cells and Electrolysis-PEMFC

Fig. 11

Par la suite la figure suivante apparaîtra

Fig. 12

En ce qui concerne parameters on fait entrer les valeurs comme suit :

Fig. 13

En ce qui concerne Anode et Cathode on sélectionne les zones comme suit :

Tableau 2 : Les noms à selectionner pour l’anode et la cathode

Anode Cathode
Current collector anode_cc cathode_cc
Flow Channel anode_flowchannel cathode_flowchannel
Porous Electrode anode_gdl cathode_gdl
TPB Layer (Catalyst) anode_cl cathode_cl

En ce qui concerne Elecrolyte :

Fig. 14

Un simple clique sur Reports :

Fig. 15

On sélectionne terminal_a et terminal_c pour Anode et Cathode respectivement puis OK.

4.2. Les matériaux

Les matériaux sont définis par défaut, en double cliquant sur Materials la figure suivante
apparaîtra.

Fig. 16

4.3. Les conditions de fonctionnement

Dans Cell Zone Conditions on fait entrer les paramètres suivants

Tableau 3 : Les différents types de zones

Zones Types
anode_cc Solid
cathode_cc Solid
anode_cl fluid
cathode_cl fluid
anode_gdl fluid
cathode_gdl fluid
membrane fluid
anode_flowchannel fluid
cathode_flowchannel fluid

Et on définit la pression dans Operating Conditions

P =101325 Pa puis on clique OK pour fermer.

Fig. 17
 4.4. Les conditions aux limites

On se limite à un modèle iso-thermique , Pour l’anode on définit inlet_a, outlet_a et terminal_a .

Inlet_a :

En ce qui concerne Velocity Magnitude 0.3 m/s

Fig. 18

Pour Thermal : T=343 K

Fig. 19

Pour Species :

Fig. 20.

Pour UDS :

En ce qui concerne Water Saturation on sélectionne Specified Value au lieu de Specified Flux.

Fig. 21

terminal_a :

Pour Thermal :

Fig. 22

Pour UDS :

Pour Electric Potential On sélectionne Specified Value au lieu de Specified Flux avec une valeur 0

Fig. 23

Pour la cathode on définit inlet_c, outlet_c et terminal_c .

Il faut suivre les mêmes instructions à l’anode sauf :

inlet_c :

Pour Velocity Magnitude(m/s)= 0.5

Fig. 24
Pour Species et Temperature:

Fig. 25

vous allez entrer la valeur 0.7 pour Electric Potenial

Fig. 26

On définit aussi les contrôles de solution dans Solution Controls. Dans l’onglet Under-Relaxation
Factors vous allez entrer les valeurs suivantes :

Tableau 4 : les quantités physiques et leurs valeurs

Quantité physique Valeur

Pressure 0.7
Density 1
Body Forces 1
Momentum 0.3
h2 1
o2 1
h2o 1
Energy 1
Electric Potential 1
Protonic Potenial 1
Water Saturation 1
Water Content 0.95

Dans l’onglet Advanced Solution Controls vous allez suivre les instructions suivantes :

Fig. 27

Dans l’onglet Monitors vous allez entrer les valeurs suivantes :

Tableau 5 : les résidus et leurs valeurs

Résidus Valeurs
continuity 0.001
x-velocity 0.001
y-velocity 0.001
z-velocity 0.001
energy 10-6
h2 0.001
h2 0.001
h2o 0.001
uds-0 0.001
uds-1 0.001
uds-2 0.001
uds-3 0.001

Fig. 28
Dans ce cas d’une simulation stationnaire on choisit le critère de convergence absolute . Nous avons
choisit 10-3 [12].

Dans Surface Monitors on clique sur Create la figure suivante apparaîtra, et on suit les instructions
suivantes :

Fig. 29

Dans l’onglet Solution Initialisation on clique sur Standard Initialisation et Relative to Cell Zone
puis Initialize et finalement OK comme vous voyez dans la figure suivante :

Fig. 30

Dans l’onglet Calculation Activities on clique sur Edit puis Data File Quantities et finalement on
sélectionne Standard Quantities et Additional Quantities.

Fig. 31

Pour Run Calculation, on fait entrer les valeurs suivantes :

Fig. 32

On clique sur Calculate.

5. Lancement du calcul et traitement des résultats

Dans cette étape on va voir les instructions à suivre pour visualiser les résultats dans Post-Traitement

5.1. La convergence

On suit l’evolution des résidus en fonction de nombre d’itérations

Fig. 33a et Fig. 33b

Vous remarquez que les calculs sont bien convergés : 100 itérations

5.2. La courbe de polarisation

On change la valeur du voltage dans terminal_c 0.7 , 0.8 et 0.9 respectivement ( Les conditions
aux limites) sur la case Electric Potential puis on note la valeur du densité de courant [13].
 Tableau 6 : Le voltage en fonction du densité de courant

La densité de courant Voltage

0.145371 0.7
La solution converge à 100 itérations
0.050094 0.8
La solution converge à 132 itérations
0.003763 0.9
La solution converge à 216 itérations

Fig.34

Vous pouvez changer la valeur de la température et faire la même démarche , pour visualiser par
exemple l’effet du changement de température sur les performances de la pile.

Fig. 35

Pour générer cette courbe il suffit de multiplier les valeurs de tensions par les valeurs de courants

5.3. La visualisation avec Post-CFD

Créer un contour :

On clique sur Graphics and Animations puis Contours et finalement on sélectionne toutes les
surfaces comme vous observez sur la figure suivante [11,14] :

Fig. 36

On choisit Species et puis Mass Fraction of o2

Fig. 37a

On choisit Mass Fraction of h2

Fig.37b

Nous allons voir par la suite quelques outils essentiels pour traiter les résultats.

On clique sur Résultats dans Workbench

Fig. 38

Cette fenêtre apparaîtra

 Fig. 39

Créer un contour :

On clique sur contour puis on écrit contour 1.

Créer un plan :

On clique sur Location pour créer un plan, et on forme trois plans Plan 1(Z=0mm) , Plan
2(Z=20mm) et Plan 3(Z=39.99 mm).

Fig. 40a

Dans contour 1 on sélectionne les trois plans sur Location puis on clique OK :

Fig. 40b , Fig. 40c

On choisit Temmperature dans Variable et Global dans Range on clique Apply

Pour visualiser la température (Temperature) on clique sur Apply et vous allez obtenir la figure
suivante :

Fig. 41

6. Conclusion

Ce tutoriel a pour but de permettre aux chercheurs qui travaillent sur les piles à combustibles de
s’initier à la modélisation par le code ANSYS-Fluent. Dans ce travail nous avons présenté les étapes
à suivre pour préparer la géométrie , discrétiser le domaine (maillage ) et finalement visualiser les
résultats tels que la courbe de polarisation , la fraction massique de l’hydrogène et l’oxygène et la
distribution de la température au sein de la pile à membrane échangeuse de protons. Les erreurs
rencontré par un débutant lors de la modélisation et l’utilisation des UDF pour cette pile seront le sujet
de notre prochain travail.

7. Références

[1] Arvay, A. ; Ahmed, A. ; Peng,A.; Kannan, A.M. , 2011. Convergence criteria establishment for
3D simulation of proton exchange membrane fuel cell. International Journal of Hydrogen Energy,
37(3), pp.2482–2489. Available at: http://dx.doi.org/10.1016/j.ijhydene.2011.11.005.

[2] Larminie, J. ; Dicks A. Fuel cell systems explained. 2nd ed. Chichester, UK: John Wiley & Sons;
2003.

[3] Cheddie, DF ; Munroe NDH. Three dimensional modeling of high temperature PEM fuel cells. J
Power Sources 2006;160:215e23.
 [4] Shah, AA ; Luo, KH ; Ralph, TR ; Walsh FC. Recent trends and developments in polymer
electrolyte membrane fuel cell modeling. Electrochim Acta 2011;56:3731e57.

[5] Zhang H, Pei P, Li P, Yuan X. The conception of in-plate adverse-flow flow field for a proton
exchange membrane fuel cell. Int J Hydrogen Energy 2010;35:9124e33.

[6] Al-Baghdadi MARS. Three-dimensional computational fluid dynamics model of a tubular-shaped

PEM fuel cell. Renewable Energy 2008;33:1334e45.

[7] ANSYS, Inc. ANSYS FLUENT 12.0 fuel cells module manual,
https://www1.ansys.com/customer/content/documentation/ 121/fluent/flfuelcells.pdf; April 2009.

[8] Siegel, C. Review of computational heat and mass transfer modeling in polymer-electrolyte-
membrane (PEM) fuel cells. Energy 2008;33:1331e52.

[9] Wang, L. ; Husa,r A. ; Zhou, T. ; Liu, H. A parametric study of PEM fuel cell performances. Int
J Hydrogen Energy 2003;28: 1263e72.

[10] Kamarajugadda, S. ; Mazumder, S. On the implementation of membrane models in

computational fluid dynamics calculations of polymer electrolyte membrane fuel cells. Comput Chem
Eng 2008;32:1650e60.

[11] Dawes, JE ; Hanspal, NS ; Family, OA ; Turan, A. Threedimensional CFD modelling of PEM

fuel cells: an investigation into the effects of water flooding. Chem Eng Sci 2009;64:2781e94.

[12] Cano-Andrade, S. ; Hernandez-Guerrero, A., Spakovsky MR, Damian-Ascencio CE, Rubio-

Arana JC. Current density and polarization curves for radial flow field patterns applied to PEMFCs.
Energy 2010;35:920e7.

[13] Min, C. Performance of a proton exchange membrane fuel cell with a stepped flow field design.
J Power Sources 2008;186: 370e6.

[14] Djilali, N. Computational modelling of polymer electrolyte membrane (PEM) fuel cells:
challenges and opportunities. Energy 2007;32:269e80.
Fig. 1. la géométrie sur DesignModeler

Fig. 2. la géométrie d’une pile à membrane échangeuse de protons

 Fig. 3. Dimensionnement des arrêts

Fig. 4. L’onglet Mettre à jour

Fig. 5. L’étape précédente est bien validée

Fig. 6. L’interface du maillage- Sélections nommées

Fig. 7. Le maillage obtenu

Fig. 8. L’interface Workbench

Fig. 9. L’interface Fluent Launcher
Fig. 10. L’interface Fluent

Fig. 11. Le chargement du modèle

Fig. 12. Les modèles donnés par Fluent

Fig. 13. Les caractéristiques d’une pile à membrane échangeuse de protons

Fig. 14. Les paramètres à définir pour l’électrolyte
Fig. 15. Interface pour les rapports

Fig. 16. Les matériaux donnés par Fluent

Fig. 17. La pression de fonctionnement

Fig. 18. La vitesse à l’entrée de l’anode

Fig. 19. La température à l’entrée de l’anode

Fig. 20. Les fractions massiques de l’eau et de l’hydrogène à l’entrée de l’anode

Fig. 21. Les paramètres UDS
Fig. 22. La température à la paroi terminal_a

Fig. 23. Les paramètres UDS

Fig. 24. La vitesse à l’entrée de la cathode

Fig. 25. La température et les fractions massiques de l’eau et de l’oxygène à l’entrée de la cathode
Fig. 26. Les paramètres UDS à la paroi terminal_c
Fig. 27. Contrôles avancés
Fig. 28. Les paramètres de convergence
Fig. 29. Les paramètres à définir pour Surface Monitors
Fig. 30. L’initialisation de la solution

Fig. 31. Les quantités physiques dans Calculation Activities

Fig. 32. Interface pour le lancement de la simulation

Fig. 33a. L’intégrale en fonction des itérations

Fig. 33b. L’évolution des résidus

Fig. 34. La courbe de polarisation

1

0,8
Voltage(V)

0,6

0,4

0,2

0
0 0,05 0,1 0,15 0,2
2
La densité de courant(A/cm )
 Fig. 35. La densité de puissance en fonction de la densité de courant

0,12

La densité de puissance
0,1
0,08
0,06
0,04
0,02
0
0 0,02 0,04 0,06 0,08 0,1 0,12 0,14 0,16
2
La desnité de courant (A/cm )

Fig. 36 .L’interface Graphics and Animations

Fig. 37a .La fraction massique de l’oxygène

 Fig. 37b .La fraction massique de l’hydrogène

Fig. 38 . Wokbench—Résultats

Fig. 39 . L’interface CFD-Post

Fig. 40a . La création des plans

Fig. 40b . L’interface Location Selector

Fig. 40c . Les détails du contour 1

 Fig. 41 . La distribution de la température

View publication stats