Volume viscosity explained

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',\mub,\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, 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.[1] Although relatively obscure in the scientific literature at large, volume viscosity is discussed in depth in many important works on fluid mechanics,[2] fluid acoustics,[3] [4] theory of liquids,[5] [6] and rheology.[7]

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,

a
-{1\over3}\sigma
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

\muv

.

Volume viscosity appears in the classic Navier-Stokes equation if it is written for compressible fluid, as described in most books on general hydrodynamics[8] [9] and acoustics.[10] [11]

\rho

Dv
Dt

=-\nablaP+\nabla\left[\mu\left(\nablav+\left(\nablav\right)T-

2
3

(\nablav)I\right)\right]+\nabla[\zeta(\nablav)I]+\rhog

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

Dv/Dt

is the material derivative. By introducing the tensors (matrices)

\boldsymbol{\epsilon}

,

\boldsymbol{\gamma}

and

eI

(where e is a scalar called dilation, and

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[12]) and compression flow (i.e. the isotropic dilation tensor), respectively,

\boldsymbol{\epsilon}=

1
2

\left(\nablav+\left(\nablav\right)T\right)

e=

1
3

\nablav

\boldsymbol{\gamma}=\boldsymbol{\epsilon}-eI

the classic Navier-Stokes equation gets a lucid form.

\rho

Dv
Dt

=-\nabla(P-3\zetae)+\nabla(2\mu\boldsymbol\gamma)+\rhog

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:

\nablav=0

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

\rho

Dv
Dt

=-\nablaP+\nabla(2\mu\boldsymbol\epsilon)+\rhog

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=

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,[9]

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[11] 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):[13] 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)[14] and Cramer.[15] 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.

References

Notes and References

  1. Lamb, H., "Hydrodynamics", Sixth Edition,Dover Publications, NY (1932)
  2. Potter, M.C., Wiggert, D.C. "Mechaniscs of Fluids", Prentics Hall, NJ (1997)
  3. Morse, P.M. and Ingard, K.U. "Theoretical Acoustics", Princeton University Press(1968)
  4. Temkin, S., "Elements of Acoustics", John Wiley and Sons, NY (1981)
  5. 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)
  6. Enskog, D. "Kungliga Svenska Vetenskapsakademiens Handlingar", 63, 4, (1922)
  7. Graves, R.E. and Argrow, B.M. "Bulk viscosity: Past to Present", Journal of Thermophysics and Heat Transfer,13, 3, 337–342 (1999)
  8. Happel, J. and Brenner, H. "Low Reynolds number hydrodynamics", Prentice-Hall, (1965)
  9. Landau, L.D. and Lifshitz, E.M. "Fluid mechanics", Pergamon Press, New York (1959)
  10. Litovitz, T.A. and Davis, C.M. In "Physical Acoustics", Ed. W.P.Mason, vol. 2, chapter 5, Academic Press, NY, (1964)
  11. Dukhin, A. S. and Goetz, P. J. Characterization of liquids, nano- and micro- particulates and porous bodies using Ultrasound, Elsevier, 2017
  12. see also Generalized Newtonian fluid
  13. Dukhin. Andrei S.. Goetz. Philip J.. Bulk viscosity and compressibility measurement using acoustic spectroscopy. The Journal of Chemical Physics. 130. 12. 2009. 124519. 0021-9606. 10.1063/1.3095471. 19334863. 2009JChPh.130l4519D .
  14. Sharma, B and Kumar, R "Estimation of bulk viscosity of dilute gases using a nonequilibrium molecular dynamics approach.", Physical Review E,100, 013309 (2019)
  15. Cramer, M.S. "Numerical estimates for the bulk viscosity of ideal gases.", Phys. Fluids,24, 066102 (2012)