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.     

Modeling two-phase flow in PEM fuel cell channels

This paper is concerned with the simultaneous flow of liquid water and gaseous reactants in mini-channels of a proton exchange membrane (PEM) fuel cell. Envisaging the mini-channels as structured and ordered porous media, we develop a continuum model of two-phase channel flow based on two-phase Darcy's law and the M² formalism, which allow estimate of the parameters key to fuel cell operation such as overall pressure drop and liquid saturation profiles along the axial flow direction. Analytical solutions of liquid water saturation and species concentrations along the channel are derived to explore the dependences of these physical variables vital to cell performance on operating parameters such as flow stoichiometric ratio and relative humility. The two-phase channel model is further implemented for three-dimensional numerical simulations of two-phase, multi-component transport in a single fuel-cell channel. Three issues critical to optimizing channel design and mitigating channel flooding in PEM fuel cells are fully discussed: liquid water buildup towards the fuel cell outlet, saturation spike in the vicinity of flow cross-sectional heterogeneity, and two-phase pressure drop. Both the two-phase model and analytical solutions presented in this paper may be applicable to more general two-phase flow phenomena through mini- and micro-channels.     

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.     

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.     

An Experimental Investigation into the Effect of Surfactants on Air-Water Two-Phase Flow in Minichannels

The complex interfacial phenomena involved in two-phase gas-liquid flow have defied mathematical simplification and modeling. However, these systems are used in heat exchangers, condensers, chemical processing plants, nuclear reactor systems, and fuel cells. The present work considers a 1 mm-square minichannel and adiabatic flows corresponding to practical PEM fuel cell conditions. Pressure drop data is collected over mass fluxes of 4.0–12.0 kg/m²s for air and 0.5–21.6 kg/m²s for water, corresponding to superficial gas and liquid velocities of 3.19–10.06 m/s and 0.0005–0.022 m/s, respectively. The experiments are repeated with water-surfactant mixtures of different concentrations in order to quantify the surface tension effects, as it is recognized that surface tension is an important parameter for two-phase flow in minichannels. The accuracy of various two-phase pressure drop models is evaluated, and a new model for laminar-laminar two-phase flow pressure drop is developed.     

Gas–liquid two-phase flow patterns in parallel channels for fuel cells

Two-phase flow in horizontal parallel channels has been experimentally investigated under fuel cell related operating conditions. Pronounced hysteresis is observed in the pressure drop versus flow characteristic curve when starting from either flooded or dry conditions. When gas is introduced into channels initially filled with water (flooded initial condition), both gas and liquid tend to flow predominantly in one channel at low gas or liquid flow velocities. As the gas flow velocity increases, even distribution of gas and liquid flow in both channels is observed, accompanied with a sudden decrease in the pressure drop. On the other hand, even gas and liquid flow distribution between both channels is found at comparatively lower gas flow velocities when starting with dry-gas flow conditions with liquid introduced into channels filled with gas (stratified flow regime). The flow regimes of this system are visualized in plots of the pressure drop against gas and liquid flow velocities. However, this phenomenon tends to vanish at high gas and liquid flow velocities, suggesting that high gas and liquid flow velocities are required to ensure even flow distribution in parallel channels. The hysteresis points appear at the same level of the pressure drop, reflecting intrinsic characteristics of the parallel channels used in this study. These results have important implications for PEM fuel cell operational strategies. In order to avoid reactant mal-distribution in parallel flow channels in the flow field in the two-phase flow regime, fuel cells should be operated at sufficiently high gas flow velocities.     

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.     

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.     

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.     

Measurement and modeling of liquid film thickness evolution in stratified two-phase microchannel flows

Polymer electrolyte membrane (PEM) fuel cells incorporating microchannels (D < 500 μm) can benefit from improved fuel delivery and convective cooling. However, this requires a better understanding of two-phase microchannel transport phenomena, particularly liquid–gas interactions and liquid clogging in cathode air-delivery channels. This paper develops optical fluorescence imaging of water films in hydrophilic channels with varying air velocity and water injection rate. Micromachined silicon test structures with optical access and distributed water injection simulate the cathode channels of a PEM fuel cell. Film thickness data vary strongly with air velocity and are consistent with stratified flow modeling. This work facilitates the study of regime transitions in two-phase microchannel flows and the effects of flow regimes on heat and mass transfer and axial pressure gradients.     

Water management studies in PEM fuel cells,Part II: Ex situ investigation of flow maldistribution,pressure drop and two-phase flow pattern in gas channels

Two-phase flow of water and reactant gases in the gas distribution channels of proton exchange membrane fuel cells (PEMFCs) plays a critical role in proper water management. In this work, the two-phase flow in PEMFC cathode parallel channels is studied over a wide range of superficial air velocity (air stoichiometry) and superficial water velocity in a specially designed ex situ experimental setup, which enables the measurement of instantaneous flow rates in individual gas channels and simultaneous visualization of the water flow structure. It is found that the two-phase flow at low superficial air velocities (air stoichiometry below 5) is dominated by slugs or semi-slugs, leading to severe flow maldistribution and large fluctuations in the pressure drop. Slug residence time, measured from the video observation and the instantaneous flow rate data, is found to be a new parameter to describe the slug flow. At higher air velocities, a water film is formed on the channel walls if they are hydrophilic. The pressure drop for the film flow is characterized by smaller but frequent fluctuations, which are found to result from the water buildup at the channel-exit manifold interface. As the superficial air velocity increases further, mist flow is obtained where little water buildup is observed. The water buildup in the gas channels at the two-phase flow is well described by the two-phase friction multiplier, defined as the ratio of the two-phase pressure drop to the single gas phase pressure drop. It is found that the two-phase friction multiplier increases with increasing water flow rate. A flow pattern map is developed using superficial water and air velocities with clearly defined transition regions.     

