Numerical investigation of parameter distributions in high-temperature PEMFCs under various cooling surface temperature gradients

Non-uniform current density distribution can lead to localized hotspots, resulting in shortened lifespans of fuel cells. The temperature distributions on the bipolar plate cooling surface significantly influence various parameters within a fuel cell. This study develops a mathematical model for a high-temperature proton exchange membrane fuel cell (HT-PEMFC) to simulate its operation in four different directions of temperature gradient on the cooling surface. The distributions of membrane temperature, proton conductivity, oxygen concentration, activation loss, and current density were obtained to clarify the interaction mechanism between cooling surface temperature and these parameters in fuel cell. The results show that the membrane temperature and proton conductivity increase and the oxygen concentrations decrease with increasing cooling surface temperature. A lower temperature in the catalyst layer regions leads to a decrease in the current density. In particular, when the cooling surface temperature distribution is characterized by a temperature decrease in the upstream area with a high reactant concentration and an increase in the downstream area with a low reactant concentration, the fuel cell has an excellent current density uniformity, which reaches 92.61 % at a voltage of 0.6 V and a cooling surface temperature difference of 10 K.