In this study, the population balance model (PBM) is coupled with computational fluid dynamics (CFD) to investigate the steady-state bubble size distribution in two types of process equipment namely, a standard Rushton turbine stirred tank reactor and a generic lab-scale flotation cell. The coupling is realized using Fluent 15.07 software, and the numerical model is validated for the stirred tank reactor. The population balance equation (PBE) is solved using the quadrature method of moments (QMOM) technique along with a correction procedure implemented to check and correct invalid moment sets. The breakage and coalescence of bubbles due to turbulence are considered. The breakage rate and daughter size distribution models proposed by Laakkonen et al. (2007) are considered. For modeling coalescence rate, models proposed by Coulaloglou and Tavlarides (1977) are considered. The interaction between the phases is handled by considering the drag model proposed by Lane et al. (2005) while ignoring the other interphase forces. The correction algorithm has been successfully implemented, and improved predictions of gas volume fraction and Sauter mean diameter (SMD, d(32)) have been observed with a good match between the predictions and experimental measurements. The local SMD predictions are compared against predictions from the past studies and the superiority of the current approach for moderate gassing rates is established. The CFD-PBM approach is then used to study and characterize different flow regimes occurring in a generic mechanical flotation cell at different aeration rates and impeller rotation speeds. Also, power numbers are calculated from torque data and are found to drop considerably with an increase in aeration rate and impeller rotation speed as the flow regime approaches recirculating flow. The predicted SMD for flotation cell indicates that smaller bubbles are concentrated near the high turbulence impeller stream, the lower recirculation region, and close to the tank walls. On the other hand, large bubbles are formed in the upper tank region and are concentrated around the shaft during the flooding, loading, and transition flow regimes. In the future, the corrected QMOM approach will be further extended by implementing kinetic models capable of predicting the flotation rate constant using local bubble size information obtained from CFD-PBM simulations. (C) 2016 Elsevier Ltd. All rights reserved.