Fine grid, two dimensional simulations of reactive gas–solid flows occurring in a fluidized bed reactor were carried out using the Eulerian multi-fluid kinetic theory of granular flow (KTGF) approach in the commercial flow solver, ANSYS FLUENT 12.1. The fuel reactor of a pilot scale Chemical Looping Combustion rig, operated in the bubbling fluidization regime at the Vienna University of Technology, was simulated. Grid dependence studies were carried out as well as sensitivity studies to the fuel inlet condition and the inclusion of gas phase turbulence. Simulations could not accurately reproduce the experimental trend for the case when highly reactive nickel oxide was used as the oxygen carrier material, but in general satisfactory quantitative agreement was observed. The failure to correctly capture the experimental trend was primarily attributed to the fine length-scales at the feed gas inlets not being adequately resolved even at the finest grid investigated. The trend quickly worsened when coarser grids were used, indicating that the generality of the model is lost when grid dependence effects are present. A number of possible dimensional effects were also discussed. Subsequently, the model was used to successfully capture another experimental trend obtained with a much less reactive ilmenite oxygen carrier material. The model captured this trend correctly because the reaction was now limited by the reaction rate and not by species transfer to the large scale gas-emulsion interfaces. Results were therefore not as sensitive to the correct hydrodynamic modelling of the interface, especially near the gas inlets, and the model retained its generality over a wide range of operating conditions.