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.     

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.     

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

In the present study, gas–liquid two-phase flow in a parallel square minichannel system oriented horizontally and at an incline is studied under operating conditions relevant to fuel cell operations. Flow mal-distribution in parallel channels occurs at low gas and liquid flow rates. In general, high superficial gas velocities are required to ensure even flow distribution, and the minimum gas flow rates required to achieve even distribution depend on the liquid flow rates, channel orientation and experimental procedures. As the inclination angle is increased, a higher gas flow rate is required to ensure even gas–liquid flow distribution while flow channels inclined downward seems to help in improving the even flow distribution. The presence of flow hysteresis phenomena indicate that multiple flow distributions exist at the same given flow conditions when the gas flow rates are varied in ascending and descending manners. Flow mal-distribution and flow hysteresis are directly linked with flow stability. More specifically, the actual gas and liquid distribution in parallel channels is determined by the stability of mathematical solutions of mass and momentum balance equations and also the flow history. For the first time, the present work investigates flow distributions in fuel cell flow fields by accounting for two-phase flow conditions. In addition, a novel approach is introduced to ensure flow distributions and their stability through contour construction of isobars where unstable flow region can be identified, which can be used in the design of parallel channel flow fields, especially for fuel cells.     

Numerical study of two-phase flow patterns in the gas channel of PEM fuel cells with tapered flow field design

In this study the air–water two-phase flow in a tapered channel of a PEMFC was numerically simulated using the volume of fluid (VOF) method. In particular, a 3D mathematical model of the fuel cell flow channel was used to obtain a reliable evaluation of the fuel cell performance for different taper angles and different temperatures and to calculate the total amount of water produced. This information was then used as boundary conditions to simulate the two-phase flow in the cell channel through a 2D VOF model. Typical operating conditions were assigned and the numerical mesh was constructed to represent the real fuel cell configuration. The results show that tapering the channel downstream enhances the water removal due to increased airflow velocity. In the rectangular channel no film formation is noted with a marked predominance of slug flow. In contrast, as the taper angle is increased the predominant two-phase flow pattern is film flow. Finally many contact angles have been used to simulate the effect of the hydrophobicity of a GDL surface on the motion of the water. As the hydrophobicity of a GDL surface is decreased the presence of film is more evident even for less tapered channels.     

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.     

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.     

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.     

Experimental analysis of the unit cell approach for two-phase flow dynamics in curved flow channels

Flow behavior of gas–liquid mixtures in thin channels has become increasingly important as a result of miniaturization of fluid and thermal systems. The present empirical study investigates the use of the unit cell or periodic boundary approach commonly used in two-phase flows. This work examines the flow patterns formed in small tube diameter (<3 mm) and curved geometry flow systems for air–water mixtures at standard conditions. Liquid and gas superficial velocities were varied from 0.1 to 7.0 (～±0.01) m/s and 0.03 to 14 (～±0.2) m/s for air and water respectively to determine the flow pattern formed in three geometries and dispersed bubble, plug, slug and annular flow patterns are reported using high-frame rate videography. Flow patterns formed were plotted on the generalized two-phase flow pattern map to interpret the effect of channel size and curvature on the flow regime boundaries. Relative to a straight a channel, it is shown that a ‘C shaped’ channel that causes a directional change in the flow induces chaotic advection and increases phase interaction to enhance gas bubble or liquid slug break-up thus altering the boundaries between the dispersed bubble and plug/slug flow regimes as well as between the annular and plug/slug flow regimes.     

Gas–liquid two-phase flow patterns in a miniature square channel with a gas permeable sidewall

Characteristics of air–water two-phase flow patterns in a miniature square channel having a gas permeable sidewall were investigated experimentally. Water was fed into the channel from its entrance, while air was injected uniformly into the channel along the permeable sidewall. This configured two-phase flow problem is encountered in direct feed methanol fuel cells. Flow patterns in both vertical upward and horizontal flows were identified using a high-speed motion analyzer. The visualization shows that the typical flow pattern encountered in the conventional co-current gas–liquid two-phase flow, such as bubbly flow, plug flow, slug flow and annular flow were also observed in the present work. However, unlike the conventional co-current gas–liquid two-phase flow in a channel with gas and liquid uniformly entering from one of its ends, for the flow configuration considered in this work, the stratified flow and wavy flow were not found in horizontal flow. And a so-called “single layer bubbly flow” was found in vertical upward flow, which is characterized by a mono small-gas-bubble layer existing adjacent to the surface of the permeable sidewall with the reminding space occupied by the liquid phase. Four transitional flow patterns such as bubbly-plug flow, bubbly-slug flow, plug–slug flow, and slug-annular flow, were found to exist between the distinct flow patterns. Finally, the flow regime maps for various liquid volumetric fluxes are presented in terms of mass quality versus the volumetric flux of gas phase.     

Effect of flow regime of circulating water on a proton exchange membrane electrolyzer

