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.     

Parameter sensitivity examination for a complete three-dimensional,two-phase,non-isothermal model of polymer electrolyte membrane fuel cell

Parameter sensitivity analysis is carried out for a complete three-dimensional, two-phase, non-isothermal model of polymer electrolyte membrane (PEM) fuel cell with a parallel flow field design. The model couples the two-phase flow of the multi-component reactants and liquid water, species transport, electrochemical reactions, proton and electron transport, and the electro-osmosis transport, back diffusion of water in the membrane, and energy transport. Twenty nine parameters, which are classified into the structural or transport parameters of porous layers (tortuosity, porosity, permeability, proton conductivity, electron conductivity, and thermal conductivity) as well as the electrochemical parameters (anodic and cathodic exchange current densities, anodic and cathodic transfer coefficients for anode and cathode reactions), are used to implement individual parameter investigation. The results show the parameters can be divided in to strongly sensitive, conditional sensitive and weak sensitive parameters according to its effect on the cell polarization curve. The optimization of parameters of cathode gas diffusion layer (GDL) and catalyst layer (CL) is more important to improve cell performance than that of anode GDL and CL because liquid water transport and removal affect significantly membrane hydration and reactant transport. Electrochemical parameters determine the activation potential and the slope of ohmic polarization hence these parameters can be used to fit experimental polarization curve more effectively than the other parameters.     

Two-phase flow in GDL and reactant channels of a proton exchange membrane fuel cell

Understanding the effect of two-phase flow in the components of proton exchange membrane fuel cells (PEMFCs) is crucial to water management and subsequently to their performance. The local water saturation in the gas diffusion layer (GDL) and reactant channels influences the hydration of the membrane which has a direct effect on the PEMFC performance. Mass transport resistance includes contributions from both the GDL and reactant channels, as well as the interface between the aforementioned components. Droplet–channel wall interaction, water area coverage ratio on the GDL, oxygen transport resistance at the GDL–channel interface, and two-phase pressure drop in the channels are interlinked. This study explores each factor individually and presents a comprehensive perspective on our current understanding of the two-phase transport characteristics in the PEMFC reactant channels.     

In situ investigation of water transport in an operating PEM fuel cell using neutron radiography: Part 2 – Transient water accumulation in an interdigitated cathode flow field

An interdigitated cathode flow field has been tested in situ with neutron radiography to measure the water transport through the porous gas diffusion layer in a PEM fuel cell. Constant current density to open circuit cycles were tested and the resulting liquid water accumulation and dissipation rates with in-plane water distributions are correlated to measured pressure differential between inlet and outlet gas streams. The effect of varying the reactant gas relative humidity on liquid water accumulation is also demonstrated. These results provide evidence that the reactant gas establishes a consistent in-plane transport path through the diffusion layer, leaving stagnant regions where liquid water accumulates. A simplified permeability model is presented and used to correlate the relative permeability to varying gas diffusion layer liquid water saturation levels.     

Modeling of water transport through the membrane electrode assembly for direct methanol fuel cells

In this work, a one-dimensional, isothermal two-phase mass transport model is developed to investigate the water transport through the membrane electrode assembly (MEA) for liquid-feed direct methanol fuel cells (DMFCs). The liquid (methanol–water solution) and gas (carbon dioxide gas, methanol vapor and water vapor) two-phase mass transport in the porous anode and cathode is formulated based on classical multiphase flow theory in porous media. In the anode and cathode catalyst layers, the simultaneous three-phase (liquid and vapor in pores as well as dissolved phase in the electrolyte) water transport is considered and the phase exchange of water is modeled with finite-rate interfacial exchanges between different phases. This model enables quantification of the water flux corresponding to each of the three water transport mechanisms through the membrane for DMFCs, such as diffusion, electro-osmotic drag, and convection. Hence, with this model, the effects of MEA design parameters on water crossover and cell performance under various operating conditions can be numerically investigated.     

Numerical investigation of liquid water distribution in the cathode side of proton exchange membrane fuel cell and its effects on cell performance

