This paper studies numerically the two-phase flow in the single-bubble electrospinning process by solving the modified Navier-Stokes equations under the influence of electric field, and the interface between the two fluids has been determined by using the volume of fluids method. A realizable k-ε model is used to model the turbulent viscosity. The numerical results offer in-depth insight into physical understanding of many complex phenomena which cannot be fully explained experimentally.