Numerical study of water management in the air flow channel of a PEM fuel cell cathode

The water management in the air flow channel of a proton exchange membrane (PEM) fuel cell cathode is numerically investigated using the FLUENT software package. By enabling the volume of fraction (VOF) model, the air–water two-phase flow can be simulated under different operating conditions. The effects of channel surface hydrophilicity, channel geometry, and air inlet velocity on water behavior, water content inside the channel, and two-phase pressure drop are discussed in detail. The results of the quasi-steady-state simulations show that: (1) the hydrophilicity of reactant flow channel surface is critical for water management in order to facilitate water transport along channel surfaces or edges; (2) hydrophilic surfaces also increase pressure drop due to liquid water spreading; (3) a sharp corner channel design could benefit water management because it facilitates water accumulation and provides paths for water transport along channel surface opposite to gas diffusion layer; (4) the two-phase pressure drop inside the air flow channel increases almost linearly with increasing air inlet velocity.     

Water and pressure effects on a single PEM fuel cell

A fuel cell is a promising energy conversion system that will eventually become the first-choice for producing power because of its clean or zero-emission nature. A steady-state, two-dimensional mathematical model with pressure and phase change effects for a single PEM fuel cell was developed to illustrate the inlet humidification and pressure effects on proton exchange membrane (PEM) fuel cell performance. This model considers the transport of species along the channel as well as water transfer through the membrane. It can be used to predict trends of the following parameters along the fuel cell channels: mole number of liquid water and water vapor, pressure, temperature, density, viscosity, velocity, saturation pressure, pressure drop, vapor mole fraction, volume flow rate, required pumping power and current density.     

A critical review of two-phase flow in gas flow channels of proton exchange membrane fuel cells

Water management in PEM fuel cells has received extensive attention due to its key role in fuel cell performance. The unavoidable water, from humidified gas streams and electrochemical reaction, leads to gas-liquid two-phase flow in the flow channels of fuel cells. The presence of two-phase flow increases the complexity in water management in PEM fuel cells, which remains a challenging hurdle in the commercialization of this technology. Unique water emergence from the gas diffusion layer, which is different from conventional gas-liquid two-phase flow where water is introduced from the inlet together with the gas, leads to different gas-liquid flow behaviors, including pressure drop, flow pattern, and liquid holdup along flow field channels. These parameters are critical in flow field design and fuel cell operation and therefore two-phase flow has received increasing attention in recent years. This review emphasizes gas-liquid two-phase flow in minichannels or microchannels related to PEM fuel cell applications. In situ and ex situ experimental setups have been utilized to visualize and quantify two-phase flow phenomena in terms of flow regime maps, flow maldistribution, and pressure drop measurements. Work should continue to make the results more relevant for operating PEM fuel cells. Numerical simulations have progressed greatly, but conditions relevant to the length scales and time scales experienced by an operating fuel cell have not been realized. Several mitigation strategies exist to deal with two-phase flow, but often at the expense of overall cell performance due to parasitic power losses. Thus, experimentation and simulation must continue to progress in order to develop a full understanding of two-phase flow phenomena so that meaningful mitigation strategies can be implemented.     

Numerical studies of liquid water behaviors in PEM fuel cell cathode considering transport across different porous layers

In this paper, a two-phase two-dimensional PEM fuel cell model, which is capable of handling liquid water transport across different porous materials, is employed for parametric studies of liquid water transport and distribution in the cathode of a PEM fuel cell. Attention is paid particularly to the coupled effects of two-phase flow and heat transfer phenomena. The effects of key operation parameters, including the outside cell boundary temperature, the cathode gas humidification condition, and the cell operation current, on the liquid water behaviors and cell performance have been examined in detail. Numerical results elucidate that increasing the fuel cell temperature would not only enhance liquid water evaporation and thus decrease the liquid saturation inside the PEM fuel cell cathode, but also change the location where liquid water is condensed or evaporated. At a cell boundary temperature of 80 °C, liquid water inside the catalyst layer and gas diffusion media under the current-collecting land would flow laterally towards the gas channel and become evaporated along an interface separating the land and channel. As the cell boundary temperature increases, the maximum current density inside the membrane would shift laterally towards the current-collecting land, a phenomenon dictated by membrane hydration. Increasing the gas humidification condition in the cathode gas channel and/or increasing the operating current of the fuel cell could offset the temperature effect on liquid water transport and distribution.     

