A semi-empirical non-isothermal model incorporating coupled momentum, heat and mass transport phenomena for predicting the performance of a proton exchange membrane (PEM) water electrolysis cell operating without flow channels is presented. Model input parameters such as electro-kinetics properties and mean pore size of the porous transport layer (PTL) were determined by rotating disc electrode and capillary flow porometry, respectively. This is the first report of a semi-empirical fully coupled model which allows one to quantify and investigate the effect of the gas phase and bubble coverage on PEM cell performance up to very high current densities of about 5 A/cm(2), The mass transport effects are discussed in terms of the operating conditions, design parameters and the microstructure of the PTL. The results show that, the operating temperature and pressure, and the inlet water flowrate and thickness of the PTL are the critical parameters for mitigating mass transport limitation at high current densities. The model presented here can serve as a tool for further development and scale-up effort in the area of PEM water electrolysis, and provide insight during the design stage. (C) 2017 Hydrogen Energy Publications LLC. Published by Elsevier Ltd. All rights reserved.