Energy stable and conservative numerical schemes for simulating two-phase flows using Cahn-Hilliard Navier Stokes equations