Three-dimension simulation of two-phase flows in a thin gas flow channel of PEM fuel cell using a volume of fluid method

This study investigates the two-phase flow in a thin gas flow channel of PEM fuel cells and wall contact angle's impact using the volume of fluid (VOF) method with tracked two-phase interface. The VOF results are compared with experimental data, theoretical solution and analytical data in terms of flow pattern, pressure drop and water fraction. Stable film flow is predicted, as observed experimentally, for the contact angle ranging from 5° to 40° including varying contact angles at different walls of a channel. The contact angle is found to have small impact on the gas pressure drop for the stratified flow regime, but it determines the meniscus of the two-phase interface, which affects the optical detection of the liquid thickness in experiment. The work is important to study of two-phase flow dynamics, multichannel design, experimental design and control of two-phase flows in thin gas flow channels for PEM fuel cells.     

Experimental research on water management in proton exchange membrane fuel cells

A simulated cathode flow channel experiment system was set up based on the gas flow rate and water flow rate in the PEM fuel cell. With the assistance of the visualization system, high-sensitivity double parallel conductance probes flow regime inspecting technique was adopted successfully in the experiment system to inspect the flow regime of the gas–liquid two-phase flow in the PEM fuel cell. The research results show that the double parallel conductance probes inspecting system and the flow regime image system for the gas–liquid two-phase flow in the PEM fuel cell simulated channel both can judge the slug flow and annular flow in it, and the double parallel conductance probes flow regime inspecting system can divide the annular flow into three subtypes. The main probes inspecting system and the assistant image system validate reciprocally, which enhances the experimental veracity. The typical flow regimes of the PEM fuel cell simulated channel include slug flow, annular flow with big water film wave, annular flow with small water film wave and annular flow without water film wave. With the increase of the liquid superficial velocity, the frequencies of liquid slug and wave of liquid film increase. The flow regime map in the flow channel of the PEM fuel cell was developed. The flow regime of the gas–liquid two-phase flow in a PEM fuel cell in different operating conditions can be forecasted with this map. With the PEM fuel cell operating condition in this study, the flow regimes of gas–liquid two-phase flow for different cases are all annular flow with small water film wave, and the liquid film waves more with bigger current density. With the location closer to the channel outlet, the liquid film waves are more for the same current density.     

An experimental and numerical investigation on the cross flow through gas diffusion layer in a PEM fuel cell with a serpentine flow channel

A serpentine flow channel is one of the most common and practical channel layouts for a polymer electrolyte membrane (PEM) fuel cell since it ensures the removal of water produced in a cell with acceptable parasitic load. During the reactant flows along the flow channel, it can also leak or cross to neighboring channel via the porous gas diffusion layer due to the high pressure gradient caused by the short distance. Such a cross flow leads to a larger effective flow area altering reactant flow in the flow channel so that the resultant pressure and flow distributions are substantially different from that without considering cross flow, even though this cross flow has largely been ignored in previous studies. In this work, a numerical and experimental study has been carried out to investigate the cross flow in a PEM fuel cell. Experimental measurements revealed that the pressure drop in a PEM fuel cell is significantly lower than that without cross flow. Three-dimensional numerical simulation has been performed for wide ranges of flow rate, permeability and thickness of gas diffusion layer to analyze the effects of those parameters on the resultant cross flow and the pressure drop of the reactant streams. Considerable amount of cross flow through gas diffusion layer has been found in flow simulation and its effect on pressure drop becomes more significant as the permeability and the thickness of gas diffusion layer are increased. The effects of this phenomenon are also crucial for effective water removal from the porous electrode structure and for estimating pumping energy requirement in a PEM fuel cell, it cannot be neglected for the analysis, simulation, design, operation and performance optimization of practical PEM fuel cells.     

