Abstract: The key idea of the phase-field methodology is to replace sharp interfaces by thin transition regions where interfacial forces are smoothly distributed. One of the main reasons for the success of the phase-field methodology is that it is based on rigorous mathematics and thermodynamics. Most phase-field models satisfy a nonlinear stability relationship, usually expressed as a time-decreasing free-energy functional, which is called energy gradient stability. To obtain accurate numerical solutions of such problems, it is desirable to use high-order approximations in space and time. Yet because of the difficulties introduced by the combination of nonlinearity and stiffness, most computations have been limited to lower-order in time, and in most cases to constant time-stepping. On the other hand, numerical simulations for many physical problems require large time integration; as a result large time-stepping methods become necessary in some time regimes. To effectively solve the relevant physical problems, the combination of higher-order and highly stable temporal discretizations becomes very useful. In this talk, we will describe some adaptive time stepping approach for phase-field problems, which inherits the energy gradient stability of the continuous model. Particular attention will be given to effectively resolve both the solution dynamics and steady state solutions.