burian
- 64
- 6
I am thrown of by the tensor notation xD. I'll learn it and reattend your reply some day.vanhees71 said:Sure, that was only the ideal-fluid case. For real fluids you have to take into account "friction", i.e., dissipation, and that's of course also part of the stress. There you need a difference in speed between the fluid in the volume and outside of it along the boundary of the volume. In the next approximation thus the corresponding stress from friction must be of first order in the gradient of ##\vec{v}##. So it should be proportional to the components ##\partial_j v_k##. Now the antisymmetric part of this tensor refers to a rotation of the fluid element as a whole and thus doesn't describe true velocity differences. So it's the symmetric part which must enter the first-order contribution to the stress. Since in a liquid or gas locally no direction is preferred there can only be two pieces, corresponding to the irreducible representations of the rotation group. For a symmetric 2nd rank tensor this are the five components of the trace-less part as well as the trace part. Thus one has
$$\sigma_{jk}^{(1)}=\eta (\partial_j v_k + \partial_k v_j - \frac{2}{3} \partial_{i} v_i \delta_{jk} )+\zeta \partial_i v_i \delta_{jk}.$$
As it must be the tensor only depends on two scalar quantities ##\eta##, and ##\zeta##, the shear and bulk viscosities.
Then you end up with the Navier-Stokes equation,
$$\rho (\partial_t \vec{v} + \vec{v} \cdot \vec{\nabla} \vec{v})=-\vec{\nabla} P + \eta \Delta \vec{v} + \left (\frac{\eta}{3}+\zeta \right) \vec{\nabla}(\vec{\nabla} \cdot \vec{v}),$$
where I have assumed the ##\eta## and ##\zeta## can be considered constant within the fluid (though they depend on temperature and density and thus are not strictly constant).