A three-dimensional unsteady two-phase model for the cathode side of proton exchange membrane fuel cell (PEMFC) consisting of gas diffusion layer (GDL) with hybrid structural model is developed to investigate liquid water behaviors under different operating and geometrical conditions and to quantitatively evaluate effects of liquid water distribution on reactant transport and current density distribution. Simulation results reveal that liquid water transport processes and distributions are significantly affected by inlet air velocity, wall wettability and water inlet position, which in turn play a prominent role on local reactant transport and cause considerable disturbances of the current density. Liquid water film spreading on the gas channel (GC) top wall is identified as the most desirable flow pattern in the GC based on overall evaluations of current density magnitude, uniformity of current density distribution and pressure drop in the GC. Modification to GDL structure is proposed to promote the formation of the desirable flow pattern.     

Humidity of reactant fuel on the cell performance of PEM fuel cell with baffle-blocked flow field designs

The objective of this work is to examine the effects of humidity of reactant fuel at the inlet on the detailed gas transport and cell performance of the PEM fuel cell with baffle-blocked flow field designs. It is expected that, due to the water management problem, the effects of inlet humidity of reactant fuel gases on both anode and cathode sides on the cell performance are considerable. In addition, the effects of baffle numbers on the detailed transport phenomena of the PEM fuel cell with baffle-blocked flow field are examined. Due to the blockage effects in the presence of the baffles, more fuel gas in the flow channel can be forced into the gas diffuser layer (GDL) and catalyst layer (CL) to enhance the chemical reactions and then augment the performance of the PEMFC systems. Effect of liquid water formation on the reactant gas transport is taken into account in the numerical modeling. Predictions show that the local transport of the reactant gas, the local current density generation and the cell performance can be enhanced by the presence of the baffles. Physical interpretation for the difference in the inlet relative humidity (RH) effects at high and low operating voltages is presented. Results reveal that, at low voltage conditions, the liquid water effect is especially significant and should be considered in the modeling. The cell performance can be enhanced at a higher inlet relative humidity, by which the occurrence of the mass transport loss can be delayed with the limiting current density raised considerably.     

Numerical investigation on the characteristics of mass transport and performance of PEMFC with baffle plates installed in the flow channel

A 3D numerical model of proton exchange membrane fuel cell (PEMFC) with the installation of baffle plates is developed. The majority of the conservation equations and physical parameters are implemented through the user defined functions (UDFs) in the FLUENT software. The characteristics of mass transport and performance of PEMFC are investigated. The results reveal that the baffle plate can enhance the mass transport efficiency and the performance of PEMFC. The baffle plate installed in the PEMFC flow channel increases the local gas velocity, which can promote the reactant gas transport and the liquid water removal in the porous electrode. As a result, the reactant gas concentration is larger in the porous electrode, which enhances the fuel cell performance for decreasing the over-potential of concentration. The fuel cell output power increases with the blockage ratio of the baffle plate. Considering the extra pumping power resulted from pressure loss caused by the baffle plate, the fuel cell with the blockage ratio of 0.8 is found to perform best in terms of the fuel cell net power generation. The fuel cell performance increases first with the baffle plate number, due to the better reactant distribution and water management, but decreases when the baffle plate number is too large, due to the excessive blockage for the reactant gas transport to the channel downstream. The PEMFC investigated with 5 baffle plates in the channel is found to be optimal. A channel design to achieve gradually increasing blockage ratios is also proposed, which exhibits better cell performance than the design with even blockage ratios.     

Numerical investigations on liquid water removal from the porous gas diffusion layer by reactant flow

Water removal from the gas diffusion layer (GDL) is crucial for the efficient operation of proton exchange membrane (PEM) fuel cell. Static pressure gradient caused by the fast reactant flow in the flow channel is one of the main mechanisms of water removal from GDL. Reactant can leak or cross directly to the neighboring channel via the porous GDL in the cells with serpentine flow channel and many of its modifications. Such cross flow plays an important role for the removal of liquid water accumulated in the GDL especially under land area. To investigate the characteristics of liquid water behavior in the GDL under pressure gradient, the fibrous porous structure of the carbon paper is modeled by three dimensional impermeable cylinders randomly distributed in the in-plane directions and unsteady two-phase simulations are conducted. It is shown that the permeability from the numerical model matches well the experimental measurements of the common GDLs in the literature. The contact angle and pressure gradient are the key parameters that determine the initiation and the process of liquid water transport in the GDL which is initially wet with stagnant liquid water. It has been observed that the larger contact angle results in faster water removal from the GDL. Numerical simulations are performed for a wide range of pressure gradient with different contact angles to determine the minimum pressure gradient that initiates the liquid water transport in the GDL. It is found that the amount of pressure gradient caused by the cross flow is sufficient and effective to get rid of the liquid water accumulated in the GDL. The simulation results are also compared with experimental data in literature showing a good agreement. The characteristics of liquid water discharging from the gas diffusion layer are also described.     

