Cauchy Momentum Equation - Stress Tensor I've been trying to understand the derivation for the Cauchy Momentum Equation for so long now, and there is one part that every derivation glides over very quickly with practically no explanation (I'm guessing they assume the reader knows it already). 
The part I'm stuck with is how they relate stress tensor, $\sigma_{ji}$, to the sum of the forces' on a infinitesimal block of volume $dV$. I'll give you a bit of the context of the situation I'm in. Here's how this part of every derivation goes.

Suppose you have a differential/infinitesimal [rectangular prism] volume of fluid $dV$, side-lengths $dx_j$, density $\rho$, and acceleration in the $i^{\textrm{th}}$ direction $a_i$. Applying Newton's second law per unit volume in the $i^{\textrm{th}}$ direction gives us
$$\rho\,a_i = \sum F_i, \,\,\,\,\textrm{$\sum F_i$ is the net body force in the $i^{\textrm{th}}$ direction}$$
Now for the net body force in the $i^{\textrm{th}}$ direction, we have the external body forces $f$ and the stress [surface really] forces, which is given as the rate of the stress variation in the $i^{\textrm{th}}$ direction, $\sum \limits_{j} \frac{\partial \sigma_{ji}}{\partial x_j}$. Thus, we have by Newton's second law per unit volume of fluid in the $i^{\textrm{th}}$ direction,
$$\rho\,a_i=f_i + \sum \limits_{j} \frac{\partial \sigma_{ji}}{\partial x_j}$$
And so the total force (all of the previous equation multiplied by the volume of the unit of fluid which it pertained to) is given as
$$\rho\,a_i\,dV = f_i\,dV + \sum \limits_{j} \frac{\partial \sigma_{ji}}{\partial x_j}dV$$

How on earth is the total stress force in a particular direction (used as a body force) given as the rate of the stress variation in that direction?! I see how the units work out, but I can't see any logic behind it. I thought that $\sigma_{ji}$ represented the stress (force per unit area) on the $dx_j$ side pointing in the $i$ direction. If that's the case, how do the relate that surface force to being a body force, especially in the way that is stated above (saying it's the rate of stress variation in that direction)?
Please help me. I've already asked this on the physics stackexchange, but nobody answered me.
 A: Consider a small volume element
$$\{(s_1,s_2,s_3)\in \mathbb{R}^3: x_1 \leq s_1 \leq x_1 + \Delta x_1,x_2 \leq s_2 \leq x_2 + \Delta x_2,x_3 \leq s_3 \leq x_3 + \Delta x_3\}.$$
The surface (stress) force acting on that element in the $i$ direction consists of contributions from each of three pairs of faces where coordinates $s_1$, $s_2$ and $s_3$ are constant, respectively.
Consider one pair of where $s_1$ is fixed: $s_1 = x_1$ and $s_1 = x_1 + \Delta x_1$. The net force in the $i$ direction is
$$F_{i,1}\approx \sum_{j=1}^{3}\sigma_{ij}n_j|_{s_1=x_1 + \Delta x_1}\Delta x_2\Delta x_3+\sum_{j=1}^{3}\sigma_{ij}n_j|_{s_1=x_1}\Delta x_2\Delta x_3$$
This is approximate because we evaluate $\sigma_{ij}(s_1,s_2,s_3)$ at representative coordinates $s_2= x_2$ and $s_3=x_3$ . In general, we assume that the stress tensor is a continuously differentiable function of spatial position.
Here we use the fact that the surface force (in terms of the stress tensor and normal vector) is the product $\mathbb{\sigma}\cdot\mathbb{n}.$ On the face $s_1 = x_1 + \Delta x_1$ we have $\mathbb{n} = (1,0,0)$, and on the face $s_1 = x_1$ we have $\mathbb{n} = (-1,0,0).$
Hence,
$$F_{i,1}\approx \sigma_{i1}|_{s_1=x_1 + \Delta x_1}\Delta x_2\Delta x_3-\sigma_{i1}|_{s_1=x_1}\Delta x_2\Delta x_3$$
The contribution to the force per unit volume is
$$\frac{F_{i,1}}{\Delta x_1 \Delta x_2 \Delta x_3}\approx \frac{\sigma_{i1}|_{s_1=x_1 + \Delta x_1}-\sigma_{i1}|_{s_1=x_1}}{\Delta x_1}.$$
Now consider the limit as the element shrinks to the point $(x_1,x_2,x_3)$ and $\Delta x_k \rightarrow 0.$
Then the contribution to the net surface force per unit volume at that point is
$$\hat{F_{i,1}}= \lim_{\Delta x_1 \rightarrow 0}\frac{\sigma_{i1}|_{s_1=x_1 + \Delta x_1}-\sigma_{i1}|_{s_1=x_1}}{\Delta x_1}=\frac{\partial \sigma_{i1}}{\partial x_1}.$$
Applying the same derivation to the other pairs of faces we obtain the total contribution to the net surface force per unit volume (in the $i$ direction) as
$$\hat{F_{i}}=\sum_{j=1}^{3}\frac{\partial \sigma_{ij}}{\partial x_j}$$
