In this paper we study the mass transfer process of ozone dissolution in a bubble plume inside a rectangular water tank, as a model problem for a water purification system. The effect of bubble diameter and plume structure on mass transfer efficiency of ozone in bubble plumes is investigated numerically. In order to capture the detailed plume structure, the interaction between liquid and bubbles is treated by a two-way coupling Euler-Lagrange method. The motion of the continuous phase (a mixture of liquid and gas bubbles) is solved using a finite difference method in an Eulerian framework. The motion of the dispersed phase (bubbles) is tracked individually in a Lagrangian approach. The ozone transfer process from bubbles to liquid is computed by modelling the mass transfer rate of individual bubbles. Our numerical results show a nonlinear dependence of the ozone dissolution efficiency on the initial bubble size. The dissolution efficiency varies rapidly when the initial bubble size reaches certain value while the change of efficiency is much slower at other bubble sizes. Therefore, for a given tank size it is not necessary to generate bubbles much smaller than the optimal size. This result is of importance for engineering since it is difficult to generate small bubbles in practice. Our results also show that the instantaneous dissolution rate of ozone could be increased by increasing the initial volumetric fraction of ozone inside bubbles even up to 20% while maintaining the dissolution efficiency. (c) 2006 Elsevier Ltd. All rights reserved.