Brevet FR3041481A1 — publié le 20170324

DETERMINATION D'UNE DISTRIBUTION SPATIALE D'UN PARAMETRE DE PRODUCTION ELECTRIQUE D'UNE 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 la production d'énergie électrique de la cellule électrochimique ainsi que celui des procédés de réalisation d'une plaque bipolaire 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. 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 3041481 2 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 distribution prend souvent la 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 à 5 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 lors de la réaction par la cellule. Le document FR2976732 illustre un exemple de cellule électrochimique réalisée de 10 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 des espèces réactives. En effet, considérant 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 à 15 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é é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 20 partielle d'oxygène. Cependant, il reste à améliorer encore l'homogénéisation de la distribution spatiale de la température effective de la cellule électrochimique, de manière à préserver la durée de vie de la cellule en limitant la cinétique des réactions de dégradations des différents éléments de la cellules et en réduisant les contraintes mécaniques d'origine thermique 25 susceptibles de fragiliser la tenue mécanique des éléments de la cellule. 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 3041481 3 distribution spatiale d'un paramètre représentatif d'une production locale d'énergie électrique d'une 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. 5 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 de la production d'énergie électrique 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 10 fonctionnement, 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 cellule en fonctionnement ; 15 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 ladite distribution spatiale de la température de consigne et de ladite distribution 20 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 choisi parmi un paramètre représentatif de la résistance électrique de la cellule électrochimique ; un paramètre représentatif d'une 25 résistance de contact entre au moins l'une des électrodes et la plaque bipolaire adjacente ; une charge d'un matériau catalyseur présent au moins dans l'une des électrodes ; et un paramètre représentatif d'une perméabilité d'au moins une électrode. De préférence, le paramètre d'intérêt est un paramètre représentatif de la résistance électrique de la cellule électrochimique, la première grandeur thermique est une 3041481 4 température locale effective mesurée au sein de la cellule en fonctionnement, et la seconde grandeur thermique est 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 : 5 une sous-étape de calcul d'une distribution spatiale d'un coefficient correcteur à 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. 10 Alternativement, l'étape v) peut comporter : une sous-étape dans laquelle on identifie au moins une zone 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 15 d'intérêt par modification de la distribution spatiale initiale dans la zone identifiée, en fonction de la valeur locale estimée de la seconde grandeur thermique dans cette zone. Par ailleurs, le paramètre d'intérêt peut être un paramètre représentatif de la résistance électrique de la cellule électrochimique, la première grandeur thermique est alors 20 représentative de l'évacuation locale de la chaleur produite par la cellule en fonctionnement, et la seconde grandeur thermique est alors représentative de la production locale de chaleur par la cellule en fonctionnement. La première grandeur thermique peut être le débit local effectif mesuré d'un fluide caloporteur circulant dans un circuit de refroidissement d'une plaque bipolaire de la 25 cellule, et la seconde grandeur thermique est le flux local de chaleur produite par la cellule en fonctionnement. L'étape v) peut comporter : 3041481 5 une première sous-étape dans laquelle on estime la distribution spatiale de la densité d'un signal électrique produit par la cellule en fonctionnement, à partir de la distribution spatiale estimée du flux de chaleur produite, et une deuxième sous-étape dans laquelle on détermine la distribution spatiale du 5 paramètre d'intérêt à partir de la densité locale du signal électrique. L'invention porte également sur un procédé de réalisation d'une plaque bipolaire de cellule é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 10 plaques bipolaires adaptées à amener des espèces réactives aux électrodes et à évacuer de la chaleur produite par la cellule en fonctionnement, la cellule présentant une résistance électrique 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 de la résistance électrique ; on réalise ladite plaque bipolaire de sorte que la résistance électrique présente la distribution spatiale déterminée. Ladite plaque bipolaire étant formée de deux tôles embouties assemblées l'une à l'autre par une pluralité de points de soudure répartis selon une distribution spatiale initiale, 20 l'étape de réalisation de la plaque bipolaire peut comporter la modification de la distribution spatiale d'un paramètre représentatif d'une résistance de contact entre les deux tôles en fonction de la distribution spatiale déterminée. La modification de la distribution spatiale du paramètre représentatif de la résistance de contact peut comporter une modification de la répartition de points de soudure assurant 25 l'assemblage des tôles l'une à l'autre en fonction de ladite distribution spatiale déterminée de la résistance électrique. La modification de la distribution spatiale du paramètre représentatif de la résistance de contact peut comporter une disposition d'un revêtement présentant une conductivité 3041481 6 électrique différente de celles des tôles dans des zones préalablement identifiées à partir de ladite distribution spatiale déterminée de la résistance électrique. BRÈVE DESCRIPTION DES DESSINS 5 D'autres aspects, buts, avantages et caractéristiques de l'invention apparaîtront 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 lesquels : la figure la est une représentation schématique en coupe d'un exemple de cellule 10 é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 cellule électrochimique en fonctionnement ; la figure 2 est un organigramme d'un procédé de détermination de la distribution spatiale 15 de la résistance électrique de la cellule électrochimique selon un premier mode de réalisation ; la figure 3 illustre des exemples de courbes de polarisation qui montrent l'évolution de la tension électrique de la cellule en fonction de la densité locale du courant électrique, pour différentes valeurs de la résistance électrique locale ; 20 la figure 4 est un organigramme d'un procédé de détermination de la distribution spatiale de la résistance électrique de la cellule électrochimique selon un deuxième mode de réalisation ; la figure 5 est un exemple de maillage du circuit de refroidissement où chaque maille comporte un terme de production locale de chaleur Qri, un terme d'évacuation local de 25 chaleur et une température Tic,i ; la figure 6a est une représentation schématique en perspective d'une partie de plaque bipolaire formée de deux tôles assemblées l'une à l'autre par un point de soudure ; la 3041481 7 figure 6b est une variante de la figure 6a où un produit de revêtement est situé entre les deux tôles ; et la figure 6c est une représentation schématique en vue de dessus d'un circuit de refroidissement où des zones localement grisées illustrent des points de soudure ou la présence du produit de revêtement. 5 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. De plus, les différents éléments ne sont pas représentés à l'échelle, de manière à privilégier la clarté des figures. 10 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 dont la cathode est alimentée en oxygène et l'anode en hydrogène. L'invention s'applique d'une manière générale à tout type de pile à combustible, en particulier à celles fonctionnant à basse température, c'est-à-dire à une température inférieure à 250°C, ainsi qu'aux électrolyseurs électrochimiques. 15 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. 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 20 disposées entre deux plaques bipolaires 40, 50 adaptées à amener des espèces réactives 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 à 25 l'autre 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 3041481 8 technologies 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 5 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 diffusion 11, 21 sont réalisées en matériau poreux permettant 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 10 jusqu'au même circuit de distribution. Les couches actives 12, 22 comportent des matériaux adaptés à permettre les réactions d'oxydation et de réduction aux interfaces respectives de l'anode et la cathode avec la membrane, par exemple un ionomère et un catalyseur. La figure lb représente schématiquement la distribution spatiale de température T de la 15 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 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 20 caloporteur dans le circuit de refroidissement. 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 l'échauffement de la cellule ne peut suffire à obtenir l'homogénéisation de la distribution de température T de la cellule. En effet, il importe de tenir compte à la fois de la présence 25 ou non d'inhomogénéités locales dans le terme Q de production de chaleur ainsi que de la présence ou non 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 3041481 9 précisément à la densité locale de courant I, comme l'exprime la relation entre les distributions spatiales respectives : Qx,y = lx,y (4H/2F - Ux,y) (1) où AH est l'enthalpie de la réaction électrochimique, F la constante de Faraday, et Ux,y la 5 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 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 10 é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 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 15 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 traduire par la formation de zones à fort débit ou à faible débit au sein du circuit de refroidissement. 20 Dans le cadre de l'invention, on cherche à adapter la distribution spatiale d'un paramètre d'intérêt représentatif de la production d'énergie électrique, et donc d'énergie thermique, pour que la distribution spatiale de la température de la cellule en fonctionnement corresponde à celle d'une température de consigne, en tenant compte de la distribution spatiale de l'évacuation de chaleur que présente la cellule 25 électrochimique. Ainsi, la cellule électrochimique dont le paramètre d'intérêt serait réparti 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 3041481 10 température 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 constante. Par homogène en gradient, on entend que le gradient local de température est 5 sensiblement constant. La cellule ne présente alors pas de zones de sur-température, ou points chauds, qui d'une part 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 10 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 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 15 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 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. Par distribution spatiale d'un paramètre, on entend la valeur locale de ce paramètre en 20 tout point de la cellule, ou plus précisément, en tout point (x,y) d'un plan parallèle à la cellule, au niveau de la zone dite active correspondant à l'étendue surfacique des couches actives des électrodes. Par température de la cellule, on entend la température de l'un ou l'autre des éléments de la cellule, par exemple la température de l'une des plaques bipolaires voire de l'une 25 des électrodes. La température de la cellule peut ainsi correspondre localement à la température du fluide caloporteur dans le circuit de refroidissement. Par paramètre représentatif de la production d'énergie électrique, on entend un paramètre dont la valeur influe localement sur la densité de courant lx,y. Il peut s'agir de la résistance électrique R de la cellule qui dépend notamment de la résistivité des 3041481 11 différents éléments de la cellule (plaques bipolaires, couches de diffusion, couches actives, membrane) et de la résistance électrique de contact entre chacun de ces éléments, en particulier entre les deux tôles de la plaque bipolaire et entre la plaque bipolaire et l'électrode. Il peut également s'agir de la charge Cx,y (ou taux de charge, ou 5 masse surfacique) en catalyseur de la couche active dans la mesure où elle influe directement sur la densité locale de courant lx,y, ou encore de la perméabilité k de la couche de diffusion des électrodes qui détermine la quantité d'espèce réactive susceptible de diffuser localement jusqu'à la couche active. Enfin, par paramètre représentatif de l'évacuation de chaleur, on entend un paramètre 10 dont la valeur traduit la capacité de la cellule à évacuer localement la chaleur produite. Il peut s'agir en particulier du débit local, massique ou volumique, du fluide caloporteur circulant dans le circuit de refroidissement. La figure 2 est un organigramme d'un procédé de détermination de la distribution 15 spatiale d'un paramètre d'intérêt représentatif de la production d'énergie électrique, selon un premier mode de réalisation. Dans cet exemple, le paramètre d'intérêt est représentatif de la résistance électrique R de la cellule électrochimique, et correspond ici à la résistance électrique R, dont la valeur influe directement sur la densité locale de courant électrique lx,y localement produit lors de la réaction électrochimique et donc sur 20 le flux local 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 Rfx,3, de la résistance électrique R à partir de l'estimation de la distribution spatiale AT'e3, d'un écart ATe entre une température effective Tr de la cellule en fonctionnement pour une résistance électrique répartie selon une distribution 25 initiale donnée, et une température de consigne r prédéterminée. On est alors en mesure de modifier la résistance électrique initiale Rix,3, de la cellule pour qu'elle présente la distribution optimisée Rf',y, de manière à ce que la température effective Tr de la cellule modifiée soit sensiblement égale à la température de consigne Tc. 3041481 12 La résistance électrique R de la cellule correspond ici à la somme des résistances des différents éléments de la cellule, à savoir principalement la résistance RGDL,a et RGDLc des couches de diffusion anodique et cathodique, la résistance Rn, de la membrane, la résistance RPB d'une plaque bipolaire, ainsi que la résistance de contact RCGDUBp entre la 5 couche de diffusion et la plaque bipolaire et la résistance de contact RcPB entre les deux tôles métalliques de la plaque bipolaire. Elle peut être préalablement estimée ou mesurée localement sur des échantillons de cellule. Lors d'une première étape 110, on définit une cellule électrochimique de référence au sein de laquelle la résistance électrique R de la cellule est répartie spatialement selon une 10 distribution initiale Rix,y. La cellule présente une structure identique ou similaire à celle décrite en référence à la figure 1. La distribution spatiale initiale Rix,3, de la résistance électrique peut être sensiblement homogène, c'est-à-dire ici qu'elle présente une valeur sensiblement constante en tout point de la zone active. A titre d'exemple, on considère que les résistances des différents éléments de la cellule sont sensiblement constantes en 15 tout point de la zone active, et que les points de soudure susceptibles de modifier localement la résistance de contact des plaques bipolaires sont répartis de manière sensiblement homogène au sein de la cellule. 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 20 une densité totale de courant Itot pour une tension Utnt 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 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, 25 par exemple à 75°C, et la température de sortie peut être estimée à partir de la puissance thermique Pth à évacuer qui peut être évaluée par la relation Pth = Itot (4H/2F - Ut0t). 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é thermique 3041481 13 cp du fluide caloporteur. On peut ensuite définir la distribution spatiale l''cy 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 5 température de consigne est sensiblement constant. 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 10 considérées dans l'étape 120. Cette distribution Txr3, n'est pas estimée mais est issue d'une 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 » 15 comportant une plaque invasive insérée entre deux plaques bipolaires et adaptée à mesurer une distribution spatiale de température. La distribution l'xily 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 20 including cooling water, J. Power Sources 162 (2006) 81-93. 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 25 refroidissement. Lors d'une étape 140, on estime la distribution spatiale d'une seconde grandeur thermique, ici une grandeur ATe 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 ATe est calculée à partir de la distribution spatiale T'cy de la température de consigne r 3041481 14 définie lors de l'étape 120 et de la distribution spatiale T,',:yde 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 5 la température effective et de la température de consigne : Ar,y = -TXy. Lors d'une étape 150, on détermine ensuite la distribution spatiale Rf',3, du paramètre d'intérêt R 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. 10 A titre d'exemple, le coefficient correcteur varie continument entre une valeur minimale, par exemple inférieure ou égale à l'unité, jusqu'à une valeur maximale, à mesure que l'écart local varie respectivement entre une valeur sensiblement nulle à une valeur maximale. En référence à la figure 3, la valeur maximale du coefficient correcteur peut être déterminée à partir de courbes de polarisation relatives à la cellule de référence 15 pour différentes valeurs locales de la résistance électrique. Dans cet exemple de courbes de polarisation issue d'un modèle électrochimique de cellule validé expérimentalement, la densité locale de courant diminue de 0,93A/m2 à 0,77A/m2 en multipliant la résistance électrique de la cellule par 3, pour une tension de cellule égale à 0,65V. Cette diminution de 17% de la densité locale de courant conduit à une diminution sensiblement identique 20 de la production locale de chaleur. On peut ensuite déterminer la distribution spatiale Rf',3, de la résistance électrique R par corrélation, par exemple une multiplication termes à termes, de la distribution spatiale initiale Rix,3, de la résistance électrique R avec la distribution spatiale du coefficient correcteur. Selon une deuxième variante, on identifie tout d'abord au moins une zone Z, de la cellule 25 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 Rf',3, de la résistance électrique R par modification de la distribution spatiale initiale Rix,3, dans la zone Z, identifiée en fonction 3041481 15 de la valeur locale de l'écart ATe dans cette zone. A titre d'exemple, la distribution spatiale initiale Rix,3, peut être localement modifiée à partir d'un coefficient correcteur dont la valeur est proportionnelle à celle de l'écart ATe dans cette zone. Comme dans la première variante, le coefficient correcteur varie continûment entre une valeur minimale 5 et une valeur maximale préalablement déterminées par des courbes de polarisation similaires ou identiques à celles de la figure 3. On obtient ainsi une distribution spatiale Rf',3, de la résistance électrique R de la cellule électrochimique. On est alors en mesure de modifier la distribution initiale Rix,3, de la résistance électrique R de la cellule de référence pour qu'elle présente la nouvelle 10 distribution déterminée lors de l'étape 150. La cellule ainsi optimisée 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 15 homogène, permettant ainsi de préserver la durée de vie de la cellule. La figure 4 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 la production d'énergie électrique, selon un deuxième mode de réalisation. Dans cet exemple, le paramètre d'intérêt est 20 représentatif de la résistance électrique de la cellule électrochimique, dont la valeur influe directement sur la densité de courant électrique localement produit lors de la réaction électrochimique. Il correspond ici à la résistance électrique R de la cellule. D'une manière générale, selon ce deuxième mode de réalisation, on détermine la distribution spatiale Rfx,3, de la résistance électrique R à partir de l'estimation de la 25 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 résistance électrique R de la cellule pour 3041481 16 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 résistance électrique R de la cellule est répartie spatialement selon une 5 distribution initiale Rix,y. La distribution spatiale initiale Rix,3, de la résistance électrique R 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. 10 Lors d'une étape 220, 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. 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 l''cy de la 15 température de consigne r en fonction de la distribution spatiale de la concentration en 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 surface 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 20 traduit par une diminution de la densité locale de courant produit par la cellule et donc de la production 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 25 température effective de la cellule en fonctionnement correspond à la température de consigne, celle-ci présentant avantageusement une distribution spatiale homogène. Pour cela, la distribution spatiale r-Qy de la température de consigne précisée 1"c peut s'écrire, par exemple : TXy= T),37+ Ki.[ max(c4) - ] (2) 3041481 17 où ci est la distribution spatiale de la concentration ci en espèce réactive i, par exemple x,y en oxygène, et Ki est une constante positive, par exemple proche de 1, qui peut ultérieurement être ajustée. La distribution spatiale cxi,3, de la concentration ci peut être estimée au premier ordre à partir du cheminement des canaux du circuit de distribution 5 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 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 (2) peuvent convenir pour préciser 10 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 1"xc,3, 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 à l'étape de 15 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 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 20 de distribution et 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é 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 25 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 : TXÿ=TXy+ [ (px,3, / (pin ] (3) où (px,3, est la distribution spatiale de l'humidité relative cp dans le circuit de distribution, (pin est l'humidité relative en entrée du circuit de distribution, et Kg) est une constante 3041481 18 positive, par exemple proche de 1, qui peut ultérieurement être ajustée. La distribution (px,3, de l'humidité relative cp 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 5 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 permet de déduire celle de l'humidité relative. D'autres relations (3) 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 t, y de la 10 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. Lors d'une étape 230, on obtient une distribution spatiale Dxr,y 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 15 caloporteur dans le circuit de refroidissement. Cette distribution Dxr,y n'est pas estimée mais est issue d'une 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 210, par exemple au moyen d'une technique de type PIV (Particle Image Velocimetry, en 20 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 DxrY 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. 25 Lors d'une étape 240, on estime la distribution spatiale Qe',y d'une seconde grandeur thermique Qe à partir de ladite distribution spatiale T'cy de la température de consigne T` définie lors de l'étape 220 et de ladite distribution spatiale K,y du débit de fluide caloporteur Dr obtenue lors de l'étape 230. La seconde grandeur thermique est 3041481 19 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. Pour cela, comme l'illustre la figure 5, on discrétise le circuit de refroidissement en un maillage dont chaque maille est un volume élémentaire (i,j) traversé par le fluide 5 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 échanges 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 10 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 : BTij -Tij- (4-1) 8Ti2i = (4-2) 15 8Ti3i = - Tic+1,1 (4-3) 8Tià = - (4-4) Les débits du fluide caloporteur au niveau des quatre facettes de la maille (i,j) peuvent être 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 20 mailles (i,j-1), (i,j) et (i,j+1) : dit = (Di5.ey + Di5+1.ey)12 (5-1) dizj = + Dir_ti.ex)12 (5-2) c1,3 + Dil (5-3) ,] = lE1i.ex)12 dit = (Di5.ey + Dri_rey)12 (5-4) 25 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. On calcule enfin la distribution spatiale (2,e,), du flux Qe de chaleur produit par la cellule à partir de la relation : 3041481 20 Qc,y Qie,j = dikj. Cp. ôTiki (6) On obtient ainsi la distribution spatiale du flux de chaleur Qe que doit produire la cellule 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 5 dans le circuit de distribution. Lors d'une étape 250, on détermine la distribution spatiale Rfx,3, de la résistance électrique R en fonction de la distribution spatiale 1:2'e,y du flux de chaleur produite Qe. Pour cela, on peut estimer tout d'abord la distribution spatiale de la densité d'un signal électrique produit par la cellule en fonctionnement, par exemple la densité de courant le, à partir de 10 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 : le e Itot = Qtot (7) où Itot est la densité de courant totale délivrée par la cellule électrochimique en 15 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 Rfx,3, de la résistance électrique R à partir de la densité locale du courant électrique Ixe,y. Pour cela, une approche consiste à déterminer les valeurs minimale Rmin et maximale Rmax de la résistance électrique de la 20 cellule au niveau de la zone active. Il peut s'agir d'une mesure expérimentale sur un échantillon de cellule présentant les mêmes propriétés que celles de la cellule de référence, ou d'une mesure par simulation numérique, par exemple à l'aide du logiciel COMSOL Multiphysics®. On calcule ensuite la distribution spatiale Rfx,3, par exemple à partir de la relation : 25 Rf = min (Rmax, Rmin xy irneax) (8) x,y où lie.nax est la valeur maximale de la densité locale de courant Ixe,y. La résistance électrique locale varie ainsi linéairement entre les valeurs minimale Rmin et maximale Rmax 3041481 21 en fonction de la valeur locale de la densité de courant le. Bien entendu, toute autre loi, par exemple polynomiale, exponentielle ou logarithmique, faisant varier la résistance électrique locale de sorte que la valeur minimale Rn-fin corresponde à une densité locale de courant maximale et inversement, peut convenir. 5 On a ainsi déterminé la distribution spatiale Rf',3, de la résistance électrique R qui tienne compte de la 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 r, 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 10 avantageusement homogène spatialement, une cellule dont la résistance électrique R est répartie suivant la distribution spatiale Rf',3, présente, lorsqu'elle est en fonctionnement au point de polarisation Utot et Itot, une température effective dont l'homogénéité de la distribution spatiale est optimisée. 15 Un procédé de réalisation d'une plaque bipolaire de la cellule électrochimique est maintenant décrit. 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 chaleur 20 produite par la cellule en fonctionnement, la cellule présentant une résistance électrique R répartie spatialement selon une distribution initiale Rix,y. On détermine, par le procédé décrit précédemment en référence aux figures 2 et 4, une distribution spatiale Rf',3, de la résistance électrique R de la cellule électrochimique. On réalise ensuite, lors d'une étape 160 (figure 2) ou 260 (figure 4), la plaque bipolaire de sorte que la résistance électrique R 25 présente la distribution spatiale Rfx,3, déterminée. Selon un mode de réalisation, l'une et/ou l'autre des plaques bipolaires est formée de deux tôles présentant des bossages obtenus par emboutissage ou moulage, assemblées l'une à l'autre. Il est possible de modifier la résistance de contact entre les deux tôles par la présence ou non de points de soudure (figure 6a) et/ou par la présence d'un 3041481 22 revêtement disposé localement entre les deux tôles (figure 6b) de manière à ce que la résistance électrique R présente la distribution spatiale Rf',y déterminée précédemment. Comme l'illustre la figure 6a, les deux tôles 42a, 42b, ici métalliques, peuvent être assemblées l'une à l'autre par une pluralité de points, ou zones, de soudure 44 répartis 5 selon une distribution spatiale initiale. La dimension des points de soudure 44 ainsi que leur densité modifient localement la résistance de contact RPB entre les deux tôles métalliques 42a, 42b et donc la résistance électrique R de la cellule. Plus précisément, un point de soudure entraîne une augmentation locale de la conductivité de la plaque bipolaire et donc une diminution de la résistance de contact de celle-ci. Dans cet exemple, 10 les tôles sont réalisées en un matériau métallique, les bossages étant de préférence réalisés par emboutissage. Comme l'illustre la figure 6b, les deux tôles 42a, 42b peuvent comporter localement un revêtement 45 au contact des deux tôles, dont le matériau et/ou l'épaisseur sont choisis en fonction de la modification de la résistance de contact désirée. Le matériau présente 15 une conductivité et une épaisseur choisies pour assurer une modification locale de la résistance de contact RBP de la plaque bipolaire et donc de la résistance électrique R de la cellule. En particulier, la conductivité électrique du revêtement est différente de celle des tôles de la plaque bipolaire. A titre illustratif, le matériau peut être choisi parmi les matériaux listés dans le tableau 2 de la publication de Herman et al., Bipolar plates for 20 PEM fuel cells: A review, International Journal of Hydrogen Energy, 30 (2005), 1297-1302. Le positionnement de la couche de revêtement entre les deux tôles présente notamment l'avantage d'éviter tout risque de mise en contact du matériau de revêtement avec les électrodes susceptible d'en altérer les propriétés électrochimiques. Dans cet exemple, les tôles peuvent être réalisées en un matériau autre que métallique, par exemple un 25 matériau composite chargé en graphite. Les tôles sont alors des plaques dont les bossages sont réalisés de préférence par moulage. Comme l'illustre la figure 6c, on réalise la plaque bipolaire en modifiant localement la distribution spatiale initiale de la résistance électrique R de la cellule par modification de la résistance de contact RBP entre les deux tôles. Sur la figure est représenté de manière 3041481 23 schématique en vue de dessus un circuit cathodique de distribution d'air 51 d'une plaque bipolaire qui s'étend entre une entrée et une sortie, et dont les différentes zones rectangulaires grisées représentent des zones où la résistance de contact RBP a été localement modifiée. Dans les zones Zi et Z2 où la température locale correspondant à 5 une répartition de la résistance électrique initialement homogène est susceptible d'être supérieure à la température de consigne r, on diminue la densité des points de soudure par rapport à la densité initiale, de sorte que la résistance électrique R soit maximale et limite la production locale d'énergie électrique et donc thermique. Dans la zone Z3 où la température locale est susceptible d'être inférieure à la température de consigne r, on 10 augmente la densité des points de soudure par rapport à la densité initiale, de sorte que la résistance électrique R soit minimale et permette une production locale d'énergie électrique et donc thermique maximale. Ainsi, la cellule électrochimique munie d'une plaque bipolaire dont la résistance de contact RBP présente une distribution ainsi modifiée fournit en fonctionnement une température effective sensiblement égale à la 15 température de consigne en tout point de la zone active, celle-ci étant avantageusement homogène. De manière équivalente, les zones grisées localement sur la figure 6c peut illustrer des zones où un revêtement présentant une conductivité électrique élevée a été disposé. On modifie ainsi la résistance de contact RBp de la plaque bipolaire pour ajuster la résistance 20 électrique R de la cellule de manière à ce que la température effective de la cellule soit sensiblement identique à la température de consigne en tout point de la zone active. Comme le montrent les figures 2 et 4, à l'issue de l'étape de réalisation de la cellule électrochimique pour que la résistance électrique R présente la distribution spatiale Rfx,3, 25 préalablement déterminée, une étape de mesure de la distribution spatiale de la température effective de la cellule en fonctionnement est effectuée, puis celle-ci est comparée à la température de consigne Tc. Les étapes 110, 120; 210, 220 peuvent être réitérées en modifiant les propriétés de la nouvelle cellule de référence et/ou en affinant la distribution spatiale de la température de consigne, notamment en ajustant les 3041481 24 coefficients K' et K. Les étapes suivantes du procédé de détermination de la distribution de la résistance électrique R puis de réalisation de la plaque bipolaire sont ensuite effectuées, puis est réalisée une nouvelle étape de comparaison entre la température effective de la cellule en fonctionnement et la température de consigne, jusqu'à ce qu'un 5 paramètre de convergence indique que cet écart est inférieur à une valeur seuil prédéterminée. Des modes de réalisation particuliers viennent d'être décrits. Différentes variantes et modifications apparaîtront à l'homme du métier. 10 Ainsi, le paramètre d'intérêt représentatif de la production d'énergie électrique d'une cellule électrochimique peut être plus particulièrement représentatif de la résistance de contact RCBp/GDL entre la couche de diffusion de l'électrode et la plaque bipolaire adjacente. La résistance de contact RCGDL/pg de l'électrode avec la plaque bipolaire peut être ajustée 15 en fonction de la distribution spatiale AT'e,y, Q,'',y d'une seconde grandeur thermique ATe, Qe selon que le procédé selon le premier mode de réalisation (figure 2) ou selon le second mode (figure 4) est utilisé. Pour cela, un procédé de réalisation d'une cellule électrochimique comporte une étape dans laquelle on ajuste localement l'épaisseur d'au moins l'une des plaques bipolaires. 20 Plus précisément, on diminue l'épaisseur de la plaque bipolaire dans une zone Z, dans laquelle la température effective est supérieure à la température de consigne, de manière proportionnelle à la valeur locale de cet écart. Ainsi, cette diminution d'épaisseur de la plaque bipolaire conduit à un effort diminué sur la couche de diffusion, ce qui augmente la résistance de contact RCpB/GDL et entraîne une diminution de la production locale 25 d'énergie électrique et donc thermique. A l'inverse, on augmente l'épaisseur de la plaque bipolaire dans une autre zone Z, dans laquelle la température effective est inférieure à la température de consigne, de manière proportionnelle à la valeur locale de cet écart. Ainsi, cette augmentation d'épaisseur de la plaque bipolaire conduit à un effort augmenté sur la couche de diffusion, ce qui diminue la résistance de contact RCpB/GDL et entraîne une 3041481 25 augmentation de la production locale d'énergie électrique et donc thermique. Il est à noter que la variation de l'épaisseur de la plaque bipolaire entraîne également une variation de la section droite des canaux de distribution des espèces réactives, qui induit également une variation à la production d'énergie électrique et donc thermique qui se 5 conjugue avec la variation de la résistance de contact. De manière alternative ou complémentaire, on place au moins une plaque de serrage au contact ou à distance de la cellule électrochimique, cette plaque étant structurée de manière à assurer l'application d'un effort de serrage sur la cellule entraînant des contraintes mécaniques se traduisant par des variations locales de la résistance de 10 contact RcPB/GDL. Il peut s'agir d'une ou de deux plaques de serrage placées de part et d'autre d'un empilement de cellules électrochimiques et/ou d'une plaque de serrage, ou insert, disposée entre deux cellules adjacentes de l'empilement. La plaque de serrage présente une épaisseur variable qui conduit à un effort variable sur la couche de diffusion, ce qui modifie la résistance de contact RcPB/GDL et entraîne une variation de la 15 production locale d'énergie électrique et donc thermique. Ainsi, l'application d'un effort inhomogène sur la cellule électrochimique se traduisant par des inhomogénéités de la distribution spatiale de la résistance de contact RcPB/GDL entre l'électrode et la plaque bipolaire revient à compenser les inhomogénéités de température de la cellule électrochimique de manière à ce que la température corresponde à la 20 température de consigne. En variante ou de manière complémentaire, le paramètre d'intérêt représentatif de la production d'énergie électrique d'une cellule électrochimique peut être la charge en catalyseur présent dans la couche active d'au moins l'une des électrodes, le catalyseur 25 pouvant être par exemple du platine situé dans la couche active de la cathode. Par charge, ou taux de charge, on entend la masse surfacique. La densité locale de courant électrique lx,y est en effet directement proportionnelle à la charge locale Cx,y en catalyseur. Ainsi, on détermine la distribution spatiale Cx,y de la charge en catalyseur en fonction de la distribution spatiale AT'ey, Qe',y de la seconde 3041481 26 grandeur thermique Ar ; Ce selon que l'on applique le procédé de détermination selon le premier mode (figure 2) ou selon le second mode de réalisation (figure 4). Un procédé de réalisation d'une cellule électrochimique comporte alors une étape de dépôt du catalyseur, par exemple à l'aide d'une encre catalytique sur la surface de la 5 membrane ou sur la couche de diffusion. Cette encre catalytique peut être composée notamment du catalyseur supporté sur du carbone (carbone platiné), un solvant et un polymère conducteur protonique (ionomère), en général de même nature que l'électrolyte. Le dépôt peut être effectué par des méthodes d'impression telles que l'enduction, la pulvérisation, le jet d'encre, ou la sérigraphie. Il est alors possible 10 d'effectuer plusieurs étapes de dépôt localisé de différentes encres catalytiques qui diffèrent entre elles par la charge en catalyseur. Il est par ailleurs avantageux de préserver la quantité totale de catalyseur sur la zone active de manière à ne pas modifier la densité total de courant électrique de la cellule en fonctionnement. 15 Ainsi, les variations en charge de catalyseur se traduisent par des inhomogénéités de la distribution spatiale de densité de courant et donc de flux de chaleur produite, qui compensent les inhomogénéités de température de la cellule électrochimique de manière à ce que la température corresponde à la température de consigne. 20 En variante ou de manière complémentaire, le paramètre d'intérêt représentatif de la production d'énergie électrique d'une cellule électrochimique peut être plus particulièrement représentatif de la perméabilité fluidique d'au moins l'une des électrodes de la cellule, et être par exemple la perméabilité k de la couche de diffusion de l'anode ou de la cathode. La perméabilité k illustre la capacité de la couche de permettre 25 la diffusion des fluides en son sein. La densité locale de courant électrique lx,y est en effet directement proportionnelle à la quantité d'espèce réactive ayant diffusée au travers de la couche de diffusion jusqu'à la couche active de l'électrode. Ainsi, une couche de diffusion ayant localement une 3041481 27 perméabilité k quasiment nulle ne permet pas à l'espèce réactive de rejoindre la couche active, ce qui se traduit par une densité locale de courant électrique sensiblement nulle. Ainsi, on détermine la distribution spatiale de la perméabilité k de la couche de diffusion en fonction de la distribution spatiale AT'',3,,1:2'e,y, de la seconde grandeur thermique ATe ; 5 Qe selon que l'on applique le procédé de détermination selon le premier mode (figure 2) ou selon le second mode de réalisation (figure 4). On réalise ensuite une électrode de la cellule électrochimique en modifiant localement la perméabilité k de la couche de diffusion de manière à former des inhomogénéités dans la distribution spatiale kx,y de la perméabilité k en vue de compenser les inhomogénéités de 10 la température de la cellule vis-à-vis d'une température de consigne. A titre d'exemple, on diminue localement la perméabilité de la couche de diffusion dans des zones où la température effective de la cellule est susceptible d'être supérieure à la température de consigne. Ainsi, la quantité locale d'espèce réactive pouvant diffuser jusqu'à la couche active est diminuée, ce qui diminue la densité locale de courant et donc le flux local de 15 chaleur produite. Ainsi, la température effective de la cellule correspond localement à la température de consigne. Pour cela, on imprègne par exemple la couche de diffusion d'un produit de colmatage dans les zones Z, identifiées. Le produit de colmatage peut comprendre, à titre d'exemple, du TiO2 ou de la silice colloïdale. La quantité de produit de colmatage imprégné 20 localement dans la couche de diffusion peut ainsi modifier la porosité de celle-ci pour ainsi obtenir une perméabilité locale k7 présentant la valeur désirée. Ainsi, les variations de perméabilité de la couche de diffusion de l'électrode se traduisent par des inhomogénéités de la distribution spatiale de densité de courant et donc de flux de chaleur produite, qui compensent les inhomogénéités de température de la cellule 25 électrochimique de manière à ce que la température corresponde à la température de consigne. REVENDICATIONS1. Procédé de détermination d'une distribution spatiale (R4) d'un paramètre d'intérêt (R) représentatif de la production d'énergie électrique 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, 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 (R) est réparti spatialement selon une distribution initiale (g,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 (1'3y ) de la température de consigne (r) et de ladite distribution spatiale (1;1(:y K,y) de la première grandeur thermique (Tr; D') ; v) on détermine (150; 250) une distribution spatiale ( Rfx,y ) du paramètre d'intérêt (R) 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 choisi parmi un paramètre représentatif de la résistance électrique (R) de la cellule électrochimique ; un paramètre représentatif d'une résistance de contact (RcBP/GDL) entre au moins l'une des électrodes et la plaque bipolaire adjacente ; une charge d'un 3041481 29 matériau catalyseur (Cpt) présent au moins dans l'une des électrodes ; et un paramètre représentatif d'une perméabilité (k) d'au moins une électrode. 3. Procédé de détermination selon la revendication 1 ou 2, dans lequel le paramètre d'intérêt est un paramètre représentatif de la résistance électrique (R) de la cellule 5 électrochimique, 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 (4r) représentative d'un écart local entre la température effective (r) et la température de consigne (r). 4. Procédé de détermination selon la revendication 3, dans lequel l'étape v) 10 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 (4r) ; b) une sous-étape dans laquelle on détermine la distribution spatiale (Rf',y) du paramètre d'intérêt (R) par corrélation de la distribution spatiale initiale (Rix,y) du 15 paramètre d'intérêt (R) avec la distribution spatiale du coefficient correcteur. 5. Procédé de détermination selon la revendication 3, 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 (4r) présente une valeur locale estimée supérieure 20 à une valeur seuil prédéterminée ; b) une sous-étape dans laquelle on détermine la distribution spatiale (Rf',y) du paramètre d'intérêt (R) par modification de la distribution spatiale initiale (Rix,y) dans la zone Z, identifiée, en fonction de la valeur locale estimée de la seconde grandeur thermique (4r) dans cette zone. 25 6. Procédé de détermination selon la revendication 1 ou 2, dans lequel le paramètre d'intérêt est un paramètre représentatif de la résistance électrique (R) de la cellule électrochimique, la première grandeur thermique est représentative de l'évacuation 3041481 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. 7. Procédé de détermination selon la revendication 6, dans lequel la première 5 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. 8. Procédé de détermination selon la revendication 7, dans lequel l'étape v) 10 comporte : a) une première sous-étape dans laquelle on estime la distribution spatiale de la densité d'un signal électrique (le) produit par la cellule en fonctionnement, à partir de la 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 15 paramètre d'intérêt (R) à partir de la densité locale du signal électrique (le). 9. Procédé de réalisation d'une plaque bipolaire de cellule électrochimique, comportant les étapes dans lesquelles : i) 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 20 bipolaires adaptées à amener des espèces réactives aux électrodes et à évacuer de la chaleur produite par la cellule en fonctionnement, la cellule présentant une résistance électrique (R) répartie spatialement selon une distribution initiale (g,y) ii) on détermine, par le procédé (100 ; 200) selon l'une quelconque des revendications 25 précédentes, une distribution spatiale (Rfx ) de la résistance électrique (R) ; iii) on réalise (160 ; 260) ladite plaque bipolaire de sorte que la résistance électrique (R) présente la distribution spatiale (Rfx,y) déterminée. 3041481 31 10. Procédé de réalisation d'une plaque bipolaire selon la revendication 9, ladite plaque bipolaire étant formée de deux tôles embouties assemblées l'une à l'autre par une pluralité de points de soudure répartis selon une distribution spatiale initiale, l'étape de réalisation de la plaque bipolaire comporte la modification de la distribution spatiale d'un 5 paramètre représentatif d'une résistance de contact (RcpB) entre les deux tôles en fonction de la distribution spatiale (Rfx,y) déterminée. 11. Procédé de réalisation d'une plaque bipolaire selon la revendication 10, dans lequel la modification de la distribution spatiale du paramètre représentatif de la résistance de contact (RcpB) comporte une modification de la répartition de points de soudure assurant 10 l'assemblage des tôles l'une à l'autre en fonction de ladite distribution spatiale déterminée (R4) de la résistance électrique (R). 12. Procédé de réalisation d'une plaque bipolaire selon la revendication 10 ou 11, dans lequel la modification de la distribution spatiale du paramètre représentatif de la résistance de contact (RcpB) comporte une disposition d'un revêtement présentant une 15 conductivité électrique différente de celles des tôles dans des zones préalablement identifiées à partir de ladite distribution spatiale déterminée (Rfx,y) de la résistance électrique (R). 20