Bulk Viscosity

Volume viscosity (also called bulk viscosity, or second viscosity or, dilatational viscosity) is a material property relevant for characterizing fluid flow. Common symbols are $\zeta, \mu', \mu_\mathrm, \kappa$ or $\xi$. It has dimensions (mass / (length × time)), and the corresponding SI unit is the pascal-second (Pa·s).

Like other material properties (e.g. density, shear viscosity, and thermal conductivity) the value of volume viscosity is specific to each fluid and depends additionally on the fluid state, particularly its temperature and pressure. Physically, volume viscosity represents the irreversible resistance, over and above the reversible resistance caused by isentropic bulk modulus, to a compression or expansion of a fluid. At the molecular level, it stems from the finite time required for energy injected in the system to be distributed among the rotational and vibrational degrees of freedom of molecular motion.

Knowledge of the volume viscosity is important for understanding a variety of fluid phenomena, including sound attenuation in polyatomic gases (e.g. Stokes's law), propagation of shock waves, and dynamics of liquids containing gas bubbles. In many fluid dynamics problems, however, its effect can be neglected. For instance, it is 0 in a monatomic gas at low density (unless the gas is moderately relativistic), whereas in an incompressible flow the volume viscosity is superfluous since it does not appear in the equation of motion.

Volume viscosity was introduced in 1879 by Sir Horace Lamb in his famous work ''Hydrodynamics''.Lamb, H., "Hydrodynamics", Sixth Edition,''Dover Publications'', NY (1932) Although relatively obscure in the scientific literature at large, volume viscosity is discussed in depth in many important works on fluid mechanics,Potter, M.C., Wiggert, D.C. "Mechaniscs of Fluids", ''Prentics Hall'', NJ (1997) fluid acoustics,Morse, P.M. and Ingard, K.U. "Theoretical Acoustics", ''Princeton University Press''(1968)Temkin, S., "Elements of Acoustics", ''John Wiley and Sons'', NY (1981) theory of liquids,Kirkwood, J.G., Buff, F.P., Green, M.S., "The statistical mechanical theory of transport processes. 3. The coefficients of shear and bulk viscosity in liquids", J. Chemical Physics, 17, 10, 988-994, (1949)Enskog, D. "Kungliga Svenska Vetenskapsakademiens Handlingar", 63, 4, (1922) rheology,Graves, R.E. and Argrow, B.M. "Bulk viscosity: Past to Present", ''Journal of Thermophysics and Heat Transfer'',13, 3, 337–342 (1999) and relativistic hydrodynamics.

Derivation and use

At thermodynamic equilibrium, the negative-one-third of the trace of the Cauchy stress tensor is often identified with the thermodynamic pressure,

:$-\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_$.

Volume viscosity appears in the classic Navier-Stokes equation if it is written for compressible fluid, as described in most books on general hydrodynamicsHappel, J. and Brenner, H. "Low Reynolds number hydrodynamics", ''Prentice-Hall'', (1965)Landau, L.D. and Lifshitz, E.M. "Fluid mechanics", ''Pergamon Press'', New York (1959) and acoustics.Litovitz, T.A. and Davis, C.M. In "Physical Acoustics", Ed. W.P.Mason, vol. 2, chapter 5, ''Academic Press'', NY, (1964)Dukhin, A. S. and Goetz, P. J. ''Characterization of liquids, nano- and micro- particulates and porous bodies using Ultrasound'', Elsevier, 2017

:\rho \frac = -\nabla P + \nabla\cdot\left mu\left(\nabla\mathbf + \left(\nabla\mathbf\right)^T - \frac (\nabla\cdot\mathbf)\mathbf\right) \right+ \nabla\cdot zeta(\nabla\cdot \mathbf)\mathbf+ \rho \mathbf

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/Dt$ is the material derivative. By introducing the tensors (matrices) $\boldsymbol$, $\boldsymbol$ and $e \mathbf$ (where ''e'' is a scalar called dilation, and $\mathbf$ 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 tensorsee also Generalized Newtonian fluid) and compression flow (i.e. the isotropic dilation tensor), respectively,

:$\boldsymbol = \frac \left( \nabla\mathbf + \left(\nabla\mathbf\right)^T \right)$

:$e = \frac \nabla \! \cdot \! \mathbf$

:$\boldsymbol = \boldsymbol - e \mathbf$

the classic Navier-Stokes equation gets a lucid form.

:$\rho \frac = -\nabla (P - 3 \zeta e) + \nabla\cdot ( 2\mu \boldsymbol \gamma) + \rho \mathbf$

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 =0$

So the incompressible Navier-Stokes equation can be simply written:

:$\rho \frac = -\nabla P + \nabla\cdot ( 2\mu \boldsymbol \epsilon) + \rho \mathbf$

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

According to Landau,

He later adds:

After an example, he concludes (with $\zeta$ used to represent volume viscosity):

Measurement

A brief review of the techniques available for measuring the volume viscosity of liquids can be found in Dukhin & Goetz and Sharma (2019). One such method is by using an acoustic rheometer.

Below are values of the volume viscosity for several Newtonian liquids at 25 °C (reported in cP):
methanol - 0.8
ethanol - 1.4
propanol - 2.7
pentanol - 2.8
acetone - 1.4
toluene - 7.6
cyclohexanone - 7.0
hexane - 2.4

Recent studies have determined the volume viscosity for a variety of gases, including carbon dioxide, methane, and nitrous oxide. These were found to have volume viscosities which were hundreds to thousands of times larger than their shear viscosities. Fluids having large volume viscosities include those used as working fluids in power systems having non-fossil fuel heat sources, wind tunnel testing, and pharmaceutical processing.

Modeling

There are many publications dedicated to numerical modeling of volume viscosity. A detailed review of these studies can be found in Sharma (2019)Sharma, B and Kumar, R "Estimation of bulk viscosity of dilute gases using a nonequilibrium molecular dynamics approach.", ''Physical Review E'',100, 013309 (2019) and Cramer.Cramer, M.S. "Numerical estimates for the bulk viscosity of ideal gases.", ''Phys. Fluids'',24, 066102 (2012) In the latter study, a number of common fluids were found to have bulk viscosities which were hundreds to thousands of times larger than their shear viscosities. For relativistic liquids and gases, bulk viscosity is conveniently modeled in terms of a mathematical duality with chemically reacting relativistic fluids.

References

*
*
{{refend

Colloidal chemistry
Fluid dynamics
Viscosity