The study of bubble size distributions in direct-contact evaporators was addressed both theoretically and experimentally. Recently developed models for calculating bubble coalescence and breakage frequencies in isothermal bubble columns were adapted to the population balance equation using the bubble mass as the internal coordinate which was discretized using an expansion of the number density function by impulse functions. A sparger model was developed based on experimental data for a non-coalescing system and using bubble formation models for isothermal and non-isothermal conditions. Bubble size distributions in a direct-contact evaporator operating in the quasi-steady-state regime for four different gas superficial velocities, including the homogeneous and heterogeneous regimes, together with the sparger model, were used for estimating the three empirical parameters from the population balance model, which were observed to be functions of the gas superficial velocity. In all cases considered, the population balance model fitted the experimental data rather well and the regressed parameters exhibit the physically expected behavior with changes in the gas superficial velocity. (C) 2004 Elsevier Ltd. All rights reserved.