Computational Fluid Dynamic Investigation of Local Flow-Field Conditions in Lab Polymer Electrolyte Membrane Fuel Cells to Identify Degradation Stressors and Performance Enhancers

Abstract

The use of polymer electrolyte membrane (PEM) fuel cells as an alternative to internal combustion engines can significantly contribute to the decarbonization of the transport sector, especially for heavy-duty applications. However, degradation is still an issue for this type of component, affecting their durability and performance. In this scenario, a detailed analysis of the anodic and cathodic distributors’ flow-field geometry may help to identify some local stressors that trigger the degradation mechanism, such as local hot spots and reactants not having a uniform distribution. A computational fluid dynamic (CFD) methodology is able to provide a volumetric description of a PEM fuel cell so it can be a useful tool to better understand the physical phenomena that govern the component operations. In this work, the open-source simulation library openFuelCell2 is adopted for a detailed analysis of two different PEM fuel cells characterized by standard distributor geometries, namely a parallel channel geometry and a serpentine configuration. The library, based on the OpenFOAM code, has been extended with a novel implementation accounting for the catalytic activity reduction due to the platinum oxide (PtOx) formation occurring under certain particular conditions. The adopted methodology is firstly validated resorting to experimental data acquired for the two different fuel cell configurations. The analysis highlights that the PtOx formation leads to a reduction in the fuel cell performance reaching up to 60–80% when operating at high voltages. Then, the effect of the distributor geometries on the component performance is investigated by resorting to in-plane and through-plane physical quantity distribution, such as reactant concentration, pressure or velocity fields. While the parallel flow channel configuration shows some diffusion losses under the rib, the serpentine channel geometry configuration can achieve some local performance peaks thanks to the convective flow in the gas diffusion layer (GDL) driven by local pressure gradients. Furthermore, the local enhancement in terms of higher current density under the rib is associated with an effective heat removal due to the high thermal capacity of the bipolar plate, avoiding the generation of local hot spots.