Optimisation d’architecture d’électrode poreuse pour pile à combustible à oxyde solide

This project is involved in the development of new green power sources. Solid Oxide Fuel Cells (SOFCs) can achieve an output power of 1kW to 2MW and an energy conversion of up to 70%. Temperatures between 700 and 1000°C are required. A typical cell is made of an electrolyte sandwiched between two porous electrodes (anode and cathode). Porous electrodes are elaborated from ceramic powders and are critical components of the whole structure. These electrodes need to be porous enough to optimize gaz diffusion and electrochemical reactions. This requirement is antagonist to the need of a good mechanical strength. This conflict could be solved using hierarchical or anisotropic electrode microstructures. The aim of this thesis is to investigate possible ways to optimize an electrode. Numerical simulations and nanotomography characterizations are used for this purpose. Electrodes are elaborated using two different protocoles leading to anisotropic and isotropic porosities. Anisotropic samples are prepared by freeze-casting from a slurry of YSZ and LSM, which are typical materials for SOFCs. Freze-casting leads to a hierarchical porosity. The overall porosity is controlled by the loading of the slurry. The microporosity decreases with sintering temperature and the macropore size is function of the freezing rate. Isotropic samples are processed using pore formers. The size and the amount of pore formers are selected to match the characteristics of the anisotropic samples. These electrodes are characterized with Archimedes technique to determine the porosity, and with scanning electron microscope (SEM) to obtain the size of macropores. Three dimensional images of the microstructures are captured using focused ion beam (FIB-SEM tomography) technique (10nm} resolution) and using X-ray nanotomography (75nm} resolution). The overpotentials in an electrode depend on different parameters: composition of YSZ/LSM, porosity, particle sizes, electronic/ionic conductivities and electrochemical resistance. These parameters are studied on numerical microstructures coupled with a resistor network. These numerical microstructures have been generated at th scale of particles, using a numerical code based on the discrete element method (DEM). Simulations can be used to determine the limiting factor on the effective conductivity. For example, we show that the composition of YSZ/LSM in a sample matters little for electrodes below a certain thickness. A new method has also been developed to compute the effective conductivity from a FIB-SEM image taking into account the electrochemical resistance at the triple point boundaries between gaz, YSZ and LSM. The mechanical response of the elaborated microstructures are tested in compression up to the fracture. In parallel, DEM simulations are performed to simulate mechanical properties based on 3D images. The mechanical behaviours of homogeneous samples (with pore formers) and anisotropic samples are compared. The yield strength and stiffness are overestimated by simulations. Qualitatively, experimental results and simulations show consistent failure mecanisms. Moreover, the yield strength and stiffness are different in the two types of sample (anisotropic and isotropic). Such an anisotropy could be used to optimize mechanical properties in one direction.Ce projet se place dans le cadre du développement des nouvelles technologies de l'énergie respectueuses de l'environnement. Les piles à combustibles à oxydes solides (SOFC) permettent, pour les applications stationnaires, la génération de puissance de 1kW à 2MW avec un rendement électrique pouvant atteindre 70%. Elles fonctionnent à très hautes températures, typiquement entre 700-1000°C. La cellule d'une SOFC est constituée d'un électrolyte dense pris en sandwich entre deux électrodes poreuses (anode et cathode). Les électrodes poreuses, élaborées à partir de poudres céramiques, représentent un élément critique de l'assemblage. En effet, elles doivent être suffisamment poreuses pour optimiser à la fois la diffusion des gaz et les réactions électrochimiques. Cette nécessité est en contradiction avec l'exigence d'une bonne tenue mécanique. Cette contradiction doit pouvoir être résolue en proposant des microstructures d'électrodes poreuses hiérarchisées ou anisotropes. L'objectif de cette thèse est de montrer différentes voies possibles pour optimiser l'électrode en s'appuyant en particulier sur des simulations numériques et sur des caractérisations tomographiques. Les électrodes sont élaborées en utilisant deux protocoles différents conduisant à des porosités isotropes et anisotropes. Les échantillons anisotropes sont préparés en utilisant la méthode de moulage par congélation à partir de poudres YSZ/LSM, typiques de matériaux d'électrode. Cette méthode de fabrication conduit à une porosité hiérarchisée. La porosité totale est définie par le taux de chargement dans la barbotine initiale. La microporosité diminue avec la température de frittage et la taille des macropores est fonction de la vitesse de solidification. Les échantillons isotropes sont préparés en utilisant des agents porogènes avec des caractéristiques identiques aux échantillons anisotropes. Ces électrodes sont caractérisées par la technique d'Archimède pour déterminer les taux de porosités (macro et micro) et par microscopie à balayage pour connaître la taille des macroporosités. Des images tridimensionnelles des microstructures sont obtenues par FIB-SEM (Focused Ion Beam, 15µm³) et par nanotomographie-X (75µm³), avec des résolutions de 10nm et 75nm, respectivement. Le rendement énergétique d'une électrode dépend de différents paramètres : composition YSZ/LSM, taux de porosité, taille des particules, conductivités électronique/ionique et résistance électrochimique. Ces paramètres sont étudiés en utilisant des microstructures numériques associées à un réseau de résistance. Les simulations permettent de déterminer les facteurs qui contrôlent la conductivité effective. Ces microstructures numériques sont élaborées à l'échelle de la taille des particules en utilisant le code dp3D basé sur la méthode des éléments discrets (DEM) et développé au sein du laboratoire SIMaP. Nous montrons par exemple qu'en dessous d'une certaine épaisseur, la composition YSZ/LSM a très peu d'influence sur la conductivité effective. Une méthode a également été développée pour calculer cette conductivité effective à partir d'une image de FIB prenant en compte la résistance électrochimique aux points triples (gaz, YSZ, LSM). La tenue mécanique des différentes microstructures est testée en compression jusqu'à la rupture. En parallèle, des calculs sur image, couplés à la DEM sont effectués pour simuler les propriétés mécaniques. Nous comparons le comportement des microstructures homogènes (obtenues avec des agents porogènes) et celui des microstructures anisotropes. Les modules et les contraintes à rupture sont surestimés par les simulations. Qualitativement, les résultats expérimentaux et de simulation montrent des mécanismes de rupture cohérents entre eux. Par ailleurs, les modules et les contraintes à rupture sont différents entre les deux types d'échantillon (anisotrope et isotrope). Cette anisotropie peut être utilisée pour optimiser les propriétés mécaniques suivant une direction