The PEFC has emerged as the most viable fuel cell type for automotive and some portable applications, and also has potential back-up power unit applications due to its low operating temperature, comparative simplicity of construction, high power density, and ease of operation. In spite of tremendous scientific advances, as well as engineering progress over the last few decades, the commercialization of PEFCs remains unrealized, owing primarily to economic viability associated with the high prices of materials and components and technical problems relating primarily to water management. The difficulty in addressing the water management issues lies mostly in the two-phase multi-component flow involving phase-change in porous media, coupled heat and mass transfer, interactions between the porous layers and gas channel (GC) and the complex relationship between water content and cell performance. Due to the low temperature of operation, water generated by the electrochemical reactions often condenses into liquid form, potentially flooding the gas diffusion layer (GDL), GC or other components. Insight into the fundamental processes of liquid water evolution and transport is still lacking, preventing further enhanced PEFC development. The aim of this paper is to give a comprehensive introduction to PEFC modeling inside GCs and GDLs, with a focus on two-phase flow and related phase-change and transport processes. Relevant momentum, mass and heat transport processes are introduced and the microstructural effects on the transport processes inside the porous GDL are extensively discussed. The selection of a computational approach, for the two-phase flow within a GDL or GC, for example, should be based on the computational resources available, concerns about time and scale (microscale, cell scale, stack scale or system scale), as well as accuracy requirements. One important feature, included in some computational approaches, is the possibility to track the front between the liquid and the gas phases. To build a PEFC model, one must make a large number of assumptions. Some assumptions have a negligible effect on the results and reliability of the model. However, other assumptions may significantly affect the result. It is strongly recommended in any modeling paper to clearly state the assumptions being implemented, for others to be able to judge the work. It is important to note that a large fraction of the expressions that presently are used to describe the transport processes inside PEFC GDLs were originally developed to describe significantly different systems, such as sand or rocks. Moreover, the flow pattern maps and pressure drop correlations of two phase flow in micro channels may not be applicable for GCs due to one side wall being porous, with the resulting interaction between the GDL and GC. (C) 2016 Elsevier Ltd. All rights reserved.