A fuel cell is most interesting new power source because it solves not only the environment problem but also natural resource exhaustion problem. CFD modeling and simulation for heat and mass transport in PEM fuel cells are being used extensively in researches and industrial applications to gain better understanding of the fundamental processes and to optimize fuel cell designs before building a prototype for engineering application. In this research, full three-dimensional, non-isothermal computational fluid dynamics model of a tubular-shaped proton exchange membrane (PEM) fuel cell has been developed. This comprehensive model accounts for the major transport phenomena such as convective and diffusive heat and mass transfer, electrode kinetics, transport and phase-change mechanism of water, and potential fields in a tubular-shaped PEM fuel cell. The model explains many interacting, complex electrochemical, and transport phenomena that cannot be studied experimentally. Three-dimensional results of the species profiles, temperature distribution, potential distribution, and local current density distribution are presented and analysed, with the focus on the physical insight and fundamental understanding.