Modélisation multiphysique des PEMFC

Une pile à combustible à membrane échangeuse de protons n'est pas gouvernée par un seul phénomène physique. Sa réponse électrique résulte de l'interaction entre l'électrochimie, le transport des gaz, la gestion de l'eau, les transferts thermiques, la dynamique des pressions et les lois de commande.

Pour cette raison, un modèle d'ingénierie utile doit être multiphysique par conception.

L'objectif n'est pas de reproduire chaque mécanisme microscopique. Il s'agit de conserver les phénomènes couplés qui influencent réellement les grandeurs d'intérêt, tout en maintenant un modèle identifiable, numériquement robuste et calculable à un coût raisonnable.

Cette page décrit la structure de modélisation que j'utilise pour mes travaux de jumeau numérique PEMFC. Pour une introduction générale à la technologie, voir Piles à combustible à membrane échangeuse de protons.

Le modèle est formulé à l'échelle du système et du stack plutôt qu'à celle des pores ou de la microstructure.

Les principaux domaines modélisés sont :

• la conversion électrochimique ;
• l'hydratation de la membrane ;
• le transport des gaz et de la vapeur ;
• l'accumulation d'eau liquide ;
• la dynamique des pressions à l'anode et à la cathode ;
• le comportement thermique ;
• les états de fonctionnement liés aux purges ;
• la sortie électrique.

Le modèle obtenu est destiné à la simulation, à l'estimation d'état, à l'analyse d'ingénierie et, à terme, aux applications orientées commande.

La tension de cellule peut être représentée schématiquement par

V_cell = E_rev - η_act - η_ohm - η_conc

La tension réversible dépend de la température et des pressions partielles des réactifs. Les trois pertes principales correspondent à :

• des pertes d'activation associées à la cinétique des réactions ;
• des pertes ohmiques dominées par le transport protonique dans la membrane et par les résistances électroniques ;
• des pertes de concentration ou de transport de matière à forte densité de courant.

Cette décomposition ne sert pas uniquement à reproduire une courbe de polarisation. Chaque terme établit un lien physique entre le comportement électrique observable et les états internes thermiques, d'hydratation et de transport.

L'hydratation de la membrane est un état interne central, car la conductivité protonique dépend fortement de la teneur en eau.

Un modèle réduit suit donc une variable d'hydratation de la membrane, typiquement écrite sous la forme d'un état de teneur en eau tel que λ.

Sa dynamique reflète plusieurs mécanismes concurrents :

• production d'eau à la cathode ;
• entraînement électro-osmotique de l'anode vers la cathode ;
• rétrodiffusion à travers la membrane ;
• échanges avec les phases gazeuses ;
• évacuation par les débits de sortie.

L'hydratation influe directement sur la résistance de la membrane et donc sur la perte de tension ohmique.

Une membrane très sèche entraîne une résistance élevée et de mauvaises performances. Un excès d'eau peut toutefois favoriser l'inondation et les limitations de transport. L'état recherché n'est donc pas simplement « davantage d'eau ».

Les pressions et compositions des réactifs évoluent dynamiquement avec le débit d'entrée, la consommation par réaction, le débit de sortie, l'humidification et les événements de purge.

Une formulation réduite par volume de contrôle peut s'écrire génériquement :

dn_i/dt = ṅ_i,in - ṅ_i,out - ṅ_i,reaction + ṅ_i,transfer

Les pressions partielles correspondantes découlent de l'inventaire gazeux, du volume, de la température et de l'équation d'état.

Ces états sont importants car le modèle électrochimique dépend des pressions partielles d'oxygène et d'hydrogène, et non simplement des débits imposés de l'extérieur.

La gestion de l'eau exige de distinguer au moins deux formes :

• la vapeur d'eau transportée dans la phase gazeuse ;
• l'eau liquide condensée.

Le passage de l'une à l'autre dépend de la température, de la pression et des conditions de saturation.

Un état réduit d'eau liquide peut représenter l'accumulation sans chercher à résoudre individuellement les gouttes ou les canaux. Son rôle est de capturer l'effet d'un excès d'eau sur le transport gazeux et le comportement en fonctionnement.

Ce point est particulièrement important pour les systèmes à cathode ouverte ou refroidis par air, dans lesquels la gestion thermique, l'alimentation en air, l'humidification et l'évacuation de l'eau sont fortement couplées.

