Investigation of two-phase flow in the compressed gas diffusion layer microstructures

This study aims to investigate how multiple parameters affect the two-phase flow in compressed gas diffusion layer (GDL). A stochastic model is adopted to reconstruct the GDL microstructures. Solid mechanics simulations on the reconstructed GDL microstructures are performed, based on the finite element method (FEM). Various pore morphologies and distributions of compressed GDLs are observed. Two-phase flow in GDL is simulated using a volume of fluid (VOF) model. Corner droplet (on the GDL surface) and water flow (emerging from GDL bottom) are considered. It is found that two-phase flow in the GDL is highly influenced by compression, fiber diameter, porosity, and GDL thickness. The results indicate that a larger fiber diameter or higher porosity contributes to the water transport due to larger average pore size. Furthermore, water removal from a thicker GDL is more difficult, whereas water transport in the lower part of a compressed thick GDL is easy.     

Numerical simulation of two-phase flow in gas diffusion layer and gas channel of proton exchange membrane fuel cells

Liquid water within the cathode Gas Diffusion Layer (GDL) and Gas Channel (GC) of Proton Exchange Membrane Fuel Cells (PEMFCs) is strongly coupled to gas transport properties, thereby affecting the electrochemical conversion rates. In this study, the GDL and GC regions are utilized as the simulation domain, which differs from previous studies that only focused on any one of them. A Volume of Fluid (VOF) method is adopted to numerically investigate the two-phase flow (gas and liquid) behavior, e.g., water transport pattern evolution, water coverage ratio as well as local and total water saturation. To obtain GDL geometries, an in-house geometry-based method is developed for GDL reconstruction. Furthermore, to study the effect of GDL carbon fiber diameter, the same procedure is used to reconstruct three GDL structures by varying the carbon fiber diameter but keeping the porosity and geometric dimensions constant. The wall wettability is introduced with static contact angles at carbon fiber surfaces and channel walls. The results show that the GDL fiber microstructure has a significant impact on the two-phase flow patterns in the cathode field. Different stages of two-phase flow pattern evolution in both cathode domains are observed. The liquid water in the GDL experiences water invasion, spreading, and rising, followed by the droplet breakthrough in the GDL/GC interface. In the GC, the water droplets randomly experience accumulation, combination, attachment, and detachment. Due to the difference in surface wettability, the water coverage of the GDL/GC interface is smaller than that of the channel side and top walls. It is also found that the water saturation inside the GDL stabilizes after the water breakthrough, while local water saturation at the interface keeps irregular oscillations. Last but not the least, a water saturation balance requirement between the GDL and GC is observed. In terms of varying fiber diameter, a larger fiber diameter would result in less water saturation in the GDL but more water in the GC, in addition to faster water movement throughout the total domain.     

Three-dimensional numerical simulations of water droplet dynamics in a PEMFC gas channel

The dynamic behavior of liquid water emerging from the gas diffusion layer (GDL) into the gas flow channel of a polymer electrolyte membrane fuel cell (PEMFC) is modeled by considering a 1000 μm long air flow microchannel with a 250 μm × 250 μm square cross section and having a pore on the GDL surface through which water emerges with prescribed flow rates. The transient three-dimensional two-phase flow is solved using Computational fluid dynamics in conjunction with a volume of fluid method. Simulations of the processes of water droplet emergence, growth, deformation and detachment are performed to explicitly track the evolution of the liquid–gas interface, and to characterize the dynamics of a water droplet subjected to air flow in the bulk of the gas channel in terms of departure diameter, flow resistance coefficient, water saturation, and water coverage ratio. Parametric simulations including the effects of air flow velocity, water injection velocity, and dimensions of the pore are performed with a particular focus on the effect of the hydrophobicity of the GDL surface while the static contact angles of the other channel walls are set to 45°. The wettability of the microchannel surface is shown to have a major impact on the dynamics of the water droplet, with a droplet splitting more readily and convecting rapidly on a hydrophobic surface, while for a hydrophilic surface there is a tendency for spreading and film flow formation. The hydrophilic side walls of the microchannel appear to provide some benefit by lifting the attached water from the GDL surface, thus freeing the GDL-flow channel interface for improved mass transfer of the reactant. Higher air inlet velocities are shown to reduce water coverage of the GDL surface. Lower water injection velocities as well as smaller pore sizes result in earlier departure of water droplets and lower water volume fraction in the microchannel.     

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

Influence of the surface microstructure of the fuel cell gas diffusion layer on the removal of liquid water

The effective removal of water on the gas diffusion layer (GDL) surface and low flow channel resistance are essential for the water management of proton exchange membrane fuel cells (PEMFCs). In this paper, a 3D two-phase volume of fluid (VOF) model is used to compare and analyze the influence of different GDL surface microstructures on the liquid hydrodynamic behavior and optimize the design of the sine wave microstructure. The results show that the surface microstructure of the GDL has a more significant impact on the water removal and flow resistance coefficient in the flow channel, and the sine wave microstructure has substantial advantages. The sine wave peak and period significantly influence the water removal and flow resistance coefficient. As the peak increases, the average relative change rate and the flow resistance coefficient also increase; the influence of the period is opposite to the peak, and the continuous decrease of the period will accelerate the water removal in the flow channel. The sine wave's height and width have little impact on water removal and the flow resistance coefficient. When the sine wave A = 75 μm, T = 1.5, H = 15 μm, and L = 25 μm, good flow channel water removal and low resistance are achieved. This work has particular guiding significance for removing liquid water on the GDL surface and obtaining low flow channel resistance.     

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.     

Liquid water distribution in hydrophobic gas-diffusion layers with interconnect rib geometry: An invasion-percolation pore network analysis

