In the present work the effect of the chemical reaction on the temperature field in an anode supported planar SOFC is numerically studied by the aid of a two-dimensional mathematical model. For the model development the mass transport phenomena, the energy conservation, the species flow governed by Darcy's law and the electrochemistry are coupled. The finite difference method is used to solve numerically the system of the equations. The temperature field within each component of the SOFC (interconnection, cathode, anode and electrolyte) is calculated via the mathematical model which is implemented in FORTRAN language. The model results demonstrate the effect of different expressions of the chemical heat source, expressed in terms of enthalpy and entropy, on the temperature field and the location of the higher temperatures that occur within the SOFC during its operation. Copyright (C) 2010, Hydrogen Energy Publications, LLC. Published by Elsevier Ltd. All rights reserved.