Computational fluid dynamics (CFD) modeling of the catalytic ozone decomposition reaction in a circulating fluidized-bed (CFB) riser, using iron-impregnated fluid catalytic cracking particles as catalyst, is carried out. The catalytic reaction is defined as a one-step reaction, and the reaction equation is modified by with respect to the particle surface area, A(p), and an empirical coefficient. The Eularian-Eularian method with the kinetic theory of granular flow is used to solve the gas-solids two-phase flow in the CFB riser. The simulation results are compared with experimental data, and the reaction rate is modified by using an empirical coefficient, to provide better simulation results than the original reaction rate. Moreover, the particle size has great effects on the reaction rate. The generality of the CFD model is further validated under different operating conditions of the riser.