FR 3041480 A1 20170324 FR 1558894 A 20150921 1 DETERMINATION D'UNE DISTRIBUTION SPATIALE DE LA RESISTANCE ELECTRIQUE DE CONTACT 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 résistance électrique de contact d'une cellule électrochimique permettant d'homogénéiser la distribution spatiale de température de la cellule en fonctionnement, ainsi que celui des procédés de réalisation d'une cellule ou d'un réacteur é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, et 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 sous forme d'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, ce qui nécessite d'apporter les espèces réactives, d'évacuer les produits et les espèces non réactives ainsi que la chaleur produite. 3041480 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 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 à 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 de diffusion et d'une couche active. 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 réactifs 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 correspondant à la couche de diffusion. L'empilement de cellules électrochimiques est habituellement disposé entre deux plaques terminales (end-plates, en anglais) qui appliquent une force de serrage, ou effort mécanique, sur l'empilement via des éléments de serrage tels que des boulons tirants, cet 25 effort permettant d'obtenir une étanchéité adéquate en bordure des cellules électrochimiques d'une part, et un contact électrique suffisant entre les cellules d'autre part. L'article de Wang et al. intitulé Experimental study on clamping pressure distribution in PEM fuel cells, Journal of Power Sources, 179 (2008) 305-309 montre que la performance d'une pile à combustible est meilleure lorsque la force de serrage induit au sein des cellules 3041480 3 une contrainte mécanique en compression répartie spatialement de manière homogène. Pour cela, l'article de Yu et al. intitulé Composite endplates with pre-curvature for PEMFC (polymer electrolyte membrane fuel cell), Composite Structures 92 (2010) 1498-1503 décrit des plaques terminales initialement convexes de sorte que, une fois assemblées à 5 l'empilement, elles présentent une courbure sensiblement nulle, entraînant ainsi la répartition sensiblement uniforme de la force de serrage. Par ailleurs, on cherche à 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 cellule en fonctionnement. En effet, les inhomogénéités spatiales de densité de courant 10 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 thermique susceptibles de fragiliser la tenue mécanique des éléments de la cellule. A titre illustratif, le document FR2976732 décrit une cellule électrochimique réalisée de 15 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 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 20 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 25 partielle d'oxygène. 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 3041480 4 distribution spatiale d'un paramètre représentatif d'une résistance électrique de contact 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 d'une résistance de contact entre une plaque bipolaire et une électrode adjacente d'une cellule électrochimique, ladite cellule comportant deux électrodes séparées l'une de l'autre par un électrolyte et disposées entre des plaques bipolaires adaptées à amener des espèces réactives aux électrodes et à 10 évacuer de la chaleur produite par la cellule en 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 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. Le paramètre d'intérêt est de préférence la résistance électrique de contact de la cellule 25 électrochimique, 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 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 : 3041480 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 5 avec la distribution spatiale du coefficient correcteur. 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 ; 10 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 identifiée, en fonction de la valeur locale estimée de la seconde grandeur thermique dans cette zone. De préférence, le paramètre d'intérêt est la résistance électrique de contact de la cellule 15 électrochimique, 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. 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 cellule, 20 et la seconde grandeur thermique peut être le flux local de chaleur produite par la cellule en fonctionnement. De préférence, l'étape v) comporte : 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 25 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 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 cellule électrochimique, comportant les étapes dans lesquelles : 3041480 6 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 des plaques bipolaires adaptées à amener des espèces réactives aux électrodes et à évacuer de la chaleur produite par la cellule en fonctionnement, les plaques bipolaires 5 présentant une distribution spatiale initiale d'épaisseur, la cellule présentant un paramètre d'intérêt représentatif de la résistance électrique de contact réparti spatialement selon une distribution initiale ; on détermine, par le procédé selon l'une quelconque des caractéristiques précédentes, une distribution spatiale du paramètre d'intérêt ; 10 on réalise ladite cellule à partir de la cellule de référence, de sorte que le paramètre d'intérêt présente la distribution spatiale déterminée. De préférence, on détermine une distribution spatiale d'épaisseur d'au moins l'une des plaques bipolaires en fonction de la distribution spatiale déterminée. De préférence, ladite plaque bipolaire présente localement une épaisseur inférieure ou 15 supérieure à une épaisseur nominale de la plaque sensiblement égale à l'épaisseur initiale, dans des zones préalablement identifiées à partir de ladite distribution spatiale déterminée de la résistance électrique de contact. L'invention porte également sur un procédé de réalisation d'un réacteur électrochimique, comportant les étapes dans lesquelles : 20 on considère un réacteur électrochimique de référence comportant un empilement de cellules électrochimiques enserrées et comprimées entre deux plaques terminales, chaque cellule comprenant deux électrodes séparées l'une de l'autre par un électrolyte et disposées entre des plaques bipolaires adaptées à amener des espèces réactives aux électrodes et à évacuer de la chaleur produite par la cellule 25 en fonctionnement, et présentant un paramètre d'intérêt représentatif de la résistance électrique de contact réparti spatialement selon une distribution initiale ; on détermine, par le procédé selon l'une quelconque des caractéristiques précédentes, une distribution spatiale du paramètre d'intérêt pour au moins l'une des cellules ; 3041480 7 - on réalise ledit réacteur électrochimique à partir du réacteur de référence, de sorte que le paramètre d'intérêt présente la distribution spatiale déterminée. De préférence, on détermine une distribution spatiale d'épaisseur d'au moins l'une des plaques terminales, et/ou une distribution spatiale d'épaisseur d'au moins une plaque 5 insérée entre deux cellules adjacentes dudit empilement, en fonction de la distribution spatiale déterminée. L'invention porte également sur un cellule électrochimique comportant deux électrodes séparées l'une de l'autre par un électrolyte et disposées entre des plaques bipolaires adaptées à amener des espèces réactives aux électrodes et à évacuer de la chaleur produite 10 par la cellule en fonctionnement, caractérisée en ce qu'au moins l'une des plaques bipolaires présente une épaisseur locale inférieure ou supérieure à une épaisseur nominale de la plaque, dans au moins une partie de sa surface. L'invention porte également sur un réacteur électrochimique comportant un empilement de cellules électrochimiques enserrés et comprimés entre deux plaques terminales, chaque 15 cellule comprenant deux électrodes séparées l'une de l'autre par un électrolyte et disposées entre des plaques bipolaires adaptées à amener des espèces réactives aux électrodes et à évacuer de la chaleur produite par la cellule en fonctionnement, caractérisée en ce qu'au moins l'une des plaques terminales présente une face en contact avec l'empilement de cellule dont la surface est localement structurée de manière à 20 présenter des protubérances ou des échancrures dans au moins une partie de sa surface, et/ou en ce qu'elle comporte une plaque insérée entre deux cellules adjacentes dont au moins l'une de ses faces en contact avec les cellules présente une surface localement structurée de manière présenter des protubérances ou des échancrures dans au moins une partie de sa surface. 25 BRÈVE DESCRIPTION DES DESSINS 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, 3041480 8 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 électrochimique, et la figure lb est une représentation schématique illustrant la relation 5 de corrélation entre la distribution spatiale de la production de chaleur Q et la distribution spatiale d'évacuation D de chaleur dont résulte la distribution spatiale de température T de la cellule électrochimique en fonctionnement ; la figure 2 est un organigramme d'un procédé de détermination de la distribution spatiale de la résistance électrique de contact de la cellule électrochimique selon un premier mode 10 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 ; la figure 4 est un organigramme d'un procédé de détermination de la distribution spatiale 15 de la résistance électrique de contact 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 chaleur et une température locale Tic,i ; 20 les figures 6a et 6b illustrent de manière schématique et partielle une plaque bipolaire située entre deux cellules électrochimiques adjacentes, la plaque bipolaire présentant une épaisseur sensiblement constante spatialement (fig.6a) ou qui varie localement (fig.6b) ; la figure 7a illustre des exemples de relation entre la résistance électrique d'une cellule électrochimique en fonction de la contrainte de compression subie par la cellule, pour 25 différents types de couche de diffusion ; et la figure 7b illustre des exemples de relation entre l'épaisseur d'une couche de diffusion de cellule électrochimique en fonction de la contrainte de compression subie par la cellule, pour différents types de couche de diffusion ; 3041480 9 la figure 8 illustre de manière schématique et partielle un réacteur électrochimique comportant un empilement de cellules comprimé entre deux plaques bipolaires à surface structurée et comportant un insert à surface structurée. 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. On décrira différents modes de réalisation et variante en référence à une pile à 10 combustible, et en particulier à une pile à hydrogène de type PEM (pour Proton Exchange 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 à 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 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 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 à l'autre, 25 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 de plaque bipolaire peuvent être utilisées, 3041480 10 par exemple les plaques réalisées en matériau composite chargé par exemple de graphite, où les bossages sont réalisés par moulage. Comme mentionné précédemment, 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 5 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 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 10 sont le lieu des réactions électrochimiques. Elles 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. 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 15 conducteur, par exemple une matrice carbonée. La cellule électrochimique présente une résistance électrique R qui dépend notamment de la somme des résistances électriques ohmiques des différents éléments de la cellule (plaques bipolaires RPB, membrane Rm, électrodes RGDL) ainsi que de la résistance de contact Rc entre chaque plaque bipolaire et la couche active de l'électrode adjacente. La résistance 20 de contact Rc dépend notamment de la surface de contact entre la plaque bipolaire et la couche de diffusion, ainsi que de la contrainte en compression issue d'une force de serrage appliquée à la cellule et notamment subie par la couche de diffusion. La résistance de contact comporte ainsi un terme lié à l'interface entre la plaque bipolaire et la couche de diffusion et un terme lié à la résistance électrique propre de la couche de diffusion. Elle 25 peut être mesurée sur des échantillons des matériaux constituant la cellule, afin de connaître son évolution en fonction de l'effort mécanique appliqué au niveau du contact. 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 3041480 11 représentative de la production de chaleur par la cellule, par exemple 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 de la demande FR2976732 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 d'éventuelles inhomogénéités locales dans le terme Q de production de chaleur ainsi que 10 de la présence d'éventuelles 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) (1) 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 3041480 12 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 représentatif de la résistance électrique de contact d'une cellule électrochimique, en particulier de la résistance de contact Rc entre la plaque bipolaire et une électrode adjacente. Le paramètre d'intérêt est donc également représentatif de la production d'énergie électrique, et donc d'énergie thermique. On cherche en particulier à déterminer la distribution spatiale du paramètre d'intérêt pour que la distribution spatiale de la 10 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. 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 s'agit ici en particulier de la 15 résistance électrique de contact Rc de la cellule entre la plaque bipolaire et l'électrode adjacente, notamment sa couche de diffusion, ou d'un paramètre équivalent tel que la conductance de contact. En effet, plus la résistance de contact Rc est petite, meilleure sera la densité de courant produit par la cellule en fonctionnement, et inversement. La valeur de la résistance de contact dépend notamment de la force de serrage exercée sur la cellule, 20 c'est-à-dire sur les deux plaques bipolaires en direction l'une de l'autre : plus la force de serrage est élevée, plus petite sera la résistance de contact, ce qui se traduit par une densité de courant élevée. Par température de la cellule, on entend la température locale, c'est-à-dire la distribution spatiale de température de l'un ou l'autre des éléments de la cellule, par exemple de l'une 25 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 circuit de refroidissement. La température effective de la cellule est la distribution spatiale de température de la cellule en fonctionnement, au point de polarisation définie par la tension de la cellule Utot et le 3041480 13 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 5 en particulier du débit local, massique ou volumique, du fluide caloporteur circulant dans le circuit de refroidissement. Enfin, 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 cellule, au niveau de la zone dite active correspondant à l'étendue surfacique des couches 10 actives des électrodes. Ainsi, la cellule électrochimique dont le paramètre représentatif de la résistance de contact est réparti spatialement selon la distribution ainsi déterminée présente une température effective, ou température lors du fonctionnement de la cellule, sensiblement identique à la température de consigne. Cette température de consigne présente avantageusement 15 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 peuvent augmenter la cinétique des réactions de dégradation 20 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 présente 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 25 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,1 à 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 3041480 14 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. La figure 2 est un organigramme d'un procédé de détermination de la distribution spatiale 5 du paramètre d'intérêt représentatif de la résistance électrique de contact de la cellule électrochimique, selon un premier mode de réalisation. Dans cet exemple, le paramètre d'intérêt est la résistance électrique de contact Rc de la cellule électrochimique entre la plaque bipolaire et l'électrode adjacente, et notamment la couche de diffusion de celle-ci, dont la valeur Rcx,y influe directement sur la densité locale de courant électrique lx,y 10 localement produit lors de la réaction électrochimique et donc sur le flux local de chaleur produite Q. D'une manière générale, selon ce premier mode de réalisation 100, on détermine une distribution spatiale optimisée Rcxf,y de la résistance de contact Rc à partir de l'estimation de la distribution spatiale AT'e3, d'un écart ATe entre une température effective r de la 15 cellule en fonctionnement dans laquelle la résistance de contact Rc est répartie spatialement selon une distribution initiale donnée, et une température de consigne r prédéterminée. On est alors en mesure de réaliser une cellule dont la résistance de contact Rc présente la distribution optimisée Rcxf,y, de manière à ce que la température effective r de la cellule modifiée soit sensiblement égale à la température de consigne r. 20 Lors d'une première étape 110, on définit une cellule électrochimique de référence au sein de laquelle la résistance de contact Rc est répartie spatialement selon une distribution initiale Rc'i,y. La cellule présente une structure identique ou similaire à celle décrite en référence à la figure 1. La distribution spatiale initiale Rc'i,3, de la résistance électrique peut être sensiblement homogène, c'est-à-dire ici qu'elle présente une valeur sensiblement 25 constante en tout point de la zone active. A titre d'exemple, on considère que les surfaces de contact entre la plaque bipolaire et la couche de diffusion sont uniformément réparties, et que la force de serrage appliquée à la cellule est sensiblement uniforme. 3041480 15 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 5 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, 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 10 délivrée lors du fonctionnement de la cellule. La puissance thermique Pth est estimée à partir de l'intégration sur la zone active du flux local de chaleur produite Clx,y obtenu à partir de la relation (1). 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é 15 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. 20 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 25 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 3041480 16 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, 5 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 les inhomogénéités locales de l'évacuation de chaleur qui dépend notamment du débit local du fluide caloporteur dans 10 le circuit de refroidissement. 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 15 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. 20 Lors d'une étape 150, on détermine ensuite la distribution spatiale Rc'f,y de la résistance de contact Rc en fonction de la distribution spatiale AT'e,3, de l'écart local 4Te. Selon une première variante, on calcule tout d'abord un coefficient correcteur, ou facteur de pondération, dont la distribution spatiale peut être proportionnelle termes à termes à celle AT'e,3, de l'écart local 4Te. A titre d'exemple, le coefficient correcteur varie 25 continument entre une valeur minimale, par exemple inférieure ou égale à l'unité, jusqu'à une valeur maximale, à mesure que l'écart local 4Te 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 3041480 17 à la cellule de référence pour différentes valeurs locales de la résistance électrique de la cellule. Dans cet exemple de courbes de polarisation obtenues à l'aide 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 5 tension de cellule égale à 0,65V. Cette diminution de 17% de la densité locale de courant conduit à une diminution sensiblement identique de la production locale de chaleur. On détermine ensuite la distribution spatiale Rc'f,y de la résistance électrique Rc par corrélation, par exemple par multiplication termes à termes, de la distribution spatiale initiale Rc'i,y de la résistance électrique Rc avec la distribution spatiale du coefficient 10 correcteur. 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 Rcxf,y de la résistance électrique Rc par modification de la distribution 15 spatiale initiale Rc'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 Rc'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 20 préalablement déterminées par des courbes de polarisation similaires ou identiques à celles de la figure 3, moyennant une étape intermédiaire de soustraction des résistances ohmiques des éléments de la cellule pour déduire de la résistance électrique totale la valeur de la résistance de contact Rc. On obtient ainsi une distribution spatiale Rc'f,y de la résistance de contact Rc de la cellule 25 électrochimique. On est alors en mesure de modifier la distribution initiale Rcxi,y de la résistance de contact Rc 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 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 3041480 18 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. 5 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 résistance électrique de contact de la cellule électrochimique, selon un deuxième mode de réalisation 200. Dans cet exemple, le paramètre d'intérêt est la résistance électrique de contact de la cellule électrochimique entre la plaque bipolaire et l'électrode adjacente, et notamment la couche de diffusion de 10 celle-ci, dont la valeur influe directement sur la densité de courant électrique localement produit 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 Rcxf,y de la résistance de contact Rc à partir de l'estimation de la distribution spatiale de la production de chaleur nécessaire pour obtenir la distribution 15 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 réaliser une cellule présentant une résistance de contact Rc ayant la distribution spatiale déterminée, de manière à ce que la température effective soit alors sensiblement égale à la température de consigne. 20 Lors d'une première étape 210, on définit une cellule électrochimique de référence au sein de laquelle la résistance de contact Rc entre une plaque bipolaire et la cellule adjacente de la cellule est répartie spatialement selon une distribution initiale Rc'i,y. La distribution spatiale initiale Rc'i,y de la résistance de contact Rc peut être sensiblement homogène en scalaire, de sorte que la valeur locale est sensiblement constante en tout point de la zone 25 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 consigne r de la cellule de référence lorsque celle-ci est en fonctionnement et produit une 3041480 19 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 température de consigne r en fonction de la distribution spatiale de la concentration en 5 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 10 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 15 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 : tc,y = Txc,v Ki.L max(cix ) - ci (2) ,y x,y où ci est la distribution spatiale de la concentration c' en espèce réactive i, par exemple x,y en oxygène, et K' est une constante positive, par exemple proche de 1, qui peut 20 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 de l'espèce réactive i 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 25 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 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 3041480 20 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. Par ailleurs, de manière facultative et éventuellement complémentaire à l'étape de précision de la température de consigne décrite précédemment, il est avantageux de 5 préciser la distribution spatiale l''cy 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 de distribution sur la pression de vapeur saturante Psat. L'humidité relative cp a un effet sur la cinétique des 10 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 fonctionnement présente une distribution spatiale homogène. Pour cela, la distribution spatiale 1"),',c3, de la 15 température de consigne précisée 7.1"'c peut s'écrire, par exemple : TXÿ=TXy+K`P.[(Px,yI(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 1(q) est une constante positive, par exemple proche de 1, qui peut ultérieurement être ajustée. La distribution 20 (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 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 25 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 y de la température de consigne 7.1"c 3041480 21 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 5 fonctionnement. La première grandeur thermique est ici le débit massique Dr de fluide 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 10 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,y 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. 15 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 représentative de la production locale de chaleur et correspond ici au flux local de chaleur 20 Ce qu'évacue le fluide caloporteur Dr pour obtenir la température de consigne Tc. Pour cela, comme l'illustre la figure 5, on discrétise un modèle du 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 ; 25 et une grandeur à déterminer : le flux local de chaleur produite Qr. On calcule ensuite les 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. 3041480 22 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 : = Tc- - Tic J +1 (4-1) 8Ti2 = Tc- - (4-2) 5 8Ti3J = Tc- - (4-3) 8T j = Tc- - (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 mailles (i,j-1), 10 (i,j) et (i,j+1) : dil = (Di5. ey + +1. ey) /2 (5-1) dit = (Di5. ex + Di 1,1.ex) /2 (5-2) dit = (Di5. ex + ex) /2 (5-3) dij = (Di5. ey + ey) /2 (5-4) 15 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 1:2'e,y du flux Qe de chaleur produit par la cellule à partir de la relation : Qc,y Qie,i = dikj. Cp. ÔT (6) 20 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 dans le circuit de distribution. Lors d'une étape 250, on détermine la distribution spatiale Rc'f,y de la résistance de contact 25 Rc en fonction de la distribution spatiale Q'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 3041480 23 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 proportionnelle à la densité de courant le, celle-ci peut être déterminée à partir de la relation : le e Itot = Qtot 5 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 Rcxf,3, de la résistance de contact Rc à partir de la densité locale du courant électrique Ixe,y. Pour cela, une approche consiste à 10 déterminer les valeurs minimale Rcmin et maximale Rcmax de la résistance de contact. 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 avec le logiciel ANSYS (cf. Z. Wu, Y. Zhou, G. Lin, S. Wang, S. J. Hu, An improved mode) for predicting electrical contact resistance between bipolar plate and gas diffusion 15 layer in proton exchange membrane fuel cells, J. Power Sources, Volume 182, Issue 1, 15 July 2008, 265-269). On calcule ensuite la distribution spatiale R43, par exemple à partir de la relation : Rcf = min ( ) Rcmax, Rc xy min e le où Ime ax est la valeur maximale de la densité locale de courant Ixe,y. La résistance locale de 20 contact varie ainsi linéairement entre les valeurs minimale Rcmjn et maximale Rcmax 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 locale de contact de sorte que la valeur minimale Rcmin corresponde à une densité locale de courant maximale et inversement, peut convenir. 25 On a ainsi déterminé la distribution spatiale Rcxf,3, de la résistance électrique de contact Rc 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 (7) (8) x,y 3041480 24 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 avantageusement homogène spatialement, une cellule dont la résistance électrique Rc est répartie suivant la distribution spatiale Rcxf,3, ainsi déterminée 5 présente, lorsqu'elle est en fonctionnement au point de polarisation Utot et Itot, une température effective dont la distribution spatiale est homogène. Un procédé de réalisation d'une cellule électrochimique est maintenant décrit. On considère une cellule électrochimique de référence identique ou similaire à la cellule de 10 référence définie lors des étapes 110 et 210. La cellule de référence présente une résistance électrique de contact Rc, entre la plaque bipolaire et la couche de diffusion de l'électrode adjacente, répartie spatialement selon une distribution initiale Rcixy. On détermine, par le procédé décrit précédemment en référence aux figures 2 ou 4, une distribution spatiale Rcxf,3, de la résistance de contact Rc 15 de la cellule électrochimique. On réalise ensuite, lors d'une étape 160 (figure 2) ou 260 (figure 4), la cellule de sorte que la résistance électrique de contact Rc présente la distribution spatiale Rcxf,3, déterminée. Pour cela, le procédé de réalisation d'une cellule électrochimique comporte une étape dans laquelle on modifie localement la contrainte mécanique subie par la cellule de référence, 20 et en particulier subie par la couche de diffusion comprimée entre la plaque bipolaire et la couche active, de manière à obtenir la valeur locale déterminée Rcxf de la résistance de contact Rc. Selon une première variante faisant référence aux figures 6a et 6b, on ajuste localement l'épaisseur de la plaque bipolaire de manière à modifier la contrainte mécanique en 25 compression subie par la couche de diffusion et en conséquence la résistance de contact Rc. Cet ajustement est réalisé de préférence sans modification de l'effort de serrage appliquée à la cellule au moyen de plaques terminales entre lesquelles est disposée la cellule électrochimique. 3041480 25 En référence à la figure 6a qui illustre une partie d'un exemple de cellules de référence dans lesquelles la résistance de contact Rc présente une distribution initiale Rcxi,3, sensiblement uniforme. Sont représentées une plaque bipolaire 60 en contact, d'une part, avec une couche de diffusion anodique 11 recouverte par une couche active et une 5 membrane électrolytique 30 d'une première cellule de référence, et d'autre part avec une couche de diffusion cathodique 21 recouverte par une couche active et une membrane électrolytique d'une seconde cellule de référence. La plaque bipolaire 60 est ici formée de l'assemblage de deux tôles 64a, 64b électriquement conductrices, par exemple métalliques, présentant des bossages 65a, 65b formant un circuit de distribution 10 d'hydrogène 61 en contact avec la couche de diffusion anodique 11 et un circuit de distribution d'oxygène 62 en contact avec la couche de diffusion cathodique 21. Un circuit de refroidissement 63 assurant l'écoulement d'un fluide caloporteur est également présent. La plaque bipolaire 60 est en contact avec les couches de diffusion 11, 21 au niveau des bossages 65a, 65b formant les conduits de refroidissement 63. 15 La résistance électrique de contact entre la plaque bipolaire 60 et chaque couche de diffusion 11, 21 est, par définition, formée dans les zones de contact mécanique entre ces éléments. Dans la mesure où la plaque bipolaire 60 présente une épaisseur nominale e 130 spatialement constante, correspondant ici à la dimension des conduits de refroidissement 63 suivant un axe sensiblement orthogonal au plan moyen de la plaque bipolaire, et sous 20 l'application d'une force de serrage spatialement uniforme au niveau des cellules, les couches de diffusion 11 et 21 sont localement comprimées entre les couches actives et la plaque bipolaire de manière homogène. Par épaisseur de la plaque bipolaire, on entend la dimension de la plaque, sensiblement orthogonale au plan moyen de la plaque, dans les zones de contact avec les couches de diffusion. Les couches de diffusion présentent donc G131'i 25 localement une épaisseur ex3, uniforme dans les zones de contact mécanique avec la plaque bipolaire, et subie une contrainte mécanique de compression a également x,y uniforme. La distribution initiale Rcxi,3, de la résistance de contact Rc est en conséquence sensiblement uniforme. On suppose ici que les cellules de référence présentent une zone 3041480 26 Zi présentant une sur-température, ou point chaud, c'est-à-dire que la valeur locale de la température effective est supérieure à une valeur seuil prédéterminée. En référence à la figure 6b qui illustre la partie de cellules similaire à celle de la figure 6a, dans laquelle la résistance de contact Rc présente une distribution déterminée Rc'f,y par le 5 procédé selon l'invention. La valeur locale R41 de la résistance de contact est supérieure à la valeur initiale RcZ1 de manière à diminuer la production locale d'énergie électrique et donc thermique. Pour cela, on diminue la contrainte mécanique de compression mécanique afzi vis-à-vis de sa valeur initiale 41 en augmentant l'épaisseur locale ezGlinf de la couche de diffusion par une modification locale de l'épaisseur eziPB,f de la plaque bipolaire, 10 plus précisément par une diminution, dans la zone Z1, de la dimension transversale des conduits de refroidissement 63 suivant un axe orthogonal au plan moyen de la cellule. Autrement dit, on peut réaliser les étapes suivantes : a. Identification d'au moins une zone Z, de sur-température dans la cellule de référence dont la résistance de contact Rc présente une distribution spatiale initiale 15 R4i b. Détermination, par le procédé selon l'invention, de la valeur locale R4i de la résistance de contact Rc dans la zone Z, pour laquelle la température effective de la cellule est sensiblement égale à une température de consigne permettant ici de réduire la sur-température ; 20 c. Estimation de la contrainte mécanique locale afzi subie par la couche de diffusion permettant d'obtenir la valeur locale R4i déterminée, pour un effort de serrage identique à celui appliqué à la cellule de référence ; d. Estimation de l'épaisseur locale ezGiDL'f de la couche de diffusion permettant d'obtenir la contrainte mécanique afzi estimée ; 25 e. Calcul de l'épaisseur locale ezPi" de la plaque bipolaire permettant d'obtenir l'épaisseur ezGiDL'f estimée puis réalisation de la cellule comportant une plaque bipolaire dont l'épaisseur locale est ainsi modifiée. 3041480 27 L'estimation de l'épaisseur locale eGDL,f de la couche de diffusion permettant d'obtenir la x,Y valeur déterminée R4 de la résistance de contact peut être effectuée à partir de mesures expérimentales reliant d'une part la résistance électrique R à la contrainte mécanique de compression o subie par la couche de diffusion, et d'autre part la contrainte de 5 compression o à l'épaisseur e' de la couche de diffusion. L'article d'Escribano et al. intitulé Characterization of PEMFCs gas diffusion loyers properties, Journal of Power Sources, 156 (2006) 8-13 montre l'évolution de la résistance électrique de la cellule et l'épaisseur de la couche de diffusion en fonction de la contrainte de compression subie. La figure 7a illustre un exemple d'évolution de la résistance électrique R de la cellule en 10 fonction de la contrainte de compression subie, pour différents types de couches de diffusion, décrit dans l'article d'Escribano 2006. Comme indiqué par les auteurs, l'évolution de la résistance électrique R est représentative de celle de la résistance de contact Rc. Aussi, connaissant les différentes résistances ohmiques de la cellule, on déduit la valeur de la résistance de contact Rc en soustrayant les résistances ohmiques de la résistance 15 électrique R. A titre d'exemple, pour augmenter localement la résistance électrique de 8mQ.cm2 à 24mQ.cm2, la contrainte de compression doit passer de 2,5MPa à 0,2MPa. La figure 7b illustre un exemple d'évolution de l'épaisseur de la couche de diffusion en fonction de la contrainte de compression subie, pour différents types de couches de diffusion, décrit dans l'article d'Escribano 2006. Connaissant la contrainte mécanique de 20 compression ax,y à générer dans la couche de diffusion, on déduit l'épaisseur exGDL,fy de la couche de diffusion à obtenir. A titre d'exemple, en générant une contrainte de compression passant de 2,5MPa à 0,2MPa, l'épaisseur de la couche de diffusion doit passer de 0,22mm à 0,32mm. La plaque bipolaire est ensuite réalisée ou localement modifiée pour que, une fois la cellule assemblée et sous l'application d'une force de serrage identique à 25 celle appliquée à la cellule de référence, l'épaisseur locale de la plaque bipolaire dans la zone Z, permette d'obtenir l'épaisseur locale estimée de la couche de diffusion. Plus généralement, on augmente l'épaisseur locale de la couche de diffusion dans une zone Z, dans laquelle la température effective de la cellule de référence est supérieure à la température de consigne, de manière proportionnelle à la valeur locale de cet écart, en 3041480 28 diminuant l'épaisseur locale de la plaque bipolaire. Cette diminution d'épaisseur de la plaque bipolaire conduit à une contrainte de compression diminuée dans la couche de diffusion, ce qui augmente la résistance locale de contact Rc et entraîne une diminution de la production locale d'énergie électrique et donc thermique. A l'inverse, on peut augmenter 5 l'épaisseur de la couche de diffusion dans une autre zone Z, dans laquelle la température effective de la cellule de référence est inférieure à la température de consigne, de manière proportionnelle à la valeur locale de cet écart, en augmentant l'épaisseur locale de la plaque bipolaire. Ainsi, cette augmentation d'épaisseur de la plaque bipolaire conduit à une contrainte de compression augmentée dans la couche de diffusion, ce qui diminue la 10 résistance de contact Rc et entraîne une 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, ce qui induit également une variation à la production d'énergie électrique et donc thermique qui se conjugue avec la variation de la résistance de contact. 15 Ainsi, une cellule électrochimique dont la structure est similaire à la cellule de référence mais équipée d'au moins une plaque bipolaire dont l'épaisseur varie localement, tel que décrit précédemment, permet en fonctionnement de fournir une température effective correspondant à la température de consigne, limitant alors avantageusement les zones de sur-température. 20 Selon une seconde variante faisant référence à la figure 8, on ajuste localement l'épaisseur d'au moins une plaque terminale enserrant l'empilement de cellules et/ou l'épaisseur d'une plaque insérée 80 à l'intérieur de l'empilement de cellule, de manière à modifier la contrainte mécanique en compression subie par la couche de diffusion des cellules 25 électrochimiques et en conséquence la résistance de contact Rc entre chaque couche active d'électrode et la plaque bipolaire adjacente. Cet ajustement est réalisé de préférence sans modification de l'effort de serrage appliqué à la cellule par les plaques terminales, ni modification de l'épaisseur locale des plaques bipolaires des cellules électrochimiques. 3041480 29 En référence à la figure 8 qui illustre de manière partielle et schématique une partie d'un empilement 2 de cellules électrochimiques dans laquelle la résistance de contact Rc présente une distribution spatiale Rcxf,y déterminée à partir du procédé décrit précédemment. Sont représentés sur la figure deux plaques terminales 70a, 70b et un 5 empilement 2 de cellules électrochimiques enserré entre les plaques terminales. Les cellules présentent une structure identique ou similaire les unes aux autres et les plaques bipolaires 60 ont ici une épaisseur sensiblement uniforme. Au moins l'une des plaques terminales, ici les deux plaques 70a, 70b, présente une épaisseur ex"ya, ex"yb qui varie localement de manière à ce que la résistance de contact des cellules présente la 10 distribution spatiale Rc'f,y déterminée. Plus précisément, chaque plaque terminale 70a, 70b comporte une surface 71a, 72b en contact avec l'empilement de cellules, structurée de manière à présenter localement des protubérances ou des échancrures, qui entraînent, une fois l'empilement de cellules enserré entre les plaques terminales, une variation locale de la contrainte de compression au niveau des couches de diffusion 11, 21 se traduisant 15 par une variation de leur épaisseur locale exGDyLf et dont de la résistance de contact Rcxf,y entre chaque couche active et la plaque bipolaire adjacente. Par protubérance et échancrure, on entend une structuration locale de la surface de contact de manière à former respectivement une zone de la surface en saillie ou en retrait vis-à-vis d'un plan principal sensiblement tangent à la surface de contact. 20 En variante ou en complément, comme représenté sur la figure 8, une plaque supplémentaire 80 peut être insérée entre deux cellules adjacentes dans l'empilement de cellules, ici entre deux plaques bipolaires. De manière similaire, cet insert 80 présente une épaisseur ex8°3, qui varie localement de manière à modifier la contrainte mécanique subie par les couches de diffusion et donc leur épaisseur locale exGDL,fy , se traduisant par une 25 variation en conséquence de la résistance de contact locale Rcxf,y. Plus généralement, on diminue l'épaisseur d'au moins une plaque terminale et/ou d'un insert par rapport à un plan principal sensiblement tangent à la surface de contact de la plaque ou de l'insert, ou on réalise localement une échancrure, dans une zone Z, dans laquelle la température effective de la cellule de référence est supérieure à la température 3041480 de consigne, de manière proportionnelle à la valeur locale de cet écart. Ainsi, on diminue la contrainte de compression subie par les couches de diffusion et on en augmente l'épaisseur locale. Cela se traduit par une augmentation locale de la résistance de contact et donc à une diminution de la production locale d'énergie électrique et donc thermique. 5 Ainsi, la sur-température dans cette zone Z, est diminuée. De manière similaire à ce qui a été décrit précédemment en référence aux figures 7a et 7b, la profondeur de l'échancrure ou la hauteur de la protubérance vis-à-vis d'un plan principal sensiblement tangent à la surface de contact de la plaque ou de l'insert est déterminée en fonction de l'épaisseur estimée des couches de diffusion permettant d'obtenir la valeur locale déterminée de la 10 résistance de contact. Ainsi, un réacteur électrochimique comportant au moins une plaque terminale en contact avec un empilement de cellules électrochimiques et/ou une plaque insérée à l'intérieur de l'empilement dont l'épaisseur varie localement tel que décrit précédemment permet, en fonctionnement, de fournir une température effective correspondant à la température de 15 consigne, limitant alors avantageusement les zones de sur-température. Comme le montre les figures 2 et 4, à l'issue de l'étape de réalisation de la cellule électrochimique ou d'un réacteur électrochimique pour que la résistance électrique de contact Rc présente la distribution spatiale Rc'f,y préalablement déterminée, une étape de 20 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 coefficients K' et K. Les étapes suivantes du procédé 25 de détermination de la distribution de la résistance électrique de contact Rc puis de réalisation de la cellule ou du réacteur électrochimique 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 paramètre de convergence indique que cet écart est inférieur à une valeur seuil prédéterminée. 30 REVENDICATIONS1. Procédé de détermination d'une distribution spatiale (Rcxf,y) d'un paramètre d'intérêt (Rc) représentatif d'une résistance de contact entre une plaque bipolaire et une électrode adjacente d'une cellule électrochimique, ladite cellule comportant deux électrodes séparées l'une de l'autre par un électrolyte et disposées entre des 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 (Rc) est réparti spatialement selon une distribution initiale (R4,3,) ; ii) on définit (120 ; 220) une distribution spatiale (1'3y) d'une température de consigne (r) au sein de la cellule en fonctionnement ; iii) on obtient (130 ; 230) une distribution spatiale (1;1(:y ; K,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 (rcly ; K,y) de la première grandeur thermique (Tr; Dr) ; v) on détermine (150 ; 250) une distribution spatiale (Rcxf,y) du paramètre d'intérêt (Rc) en fonction de la distribution spatiale (AT'e3, ; Qe',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 résistance électrique de contact (Rc) de la cellule é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 (ATe) représentative d'un écart local entre la température effective (r) et la température de consigne (r). 3041480 32 3. Procédé de détermination selon la revendication 2, 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 (Rcxf,y) du 5 paramètre d'intérêt (Rc) par corrélation de la distribution spatiale initiale (R4,3,) du paramètre d'intérêt (Rc) avec la distribution spatiale du coefficient correcteur. 4. Procédé de détermination selon la revendication 2, 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 10 à une valeur seuil prédéterminée ; b) une sous-étape dans laquelle on détermine la distribution spatiale (Rcxf,y) du paramètre d'intérêt (Rc) par modification de la distribution spatiale initiale (R4,3,) dans la zone Z, identifiée, en fonction de la valeur locale estimée de la seconde grandeur thermique (4Te) dans cette zone. 15 5. Procédé de détermination selon la revendication 1, dans lequel le paramètre d'intérêt est la résistance électrique de contact (Rc) de la cellule électrochimique, 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. 20 6. Procédé de détermination selon la revendication 5, 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. 7. Procédé de détermination selon la revendication 6, dans lequel l'étape v) comporte : 25 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 3041480 33 b) une deuxième sous-étape dans laquelle on détermine la distribution spatiale du paramètre d'intérêt (Rc) à partir de la densité locale du signal électrique (le). 8. Procédé de réalisation d'une cellule électrochimique, comportant les étapes dans lesquelles : 5 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 des plaques bipolaires adaptées à amener des espèces réactives aux électrodes et à évacuer de la chaleur produite par la cellule en fonctionnement, les plaques bipolaires présentant une distribution spatiale initiale (e),ÿ') d'épaisseur, la cellule présentant 10 un paramètre d'intérêt représentatif de la résistance électrique de contact (Rc) réparti spatialement selon une distribution initiale (Rcxi,y) ; ii) on détermine, par le procédé (100 ; 200) selon l'une quelconque des revendications précédentes, une distribution spatiale (Rcxf,y) du paramètre d'intérêt (Rc) ; iii) on réalise (160 ; 260) ladite cellule à partir de la cellule de référence, de sorte que le 15 paramètre d'intérêt (Rc) présente la distribution spatiale (Rcxf,y) déterminée. 9. Procédé de réalisation d'une cellule électrochimique selon la revendication 8, dans laquelle on détermine une distribution spatiale (e),,By'f) d'épaisseur d'au moins l'une des plaques bipolaires en fonction de la distribution spatiale (Rcf ) déterminée. x,y 10. Procédé de réalisation d'une cellule électrochimique selon la revendication 9, dans 20 lequel ladite plaque bipolaire présente localement une épaisseur inférieure ou supérieure PB f (ey ' ) à une épaisseur nominale de la plaque sensiblement égale à l'épaisseur initiale x, (exP,Bii) dans des zones préalablement identifiées à partir de ladite distribution spatiale déterminée (Rcfx, ) de la résistance électrique de contact (Rc). 11. Procédé de réalisation d'un réacteur électrochimique, comportant les étapes dans 25 lesquelles : i. on considère un réacteur électrochimique de référence comportant un empilement de cellules électrochimiques enserrées et comprimées entre deux plaques 3041480 34 terminales, chaque cellule comprenant deux électrodes séparées l'une de l'autre par un électrolyte et disposées entre des plaques bipolaires adaptées à amener des espèces réactives aux électrodes et à évacuer de la chaleur produite par la cellule en fonctionnement, et présentant un paramètre d'intérêt représentatif de la 5 résistance électrique de contact (Rc) réparti spatialement selon une distribution initiale (R4,3,) ; ii. on détermine, par le procédé (100 ; 200) selon l'une quelconque des revendications 1 à 7, une distribution spatiale (Rcxf,y) du paramètre d'intérêt (Rc) pour au moins l'une des cellules ; 10 iii. on réalise (160 ; 260) ledit réacteur électrochimique à partir du réacteur de référence, de sorte que le paramètre d'intérêt (Rc) présente la distribution spatiale (Rcly) déterminée. 12. Procédé de réalisation d'un réacteur électrochimique selon la revendication 11, dans laquelle on détermine une distribution spatiale ( x e7Oya,70b) d'épaisseur d'au moins l'une des 15 plaques terminales, et/ou une distribution spatiale (en) d'épaisseur d'au moins une plaque insérée entre deux cellules adjacentes dudit empilement, en fonction de la distribution spatiale (Rc,f,y) déterminée. 13. Cellule électrochimique comportant deux électrodes (10, 20) séparées l'une de l'autre par un électrolyte (30) et disposées entre des plaques bipolaires adaptées à amener 20 des espèces réactives aux électrodes et à évacuer de la chaleur produite par la cellule en fonctionnement, caractérisée en ce qu'au moins l'une des plaques bipolaires (60') présente une épaisseur locale (ezPiB'f) inférieure ou supérieure à une épaisseur nominale (e2,,) de la plaque, dans au moins une partie (4) de sa surface. 14. Réacteur électrochimique comportant un empilement (2) de cellules 25 électrochimiques enserrés et comprimés entre deux plaques terminales (70a, 70b), chaque cellule comprenant deux électrodes (10, 20) séparées l'une de l'autre par un électrolyte (30) et disposées entre des plaques bipolaires (60) adaptées à amener des espèces réactives aux électrodes et à évacuer de la chaleur produite par la cellule en 3041480 35 fonctionnement, caractérisée en ce qu'au moins l'une des plaques terminales (70a, 70b) présente une face en contact avec l'empilement de cellule dont la surface (71a, 72b) est localement structurée de manière à présenter des protubérances ou des échancrures dans au moins une partie (Z,) de sa surface, et/ou en ce qu'elle comporte une plaque (80) insérée 5 entre deux cellules adjacentes dont au moins l'une de ses faces en contact avec les cellules présente une surface (81, 82) localement structurée de manière présenter des protubérances ou des échancrures dans au moins une partie (Z,) de sa surface. 10