Lors de l'importation des données SEM des électrodes dans COMSOL, la structure poreuse ressemblait à un paysage lunaire complexe. L'approche initiale consistant à importer directement les données de niveaux de gris pour générer un modèle de pores réels a provoqué une explosion du nombre de mailles, rendant les calculs prohibitifs. Une solution de compromis a été adoptée en utilisant une caractérisation des paramètres réels - porosité de 0.32 et facteur de tortuosité de 1.8 - intégrés dans l'équation de Brinkman :
% Paramètres de la couche de diffusion du cathode
cathode.porosite = 0.32;
cathode.permeabilite = 1e-12*[0.7,0.2]; // Permeabilité anisotrope
Une observation intéressante a été que lorsque le tenseur de perméabilité était défini comme anisotrope, le champ de débit de l'hydrogène du côté de l'anode devenait remarquablement ordonné, expliquant les caractéristiques subtiles des bandes observées dans les mesures expérimentales de distribution des gaz.
Lors du couplage des champs de température, un phénomène inhabituel a été détecté : à une température de fonctionnement de 800°C, des zones de concentration de contraintes en forme de papillon sont apparues aux bords de la couche d'électrolyte. La vérification de la bibliothèque des matériaux a révélé que le coefficient de dilatation thermique par défaut pour l'oxyde de zirconium stabilisé par l'oxyde d'yttrium (YSZ) était de 10.5e-6/K, tandis que les mesures récentes de la littérature [1] indiquaient en réalité 11.3e-6/K. Après ajustement de ce paramètre et recalcul, le pic de contrainte est passé de 358 MPa à 289 MPa, correspondant aux résultats d'observations micro-CT publiés.
Lors de l'analyse des courbes de polarisation, une zone plate inattendue a été observée autour de 0.6V, rappelant une onde cardiaque irrégulière. L'examen détaillé de l'équation de conservation de la charge a révélé que le coefficient de transfert de charge aux frontières triphasées de l'anode était responsable de ce phénomène. L'ajustement des facteurs de symétrie dans l'équation de Butler-Volmer a conduit à une nouvelle compréhension :
// Réactions électrochimiques dans les champs physiques
eta_act = R*T/(alpha*n*F)*asinh(j/(2*j0));
Cette fonction sinus hyperbolique se comporte comme un régleur subtil, et lorsque la densité de courant d'échange j0 est définie comme fonction de la température, toute la courbe de polarisation commence à "danser" avec la distribution de température. Cela explique pourquoi la littérature [2] insiste sur la nécessité d'acquérir simultanément des données thermiques pendant les essais.
Lors de l'analyse post-traitement avec Paraview, des motifs de dissipation similaires à des structures fractales ont été observés dans le champ de concentration d'oxygène du côté du cathode. L'examen de la matrice des coefficients de diffusion croisée pour le transfert de masse a révélé que les termes de couplage entre la diffusion de Knudsen et la diffusion ordinaire étaient définis en mode non linéaire. Cette observation correspond aux expériences de traçage par fluorescence présentées lors des réunions de groupe, où les chemins sinueux des gaz à l'intérieur des électrodes étaient déjà prédits par le modèle mathématique.
Lorsque la dernière carte de contraintes thermiques a montré une erreur inférieure à 5% par rapport aux données expérimentales, l'analyse était considérée comme réussie. L'examen des journaux d'itération a révélé 187 tentatives de résolution infructueuses, soulignant que chaque résultat convergent représente une miniaturisation奇迹 - tout comme les ions d'oxygène dans la pile à combustible, cherchant avec persistance leurs partenaires électroniques dans le labyrinthe céramique.