A numerical model is developed to study the dissolution of droplets in a binary mixture with miscibility gap. The moving boundary problem is solved with a modified volume of fraction method to compute the time evolution of the average drop radius and the velocity and concentration distributions around the dissolving drop. The modeling results are presented to explain experimental findings that show stable and oscillatory plumes rising from droplets of methanol that dissolve into a cyclohexane liquid matrix. (C) 2003 Elsevier Ltd. All rights reserved.