In the present study, we develop a two-phase, non-isothermal, three-dimensional model of PEMFC to simulate cell performance. The model is based on five conservation equations(mass, momentum, species, energy, charge) and the liquid water equation. In this model, it can describe the phase change of water by evaporation/condensation effect, and the water transportation from electro-osmotic
drag effect.
In the practical application of the PEM fuel cell, liquid water would exist at high current density, it reduce the diffusivity of the gas phase species. If the pores in
a porous media are occupied by liquid water, the gas phase species can not reach the catalyst layer. In past studies, effective medium theory is commonly adopted to calculated effective diffusivity. Here, we modify the effective diffusivity to simulate the area blockage effect due to the liquid water effect, validated by the network model
proposed by Nam and Kaviany[7], the effective medium theory always overestimates the cell performance at high current density. However in our model the novel formulation can predict the cell performance precise validated by experiment data. Finally, the "one-dimensional two phase and non-isothermal" model of MEA with MPL has been adopted. The micro porous layer will reduce the liquid water at gas
diffusion media and also improve the cell performance.

[1] T. E. Springer, T. Rockward, T. A. Zawodzinski, S. Gottesfeld, "Model for Polymer Electrolyte Fuel Cell Operation on Reformate Feed," J. Electrochem. Soc. 148 (1), A11 (2001)
[2] C. Marr, X. Li, "Composition and performance modeling of catalyst layer in a proton exchange membrane fuel cell," J. Power Source 77, 17 (1999)
[3] ZN. Farhat, "Modeling of catalyst layer microstructure refinement and catalyst utilization in a PEM fuel cell," J. Power Source 138, 68 (2004)
[4] L. You, H. T. Liu, "A parametric study of the cathode catalyst layer of PEM fuel cells using a pseudo-homogeneous model," Int. J. Hydrogen Energy 26, 991 (2001)
[5] H. S. Chu, C. Yeh, F. Chen, "Effects of porosity change of gas diffuser on performance of proton exchange membrane fuel cell," J. Power Source 123, 1 (2003)
[6] V. Gurau, F. Barbir, H. T. Liu, "An Analytical Solution of a Half-Cell Model for PEM Fuel Cells," J. Electrochem. Soc. 147 (7), 2468 (2000)
[7] J. H. Nam, M. Kaviany, "Effective diffusivity and water-saturation distribution in single- and two-layer PEMFC diffusion medium," Int. J. heat and mass transfer 46, 4595 (2003)
[8] Z. H. Wang, C. Y. Wang, K. S. Chen, "Two-phase °ow and transport in the air cathode of proton exchange membrane fuel cells," J. Power Source 94, 40 (2001)
[9] G. Lin, W. He, T. V. Nguyen, "Modeling Liquid Water Effects in the Gas Diffusion and Catalyst layers of the Cathode of a PEM Fuel Cell," J.
Electrochem. Soc. 151 (12), A1999, (2004)
[10] T. F. Fuller, J. Newman, "Water and Heat Management Model for Solid-Polymer-Electrolyte Fuel cells," J. Electrochem. Soc. 140 (5), 1218 (1993)
[11] H. Ju, H. Meng, C. Y. Wang, "A single-phase, non-isothermal model for PEM fuel cells," I. J. Heat and Mass Transfer 48, 1303 (2005)
[12] Hua Meng, Chao-Yang Wang, "Large-scale simulation of polymer electrolyte fuel cells by parallel computing," Chemical Engineering Science 59, 3331 (2004)
[13] Sukkee Um, C. Y. Wang, "Three-dimensional analysis of transport and electrochemical reactions in polymer electrolyte fuel cells," J. Power Source 125, 40 (2004)
[14] Hua Meng and Chao-Yang Wang, "Model of Two-Phase Flow and Flooding Dynamics in Polymer Electrolyte Fuel Cells," J. Electrochem. Soc. 152 (9), A1733 (2005)
[15] Qiang Ye, T V Nguyen, "Three-Dimensional Simulation of Liquid Water Distribution in a PEMFC with Experimentally Measured Capillary Functions"J. Electrochem. Soc. 154 (12), B1242 (2007)
[16] S.Mazumder, and J. V. Cole, "Rigorous 3-D Mathematical Modeling of PEM fuel cells I.Model Predictions without LiquidWater Transport" J. Electrochem. Soc. 150 (11), A1503 (2003)
[17] S.Mazumder, and J. V. Cole, "Rigorous 3-D Mathematical Modeling of PEM fuel cells II.Model Predictions with Liquid Water Transport" J. Electrochem. Soc. 150 (11), A1510 (2003)
[18] P. C. Sui, S. Kumar, N. Djilali, "Advanced computational tools for PEM fuel cell design Part 1.Development and base case simulations" J. Power Source
180, 410 (2008)
[19] P. C. Sui, S. Kumar, N. Djilali, "Advanced computational tools for PEM fuel cell design Part2.Detialed experimental validation and parametric study" J.Power Source 180, 423 (2008)
[20] E. A. Ticianelli, C. R. Derouin, and S. Srinivasan "Localization of platinum in low catalyst loading electrodes to to attain high power densities in SPE fuel cells" J. Electroanal. Chem. Interfacial Electrochem, 251,275(1988)
[21] E. A. Ticianelli, C. R. Derouin, A. Redondo, and S. Srinivasan "Methods to advance technology of proton-exchange membrane fuel-cells" J. Electrochem. Soc, 135,2209(1988)
[22] U. Pasaogullari, C.Y.Wang,"Two-phase transport and the role of micro-porous layer in polymer electrolyte fuel cells"Electrochumica Acta. 49,4359(2004)
[23] Z. Zhan, J. Xiao, D. Li, M. Pan, R. Yuan,"E®ects of porosity distribution variation on the liquid water °ux through gas diffusion layers of PEM fuel cells"J. power sources 160,1041(2006)
[24] J. Newman, Electrochemical Systems, Prentice-Hall, Englewood Cliffs, NJ1911,pp. 290-291.
[25] R. B. Bird, W. Stewart, and E. N. Lightfoot, Transport Phenomena, Wiley, New york(1960)
[26] M.M. Mezedur, M.Kaviany, W.Moore, "Effect of pore structure, randomness and size on effective mass diffusivity" AIChE Journal 48(1), 15 (2002)
[27] M. M. Tomadakis, S. V. Sotirchos, "Ordinary and transition regime in random fiber structure" AIChE Journal 39(3),397 (1993)
[28] J.J. Lin and C.A. Lin, "A non-isothermal two-phase flow MEA model for proton exchange membrane fuel cell," National Tsing Hua University, Taiwan (2008)
[29] C. Y. Wang and P. Cheng, Adv. Heat Transfer,30,93(1997)
[30] G. Dagan, Flow and Transport in Porous Formations, Springer-Verlag, Berlin(1989)
[31] K. S. Udell, "Heat transfer in porous media considering phase change and capillary-the heat pipe effect" Int. J. Heat mass transfer 28,485-495(1985)
[32] CFD-ACE V2007 User Manual
[33] CFD-ACE V2007 Modules Manual V1