We consider the numerical solution of the free boundary Beoulli problem by employing level set formulations.Using a perturbation technique,we derive a second order method that leads to a fast iteration solver.The iteration procedure is adapted in order to work in the case of topology changes.Various numerical experiments confirm the efficiency of the derived numerical method.