Accurate time stepping of vesicle suspensions is crucial for obtaining stable and accurate simulations. I will describe two algorithms that can be used to accomplish this task. First, a deferred correction method uses a low-order semi-implicit time integrator to construct high-order results. Second, physically conserved quantities are used to adaptively select appropriate time step sizes. The methodology will be demonstrated for simulations involving stenosis, high concentrations, and adhesion.