A two-phase flow model was developed for liquid-feed methanol fuel cells (DMFC) to evaluate the effects of various operating parameters on the DMFC performance. In this study, a general homogenous two-dimensional model is described in details for both porous layers and fluid channels. This two-dimensional general model accounts for fluid flow, electrochemical kinetics, current density distribution, hydrodynamics, multi-component transport, and methanol crossover. It starts from basic transport equations including mass conservation, momentum transport, energy balance, and species concentration conservation in different elements of the fuel cell sandwich, as well as the equations for the phase potential in the membrane and the catalyst layers. These governing equations are coupled with chemical reaction kinetics by introducing various source terms. It is found that all these equations are in a very similar form except the source terms. Based on this observation, all the governing equations can be solved using the same numerical formulation in the single domain without prescribing the boundary conditions at the various interfaces between the different elements of the fuel cell. Detailed numerical formulations are presented in this paper. The numerical simulation results, such as velocity field, local current density distribution, and species concentration variation along the flow channel, under various operation conditions are computed. The performance of the DMFC affected by various parameters such as temperature, pressure, and methanol concentration is investigated in this paper. The numerical results are further validated with available experimental data from the published literatures.