3D Eulerian-Eulerian model was applied to simulate the gas-liquid two-phase flow in a stirred tank of dual Rushton turbines using computational fluid dynamics (CFD). The effects of two different bubble treatment methods (constant bubble size vs. population balance model, PBM) and two different coalescence models (Luo model vs. Zaichik model) on the prediction of liquid flow field, local gas holdup or bubble size distribution were studied. The results indicate that there is less difference between the predictions of liquid flow field and gas holdup using the above models, and the use of PBM did not show any advantage over the constant bubble size model under lower gas holdup. However, bubble treatment methods have great influence on the local gas holdup under larger gas flow rate. All the models could reasonably predict the gas holdup distribution in the tank operated at a low aeration rate except the region far from the shaft. Different coalescence models have great influence on the prediction of bubble size distribution (BSD). Both the Luo model and Zaichik model could qualitatively estimate the BSD, showing the turning points near the impellers along the height, but the quantitative agreement with experiments is not achieved. The former over-predicts the BSD and the latter under-predicts, showing that the existing PBM models need to be further developed to incorporate more physics.