In this work, mathematical modeling of SO2 removal reaction with CuO was accomplished by the random pore model. The partial differential equations, describing the reaction of a gaseous reactant with a single pellet and also a packed bed of solid reactant, were solved by the finite element method. The results of modeling consist of CuO conversion-time profiles at different temperatures and SO2 concentrations, and also break through curves which were compared with the literature experimental data in a good accuracy. The rate constants were estimated from the initial slope of the conversion-time curves, and the product layer diffusivities were evaluated from the whole conversion-time profiles. (C) 2013 The Institution of Chemical Engineers. Published by Elsevier B.V. All rights reserved.