The motion of noncolloidal particles convected by a nonhomogeneous and nonstationary viscous fluid flow is investigated, assuming that inertial effects can be neglected. It appears that the particle volumetric flux is the sum of a convective part, V(E)phi, and a diffusive term, -D-self.delphi, where V-E is the Eulerian mean velocity of a test particle, phi is the particle volume fraction, and D-self is the coefficient of self-diffusion. The latter measures the local temporal growth of the mean square displacement of a tracer particle from its average position and can be written as D-self(r) = D(r,r), where the cross-diffusion tensor D(r(1),r(2)) is the time integral of the velocity cross-correlation function. On the other hand, the Eulerian mean velocity V-E is the sum of the coarse-grained average particle velocity, (V) over bar, and a drift velocity, V-d(r) = -[(partial derivative/partial derivativer(2)).D-T(r,r(2))](r2) = r. This last term, which is identically zero when the suspended particles are passive tracers, indicates that the suspended particles tend to move toward regions with smaller diffusivities. This result demonstrates that the motion of each suspended particle is a random process satisfying a generalized nonlinear Langevin equation, where the fluctuating term is described through the cross diffusivity D. (C) 2003 American Institute of Physics.