Abstract: Reasonably describing gas-liquid interaction force and bubble size distribution is key to accurately simulate the hydrodynamic parameters of bubble column. For the gas-liquid system, a new bubble breakup model was implemented into the population balance model (PBM) and coupled with the Euler two-fluid model. The mechanisms of bubble coalescence included coalescence due to turbulent eddies, different bubble rise velocities and bubble wake entrainment, and bubble breakup model, energy constraint and pressure constraint were considered. The gas-liquid two-phase flow in a bubble column with an inner diameter of 0.44 m and a height of 2.43 m was simulated. Firstly, the predicted values of bubble breakup rate and daughter bubble size distribution of the breakup model were verified. Then, the gas holdup and axial liquid velocity predicted by coupled model at different superficial gas velocities were compared with experimental values. Moreover, the effect of lift force to model (prediction result) was investigated. The results show that both gas holdup and axial liquid velocity are improved in radial distribution when the lift force is included. The simulation value of the CFD-PBM coupled model agrees well with the experimental data at low and high superficial velocity. Both the turbulent dissipation rate in the reactor and the interaction between the bubbles rise with an increase in superficial gas velocity, which make the distribution range of bubble size widened.