Modeling the hydrodynamics of dense-solid gas flows is strongly affected by the wall boundary condition and in particular, the specularity coefficient φ which characterizes the tangential momentum transfer from the particles to the wall. The focus of this study is to investigate the impact of φ on the fluidization hydrodynamics using a fully Eulerian description of the solid and gas phases in 3D cylindrical coordinates. In order to quantify this impact, tools for characterizing the bubbling dynamics and solids circulation are developed and applied to both lab-scale (diameters 10 cm and 14.5 cm) and pilot-scale (diameter 30 cm) cylindrical beds. Comparison of simulation predictions with experimental data for different fluidization regimes and particle properties suggests that values of φ in the range [0.01,0.3] are suitable for simulating most dense solid-gas flows of practical interest. It is also shown that for this range of φ, the fluidization hydrodynamics are not significantly dependent on the choice of φ especially as the bed diameter is increased. Additionally, 3D validation of the variable φ model by Li and Benyahia [1] shows the bubble diameter predictions to be in excellent agreement with experiment and the average value of φ predicted within the range [0.01,0.3]. Quantifying the impact of φ and establishing an appropriate range is not only important for accurate simulations at both lab and pilot scales but also validation of models and sub-models for a better understanding of the fluidization phenomenon. Finally, a comparison of the Gidaspow and Syamlal-O'Brien gas-solids drag model shows that the former is more applicable to homogeneous bubbling fluidization (U/U mf <4) while the latter is only suitable for high velocities (U/U mf >4) associated with larger bubbles and slugs.