Designing the jet ejector optimally is a challenging task and has a great impact on industrial applications. Three different sets of nozzles (namely, 1, 3, and 5) inside the jet ejector are compared in this study by using numerical simulations. More precisely, dynamics of bubble coalescence and breakup in the multinozzle jet ejectors are studied by means of Computational Fluid Dynamics (CFD). The population balance approach is used for the gas phase such that different bubble size groups are included in CFD and the number densities of each of them are predicted in CFD simulations. Here, commercial CFD softwareANSYS Fluent 14.0is used. The realizablek-εturbulence model is used in CFD code in three-dimensional computational domains. It is clear that Reynolds-Averaged Navier-Stokes (RANS) models have their limitations, but on the other hand, turbulence modeling is not the key issue in this study and we can assume that the RANS models can predict turbulence of the carrying phase accurately enough. In order to validate our numerical predictions, results of one, three, and five nozzles are compared to laboratory experiments data for Cl2-NaOH system. Predicted gas volume fractions, bubble size distributions, and resulting number densities of the different bubble size groups as well as the interfacial area concentrations are in good agreement with experimental results.