In this paper, the effects of dense vertical internals (rods) on gas holdup and local gas and liquid velocities were investigated by using the Eulerian-Eulerian model coupled with the population balance model. The inclusion of lift and wall lubrication forces was studied by applying different models. The results indicated that just by choosing the appropriate interfacial forces, the numerical model agrees well with the experimental data. A sharper gas holdup, a stronger gas velocity gradient, and a more intense liquid recirculation were observed as the important impacts of the internals. Moreover, three circular internals' arrangements were considered to study the effect of wall and core clearance distances on the bubble column hydrodynamics. The results revealed that by increasing the wall clearance distance, flatter gas holdup and velocity distributions could be achieved. Also, the turbulence parameters were used to evaluate the capability of the model in the prediction of the bubble size distribution. A modification factor for the breakage and coalescence kernels was proposed; the factor enabled the model to reach equilibrium mean diameters for bubbles. The comparison of the probability distribution function (PDF) of bubble sizes in the bubble column with and without internals showed a narrower bubble size distribution with a smaller mean diameter in the presence of internals. (C) 2018 Elsevier Ltd. All rights reserved.