The role of the gas diffusion layer on slug formation in gas flow channels of fuel cells

Water drops emerge from large pores of the hydrophobic Gas Diffusion Layers (GDL) into the cathode gas flow channel of Polymer Electrolyte Membrane (PEM) Fuel Cells. The drops grow into slugs that span the cross-section of the flow channels. The slugs detach and are forced out the gas flow channel by the air flow. An acrylic micro-fluidic flow cell with a 1.6 mm gas flow channel and a 100 μm liquid pore through a carbon paper GDL has been used to quantitatively determine slug volumes, velocity of slug motion, and the force required to move slugs as functions of the gas and liquid flow rates. In a channel with 4 acrylic walls, slugs detach immediately upon formation. A porous GDL wall allows gas flow to bypass the slugs, thus allowing slugs to continue to grow after spanning the open area of the channel. The differential pressure to detach and move slugs is equal to the dynamic interfacial force on a slug normalized by the cross-sectional area of the channel. The dynamic interfacial force is equal to the difference between the downstream (advancing) and upstream (receding) contact lines of the water with the channel walls. Slugs will stop moving if the differential pressure drop for gas flow to bypass the slug and flow through the GDL under the rib separating the channels is less than the differential pressure required to move the slug. The results improve our physical insight into the state of water hold up in PEM fuel cells.     

Effects of needle orientation and gas velocity on water transport and removal in a modified PEMFC gas flow channel having a hydrophilic needle

Water management is critical to the performance and operation of the proton exchange membrane fuel cell (PEMFC). Effective water removal from the gas diffusion layer (GDL) surface exposed to the gas flow channel in PEMFC mitigates the water flooding of and improves the reactants transport into the GDL, hence benefiting the PEMFC performance. In this study, a 3D numerical investigation of water removal from the GDL surface in a modified PEMFC gas flow channel having a hydrophilic needle is carried out. The effects of the needle orientation (inclination angle) and gas velocity on the water transport and removal are investigated. The results show that the water is removed from the GDL surface in the channel for a large range of the needle inclination angle and gas velocity. The water is removed more effectively, and the pressure drop for the flow in the channel is smaller for a smaller needle inclination angle. It is also found that the modified channel is more effective and viable for water removal in fuel cells operated at smaller gas velocity.     

Two-phase flow pressure drop hysteresis in parallel channels of a proton exchange membrane fuel cell

Two-phase flow pressure drop hysteresis was studied in a non-operational PEM fuel cell to understand the effect of stoichiometry, GDL characteristics, operating range, and initial conditions (dry vs. flooded) for flow conditions typical of an operating fuel cell. This hysteresis is noted when the air and water flow rates are increased and then decreased along the same path, exhibiting different pressure drops. When starting from dry conditions, the descending pressure drop tended to be higher than the ascending pressure drop at lower simulated current densities. The hysteresis effect was noted for stoichiometries of 1-4 and was eliminated at a stoichiometry of 5. It was found that the hysteresis was greater when water breakthrough occurred at higher simulated current densities, which is a function of GDL properties. The operating range had to reach a critical simulated current density (800 mA cm⁻² in this case) between the ascending and descending approach to create a pressure drop hysteresis zone. The descending step size does not change the size of the hysteresis effect, but a larger step size leads to lower fluctuations in the pressure drop signal. An initially flooded condition also showed hysteresis, but the ascending approach tended to have a higher pressure drop than the descending approach.     

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.     

Experimental investigations on liquid water removal from the gas diffusion layer by reactant flow in a PEM fuel cell

The cross flow from channel to channel through gas diffusion layer (GDL) under the land could play an important role for water removal in proton exchange membrane (PEM) fuel cells. In this study, characteristics of liquid water removal from GDL have been investigated experimentally, through measuring unsteady pressure drop in a cell which has the GDL initially wet with liquid water. The thickness of GDL is carefully controlled by inserting various thicknesses of metal shims between the plates. It has been found that severe compression of GDL could result in excessive pressure drop from channel inlet to channel outlet. Removing liquid water from GDL by cross flow is difficult for GDL with high compression levels and for low inlet air flow rates. However, effective water removal can still be achieved at high compression levels of GDL if the inlet air flow rate is high. Based on different compressed GDL thicknesses, different GDL porosities and permeabilities were calculated and their effects on the characteristics of liquid water removal from GDL were evaluated. Visualization of liquid water transport has been conducted by using transparent flow channel, and liquid water removal from GDL under the land was observed for all the tested inlet air flow rates, which confirms that cross flow is practically effective to remove the liquid water accumulated in GDL under the land area.     

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.     

