A recently developed model for coupled heat and mass transfer in binary systems during the formation and ascension of superheated bubbles was extended to a multicomponent system comprising N volatile species. The model allows variable properties and bubble radius changes, assuming diffusive mass fluxes to be properly described by Fick's law. Experimental direct-contact evaporation tests were conducted with ethyl acetate aqueous solutions to provide data for assessing the developed model. In addition, the model was tested against available literature data for an air-stripper. In both cases, a good agreement between simulation and experimental results was verified. (C) 2004 Elsevier Ltd. All rights reserved.