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.     

Three dimensional numerical simulation of gas-liquid two-phase flow patterns in a polymer-electrolyte membrane fuel cells gas flow channel

Water management in polymer-electrolyte membrane fuel cells (PEMFCs) has a major impact on fuel cell performance and durability. To investigate the two-phase flow patterns in PEMFC gas flow channels, the volume of fluid (VOF) method was employed to simulate the air-water flow in a 3D cuboid channel with a 1.0 mm × 1.0 mm square cross section and a 100 mm in length. The microstructure of gas diffusion layers (GDLs) was simplified by a number of representative opening pores on the 2D GDL surface. Water was injected from those pores to simulate water generation by the electrochemical reaction at the cathode side. Operating conditions and material properties were selected according to realistic fuel cell operating conditions. The water injection rate was also amplified 10 times, 100 times and 1000 times to study the flow pattern formation and transition in the channel. Simulation results show that, as the flow develops, the flow pattern evolves from corner droplet flow to top wall film flow, then annular flow, and finally slug flow. The total pressure drop increases exponentially with the increase in water volume fraction, which suggests that water accumulation should be avoided to reduce parasitic energy loss. The effect of material wettability was also studied by changing the contact angle of the GDL surface and channel walls, separately. It is shown that using a more hydrophobic GDL surface is helpful to expel water from the GDL surface, but increases the pressure drop. Using a more hydrophilic channel wall reduces the pressure drop, but increases the water residence time and water coverage of the GDL surface.     

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.     

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.     

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.     

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.     

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.     

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.     

Interface resolving two-phase flow simulations in gas channels relevant for polymer electrolyte fuel cells using the volume of fluid approach

With the increased concern about energy security, air pollution and global warming, the possibility of using polymer electrolyte fuel cells (PEFCs) in future sustainable and renewable energy systems has achieved considerable momentum. A computational fluid dynamic model describing a straight channel, relevant for water removal inside a PEFC, is devised. A volume of fluid (VOF) approach is employed to investigate the interface resolved two-phase flow behavior inside the gas channel including the gas diffusion layer (GDL) surface. From this study, it is clear that the impact on the two-phase flow pattern for different hydrophobic/hydrophilic characteristics, i.e., contact angles, at the walls and at the GDL surface is significant, compared to a situation where the walls and the interface are neither hydrophobic nor hydrophilic (i.e., 90° contact angle at the walls and also at the GDL surface). A location of the GDL surface liquid inlet in the middle of the gas channel gives droplet formation, while a location at the side of the channel gives corner flow with a convex surface shape (having hydrophilic walls and a hydrophobic GDL interface). Droplet formation only observed when the GDL surface liquid inlet is located in the middle of the channel. The droplet detachment location (along the main flow direction) and the shape of the droplet until detachment are strongly dependent on the size of the liquid inlet at the GDL surface. A smaller liquid inlet at the GDL surface (keeping the mass flow rates constant) gives smaller droplets.     

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

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

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.     

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.     

VOF modelling of gas–liquid flow in PEM water electrolysis cell micro-channels

In this study, the gas–liquid flow through an interdigitated anode flow field of a PEM water electrolysis cell (PEMEC) is analysed using a three-dimensional, transient, computational fluid dynamics (CFD) model. To account for two-phase flow, the volume of fluid (VOF) method in ANSYS Fluent 17.2 is used. The modelled geometry consists of the anode channels and the anode transport layer (ATL). To reduce the complexity of the phenomena governing PEMEC operation, the dependence upon electro-chemistry is disregarded. Instead, a fixed source of the gas is applied at the interface between the ATL and the catalyst layer. An important phenomenon that the model is able to capture is the gas–liquid contact angle on both the channel wall and ATL-channel interface. Particularly, the latter interface is crucial in capturing bubble entrainment into the channel. To validate the numerical simulation, photos taken of the gas–liquid flow in a transparent micro-channel, are qualitative compared against the simulation results. The experimental observations confirm the models prediction of long Taylor bubbles with small bubbles in between. From the simulation results, further intriguing details of the flow are revealed. From the bottom to the top of the outgoing channel, the film thickness gradually increases from zero to 200 μm. This increase in the film thickness is due to the particular superficial velocity field that develops in an interdigitated flow. Here both the superficial velocities change along the length of the channel. The model is capable of revealing effect of different bubble shapes/lengths in the outgoing channel. Shape and the sequence of the bubbles affect the water flow distribution in the ATL. The model presented in this work is the first step in the development of a comprehensive CFD model that comprises multiphase flow in porous media and micro-channel, electro-chemistry in catalyst layers, ion transport in membrane, hydrogen evolution, etc. The model can aid in the study of gas–liquid flow and its impact on the performance of a PEMEC.     

