In this work, the influence of the liquid water inlet boundary conditions at the gas diffusion layer (GDL)/catalyst layer interface on the spatial distribution of liquid water within the GDL was studied. We used pore network modeling with invasion percolation to simulate liquid water transport in a commercially available GDL, where the detailed, 3D microstructure of the GDL was obtained through X-ray imaging. Three inlet boundary conditions were studied: uniform pressure (single reservoir), uniform flux (completely discretized reservoirs), and distributed uniform pressure (random spatial-and size-distributions of reservoirs). We presented the distributed uniform pressure boundary condition as a more realistic inlet, where inlets are randomly distributed reservoirs that are connected to multiple inlet pores. It was found that the overall saturation ranged from 6% to 28% when the number of inlets ranged from 20 to 300; however, the GDL/catalyst layer delamination dominated water transport behavior. Proton exchange membrane fuel cells (PEMFCs) are promising electrochemical energy conversion devices; however, at high current densities they are prone to the accumulation of excess liquid water (flooding) in the cathode gas diffusion layer (GDL). This excess water leads to problems such as reduced oxygen transport pathways to the catalyst layer (CL), performance degradation, and reduced efficiency.Pore network modeling is a simulation tool for studying transport phenomena in porous media which is utilized here to investigate the steady state liquid water configurations in GDLs. A pore network model is a network diagram of pores (nodes) and throats (bonds) that represent a porous media, where pores are locations of large void spaces and throats are the local constrictions that connect adjacent pores. Several researchers have developed and applied pore network models to simulate the distribution of liquid water in the GDL.
1-11The introduction of water into the GDL occurs through a variety of mechanisms: electro-osmotic water transport across the membrane, water production from the oxygen reduction reaction, water vapor condensation within the bulk of the GDL, and condensation near the low temperature ribs of the flow field. The spatial distribution of liquid water and saturation (fraction of pore volume occupied by liquid water) within the bulk GDL depends on how liquid water is introduced to this porous layer, specifically at the interface between the GDL and the CL. However, the precise configuration of liquid water inlets at this interface is not well understood.Many authors have provided valuable insight into how inlet conditions heavily influence the spatial distribution of liquid water in the material. 2,4,5,[12][13][14][15] Two types of boundary conditions have been generally considered in isolation for studying liquid water transport in the GDL: uniform pressure boundary condition 2,15 and uniform flux boundary condition. 4 For the uniform pressure boundary condition, all inlet pores along the GDL/CL interface a...