Two-phase flow pressure drop hysteresis in an operating proton exchange membrane fuel cell

Two-phase flow pressure drop hysteresis was studied in an operating PEM fuel cell. The variables studied include air stoichiometry (1.5, 2, 3, 4), temperature (50, 75, 90 °C), and the inclusion of a microporous layer. The cathode channel pressure drops can differ in PEM fuel cells when the current density is increased along a path and then decreased along the same path (pressure drop hysteresis). Generally, the descending pressure drop is greater than the ascending pressure drop at low current densities (<200 mA cm⁻²), and the effect is worse at low stoichiometries and low temperatures. The results show that the hysteresis occurs with or without the inclusion of a microporous layer. Initial results show a modified Lockhart-Martinelli approach seems to be able to predict the two-phase flow pressure drop during the ascending path. The results compare well with photographs taken from the cathode flow field channel of a visualization cell.     

An analytical analysis on the cross flow in a PEM fuel cell with serpentine flow channel

A serpentine flow channel can be considered as neighboring channels connected in series, and is one of the most common and practical channel layouts for polymer electrolyte membrane (PEM) fuel cells, as it ensures the removal of liquid water produced in a cell with good performance and acceptable parasitic load. During the reactant flows along the flow channel, it can also leak or cross directly to the neighboring channel via the porous gas diffusion layer (GDL) due to the high‐pressure gradient caused by the short distance. Such a cross flow leads to a larger effective flow area resulting in a substantially lower amount of pressure drop in an actual PEM fuel cell compared with the case without cross flow. In this study, an analytical solution is obtained for the cross flow in a PEM fuel cell with a serpentine flow channel based on the assumption that the velocity of cross flow is linearly distributed in the GDL between two successive U‐turns. The analytical solution predicts the amount of pressure drop and the average volume flow rate in the flow channel and the GDL. The solution is validated over a wide range of the thickness and permeability of the GDL by comparing the results with experimental measurements and 3‐D numerical simulations in literature. Excellent agreement is obtained for the permeability less than 10^?9 m², which covers the typical permeability values of the GDLs in actual PEM fuel cells. The solution presents an accurate and efficient estimation for cross flow providing a useful tool for the design and optimization of PEM fuel cells with serpentine flow channels. Copyright © 2010 John Wiley & Sons, Ltd.     

Determination of the optimal active area for proton exchange membrane fuel cells with parallel,interdigitated or serpentine designs

Effects of active area size on steady-state characteristics of a working PEM fuel cell, including local current densities, local oxygen transport rates, and liquid water transport were studied by applying a three-dimensional, two-phase PEM fuel cell model. The PEM fuel cells were with parallel, interdigitated, and serpentine flow channel design. At high operating voltages, the size effects on cell performance are not noticeable owing to the occurrence of oxygen supply limit. The electrochemical reaction rates are high at low operating voltages, producing large quantity of water, whose removal capability is significantly affected by flow channel design. The cells with long parallel flow field experience easy water accumulation, thereby presenting low oxygen transport rate and low current density. The cells with interdigitated and serpentine flow fields generate forced convection stream to improve reactant transport and liquid water removal, thereby leading to enhanced cell performance and different size effect from the parallel flow cells. Increase in active area significantly improves performance for serpentine cells, but only has limited effect on that of interdigitated cells. Size effects of pressure drop over the PEM cells were also discussed.     

Treatment of two-phase flow in cathode gas channel for an improved one-dimensional proton exchange membrane fuel cell model