Two-phase flow and maldistribution in gas channels of a polymer electrolyte fuel cell

Liquid water transport in a polymer electrolyte fuel cell (PEFC) is a major issue for automotive applications. Mist flow with tiny droplets suspended in gas has been commonly assumed for channel flow while two-phase flow has been modeled in other cell components. However, experimental studies have found that two-phase flow in the channels has a profound effect on PEFC performance, stability and durability. Therefore, a complete two-phase flow model is developed in this work for PEFC including two-phase flow in both anode and cathode channels. The model is validated against experimental data of the wetted area ratio and pressure drop in the cathode side. Due to the intrusion of soft gas diffusion layer (GDL) material in the channels, flow resistance is higher in some channels than in others. The resulting flow maldistribution among PEFC channels is of great concern because non-uniform distributions of fuel and oxidizer result in non-uniform reaction rates and thus adversely affect PEFC performance and durability. The two-phase flow maldistribution among the parallel channels in an operating PEFC is explored in detail.     

Water management studies in PEM fuel cells, part IV: Effects of channel surface wettability, geometry and orientation on the two-phase flow in parallel gas channels

In this study, the effects of channel surface wettability, cross-sectional geometry and orientation on the two-phase flow in parallel gas channels of proton exchange membrane fuel cells (PEMFCs) are investigated. Ex situ experiments were conducted in flow channels with three different surface wettability (hydrophilically coated, uncoated, and hydrophobically coated), three cross-sectional geometries (rectangular, sinusoidal and trapezoidal), and two orientations (vertical and horizontal). Flow pattern map, individual channel flow variation due to maldistribution, pressure drop and flow visualization images were used to analyze the two-phase flow characteristics. It is found that hydrophilically coated gas channels are advantageous over uncoated or slightly hydrophobic channels regarding uniform water and gas flow distribution and favoring film flow, the most desirable two-phase flow pattern in PEMFC gas channels. Sinusoidal channels favor film flow and have lower pressure drop than rectangular and trapezoidal channels, while the rectangular and trapezoidal channels behave similarly to each other. Vertical channel orientation is advantageous over horizontal orientation because the latter is more prone to slug flow, nonuniform liquid water distribution and instable operation.     

Dominant frequency of pressure drop signal as a novel diagnostic tool for the water removal in proton exchange membrane fuel cell flow channel

In this work, a transparent assembly was self-designed and manufactured to perform ex situ experimental study on the liquid water removal characteristics in PEM fuel cell parallel flow channels. It was found that the dominant frequency of the pressure drop across the flow channels may be utilized as an effective diagnostic tool for water removal. Peaks higher than 1 Hz in dominant frequency profile indicated water droplet removals at the outlet, whereas relatively lower peaks (between 0.3 and 0.8 Hz) corresponded to water stream removals. The pressure drop signal, although correlated with the water removal at the outlet, was readily influenced by the two phase flow transport in channel, particularly at high air flow rates. The real-time visualization images were presented to show a typical water droplet removal process. The findings suggest that dominant frequency of pressure drop signal may substitute pressure drop as a more effective and reliable diagnostic tool for water removal in PEM fuel cell flow channels.     

Liquid water transport in parallel serpentine channels with manifolds on cathode side of a PEM fuel cell stack

Water management in a proton exchange membrane (PEM) fuel cell stack has been a challenging issue on the road to commercialization. This paper presents a numerical investigation of air–water flow in parallel serpentine channels on cathode side of a PEM fuel cell stack by use of the commercial Computational Fluid Dynamics (CFD) software package FLUENT. Different air–water flow behaviours inside the serpentine flow channels with inlet and outlet manifolds were discussed. The results showed that there were significant variations of water distribution and pressure drop in different cells at different times. The “collecting-and-separating effect” due to the serpentine shape of the gas flow channels, the pressure drop change due to the water distribution inside the inlet and outlet manifolds were observed. Several gas flow problems of this type of parallel serpentine channels were identified and useful suggestions were given through investigating the flow patterns inside the channels and manifolds.     

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.     

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.     

Numerical analysis on performances of polymer electrolyte membrane fuel cells with various cathode flow channel geometries

This paper investigated numerically the effect of cathode channel shapes on the local transport characteristics and cell performance by using a three-dimensional, two-phase, and non-isothermal polymer electrolyte membrane (PEM) fuel cell model. The cells with triangle, trapezoid, and semicircle channels were examined using that with rectangular channel as comparison basis. At high operating voltages, the cells with various channel shapes would have similar performance. However, at low operating voltages, the fuel cell performance would follow: triangle > semicircle > trapezoid > rectangular channel. Analyses of the local transport phenomena in the cell indicate that triangle, trapezoid, and semicircle channel designs increase remarkably flow velocity of reactant, enhancing liquid water removal and oxygen utilization. Thus, these designs increase the limiting current density and improve the cell performance relative to rectangular channel design.     

A two-phase flow and transport model for the cathode of PEM fuel cells

A unified two-phase flow mixture model has been developed to describe the flow and transport in the cathode for PEM fuel cells. The boundary condition at the gas diffuser/catalyst layer interface couples the flow, transport, electrical potential and current density in the anode, cathode catalyst layer and membrane. Fuel cell performance predicted by this model is compared with experimental results and reasonable agreements are achieved. Typical two-phase flow distributions in the cathode gas diffuser and gas channel are presented. The main parameters influencing water transport across the membrane are also discussed. By studying the influences of water and thermal management on two-phase flow, it is found that two-phase flow characteristics in the cathode depend on the current density, operating temperature, and cathode and anode humidification temperatures.