Accurate three-dimensional CFD simulation is of great importance to the flow and heat transfer of flow boiling in the steam generator in nuclear power plant, in which bubbly flow is the main flow pattern. One of the biggest challenges for the CFD simulation is the coalescence and breakup model for bubbles dispersed in the continuous fluid. Recently, several new bubble coalescence and breakup models have been developed combined with the inhomogeneous Multi-Size-Group (MUSIG) model in the Eulerian framework in addition to the standard closure model (Luo and Svendsen (1996), and Prince and Blanch (1990)), which has been widely adopted in the CFD simulations in bubbly flow using the Eulerian method. For example, Liao had developed a comprehensive closure model (2015), in which five collision mechanisms and four breakup mechanisms have been included in the model. However, proper coefficients for the Liao model must be determined to properly address the proportion for various mechanism in the coalescence and breakup models. In this paper, we compared the predictions using the standard closure model and the Liao model (2015) combined with the iMUSIG model for the case of adiabatic upward bubbly flow in vertical pipe, as well as the corresponding experimental results reported by Lucas et al (2008) which were obtained at TOPFLOW test facilities. Six simulated working conditions are chosen. The significance of various mechanisms on the bubble coalescence and breakup were analyzed based on the bubble size distributions in various heights obtained in the CFD predictions and experiments. Finally, a set of proper coefficients for various mechanisms in the coalescence and breakup models of Liao (2015) was developed for the simulated conditions.