Numerical simulations of unsteady gas flows are studied on the basis of Gas-Kinetic Unified Algorithm (GKUA) from rarefied transition to continuum flow regimes. Several typical examples are adopted. An unsteady flow solver is developed by solving the Boltzmann model equations, including the Shakhov model and the Rykov model etc. The Rykov kinetic equation involving the effect of rotational energy can be transformed into two kinetic governing equations with inelastic and elastic collisions by integrating the molecular velocity distribution function with the weight factor on the energy of rotational motion. Then, the reduced velocity distribution functions are devised to further simplify the governing equation for one- and twodimensional flows. The simultaneous equations are numerically solved by the discrete velocity ordinate (DVO) method in velocity space and the finite-difference schemes in physical space. The time-explicit operator-splitting scheme is constructed, and numerical stability conditions to ascertain the time step are discussed. As the application of the newly developed GKUA, several unsteady varying processes of one- and twodimensional flows with different Knudsen number are simulated, and the unsteady transport phenomena and rarefied effects are revealed and analyzed. It is validated that the GKUA solver is competent for simulations of unsteady gas dynamics covering various flow regimes.