La température influe sur presque tous les sous-systèmes :

• cinétique des réactions ;
• tension réversible ;
• conductivité de la membrane ;
• pression de saturation ;
• masse volumique des gaz ;
• condensation et évaporation de l'eau ;
• évacuation de la chaleur.

Un bilan thermique global peut s'écrire :

C_th dT/dt = Q̇_gen - Q̇_loss

La chaleur produite comprend les pertes électrochimiques irréversibles et les autres dissipations internes. Le rejet thermique peut inclure convection, flux d'air forcé, conduction et échanges thermiques avec les structures environnantes.

Le modèle thermique ferme donc plusieurs boucles de rétroaction importantes plutôt que de constituer un calcul de post-traitement indépendant.

Le fonctionnement d'une anode en cul-de-sac ou purgée périodiquement introduit une échelle de temps de fonctionnement spécifique.

Entre deux purges, des espèces inertes et de l'eau peuvent s'accumuler. Une purge modifie brutalement l'inventaire gazeux et affecte donc :

• la concentration d'hydrogène ;
• la pression ;
• l'humidité ;
• l'eau accumulée ;
• le comportement en tension qui suit.

Un modèle pratique peut représenter ce phénomène par un état continu décrivant l'accumulation depuis la dernière purge, associé à des événements discrets de purge.

On obtient naturellement un modèle dynamique hybride contenant à la fois des états continus et des transitions déclenchées par événements.

Le point important n'est pas l'existence des sous-modèles pris séparément, mais leur couplage.

Par exemple :

1. le courant détermine la consommation de réactifs et la production d'eau ;
2. la production d'eau modifie l'hydratation de la membrane et l'humidité des gaz ;
3. l'hydratation modifie la conductivité de la membrane ;
4. la conductivité modifie la perte de tension ohmique ;
5. les pertes génèrent de la chaleur ;
6. la température modifie la pression de saturation et la cinétique des réactions ;
7. la condensation modifie l'accumulation d'eau liquide ;
8. l'eau liquide accumulée peut limiter le transport de matière ;
9. la perte de concentration qui en résulte modifie la tension aux bornes.

Un modèle qui traite ces effets indépendamment peut reproduire des courbes isolées tout en échouant sur la dynamique.

Un vecteur d'état réduit peut contenir des grandeurs telles que

x = [T  λ  p_H₂  p_O₂  p_H₂O  s_l  q_p]ᵀ

Les variables d'état sont :

• T : température de la cellule, qui agit sur la cinétique électrochimique, la conductivité protonique de la membrane et la gestion de l'eau ;
• λ : teneur en eau de la membrane, généralement définie comme le nombre de molécules d'eau par groupe acide sulfonique (SO₃H). Elle influence fortement la conductivité protonique de la membrane ;
• p_H₂ : pression partielle d'hydrogène à l'anode, représentant la disponibilité locale du combustible pour la réaction d'oxydation de l'hydrogène ;
• p_O₂ : pression partielle d'oxygène à la cathode, représentant la disponibilité locale de l'oxydant pour la réaction de réduction de l'oxygène ;
• p_H₂O : pression partielle de vapeur d'eau, décrivant la quantité d'eau présente dans la phase gazeuse et donc les conditions d'humidification de la cellule ;
• s_l : saturation en eau liquide dans les milieux poreux, représentant la fraction du volume poreux occupée par l'eau liquide. Elle est particulièrement pertinente pour décrire l'inondation et le transport diphasique ;
• q_p : état d'accumulation lié à la purge, représentant l'accumulation progressive d'espèces ou d'eau entre deux purges successives de l'anode.

Ensemble, ces variables décrivent les principaux états dynamiques thermiques, d'hydratation, de phase gazeuse, d'eau liquide et électrochimiques du modèle PEMFC.

Toutes les grandeurs physiques ne doivent pas nécessairement devenir des états estimés.

Le choix des états est un compromis entre :

• pertinence physique ;
• observabilité ;
• identifiabilité des paramètres ;
• coût de calcul ;
• conditionnement numérique ;
• mesures disponibles.

Dans les travaux actuels d'estimation d'état, le vecteur estimé est volontairement plus petit que l'état physique complet.

