In the present study a two-dimensional model of a tubular solid oxide fuel cell operating in a stack is presented. The model analyzes electrochemistry, momentum, heat and mass transfers inside the cell. Internal steam reforming of the reformed natural gas is considered for hydrogen production and Gibbs energy minimization method is used to calculate the fuel equilibrium species concentrations. The conservation equations for energy, mass, momentum and voltage are solved simultaneously using appropriate numerical techniques. The heat radiation between the preheater and cathode surface is incorporated into the model and local heat transfer coefficients are determined throughout the anode and cathode channels. The developed model has been compared with the experimental and numerical data available in literature. The model is used to study the effect of various operating parameters such as excess air, operating pressure and air inlet temperature and the results are discussed in detail. The results show that a more uniform temperature distribution can be achieved along the cell at higher air-flow rates and operating pressures and the cell output voltage is enhanced. It is expected that the proposed model can be used as a design tool for SOFC stack in practical applications. Copyright © 2009 John Wiley & Sons, Ltd.