It has been reported recently that water flooding in the cathode gas channel has significant effects on the characteristics of a proton exchange membrane fuel cell. A better understanding of this phenomenon with the aid of an accurate model is necessary for improving the water management and performance of fuel cell. However, this phenomenon is often not considered in the previous one-dimensional models where zero or a constant liquid water saturation level is assumed at the interface between gas diffusion layer and gas channel. In view of this, a one-dimensional fuel cell model that includes the effects of two-phase flow in the gas channel is proposed. The liquid water saturation along the cathode gas channel is estimated by adopting Darcy’s law to describe the convective flow of liquid water under various inlet conditions, i.e. air pressure, relative humidity and air stoichiometry. The averaged capillary pressure of gas channel calculated from the liquid water saturation is used as the boundary value at the interface to couple the cathode gas channel model to the membrane electrode assembly model. Through the coupling of the two modeling domains, the water distribution inside the membrane electrode assembly is associated with the inlet conditions. The simulation results, which are verified against experimental data and simulation results from a published computational fluid dynamics model, indicate that the effects of relative humidity and stoichiometry of inlet air are crucial to the overall fuel cell performance. The proposed model gives a more accurate treatment of the water transport in the cathode region, which enables an improved water management through an understanding of the effects of inlet conditions on the fuel cell performance.     

Impact of channel geometry on two-phase flow in fuel cell microchannels

An important function of the gas delivery channels in PEM fuel cells is the evacuation of water at the cathode. The resulting two-phase flow impedes reactant transport and causes parasitic losses. There is a need for research on two-phase flow in channels in which the phase fraction varies along the flow direction as in operating fuel cells. This work studies two-phase flow in 60 cm long channels with distributed water injection through a porous GDL wall to examine the physics of flows relevant to fuel cells. Flow regime maps based on local gas and liquid flow rates are constructed for experimental conditions corresponding to current densities between 0.5 and 2 A cm⁻² and stoichiometric coefficients from 1 to 4. Flow structures transition along the length of the channel. Stratified flow occurs at high liquid flow rates, while intermittent slug flow occurs at low liquid flow rates. The prevalence of stratified flow in these serpentine channels is discussed in relation to water removal mechanisms in the cathode channels of PEM fuel cells. Corners facilitate formation of liquid films in the channel, but may reduce the water-evacuation capability. This analysis informs design guidelines for gas delivery microchannels for fuel cells.     

A flow channel design procedure for PEM fuel cells with effective water removal

Proper water management in polymer electrolyte membrane (PEM) fuel cells is critical to achieve the potential of PEM fuel cells. Membrane electrolyte requires full hydration in order to function as proton conductor, often achieved by fully humidifying the anode and cathode reactant gas streams. On the other hand, water is also produced in the cell due to electrochemical reaction. The combined effect is that liquid water forms in the cell structure and water flooding deteriorates the cell performance significantly. In the present study, a design procedure has been developed for flow channels on bipolar plates that can effectively remove water from the PEM fuel cells. The main design philosophy is based on the determination of an appropriate pressure drop along the flow channel so that all the liquid water in the cell is evaporated and removed from, or carried out of, the cell by the gas stream in the flow channel. At the same time, the gas stream in the flow channel is maintained fully saturated in order to prevent membrane electrolyte dehydration. Sample flow channels have been designed, manufactured and tested for five different cell sizes of 50, 100, 200, 300 and 441 cm². Similar cell performance has been measured for these five significantly different cell sizes, indicating that scaling of the PEM fuel cells is possible if liquid water flooding or membrane dehydration can be avoided during the cell operation. It is observed that no liquid water flows out of the cell at the anode and cathode channel exits for the present designed cells during the performance tests, and virtually no liquid water content in the cell structure has been measured by the neutron imaging technique. These measurements indicate that the present design procedure can provide flow channels that can effectively remove water in the PEM fuel cell structure.     

