Equations of fluid dynamics

Appendix J Derivation of the stress tensor for a newtonian fluid

It is generally assumed, that friction between fluid elements is proportional to the area of their surfaces. So in general the frictional or vicous force on a fluid element can be expressed like

Fvisc,i=Aσijnj𝑑A=Vrjσij𝑑V. (J.1)

This force leads to an irreversible rise of temperature in the fluid or an irreversible decrease of kinetic energy expressed by the equation for the dissipation

tEkin,visc=-Vσijrjvi𝑑V (J.2)

For a motionless fluid (vi=0) and for a fluid with constant velocity (virj=0) this integral is zero. But also a rotating observer of a motionless fluid should not see a rise in the temperature of a fluid 3131We do not consider here the a rigidly rotating fluid as it os often done in the literature, because a rigidly rotating fluid is unphysical. that means

Vσijrjvi𝑑V=0 (J.3)

A rotating observer of a motionless fluid sees a velocity field of the form

vi=ϵijkωjrk (J.4)

where ωj is the angular velocity vector and rk is the position vector. It can be shown, that for such a velocity field the jacobian is antisymmetric, that means

virj=-vjri (J.5)

Using this and equation (H.5) in equation (J.3) we get

V12(σij-σji)rjvi𝑑V=0 (J.6)

This relation can only be fulfilled if the stress tensor σij is symmetric

σij=σji (J.7)

For a newtonian fluid is it assumed that the stress tensor is proportional only to the first derivatives of the velocity field. Together with the requirement of symmetry the most general for of such a tensor is

σij=a(vjri+virj)+bδijvkrk (J.8)

Usually the trace is split off the first term and added the second term so

σij=a(vjri+virj-23δijvkrk)+(2a3+b)δijvkrk (J.9)

Using the definitions 2a=η and 2a3+b=ζ we get the form most common in literature

σij=η[12(virj+vjri)-13δijvkrk]+ζδijvkrk (J.10)