In this research, a mathematical model to simulate a two-phase fluidized bed reactor for oxidative dehydrogenation of propane to propylene was developed. The developed model was solved numerically using the MATLAB software. Moreover, effects of the reaction temperature, operating pressure, feed composition, particle's size and gas superficial velocity on the reactant conversion, product selectivity and desired product's yield were understudied. The model was validated using previously published experimental data for simulation of the case at the exit of the circulating fluidized bed reactor. The results revealed that, as the temperature raised from 783 to 823 K, the propane conversion enhanced from 4.47% to 13.01%, while the propylene selectivity lowered from 85.35% to 68.03%. In addition, it was shown that doubling of the propane concentration increased the propane conversion from 4.47% to 4.73%. On one hand, doubling of the oxygen concentration enhanced the propane conversion from 4.41% to 4.47%, while on the other, tripling of the pressure led to increasing of the propane conversion from 4.47% to 5.54%. Nonetheless, this pressure enhancement lowered the propylene selectivity from 85.35% to 77.90%. Moreover, changing of the particle's type from Geldart A to B resulted in decreasing of the propane conversion from 4.56% to 1.40% and increasing of the propylene selectivity from 84.83% to 89.69%. The model further revealed that effect of the superficial velocity was not very significant. Ultimately, to predict the maximum propylene selectivity, different kinetic expressions were understudied.