To develop operating strategies in polymer electrolyte membrane (PEM) fuel cell-powered applications, precise computationally efficient models of the fuel cell stack voltage are required. Models are needed for all operating conditions, including transients. In this work, transient evolutions of voltage, in response to load changes, are modeled with a sum of three exponential decay functions. Amplitude factors are correlated to steady-state operating data (temperature, humidity, average current, resistance, and voltage). The obtained time constants reflect known processes of the membrane heat/water transport. These model parameters can form the basis for the prediction of voltage overshoot/undershoot used in computational-based control systems, used in real-time simulation. Furthermore, the results provide an empirical basis for the estimation of the magnitude of temporary voltage loss to be expected with sudden load changes, as well as a systematic method for the analysis of experimental data. Its applicability is currently limited to thin membranes with low to moderate humidity gases, and with adequately high reactant-gas stoichiometry.

Modeling works which simulate the proton-exchange membrane fuel cell with the computational fluid dynamics approach involve the simultaneous solution of multiple, interconnected physics equations for fluid flows, heat transport, electrochemical reactions, and both protonic and electronic conduction. Modeling efforts vary by how they treat the physics within and adjacent to the membrane-electrode assembly (MEA). Certain approaches treat the MEA not as part of the computational domain, but rather an interface connecting the anode and cathode computational domains. These approaches may lack the ability to consistently model catalyst layer losses and MEA ohmic resistance. This work presents an upgraded interface formulation where coupled water, heat, and current transport behaviors of the MEA are modeled analytically. Improving upon the previous work, catalyst layer losses can now be modeled accurately without ad-hoc selection of model kinetic parameters. Key to the formulation is the incorporation of water absorption/desorption resistances. The interface model is developed with the consideration of only thru-plane variation, based upon varied fundamental research into each of the relevant physics. The model is validated against differential cell data with high-and low-humidity reactants. The agreement is very good in each case.

No abstract

No abstract

No abstract