^{1}

^{*}

^{2}

^{1}

This paper presents an experimental testing and validation results for a zero-dimensional self-humidifying PEM (Proton Exchange Membrane) fuel cell stack. The model incorporates major electric and thermodynamic variables and parameters involved in the operation of the PEM fuel cell under different operational conditions. The mathematical equations are modelled by using Matlab-Simulink tools in order to simulate the operation of the developed model with a commercially available 1 kW Horizon (H-1000) PEM fuel cell stack, which is used for the purposes of model validation and tuning of the developed model. The model is mathematically modelled and presented in the recent published work of authors. The observations from model simulations provide sufficient evidence and support to the results and observations obtained from testing 1 kW Horizon (H-1000) PEM fuel cell stack used in this research. The developed model can be used as a generic model and simulation platform for a self-humidifying PEM fuel cell with an output power varying from 50 W to 1 kW, with extrapolation to higher powers is also possible.

A fuel cell is a device which directly converts the energy in the reactants into electricity. The efficiency of energy conversion-production for the fuel cell is comparatively higher than internal combustion engines because there is no intermediate thermal conversion process similar to the internal combustion engines or gas turbines [

PEM fuel cells have several features comparing to other types of fuel cells, such as: low corrosion, low weight and compact size, which make them suitable for mobile and stationary applications. The operational temperatures of the PEM fuel cell range between 30˚C - 100˚C, with a dynamical response faster than solid oxide fuel cell which operates at temperatures over 700˚C [

A one dimensional isothermal steady-state model for a PEM fuel cell with Nafion117 membrane has been developed [

A non-linear dynamic model of the PEM fuel cells using electrochemical-thermodynamic and zero-dimensional fluid mechanics principles has been developed [

A simple dynamic electrical model of a PEM fuel cell has been developed by [

A new dynamic model of 20 cells stack has been proposed [

A three-dimensional multi-phase fuel cell model has been developed [

Although many of the fuel cell models are available in literatures, but most of these models estimate the output voltage of the fuel cell for a particular set of operating conditions and some often for a limited dynamic variations. The proposed model in this paper presents a simplified zero-dimensional mathematical model for a self-humidifying 1 kW PEM fuel cell developed by modelling the major electric and thermodynamic variables and parameters involved in the operation of a PEM fuel cell.

A simplified mathematical model of the self-humidifying PEM fuel cell is developed by modelling the major electric and thermodynamic variables and parameters, as it was mathematically modelled and presented in the recent published work of the authors in [

Matlab-Simulink tools is used in order to simulate the operation of the developed model with a commercially available 1 kW Horizon (H-1000) PEM fuel cell stack which is used as an experimental device in running the test physically in order to validate and tune the developed model. The electric and thermodynamic governing equations of the developed model were mathematically modelled and used to determine the output voltage of the PEM fuel cell stack as a function to the changes of: fuel cell temperature, water content in the membrane, dimensions and volumes of the fuel cell and membrane, drawn current from the cell, and pressure and the flow rates of the hydrogen and oxygen inside the fuel cell, as presented in the previously published work of the authors in [

The output voltage of the PEM fuel cell V_{fc} is subject to various losses, and could be expressed by the following equations.

where, E_{oc} represents open circuit voltage of the PEM fuel cell as it has been derived previously for a liquid state of water produced by the fuel cell and given by [_{act}, V_{ohm}, and V_{con}, represent the activation overvoltage, ohmic overvoltage, and concentration overvoltage of the PEM fuel cell, respectively. n is the number of cells connected in series in the stack. T is stack temperature in Kelvin which is nearly equal to cell temperature. I is the drawn current in Ampere, and ζ_{n} represents the empirical parametric coefficient based on the experimental data, which may vary from one stack to another or cell to another depends on the geometrical design and the materials used in the construction of the PEM fuel cell. A_{fc} is the membrane active area (cm^{2}), and l is the thickness of membrane (cm). R is the universal gas constant, and F is Faraday’s constant.

The mass flow rate of air flows between the exit of supply manifold W_{sm}_{,out,ca} and the cathode is determined as given in the equation below [

where, P_{ca}_{,in} is the pressure of air enter the cathode, and W_{sm}_{,out,ca} is assumed to be equal to W_{sm}_{,in,ca} under condition of steady flow. k_{sm}_{,}_{out}_{,ca} is the nozzle constant of supply manifold outlet (kg/s・kPa) which represents the ratio of mass flow rate of air to the pressure. While, the pressure of air enter and exit the cathode can be determined as given in equation below [

where, m_{w}_{,gen}, m_{w}_{,mbr}, and m_{O}_{2,rct} represent mass (kg) of produced water as a result of electrochemical reaction, the mass of water vapour across the electrolyte membrane, and the mass of reacted oxygen in the cathode, respectively, V_{ca} is the volume of cathode (m^{3}), and T_{st} is the stack temperature. M_{w} and M_{O}_{2} represents the molar mass of water and oxygen (kg/mol), respectively. The mass flow rate of hydrogen at the exit of the supply manifold of anode W_{sm}_{,out,an} is determined by the equation below [

where, P_{an}_{,in} is the pressure of hydrogen enter the anode, and W_{sm}_{,out,an} is assumed to be equal to W_{sm}_{,in,an} under condition of steady flow. k_{sm}_{,}_{out}_{,an} represents the nozzle constant of supply manifold outlet of anode (kg/s・kPa), which represents the ratio of mass flow rate of hydrogen to the pressure. While, the pressure of hydrogen enter and exit the anode can be determined as given in equation below [

where, V_{an} is the volume of anode (m^{3}), m_{H}_{2,rct} represents mass of the reacted hydrogen (kg), and M_{H}_{2} is the molar mass of pure hydrogen (kg/mol), and W_{H}_{2,rct} is the mass flow rate of the reacted hydrogen as a result of electrochemical reaction.

In this research, Horizon (H-1000) fuel cell stack system is adopted as an experimental prototype, which is designed by the manufacturer to be self-humidified fuel cell stack [

Horizon (H-1000) 1 kW PEM fuel cell stack is used as an experimental device in running the test physically in order to validate and tune the developed mathematical model of PEM fuel cell. The developed model is simulated under operating conditions of room temperature, stack temperature, and drawn current similar to the real operating condition of the tested Horizon stack in the laboratory.

Different values of tuning are adopted for the developed model, in order to find the best tuning value, which enables the developed model to perform and produce a steady state output voltages and currents close to the output voltages and currents of the tested stack. Moreover, in order to get the most accurate real data outputs, the test of the PEM fuel cell stack is repeated several times on different occasions, in order to achieve a minimum deviation between the output voltages of the tested stack and output voltages of the developed model, also a closeness in the behaviour and steady state performance between the tested stack and the developed model.

Horizon fuel cell stack (H-1000) is designed to have four fans installed at the exit outlet of the return manifolds of cathodes. Hence, fans are configured as a ducted inlet-free outlet, working as suction devices at the outlet of the return manifolds of cathodes. The rotational speed of the four axial fans and the frequency of purging for the outlet hydrogen valve every 10 seconds are controlled by an inbuilt stack’s controller, in order to maintain enough operational pressure at the cathode and anode, and also to provide a sufficient flow of air and hydrogen leading to maintaining a certain level of stack temperature and retaining the continuity of the electrochemical reaction. The controller provides a protection shutdown for the stack at 30 A over current and 36 V low voltage operation [

An AttoPilot 50 V/90 A, DC voltage-current senor with 2 analog outputs is used to measure the output voltages and currents from the stack via stepping down the detected voltage and current at a ratio of (63.69 mV/1V), and (36.6 mV/1A). A USB interface data acquisition (NI USB-6008) with 8 analog inputs and 2 analog outputs ports is used to capture the real time currents and voltages, while LabVIEW v13.0 software is used for the purpose of data recording and analysis, with sampling frequency set on 1 Hz. SkyRC i-Meter operates for a maximum input voltage and current of 60 V and 100 A is used to manually measure the output voltages and currents from the fuel cell stack. DC-DC Convertor (Mean Well SD-1000L-24) with an input range of 19 V to 72 V and output of 24 V, with 40 A maximum current, is installed between the fuel cell stack and the BLDC motor in order to stabilise and maintain the output voltage.

A three phases, 14 poles, brushless DC Motor BLDC (KMS Quantum 4130/07) is used as an attached load to the Horizon fuel cell stack. The motor operates at supply voltage of 18.5 - 29.6 V, with a maximum continuous current of 40 A and a maximum efficiency higher than 90%, and the ratio of rotational speed (rpm) to the supply voltage is 360/V. Al’s Hobbies Professional analogue-digital servo tester is used to adjust the rotational speed of the BLDC motor and its propeller via generating a PWM signal to the electronic speed controller (ESC) in order to increase-decrease the current drawn from the fuel cell stack.

A pressure reduce valve (Swagelok) is used to maintain the supply pressure of hydrogen to the stack at 55 kPa. A temperature and humidity data-logger (KTH-300 Kistock), integrated with thermo-hygrometry (TH) probe sensor, is used to measure the temperature of the exit air from the stack. TH sensor is mounted at each fan outlet in an attempt to obtain an accurate estimate of temperature. The stack temperature is determined by taking the average of the sums of temperature readings for these four TH sensors.

during the testing process carried out in the laboratory. The defined values of variables and constants for the operational parameters involved in the simulation of the developed mathematical model of PEM fuel cell stack are presented in

Horizon (H-1000) PEM fuel cell stack is tested under atmospheric pressure and 23˚C room temperature. The stack is operated for about 25 minutes under different levels of load current varied from 1 to 17 A in step of 1 A, with approximate holding time of one and a half minute. One TH sensor is installed at the exit outlet of the fan in order to measure the temperature of the stack.

Different values of stack efficiency (83%, 84%, 85%, and 86%) are adopted for the developed model of PEM fuel cell stack in order to find the best tuning value between the tested stack and the developed model. It has been found that the best value of efficiency for the developed model that enables the model to perform and produce output voltages close to the stack output voltages was at 84% with 0.78 V average of deviation between the real test output voltages and output voltages of the developed model.

The output voltages for the tested Horizon (H-1000) fuel cell stack and the developed model are shown in

time of one and a half minute. The red line represents output voltages of mathematical developed model based on 84% model tuning efficiency.

It is clear that the performance and output voltages of the mathematical model developed for a PEM fuel cell is fairly close to the performance and output voltages of the tested Horizon stack operated under varied levels of load current and stack temperature, which provides initial satisfaction about the validity of the developed model.

Horizon (H-1000) PEM fuel cell stack is tested under atmospheric pressure and 18˚C room temperature. The stack is operated for about 38 minutes under different levels of current load varied from 0 to 25 A, with approximate holding time of two minutes. Four TH sensors are installed at the exit outlets of the four fans in order to have more accurate temperature measurements. The impact of drawn load current and the time of operation upon the temperature of Horizon (H-1000) fuel cell stack are measured and presented in Figures 4-6, respectively.

Horizon (H-1000) PEM fuel cell stack is tested under atmospheric pressure and 20˚C room temperature. The stack is operated for about 43 minutes under different levels of current load varied from 1 to 20 A, with approximate holding time of three minutes. Data acquisition is set to capture and record the output voltages and currents from the stack at every 30 seconds. Four TH sensors are installed at the exit outlets of the four fans in order to have more accurate

temperature measurements.

Different values of efficiency (83%, 84%, 85%, and 86%) for developed PEM fuel cell are adopted in order to find the best tuning value between the tested fuel cell stack and the developed model. It is noticed that the best value of efficiency for the developed model that enables the model to perform and produce output

voltages close to the real outputs of tested stack was at 85% with 1.44 V average of deviation between the stack test output voltages and output voltages of the developed model. Output voltages for the tested Horizon stack and for the developed model under varied stack temperatures and load currents, and the output voltages of the developed model at particularly 85% tuning efficiency and constant average stack Temperature (27.73˚C) are presented in

It can be noticed from

stack’s temperature will be adopted for any further simulations for the developed mathematical model of PEM fuel cell in this research.

Horizon (H-1000) PEM fuel cell stack is tested under atmospheric pressure and 21.5˚C room temperature. Test is started after leaving the fuel cell stack for about one hour from the last running test in order to refresh and rest the stack. The stack is operated for about 32 minutes under constant load current (9.08 A). Four TH sensors are installed at the exit outlets of the four fans in order to have more accurate temperature measurements, and TH sensors readings are captured every two minutes. The impact of constant drawn load current and time of operation upon the temperature of Horizon fuel cell stack are measured and drawn in

It is clear from the results presented for the tested Horizon (H-1000) fuel cell stack that the stack’s temperature has intendancy to increase with operating time and drawn current. The temperature of the stack is maintained around the level of (30˚C) even under further extend in operating time or further increase in the drawn current. This is returned to the potential role of stack’s controller which works to suppress any further increase in the stack’s temperature above (30˚C) by pumping more air to the cathode in order to maintain the operating temperature of the stack around certain level of operating temperature.

[

Furthermore, the observations from model simulations provide sufficient evidence and support to the results and observations obtained from testing 1 kW Horizon (H-1000) PEM fuel cell stack used in this research. The developed model can be used as a generic model and simulation platform for a self-humidifying PEM fuel cell with an output power varying from 50 W to 1 kW, with extrapolation to higher powers is also possible.

This paper presents a validation of the zero dimensional mathematical model for a self-humidifying PEM fuel cell stack. The major electric and thermodynamic variables and parameters involved in the operation of the PEM fuel cell with the association of the influence of the environment and conditions of fuel cell operation are considered. A 1 kW Horizon (H-1000) PEM fuel cell stack is used as a prototype device in order to validate and tune the developed mathematical model with the output results of the test.

The results of testing 1 kW Horizon stack and the simulations of the developed mathematical model are presented for different changes in the operational variables in order to improve the accuracy of the model. The observations from model simulations provide sufficient evidence and support to the results and observations obtained from testing 1 kW Horizon PEM fuel cell stack used in this research. Moreover, the developed model can be used as a generic model and simulation platform for a self-humidifying PEM fuel cell with an output power varying from 50 W to 1 kW, with extrapolation to higher powers is also possible.

Further works need to be carried out in the future to extend the developed model to implicit the impact of produced heat as a result of electrochemical reaction and the impact of the friction of the flowing reactants inside the fuel cell on the performance and power output of PEM fuel cell stack.

The authors would like to thank Dr. George Pissanidis and Leonine Kunzwa from School of Engineering & Technology at the University of Hertfordshire for their assistance and technical support during the field of experiments.

Saleh, I.M. Ali, R. and Zhang, H.W. (2018) Experimental Testing and Validation of the Mathematical Model for a Self-Humidifying PEM Fuel Cell. Journal of Materials Science and Chemical Engineering, 6, 202-218. https://doi.org/10.4236/msce.2018.64019

Variables and Constants | Defined Values |
---|---|

Number of cells in the stack | 72 |

Active area of electrolyte membrane | 80 cm^{2} |

Water content in electrolyte membrane (λ) | 7 |

Membrane thickness (L) | 25 × 10^{−}^{4} cm |

Maximum current density (i_{m}) | 500 mA/cm^{2} |

Maximum output voltage | 43.2 V at 24 A |

Maximum stack operating temperature | 65˚C |

Temperature of supply air and hydrogen | Room temperature |

Ambient pressure | 101.325 kPa |

Ambient air density | 1.225 kg/m^{3} |

Pressure of supply Hydrogen | 55 kPa |

Density of supply Hydrogen | 0.899 m^{3} |