Numerical simulations of two dimensional gas-solid bubbling fluidized beds with and without immersed horizontal tubes were performed using Eulerian - Eulerian Two Fluid Model (TFM). The influences of immersed horizontal tubes and different data analyses techniques on the bed and bubble hydrodynamics were investigated. The results were compared with experimental data and correlations available in the literature.Different ways for extracting and defining hydrodynamic properties, such as bed expansion ratio and bubble properties, were found to influence the simulation results. Furthermore, the time-averaged values showed greater sensitivity to the length of averaging time in the first few seconds. With regard to tube influence, immersed tubes were found to be the main source of bubble breakup. Thus, the calculated mean bubble diameters and rise velocities were found to be lower with tubes than without for the same bed geometry and superficial velocity. The bubble shapes were observed to elongate in the vertical direction in the tube bank region compared to the bed region below and above the tube bank. In addition, the TFM was found to successfully predict the overall time-averaged solid motion and distribution. For beds with immersed tubes, defluidized regions were observed at the upper part of the tubes where solid particles rested without moving. On the other hand, the lower parts of the tubes were usually covered with gas pockets. These effects were seen to reduce with increasing superficial velocity.