  We derive the hydrodynamical equations of r-mode oscillations in neutron
stars in presence of a novel damping mechanism related to particle number
changing processes. The change in the number densities of the various species
leads to new dissipative terms in the equations which are responsible of the
{\it rocket effect}. We employ a two-fluid model, with one fluid consisting of
the charged components, while the second fluid consists of superfluid neutrons.
We consider two different kind of r-mode oscillations, one associated with
comoving displacements, and the second one associated with countermoving, out
of phase, displacements.
