A comprehensive 3D (three-dimensional) multiphase model of PEMFC (proton exchange membrane fuel cell) is developed, in which the gas and liquid two-phase flow in channel and porous electrodes are investigated in detail. In the simulation of gas and liquid two-phase flow in channels, the effect of surface tension, wall adhesion and gravity is taken into account, including the influence of pressure difference between the inlet and outlet on inlet reactant gas concentration; while in porous electrodes, the anisotropy of GDL (gas diffusion layer) and liquid saturation jump at the interface of two different porous layers (e.g. GDL and MPL (micro-porous layer)) are also considered in this model. It is found that the amount of liquid water in channels increases with the increment of current density. In addition, increasing the contact angle at GDL/channel interface is found to be able to improve the performance of PEMFC by facilitating the water removal process in channels. Moreover, it can be concluded that adding baffles in cathode channel not only increases the oxygen concentration in porous electrodes but also facilitates the water removal process, both of which prevent PEMFC from concentration loss effectively. (C) 2017 Elsevier Ltd. All rights reserved.