Les entrées typiques du modèle comprennent :

• le courant de charge ;
• les conditions des gaz à l'entrée ;
• les conditions ambiantes ;
• la commande de refroidissement ou du ventilateur ;
• la commande de purge ;
• les consignes de fonctionnement.

Les sorties mesurables typiques comprennent :

• la tension du stack ou de la cellule ;
• la température ;
• la pression ;
• les mesures de débit lorsqu'elles sont disponibles.

Les autres grandeurs internes sont des états du modèle et non des signaux directement mesurés.

L'implémentation est développée en Python.

Les exigences principales sont :

• définitions explicites des états ;
• évaluation déterministe du modèle ;
• unités physiques clairement définies ;
• intégration numérique stable ;
• états bornés ou physiquement contraints ;
• tests de régression ;
• jeux de paramètres reproductibles ;
• séparation entre les équations du modèle et le code d'analyse.

La qualité du logiciel scientifique est directement importante ici. Un modèle numériquement fragile peut produire des courbes physiquement plausibles tout en restant inadapté à l'estimation d'état, à l'optimisation ou à la commande.

Un modèle multiphysique contient des paramètres de natures très différentes :

• paramètres géométriques ;
• propriétés thermophysiques ;
• coefficients électrochimiques ;
• coefficients de transport ;
• termes correctifs empiriques ;
• paramètres d'actionneurs et de capteurs.

Ils ne doivent pas être traités de la même manière.

Certains sont connus par la géométrie ou les propriétés des matériaux. D'autres doivent être identifiés à partir de données. Certains ne sont que des paramètres effectifs compensant la résolution spatiale réduite du modèle.

Cette distinction doit rester explicite, car l'ajustement de paramètres ne transforme pas une mauvaise structure de modèle en modèle physiquement valide.

Un modèle physique complet n'est pas automatiquement un meilleur modèle d'ingénierie.

Une résolution spatiale élevée peut être utile pour une conception détaillée, mais contre-productive pour l'estimation temps réel ou la commande.

La réduction vise donc à préserver le comportement entrée-sortie dominant et les états internes physiquement significatifs tout en supprimant les détails spatiaux inutiles.

Le modèle obtenu doit être assez complexe pour représenter les couplages pertinents, mais assez simple pour être analysé et exécuté de manière répétée.

La validation expérimentale n'est pas encore achevée pour le jumeau numérique actuel.

Avant toute comparaison aux données expérimentales, l'implémentation peut néanmoins être contrôlée systématiquement au moyen de :

• la cohérence dimensionnelle ;
• les lois de conservation ;
• les cas limites ;
• les contrôles de monotonie lorsque la physique l'exige ;
• la sensibilité aux paramètres ;
• la convergence numérique ;
• les tests de régression ;
• la comparaison de formulations équivalentes ;
• la cohérence interne des bilans couplés.

Ces contrôles relèvent de la vérification, et non de la validation. Ils établissent que le modèle est implémenté de manière cohérente avec ses équations et ses hypothèses ; ils ne prouvent pas qu'il reproduit le système réel.

Le jumeau numérique ne se réduit pas au modèle multiphysique direct.

Une architecture complète exige aussi :

• des interfaces de mesure ;
• une gestion des paramètres ;
• une estimation d'état ;
• une gestion de l'incertitude ;
• la comparaison entre modèle et comportement du système ;
• un état de fonctionnement persistant ;
• à terme, une logique de diagnostic ou de commande.

Le modèle multiphysique constitue le cœur physique de cette architecture.

L'étape suivante est la couche d'estimation, où des mesures incomplètes sont combinées aux prédictions du modèle pour reconstruire des états internes comme la température, l'hydratation de la membrane et l'accumulation liée aux purges.

Voir aussi : Estimation d'état et jumeaux numériques.

Les principaux principes que j'applique sont :

• conserver la signification physique chaque fois que possible ;
• rendre les hypothèses explicites ;
• séparer les grandeurs mesurées, les états et les paramètres ;
• éviter d'ajuster des paramètres non identifiables ;
• tester les lois de conservation et les comportements limites ;
• préférer des jeux de paramètres reproductibles au réglage manuel ;
• distinguer vérification et validation expérimentale ;
• conserver un modèle assez simple pour être compris et inspecté.