Numerical Simulation of Full-Scale Cell and Stack at Very High Current Density: Influence of Heat and Liquid Water Saturation

In order to further improve power density of PEFC stack, the operations at high current density and high voltage are required in addition to downsizing. To achieve 6.0kW/L around 2030 formulated by NEDO in Japan, the extreme targets such as 0.7V@3.0A/cm 2 and 0.66V@3.8A/cm 2 being proposed. However, such operations cannot be achieved by the currently common MEA and GDL, and it is difficult to experimentally investigate the issues of internal states at such very high current density. In the present study, we carried out simulation of fuel cell performance by our originally developed software P-Stack. The results indicated that improving GDL properties is very important in addition to enhance performances of PEM and catalyst layer. The most important feature of this software is the employed macroscopic model comprehensively dealing with all the relevant transport phenomena coupled with electrochemical reactions. In anode and cathode channels and GDLs, mass transport of gas and liquid phase are solved by considering the two-phase fluid dynamics with taking into account the effects of phase changing. Here mass transport is coupled with heat and chemical species transport equations. Those all transport phenomena are coupled with electrochemical reactions in the MEA. Electrochemical reactions are modeled by Butler-Volmer equation. The parameters in the equation and other models need to be fitted to reproduce the fuel cell performance with target MEA and GDL for different operating conditions. First, we assembled 1cm 2 cell with a typical MEA (PEM: Nafion NR211, Ca/An CL: TKK TEC10E50E 0.3mgPt/cm 2 , Ca/An GDL: SGL28BC), and then measured I-V and I-R curves under various operating conditions and transient behavior when current density is rapidly increased assuming GDL flooding. These experiments were carried out with massive flow rate and constant temperature, and therefore it can be assumed the in-plane distributions are uniform. In order to reproduce the measured results of the 1cm 2 cell by simulation, we second fitted the various simulation parameters including electrochemical parameters and mass transport resistance of cathode catalyst layer and the parameters related to GDL flooding. As a validation of the fitted parameters, we compared measured and calculated I-V and I-R curves for 25cm 2 dual serpentine cell with the reference MEA. These calculated results are in good agreement with the measured values. 300cm 2 cell with reference MEA were calculated unde