Numerical simulations of gas–solid fluidized beds based on the kinetic theory of granular flow exhibit a significant dependence on domain discretization. Bubble formation, bubble size and shape all vary greatly with the discretization, and the use of an inappropriate scale resolution leads to inaccurate predictions of fluidization hydrodynamics. In this study, grid-independent solutions of the two fluid model were examined by comparing the bed expansions obtained from numerical simulations with experimental results and empirical predictions, based on bubbling fluidized beds of Geldart B particles. Grid independence was achieved with a grid resolution equal to 18 times the particle diameter. The simulation results were compared with previously published data for verification purposes. The results of this work should provide a guideline for choosing the appropriate grid size and thereby minimize the time and expense associated with large simulations.