In this work, an efficient model for simulating bubble dispersion and coalescence due to turbulence is developed in the Euler-Lagrange framework. The primary liquid phase is solved on the Euler grid with the RANS turbulence model. Bubble motion is computed with the force balance equations. One-way coupling between two phases is assumed and the framework is designed for the computation of disperse bubbly flows at low Eötvös number. The turbulent dispersion of the dispersed phase is reconstructed with the continuous random walk (CRW) model. Bubble-bubble collisions and coalescence are accounted for deterministically. To accelerate the time-consuming search for potential collision partners in dense bubbly flows, the sweep and prune algorithm is employed, which can be utilized in arbitrary mesh types and sizes. Validation against experiments of turbulent pipe flows demonstrates that the one-way coupled EL-CRW dispersion model can well reproduce the bubble distribution in a typical dense bubbly pipe flow. Good agreement of the bubble size distribution at the pipe outlet between the simulation and the experiment is obtained.