We propose a free boundary model described by coupled Navier-Stokes and chemical reaction equations with discontinuous coefRcients to simulate the chemical re- ¿ctions in viscous drops moving in a viscous incompressible ûuid. Approximation of the solution by a special ñnite element method (FEM) with a method of mapping is discussed. Several numerical resulùs åre presented.