In the current work, a modeling simulation based on the Computational Fluids Dynamic (CFD) method was done for two phase flow (gas-liquid) in two bubble columns. The effect of geometrical (two column diameters) and different gas velocity on the system hydrodynamics (gas hold up, mixing time) was investigated. Furthermore, the liquid circulation velocity in the columns was predicted. The simulated results were agreement with the experimental data reported in the literature. It was concluded that gas hold up was decreased while bubble column diameter increased from 15cm to 30cm in both experimental and CFD data's. But gas hold up was increased by increasing superficial gas velocity. In addition, it was seen that increasing of bubble column diameter from 15 to 30 cm leads to decreasing of mixing time. also mixing time reduced when superficial gas velocity was increased from 1 to 10cm/s. Moreover, for more study on CFD simulation of bubble columns with different diameters, the effect of different superficial gas velocity on the liquid circulation velocity was predicted for two reactors. This prediction has shown increasing of liquid circulation velocity versus of increasing superficial gas velocity but it reduced when reactor diameter was increased.