A major challenge in the application of proton exchange membrane fuel cells (PEMFCs) is water management, with the flooding of electrodes as the main issue. The Lattice-Boltzmann method (LBM) is a relatively new technique that is superior in modeling the dynamic interface of multiphase fluid flow in complex microstructures such as non-homogeneous and anisotropic porous media of PEMFC electrodes. In this study, the dynamic behavior of a water droplet during removal from gas diffusion layer (GDL) of a PEMFC electrode with interdigitated flow field is simulated using LBM. The effects of GDL wettability and its spanwise and transverse gradients on the removal process are investigated. The results demonstrate great influence of wettability and its spanwise and transverse gradients on the dynamic behavior of droplets during the removal process. Although increasing the hydrophobicity of GDL results in better droplet removal,
its increase beyond a critical value does not show a significant effect.