Effects of gas diffusion layer deformation on the transport phenomena and performance of PEM fuel cells with interdigitated flow fields

In this study, a three-dimensional, non-isothermal, two-phase flow mathematical model is developed and applied to investigate the effect of the GDL deformation on transport phenomena and performance of proton exchange membrane (PEM) fuel cells with interdigitated flow fields. The thickness and porosity of the GDL is decreased after compression, and the corresponding transport parameters (permeability, mass diffusivity, thermal conductivity and electrical conductivity) are affected significantly. The alterations in geometry and transport parameters of the GDL are considered in the mathematical model. The oxygen concentration, temperature, liquid water saturation and volumetric current density distributions of PEM fuel cells without compression are investigated and then compared to the PEM fuel cells with various assembly forces. The numerical results show that the cell performance is considerably improved with increasing assembly forces. However, the pressure drops in the gas flow channels are also substantially increased. It is concluded that the assembly force should be as small as possible to decrease the parasitic losses with consideration of gas sealing concern.     

Numerical investigation of the impact of two-phase flow maldistribution on PEM fuel cell performance

Flow maldistribution usually happens in PEM fuel cells when using common inlet and exit headers to supply reactant gases to multiple channels. As a result, some channels are flooded with more water and have less air flow while other channels are filled with less water but have excessive air flow. To investigate the impact of two-phase flow maldistribution on PEM fuel cell performance, a Volume of Fluid (VOF) model coupled with a 1D MEA model was employed to simulate two parallel channels. The slug flow pattern is mainly observed in the flow channels under different flow maldistribution conditions, and it significantly increases the gas diffusion layer (GDL) surface water coverage over the whole range of simulated current densities, which directly leads to poor fuel cell performance. Therefore, it is recommended that liquid and gas flow maldistribution in parallel channels should be avoided if possible over the whole range of operation. Increasing the gas stoichiometric flow ratio is not an effective method to mitigate the gas flow maldistribution, but adding a gas inlet resistance to the flow channel is effective in mitigating maldistribution. With a carefully selected value of the flow resistance coefficient, both the fuel cell performance and the gas flow distribution can be significantly improved without causing too much extra pressure drop.     

Numerical simulation of two-phase flow in a multi-gas channel of a proton exchange membrane fuel cell

Understanding the two-phase distribution characteristics within the multi-gas channel of a fuel cell is important for improving fuel cell performance. In the paper, the volume of fluid model is used to predict the dynamic behaviour of water in the multi-gas channel, analyze the pressure drop, velocity distribution, and flow resistance coefficient between different channels, and investigate the influence of operating conditions, surface wettability and channel structure on the two-phase distribution characteristics in the channel. The results show that water undergoes the processes of growth, separation, single droplet transport, wall impact, droplet collision, liquid film formation, and liquid film transport in the multi-gas channel. Inlet velocity and surface wettability significantly affect the pressure drop, water saturation, and surface water coverage. As the inlet velocity and gas diffusion layer surface wettability increase, the flow resistance coefficient and unevenness of the distribution decrease, indicating that the in-channel flow distribution homogeneity is enhanced. The rectangular channel has better water removal and flow distribution uniformity than the tapered channel, and the unevenness of distribution decreases significantly with decreasing rectangular width, from 0.15715 to 0.00315. The research work is a guide to understanding water transport in multi-gas channels, accelerating water removal, and improving inter-channel flow distribution uniformity.     

Investigation of liquid water heterogeneities in large area proton exchange membrane fuel cells using a Darcy two-phase flow model in a multiphysics code

