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.     

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.     

Ex situ and modeling study of two-phase flow in a single channel of polymer electrolyte membrane fuel cells

In this study, we investigate the air-water two-phase flow in a single flow channel of polymer electrolyte membrane (PEM) fuel cells. In the ex situ study, both straight and serpentine channels with various gas diffusion layer (GDL) surfaces are studied. Focus is placed on the two-phase flow patterns, which are optically characterized using a microscope with a high-resolution camera, and the two-phase pressure amplifiers. We find that the GDL surface properties slightly affect the flow pattern and two-phase pressure amplifier in the flow field configuration. Flow pattern transition occurs at the superficial gas velocity of around 1 m s⁻¹, and the pressure amplifier can reach as high as 10. A two-fluid model is also presented together with one dimensional (1-D) analytical solution, and acceptable agreement is achieved between the model prediction and experimental data at high gas flow rates.     

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.     

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.     

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.     

A PEM fuel cell model for cold-start simulations

In this paper, a transient multiphase multi-dimensional PEM fuel cell model has been developed in the mixed-domain framework for elucidating the fundamental physics of fuel cell cold start. Cold-start operations of a PEM fuel cell at a subfreezing boundary temperature of −20 °C under both constant current and constant cell voltage conditions have been numerically examined. Numerical results indicate that the water vapor concentration inside the cathode gas channel affects ice formation in the cathode catalyst layer and thus the cold-start process of the fuel cell. This conclusion demonstrates that high gas flow rates in the cathode gas channel could increase fuel cell cold-start time and benefit the cold-start process. It is shown that the membrane plays a significant role during the cold-start process of a PEM fuel cell by absorbing the product water and becoming hydrated. The time evolutions of ice formation, current density and water content distributions during fuel cell cold-start processes have also been discussed in detail.     

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.     

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.     

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.     

An inverse geometry design problem for optimization of single serpentine flow field of PEM fuel cell

The cathode flow-field design of a proton exchange membrane fuel cell (PEMFC) determines its reactant transport rates to the catalyst layer and removal rates of liquid water from the cell. This study optimizes the cathode flow field for a single serpentine PEM fuel cell with 5 channels using the heights of channels 2–5 as search parameters. This work describes an optimization approach that integrates the simplified conjugated-gradient scheme and a three-dimensional, two-phase, non-isothermal fuel cell model. The proposed optimal serpentine design, which is composed of three tapered channels (channels 2–4) and a final diverging channel (channel 5), increases cell output power by 11.9% over that of a cell with straight channels. These tapered channels enhance main channel flow and sub-rib convection, both increasing the local oxygen transport rate and, hence, local electrical current density. A diverging, final channel is preferred, conversely, to minimize reactant leakage to the outlet. The proposed combined approach is effective in optimizing the cathode flow-field design for a single serpentine PEMFC. The role of sub-rib convection on cell performance is demonstrated.     

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.     

Innovative gas diffusion layers and their water removal characteristics in PEM fuel cell cathode

Liquid water transport is one of the key challenges regarding the water management in a proton exchange membrane (PEM) fuel cell. Conventional gas diffusion layers (GDLs) do not allow a well-organized liquid water flow from catalyst layer to gas flow channels. In this paper, three innovative GDLs with different micro-flow channels were proposed to solve liquid water flooding problems that conventional GDLs have. This paper also presents numerical investigations of air–water flow across the proposed innovative GDLs together with a serpentine gas flow channel on PEM fuel cell cathode by use of a commercial computational fluid dynamics (CFD) software package FLUENT. The results showed that different designs of GDLs will affect the liquid water flow patterns significantly, thus influencing the performance of PEM fuel cells. The detailed flow patterns of liquid water were shown. Several gas flow problems for the proposed different kinds of innovative GDLs were observed, and some useful suggestions were given through investigating the flow patterns inside the proposed GDLs.     

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.     

Water flooding and two-phase flow in cathode channels of proton exchange membrane fuel cells

The water flooding and two-phase flow of reactants and products in cathode flow channels (0.8 mm in width, 1.0 mm in depth) were studied by means of transparent proton exchange membrane fuel cells. Three transparent proton exchange membrane fuel cells with different flow fields including parallel flow field, interdigitated flow field and cascade flow field were used. The effects of flow field, cell temperature, cathode gas flow rate and operation time on water build-up and cell performance were studied, respectively. Experimental results indicate that the liquid water columns accumulating in the cathode flow channels can reduce the effective electrochemical reaction area; it makes mass transfer limitation resulting in the cell performance loss. The water in flow channels at high temperature is much less than that at low temperature. When the water flooding appears, increasing cathode flow rate can remove excess water and lead to good cell performance. The water and gas transfer can be enhanced and the water removal is easier in the interdigitated channels and cascade channels than in the parallel channels. The cell performances of the fuel cells that installed interdigitated flow field or cascade flow field are better than that installed with parallel flow field. The images of liquid water in the cathode channels at different operating time were recorded. The evolution of liquid water removing out of channels was also recorded by high-speed video.     

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.     

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.     

Microscale and Macroscale Aspects of Water Management Challenges in PEM Fuel Cells

Water management is critical to the successful implementation of proton exchange membrane (PEM) fuel cells in the automotive transportation sector. Liquid water appears in the fuel cells not only from the water generated at the cathode catalyst layer but also as a result of condensation of water vapor within the humidified gases. Topics of intense current interest include the microscopic flow of liquid water through the membrane, catalyst layers, and gas diffusion layers; the macroscopic interaction between water and the gas flow at the gas diffusion layer interface; and the two-phase, multicomponent flow through the gas channels. Recent work published in this area is reviewed and recommendations for future work are outlined.     

Application of a thermally conductive pyrolytic graphite sheet to thermal management of a PEM fuel cell

This work experimentally investigates the thermal performance of a pyrolytic graphite sheet (PGS) in a single proton exchange membrane fuel cell (PEMFC). This PGS with high thermal conductivity serves as a heat spreader, reduces the volume and weight of cooling systems, and reduces and homogenizes the temperature in the reaction area of the fuel cells. A transparent PEMFC is constructed with PGS of thickness 0.1 mm cut into the shape of a flow channel and bound with the cathode gas channel plate. Eleven thermocouples are embedded at different positions on the cathode gas channel plate to measure the temperature distribution. The water and water flooding inside the cathode gas channels, with and without PGS, were successfully visualized. The locations of liquid water are correlated with the temperature measurement. PGS reduces the maximum cell temperature and improves cell performance at high cathode flow rates. The temperature distribution is also more uniform in the cell with PGS than in the one without PGS. Results of this study demonstrate the promising application of PGS to the thermal management of a fuel cell system.     

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.