This paper presents a novel approach to model the transient behavior of solid-oxide fuel cell (SOFC) stacks in two and three dimensions. A hierarchical model is developed by decoupling the temperature of the solid phase from the fluid phase. The solution of the temperature field is considered as an elliptic problem, while each channel within the stack is modeled as a marching problem. This paper presents the numerical model and cluster algorithm for coupling between the solid phase and fluid phase. For demonstration purposes, results are presented for a stack operated on pre-reformed hydrocarbon fuel. Transient response to load changes is studied by introducing step changes in cell potential and current. Furthermore, the effect of boundary conditions and stack materials on response time and internal temperature distribution is investigated. (C) 2012 Elsevier B.V. All rights reserved.