In this work, the experimental investigation of SO2 removal reaction by CuO was accomplished in a packed bed reactor between 400 and 600 degrees C at 12505000 ppm of inlet SO2. Moreover, the mathematical modeling of this reaction was performed by the random pore model. The partial differential equations in a packed bed reactor were solved using finite element technique. The kinetic parameters were previously determined from thermogravimetry tests on a single-pellet CuO reaction with SO2 at various concentrations and temperatures. The results of modeling consist of packed bed breakthrough curves at different temperatures and SO2 concentrations, which were compared with the experimental results successfully. By increasing the reactor temperature from 400 to 600 degrees C, the appearance of the breakthrough curve shifts from 200 to 550 min. Similarly, by decreasing the inlet SO2 concentration from 5000 to 1250 ppm, the breakthrough profile appears at 900 min instead of 200 min.