At thermodynamic equilibrium, the negative-one-third of the
trace of the
Cauchy stress tensor is often identified with the thermodynamic
pressure, :-{1\over3}\sigma_a^a = P, which depends only on equilibrium state variables like temperature and density (
equation of state). In general, the trace of the stress tensor is the sum of thermodynamic pressure contribution and another contribution which is proportional to the
divergence of the velocity field. This coefficient of proportionality is called volume viscosity. Common symbols for volume viscosity are \zeta and \mu_{v}. Volume viscosity appears in the classic
Navier-Stokes equation if it is written for
compressible fluid, as described in most books on general hydrodynamics and acoustics. :\rho \frac{D \mathbf{v}}{Dt} = -\nabla P + \nabla\cdot\left[\mu\left(\nabla\mathbf{v} + \left(\nabla\mathbf{v}\right)^T - \frac{2}{3} (\nabla\cdot\mathbf{v})\mathbf{I}\right) \right] + \nabla\cdot[\zeta(\nabla\cdot \mathbf{v})\mathbf{I}] + \rho \mathbf{g} where \mu is the
shear viscosity coefficient and \zeta is the volume viscosity coefficient. The parameters \mu and \zeta were originally called the first and bulk viscosity coefficients, respectively. The operator D\mathbf{v}/Dt is
the material derivative. By introducing the tensors (matrices) \boldsymbol{\epsilon} , \boldsymbol{\gamma} and e \mathbf{I} (where
e is a scalar called
dilation, and \mathbf{I} is the
identity tensor), which describes crude shear flow (i.e. the
strain rate tensor), pure shear flow (i.e. the
deviatoric part of the strain rate tensor, i.e. the
shear rate tensor) and compression flow (i.e. the isotropic dilation tensor), respectively, : \boldsymbol{\epsilon} = \frac{1}{2} \left( \nabla\mathbf{v} + \left(\nabla\mathbf{v}\right)^T \right) : e = \frac{1}{3} \nabla \! \cdot \! \mathbf{v} : \boldsymbol{\gamma} = \boldsymbol{\epsilon} - e \mathbf{I} the classic Navier-Stokes equation gets a lucid form. :\rho \frac{D \mathbf{v}}{Dt} = -\nabla (P - 3 \zeta e) + \nabla\cdot ( 2\mu \boldsymbol \gamma) + \rho \mathbf{g} Note that the term in the momentum equation that contains the volume viscosity disappears for an
incompressible flow because there is no
divergence of the flow, and so also no flow dilation
e to which is proportional: : \nabla \! \cdot \! \mathbf{v} =0 So the incompressible Navier-Stokes equation can be simply written: :\rho \frac{D \mathbf{v}}{Dt} = -\nabla P + \nabla\cdot ( 2\mu \boldsymbol \epsilon) + \rho \mathbf{g} In fact, note that for the incompressible flow the strain rate is purely deviatoric since there is no dilation (
e=0). In other words, for an incompressible flow the isotropic stress component is simply the pressure: :p= \frac 1 3 Tr(\boldsymbol \sigma) and the deviatoric (
shear) stress is simply twice the product between the shear viscosity and the strain rate (
Newton's constitutive law): :\boldsymbol \tau = 2 \mu \boldsymbol \epsilon Therefore, in the incompressible flow the volume viscosity plays no role in the fluid dynamics. However, in a compressible flow there are cases where \zeta\gg\mu, which are explained below. In general, moreover, \zeta is not just a property of the fluid in the classic thermodynamic sense, but also depends on the process, for example the compression/expansion rate. The same goes for shear viscosity. For a
Newtonian fluid the shear viscosity is a pure fluid property, but for a
non-Newtonian fluid it is not a pure fluid property due to its dependence on the velocity gradient. Neither shear nor volume viscosity are equilibrium parameters or properties, but transport properties. The velocity gradient and/or compression rate are therefore independent variables together with pressure, temperature, and other
state variables. ==Landau's explanation==