The interfacial dynamics of immiscible and incompressible two-phase fluids is of great interest but computationally not easy to resolve. In this talk, we will focus on a phase field description of the interfacial dynamics. The model is a set of coupling equations, which consists of the Navier-Stokes equations and the Cahn-Hilliard equation. The talk starts with a review of the existing methods for numerical solutions of the Navier-Stokes-Cahn-Hilliard coupling equations. Then we propose and analyze a class of efficient time-stepping schemes for the model. A detailed comparison with existing schemes will be presented, and the advantage of the new schemes are highlighted.