The flow characteristics of circulating water in a proton exchange membrane (PEM) electrolyzer were experimentally evaluated using a small cell and two-phase flow theory. Results revealed that when a two-phase flow of circulating water at the anode is either slug or annular, then mass transport of the water for the anode reaction is degraded, and that the concentration overvoltage increases at higher current density compared to that when the flow is bubbly. In a serpentine-dual flow field, when both phases of the two-phase flow are assumed laminar, then the increase in pressure drop caused by the increase in gas production can be explained relatively well using the Lockhart–Martinelli method with the Chisholm parameter. The optimal flow rate of circulating water was also discussed based on mass balance analysis.     

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.     

Effects of electrode wettabilities on liquid water behaviours in PEM fuel cell cathode

Liquid water transport is one of the key challenges for water management in a proton exchange membrane (PEM) fuel cell. Investigation of the air–water flow patterns inside fuel cell gas flow channels with gas diffusion layer (GDL) would provide valuable information that could be used in fuel cell design and optimization. This paper presents numerical investigations of air–water flow across an innovative GDL with catalyst layer and serpentine channel on PEM fuel cell cathode by use of a commercial Computational Fluid Dynamics (CFD) software package FLUENT. Different static contact angles (hydrophilic or hydrophobic) were applied to the electrode (GDL and catalyst layer). The results showed that different wettabilities of cathode electrode could affect liquid water flow patterns significantly, thus influencing on the performance of PEM fuel cells. The detailed flow patterns of liquid water were shown, several gas flow problems were observed, and some useful suggestions were given through investigating the flow patterns.     

Effects of temperature on the location of the gas–liquid interface in a PEM fuel cell

The objective of this study is to investigate the location of the gas–liquid interface at various temperatures in a polymer electrolyte membrane fuel cell under non-isothermal conditions. A mathematical model, coupled with the electrochemical process, two-phase flows, species transfer, and heat transfer is employed. A finite volume-based CFD approach is applied to investigate the species transport behavior in a fuel cell. The effects of two model parameters, namely cell temperature (T_cell) and humidification temperature (T_h), on the gas–liquid interface and cell performance are presented. Simulation results indicate that variations of these two parameters influence the location of the gas–liquid interface, the cell performance, and the distribution of liquid water saturation. At lower cell temperatures, the gas–liquid interface moves toward the inlet port of the channel when the humidification temperature is greater than the cell temperature. Therefore, the cell performance decreases as the liquid water clogs the passage for the transport of oxygen. Furthermore, these two factors are closely related to the membrane temperature distribution. Obvious variations in magnitude are seen at a cell temperature of 323 K and a humidification temperature of 343 K.     

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.     

A hybrid model of cathode of PEM fuel cell using the interdigitated gas distributor

A two-dimensional (2D), single- and two-phase, hybrid multi-component transport model is developed for the cathode of PEM fuel cell using interdigitated gas distributor. The continuity equation and Darcy's law are used to describe the flow of the reactant gas and production water. The production water is treated as vapor when the current density is small, and as two-phase while the current density is greater than the critical current density. The advection–diffusion equations are utilized to study species transport of multi-component mixture gas. The Butler–Volmer equation is prescribed for the domain in the catalyst layer. The predicted results of the hybrid model agree well with the available experimental data. The model is used to investigate the effects of operating conditions and the cathode structure parameters on the performance of the PEM fuel cell. It is observed that liquid water appears originally in the cathodic catalyst layer over outlet channel under intermediate current and tends to be distributed uniformly by the capillary force with the increase of the current. It is found that reduction of the width of outlet channel can enhance the performance of PEM fuel cell via the increase of the current density over this region, which has, seemingly, not been discussed in previous literatures.     

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.     

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.     

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.     

Numerical study of bubble generation and transport in a serpentine channel with a T-junction

The two-phase flow problem is one of the most critical challenges for direct methanol fuel cells (DMFC). With carbon dioxide generated in the anode catalyst layer, a bubble that could not be removed from the flow channel would hinder the oxidation reaction. In this paper, bubble generation and flow in a micro-serpentine channel is modeled using a multi-phase three-dimensional Navier–Stokes plus volume of fluid (VOF) method. The bubble generation process, which results from the gas–liquid flow and the effects of surface tension and viscosity, are discussed and validated by a comparison with a serious of experimental results. The results prove that the VOF method is an effective method to simulate flow fields where the bubble flow phenomena exist, such as the two-phase flow fields in DMFCs.     

Effect of air flow on liquid water transport through a hydrophobic gas diffusion layer of a polymer electrolyte membrane fuel cell

The transport of liquid water through an idealized 2-D reconstructed gas diffusion layer (GDL) of a polymer electrolyte membrane (PEM) fuel cell is computed subject to hydrophobic boundary condition at the fibre–fluid interface. The effect of air flow, as would occur in parallel/serpentine/interdigitated type of flow fields, on the liquid water transport through the GDL, ejection into the channel in the form of water droplets and subsequent removal of the droplets has been simulated. Results show that typically water flow through the fibrous GDL occurs through a fingering and channelling type of mechanism. The presence of cross-flow of air has an effect both on the path created within the GDL and on the ejection of water into the channel in the form of droplets. A faster rate of liquid water evacuation through the GDL (i.e., more frequent ejection of water droplets) as well as less flooding of the void space results from the presence of cross-flow. These results agree qualitatively with experimental observations reported in the literature.