Numerical simulation of two‐phase cross flow in the gas diffusion layer microstructure of proton exchange membrane fuel cells

The cross flow in the under‐land gas diffusion layer (GDL) between 2 adjacent channels plays an important role on water transport in proton exchange membrane fuel cell. A 3‐dimensional (3D) two‐phase model that is based on volume of fluid is developed to study the liquid water‐air cross flow within the GDL between 2 adjacent channels. By considering the detailed GDL microstructures, various types of air‐water cross flows are investigated by 3D numerical simulation. Liquid water at 4 locations is studied, including droplets at the GDL surface and liquid at the GDL‐catalyst layer interface. It is found that the water droplet at the higher‐pressure channel corner is easier to be removed by cross flow compared with droplets at other locations. Large pressure difference Δp facilitates the faster water removal from the higher‐pressure channel. The contact angle of the GDL fiber is the key parameter that determines the cross flow of the droplet in the higher‐pressure channel. It is observed that the droplet in the higher‐pressure channel is difficult to flow through the hydrophobic GDL. Numerical simulations are also performed to investigate the water emerging process from different pores of the GDL bottom. It is found that the amount of liquid water removed by cross flow mainly depends on the pore's location, and the water under the land is removed entirely into the lower‐pressure channel by cross flow.     

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.     

Numerical investigation of drop dynamics in presence of wettability gradient inside a serpentine channel of proton exchange membrane fuel cell

Water management significantly affects the performance of a proton exchange membrane fuel cell (PEMFC). Therefore, interest is felt to numerically investigate water droplet movement and slug formation inside the microchannel (gas) of PEMFC. Two important parameters—water coverage ratio and pressure drop have been studied in detail. A U-shaped geometry with a round corner is used for this purpose. 3D unsteady-state models are used to study the drop dynamics using commercial CFD software ANSYS FLUENT 18. For tracking of water drop dynamics, the volume of fluid model is used. Two different situations are simulated. In the first case, the investigation of hydrodynamics of the 0.4 mm drop adhered to the surface of the gas diffusion layer (GDL) has been done. In the second case, simulation of air-water slug flow has been done. GDL surfaces at upstream and downstream of bend are modified using user-defined functions, such that the GDL surface has a dynamic contact angle with respect to the direction of flow. This makes it a continuously hydrophilic surface at upstream and continuously hydrophobic surface at downstream with respect to the direction of flow. The impact of GDL wettability on water retention and removal has been discussed. It is noted that the presence of a gradient facilitates the removal of water drop adhered to the GDL surface. For the case of a suspended drop with an increase of 1°/mm in the magnitude of the gradient, a decrease of 30% is observed in water coverage ratio and pressure drop observed in the channel. Such modified surfaces aids in the conversion of slugs to film at the downstream of bend that reduces maldistribution. The pressure fluctuations and average pressure drop are reduced by 66% when subjected to the aforementioned hybrid gradient.     

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.     

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.     

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.     

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.     

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 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.     

Modeling gas flow in PEMFC channels: Part I – Flow pattern transitions and pressure drop in a simulated ex situ channel with uniform water injection through the GDL