Cross-leakage flow between adjacent flow channels in PEM fuel cells

In polymer electrolyte membrane (PEM) fuel cells, serpentine flow channels are used conventionally for effective water removal. The reactant flows along the flow channel with pressure decrease due to the frictional and minor losses as well as the reactant depletion because of electrochemical reactions in the cells. Because of the short distance between the adjacent flow channels, often in the order of 1 mm or smaller, the pressure gradient between the adjacent flow channels is very large, driving part of reactant to flow through the porous electrode backing layer (or the so-called gas diffusion layer)—this cross-leakage flow between adjacent flow channels in PEM fuel cells has been largely ignored in previous studies. In this study, the effect of cross-flow in an electrode backing layer has been investigated numerically by considering bipolar plates with single-channel serpentine flow field for both the anode and cathode side. It is found that a significant amount of reactant gas flows through the porous electrode structure, due to the pressure difference, and enters the next flow channel, in addition to a portion entering the catalyst layer for reaction. Therefore, mixing occurs between the relatively high concentration reactant stream following the flow channel and the relatively low reactant concentration stream going through the electrode. It is observed that the cross-leakage flow influences the reactant concentration at the interface between the electrode and the catalyst layer, hence the distribution of reaction rate or current density generated. In practice, this cross-leakage flow in the cathode helps drive the liquid water out of the electrode structure for effective water management, partially responsible for the good PEM fuel cell performance using the serpentine flow 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.     

Numerical simulation of three-dimensional gas/liquid two-phase flow in a proton exchange membrane fuel cell

Investigation into the formation and transport of liquid water in proton exchange membrane fuel cells (PEMFCs) is the key to fuel cell water management. A three-dimensional gas/liquid two-phase flow and heat transfer model is developed based on the multiphase mixture theory. The reactant gas flow, diffusion, and chemical reaction as well as the liquid water transport and phase change process are modeled. Numerical simulations on liquid water distribution and its effects on the performance of a PEMFC are conducted. Results show that liquid water distributes mostly in the cathode, and predicted cell performance decreases quickly at high current density due to the obstruction of liquid water to oxygen diffusion. The simulation results agree well with experimental data. Translated from J Tsinghua Univ (Sci & Tech), 2006, 46(2): 252–256 [译自: 清华大学学报]     

Effects of metal foam properties on flow and water distribution in polymer electrolyte fuel cells (PEFCs)

Using a two-phase polymer electrolyte fuel cell (PEFC) model, we numerically investigated the influence of metal foam porous properties and wettability on key species and current distributions within a PEFC. Three-dimensional simulations were conducted under practical low humidity inlet hydrogen and air gases, and numerical comparisons were made for different metal foam design variables. These simulations were conducted to elucidate the detailed operating characteristics of PEFCs using metal foams as a flow distributor, and the simulation results showed that two-phase transport and the resulting flooding behavior in a PEFC are influenced by both the metal foam porous properties and the porous properties of an adjacent layer (e.g., the gas diffusion layer). This paper offers basic directions to design metal foams for optimal water management of PEFCs.     

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.     

An approach to modeling two-phase transport in the gas diffusion layer of a proton exchange membrane fuel cell

A series of analysis methods is proposed to simulate the liquid–gas two-phase and multi-component transport phenomena in the gas diffusion layer (GDL) of a proton exchange membrane fuel cell (PEMFC). These methods involve measuring and predicting the two-phase flow properties of a GDL, and simulating the two-phase multi-component transport in the GDL. The capillary pressure is measured by the porous diaphragm method and predicted by the pore network model. The relative permeability is measured by the steady-state method and predicted by a combination of the single-phase and the two-phase lattice Boltzmann method. And the simulation of the liquid–gas two-phase transport is done using the multi-phase mixture model. The methods are applied to a carbon-fiber paper GDL to identify the two-phase multi-component transport in the GDL.     

