Dynamic behaviour of liquid water emerging from a GDL pore into a PEMFC gas flow channel

A numerical investigation of the dynamic behaviour of liquid water entering a polymer electrolyte membrane fuel cell (PEMFC) channel through a GDL pore is reported. Two-dimensional, transient simulations employing the volume of fluid (VOF) method are performed to explicitly track the liquid–gas interface, and to gain understanding into the dynamics of a water droplet subjected to air flow in the bulk of the gas channel. The modeled domain consists of a straight channel with air flowing from one side and water entering the domain from a pore at the bottom wall of the channel. The channel dimensions, flow conditions and surface properties are chosen to be representative of typical conditions in a PEMFC. A series of parametric studies, including the effects of channel size, pore size, and the coalescence of droplets are performed with a particular focus on the effect of geometrical structure. The simulation results and analysis of the time evolution of flow patterns show that the height of the channel as well as the width of the pore have significant impacts on the deformation and detachment of the water droplet. Simulations performed for droplets emerging from two pores with the same size into the channel show that coalescence of two water droplets can accelerate the deformation rate and motion of the droplets in the microchannel. Accounting for the initial connection of a droplet to a pore was found to yield critical air inlet velocities for droplet detachment that are significantly different from previous studies that considered an initially stagnant droplet sitting on the surface. The predicted critical air velocity is found to be sensitive to the geometry of the pore, with higher values obtained when the curvature associated with the GDL fibres is taken into account. The critical velocity is also found to decrease with increasing droplet size and decreasing GDL pore diameter.     

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.     

Effects of hydrophilic/hydrophobic properties of gas flow channels on liquid water transport in a serpentine polymer electrolyte membrane fuel cell

The droplet dynamics in the serpentine flow channel of a hydrogen fuel cell has been numerically investigated to obtain ideas for designing a serpentine channel with the aim of effectively preventing flooding. Three-dimensional two-phase flow simulations employing the volume of fluid (VOF) method have been performed. Liquid droplets emerging from four adjacent pores at the hydrophobic bottom wall are subjected to airflow in the bulk of the serpentine flow channel. The effects of contact angle variation of the channel walls on liquid water removal have been tested in terms of liquid water saturation and coverage of liquid water on the gas diffusion layer (GDL) surface. The numerical results show that the hybrid case, which consists of hydrophilic channel walls at the straight part and hydrophobic walls at the turning part of the serpentine flow channels, enhances water removal compared with two other cases in which the channel wall is homogeneously hydrophilic or hydrophobic. The three-dimensional visualization of liquid water droplets reveals that the hydrophobic wall at the turning part reduces the water saturation in the channel and the hydrophilic wall at the straight part prevents the liquid water from covering the GDL surface.     

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.     

Numerical simulation of liquid water emerging and transport in the flow channel of PEMFC using the volume of fluid method

Three-dimensional numerical simulation of liquid water emerging from the gas diffusion layer (GDL) surface to the gas flow channel in the proton exchange membrane (PEM) fuel cell (PEMFC) is carried out using the volume of fluid (VOF) method. The effects of the water velocity in the GDL hole, the airflow velocity and the wettability of the channel surfaces on the water emerging process and transport in the flow channel are investigated. It is found that at low water velocity, the water detaches from the water hole, forming discrete water droplets on the GDL surface, and is transported downstream on the GDL surface until removed from the GDL surface by the U-turn part of the flow channel; whereas at high water velocity, the continuous water column impinges the hydrophilic channel surface counter to the GDL surface, being directly removed from the GDL surface. The airflow velocity affects water detachment and impact process in the channel corner, and water droplet breakup is observed under high airflow velocity. The channel surface wettability influences water droplet shape and its transport in the channel. Rather than forming corner water films at the U-turn for hydrophilic channel surface, water maintains the droplet shape and smoothly passes through the U-turn for hydrophobic channel surface. The importance of the U-turn to the water removal is also discussed. The U-turn promotes water removal from the GDL surface at low water velocity and water breakup at high airflow velocity.     

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.     

Comparative analysis of two-phase flow in sinusoidal channel of different geometric configurations with application to PEMFC