The two-phase flow in the gas channels of a proton exchange membrane fuel cell (PEMFC) is studied with an ex situ setup using a gas diffusion layer (GDL) as the sidewall of the channels. Air is supplied at the channel inlet manifold and water is supplied continuously and uniformly through the GDL along the length of the channel. This is different from the simultaneous air and water introduction at the inlet of the channel as studied by previous two-phase flow researchers. The GDL is compressed between the gas channels and the water chambers to simulate PEMFC conditions. The superficial velocity for air and water ranged from 0.25 to 34.5 m/s and 1.54 × 10⁻⁵ to 1.54 × 10⁻⁴ m/s, respectively. The ex situ setup was run in both vertical and horizontal orientations with two GDLs, – Baseline (Mitsubishi Rayon Co. MRC 105 with 5 wt.% PTFE and coated with an in-house MPL by General Motors) and SGL 25 BC – and three channel treatments – hydrophobic, hydrophilic, and untreated Lexan, with contact angles of 116°, 11° and 86°, respectively. No appreciable effect was noted because of the orientation, GDL type or channel coatings. The flow regime is observed at different locations along the channel and is expressed as a function of the superficial air and water velocities. Flow regime criteria are developed and validated against the range of ex situ data observations. A new variable water flow rate pressure drop model is developed in order to account for the variation of water entering the channel at multiple locations along the flow length. Pressure drop models are developed for specific flow regimes and validated against experimental data. The models are able to predict the experimental pressure drop data with a mean error of less than 14%.     

Dynamic water transport and droplet emergence in PEMFC gas diffusion layers

The dynamics of liquid water transport through the gas diffusion layer (GDL) and into a gas flow channel are investigated with an ex situ experimental setup. Liquid water is injected through the bottom surface of the GDL, and the through-plane liquid pressure drop, droplet emergence and droplet detachment are studied. The dynamic behaviour of water transport in and on the surface of the GDL is observed through fluorescence microscopy, and the through-plane liquid pressure drop is measured with a pressure transducer. With an initially dry GDL, the initial breakthrough of liquid water in the GDL is preceded by a substantial growth of liquid water pressure. Post-breakthrough, droplets emerge with a high frequency, until a quasi-equilibrium liquid water pressure is achieved. The droplet emergence/detachment regime is followed by a transition into a slug formation regime. During the slug formation regime, droplets tend to pin near the breakthrough location, and the overall channel water content increases due to pinning and the formation of water slugs. Droplets emerge from the GDL at preferential breakthrough locations; however, these breakthrough locations change intermittently, suggesting a dynamic interconnection of water pathways within the GDL. The experiments are complemented by computational fluid dynamics (CFD) simulations using the volume of fluid method to illustrate the dynamic eruption mechanism.     

Effects of internal flow modification on the cell performance enhancement of a PEM fuel cell

This study presents a numerical investigation on the cell performance enhancement of a proton exchange membrane fuel cell (PEMFC) using the finite element method (FEM). The cell performance enhancement in this study has been accomplished by the transverse installation of a baffle plate and a rectangular block for the modification of flow pattern in the flow channel of the fuel cell. The baffle plates (various gap ratios, λ = 0.005–10) and the rectangular block (constant gap ratio, λ = 0.2) are installed along the same gas diffusion layer (GDL) in the channel at constant Reynolds number for the purpose of investigating the cell performance. The results show that the transverse installation of a baffle plate and a rectangular block in the fuel flow channel can effectively enhance the local cell performance of a PEMFC. Besides, the effect of a rectangular block on the overall cell performance is more obvious than a baffle plate.     

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.     

Effect of GDL compression on pressure drop and pressure distribution in PEMFC flow field

Improving reactant distribution is an important technological challenge in the design of a PEMFC. Flow field and the Gas Diffusion Layer (GDL) distribute the reactant over the catalyst area in a cell. Hence it is necessary to consider flow field and GDL together to improve their combined effectiveness. This paper describes a simple and unique off-cell experimental setup developed to determine pressure as a function of position in the active area, due to reactant flow in a fuel cell flow field. By virtue of the experimental setup being off-cell, reactant consumption, heat production, and water generation, are not accounted as experienced in a real fuel cell. A parallel channel flow field and a single serpentine flow field have been tested as flow distributors in the experimental setup developed. In addition, the interaction of gas diffusion layer with the flow distributor has also been studied. The gas diffusion layer was compressed to two different thicknesses and the impact of GDL compression on overall pressure drop and pressure distribution over the active area was obtained using the developed experimental setup. The results indicate that interaction of GDL with the flow field and the effect of GDL compression on overall pressure drop and pressure distribution is more significant for a serpentine flow field relative to a parallel channel flow field.