A two-fluid model is used to study the time evolution of stellar winds including the dynamical effect of cosmic rays. Neglecting the diffusion of cosmic rays, we seek self-similar solutions to spherically symmetric winds with a termination shock. The velo
