The distribution and migration of liquids in various layers of a PEFC is commonly modelled by the 3D flow equations. Given the fact these layers are very thin, there are major problems with such models, including heavy computational efforts and doubts in their applicability to the gas diffusion layer (GDL). Recently, a new approach for modelling multiphase flow through a stack of thin porous layers has been developed [Qin and Hassanizadeh, Int. J. Heat Mass Transfer 70 (2014) 693-708]. In this approach, which is called "reduced continua model", each layer is modelled as a 2D domain with governing equations formulated in terms of thickness-averaged properties. The mass exchange between layers is prescribed by a new constitutive equation. The aim of this paper is to illustrate the implementation of the reduced continua model and show its advantages in modelling liquid water flooding in the GDL and micro porous layer (MPL) of a PEFC. We find that, in comparison to the Richards model, the reduced continua model predicts quite similar water dynamics in the MPL, but a lower steady-state water saturation in the GDL, particularly under the channel area. We provide a quantitative indication of the enormous computational efficiency of the reduced continua model as compared to the Richards model. Finally, the sensitivity studies of major material parameters of the reduced continua model have been provided. Copyright (C) 2015, Hydrogen Energy Publications, LLC. Published by Elsevier Ltd. All rights reserved.