FR 3041479 A1 20170324 FR 1558898 A 20150921 1 DETERMINATION D'UNE DISTRIBUTION SPATIALE DE L'ACTIVITE CATALYTIQUE D'UNE ELECTRODE DE CELLULE ELECTROCHIMIQUE DOMAINE TECHNIQUE Le domaine de l'invention est celui des réacteurs électrochimiques comportant un empilement de cellules électrochimiques, tels que les piles à combustible et les électrolyseurs, et plus précisément celui des procédés de détermination d'un paramètre représentatif de l'activité catalytique d'au moins une électrode de la cellule électrochimique ainsi que celui des procédés de réalisation d'une électrode de cellule électrochimique. ÉTAT DE LA TECHNIQUE ANTÉRIEURE Un réacteur électrochimique tel qu'une pile à combustible ou un électrolyseur comporte habituellement un empilement de cellules électrochimiques qui comprennent chacune une anode et une cathode séparées électriquement l'une de l'autre par un électrolyte, dans laquelle a lieu une réaction électrochimique entre deux réactifs qui sont introduits de manière continue. Dans le cas d'une pile à hydrogène, le combustible (l'hydrogène) est apporté au contact de l'anode, alors que le comburant (l'oxygène), par exemple contenu dans de l'air, est amené au contact de la cathode. La réaction électrochimique est subdivisée en deux demi-réactions, une réaction d'oxydation et une réaction de réduction, qui ont lieu respectivement à l'interface anode/électrolyte et à l'interface cathode/électrolyte. Pour avoir lieu, la réaction électrochimique requiert la présence d'un conducteur ionique entre les deux électrodes, à savoir l'électrolyte éventuellement contenu dans une membrane polymère, et un conducteur électronique formé par le circuit électrique externe. L'empilement de cellules est ainsi le lieu de la réaction électrochimique : les réactifs doivent y être apportés, les produits et les espèces non réactives doivent en être évacués, tout comme la chaleur produite lors de la réaction. 3041479 2 Les cellules sont habituellement séparées les unes des autres par des plaques bipolaires qui assurent l'interconnexion électrique des cellules. Les plaques comportent un circuit de distribution du carburant formé sur une face anodique et un circuit de distribution du comburant formé sur une face cathodique opposée à la face anodique. Chaque circuit de 5 distribution prend souvent le forme d'un réseau de canaux agencés par exemple en parallèle ou en serpentin, adaptés à amener l'espèce réactive de manière homogène à l'électrode correspondante. Les plaques bipolaires peuvent également comporter un circuit de refroidissement formé d'un réseau de conduits internes qui assurent l'écoulement d'un fluide caloporteur permettant d'évacuer la chaleur produite localement 10 lors de la réaction par la cellule. Chaque électrode est habituellement formée d'une couche active et d'une couche de diffusion. La couche de diffusion est située entre la couche active et la plaque bipolaire correspondante. Elle est réalisée en un matériau poreux permettant la diffusion des espèces réactives gazeuses jusqu'à la couche active ainsi que la diffusion des produits issus 15 de la réaction électrochimique. La couche active, située entre la membrane électrolytique et la couche de diffusion, est le lieu de la réaction électrochimique. Elle comporte généralement un conducteur électronique (par exemple du carbone), un conducteur protonique (l'électrolyte) et un catalyseur, par exemple des particules de platine. Elle peut être formée d'un mélange de poudre de carbone platinée, de particules de PTFE et 20 d'électrolyte polymère (par exemple du Nafion), disposé sur un tissu de carbone. La couche active se caractérise par une activité catalytique qui traduit la capacité du catalyseur à catalyser une réaction (ici la réduction de l'oxygène). Elle peut être exprimée en termes d'activité spécifique et d'activité massique du catalyseur, qui correspondent à la densité de courant échangée par la cellule en fonctionnement normalisé respectivement 25 par la surface active de catalyseur et par la masse de catalyseur. L'activité catalytique dépend notamment de la surface active du catalyseur et de sa charge (également appelée taux de chargement ou teneur). On cherche habituellement à préserver la durée de vie des cellules électrochimiques, par exemple en uniformisant la distribution spatiale de la densité de courant produit par la 3041479 3 cellule en fonctionnement. En effet, les inhomogénéités spatiales de densité de courant entraînent des inhomogénéités spatiales de température de la cellule, qui peuvent conduire, d'une part à augmenter la cinétique des réactions de dégradations des différents éléments de la cellules et d'autre part à générer des contraintes mécaniques d'origine 5 thermique susceptibles de fragiliser la tenue mécanique des éléments de la cellule. Ainsi, le document FR2976732 décrit une cellule électrochimique réalisée de manière à obtenir un échauffement local homogène au sein de la cellule en fonctionnement. L'échauffement dépend de la densité de courant électrique en chaque point de la cellule qui est elle-même fonction de la pression partielle d'espèce réactive. En effet, considérant 10 ici le côté cathodique de la cellule, la quantité d'oxygène contenu dans le gaz circulant dans le circuit de distribution diminue progressivement à mesure que l'oxygène est consommé par la cellule, ce qui entraîne une variation spatiale de la densité de courant électrique produit par la cellule, et donc une variation spatiale de l'échauffement de la cellule. Pour éviter cette inhomogénéité spatiale de l'échauffement de la cellule, la conductivité 15 électrique entre la plaque bipolaire amenant l'oxygène et la cellule est ajustée localement de manière à compenser la diminution de la pression partielle d'oxygène. Dans le document W02007/032903, la durée de vie d'une cellule électrochimique est préservée par l'utilisation d'une couche active dont la répartition spatiale de l'activité catalytique permet d'obtenir une distribution spatiale uniforme de la densité de courant, 20 qui, en retour, se traduit par une distribution spatiale uniforme de la température de la cellule en fonctionnement. En particulier, ce document décrit une répartition spatiale d'activité catalytique qui varie entre l'entrée et la sortie des circuits de distribution des espèces réactives. 25 EXPOSÉ DE L'INVENTION L'invention a pour objectif de remédier au moins en partie aux inconvénients de l'art antérieur, et plus particulièrement de proposer un procédé de détermination de la distribution spatiale d'un paramètre représentatif de l'activité catalytique d'une électrode 3041479 4 de cellule électrochimique permettant d'homogénéiser la température locale de la cellule électrochimique en fonctionnement et ainsi de préserver la durée de vie de cette dernière. L'invention propose à cet effet un procédé de détermination d'une distribution spatiale d'un paramètre d'intérêt représentatif d'une activité catalytique d'une couche active d'au 5 moins une électrode d'une cellule électrochimique, ladite cellule comportant deux électrodes séparées l'une de l'autre par un électrolyte et disposées entre deux plaques bipolaires adaptées à amener des espèces réactives aux électrodes et à évacuer de la chaleur produite par la cellule en fonctionnement, au moins une électrode comportant une couche active comportant au moins un catalyseur, la couche active étant adaptée à générer 10 une réaction électrochimique à partir de l'une des espèces réactives, ledit procédé comprenant les étapes dans lesquelles : on définit une cellule électrochimique au sein de laquelle le paramètre d'intérêt est réparti spatialement selon une distribution initiale ; on définit une distribution spatiale d'une température de consigne au sein de la 15 cellule en fonctionnement ; on obtient une distribution spatiale d'une première grandeur thermique représentative de la température ou de l'évacuation de la chaleur au sein de la cellule en fonctionnement ; on estime la distribution spatiale d'une seconde grandeur thermique à partir de 20 ladite distribution spatiale de la température de consigne et de ladite distribution spatiale de la première grandeur thermique ; on détermine une distribution spatiale du paramètre d'intérêt en fonction de la distribution spatiale de la seconde grandeur thermique. De préférence, le paramètre d'intérêt est la charge en catalyseur de la couche active. 25 De préférence, la distribution spatiale du paramètre d'intérêt est déterminée de sorte qu'elle présente une valeur moyenne sensiblement identique à une valeur moyenne de la distribution initiale. De préférence, la première grandeur thermique est une température locale effective mesurée au sein de la cellule en fonctionnement, et la seconde grandeur thermique est 3041479 5 une grandeur représentative d'un écart local entre la température effective et la température de consigne. De préférence, l'étape v) comporte : une sous-étape de calcul d'une distribution spatiale d'un coefficient correcteur à 5 partir de la distribution spatiale de la seconde grandeur thermique ; une sous-étape dans laquelle on détermine la distribution spatiale du paramètre d'intérêt par corrélation de la distribution spatiale initiale du paramètre d'intérêt avec la distribution spatiale du coefficient correcteur. De préférence, l'étape v) comporte : 10 une sous-étape dans laquelle on identifie au moins une zone Zi de la cellule où la seconde grandeur thermique présente une valeur locale estimée supérieure à une valeur seuil prédéterminée ; une sous-étape dans laquelle on détermine la distribution spatiale du paramètre d'intérêt par modification de la distribution spatiale initiale dans la zone Zi 15 identifiée, en fonction de la valeur locale estimée de la seconde grandeur thermique dans cette zone. De préférence, la première grandeur thermique est représentative de l'évacuation locale de la chaleur produite par la cellule en fonctionnement, et la seconde grandeur thermique est représentative de la production locale de chaleur par la cellule en fonctionnement. 20 De préférence, la première grandeur thermique est le débit local effectif mesuré d'un fluide caloporteur circulant dans un circuit de refroidissement d'une plaque bipolaire de la cellule, et la seconde grandeur thermique est le flux local de chaleur produite par la cellule en fonctionnement. L'étape v) peut comporter : 25 une première sous-étape dans laquelle on estime une distribution spatiale de la densité d'un signal électrique produit par la cellule en fonctionnement, à partir de ladite distribution spatiale estimée du flux de chaleur produite, et 3041479 6 - une deuxième sous-étape dans laquelle on détermine la distribution spatiale du paramètre d'intérêt à partir de ladite distribution spatiale de la densité du signal électrique. L'invention porte également sur un procédé de réalisation d'une électrode de cellule 5 électrochimique, comportant les étapes dans lesquelles : - on considère une cellule électrochimique de référence comportant deux électrodes séparées l'une de l'autre par un électrolyte et disposées entre deux plaques bipolaires adaptées à amener des espèces réactives aux électrodes et à évacuer de la chaleur produite par la cellule en fonctionnement, au moins une électrode 10 comportant une couche active adaptée à générer une réaction électrochimique à partir de l'une des espèces réactives, la couche active présentant un paramètre d'intérêt représentatif d'une activité catalytique répartie spatialement selon une distribution initiale ; - on détermine, par le procédé selon l'une quelconque des caractéristiques 15 précédentes, une distribution spatiale du paramètre représentatif de l'activité catalytique ; on réalise ladite électrode de sorte que le paramètre représentatif de l'activité catalytique de la couche active de celle-ci présente la distribution spatiale déterminée. 20 L'étape de réalisation de ladite électrode peut comporter le dépôt, dans au moins une zone de la couche active préalablement identifiée à partir de ladite distribution spatiale déterminée du paramètre représentatif de l'activité catalytique, d'une encre catalytique présentant une charge catalytique inférieure à une charge catalytique moyenne hors de ladite zone. 25 L'encre catalytique peut être déposée par sérigraphie au moyen d'un écran de sérigraphie comportant au moins une ouverture obturée destinée à recouvrir ladite zone de la couche active. L'invention porte également sur une électrode de cellule électrochimique comportant une couche active adaptée à générer une réaction électrochimique au sein de la cellule, 3041479 7 caractérisée en ce qu'elle présente un paramètre représentatif d'une activité catalytique de la couche active, de préférence une charge catalytique, qui présente, dans au moins une zone de la couche active, une valeur locale inférieure à une valeur moyenne hors de ladite zone. 5 BRÈVE DESCRIPTION DES DESSINS D'autres aspects, buts, avantages et caractéristiques de l'invention appara Iront mieux à la lecture de la description détaillée suivante de formes de réalisation préférées de celle-ci, donnée à titre d'exemple non limitatif, et faite en référence aux dessins annexés sur 10 lesquels : la figure la est une représentation schématique en coupe d'un exemple de cellule électrochimique, et la figure lb est une représentation schématique illustrant la relation de corrélation entre la distribution spatiale de la production de chaleur et la distribution spatiale d'évacuation de chaleur dont résulte la distribution spatiale de température de la 15 cellule électrochimique en fonctionnement ; la figure 2 est un organigramme d'un procédé de détermination de la distribution spatiale de la charge catalytique d'une cellule électrochimique selon un premier mode de réalisation ; la figure 3 est un organigramme d'un procédé de détermination de la distribution spatiale 20 de la charge catalytique d'une cellule électrochimique selon un deuxième mode de réalisation ; la figure 4 est un exemple de maillage du circuit de refroidissement où chaque maille comporte un terme de production locale de chaleur In, un terme d'évacuation local de chaleur Dfi, et une température locale T15 ; 25 la figure 5 est un exemple de relation entre la densité de puissance d'une pile à combustible en fonction de la charge catalytique des cellules électrochimiques ; 3041479 8 la figure 6a est un exemple de distribution spatiale de la température effective d'une cellule électrochimique de référence en fonctionnement pour laquelle la charge catalytique des couches actives des électrodes est spatialement uniforme, et la figure 6b est un exemple de distribution spatiale de la température effective de la cellule électrochimique pour 5 laquelle la charge catalytique des couches actives présente une distribution spatiale déterminée selon le procédé du premier mode de réalisation ; les figures 7a à 7c illustrent trois étapes du dépôt d'une encre catalytique sur la surface d'une couche active de cellule électrochimique au moyen d'un écran de sérigraphie. 10 EXPOSÉ DÉTAILLÉ DE MODES DE RÉALISATION PARTICULIERS Sur les figures et dans la suite de la description, les mêmes références représentent des éléments identiques ou similaires. On décrira différents modes de réalisation et variante en référence à une pile à combustible, et en particulier à une pile à hydrogène de type PEM (pour Proton Exchange 15 Membrane, en anglais) dont la cathode est alimentée en oxygène et l'anode en hydrogène. L'invention s'applique cependant à tout type de pile à combustible, en particulier à celles fonctionnant à basse température, c'est-à-dire à une température inférieure à 200°C, ainsi qu'aux électrolyseurs électrochimiques. 20 La figure la illustre de manière partielle et schématique un exemple de cellule électrochimique 1 appartenant à un empilement de cellules d'une pile à combustible de type PEM. La cellule 1 comporte une anode 10 et une cathode 20 séparées l'une de l'autre par un électrolyte comprenant ici une membrane polymère 30, les électrodes 10, 20 étant disposées entre deux plaques bipolaires 40, 50 adaptées à amener des espèces réactives 25 aux électrodes et à évacuer la chaleur produite lors de la réaction électrochimique. Les plaques bipolaires comportent un circuit de distribution d'hydrogène 41 situé sur une face anodique et un circuit de distribution d'oxygène 51 situé sur une face cathodique. Elles sont ici formées de deux tôles métalliques 42a, 42b ; 52a, 52b, assemblées l'une à l'autre 3041479 9 par des zones ou points de soudure, et sont embouties de manière à former les circuits de distribution. L'agencement des bossages permet de réaliser également un circuit de refroidissement 43, 53 à l'intérieur des plaques assurant l'écoulement d'un fluide caloporteur sans que celui-ci ne vienne au contact des électrodes. D'autres technologies 5 de plaque bipolaire peuvent être utilisées, par exemple les plaques réalisées en matériau composite chargé par exemple de graphite, dont les bossages sont réalisés par moulage. Chaque électrode 10, 20 comporte une couche de diffusion 11, 21 (GDL, pour Gas Diffusion Layer, en anglais) disposée au contact de l'une des plaques bipolaires 40, 50 et une couche active 12, 22 située entre la membrane 30 et la couche de diffusion 11, 21. Les couches de 10 diffusion 11, 21 sont réalisées en matériau poreux autorisant la diffusion des espèces réactives à partir du circuit de distribution des plaques bipolaires jusqu'aux couches actives, ainsi que la diffusion des produits issus de la réaction électrochimique jusqu'au même circuit de distribution. Les couches actives 12, 22 sont le lieu des réactions électrochimiques. Elles comportent des matériaux adaptés pour permettre les réactions 15 d'oxydation et de réduction aux interfaces respectives de l'anode et la cathode avec la membrane. Plus précisément, elles comportent chacune un ionomère assurant la conductivité protonique, par exemple du Nafion, un catalyseur pour générer la réaction électrochimique, par exemple du platine, et un support électriquement conducteur, par exemple une matrice carbonée. 20 Comme mentionné précédemment, l'activité catalytique de la couche active dépend notamment de la surface active du catalyseur, c'est-à-dire de la surface du catalyseur susceptible d'être en contact avec l'espèce réactive, exprimée en centimètre carré ou parfois en centimètre carré par unité de masse de catalyseur, ainsi que de la charge en catalyseur (ou taux de chargement ou teneur), exprimée en gramme par unité de surface 25 de la cellule. La figure lb représente schématiquement la distribution spatiale de température T de la cellule électrochimique comme résultant d'une relation de corrélation entre la distribution spatiale d'un terme source de production de chaleur, autrement dit d'une grandeur 3041479 10 représentative de la production de chaleur par la cellule, tel que le flux de chaleur produite Q, et la distribution spatiale d'une grandeur représentative de l'évacuation de la chaleur produite, par exemple le débit massique D du fluide caloporteur dans le circuit de refroidissement. 5 Ainsi, contrairement à l'enseignement du document relatif à l'art antérieur cité précédemment, homogénéiser la distribution de production de chaleur Q et donc celle de l'échauffement de la cellule ne peut suffire à obtenir l'homogénéisation de la distribution de température T effective de la cellule. En effet, il importe de tenir compte à la fois de la présence éventuelle d'inhomogénéités locales dans le terme Q de production de chaleur 10 ainsi que de la présence éventuelle d'inhomogénéités locales dans le terme D d'évacuation de chaleur. La production locale de chaleur, ou plus précisément le flux local de chaleur produite Q, est directement proportionnelle à la production locale d'énergie électrique, ou plus précisément à la densité locale de courant I, comme l'exprime la relation entre les 15 distributions spatiales respectives : Qx,y = lx,y (4H/2F - Ux,y) (2) où AH est l'enthalpie de la réaction électrochimique, F la constante de Faraday, et Ux,y la distribution spatiale de la tension locale de la cellule, l'enthalpie et la tension pouvant être considérées comme quasiment homogènes en tout point de la cellule. Ainsi, la production 20 de chaleur est impactée par toute inhomogénéité qui provient de paramètres fluidiques (dimensionnement des circuits de distribution des espèces réactives...), électrochimiques (propriétés locales des électrodes et de la membrane...), mais également électriques (résistances électriques des différents éléments de la cellule, par exemple résistivités des matériaux et résistances de contact entre les éléments de la cellule...), qui tous influent sur 25 la distribution de densité de courant. L'évacuation de chaleur par l'écoulement du fluide caloporteur peut également présenter des inhomogénéités locales qui proviennent notamment de pertes de charge singulières dans le circuit de refroidissement. Ces pertes de charge sont issues du dimensionnement du circuit de refroidissement lors de la réalisation des plaques bipolaires, et peuvent se 3041479 11 traduire par la formation de zones à fort débit ou à faible débit au sein du circuit de refroidissement. Dans le cadre de l'invention, on cherche à adapter la distribution spatiale d'un paramètre 5 d'intérêt influençant la production d'énergie électrique, et donc d'énergie thermique, pour que la distribution spatiale de la température effective de la cellule en fonctionnement corresponde à celle d'une température de consigne, en tenant compte de la distribution spatiale de l'évacuation effective de chaleur que présente la cellule électrochimique. Le paramètre influençant la production d'énergie électrique est ici un paramètre 10 représentatif de l'activité catalytique de la couche active d'au moins une électrode de la cellule électrochimique. La valeur de ce paramètre d'intérêt influe localement sur la densité de courant lx,y dans la mesure où elle détermine la quantité d'électrons recombinés ou de protons générés respectivement à la cathode et à l'anode lors des réactions électrochimiques de réduction et d'oxydation. Le paramètre d'intérêt est ici la charge Wc 15 du catalyseur, ou charge catalytique. Il peut également s'agir de tout paramètre équivalent, telle que la surface active du catalyseur. On comprend ici que, plus la charge locale Wcx,y en catalyseur est faible, plus faible sera la densité locale lx,y de courant produit par la cellule et donc le flux local Clx,y de chaleur produite. Par température de la cellule, on entend la température locale, c'est-à-dire la distribution 20 spatiale de température de l'un ou l'autre des éléments de la cellule, par exemple de l'une des plaques bipolaires voire de l'une des électrodes. La température de la cellule peut ainsi correspondre à la distribution spatiale de température du fluide caloporteur dans le circuit de refroidissement. La température effective de la cellule est la distribution spatiale de température lors du fonctionnement de la cellule au point de polarisation définie par la 25 tension de la cellule Umt et la densité totale de courant Itot, c'est-à-dire la densité locale de courant lx,y intégrée sur toute la surface de la cellule. Par paramètre représentatif de l'évacuation de chaleur, on entend un paramètre dont la valeur traduit la capacité de la cellule à évacuer localement la chaleur produite. Il peut s'agir 3041479 12 en particulier du débit local, massique ou volumique, du fluide caloporteur circulant dans le circuit de refroidissement. Par distribution spatiale d'un paramètre, on entend la valeur locale de ce paramètre en tout point de la cellule, ou plus précisément, en tout point (x,y) d'un plan parallèle à la 5 cellule, au niveau de la zone dite active correspondant à l'étendue surfacique des couches actives des électrodes. Ainsi, la cellule électrochimique dont la charge en catalyseur serait répartie spatialement selon la distribution ainsi déterminée présenterait une température effective, ou température lors du fonctionnement de la cellule, sensiblement identique à la température 10 de consigne. Cette température de consigne présente avantageusement une distribution spatiale sensiblement homogène, en scalaire ou en gradient. Par homogène en scalaire, on entend que la valeur locale de la température est sensiblement constant. Par homogène en gradient, on entend que le gradient local de température est sensiblement constant. La cellule ne présente alors pas de zones de sur-température, ou points chauds, qui d'une part 15 peuvent augmenter la cinétique des réactions de dégradation des éléments de la cellule et d'autre part peuvent générer des contraintes mécaniques susceptibles de dégrader la tenue mécanique des éléments de la cellule. La durée de vie de la cellule électrochimique est alors préservée. Par point chaud, on entend par exemple une zone de la cellule qui prend un pic de température ou un pic de gradient de température. Plus précisément, un 20 point chaud peut être une zone où l'écart entre la température locale et la température d'entrée du circuit de refroidissement est supérieure au produit d'un coefficient et de l'écart de température entre l'entrée et la sortie du circuit de refroidissement, le coefficient pouvant être de l'ordre de 1.2 à 3 voire davantage, et est de préférence de l'ordre de 1.5. A titre d'exemple, pour une température de 77°C à l'entrée du circuit de refroidissement 25 et de 80°C à la sortie du circuit, et pour un coefficient égal à 1.5, un point chaud est une zone de la cellule où la température locale excède 81.5°C. La figure 2 est un organigramme d'un procédé de détermination de la distribution spatiale du paramètre d'intérêt représentatif de l'activité catalytique d'une électrode, selon un 3041479 13 premier mode de réalisation. Dans cet exemple, le paramètre d'intérêt est la charge Wc en catalyseur d'une couche active d'une électrode de la cellule, par exemple la couche active cathodique, dont la valeur Wcx,y influe directement sur la densité locale de courant électrique lx,y localement produit lors de la réaction électrochimique et donc sur le flux local 5 de chaleur produite Q. D'une manière générale, selon ce premier mode de réalisation, on détermine une distribution spatiale optimisée Wcxf,y de la charge Wc en catalyseur à partir de l'estimation de la distribution spatiale AT'e3, d'un écart 4Te entre une température effective r de la cellule en fonctionnement dans laquelle la charge catalytique est répartie spatialement 10 selon une distribution initiale donnée, et une température de consigne r prédéterminée. On est alors en mesure de modifier localement la charge catalytique Wc de la couche active pour qu'elle présente la distribution optimisée WcxfY' de manière à ce que la température effective r de la cellule modifiée soit sensiblement égale à la température de consigne r. Lors d'une première étape 110, on définit une cellule électrochimique de référence au sein 15 de laquelle la charge catalytique Wc de la couche active cathodique est répartie spatialement selon une distribution initiale Wcxi,y. La cellule présente une structure identique ou similaire à celle décrite en référence à la figure 1. La distribution spatiale initiale Wc'i,3, de la charge catalytique Wc peut être sensiblement homogène, c'est-à-dire ici qu'elle présente une valeur sensiblement constante en tout point de la zone active. 20 Lors d'une étape 120, on définit une distribution spatiale l'xcy d'une température de consigne r de la cellule de référence lorsque celle-ci est en fonctionnement et produit une densité totale de courant Itot pour une tension Utot donnée. Au premier ordre, la température de consigne r de la cellule peut correspondre à une température du fluide caloporteur au sein du circuit de refroidissement, la distribution de cette température 25 dépendant alors notamment des valeurs en entrée Te et en sortie Tsc du circuit de refroidissement. A titre illustratif, la température d'entrée peut être fixée préalablement, par exemple à 75°C, et la température de sortie peut être estimée à partir de la puissance thermique Pth à évacuer, celle-ci correspondant à la puissance électrique Pe = Itot.Utot délivrée lors du fonctionnement de la cellule. La puissance thermique Pth est estimée à 3041479 14 partir de l'intégration sur la zone active du flux local de chaleur produite Clx,y obtenu à partir de la relation (2). L'estimation de la température de sortie Tsc est ensuite obtenue par corrélation de la puissance thermique Pth estimée précédemment avec le débit total moyen du liquide caloporteur dans le circuit de refroidissement, moyennant la capacité 5 thermique cp du fluide caloporteur. On peut ensuite définir la distribution spatiale l'xcy de la température de consigne r à partir des valeurs de température du fluide caloporteur en entrée Te et en sortie Tsc du circuit de refroidissement, la distribution l'xcy étant avantageusement homogène en gradient, c'est-à-dire que le gradient local de température de consigne est sensiblement constant. 10 Lors d'une étape 130, on obtient une distribution spatiale Txr3, d'une première grandeur thermique représentative de la température de la cellule en fonctionnement. La première grandeur thermique est ici la température effective r de la cellule électrochimique lorsqu'elle est en fonctionnement dans les mêmes conditions opératoires que celles considérées dans l'étape 120. Cette distribution Txr3, n'est pas estimée mais est issue d'une 15 mesure par des moyens expérimentaux ou numériques. Elle peut ainsi être obtenue par mesure expérimentale au sein d'une cellule électrochimique présentant les mêmes propriétés que la cellule de référence définie lors de l'étape 110, par exemple au moyen d'une carte S++ commercialisée par la société « S++ Simulation Services » comportant une plaque invasive insérée entre deux plaques bipolaires et adaptée à mesurer une 20 distribution spatiale de température. La distribution Txr3, de température effective peut également être obtenue par simulation numérique à partir d'un modèle électrochimique de cellule, par exemple celui décrit dans la publication d'Inoue et al, Numerical analysis of relative humidity distribution in polymer electrolyte fuel cell stock including cooling water, J. Power Sources 162 (2006) 81-93. 25 La distribution Txr3, de la température effective r obtenue par mesure expérimentale ou numérique prend ainsi en compte les inhomogénéités locales du flux de chaleur produite qui dépend de la densité locale de courant, ainsi que celles de l'évacuation de chaleur qui dépend notamment du débit local du fluide caloporteur dans le circuit de refroidissement. 3041479 15 Lors d'une étape 140, on estime la distribution spatiale d'une seconde grandeur thermique, ici une grandeur 4Te représentative d'un écart local entre la température effective r et la température de consigne r. L'estimation de cette grandeur d'écart local 4Te est calculée à partir de la distribution spatiale l''cy de la température de consigne r définie lors de l'étape 5 120 et de la distribution spatiale Txryde la température effective r mesurée lors de l'étape 130. Il peut s'agir de la différence entre la valeur locale de la température effective et celle de la température de consigne, ou d'un rapport entre ces valeurs, ou autre. On considère ici la différence termes à termes entre les distributions de la température effective et de la température de consigne : Ar,y = T,',:3, - 'I',y. 10 Lors d'une étape 150, on détermine ensuite la distribution spatiale Wc'f,y de la charge catalytique Wc en fonction de la distribution spatiale Ar,y de l'écart local 4Te. Selon une première variante, on calcule tout d'abord un coefficient correcteur dont la distribution spatiale est proportionnelle termes à termes à celle AT'e,3, de l'écart local 4Te. A titre d'exemple, le coefficient correcteur varie continument entre une valeur minimale 15 et une valeur maximale à mesure que l'écart local varie respectivement entre une valeur sensiblement nulle à une valeur maximale. La valeur minimale du coefficient correcteur peut être ici sensiblement comprise entre 0 et 0.75, voire entre 0 et 0.5, et par exemple être sensiblement égale à 0, 0.1, 0.2, 0.3, et la valeur maximale peut être sensiblement égale à l'unité. On peut ensuite déterminer la distribution spatiale Wc'f,y de la charge 20 catalytique Wc par corrélation, par exemple une multiplication termes à termes, de la distribution spatiale initiale Wc'i,3, de la charge catalytique Wc avec la distribution spatiale du coefficient correcteur. Ainsi, dans les zones de la cellule où l'écart 4Te entre la température effective et la température de consigne présente une valeur maximale, c'est-à-dire dans les zones dites de points chauds, la valeur locale initiale Wc'i,3, de la charge 25 catalytique Wc est multipliée par la valeur locale du coefficient correcteur, par exemple égale à 0.25, voire moins. Ainsi, la charge catalytique Wc présente une nouvelle valeur locale diminuée par rapport à la valeur locale initiale, qui réduit la densité locale de courant et donc le flux de chaleur produite. 3041479 16 Selon une deuxième variante, on identifie tout d'abord au moins une zone Z, de la cellule dans laquelle l'écart 4Te présente une valeur supérieure à une valeur seuil prédéterminée, la valeur seuil étant par exemple représentative d'un point chaud. On détermine ensuite la distribution spatiale Wcxf,y de la charge catalytique Wc par modification de la distribution 5 spatiale initiale Wc'i,y dans la zone Z, identifiée en fonction de la valeur locale de l'écart 4Te dans cette zone. A titre d'exemple, la distribution spatiale initiale Wc'i,y peut être localement modifiée à partir d'un coefficient correcteur dont la valeur est proportionnelle à celle de l'écart 4Te dans cette zone. Comme dans la première variante, le coefficient correcteur varie continûment entre une valeur minimale et une valeur maximale, par 10 exemple entre 0 et 1. On obtient ainsi une distribution spatiale Wc'f,y de la charge catalytique Wc de la couche active. On est alors en mesure de modifier la distribution initiale Wc de la charge catalytique Wc de la couche active cathodique de la cellule de référence pour qu'elle présente la nouvelle distribution déterminée lors de l'étape 150. La cellule ainsi optimisée 15 présente alors, en fonctionnement, une température effective dont la distribution spatiale est sensiblement égale à celle de la température de consigne. Dans la mesure où la distribution de la température de consigne est avantageusement homogène, la cellule en fonctionnement présente une température effective dont la distribution est également sensiblement homogène, permettant ainsi de préserver la durée de vie de la cellule. 20 La figure 3 est un organigramme d'un procédé de détermination de la distribution spatiale d'un paramètre d'intérêt représentatif de l'activité catalytique d'une couche active d'au moins une électrode d'une cellule électrochimique, selon un deuxième mode de réalisation. Dans cet exemple, le paramètre d'intérêt est la charge catalytique Wc de la 25 couche active, ici de la cathode, dont la valeur influe directement sur la densité de courant électrique localement échangée lors de la réaction électrochimique. D'une manière générale, selon ce deuxième mode de réalisation, on détermine la distribution spatiale Wc'f,y de la charge catalytique Wc à partir de l'estimation de la 3041479 17 distribution spatiale de la production de chaleur nécessaire pour obtenir la distribution spatiale d'une température de consigne, en tenant compte de la distribution spatiale d'une grandeur thermique représentative de l'évacuation effective de chaleur dans la cellule. On est alors en mesure de modifier la distribution initiale de la charge catalytique Wc de la 5 couche active pour qu'elle présente une distribution spatiale optimisée, de manière à ce que la température effective soit alors sensiblement égale à la température de consigne. Lors d'une première étape 210, on définit une cellule électrochimique de référence au sein de laquelle la charge catalytique Wc de la couche active cathodique est répartie spatialement selon une distribution initiale Wc'i,y. La distribution spatiale initiale Wc'i,3, de 10 la charge catalytique Wc peut être sensiblement homogène en scalaire, de sorte que la valeur locale est sensiblement constante en tout point de la zone active. La cellule présente une structure identique ou similaire à celle décrite en référence à la figure 1 et cette étape est similaire ou identique à l'étape 110 décrite précédemment. Lors d'une étape 220, on définit une distribution spatiale l'xcy d'une température de 15 consigne r de la cellule de référence lorsque celle-ci est en fonctionnement et produit une densité totale de courant Itot pour une tension Utot donnée. Cette étape est similaire ou identique à l'étape 120 décrite précédemment. De manière facultative, il est avantageux de préciser la distribution spatiale T'cy de la température de consigne r en fonction de la distribution spatiale de la concentration en 20 espèce réactive au niveau de la zone active, entre l'entrée et la sortie du circuit de distribution correspondant. En effet, la consommation d'espèce réactive au sein de la zone active de la cellule conduit à une diminution progressive de la concentration d'espèce réactive le long du circuit de distribution. Cette diminution progressive se traduit par une diminution de la densité locale de courant produit par la cellule et donc de la production 25 locale de chaleur, ce qui peut conduire à la formation d'inhomogénéités de la température de la cellule. Pour compenser cette diminution progressive de la production de chaleur, il est avantageux de définir une température de consigne qui tienne compte de la diminution de la concentration d'espèce réactive, de sorte que la température effective de la cellule en fonctionnement correspond à la température de consigne, celle-ci présentant 3041479 18 avantageusement une distribution spatiale homogène. Pour cela, la distribution spatiale 1"'c3, de la température de consigne précisée 7.1"c peut s'écrire, par exemple : TXy= T, Ki.[ max(c4) - ] (3) où ci est la distribution spatiale de la concentration c' en espèce réactive i, par exemple x,y 5 en oxygène, et K' est une constante positive, par exemple proche de 1, qui peut ultérieurement être ajustée. La distribution spatiale cxi,3, de la concentration c' peut être estimée au premier ordre à partir du cheminement des canaux du circuit de distribution de l'espèce réactive considérée et en supposant une consommation homogène sur toute la zone active. Elle peut également être déterminée plus précisément par mesure 10 expérimentale ou numérique de la distribution spatiale de la densité de courant dans une cellule similaire ou identique à la cellule de référence, qui permet de déduire celle de la concentration de l'espèce réactive. D'autres relations (3) peuvent convenir pour préciser la distribution spatiale de la température de consigne en tenant compte de l'évolution spatiale de la concentration en espèce réactive. On obtient ainsi une distribution spatiale 15 1"xc3, de la température de consigne 7.1"c qui permet ainsi d'obtenir une distribution de la température effective de la cellule dont l'homogénéité est améliorée. Par ailleurs, de manière facultative et éventuellement complémentaire de l'étape de précision de la température de consigne décrite précédemment, il est avantageux de préciser la distribution spatiale l'xcy de la température de consigne r en fonction de la 20 distribution spatiale cpx,y de l'humidité relative cp dans les circuits de distribution. L'humidité relative cp est définie classiquement comme le rapport entre la pression partielle PH2O de la vapeur d'eau contenue localement dans le gaz circulant dans le circuit de distribution sur la pression de vapeur saturante Psat. L'humidité relative cp a un effet sur la cinétique des réactions électrochimiques. Aussi, pour compenser la variation locale de l'humidité 25 relative, il est avantageux de définir une température de consigne qui compense cette variation locale, par exemple l'engorgement ou l'assèchement local des circuits de distribution, de sorte que la température effective de la cellule en fonctionnement présente une distribution spatiale homogène. Pour cela, la distribution spatiale 1';',c3, de la température de consigne précisée 7.1"'c peut s'écrire, par exemple : 3041479 19 = T),37 WP.[ (Px,y (pin] (4) où (px,3, est la distribution spatiale de l'humidité relative cf) dans le circuit de distribution, (pin est l'humidité relative en entrée du circuit de distribution, et 1(q) est une constante positive, par exemple proche de 1, qui peut ultérieurement être ajustée. La distribution 5 (px,3, de l'humidité relative cf) d'eau total peut être estimée au premier ordre à partir du cheminement des canaux du circuit de distribution considéré et en supposant une densité de courant homogène sur toute la zone active. Elle peut également être déterminée plus précisément par mesure expérimentale ou numérique de la distribution spatiale de la densité de courant dans une cellule similaire ou identique à la cellule de référence, qui 10 permet de déduire celle de l'humidité relative. D'autres relations (4) peuvent convenir pour préciser la distribution spatiale de la température de consigne à partir de l'évolution spatiale de l'humidité relative. On obtient ainsi une distribution spatiale r-Qy de la température de consigne 7.1"c qui permet ainsi d'obtenir une distribution de la température effective de la cellule dont l'homogénéité est améliorée. 15 Lors d'une étape 230, on obtient une distribution spatiale Dxr,3, d'une première grandeur thermique représentative de l'évacuation de la chaleur Dr au sein de la cellule en fonctionnement. La première grandeur thermique est ici le débit massique Dr de fluide caloporteur dans le circuit de refroidissement. Cette distribution Dxr,3, n'est pas estimée mais est issue d'une mesure par des moyens expérimentaux ou numériques. Elle peut ainsi 20 être obtenue par mesure expérimentale au sein d'une cellule électrochimique présentant les mêmes propriétés que la cellule de référence définie lors de l'étape 210, par exemple au moyen d'une technique de type PIV (Particle Image Velocimetry, en anglais) ou de toute autre technique adaptée, réalisée dans un circuit de refroidissement présentant les mêmes caractéristiques dimensionnelles que celui de la cellule de référence. La distribution Dxr,3, 25 de débit massique Dr peut également être obtenue par simulation numérique à partir d'un logiciel de simulation d'écoulement tels que, par exemple, FLUENT ou COMSOL. Lors d'une étape 240, on estime la distribution spatiale Qex,y d'une seconde grandeur thermique Qe à partir de ladite distribution spatiale l'xcy de la température de consigne T` 3041479 20 définie lors de l'étape 220 et de ladite distribution spatiale K,y du débit de fluide caloporteur D' obtenue lors de l'étape 230. La seconde grandeur thermique est représentative de la production locale de chaleur et correspond ici au flux local de chaleur Ce qu'évacue le fluide caloporteur D' pour obtenir la température de consigne r. 5 Pour cela, comme l'illustre la figure 4, on discrétise le circuit de refroidissement en un maillage dont chaque maille est un volume élémentaire (i,j) traversé par le fluide caloporteur. Ainsi, chaque maille (i,j) du circuit de distribution présente deux grandeurs connues : la température locale Tici de consigne et le débit local Dfi de fluide caloporteur ; et une grandeur à déterminer : le flux local de chaleur produite Qr. On calcule ensuite les 10 transferts thermiques et fluidiques entre la maille considérée et les mailles adjacentes en déterminant d'une part les différences de température, et d'autre part les débits du fluide caloporteur au niveau des quatre facettes de la maille considérée. Les différences de température au niveau des quatre facettes de la maille (i,j) peuvent être calculées de la manière suivante : 15 (Vil" = Tij - Tici+1 (5-1) 8Ti2i = Tl - (5-2) 8Ti3i = - Tic+1,1 (5-3) 8Tii = - (5-4) Les débits du fluide caloporteur au niveau des quatre facettes de la maille (i,j) peuvent être 20 obtenues par projection du débit massique Dfi (ici une donnée vectorielle) sur les vecteurs ex et ey passant respectivement par les mailles (i-1,j), (i,j) et (i+1,j), et par les mailles (i,j-1), (i,j) et (i,j+1) : dit = (Di5. ey + D +1. ey) 12 (6-1) di2i = (Di5. ex + Di 1,1.ex) 12 (6-2) 25 di3i = (Di5. ex + Dr+1 J. e ,) 12 (6-3) d = (Di5. ey + ey) 12 (6-4) Par convention, on considère que le débit local d,,, est positif lorsque le fluide entre dans la maille (i,j) et négatif lorsqu'il en sort. 3041479 21 On calcule enfin la distribution spatiale (2,e,), du flux Qe de chaleur produit par la cellule à partir de la relation : (2c,y (2,; - 8T15 (7) On obtient ainsi la distribution spatiale du flux de chaleur Qe que doit produire la cellule 5 pour que la distribution de température effective corresponde à celle de la température de consigne, connaissant la distribution du débit massique effectif du fluide caloporteur dans le circuit de distribution. Lors d'une étape 250, on détermine la distribution spatiale Wcxf,3, de la charge catalytique Wc en fonction de la distribution spatiale 1:2'e,y du flux de chaleur produite Qe. Pour cela, on 10 peut estimer tout d'abord la distribution spatiale de la densité d'un signal électrique produit par la cellule en fonctionnement, par exemple une densité de courant le, à partir de la distribution spatiale 1:2'e,y estimée du flux de chaleur produite Qe. Dans la mesure où le flux de chaleur produite Qe est approximativement proportionnel à la densité de courant le, celle-ci peut être déterminée à partir de la relation : 15 le e Itot (8) = Qtot où Itot est la densité de courant totale délivrée par la cellule électrochimique en fonctionnement, et (lot est le flux de chaleur produite totale obtenue par intégrale de la distribution spatiale Qex,3, sur toute la surface active. Ensuite, on détermine la nouvelle distribution spatiale Wcxf,3, de la charge catalytique Wc à 20 partir de la densité locale du courant électrique Ixe,y. Pour cela, une approche consiste à déterminer les valeurs minimale Wcmin et maximale Wcmax de la charge Wc de la couche active. Il peut s'agir de valeurs couramment admises dans le cadre des piles de type PEM. On calcule ensuite la distribution spatiale Wcxf,3, par exemple à partir de la relation : le (9) WCf y = max (WCmax, WCmin ee lx,y 25 où Ime ax est la valeur maximale de la densité locale de courant Ixe,y. La charge Wc locale varie ainsi linéairement entre les valeurs minimale Wcn-fin et maximale Wcmax en fonction de la 3041479 22 valeur locale de la densité de courant le. Bien entendu, toute autre loi, par exemple polynomiale, exponentielle ou logarithmique, faisant varier la charge Wc locale de sorte que la valeur maximale Wcmax corresponde à une densité locale de courant maximale et inversement, peut convenir. 5 De manière alternative à la relation (9), on détermine la distribution spatiale Wcxf,3, à partir d'une courbe donnant la relation entre la densité de courant I et la charge Wc. La figure 5 est un exemple d'une telle relation, qui relie ici de manière équivalente la densité de puissance d'une pile à combustible en fonctionnement à la charge catalytique en platine des cellules électrochimique, pour un rendement de 50%. En fonction de la distribution 10 spatiale Ixe,y, on détermine la distribution spatiale Wcxf,y. On remarque ici que la charge catalytique peut varier d'une valeur minimale de l'ordre de 0,05mg/cm2 à une valeur maximale de l'ordre de 0,4mg/cm2 à 0,5mg/cm2, ce qui entraîne une variation de la densité de courant et donc de puissance allant de 100mW/cm2 à 400mW/cm2. On a ainsi déterminé la distribution spatiale Wcxf,3, de la charge Wc qui tienne compte de la 15 distribution de production d'énergie électrique le et donc thermique Qe qu'il faudrait pour que la température effective de la cellule en fonctionnement corresponde à la température de consigne Tc, en considérant l'évacuation effective Dr de la chaleur par le circuit de refroidissement. Dans la mesure où la température de consigne est avantageusement homogène spatialement, une cellule dont au moins l'une des électrodes comporte une 20 couche active dont la charge catalytique Wc est répartie suivant la distribution spatiale ainsi déterminé Wcfx, présente, lorsqu'elle est en fonctionnement au point de polarisation U tot et Itot, une température effective dont la distribution spatiale est homogène. La figure 6a illustre un exemple de distribution spatiale de température effective d'une cellule électrochimique en fonctionnement, dans le cas où les couches actives anodique et 25 cathodique présentent une charge catalytique sensiblement uniforme. Trois zones de sur- température, ou zones de points chauds, sont présentes, référenciées Zi à Z3. Les conditions opératoires dans ce cas sont une pression de 1,5bar, une température de 80°C en sortie du circuit de refroidissement, une humidité relative en entrée de 50% à l'anode et à la cathode, une stoechiométrie de 1,8 pour l'air et de 1,5 pour l'hydrogène. 3041479 23 La figure 6b illustre la distribution spatiale de température effective de la même cellule électrochimique en fonctionnement, dans le cas où au moins une couche active présente une charge catalytique présentant la distribution spatiale modifiée selon le premier mode de réalisation du procédé (figure 2) uniquement dans une partie de la cellule indiquée par 5 les traits pointillés. On peut remarquer que, dans cette partie où un point chaud Z2 était présent, la distribution de température est sensiblement homogène et ne présente pas de zones de sur-température. Dans les autres parties de la cellule où la charge catalytique n'a pas été modifiée, les points chauds Zi et Z3 restent présents. 10 Un procédé de réalisation d'une électrode de cellule électrochimique est maintenant décrit, ici la cathode de la cellule. On considère une cellule électrochimique identique ou similaire à la cellule de référence définie lors des étapes 110 et 210. Elle comporte deux électrodes séparées l'une de l'autre par un électrolyte et disposées entre deux plaques bipolaires adaptées à amener des espèces réactives aux électrodes et à évacuer de la 15 chaleur produite par la cellule en fonctionnement. Les électrodes comportent chacune une couche de diffusion et une couche active. La cathode présente ici une couche active dont la charge catalytique Wc est répartie spatialement selon une distribution initiale Wcxi,y. On détermine, par le procédé décrit précédemment en référence aux figures 2 ou 3, une distribution spatiale Wc'f,y de la charge catalytique de la couche active cathodique. On 20 réalise ensuite, lors d'une étape 160 (figure 2) ou 260 (figure 3), la cathode de sorte que la charge catalytique Wc de la couche active présente la distribution spatiale Wc'f,y déterminée. Ainsi, on modifie localement la charge catalytique Wc de la couche active de manière à former des inhomogénéités dans la distribution spatiale Wcx,y en vue de compenser les 25 inhomogénéités de la température effective de la cellule en fonctionnement vis-à-vis d'une température de consigne. A titre d'exemple, on diminue localement la charge Wc de la couche active dans des zones où la température effective de la cellule est susceptible d'être supérieure à la température de consigne. Ainsi, la densité locale de courant est diminuée, 3041479 24 ce qui diminue en conséquence le flux local de chaleur produite. Ainsi, la température effective de la cellule correspond localement à la température de consigne. Pour cela, on dépose sur un support de la couche active, par exemple carboné, une encre catalytique comportant, à titre illustratif, un catalyseur, par exemple 4g de Pt3Co, un 5 ionomère, par exemple 6,4g de Nafion, de l'eau, par exemple à hauteur de 18,5g, et un solvant, par exemple 1,5g d'éthanol. L'encre catalytique peut être déposée à la surface du support de la couche active de la cathode, dans au moins une zone Z, de sur-température identifiée, voire sur toute la surface de la couche active, par toutes les techniques de dépôts connues de l'homme du métier, 10 par exemple par sérigraphie. En référence aux figures 7a à 7c, dans le cas où un dépôt est réalisé par sérigraphie dans certaines zones de la couche active, un écran de sérigraphie 60 est réalisé qui forme un maillage d'ouvertures traversantes 61, par exemple de quelques centimètres carrés de surface. Certaines ouvertures 62 sont obturées pour empêcher le dépôt de l'encre 15 catalytique dans les zones Z, identifiées (figure 7a). On dépose ensuite l'écran 60 sur le support de la couche active, puis on dépose l'encre catalytique 63 sur la couche active au travers de l'écran 60 (figure 7b). Après retrait de l'écran, la couche active comporte l'encre catalytique 63 sur sensiblement toute sa surface hormis les zones Z, (figure 7c). La quantité de catalyseur peut être ajustée dans les zones Z, en effectuant éventuellement un ou 20 plusieurs dépôts successifs dans ces zones avec une encre catalytique présentant une quantité de catalyseur différente de celle utilisée pour les zones adjacentes aux zones Z. Ainsi, la couche active comporte une charge catalytique dont la valeur varie localement. Dans au moins une zone Z,, la valeur locale de la charge catalytique est inférieure à la valeur moyenne de la charge hors de cette zone Z. 25 A titre illustratif, les zones Z,, dans lesquelles la charge catalytique Wc est modifiée selon la distribution spatiale Wcxf,y déterminée par le procédé de détermination décrit en référence aux figures 2 ou 3, peuvent être de l'ordre de quelques centaines de microns carrés à quelques centimètres carrés, par exemple de 400urn2 qui est la taille typique d'une maille d'un écran de sérigraphie classiquement utilisé, à 110cm2 qui est ici la moitié de la surface 3041479 25 totale d'une cellule de 220cm2. Bien entendu, ces données sont ici purement illustratives. Plus précisément, une telle zone Z, présente une surface allant de 0,5cm2 à 40cm2, où 0,5cm2 est ici la taille minimale typique d'une maille d'un dispositif de mesure permettant de mesurer les inhomogénéités de température, et où 40cm2 est une taille typique des 5 points chauds mesurés. La surface des zones Z, peut ainsi être de l'ordre de 6cm2. Il a été observé qu'une diminution de 100% de catalyseur sur une bande de 2x16cm2 permet de diminuer la densité locale de courant de 0,3A/cm2, ce qui entraîne une baisse de la température effective locale de 2°C en moyenne. Ainsi, par extrapolation, on peut supposer qu'une diminution de 0,1mg/cm2 de platine sur au moins 1cm2 entraîne une 10 baisse de la température effective dans cette zone de 0,5°C. La charge catalytique peut être comprise entre 0 et 0,5mg/cm2 dans le cas du platine pour une cellule électrochimique de pile à combustible. La charge catalytique présente une valeur locale maximale hors des zones Z, et une valeur locale dans les zones Z, inférieure à la valeur nominale. Une valeur moyenne nominale de la charge est comprise entre les 15 valeurs locales maximale et minimale, et correspond à la moyenne de la charge locale sur toute la surface de la couche active. La valeur moyenne nominale dépend notamment de la densité totale de puissance souhaitée. Elle peut ainsi être de l'ordre de 0,4mg/cm2 de platine. La teneur en ionomère, par exemple en Nafion, peut être comprise entre 2% et 30% de la charge catalytique, de préférence comprise entre 7% et 20%, par exemple de 20 l'ordre de 15%. Il peut être avantageux que la distribution spatiale Wcxf,3, déterminée de la charge Wc présente une valeur moyenne nominale sensiblement identique à celle de la distribution initiale Wc'i,y, de manière à ne pas modifier la densité totale de puissance de la cellule en fonctionnement. Ainsi, alors que, dans les zones Z, de points chauds, la distribution spatiale 25 x Wcfx, présente une valeur locale inférieure à celle de la distribution initiale Wci y, on augmente la valeur locale Wcxf,3, hors des zones Z. A titre d'exemple, pour maintenir une valeur moyenne de 0,4mg/cm2 de platine sur l'ensemble de la couche active, une diminution de 0,1mg/cm2 sur une zone Z, de 6cm2 conduit à répartir 0,6mg sur le reste de 3041479 26 la surface de la couche active, par exemple sur 214cm2 pour une cellule de 220cm2 de couche active. Des modes de réalisation particuliers viennent d'être décrits. Différentes variantes et modifications apparaîtront à l'homme du métier. 5 REVENDICATIONS1. Procédé de détermination d'une distribution spatiale (Wc'f,y) d'un paramètre d'intérêt (Wc) représentatif d'une activité catalytique d'une couche active d'au moins une électrode d'une cellule électrochimique, ladite cellule comportant deux électrodes séparées l'une de l'autre par un électrolyte et disposées entre deux plaques bipolaires adaptées à amener des espèces réactives aux électrodes et à évacuer de la chaleur produite par la cellule en fonctionnement, au moins une électrode comportant une couche active comportant au moins un catalyseur, la couche active étant adaptée à générer une réaction électrochimique à partir de l'une des espèces réactives, ledit procédé comprenant les étapes dans lesquelles : i) on définit (110; 210) une cellule électrochimique au sein de laquelle le paramètre d'intérêt (Wc) est réparti spatialement selon une distribution initiale (Wcxi,y) ; ii) on définit (120 ; 220) une distribution spatiale (1',y) d'une température de consigne (r) au sein de la cellule en fonctionnement ; iii) on obtient (130 ; 230) une distribution spatiale (ny ; Dxr,y) d'une première grandeur thermique représentative de la température (r) ou de l'évacuation de la chaleur (Dr) au sein de la cellule en fonctionnement ; iv) on estime (140; 240) la distribution spatiale (Ar,3, ; 1:gc,y) d'une seconde grandeur thermique (Ar ; Qe) à partir de ladite distribution spatiale (ny) de la température de consigne (r) et de ladite distribution spatiale (rcly ; K,y) de la première grandeur thermique (Tr; Dr) ; v) on détermine (150 ; 250) une distribution spatiale (Wcxf,y) du paramètre d'intérêt (Wc) en fonction de la distribution spatiale (Ar,3, ; 1:gc,y) de la seconde grandeur thermique (Ar ; Qe). 2. Procédé de détermination selon la revendication 1, dans lequel le paramètre d'intérêt est la charge en catalyseur (Wc) de la couche active. 3041479 28 3. Procédé de détermination selon la revendication 1 ou 2, dans lequel la distribution spatiale (Wc'f,y) du paramètre d'intérêt (Wc) est déterminée de sorte qu'elle présente une valeur moyenne sensiblement identique à une valeur moyenne de la distribution initiale (Wc4). 5 4. Procédé de détermination selon l'une quelconque des revendications 1 à 3, dans lequel la première grandeur thermique est une température locale effective (r) mesurée au sein de la cellule en fonctionnement, et la seconde grandeur thermique est une grandeur (4Te) représentative d'un écart local entre la température effective (r) et la température de consigne (r). 10 5. Procédé de détermination selon la revendication 4, dans lequel l'étape v) comporte : a) une sous-étape de calcul d'une distribution spatiale d'un coefficient correcteur à partir de la distribution spatiale (AT'e,y) de la seconde grandeur thermique (4Te) ; b) une sous-étape dans laquelle on détermine la distribution spatiale (Wc'f,y) du paramètre d'intérêt (Wc) par corrélation de la distribution spatiale initiale (Wcxi,y) 15 du paramètre d'intérêt (Wc) avec la distribution spatiale du coefficient correcteur. 6. Procédé de détermination selon la revendication 4, dans lequel l'étape v) comporte : a) une sous-étape dans laquelle on identifie au moins une zone Z, de la cellule où la seconde grandeur thermique (4Te) présente une valeur locale estimée supérieure à une valeur seuil prédéterminée ; 20 b) une sous-étape dans laquelle on détermine la distribution spatiale (Wc'f,y) du paramètre d'intérêt (Wc) par modification de la distribution spatiale initiale (Wcxi,y) dans la zone Z, identifiée, en fonction de la valeur locale estimée de la seconde grandeur thermique (4Te) dans cette zone. 7. Procédé de détermination selon l'une quelconque des revendications 1 à 3, dans 25 lequel la première grandeur thermique est représentative de l'évacuation locale (Dr) de la chaleur produite par la cellule en fonctionnement, et la seconde grandeur thermique est représentative de la production locale de chaleur (Qe) par la cellule en fonctionnement. 3041479 29 8. Procédé de détermination selon la revendication 7, dans lequel la première grandeur thermique est le débit local effectif (Dr) mesuré d'un fluide caloporteur circulant dans un circuit de refroidissement d'une plaque bipolaire de la cellule, et la seconde grandeur thermique est le flux local de chaleur (Qe) produite par la cellule en fonctionnement. 5 9. Procédé de détermination selon la revendication 8, dans lequel l'étape v) comporte : a) une première sous-étape dans laquelle on estime une distribution spatiale de la densité d'un signal électrique (le) produit par la cellule en fonctionnement, à partir de ladite distribution spatiale estimée du flux de chaleur produite (Qe), et b) une deuxième sous-étape dans laquelle on détermine la distribution spatiale du 10 paramètre d'intérêt (Wc) à partir de ladite distribution spatiale de la densité du signal électrique (le). 10. Procédé de réalisation d'une électrode de cellule électrochimique, comportant les étapes dans lesquelles : i) on considère une cellule électrochimique de référence comportant deux électrodes 15 séparées l'une de l'autre par un électrolyte et disposées entre deux plaques bipolaires adaptées à amener des espèces réactives aux électrodes et à évacuer de la chaleur produite par la cellule en fonctionnement, au moins une électrode comportant une couche active adaptée à générer une réaction électrochimique à partir de l'une des espèces réactives, la couche active présentant un paramètre 20 d'intérêt représentatif d'une activité catalytique (Wc) répartie spatialement selon une distribution initiale (Wci ) - x,y / ii) on détermine, par le procédé (100 ; 200) selon l'une quelconque des revendications précédentes, une distribution spatiale (Wcxf,y) du paramètre représentatif de l'activité catalytique (Wc) ; 25 iii) on réalise (160 ; 260) ladite électrode de sorte que le paramètre représentatif de l'activité catalytique (Wc) de la couche active de celle-ci présente la distribution spatiale (Wcxf,y) déterminée. 3041479 11. Procédé de réalisation selon la revendication 10, dans lequel l'étape de réalisation de ladite électrode comporte le dépôt, dans au moins une zone (Z,) de la couche active préalablement identifiée à partir de ladite distribution spatiale déterminée (Wfxf,y) du paramètre représentatif de l'activité catalytique (Wc), d'une encre catalytique présentant 5 une charge catalytique inférieure à une charge catalytique moyenne hors de ladite zone (Z,). 12. Procédé de réalisation selon la revendication 11, dans lequel l'encre catalytique est déposée par sérigraphie au moyen d'un écran de sérigraphie comportant au moins une ouverture obturée destinée à recouvrir ladite zone (Z,) de la couche active. 10 13. Electrode de cellule électrochimique comportant une couche active adaptée à générer une réaction électrochimique au sein de la cellule, caractérisée en ce qu'elle présente un paramètre représentatif d'une activité catalytique (Wc) de la couche active, de préférence une charge catalytique, qui présente, dans au moins une zone (Z,) de la couche active, une valeur locale inférieure à une valeur moyenne hors de ladite zone (Z,). 15