Numerical investigation of liquid water dynamics in wave‐like gas channels of PEMFCs

For the air feed in proton exchange membrane fuel cells (PEMFCs), the wave‐like gas channel (GC) shows obvious advantages over the straight GC because the former enhances collision of secondary flow and diffusion in the gas diffusion layer (GDL). However, it is prone to water flooding, which brings greater pressure drop, larger pressure oscillation, and blocking of reaction area. In the present study, numerical models of the water dynamic processes, including water droplets emerging from micropores on the GDL surface and removing through the GC, are established based on the volume of fluid (VOF) method. Water coverage ratio and pressure drop are calculated to evaluate the water flooding. The effects of the dimensional parameters of wave‐like GC and contact angle of channel walls on the water accumulation are studied. The emergence and removal of liquid water is a quasiperiodic and oscillating process. Multicycle simulations show that channel pressure drop increases linearly with greater growth rate than channel length. The equilibrium position of water droplet is strongly dependent on the relative wettability of the GDL and bipolar plate (BPP) surfaces. And the geometric parameters of GC have a significant impact on the pressure, water removal behavior and detachment time. Smaller bent angle brings bigger pressure drop, and larger cycle length is helpful for relieving the oscillation of pressure.     

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.     

Gas flow rate distributions in parallel minichannels for polymer electrolyte membrane fuel cells: Experiments and theoretical analysis

In the present work, instantaneous gas flow rates in each of two parallel channels of gas-liquid two-phase flow systems were investigated through measurements of the pressure drop across the entrance region. Liquid flow rates in two branches were pre-determined through liquid injection independently into each channel. Experiments were conducted in two different manners, i.e., the gas flow rate was varied in both ascending and descending paths. Flow hysteresis was observed in both gas flow rate distributions and the overall pressure drop of two-phase flow systems. Effects of liquid flow rates on gas flow distributions were examined experimentally. The presence of flow hysteresis was found to be associated with different flow patterns at different combinations of gas and liquid flow rates and flow instability conditions. A new and simple method was developed to predict gas flow distributions based on flow regime-specific pressure drop models for different experimental approaches and flow patterns. In particular, two different two-phase pressure drop models were used for slug flow and annular flow, separately. Good agreement was achieved between theoretical predictions and our experimental data. The developed new method can be potentially applied to predict gas flow distributions in parallel channels for fuel cells.     

Dynamics of liquid water in the anode flow channels of PEM fuel cells: A numerical parametric study

A 3D volume of fluid (VOF) model for an anode channel in a PEM fuel cell has been built. The effects of the initial position of the water droplet, its size as well as the wettability of the gas diffusion layer (GDL) are investigated under different operating conditions. It is found that the initial position of the relatively small water droplet in the channel has almost no effect on the pressure drop and the time taken for the liquid water to move out from the channel; however, such effects become more profound as the size of the water droplet increases. Also, when the droplet is placed at the side wall of the channel, then it develops into pockets of water that are mainly located at the upper corners of the channel, thus causing a smaller pressure drop compared to the cases in which the water droplet is placed either on the surface of the GDL or on the top wall of the channel. Furthermore, the hydrogen velocity is found to have a negligible effect on the dynamics of liquid water; however, the pressure drop and removal time are significantly influenced by the hydrogen velocity. Moreover, as the size of the water droplet increases, the pressure drop increases and the time required for the liquid water to move out of the channel decreases. Finally, the pressure drop in the channel decreases and the removal time of the liquid water increases as the contact angle of the GDL decreases.     

