An investigation of electrochemical behavior of PEMFC (proton exchange membrane fuel cell) is performed by using a single-phase two-dimensional finite element analysis. Equations of current balance, mass balance, and momentum balance are implemented to simulate the behavior of PEMFC. The analysis results for the co-flow and counter flow mode of gas flow direction are examined in detail in order to compare how the gas flow direction affects quantitatively. The characteristics of internal properties, such as gas velocity distribution, mass fraction of the reactants, traction of water and current density distribution in PEMFC are illustrated in the electrode and GDL (gas diffusion layer). It is found that the dry reactant gases can be well internally humidified and maintain high performance in the case of the counter-flow mode without external humidification while it is not advantageous for highly humidified or saturated reactant gases. It is also found that the co-flow mode improves the current density distribution with humidified normal condition compared to the counter-flow mode.