Proper management of the liquid water and heat produced in proton exchange membrane (PEM) fuel cells remains crucial to increase both its performance and durability. In this study, a two-phase flow and multicomponent model, called two-fluid model, is developed in the commercial COMSOL Multiphysics® software to investigate the liquid water heterogeneities in large area PEM fuel cells, considering the real flow fields in the bipolar plate. A macroscopic pseudo-3D multi-layers approach has been chosen and generalized Darcy's relation is used both in the membrane-electrode assembly (MEA) and in the channel. The model considers two-phase flow and gas convection and diffusion coupled with electrochemistry and water transport through the membrane. The numerical results are compared to one-fluid model results and liquid water measurements obtained by neutron imaging for several operating conditions. Finally, according to the good agreement between the two-fluid and experimentation results, the numerical water distribution is examined in each component of the cell, exhibiting very heterogeneous water thickness over the cell surface.     

Modelling and optimization of reactant gas transport in a PEM fuel cell with a transverse pin fin insert in channel flow

A proton exchange membrane (PEM) fuel cell has many distinctive features which make it an attractive alternative clean energy source. Some of those features are low start-up, high power density, high efficiency and remote applications. In the present study, a numerical investigation was conducted to analyse the flow field and reactant gas distribution in a PEM fuel cell channel with transversely inserted pin fins in the channel flow aimed at improving reactant gas distribution. A fin configuration of small hydraulic diameter was employed to minimise the additional pressure drop. The influence of the pin fin parameters, the flow Reynolds number, the gas diffusion layer (GDL) porosity on the reactant gas transport and the pressure drop across the channel length were explored. The parameters examined were optimized using a mathematical optimization code integrated with a commercial computational fluid dynamics code. The results obtained indicate that a pin fin insert in the channel flow considerably improves fuel cell performance and that optimal pin fin geometries exist for minimized pressure drop along the fuel channel for the fuel cell model considered. The results obtained provide a novel approach for improving the design of fuel cells for optimal performance.     

Voltage loss and fluctuation in proton exchange membrane fuel cells: The role of cathode channel plurality and air stoichiometric ratio

Water management is a critical issue in the development of proton exchange membrane (PEM) fuel cells with robust operation. Liquid water can accumulate and flood the gas delivery microchannels and the porous electrodes within PEM fuel cells and deteriorate performance. Since the liquid distribution fluctuates in time for two-phase flow, the rate of oxygen transport to the cathode catalyst layer also fluctuates, resulting in unstable power density and efficiency. This paper reports experimental data on the mean voltage loss and the voltage fluctuations during constant current operation as a function of both the number of parallel microchannels and the air flow rate stoichiometric ratio. We define channel plurality as a flow field design parameter to describe the number of channels per unit of active area. The voltage loss was found to scale proportionally to channel plurality divided by the air stoichiometric ratio. The amplitude of the voltage fluctuations was found to be linearly proportional to channel plurality and inversely proportional to the air stoichiometric ratio squared. By analyzing pressure drop data and power spectra, we conclude that the voltage fluctuations are well-correlated with two-phase flow instabilities in the cathode's parallel microchannels. Finally, a scaling analysis is presented for generalizing the results for fuel cells having different active area and channel cross-section.     

Numerical study of a new cathode flow-field design with a sub-channel for a parallel flow-field polymer electrolyte membrane fuel cell

The cathode flow-field design of a polymer electrolyte membrane (PEM) fuel cell is crucial to its performance, because it determines the distribution of reactants and the removal of liquid water from the fuel cell. In this study, the cathode flow-field of a parallel flow-field PEM fuel cell was optimized using a sub-channel. The main-channel was fed with moist air, whereas the sub-channel was fed with dry air. The influences of the sub-channel flow rate (SFR, the amount of air from the sub-channel inlet as a percentage of the total cathode flow rate) and the inlet positions (SIP, where the sub-channel inlets were placed along the cathode channel) on fuel cell performance were numerically evaluated using a three-dimensional, two-phase fuel cell model. The results indicated that the SFR and SIP had significant impacts on the distribution of the feed air, removal of liquid water, and fuel cell performance. It was found that when the SIP was located at about 30% along the length of the channel from main-channel inlet and the SFR was about 70%, the PEM fuel cell exhibited much better performance than seen with a conventional design.