Elucidating non-uniform assembling effect in large-scale PEM fuel cell by coupling mechanics and performance models

The assembly of proton exchange membrane (PEM) fuel cell has a significant influence on contact pressure of BP and GDL, which affects key parameters of GDL including contact resistance, porosity, and permeability. In this study, a large-scale finite element method (FEM) PEM fuel cell model is developed considering realistic assembly, which is assembled by 14 bolts. The contact pressure at the interface of GDL and BP is converted into non-uniform contact resistance incorporated into a self-developed “3D + 1D” full PEM fuel cell multi-phase model with the active area of 108 cm2. The “3D + 1D” model simplifies part of components along through-plane direction into 1D model to improve calculation efficiency and stability. The results show that the performance reduces slightly when considering non-uniform contact resistance compared with that considering uniform contact resistance. The uniformity of current density and temperature distributions reduces evidently when considering non-uniform contact resistance affecting the performance and durability, which demonstrates the necessity of treatment of non-uniform contact resistance. Furthermore, the effects of porosity and permeability of deformed GDL under different preloading torques on performance are investigated. It is found that the contact resistance reduces first and then flats as the preloading torques increase, which results in the same trend of electronic ohmic loss, while the concentration loss nearly linearly increases. The simulation shows that the PEM fuel cell shows the best performance when the preloading torque is 3 Nm. This study provides some meaningful guidance for PEM fuel cell stack assembly.