Three‐dimensional simulation of water droplet movement in PEM fuel cell flow channels with hydrophilic surfaces

Water management is one of the critical issues in proton exchange membrane fuel cells, and proper water management requires effective removal of liquid water generated in the cathode catalyst layer, typically in the form of droplets through cathode gas stream in the cathode flow channel. It has been reported that a hydrophilic channel sidewall with a hydrophobic membrane electrode assembly (MEA) surface would have less chance for water accumulation on the MEA surface. Therefore, a comprehensive study on the effect of surface wettability properties on water droplet movement in flow channels has been conducted numerically. In this study, the water droplet movements in a straight flow channel with a wide range of hydrophilic surface properties and effects of inlet air velocities are analyzed by using three‐dimensional computational fluid dynamics method coupled with the volume‐of‐fluid (VOF) method for liquid–gas interface tracking. The results show that the water droplet movement is greatly affected by the channel surface wettability and air flow conditions. With low contact angle, droplet motion is slow due to more liquid–wall contact area. With high air flow velocities, increasing the contact angle of the channel surface results in faster liquid water removal due to lesser liquid–wall contact area. Copyright © 2010 John Wiley & Sons, Ltd.     

Three‐dimensional optimisation of a fuel gas channel of a proton exchange membrane fuel cell for maximum current density

Proton exchange membrane (PEM) fuel cells operated with hydrogen and air offer promising alternative to conventional fossil fuel sources for transport and stationary applications because of its high efficiency, low‐temperature operation, high power density, fast start‐up and potable power for mobile application. Power levels derivable from this class of fuel cell depend on the operating parameters. In this study, a three‐dimensional numerical optimisation of the effect of operating and design parameters of PEM fuel cell performance was developed. The model computational domain includes an anode flow channel, membrane electrode assembly and a cathode flow channel. The continuity, momentum, energy and species conservation equations describing the flow and species transport of the gas mixture in the coupled gas channels and the electrodes were numerically solved using a computational fluid dynamics code. The effects of several key parameters, including channel geometries (width and depth), flow orientation and gas diffusion layer (GDL) porosity on performance and species distribution in a typical fuel cell system have been studied. Numerical results of the effect of flow rate and GDL porosity on the flow channel optimal configurations for PEM fuel cell are reported. Simulations were carried out ranging from 0.6 to 1.6 mm for channel width, 0.5 to 3.0 mm for channel depth and 0.1 to 0.7 for the GDL porosity. Results were evaluated at 0.3 V operating cell voltage of the PEM fuel cell. The optimisation results show that the optimum dimension values for channel depth and channel width are 2.0 and 1.2 mm, respectively. In addition, the results indicate that effective design of fuel gas channel in combination with the reactant species flow rate and GDL porosity enhances the performance of the fuel cell. The numerical results computed agree well with experimental data in the literature. Consequently, the results obtained provide useful information for improving the design of fuel cells. Copyright © 2011 John Wiley & Sons, Ltd.     

A two-fluid model for two-phase flow in PEMFCs

In this study, a two-fluid (TF) model is developed for two-phase flows in proton exchange membrane fuel cells (PEMFCs). The drag force and lift force between gas and liquid phase are considered in N-S equations. In addition, a simplified model is introduced to obtain the liquid water droplet detachment diameter on the gas diffusion layer (GDL)/channel interface which involves the properties of the GDL/channel interface (contact angle and surface tension). The TF model and the simplified model for the prediction of water droplet detachment diameter on GDL/channel interface are validated by the comparison between the experimental data and the model results, respectively. The effect of the properties of GDL/channel interface (contact angle and surface tension) on two-phase behavior in PEMFCs is investigated, The results show that a high contact angle and a low surface tension are advantageous for liquid water removal in the gas channel and the GDL even though a low surface tension will lead to a low capillary force in the GDL.