For gas-liquid bubble column with dense vertical bundle internals
the resistance force of bundles on gas and liquid phases was taken into account in the two-phase Euler axisymmetric k-ε model. By introduction of corresponding momentum source
turbulence source and dissipation source
a 2-D CFD model of bubble column with bundle internals was established. The model could clearly and accurately demonstrate the features of gas-liquid flow in bubble column with bundles: "funnel effect" and extension of sparger zone. The calculated two-dimensional distributions of hold-up and liquid velocity agreed with experimental data within a wide range of superficial gas velocity (0.12—0.62 m·s-1).