Water distribution in gas diffusion layers (GDLs) of polymer electrolyte membrane fuel cells (PEMFCs) is determined by the pore morphology of the GDL as well as the flow conditions between the GDL and the gas flow field, where interconnect ribs and gas channels are placed side-by-side. The present study employs a steady state pore network model based on the invasion-percolation (IP) process to investigate the water transport in the under-rib region, in the under-channel region, and in between those regions inside the GDL. The interconnect rib partially blocks the outlet surface of the GDL, which forces water transport paths from the under-rib region to grow towards the gas channel through an extra IP process. The pore network model predicts spatially non-uniform water distributions inside the GDL due to the interconnect ribs, especially with an increased saturation level in the under-rib region. Parametric studies are also conducted to investigate the effects of several geometrical factors, such as width of the rib and the channel, thickness of the GDL, and water intruding condition at the inlet surface of the GDL.     

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.     

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.     

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

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

The gas diffusion layer in polymer electrolyte membrane fuel cells: A process model of the two-phase flow

A two-phase flow process model for the gas diffusion layer (GDL) of a polymer electrolyte membrane fuel cell, considering also the cathode catalyst layer (CL), is presented. For this purpose, a systematic analysis of the factors affecting flooding and drying, including the liquid accumulation in the gas channel (CH), was performed using a one-dimensional reference model for the GDL and a compact channel model. The treatment proposed for the CH-GDL interface was compared with other boundary conditions in the literature. It was concluded that the liquid accumulation in the channel is determinant for estimating the steady state and transient GDL flooding, but that predicting the saturation level in the CL can help for determining operation policies for precluding flooding in the GDL-CL composite, in the absence of an adequate channel model. Bifurcation behavior, associated with the water phase change, was identified by means of the compact model.     

Development of simulated gas diffusion layer of polymer electrolyte fuel cells and evaluation of its structure

In polymer electrolyte fuel cell (PEFC), it is important to understand the behavior of liquid water in gas diffusion layer (GDL) which is one of the constructional elements so as to improve the output performance and the durability. As this behavior of liquid water is attributed to not only the hydrophilicity but also inhomogeneous structure, it is needed to examine in consideration of an actual GDL structure. In this study, as the basic examination of two-phase flow analysis in an actual GDL, a simulated GDL was made by numerical analysis considering the fiber placement. Furthermore, the prediction methods for pore size distribution, permeability and tortuosity of this simulated GDL were developed with the numerical analysis. These parameters of flow and mass transfer were compared with other studies, and the validity of this simulated GDL was confirmed. In addition, effective diffusion coefficient was calculated from tortuosity in simulated GDL, and PEFC output performance was evaluated by a simple model. Moreover, the optimal GDL was examined in consideration of the effect of porosity and fiber diameter at the fiber level.     

PEMFC碳纸气体扩散层内气液两相流格子Boltzmann模拟

为了提高质子交换膜燃料电池（PEMFC）水管理,本文借助多相流格子Boltzmann模型（LBM）模拟分析了PEMFC碳纸气体扩散层（GDL）内的气液两相输运过程,主要研究了GDL疏水性对气液两相流的影响。结果表明：液态水流路径不仅受到GDL结构形态的影响,而且受到材料疏水性影响。液态水在疏水性弱的GDL中不仅容易沁入,而且容易在孔隙中达到饱和;相反,在疏水性较强的GDL中,液态水很难突破沁入小尺寸孔隙,而从孔径较大的孔隙流通,从而形成毛细力主导的指进流动。     

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.     

Dynamic behavior of droplet transport on realistic gas diffusion layer with inertial effect via a unified lattice Boltzmann method

The dynamic behavior of liquid droplets on a reconstructed real gas diffusion layer (GDL) surface with the inertial effect produced by the three dimensional (3D) flow channel is investigated using an improved pseudopotential multiphase model within the unified lattice Boltzmann model (ULBM) framework, which can realize thermodynamic consistency and tunable surface tension. The microstructure of the GDL (Toray-090) including carbon fibers and polytetrafluoroethylene (PTFE) is reconstructed by a stochastic and mixed-wettability model. The critical force formulation for the Cassie-Wenzel transition of a droplet on GDL surface is derived. The effects of inertia and contact angles on the liquid droplet transport process on a reconstructed real GDL surface with a 3D flow channel are investigated. The results show the normalized center-of-mass coordinate X may enter the channel wall area or fluctuate around the initial position. With increased inertia applied on the droplet, the normalized center-of-mass coordinate Y grows faster and the normalized center-of-mass coordinate Z decreases. It is found by the ULBM for the first time that the liquid droplet is pushed back into the GDL by inertial effect. With the increase of inertia and the decrease of contact angle of GDL, both the droplet penetration depth in GDL and the droplet invasion fraction increase. The droplet invasion fraction in GDL is up to 30%.     

Scale effect and two-phase flow in a thin hydrophobic porous layer. Application to water transport in gas diffusion layers of proton exchange membrane fuel cells

Pore network simulations are performed to study water transport in gas diffusion layers (GDLs) of polymer electrolyte membrane fuel cells (PEMFCs). The transport and equilibrium properties are shown to be scale dependent in a thin system like a GDL. A distinguishing feature of such a thin system is the lack of length scale separation between the system size and the size of the representative elementary volume (REV) over which are supposed to be defined the macroscopic properties within the framework of the continuum approach to porous media. Owing to the lack of length scale separation, two-phase flow traditional continuum models are expected to offer poor predictions of water distribution in a GDL. This is illustrated through comparisons with results from the pore network model. The influence of inlet boundary conditions on invasion patterns is studied and shown to affect greatly the saturation profiles. The effects of GDL differential compression and partial coverage of outlet surface are also investigated.     

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.     

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.