TITLE:
Numerical Analysis on Temperature Distribution in a Single Cell of PEFC Operated at Higher Temperature by1D Heat Transfer Model and 3D Multi-Physics Simulation Model
AUTHORS:
Akira Nishimura, Kyohei Toyoda, Daiki Mishima, Eric Hu
KEYWORDS:
PEFC, Heat Transfer Model, Temperature Distribution, Numerical Simulation, High Temperature Operation
JOURNAL NAME:
Energy and Power Engineering,
Vol.15 No.5,
May
31,
2023
ABSTRACT: This study is to understand the impact of operating conditions, especially initial operation temperature (Tini) which is set in a high temperature range, on the temperature profile of the interface between the polymer electrolyte membrane (PEM) and the catalyst layer at the cathode (i.e., the reaction surface) in a single cell of polymer electrolyte fuel cell (PEFC). A 1D multi-plate heat transfer model based on the temperature data of the separator measured using the thermograph in a power generation experiment was developed to evaluate the reaction surface temperature (Treact). In addition, to validate the proposed heat transfer model, Treact obtained from the model was compared with that from the 3D numerical simulation using CFD software COMSOL Multiphysics which solves the continuity equation, Brinkman equation, Maxwell-Stefan equation, Butler-Volmer equation as well as heat transfer equation. As a result, the temperature gap between the results obtained by 1D heat transfer model and those obtained by 3D numerical simulation is below approximately 0.5 K. The simulation results show the change in the molar concentration of O2 and H2O from the inlet to the outlet is more even with the increase in Tini due to the lower performance of O2 reduction reaction. The change in the current density from the inlet to the outlet is more even with the increase in Tini and the value of current density is smaller with the increase in Tini due to the increase in ohmic over-potential and concentration over-potential. It is revealed that the change in Treact from the inlet to the outlet is more even with the increase in Tini irrespective of heat transfer model. This is because the generated heat from the power generation is lower with the increase in Tini due to the lower performance of O2 reduction reaction.