Brevet FR3131664A1 — publié le 20230707

Méthode de détermination d’un indicateur de dégradation d’une cellule d’une pile à combustible (PEM) comprenant plusieurs cellules ayant chacune une surface active, la méthode comprenant les étapes de : recevoir des mesures de grandeurs physiques caractéristiques de conditions opératoires de la pile à combustible, recevoir une mesure expérimentale de la tension aux bornes d’une cellule, appliquer un algorithme Bayésien non linéaire (UKF) à partir du vecteur de commande u constitué des mesures de grandeurs physiques et de la mesure expérimentale, pour déterminer une estimée, à l’instant t+1, d’un vecteur d’état de la cellule défini par une première variable d’état caractéristique d’un taux de dissolution de platine et d’au moins une deuxième variable d’état caractéristique d’un taux de couverture d’oxydation de platine, l’algorithme Bayesien non linéaire étant en outre défini par un modèle d’observation électrochimique de la tension aux bornes d’une cellule fonction desdites variables d’état, déterminer un indicateur de dégradation de la cellule Figure 3 Méthode de détermination d’un indicateur de dégradation d’une pile à combustible, basée sur un observateur Bayésien hybride L’invention concerne le domaine des piles à combustible, en particulier les piles à combustible pour alimenter des véhicules électriques, notamment les piles à combustible à membrane d’échange de protons ou piles PEMFC (Proton Exchange Membrane Fuel Cell). Les piles à combustible sont des systèmes qui permettent de convertir de l’énergie chimique contenue dans une réaction d’oxydo-réduction en énergie électrique. Ce sont des systèmes ouverts qui échangent de la matière (le combustible et le comburant) et de l’énergie sous forme de chaleur et d’’électricité avec l’extérieur. L’invention concerne une méthode permettant de déterminer un indicateur de dégradation d’une pile à combustible, la méthode étant basée sur un observateur bayésien hybride, couplant un modèle de tension avec des données issues de capteurs, et prenant en compte à la fois des phénomènes de dégradation irréversibles et réversibles. Plusieurs mécanismes de dégradation réversibles et irréversibles interviennent dans une cellule d’une pile à combustible. Ces dégradations affectent la plupart des composants d’une cellule. En particulier, le phénomène de dégradation impacte la couche active d’une cellule de sorte que la surface active peut diminuer au cours du temps, impactant le rendement global du dispositif. La couche active d’une cellule englobe plusieurs éléments. Une couche active d’une cellule est organisée en un réseau microporeux constitué d’un ionomère, d’un support de catalyseur et d’un catalyseur. L’ionomère assure la conductivité protonique. Le support de catalyseur peut être du noir de carbone ou plus généralement du graphite poreux. Il assure la conductivité électronique et constitue un milieu poreux pour le transport des gaz et de l’eau liquide. Le catalyseur est le plus souvent à base de platine. Il est déposé sous forme de nanoparticules et permet d’augmenter la cinétique de la réaction électrochimique. Il y a deux couches actives de part et d’autre de la membrane d’une cellule, l’une constituant l’anode et l’autre la cathode. La notion de surface active se distingue de la surface géométrique en ce qu’elle définit la surface développée totale qui participe aux réactions chimiques. Ainsi, une dégradation de la surface active du catalyseur intervient lorsqu’on observe une perte de surface active de platine ou une altération des propriétés du catalyseur. La dégradation de la surface active du catalyseur peut être réversible, par exemple lorsqu’il s’agit d’une oxydation du platine, ou bien irréversible, par exemple lorsqu’il s’agit de phénomènes de dissolution-redéposition de platine. Un problème à résoudre dans ce contexte consiste à pouvoir estimer et prédire le niveau de dégradation de la surface active d’une cellule à partir de mesures accessibles sur la pile et de modèles de dégradation. L’art antérieur comporte différentes méthodes permettant d’estimer l’évolution du vieillissement d’une pile à combustible. La publication scientifique ”Yunjin Ao, Salah Laghrouche, Daniel Depernet, and Kui Chen. Lifetime prediction for proton exchange membrane fuel cell under real driving cycles based on platinum particle dissolve model. International Journal of Hydrogen Energy,45(56):32388–32401, 2020” fait état de travaux lors desquels les auteurs sélectionnent la surface active de l’électrode (ECSA) en tant qu’indicateur de santé. Il est établi dans cette publication qu’une pile atteint sa fin de vie lorsque sa tension diminue de 10% par rapport à sa valeur initiale. Donc, les pertes de tension (puissance sortie) et les pertes de surface active sont liées pour déterminer la durée de vie restante. L’inconvénient majeur de cette approche est la nécessité d’avoir des données de surface active mesurées pour définir le modèle de dégradation à partir de courbes de calibration. Les deux documents suivants proposent un autre type de méthode : ”Kui Chen, Salah Laghrouche, and Abdesslem Djerdir. Fuel cell health prognosis using unscented kalman filter: Postal fuel cell electric vehicles case study. International Journal of Hydrogen Energy, 44(3):1930–1939, 2019” ”Mathieu Bressel, Mickael Hilairet, Daniel Hissel, and Belkacem Ould Boua-mama. Fuel cell remaining useful life prediction and uncertainty quantification under an automotive profile. In IECON 2016-42nd Annual Conference of the IEEE Industrial Electronics Society, pages 5477–5482. IEEE, 2016” Dans ces travaux, les auteurs utilisent un filtre de Kalman non linéaire du type « Unscented Kalman Filter » couplé à un modèle empirique de dégradation du niveau de tension. L’approche est validée à l’aide de données expérimentales mesurées sur une pile destinée à l’automobile. Deux inconvénients majeurs se dégagent de ces travaux: la tension est la seule variable sélectionnée comme indicateur de santé, ce qui limite la robustesse et la précision de l’estimateur le modèle de dégradation utilisé est un modèle purement empirique et lors de l’étape de pronostic (prédiction du vieillissement), l’approche présente des pertes de performance pour suivre la dynamique de vieillissement de la tension. La demande de brevet FR2928777 propose un dispositif et un procédé d'estimation du vieillissement de pile à combustible dû à l'empoisonnement au monoxyde de carbone, grâce à l’utilisation d’un observateur d’état de type Kalman. Ce système est basé sur une mesure de surtension anodique, et le contrôle de la purge des monoxydes de carbone stockés. L’idée proposée ici affiche des limitations, car elle ne permet pas de dissocier précisément dégradations réversibles et irréversibles. L’invention propose une approche hybride basée à la fois sur des mesures et sur un modèle prenant en compte à la fois les dégradations réversibles et irréversibles. L’algorithme proposé est de type Bayésien non linéaire et repose sur plusieurs modèles physiques de dissolution et d’oxydation du platine dont est composée la surface active de la cellule. Un observateur d’état est conçu pour prédire des variables d’états particulières liées au vieillissement de la pile sur la base d’un modèle d’état, puis pour corriger ces prédictions en utilisant de manière optimale les mesures réalisées sur la pile. L’étape de prédiction utilise un modèle physique d’observation qui relie les variables mesurées sur la pile aux variables d’état du modèle d’observateur. Ce cycle prédiction-correction permet d’obtenir une solution d’estimation précise et robuste. La solution proposée permet de déterminer l’état de santé courant (diagnostic) et l’état de santé futur (pronostic) en temps réel des systèmes de pile à combustible à membrane d’échange de protons. Contrairement aux solutions connues, l’invention utilise un modèle robuste capable de prendre en compte les phénomènes conduisant à une perte de performances liées à la dégradation réelle des matériaux du cœur de pile mais aussi les pertes de performances induites par une modification temporaire (réversible) des propriétés physiques de ces matériaux. Cette approche présente l’avantage d’une meilleure prédiction de la durée de vie réelle de la pile. Les modèles physiques utilisés comprennent un modèle de prédiction basé sur des lois de dégradation irréversibles du catalyseur de la pile comme la dissolution du platine mais aussi des lois de dégradation réversibles comme l’oxydation du platine. Le modèle d’observation est un modèle de la tension aux bornes d’une cellule qui permet de relier les données mesurées aux variables d’état prédites. L’observateur d’état réalise la correction optimale du modèle physique par traitement des données capteurs. L’observateur suit l’état de santé de la pile via deux indicateurs : la tension cellule et la surface active du catalyseur. Cette information estimée par l’algorithme permet d’estimer la durée de vie restante de la pile et de gérer de manière optimale les conditions opératoires du système complet. L’invention a pour objet une méthode, mise en œuvre par ordinateur, de détermination d’un indicateur de dégradation d’une cellule d’une pile à combustible comprenant une pluralité de cellules ayant chacune une surface active, la méthode comprenant les étapes de, à chaque nouvel instant t : recevoir un ensemble de mesures de grandeurs physiques caractéristiques de conditions opératoires de la pile à combustible, cet ensemble de mesures formant un vecteur de commande u, recevoir une mesure expérimentale de la tension aux bornes d’une cellule, appliquer un algorithme Bayésien non linéaire à partir du vecteur de commande u et de la mesure expérimentale, pour déterminer une estimée, à l’instant t+1, d’un vecteur d’état de la cellule défini par une première variable d’état caractéristique d’un taux de dissolution de platine et d’au moins une deuxième variable d’état caractéristique d’un taux de couverture d’oxydation de platine, l’algorithme Bayesien non linéaire étant en outre défini par un modèle d’observation électrochimique de la tension aux bornes d’une cellule fonction desdites variables d’état, déterminer un indicateur de dégradation de la cellule à partir de l’estimée du vecteur d’état de la cellule. Selon un aspect particulier de l’invention, l’algorithme Bayesien non linéaire est un filtre de Kalman UKF. Selon un aspect particulier de l’invention, le vecteur d’état est composé : d’une première variable d’état caractéristique d’un taux de dissolution de platine, d’une deuxième variable d’état caractéristique du taux de couverture d’oxydation de platine lié à la formation d’oxyde de platine, d’une troisième variable d’état caractéristique du taux de couverture d’oxydation de platine lié à la formation d’hydroxyde de platine Selon un aspect particulier de l’invention, le modèle d’observation de la tension aux bornes d’une cellule dépend d’un modèle de surtension fonction d’au moins les variables d’état. Selon un aspect particulier de l’invention, les variables d’état sont déterminées récursivement à partir du vecteur de commande et de modèles de prédiction chimiques de l’effet de dissolution du platine et de l’effet d’oxydation du platine. Selon un aspect particulier de l’invention, l’indicateur de dégradation de la cellule est un indicateur de dégradation de la surface active électro-chimique de la cellule qui dépend de la première variable d’état. Selon un aspect particulier de l’invention, les mesures formant le vecteur de commande u sont prises parmi des mesures de courant, de température de refroidissement, de stochiométrie, d’humidité relative, de pression. Selon un aspect particulier de l’invention, le filtre de Kalman réalise au moins les opérations suivantes : Calculer la différence entre la mesure expérimentale de la tension et la valeur estimée de la tension via le modèle d’observation, Calculer une innovation comme le résultat de la différence pondérée par un gain de Kalman propre à chaque variable d’état, Corriger l’estimation du modèle de prédiction des variables d’état par l’innovation. Selon un aspect particulier de l’invention, l’indicateur de dégradation est déterminé pour l’ensemble des cellules de la pile à combustible. L’invention a aussi pour objet un programme d'ordinateur comportant des instructions pour l'exécution de la méthode de détermination d’un indicateur de dégradation d’au moins une cellule d’une pile à combustible selon l’invention, lorsque le programme est exécuté par un processeur. L’invention a aussi pour objet un support d'enregistrement lisible par un processeur sur lequel est enregistré un programme comportant des instructions pour l'exécution de la méthode de détermination d’un indicateur de dégradation d’au moins une cellule d’une pile à combustible selon l’invention, lorsque le programme est exécuté par un processeur. L’invention a aussi pour objet un système pour déterminer un indicateur de dégradation d’au moins une cellule d’une pile à combustible comprenant une pluralité de cellules, le système comprenant plusieurs capteurs pour mesurer différentes grandeurs physiques caractéristiques du fonctionnement de la pile à combustible et un calculateur pour exécuter la méthode de détermination d’un indicateur de dégradation selon l’invention. Selon un aspect particulier de l’invention, le calculateur est apte à déterminer un état de santé de la pile à combustible à partir de l’indicateur de dégradation déterminé pour l’ensemble des cellules de la pile à combustible. D’autres caractéristiques et avantages de la présente invention apparaîtront mieux à la lecture de la description qui suit en relation aux dessins annexés suivants. la représente un schéma d’une pile à combustible, la représente un schéma général de l’approche hybride proposée par l’invention, la représente un schéma de principe du fonctionnement de l’invention basé sur un filtre de Kalman non linéaire UKF, La représente, sur un schéma, une pile à combustible PAC composée de plusieurs cellules C montées en série. Chaque cellule comprend une membrane MEM disposée entre deux plaques bipolaires PB1, PB2. A chaque extrémité opposée de la pile à combustible, une plaque monopolaire PM1, PM2 est pourvue de canaux CD de distribution d’hydrogène et d’oxygène pour alimenter la pile. Sur la droite de la , on a représenté une seule cellule élémentaire C de la pile à combustible PAC. Une cellule C comprend les éléments suivants agencés en série : une première plaque bipolaire PB1, un premier support S1, un premier catalyseur CAT1, une membrane MEM, un second catalyseur CAT2, un second support S2, une seconde plaque bipolaire PB2. Une pile à combustible PAC est soumise à une dégradation au cours du temps principalement du fait de la dégradation des membranes MEM. L’ensemble des cellules de la pile à combustible PAC est traversée par un même courant, cependant la densité locale de courant peut ne pas être uniforme dans une cellule et peut avoir une répartition spatiale différente d’une cellule à l’autre du fait de la dégradation potentielle des membranes MEM au cours du temps. La représente un autre schéma d’un exemple de pile à combustible et de son intégration dans un système hybride d’évaluation de l’état de santé de la pile selon un mode de réalisation de l’invention. Sur la gauche de la , on a représenté un autre schéma d’une cellule d’une pile à combustible PEM constituée des éléments déjà décrits à la à savoir : une membrane MEM entourée par deux électrodes EL1,EL2 de sorte à réaliser un assemblage membrane électrodes AME. De chaque côté de cet assemblage est placé une couche de diffusion poreuse S1,S2. L’ensemble est placé entre deux plaques d’alimentation bipolaires PB1,PB2. L’invention a notamment pour objectif d’estimer la dégradation du catalyseur à la surface des électrodes et en particulier les pertes de surface active des électrodes notamment les pertes de platine. Sur la droite de la , on a représenté schématiquement le fonctionnement de la méthode selon l’invention qui est basée sur un modèle d’observation 200 définit d’une part par un vecteur d’état x qui suit une évolution modélisée par une fonction f, qui dépend d’un vecteur de commande u et d’autre part une variable d’observation y qui suit une évolution définie par la fonction h qui dépend du vecteur d’état et du vecteur de commande. La prédiction y’ de la variable d’observation fournie par le modèle d’observation 200 est comparée à une mesure expérimentale de la même grandeur y sur la pile PEM et l’erreur e(y-y’), appelée innovation, calculée est utilisée pour corriger le vecteur d’état. Le vecteur de commande u est un vecteur de mesures de différentes paramètres de fonctionnement de la pile PEM. Le modèle d’observation 200 est résolu au moyen d’un algorithme Bayésien non linéaire, par exemple un filtre de Kalman non linéaire du type SR-UKF (Square Root Unscented Kalman Filter). Dans un mode de réalisation, la méthode selon l’invention vise à estimer la perte de surface active électro-chimique d’une électrode. Pour cela, l’invention repose en partie sur le choix de trois variables d’état à estimer qui sont : la dissolution directe du rayon des particules de platine ( r dissPt ), la formation d’oxydes de platine PtO ( θ PtO ) et la formation d’hydroxydes de platine PtOH ( θ PtOH ). La dissolution du platine est un phénomène irréversible tandis que l’oxydation est un phénomène de dégradation réversible. La prise en compte de ces deux phénomènes dans le modèle d’observation permet de mieux estimer les variations de surface active et leur impact sur la surtension aux bornes de la cellule. En effet, les électrodes EL1,EL2 sont composées de particules de platine Pt supportées par une surface en carbone. Les particules de platine Pt peuvent être dissoutes en présence de fortes densités de courant qui traversent la pile notamment durant les cycles de recharge de la pile, notamment pour des batteries automobiles. Cette dégradation irréversible augmente les pertes d’activation et de concentration, ce qui engendre une dégradation du fonctionnement de la pile. L’oxydation du platine PtO affecte la dégradation catalytique et donc directement le phénomène de dissolution du platine. La formation d’oxyde de platine PtO apparait généralement pour des valeurs de tension aux bornes d’une cellule supérieures à 0.8 V. La dissolution directe des particules de platine est décrite par l’équation chimique (a) qui correspond à une dégradation irréversible. L’oxydation du platine Pt peut être décrite selon un mécanisme en quatre étapes décrit par les équations chimiques (b),(c),(d),(e). L’équation (b) correspond à une oxydation réversible, l’équation (c) correspond à une oxydation irréversible, l’équation (d) correspond à une réduction irréversible et enfin l’équation (e) correspond à une adsorption irréversible de l’oxygène. (a) (b) (c) (d) (e) Ces trois variables d’état constituent ensemble un vecteur d’état x(t) et leurs évolutions sont modélisées à partir de modèles physico-chimiques qui décrivent la cinétique des réactions chimiques (a)-(e) ci-dessus. La dissolution du platine est modélisée par l’équation (1) suivante : (1) V diss est la vitesse de dissolution du platine qui dépend d’un certain nombre de paramètres PARAM relatives aux conditions de fonctionnement de la pile. M Pt est la masse molaire du platine ρ Pt est la densité du platine L’équation (1) est issue de la publication de ”Andres Jacome, Daniel Hissel, Vincent Heiries, Mathias Gérard, and Sébastien Rosini. Pemfc state-of-health estimation using a model-based state bayesian observer under an automotive load profile. In 2020 IEEE Vehicle Power and Propulsion Conference (VPPC), pages 1–4. IEEE, 2020”. La formation d’hydroxyde de platine PtOH est modélisée par l’équation (2) suivante : (2) La formation d’oxyde de platine PtO est modélisée par l’équation (3) suivante : (3) r 1 , r 2 et r 3 sont des cinétiques de réaction qui dépendent du potentiel local. Le taux de couverture du platine , autrement dit le pourcentage de la surface active occupé par le platine, est donné par les relations suivantes : avec Les équations (2) et (3) sont issues de la publication ”Thomas Jahnke, Georg A Futter, Andrea Baricci, Claudio Rabissi, and Andrea Casalegno. Physical modeling of catalyst degradation in low temperature fuel cells: Platinum oxidation, dissolution, particle growth and platinum band formation. Journal of The Electrochemical Society”. Un modèle d’observation (4) de tension cellule dérivé de l’équation de Butler-Volmer est utilisé pour corriger l’estimation. Ce modèle est décrit dans la publication ”C Robin, M Gérard, M Quinaud, J d’Arbigny, and Y Bultel. Proton exchange membrane fuel cell model for aging predictions: Simulated equivalent active surface area loss and comparisons with durability tests. Journal of Power Sources, 326:417–427, 2016”. Le modèle d’observation est donné par la relation (4) : (4) U c ,est est une estimée de la tension aux bornes de la cellule. est une tension réversible qui peut être déterminée à partir d’un bilan thermochimique à partir des paramètres de résistance R, température T et de pression P. R ohm est une résistance totale déterminée à partir d’une résistance de la membrane et d’une valeur fixe. I est le courant traversant la cellule. est une surtension qui dépend des variables d’état x(t) et de paramètres relatifs aux conditions opératoires de la pile, température T, courant I, pression P notamment. La représente, sur un schéma fonctionnel, le principe de mise en œuvre de la méthode de détermination d’un indicateur de dégradation de la pile selon un mode de réalisation de l’invention. La méthode est basée sur un algorithme Bayésien non linéaire, par exemple un filtre de Kalman non linéaire du type « SR-UKF » tel que décrit dans la publication ” Shichun Yang, Song Feng, Kangfeng Sun, Shuai Wang, and Yaoguang Cao. Square-root unscented kalman filter for ammonia coverage ratio and input ammonia estimations in diesel-engine urea-scr system. ISA transactions, 96:299–308, 2020 ”. Le fonctionnement de la méthode est itératif. A un instant t, le filtre reçoit en entrée un vecteur de commande u constituée de mesures de paramètres physiques de la pile représentant ses conditions opératoires. Le vecteur de commande u comporte, par exemple, une mesure du courant I traversant une cellule de la pile, une mesure de température de refroidissement en entrée de la pile T cool,in , une mesure de température de refroidissement en sortie de la pile T cool,out , une mesure de stochiométrie de l’anode Sta, une mesure de stochiométrie de la cathode Stc, une mesure d’humidité relative de l’anode RHa, une mesure d’humidité relative de la cathode RHc, une mesure de pression. D’autres mesures de paramètres influant le fonctionnement de la pile peuvent être ajoutés au vecteur de commande u. D’autre part, le filtre reçoit également en entrée une mesure y=U cell de la tension aux bornes d’une cellule. Le vecteur de commande est exploité pour construire les équations d’état f(x) données par les relations (1), (2) et (3) qui modélisent les transitions du vecteur d’état entre un instant k et un instant k+1. Un terme de bruit Wc est pris en compte dans les équations d’état pour caractériser les incertitudes. Les variables d’état estimées a priori sont sommées avec une variable d’innovation multipliée par un gain de Kalman K non constant qui évolue avec le temps, pour obtenir le vecteur d’état x(t) a posteriori. L’innovation est obtenue en calculant l’erreur entre la mesure y de la tension aux bornes de la cellule et une estimée de cette tension calculée à partir de l’équation d’observation h(x) donnée par la relation (4). Un terme de bruit Wm est pris en compte dans l’équation d’observation pour caractériser les incertitudes. Un indicateur de dégradation de la surface active électrochimique ECSA(t) est ensuite calculé à partir du vecteur d’état x() a posteriori, par exemple à l’aide de la formule suivante ; (5) S 0 est la surface active initiale à l’instant k=0, N part est le nombre initial de particules de platine sur la surface active. N part peut être déterminé à partir de la valeur initiale ou nominale de la surface active de l’électrode S(0) et de la valeur initiale du rayon des particules de platine . L’indicateur de dégradation de la surface active ECSA(t) dépend essentiellement du rayon des particules de platine qui sont affectées par la dissolution du platine. Cependant, la prise en compte de l’oxydation du platine via les deux autres variables d’état permet de prendre en compte l’erreur liée à l’oxydation sur la surtension . En effet la surtension peut être modélisée en fonction de différents paramètres dont le taux d’oxydation du platine comme indiqué par la relation (6) (6) Les coefficients β 1 , β 2 , β 3 , β 4 , β 5 , β 6 ,β 7 sont des paramètres du modèle. , sont des pressions moyennes de l’air, de l’hydrogène et de l’eau est la conductivité protonique de la couche active au niveau de l’électrode. Les différents paramètres intervenant dans les équations d’état et dans l’équation d’observation, ainsi que les paramètres du filtre de Kalman doivent être calibrés dans une phase d’initialisation préalable de l’algorithme. Le filtre de Kalman peut fonctionner en boucle fermée pour fournir un diagnostic en temps réel de l’état de santé de la pile. Le calcul de l’innovation via l’équation d’observation, autrement dit l’erreur entre la tension mesurée et la tension estimée, permet de corriger les variables d’état estimées à priori. Après une phase de convergence, le filtre de Kalman peut également fonctionner en boucle ouverte. Dans ce cas, la tension estimée via l’équation d’observation n’est plus calculée et les variables d’état ne sont plus corrigées par l’innovation pondérée par le gain de Kalman. Le filtre de Kalman réalisé alors un pronostic via une estimation a priori de la perte de surface active ECSA. L’invention peut être mise en œuvre en tant que programme d’ordinateur comportant des instructions pour son exécution. Le programme d’ordinateur peut être enregistré sur un support d’enregistrement lisible par un processeur. La référence à un programme d'ordinateur qui, lorsqu'il est exécuté, effectue l'une quelconque des fonctions décrites précédemment, ne se limite pas à un programme d'application s'exécutant sur un ordinateur hôte unique. Au contraire, les termes programme d'ordinateur et logiciel sont utilisés ici dans un sens général pour faire référence à tout type de code informatique (par exemple, un logiciel d'application, un micro logiciel, un microcode, ou toute autre forme d'instruction d'ordinateur) qui peut être utilisé pour programmer un ou plusieurs processeurs pour mettre en œuvre des aspects des techniques décrites ici. Les moyens ou ressources informatiques peuvent notamment être distribués (" Cloud computing "), éventuellement selon des technologies de pair-à-pair. Le code logiciel peut être exécuté sur n'importe quel processeur approprié (par exemple, un microprocesseur) ou cœur de processeur ou un ensemble de processeurs, qu'ils soient prévus dans un dispositif de calcul unique ou répartis entre plusieurs dispositifs de calcul (par exemple tels qu’éventuellement accessibles dans l’environnement du dispositif). Le code exécutable de chaque programme permettant au dispositif programmable de mettre en œuvre les processus selon l'invention, peut être stocké, par exemple, dans le disque dur ou en mémoire morte. De manière générale, le ou les programmes pourront être chargés dans un des moyens de stockage du dispositif avant d'être exécutés. L'unité centrale peut commander et diriger l'exécution des instructions ou portions de code logiciel du ou des programmes selon l'invention, instructions qui sont stockées dans le disque dur ou dans la mémoire morte ou bien dans les autres éléments de stockage précités. Le code exécutable peut également être téléchargeable depuis un serveur distant. Le programme d’ordinateur peut comprendre un code source, un code objet, un code source intermédiaire ou un code objet partiellement compilé ou toute autre forme d’instructions de code de programme adaptées pour mettre en œuvre l’invention sous la forme d’un programme d’ordinateur. Un tel programme peut présenter diverses architectures fonctionnelles. Par exemple, un programme d’ordinateur selon l’invention peut être décomposé en une ou plusieurs routines qui peuvent être adaptées à exécuter une ou plusieurs fonctions de l’invention telles que décrites précédemment. Les routines peuvent être enregistrées ensemble dans un même fichier exécutable mais peuvent également être sauvegardées dans un ou plusieurs fichiers externes sous la forme de librairies qui sont associées à un programme principal de façon statique ou dynamique. Les routines peuvent être appelées depuis le programme principal mais peuvent également comprendre des appels à d’autres routines ou sous-routines. Tous les procédés ou étapes de procédés, programmes ou sous-programmes décrits sous la forme d’organigrammes doivent être interprétés comme correspondant à des modules, segments ou portions de code de programme qui incluent une ou plusieurs instructions de code pour implémenter les fonctions logiques et les étapes de l’invention décrites. Alternativement, l’invention peut aussi être mise en œuvre au moyen d’un processeur qui peut être un processeur générique, un processeur spécifique, un circuit intégré propre à une application (connu aussi sous le nom anglais d’ASIC pour « Application-Specific Integrated Circuit ») ou un réseau de portes programmables in situ (connu aussi sous le nom anglais de FPGA pour « Field-Programmable Gate Array »). La technique de l'invention peut se réaliser sur une machine de calcul reprogrammable (un processeur ou un micro contrôleur par exemple) exécutant un programme comprenant une séquence d'instructions, ou sur une machine de calcul dédiée (par exemple un ensemble de portes logiques comme un FPGA ou un ASIC, ou tout autre module matériel). Les différents indicateurs de dégradation calculés au cours du temps peuvent être restitués à un utilisateur via un écran d’affichage ou toute autre interface graphique. L’invention peut aussi être implémentée dans un organe d’aide à la maintenance de manière à réaliser un outil permettant d’automatiser la maintenance d’une pile à combustible à bord d’un véhicule. Une alerte de maintenance peut être déclenchée lorsqu’un indicateur de dégradation dépasse ou diminue en dessous d’un seuil d’alerte prédéterminé. L’invention permet également de déterminer un indicateur global de l’état de santé de la pile à partir de l’ensemble des indicateurs de dégradation de chaque cellule. Méthode, mise en œuvre par ordinateur, de détermination d’un indicateur de dégradation d’une cellule d’une pile à combustible (PEM) comprenant une pluralité de cellules ayant chacune une surface active, la méthode comprenant les étapes de, à chaque nouvel instant t : recevoir un ensemble de mesures de grandeurs physiques caractéristiques de conditions opératoires de la pile à combustible, cet ensemble de mesures formant un vecteur de commande u, recevoir une mesure expérimentale de la tension aux bornes d’une cellule, appliquer un algorithme Bayésien non linéaire (UKF) à partir du vecteur de commande u et de la mesure expérimentale, pour déterminer une estimée, à l’instant t+1, d’un vecteur d’état de la cellule défini par une première variable d’état caractéristique d’un taux de dissolution de platine et d’au moins une deuxième variable d’état caractéristique d’un taux de couverture d’oxydation de platine, l’algorithme Bayesien non linéaire étant en outre défini par un modèle d’observation électrochimique de la tension aux bornes d’une cellule fonction desdites variables d’état, déterminer un indicateur de dégradation de la cellule à partir de l’estimée du vecteur d’état de la cellule. Méthode de détermination d’un indicateur de dégradation selon la revendication 1 dans laquelle l’algorithme Bayesien non linéaire est un filtre de Kalman UKF. Méthode de détermination d’un indicateur de dégradation selon l’une quelconque des revendications précédentes dans laquelle le vecteur d’état est composé : d’une première variable d’état caractéristique d’un taux de dissolution de platine, d’une deuxième variable d’état caractéristique du taux de couverture d’oxydation de platine lié à la formation d’oxyde de platine, d’une troisième variable d’état caractéristique du taux de couverture d’oxydation de platine lié à la formation d’hydroxyde de platine Méthode de détermination d’un indicateur de dégradation selon l’une quelconque des revendications précédentes dans laquelle le modèle d’observation de la tension aux bornes d’une cellule dépend d’un modèle de surtension fonction d’au moins les variables d’état. Méthode de détermination d’un indicateur de dégradation selon l’une quelconque des revendications précédentes dans laquelle les variables d’état sont déterminées récursivement à partir du vecteur de commande et de modèles de prédiction chimiques de l’effet de dissolution du platine et de l’effet d’oxydation du platine. Méthode de détermination d’un indicateur de dégradation selon l’une quelconque des revendications précédentes dans laquelle l’indicateur de dégradation de la cellule est un indicateur de dégradation de la surface active électro-chimique de la cellule qui dépend de la première variable d’état. Méthode de détermination d’un indicateur de dégradation selon l’une quelconque des revendications précédentes dans laquelle les mesures formant le vecteur de commande u sont prises parmi des mesures de courant, de température de refroidissement, de stochiométrie, d’humidité relative, de pression. Méthode de détermination d’un indicateur de dégradation selon l’une quelconque des revendications précédentes dans laquelle le filtre de Kalman réalise au moins les opérations suivantes : Calculer la différence entre la mesure expérimentale de la tension et la valeur estimée de la tension via le modèle d’observation, Calculer une innovation comme le résultat de la différence pondérée par un gain de Kalman propre à chaque variable d’état, Corriger l’estimation du modèle de prédiction des variables d’état par l’innovation. Méthode de détermination d’un indicateur de dégradation selon l’une quelconque des revendications précédentes dans laquelle l’indicateur de dégradation est déterminé pour l’ensemble des cellules de la pile à combustible. Programme d'ordinateur comportant des instructions pour l'exécution de la méthode de détermination d’un indicateur de dégradation d’au moins une cellule d’une pile à combustible selon l’une quelconque des revendications précédentes, lorsque le programme est exécuté par un processeur. Support d'enregistrement lisible par un processeur sur lequel est enregistré un programme comportant des instructions pour l'exécution de la méthode de détermination d’un indicateur de dégradation d’au moins une cellule d’une pile à combustible selon l’une quelconque des revendications 1 à 9, lorsque le programme est exécuté par un processeur. Système pour déterminer un indicateur de dégradation d’au moins une cellule d’une pile à combustible comprenant une pluralité de cellules, le système comprenant plusieurs capteurs pour mesurer différentes grandeurs physiques caractéristiques du fonctionnement de la pile à combustible et un calculateur pour exécuter la méthode de détermination d’un indicateur de dégradation selon l’une quelconque des revendications 1 à 8. Système selon la revendication 12 dans lequel le calculateur est apte à déterminer un état de santé de la pile à combustible à partir de l’indicateur de dégradation déterminé pour l’ensemble des cellules de la pile à combustible.