The droplet dynamics inside a sinusoidal channel for PEMFC (polymer electrolyte membrane fuel cell) are investigated numerically using the VOF (volume of fluid) method. This study is done for three geometrically different channels corresponding to various non-dimensional sinusoidal distances (50, 25, 12.5, 16.7 and 8.3). The effects of key parameters like sinusoidal distance (pitch-amplitude ratio), radius of curvature and wall contact angle on the droplet removal in the flow channel are investigated. The performance of the sinusoidal as compared to the conventional channel is studied based on droplet removal rate and GDL (gas diffusion layer) surface water coverage. It is found that the droplet removal rate increases with increasing sinusoidal distance and wall contact angle. In addition, decrease in the sinusoidal distance results in a significant reduction in the average droplet speed and gas diffusion layer surface water coverage. It was also observed that broken bits of the droplet stuck on the wall corners accrued with a reduction in the wall contact angle. The curvy nature of the side walls generally induces a secondary flow effect which would be most beneficial in enhanced reactant diffusion and cell performance. It is suggested that the sinusoidal distance and wall contact angle effect on two-phase flow in a channel is highly significant. As such, needs to be considered for water management in sinusoidal channels.     

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.     

Lattice Boltzmann simulations of liquid droplet dynamic behavior on a hydrophobic surface of a gas flow channel

Using the multiphase free-energy lattice Boltzmann method (LBM), the formation of a water droplet emerging through a micro-pore on the hydrophobic gas diffusion layer (GDL) surface in a proton exchange membrane fuel cell (PEMFC) and its subsequent movement on the GDL surface under the action of gas shear are simulated. The dynamic behavior of the water droplet emergence, growth, detachment and movement in the gas flow channel is presented. The size of the detached droplet and the time of the droplet removing out of the channel under the influence of gas flow velocity and GDL surface wettability are investigated. The results show that water droplet removal is facilitated by a high gas flow velocity on a more hydrophobic GDL surface. A highly hydrophobic surface is shown to be capable of lifting the water droplet from the GDL surface, resulting in more GDL surface available for gas reactant transport. Furthermore, an analytical model based on force balance is presented to predict the droplet detachment size, and the predicted results are in good agreement with the simulation results. It is shown that the LBM approach is an effective tool to investigate water transport phenomena in the gas flow channel of PEMFCs with surface wettability taken into consideration.     

Numerical simulations of two-phase flow in a proton-exchange membrane fuel cell based on the generalized design method

The water management of proton-exchange membrane fuel cell (PEMFC) has a major impact on the performance of the cell system. In order to investigate the influence of air velocity and wettability on the whole process during penetration of liquid water, a generalized two-dimensional model in conjunction with the volume of fluid (VOF) method was used to simulate the whole processes from gas diffusion layer (GDL) to gas channel (GC). The results show that the wettability of the medium plays a significant role than flow rate for the penetration of liquid water in the GDL. It is shown that favorable hydrophobicity and high air velocity in GC is helpful to remove liquid droplets on the GDL surface. By contrast, the stable droplets spacing on GDL surface is more concentrated and the percentage of liquid area is more extensive under the hydrophilic and low-velocity case, which would aggravate the liquid water and hard to remove from the GDL surface.     

Experimental investigation of water droplet–air flow interaction in a non-reacting PEM fuel cell channel

It has been well documented that water production in PEM fuel cells occurs in discrete locations, resulting in the formation and growth of discrete droplets on the gas diffusion layer (GDL) surface within the gas flow channels (GFCs). This research uses a simulated fuel cell GFC with three transparent walls in conjunction with a high speed fluorescence photometry system to capture videos of dynamically deforming droplets. Such videos clearly show that the droplets undergo oscillatory deformation patterns. Although many authors have previously investigated the air flow induced droplet detachment, none of them have studied these oscillatory modes. The novelty of this work is to process and analyze the recorded videos to gather information on the droplets induced oscillation. Plots are formulated to indicate the dominant horizontal and vertical deformation frequency components over the range of sizes of droplets from formation to detachment. The system is also used to characterize droplet detachment size at a variety of channel air velocities. A simplified model to explain the droplet oscillation mechanism is provided as well.     

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.     

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.     

3D lattice Boltzmann modeling of droplet motion in PEM fuel cell channel with realistic GDL microstructure and fluid properties

A 3D multi-component multi-phase lattice Boltzmann model is developed to study the droplet motion in the flow channel of proton exchange membrane fuel cell. The model is capable of reaching realistic density and viscosity ratio, tunable surface tension with low spurious velocity and is also validated by various benchmark tests in both static and dynamic states. For the first time, the effect of realistic microstructure of gas diffusion layer (GDL) on droplet dynamic behavior is comprehensively studied in terms of comparison with smooth channel, contact angle and droplet size with the motion processes clearly illustrated. The simulation results show the GDL microstructure can amplify the material wettability, affect the motion direction and impede the droplet motion. More hydrophobic GDL can effectively accelerate the transport. However, it is observed the droplet may reach the sidewall due to the presence of GDL and the motion is therefore severely impeded regardless of the GDL contact angle or droplet size, which is hard to avoid but deadly for the water management. For this problem, a novel water management strategy is proposed and the results show the hydrophilic side & top wall can effectively remove the liquid water from the GDL surface, decrease pressure drop and prevent reactant maldistribution.     

