Temperature dependence of liquid viscosity

The temperature dependence of liquid viscosity is the phenomenon by which liquid viscosity tends to decrease (or, alternatively, its fluidity tends to increase) as its temperature increases. This can be observed, for example, by watching how cooking oil appears to move more fluidly upon a frying pan after being heated by a stove.

Physical causes

A molecular view of liquids can be used for a qualitative picture of the process of decrease in the shear (or bulk) viscosity of a simple fluid with temperature. As the temperature increases, the time of interaction between neighbouring molecules of a liquid decreases because of the increased velocities of individual molecules. The macroscopic effect is that the intermolecular force appears to decrease and so does the bulk (or shear) viscosity. The actual process can be quite complex and is typically represented by simplified mathematical or empirical models, some of which are discussed below.[1] The models are valid over limited temperature ranges and for selected materials.

Models for shear viscosity

Exponential Model

An exponential model for the temperature-dependence of shear viscosity (μ) was first proposed by Reynolds in 1886.[2]

\mu(T)\,=\,\mu_0 \exp(-bT)

where T is temperature and \mu_0 and b are coefficients. See first-order fluid and second-order fluid. This is an empirical model that usually works for a limited range of temperatures.

Arrhenius model

The model is based on the assumption that the fluid flow obeys the Arrhenius equation for molecular kinetics:

\mu(T)\,=\,\mu_0 \exp\left( \frac {E}{RT} \right)

where T is temperature, \mu_0 is a coefficient, E is the activation energy and R is the universal gas constant. A first-order fluid is another name for a power-law fluid with exponential dependence of viscosity on temperature.

Williams-Landel-Ferry model

The Williams-Landel-Ferry model, or WLF for short, is usually used for polymer melts or other fluids that have a glass transition temperature.

The model is:

\mu(T)\,=\,\mu_0 \exp \left( \frac {-C_1 (T-T_r)} {C_2+ T -T_r} \right)

where T-temperature, C_1, C_2, T_r and \mu_0 are empiric parameters (only three of them are independent from each other).

If one selects the parameter T_r based on the glass transition temperature, then the parameters C_1, C_2 become very similar for the wide class of polymers. Typically, if T_r is set to match the glass transition temperature T_g, we get

C_1 \approx17.44

and

C_2 \approx 51.6 K.

Van Krevelen recommends to choose

T_r\,=\,T_g+43 K, then
C_1 \approx 8.86

and

C_2 \approx101.6 K.

Using such universal parameters allows one to guess the temperature dependence of a polymer by knowing the viscosity at a single temperature.

In reality the universal parameters are not that universal, and it is much better to fit the WLF parameters from the experimental data.

Masuko and Magill model

The model is usually used for polymer melts or other fluids that have a glass transition temperature as well as the WLF model. Ordinarily, The WLF model is limited to the temperature interval between Tg and Tg+ 100 K, But this model can be applied to more wide temperature range. [3]

The model is:

\log (\eta/\eta_g) = A\left[\exp \left\{\frac{B(T_g-T)}{T}\right\}-1\right]

The A and B are empirical parameters that does not depend on the materials. The average values are:

A = 14.25 to 16.24,

and

B = 5.34 to 7.60.

Viscosity of water

Viscosity of water equation accurate to within 2.5% from 0 °C to 370 °C:


 \mu(T) = 2.414\times10^{-5} \times 10^{247.8/(T - 140)} 
[4]

where T has units of Kelvin, and μ has units of Kg·s/m².

Models for kinematic viscosity

The effect of temperature on the kinematic viscosity (ν) has also been described by a number of empirical equations.

Walther formula

The Walther formula[1] is typically written in the form


  \log_{10}[\log_{10}(\nu+\lambda)] = A - B\,\log_{10}(T)

where λ is a shift constant, and A, B are empirical parameters.

Wright model

The Wright model[1] has the form


  \log_{10}[\log_{10}[\nu + \lambda + f(\nu)]] = A - B\,\log_{10}(T)

where an addition function f(ν), often a polynomial fit to experimental data, has been added to the Walther formula.

Seeton model

The Seeton model[1]is based on curve fitting the viscosity dependence of many liquids (refrigerants, hydrocarbons and lubricants) versus temperature and applies over a large temperature and viscosity range:

\ln \left( {\ln \left( {\nu  + 0.7 + e^{ - \nu } K_0 \left( {\nu  + 1.244067} \right)} \right)} \right) = A - B*\ln \left( T \right)

where T is absolute temperature in kelvins, \nu is the kinematic viscosity in centistokes, K_0 is the zero order modified Bessel function of the second kind, and A and B are liquid specific values. This form should not be applied to ammonia or water viscosity over a large temperature range.

For liquid metal viscosity as a function of temperature, Seeton proposed:

\ln \left( {\ln \left( {\nu  + 0.7 + e^{ - \nu } K_0 \left( {\nu  + 1.244067} \right)} \right)} \right) = A - {B \over T}

Notes

  1. 1.0 1.1 1.2 1.3 Seeton, Christopher J. (2006), "Viscosity-temperature correlation for liquids", Tribology Letters, doi:10.1007/s11249-006-9071-2
  2. Reynolds O. (1886). Phil Trans Royal Soc London, v. 177, p.157.
  3. Toru Masuko, Joseph H. Magill (1988), A comprehensive expression for temperature dependence of liquid viscosity., Journal of the Society of Rheology Japan, Vol.16
  4. Al-Shemmeri, Tarik (2012). Engineering Fluid Mechnanics. Ventus Publishing ApS. pp. 17–18. ISBN 978-87-403-0114-4.

See also