A two-dimensional mathematical model for a planar SOFC (solid oxide fuel cell) is constructed. The distribution of the chemical species, the temperature and the performance (power) and the current density were calculated using a single-unit model with double channels of co-flow pattern. The finite volume method was employed for the calculation. The method was based on the fundamental conservation laws of continuity, momentum, energy and mass. The equations are implemented in FORTRAN language. The effects of several heat sources and flow rates on the calculated results were also investigated. The reference SOFC polarization curve has been calculated by imposing a uniform temperature of 800 K, a pressure equal to 1 atm; H(2) and O(2) molar fractions equal to 0.97 and 0.21 respectively. Results of temperature, chemical species distributions, performance and efficiency under several heat sources are shown and discussed. At a current density of about 23500 A/m(2), the power densities under all sources and chemical sources reached their maximums of 12965 W/m(2) and 16209 W/m(2) (i.e. 25% lower) respectively. However the temperature increment in the anode is analyzed toward all sources and chemical reaction. The temperature maximum values for each heat source type reached 1005.81 K and 984.69 K respectively. Copyright (C) 2010, Hydrogen Energy Publications, LLC. Published by Elsevier Ltd. All rights reserved.