We present a fully practical residual-based adaptive simulation framework for two-phase flows with variable densities governed by a Cahn-Hilliard Navier-Stokes model with double obstacle potential.The method is based on a new stable time integration scheme,which conserves the energy decay.In the adaptive framework the generation of discrete energy contributions stemming from coarsening is avoided by an appropriate modification of the coarsening strategy.