A generalized numerical model for liquid water in a proton exchange membrane fuel cell with interdigitated design

A new approach to numerical simulation of liquid water distribution in channels and porous media including gas diffusion layers (GDLs), catalyst layers, and the membrane of a proton exchange membrane fuel cell (PEMFC) was introduced in this study. The three-dimensional, PEMFC model with detailed thermo-electrochemistry, multi-species, and two-phase interactions. Explicit gas-liquid interface tracking was performed by using Computational Fluid Dynamics (CFD) software package FLUENT^® v6.2, with its User-Defined Functions (UDF) combined with volume-of-fluid (VOF) algorithm. The liquid water transport on a PEMFC with interdigitated design was investigated. The behavior of liquid water was understood by presenting the motion of liquid water droplet in the channels and the porous media at different time instants. The numerical results show that removal of liquid water strongly depends on the magnitude of the flow field. Due to the blockage of liquid water, the gas flow is unevenly distributed, the high pressure regions takes place at the locations where water liquid appears. In addition, mass transport of the species and the current density distribution is significantly degraded by the presence of liquid water.     

Novel serpentine-baffle flow field design for proton exchange membrane fuel cells

An appropriate flow field in the bipolar plates of a fuel cell can effectively enhance the reactant transport rates and liquid water removal efficiency, improving cell performance. This paper proposes a novel serpentine-baffle flow field (SBFF) design to improve the cell performance compared to that for a conventional serpentine flow field (SFF). A three-dimensional model is used to analyze the reactant and product transport and the electrochemical reactions in the cell. The results show that at high operating voltages, the conventional design and the baffled design have the same performance, because the electrochemical rate is low and only a small amount of oxygen is consumed, so the oxygen transport rates for both designs are sufficient to maintain the reaction rates. However, at low operating voltages, the baffled design shows better performance than the conventional design. Analyses of the local transport phenomena in the cell indicate that the baffled design induces larger pressure differences between adjacent flow channels over the entire electrode surface than does the conventional design, enhancing under-rib convection through the electrode porous layer. The under-rib convection increases the mass transport rates of the reactants and products to and from the catalyst layer and reduces the amount of liquid water trapped in the porous electrode. The baffled design increases the limiting current density and improves the cell performance relative to conventional design.     

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.     

Three-dimensional two-phase flow model of proton exchange membrane fuel cell with parallel gas distributors

A non-isothermal, steady-state, three-dimensional (3D), two-phase, multicomponent transport model is developed for proton exchange membrane (PEM) fuel cell with parallel gas distributors. A key feature of this work is that a detailed membrane model is developed for the liquid water transport with a two-mode water transfer condition, accounting for the non-equilibrium humidification of membrane with the replacement of an equilibrium assumption. Another key feature is that water transport processes inside electrodes are coupled and the balance of water flux is insured between anode and cathode during the modeling. The model is validated by the comparison of predicted cell polarization curve with experimental data. The simulation is performed for water vapor concentration field of reactant gases, water content distribution in the membrane, liquid water velocity field and liquid water saturation distribution inside the cathode. The net water flux and net water transport coefficient values are obtained at different current densities in this work, which are seldom discussed in other modeling works. The temperature distribution inside the cell is also simulated by this model.     

Liquid water transport in PEMFC cathode with symmetrical biomimetic flow field design based on Murray's law

Water management in the flow field as well as the flooding process in the gas diffusion and catalyst layers enormously influence proton exchange membrane fuel cells (PEMFCs) performance and reliability. Researchers have developed many different designs for flow channels that can be used to distribute fuel or oxidant in PEMFCs (proton exchange membrane fuel cells). Among these designs, novel biomimetic designs have captured special attentions from researchers due to their capability of distributing fluids effectively. This study presents an investigation of the liquid water transport within a porous layer and a symmetrical biomimetic flow field based on Murray's law. The volume of fluid (VOF) method is employed, and the dynamic contact angle (DCA) effects are also considered for better prediction of water distribution. The water transport process and water distribution inside the porous layer and flow field are obtained from the simulation results. Recommendations are given for this type of flow field design based on the behaviors of liquid water in the porous layer and flow field.