The effects of internal manifold designs on the reactants feed-stream in Polymer Electrolyte Fuel Cells [PEFCs] is studied to figure out flow and thermal distribution patterns over an entire fuel cell stack. Reactants flows are modeled either laminar of turbulent depending on regions and the open channels in the bipolar plates are simulated by porous media where permeability should be pre-deter-mined for computational analysis. In this work, numerical models for reactants feed-stream In the PEFC manifolds are classified Into two major flow patterns: Z-shape and U-shape. Several types of manifold geometries are analyzed to find the optimal manifold configurations. The effect of heat generation in PEFC on the flow distribution is also Investigated applying a simplified heat transfer model in the stack level (i.e. multi-cell electrochemical power-generation unit). This modeling technique Is well suited for many large scale problems and this scheme can be used not only to account for the manifold flow pattern but also to obtain Information on the optimal design and operation of PEFC systems.