  Dynamic critical behavior in superfluid systems is considered in a presence
of external stirring and advecting processes. The latter are generated by means
of the Gaussian random velocity ensemble with white-noise character in time
variable and self-similar spatial dependence. The main focus of this work is to
analyze an effect of compressible modes on the critical behavior. The model is
formulated through stochastic Langevin equations, which are then recast into
Janssen-De Dominicis response formalism. Employing the field-theoretic
perturbative renormalization group method we analyze large-scale properties of
the model. Explicit calculations are performed to the leading one-loop
approximation in the double $(\varepsilon, y)$ expansion scheme, where
$\varepsilon$ is a deviation from the upper critical dimension $d_c = 4$ and
$y$ describes a scaling properties of the velocity ensemble. Altogether five
distinct universality classes are expected to be macroscopically observable. In
contrast to the incompressible case, we found that compressibility leads to an
enhancement and stabilization of non-trivial asymptotic regimes.
