naive question on Boltzmann equation and conservation laws The Boltzmann equation in absence of external force reads:
$\frac{\partial f}{\partial t} + \vec{v} \cdot \frac{\partial f}{\partial \vec{r}} = \left( \frac{\partial f}{\partial t}\right)_{coll}$
Where the r.h.s. stands for the change in the distribution function of the velocities owing to collisions. I won't specify here a particular collision kernel as it is not necessary for my question...I think.
It is common to introduce 


*

*The particle density $n(\vec{r},t) \equiv \int d^3v \:f(\vec{r}, \vec{v},t)$

*The mean velocity field $\vec{u}(\vec{r},t) \equiv \frac{1}{n(\vec{r},t)}\int d^3 v \:\vec{v} f(\vec{r}, \vec{v},t)$
It is then easy to integrate over the velocities the force-free Boltzmann equation above and find the famous conservation law:
$\frac{\partial n(\vec{r},t)}{\partial t} + \nabla_{\vec{r}}\cdot(n(\vec{r},t)\vec{u}(\vec{r},t)) = 0$
That is because, if the collisions conserve the number of particles, then there should not be any 'source' or 'sink' on the r.h.s.
My naive question is the following:
The above conservation law for matter can be recast as $\frac{\partial n(\vec{r},t)}{\partial t}+ \nabla_{\vec{r}}\cdot \vec{J}(\vec{r},t) = 0$ with the particle flux $\vec{J}(\vec{r},t) = n(\vec{r},t)\vec{u}(\vec{r},t)$.
This is fine but it does not seem very general. For instance, it does not seem to encompass a diffusive flux does it?
Hence, although this equation is very general, it does not seem to lead to Fick's laws of diffusion and yet I kind of remember that one can get the latter from some kind of expansion of at least a linearized Boltzmann equation...
Could somene unlighten me on this apparent paradox? 
 A: The confusion comes from the average flow velocity, the flux in general can be written as:

The previous expression is obviously for the flux, the first term which is the product of density and flow velocity is the most general definition of flux. Yet it can be expanded into two components (as in the left most part of my previous equation), a diffusive component described by Ficks law and a convective component.
The natural question now is, what is the difference between the flow velocity u and the convection velocity v. Physically speaking, the flow can be thought of being a result of two mechanisms, the first is Diffusion which is the spread of particles from higher density region to lower density region (driven by the internal energy of the system), have a look here. The second is Convection where particles are driven mechanically into new regions (driven by external forces exerted on the system). So the splitting into two components has a physical background as you can see. 
Mathematically speaking, the difference depends on your approach of solving for the total flux. The general way (solving Boltzmann equation is more general but physically impractical) is to solve for u by multiplying Boltzmann equation by momentum and integrating to obtain the 1st moment (the zero moment being the density equation you wrote in your question). If you did that you get an equation for u and you don’t need any further expansion of convection diffusion components. The 1st moment of Boltzmann equation yields Navier-Stokes equation. Not everyone is welling to solve it.
The alternative simpler less general approach is to determine v by some very simplified means (without solving Navier-Stokes equation) and represent the general flux as summation of a diffusive component and a convective component. A good example here is plasma, where the force driving convection is electrical, so convection velocity is determined by:

The convection velocity in this case is called the drift velocity, mu is the mobility of electrons and E is the electric field. So as you can see, in such a scenario one can solve for the density of a species without a need to solve Navier-Stokes equation for the average flow velocity u of that species. After solving, one can compute the average flow velocity by dividing the flux by the density.
(Addition)
Let us have a look at an infinitesimal volume element as described by Boltzmann equation. To make everything easier to explain for me I am going to assume there is only one dimension in space and velocity space. Since there is no force acting, the velocity distribution is isotropic.
Let us have case (a) as shown below where two volume elements have the same number of particles (uniform density). In each one of the two elements, the number of particles going to the left is equal to the number of particles going to the right (because the velocity distribution is isotropic). If you look at the interface between the particles and try to count the average flux between the two elements, it will be zero because 3 particles from the left box go the right and 3 particles from the right box go to the left, so the net momentum that crossed the interface is zero. Since the flux is zero and density is not zero, that means mathematically the average flow velocity is zero.

Compare that to case (b) where the number of particles is different (non-uniform density). On the interface there are 3 particles leaving from left box to right and 1.5 particles (I should have chosen an even number, my bad!!)  go from the right box to the left. The net flux crossing the interface is directed from left to right, i.e. from higher density to lower density region. If you divide the net flux by density you get an average flow velocity to the right.

Returning back to the first argument of splitting flux into diffusive (driven by internal energy) and convective (driven by external force) you see clearly it is a fully diffusive flux in case (b) because there is no external force applied to the system.  This net motion is driven by the thermal motion of particles, which is directly related to their internal energy.
A: In order to talk about diffusion you need at least two particle species (something has to diffuse relative to something else). Then the continuity equation refers to the total mass density of the fluid
$$
\rho = \int d\Gamma\, m_1f_1 + \int d\Gamma\, m_2 f_2
$$ 
and the fluid velocity is defined by
$$
\rho\vec{u} = \int d\Gamma\, m_1v_1f_1 + \int d\Gamma\, m_2v_2 f_2\, . 
$$ 
Mass conservation in the collision term ensures that $\partial_t\rho+\vec{\nabla}(\rho \vec{u})=0$ and no dissipative fluxes appear. Now we can define a concentration, e.g.
$$
c =  \int d\Gamma\, f_1 / \int d\Gamma\, f_2\, . 
$$
The quantity $\rho c$ satisfies a more complicated continuity equation, with a diffusive current $\vec{\jmath}$ on the right hand side. This current satisfies (to leading order in gradients) Fick's law,
$$
\vec{\jmath} = -\rho D\vec{\nabla} c + \ldots
$$
where the $\ldots$ are thermal diffusion terms, and higher order gradients.