Liquid water transport between graphite paper and a solid surface

We studied the interaction of a water droplet with a solid wall on a hydrophobic gas diffusion layer (GDL). Of particular interest is the stability of the droplet as a function of plate wetting properties and the potential for liquid entrapment in the GDL/land contact area. Such transport is of relevance to breakthrough dynamics and convective liquid droplet transport in polymer electrolyte membrane (PEM) fuel cell cathode gas channels. While a variety of complex coupled transport phenomena are present in the PEM fuel cell gas channel, we utilize a very simplified experimental model of the system where a droplet originally placed on a hydrophobic GDL is translated quasistatically across the GDL surface by a solid surface. Transport and entrapment are imaged using fluorescence microscopy. This work provides new insights into droplet behaviour at the GDL/land interface in a PEM fuel cell and suggests that hydrophobic land areas are preferable for mitigating the accumulation of liquid water under the land area of the gas flow channels.     

Effect of hydrophilic pipe structure of proton exchange membrane fuel cell on water removal from the gas diffusion layer surface

In a proton exchange membrane fuel cell (PEMFC), effective GDL surface water elimination is significant to water management. This paper used the volume-of-fluid method (VOF) method to carry out simulation research on transferring liquid water in the flow channel with a hydrophilic pipe. The findings indicated that compared with a straight channel, a hydrophilic pipe structure could effectively remove water from the gas diffusion surface (GDL) and reduce the surface water coverage of the GDL. With the increase in the diameter and height of the pipe structure, the GDL surface's water coverage first increased and then decreased, and it was less with the pipe structure than with the direct flow channel. The removal rate of water on the GDL surface was accelerated. The spacing of hydrophilic pipes has a significant impact on the transportation of water. As the spacing increases, the removal rate of water on the GDL surface slowed. A hydrophilic pipe structure with a diameter of 75 μm, a height of 400 μm, and spacing of 300 μm has good water removal performance on the GDL surface. This research work proposes a new internal structure design of the flow channel, which has specific implications for removing water on the GDL surface.     

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.     

Three-dimensional numerical simulation of water droplet emerging from a gas diffusion layer surface in micro-channels

The dynamic formation of water droplets emerging from a gas diffusion layer (GDL) surface in micro-channels was simulated using the volume of fluid (VOF) method. The influence of GDL surface microstructure was investigated by changing the pore diameter and the number of pore openings on the GDL surface. Simulation results show that the microstructure of the GDL surface has a significant impact on the two-phase flow patterns in gas flow channels. For a non-uniform GDL surface, three stages were identified, namely emergence and merging on the GDL surface, accumulation on the channel sidewalls and detachment from the top wall. It was also found that if the pore size is small enough, the flow pattern in the channel does not change with further reduction in the pore diameter. However, the two-phase flow patterns change significantly with the wettability of the GDL surface and sidewalls, but remain the same when the liquid flow rate is reduced by two orders of magnitude from the reference case.     

Roughness effects of gas diffusion layers on droplet dynamics in PEMFC flow channels

Water management remains one of the major challenges in optimising the performance of PEMFCs, in which liquid accumulation and removal in gas diffusion layers (GDLs) and flow channels should be addressed. Here, effects of GDL surface roughness on the water droplet removal inside a PEMFC flow channel have been investigated using the Volume of Fluid method. Rough surfaces are generated according to realistic GDL properties by incorporating RMS roughness and roughness wavelength as the main characteristic parameters. Droplet dynamics including emergence, growth, detachment, and removal in flow channels with various airflow rates are simulated on rough substrates. The influences of airflow rate on droplet dynamics are also discussed by comparing the detachment time and droplet morphology. The liquid removal efficiency subject to different surface roughness parameters is evaluated by droplet detachment time and elongation, and regimes of detachment modes are identified based on the droplet breakup location and detachment ratio. The results suggest that rough surfaces with higher RMS roughness can facilitate the removal of liquid inside flow channel. Whilst surface roughness wavelength is found less significant to the liquid removal efficiency. The results here provide qualitative assessments on identifying the key surface characteristics controlling droplet motion in PEMFC channels.     

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.