Agfdhyk

Viscosity
A simulation of liquids with different viscosities. The liquid on the right has higher viscosity than the liquid on the left.
Common symbols	η, μ
Derivations from other quantities	μ = G·t

The viscosity of a fluid is the measure of its resistance to gradual deformation by shear stress or tensile stress.^[1] For liquids, it corresponds to the informal concept of "thickness": for example, syrup has a higher viscosity than water.^[2]

Viscosity is the property of a fluid which opposes the relative motion between two surfaces of the fluid that are moving at different velocities. In simple terms, viscosity means friction between the molecules of fluid. When the fluid is forced through a tube, the particles which compose the fluid generally move more quickly near the tube's axis and more slowly near its walls; therefore some stress (such as a pressure difference between the two ends of the tube) is needed to overcome the friction between particle layers to keep the fluid moving. For a given velocity pattern, the stress required is proportional to the fluid's viscosity.

A fluid that has no resistance to shear stress is known as an ideal or inviscid fluid. Zero viscosity is observed only at very low temperatures in superfluids. Otherwise, all fluids have positive viscosity and are technically said to be viscous or viscid. A fluid with a relatively high viscosity, such as pitch, may appear to be a solid.

Etymology

The word "viscosity" is derived from the Latin "viscum", meaning mistletoe and also a viscous glue made from mistletoe berries.^[3]

Definition

Simple definition

Illustration of a planar Couette flow. Since the shearing flow is opposed by friction between adjacent layers of fluid (which are in relative motion), a force is required to sustain the motion of the upper plate. The relative strength of this force is a measure of the fluid's viscosity.

In a general parallel flow, the shear stress is proportional to the gradient of the velocity.

The viscosity of a fluid expresses its resistance to shearing flows, in which adjacent fluid layers are in relative motion. A simple example of such a shearing flow is a planar Couette flow, where a fluid is trapped between two infinitely large plates, one fixed and one in parallel motion at constant speed u{displaystyle u} (see illustration to the right). Although viscosity applies to general flows, it is easy to define and visualize in a Couette flow.

If the speed of the top plate is low enough, then in steady state the fluid particles move parallel to it, and their speed varies from 0{displaystyle 0} at the bottom to u at the top. Each layer of fluid moves faster than the one just below it, and friction between them gives rise to a force resisting their relative motion. In particular, the fluid applies on the top plate a force in the direction opposite to its motion, and an equal but opposite force on the bottom plate. An external force is therefore required in order to keep the top plate moving at constant speed.

In many fluids, the flow velocity is observed to vary linearly from zero at the bottom to u{displaystyle u} at the top. Moreover, the magnitude F of the force acting on the top plate is found to be proportional to the speed u and the area A of each plate, and inversely proportional to their separation y:

F=μAuy.{displaystyle F=mu A{frac {u}{y}}.}

The proportionality factor μ is the viscosity of the fluid, with units of Pa⋅s{displaystyle {text{Pa}}cdot {text{s}}} (pascal-second). The ratio u/y is called the rate of shear deformation or shear velocity, and is the derivative of the fluid speed in the direction perpendicular to the plates (see illustrations to the right). If the velocity does not vary linearly with y, then the appropriate generalization is

τ=μ∂u∂y,{displaystyle tau =mu {frac {partial u}{partial y}},}

where τ = F/A, and ∂u/∂y is the local shear velocity. This expression is referred to as Newton's law of viscosity. In shearing flows with planar symmetry, it is what defines μ{displaystyle mu }. It is a special case of the general definition of viscosity (see below), which can be expressed in coordinate-free form.

Use of the Greek letter mu (μ) for the viscosity is common among mechanical and chemical engineers, as well as physicists.^[4][5][6] However, the Greek letter eta (η) is also used by chemists, physicists, and the IUPAC.^[7] The viscosity μ{displaystyle mu } is sometimes also referred to as the shear viscosity. However, at least one author discourages the use of this terminology, noting that μ{displaystyle mu } can be appear in nonshearing flows in addition to shearing flows.^[8]

General definition

In general, the stresses within a flow can be attributed partly to the deformation of the material from some rest state (elastic stress), and partly to the rate of change of the deformation over time (viscous stress). In a fluid, by definition, the elastic stress includes only the hydrostatic pressure. In very general terms, the fluid's viscosity is the relation between the strain rate and the viscous stress. In the Newtonian fluid model the relationship is by definition a linear map, described by a viscosity tensor that, when multiplied by the strain rate tensor (which is the gradient of the flow's velocity), gives the viscous stress tensor. In Cartesian coordinates, this gives^[9]

τij=∑k∑lμijkl∂vk∂rl.{displaystyle tau _{ij}=sum _{k}sum _{l}mu _{ijkl}{frac {partial v_{k}}{partial r_{l}}}.}

Since the indices in the above expression can vary from 1 to 3, there are 81 "viscosity coefficients" μijkl{displaystyle mu _{ijkl}} in total. However, due to spatial symmetries these coefficients are not all independent. For instance, for isotropic Newtonian fluids, the 81 coefficients can be reduced to 2 independent parameters. The most usual decomposition yields the standard (scalar) viscosity μ and the bulk viscosity κ{displaystyle kappa }:^[10]

τ=μ[∇v+(∇v)†]−(23μ−κ)(∇⋅v)δ,{displaystyle mathbf {tau } =mu left[nabla mathbf {v} +(nabla mathbf {v} )^{dagger }right]-left({frac {2}{3}}mu -kappa right)(nabla cdot mathbf {v} )mathbf {delta } ,}

where δ{displaystyle mathbf {delta } } is the unit tensor, and the dagger †{displaystyle dagger } denotes the transpose. This equation can be thought of as a generalized form of Newton's law of viscosity.

The bulk viscosity expresses a type of internal friction that resists the shearless compression or expansion of a fluid. Knowledge of κ{displaystyle kappa } is frequently not necessary in fluid dynamics problems. For example, incompressible liquids satisfy ∇⋅v=0{displaystyle nabla cdot mathbf {v} =0} and so the term containing κ{displaystyle kappa } is absent. Moreover, κ{displaystyle kappa } is often assumed to be negligible for gases since it is 0{displaystyle 0} in a monoatomic ideal gas.^[10] One situation in which κ{displaystyle kappa } can be important is the calculation of energy loss in sound and shock waves, described by Stokes' law of sound attenuation, since these phenomena involve rapid expansions and compressions.

Dynamic and kinematic viscosity

In fluid dynamics, it is common to work in terms of the kinematic viscosity (also called "momentum diffusivity"), defined as the ratio of the viscosity μ to the density of the fluid ρ. It is usually denoted by the Greek letter nu (ν) and has units (length)2/time{displaystyle mathrm {(length)^{2}/time} }:

ν=μρ{displaystyle nu ={frac {mu }{rho }}}.

Consistent with this nomenclature, the viscosity μ{displaystyle mu } is frequently called the dynamic viscosity.

Momentum transport

Transport theory provides an alternate interpretation of viscosity in terms of momentum transport: viscosity is the material property which characterizes momentum transport within a fluid, just as thermal conductivity characterizes heat transport, and (mass) diffusivity characterizes mass transport.^[11] To see this, note that in Newton's law of viscosity, τ=μ(∂u/∂y){displaystyle tau =mu (partial u/partial y)}, the shear stress τ{displaystyle tau } has units equivalent to a momentum flux, i.e. momentum per unit time per unit area. Thus, τ{displaystyle tau } can be interpreted as specifying the flow of momentum in the y{displaystyle y} direction from one fluid layer to the next. Per Newton's law of viscosity, this momentum flow occurs across a velocity gradient, and the magnitude of the corresponding momentum flux is determined by the viscosity.

The analogy with heat and mass transfer can be made explicit. Just as heat flows from high temperature to low temperature and mass flows from high density to low density, momentum flows from high velocity to low velocity. These behaviors are all described by compact expressions, called constitutive relations, whose one-dimensional forms are given here:

J=−D∂ρ∂x(Fick's law of diffusion)q=−kt∂T∂x(Fourier's law of heat conduction)τ=μ∂u∂y(Newton's law of viscosity){displaystyle {begin{aligned}mathbf {J} &=-D{frac {partial rho }{partial x}}qquad ;;;,{text{(Fick's law of diffusion)}}\mathbf {q} &=-k_{t}{frac {partial T}{partial x}}qquad ;;,{text{(Fourier's law of heat conduction)}}\tau &=mu {frac {partial u}{partial y}}qquad qquad {text{(Newton's law of viscosity)}}\end{aligned}}}

where ρ{displaystyle rho } is the density, J{displaystyle mathbf {J} } and q{displaystyle mathbf {q} } are the mass and heat fluxes, and D{displaystyle D} and kt{displaystyle k_{t}} are the mass diffusivity and thermal conductivity.^[12]

The fact that mass, momentum, and energy (heat) transport are among the most relevant processes in continuum mechanics is not a coincidence: these are among the few physical quantities that are conserved at the microscopic level in interparticle collisions. Thus, rather than being dictated by the fast and complex microscopic interaction timescale, their dynamics occurs on macroscopic timescales, as described by the various equations of transport theory and hydrodynamics.

Newtonian and non-Newtonian fluids

Viscosity, the slope of each line, varies among materials.

Newton's law of viscosity is a constitutive equation (like Hooke's law, Fick's law, and Ohm's law): it is not a fundamental law of nature but an approximation that holds in some materials and fails in others.

A fluid that behaves according to Newton's law, with a viscosity μ that is independent of the stress, is said to be Newtonian. Gases, water, and many common liquids can be considered Newtonian in ordinary conditions and contexts. There are many non-Newtonian fluids that significantly deviate from that law in some way or other. For example:

• 
  Shear-thickening liquids, whose viscosity increases with the rate of shear strain.


• 
  Shear-thinning liquids, whose viscosity decreases with the rate of shear strain.


• 
  Thixotropic liquids, that become less viscous over time when shaken, agitated, or otherwise stressed.


• 
  Rheopectic (dilatant) liquids, that become more viscous over time when shaken, agitated, or otherwise stressed.


• 
  Bingham plastics that behave as a solid at low stresses but flow as a viscous fluid at high stresses.

Shear-thinning liquids are very commonly, but misleadingly, described as thixotropic.

Even for a Newtonian fluid, the viscosity usually depends on its composition and temperature. For gases and other compressible fluids, it depends on temperature and varies very slowly with pressure.

The viscosity of some fluids may depend on other factors. A magnetorheological fluid, for example, becomes thicker when subjected to a magnetic field, possibly to the point of behaving like a solid.

In solids

The viscous forces that arise during fluid flow must not be confused with the elastic forces that arise in a solid in response to shear, compression or extension stresses. While in the latter the stress is proportional to the amount of shear deformation, in a fluid it is proportional to the rate of deformation over time. (For this reason, Maxwell used the term fugitive elasticity for fluid viscosity.)

However, many liquids (including water) will briefly react like elastic solids when subjected to sudden stress. Conversely, many "solids" (even granite) will flow like liquids, albeit very slowly, even under arbitrarily small stress.^[13] Such materials are therefore best described as possessing both elasticity (reaction to deformation) and viscosity (reaction to rate of deformation); that is, being viscoelastic.

Indeed, some authors have claimed that amorphous solids, such as glass and many polymers, are actually liquids with a very high viscosity (greater than 10¹² Pa·s).
^[14] However, other authors dispute this hypothesis, claiming instead that there is some threshold for the stress, below which most solids will not flow at all,^[15] and that alleged instances of glass flow in window panes of old buildings are due to the crude manufacturing process of older eras rather than to the viscosity of glass.^[16]

Viscoelastic solids may exhibit both shear viscosity and bulk viscosity. The extensional viscosity is a linear combination of the shear and bulk viscosities that describes the reaction of a solid elastic material to elongation. It is widely used for characterizing polymers.

In geology, earth materials that exhibit viscous deformation at least three orders of magnitude greater than their elastic deformation are sometimes called rheids.^[17]

Measurement

Viscosity is measured with various types of viscometers and rheometers. A rheometer is used for those fluids that cannot be defined by a single value of viscosity and therefore require more parameters to be set and measured than is the case for a viscometer. Close temperature control of the fluid is essential to acquire accurate measurements, particularly in materials like lubricants, whose viscosity can double with a change of only 5 °C.

For some fluids, the viscosity is constant over a wide range of shear rates (Newtonian fluids). The fluids without a constant viscosity (non-Newtonian fluids) cannot be described by a single number. Non-Newtonian fluids exhibit a variety of different correlations between shear stress and shear rate.

One of the most common instruments for measuring kinematic viscosity is the glass capillary viscometer.

In coating industries, viscosity may be measured with a cup in which the efflux time is measured. There are several sorts of cup – such as the Zahn cup and the Ford viscosity cup – with the usage of each type varying mainly according to the industry. The efflux time can also be converted to kinematic viscosities (centistokes, cSt) through the conversion equations.^[18]

Also used in coatings, a Stormer viscometer uses load-based rotation in order to determine viscosity. The viscosity is reported in Krebs units (KU), which are unique to Stormer viscometers.

Vibrating viscometers can also be used to measure viscosity. Resonant, or vibrational viscometers work by creating shear waves within the liquid. In this method, the sensor is submerged in the fluid and is made to resonate at a specific frequency. As the surface of the sensor shears through the liquid, energy is lost due to its viscosity. This dissipated energy is then measured and converted into a viscosity reading. A higher viscosity causes a greater loss of energy.^{[citation needed]}

Extensional viscosity can be measured with various rheometers that apply extensional stress.

Volume viscosity can be measured with an acoustic rheometer.

Apparent viscosity is a calculation derived from tests performed on drilling fluid used in oil or gas well development. These calculations and tests help engineers develop and maintain the properties of the drilling fluid to the specifications required.

Units

The SI unit of dynamic viscosity is Pa·s or kg·m⁻¹·s⁻¹. The cgs unit is called the poise^[19] (P), named after Jean Léonard Marie Poiseuille. It is commonly expressed, particularly in ASTM standards, as centipoise (cP) since the latter is equal to the SI multiple millipascal seconds (mPa·s).

The SI unit of kinematic viscosity is m²/s, whereas the cgs unit for kinematic viscosity is the stokes (St), named after Sir George Gabriel Stokes.^[20] It is sometimes expressed in terms of centistokes (cSt). In U.S. usage, stoke is sometimes used as the singular form.

The reciprocal of viscosity is fluidity, usually symbolized by φ = 1/μ or F = 1/μ, depending on the convention used, measured in reciprocal poise (P⁻¹, or cm·s·g⁻¹), sometimes called the rhe. Fluidity is seldom used in engineering practice.

Nonstandard units include the reyn, a British unit of dynamic viscosity.^{[citation needed]} In the automotive industry the viscosity index is used to describe the change of viscosity with temperature.

At one time the petroleum industry relied on measuring kinematic viscosity by means of the Saybolt viscometer, and expressing kinematic viscosity in units of Saybolt universal seconds (SUS).^[21] Other abbreviations such as SSU (Saybolt seconds universal) or SUV (Saybolt universal viscosity) are sometimes used. Kinematic viscosity in centistokes can be converted from SUS according to the arithmetic and the reference table provided in ASTM D 2161.^[22]

Molecular origins

In general, the viscosity of a system depends in detail on how the molecules constituting the system interact. There are no simple but correct expressions for the viscosity of a fluid. The simplest exact expressions are the Green–Kubo relations for the linear shear viscosity or the transient time correlation function expressions derived by Evans and Morriss in 1985.^[23] Although these expressions are each exact, calculating the viscosity of a dense fluid using these relations currently requires the use of molecular dynamics computer simulations. On the other hand, much more progress can be made for a dilute gas. Even elementary assumptions about how gas molecules move and interact lead to a basic understanding of the molecular origins of viscosity. More sophisticated treatments can be constructed by systematically coarse-graining the equations of motion of the gas molecules. An example of such a treatment is Chapman–Enskog theory, which derives expressions for the viscosity of a dilute gas from the Boltzmann equation.^[24]

Momentum transport in gases is generally mediated by discrete molecular collisions, and in liquids by attractive forces which bind molecules close together.^[25] Because of this, the dynamic viscosities of liquids are typically much larger than those of gases.

Gases

Elementary calculation of viscosity for a dilute gas

Consider a dilute gas moving parallel to the x{displaystyle x}-axis with velocity u(y){displaystyle u(y)} that depends only on the y{displaystyle y} coordinate. To simplify the discussion, the gas is assumed to have uniform temperature and density. Under these assumptions, the x{displaystyle x} velocity of a molecule passing through y=0{displaystyle y=0} is equal to whatever velocity that molecule had when its mean free path λ{displaystyle lambda } began. Because λ{displaystyle lambda } is typically small compared with macroscopic scales, the average x{displaystyle x} velocity of such a molecule has the form u(0)±αλdudy(0){displaystyle u(0)pm alpha lambda {frac {du}{dy}}(0)}, where α{displaystyle alpha } is a numerical constant on the order of 1{displaystyle 1}. (Some authors estimate α=2/3{displaystyle alpha =2/3};[25][26] on the other hand, a more careful calculation for rigid elastic spheres gives α≃0.998{displaystyle alpha simeq 0.998}.) Now, because half the molecules on either side are moving towards y=0{displaystyle y=0}, and doing so on average with half the average moleculer speed (8kBT/πm)1/2{displaystyle (8k_{text{B}}T/pi m)^{1/2}}, the momentum flux from either side is 14ρ⋅8kBTπm⋅(u(0)±αλdudy(0)).{displaystyle {frac {1}{4}}rho cdot {sqrt {frac {8k_{text{B}}T}{pi m}}}cdot left(u(0)pm alpha lambda {frac {du}{dy}}(0)right).} The net momentum flux at y=0{displaystyle y=0} is the difference of the two: −12ρ⋅8kBTπm⋅αλdudy(0).{displaystyle -{frac {1}{2}}rho cdot {sqrt {frac {8k_{text{B}}T}{pi m}}}cdot alpha lambda {frac {du}{dy}}(0).} According to the definition of viscosity, this momentum flux should be equal to −μdudy(0){displaystyle -mu {frac {du}{dy}}(0)}, which leads to μ=αρλ2kBTπm.{displaystyle mu =alpha rho lambda {sqrt {frac {2k_{text{B}}T}{pi m}}}.}

Viscosity in gases arises principally from the molecular diffusion that transports momentum between layers of flow. An elementary calculation for a dilute gas at temperature T{displaystyle T} and density ρ{displaystyle rho } gives

μ=αρλ2kBTπm,{displaystyle mu =alpha rho lambda {sqrt {frac {2k_{text{B}}T}{pi m}}},}

where kB{displaystyle k_{text{B}}} is the Boltzmann constant, m{displaystyle m} the molecular mass, and α{displaystyle alpha } a numerical constant on the order of 1{displaystyle 1}. The quantity λ{displaystyle lambda }, the mean free path, measures the average distance a molecule travels between collisions. Even without a priori knowledge of α{displaystyle alpha }, this expression has interesting implications. In particular, since λ{displaystyle lambda } is typically inversely proportional to density and increases with temperature, μ{displaystyle mu } itself should increase with temperature and be independent of density. In fact, both of these predictions persist in more sophisticated treatments, and accurately describe experimental observations. Note that this behavior runs counter to common intuition regarding liquids, for which viscosity typically decreases with temperature.^[25][26]

For rigid elastic spheres of diameter σ{displaystyle sigma }, λ{displaystyle lambda } can be computed, giving

μ=απ3/2kBmTσ2.{displaystyle mu ={frac {alpha }{pi ^{3/2}}}{frac {sqrt {k_{text{B}}mT}}{sigma ^{2}}}.}

In this case λ{displaystyle lambda } is independent of temperature, so μ∝T1/2{displaystyle mu propto T^{1/2}}. For more complicated molecular models, however, λ{displaystyle lambda } depends on temperature in a non-trivial way, and simple kinetic arguments as used here are inadequate. More fundamentally, the notion of a mean free path becomes imprecise for particles that interact over a finite range, which limits the usefulness of the concept for describing real-world gases.^[27]

Chapman–Enskog theory

A technique developed by Sydney Chapman and David Enskog in the early 1900s allows a more refined calculation of μ{displaystyle mu }.^[24] It is based on the Boltzmann equation, which provides a systematic statistical description of a dilute gas in terms of intermolecular interactions.^[28] As such, their technique allows accurate calculation of μ{displaystyle mu } for more realistic molecular models, such as those incorporating intermolecular attraction rather than just hard-core repulsion.

It turns out that a more realistic modeling of interactions is essential for accurate prediction of the temperature dependence of μ{displaystyle mu }, which experiments show increases more rapidly than the T1/2{displaystyle T^{1/2}} trend predicted for rigid elastic spheres.^[25] Indeed, the Chapman–Enskog analysis shows that the predicted temperature dependence can be tuned by varying the parameters in various molecular models. A simple example is the Sutherland model,^[29] which describes rigid elastic spheres with weak mutual attraction. In such a case, the attractive force can be treated perturbatively, which leads to a particularly simple expression for μ{displaystyle mu }:

μ=516σ2(kBmTπ)1/2(1+ST)−1,{displaystyle mu ={frac {5}{16sigma ^{2}}}left({frac {k_{text{B}}mT}{pi }}right)^{1/2}left(1+{frac {S}{T}}right)^{-1},}

where S{displaystyle S} is independent of temperature, being determined only by the parameters of the intermolecular attraction. To connect with experiment, it is convenient to rewrite as

μ=μ0(TT0)3/2T0+ST+S,{displaystyle mu =mu _{0}left({frac {T}{T_{0}}}right)^{3/2}{frac {T_{0}+S}{T+S}},}

where μ0{displaystyle mu _{0}} is the viscosity at temperature T0{displaystyle T_{0}}. If μ{displaystyle mu } is known from experiments at T=T0{displaystyle T=T_{0}} and at least one other temperature, then S{displaystyle S} can be calculated. It turns out that expressions for μ{displaystyle mu } obtained in this way are accurate for a number of gases over a sizable range of temperatures. On the other hand, Chapman and Cowling^[24] argue that this success does not imply that molecules actually interact according to the Sutherland model. Rather, they interpret the prediction for μ{displaystyle mu } as a simple interpolation which is valid for some gases over fixed ranges of temperature, but otherwise does not provide a picture of intermolecular interactions which is fundamentally correct and general. Slightly more sophisticated models, such as the Lennard–Jones potential, may provide a better picture, but only at the cost of a more opaque dependence on temperature. In some systems the assumption of spherical symmetry must be abandoned as well, as is the case for vapors with highly polar molecules like H₂O.^[30][31]

Liquids

Play media

Video showing three liquids with different viscosities

Play media

Experiment showing the behavior of a viscous fluid with blue dye for visibility

In contrast with gases, there is no simple yet accurate picture for the molecular origins of viscosity in liquids.

At the simplest level of description, the relative motion of adjacent layers in a liquid is opposed primarily by attractive molecular forces
acting across the layer boundary. In this picture, one (correctly) expects viscosity to decrease with increasing temperature. This is because
increasing temperature increases the random thermal motion of the molecules, which makes it easier for them to overcome their attractive interactions.^[32]

Building on this visualization, a simple theory can be constructed in analogy with the discrete structure of a solid: groups of molecules in a liquid
are visualized as forming "cages" which surround and enclose single molecules.^[33] These cages can be occupied or unoccupied, and
stronger molecular attraction corresponds to stronger cages.
Due to random thermal motion, a molecule "hops" between cages at a rate which varies inversely with the strength of molecular attractions. In equilibrium these "hops" are not biased in any direction.
On the other hand, in order for two adjacent layers to move relative to each other, the "hops" must be biased in the direction
of the relative motion. The force required to sustain this directed motion can be estimated for a given shear rate, leading to

μ≈NAhVexp⁡(3.8TbT),{displaystyle mu approx {frac {N_{A}h}{V}}operatorname {exp} left(3.8{frac {T_{b}}{T}}right),}

(1)

where NA{displaystyle N_{A}} is the Avogadro constant, h{displaystyle h} is the Planck constant, V{displaystyle V} is the volume of a mole of liquid, and Tb{displaystyle T_{b}} is the normal boiling point. This result has the same form as the widespread and accurate empirical relation

μ=AeBT,{displaystyle mu =Ae^{frac {B}{T}},}

(2)

where A{displaystyle A} and B{displaystyle B} are constants fit from data.^[33][34] One the other hand, several authors express caution with respect to this model.
Errors as large as 30% can be encountered using equation (1), compared with fitting equation (2) to experimental data.^[33] More fundamentally, the physical assumptions underlying equation (1) have been extensively criticized.^[35] It has also been argued that the exponential dependence in equation (1) does not necessarily describe experimental observations more accurately than simpler, non-exponential expressions.^[36][37]

In light of these shortcomings, the development of a less ad-hoc model is a matter of practical interest.
Foregoing simplicity in favor of precision, it is possible to write rigorous expressions for viscosity starting from the fundamental equations of motion for molecules. A classic example
of this approach is Irving-Kirkwood theory.^[38] On the other hand, such expressions
are given as averages over multiparticle correlation functions and are therefore difficult to apply in practice.

In general, empirically derived expressions (based on existing viscosity measurements) appear to be the only consistently reliable means of calculating viscosity in liquids.^[39]

Mixtures, blends, and suspensions

Gaseous mixtures

The same molecular-kinetic picture of a single component gas can also be applied to a gaseous mixture. For instance, in the Chapman-Enskog approach the viscosity μmix{displaystyle mu _{text{mix}}} of a binary mixture of gases can be written in terms of the individual component viscosities μ1,2{displaystyle mu _{1,2}}, their respective volume fractions, and the intermolecular interactions.^[40] As for the single-component gas, the dependence of μmix{displaystyle mu _{text{mix}}} on the parameters of the intermolecular interactions enters through various collisional integrals which may not be expressible in terms of elementary functions. To obtain usable expressions for μmix{displaystyle mu _{text{mix}}} which reasonably match experimental data, the collisional integrals typically must be evaluated using some combination of analytic calculation and empirical fitting. An example of such a procedure is the Sutherland approach for the single-component gas, discussed above.

Blends of liquids

As for pure liquids, the viscosity of a blend of liquids is difficult to predict from molecular principles. One method is to extend the molecular "cage" theory presented above for a pure liquid. This can be done with varying levels of sophistication. One useful expression resulting from such an analysis is the Lederer-Roegiers equation for a binary mixture:

μblend=x1x1+αx2ln⁡μ1+αx2x1+αx2ln⁡μ2,{displaystyle mu _{text{blend}}={frac {x_{1}}{x_{1}+alpha x_{2}}}ln mu _{1}+{frac {alpha x_{2}}{x_{1}+alpha x_{2}}}ln mu _{2},}

where α{displaystyle alpha } is an empirical parameter, and x1,2{displaystyle x_{1,2}} and μ1,2{displaystyle mu _{1,2}} are the respective mole fractions and viscosities of the component liquids.^[41]

Since blending is an important process in the lubricating and oil industries, a variety of empirical or propriety equations exist for predicting the viscosity of a blend, besides those stemming directly from molecular theory.^[41]

Suspensions

In a suspension of solid particles (e.g. micron-size spheres suspended in oil), an effective viscosity μeff{displaystyle mu _{text{eff}}} can be defined in terms of stress and strain components which are averaged over a volume large compared with the distance between the suspended particles, but small with respect to macroscopic dimensions.^[42] Such suspensions generally exhibit non-Newtonian behavior. However, for dilute systems in steady flows, the behavior is Newtonian and expressions for μeff{displaystyle mu _{text{eff}}} can be derived directly from the particle dynamics. In a very dilute system, with volume fraction ϕ≲0.02{displaystyle phi lesssim 0.02}, interactions between the suspended particles can be ignored. In such a case one can explicitly calculate the flow field around each particle independently, and combine the results to obtain μeff{displaystyle mu _{text{eff}}}. For spheres, this results in the Einstein equation:

μeff=μ0(1+52ϕ),{displaystyle mu _{text{eff}}=mu _{0}left(1+{frac {5}{2}}phi right),}

where μ0{displaystyle mu _{0}} is the viscosity of the suspending liquid. The linear dependence on ϕ{displaystyle phi } is a direct consequence of neglecting interparticle interactions; in general, one will have:

μeff=μ0(1+Bϕ),{displaystyle mu _{text{eff}}=mu _{0}left(1+Bphi right),}

where the coefficient B{displaystyle B} may depend on the particle shape (e.g. spheres, rods, disks).^[43] Perhaps surprisingly, in the century since its introduction the coefficient for spheres, B=5/2{displaystyle B=5/2}, has not been conclusively pinned down by experiments, with various experiments finding values in the range 1.5≲B≲5{displaystyle 1.5lesssim Blesssim 5}.^[44] This deficiency has been attributed to difficulty in controlling experimental conditions.^[45]

In denser suspensions, μeff{displaystyle mu _{text{eff}}} acquires a nonlinear dependence on ϕ{displaystyle phi }, which indicates the importance of interparticle interactions. Various analytical and semi-empirical schemes exist for capturing this regime. At the most basic level, a term quadratic in ϕ{displaystyle phi } is added to μeff{displaystyle mu _{text{eff}}}:

μeff=μ0(1+Bϕ+B1ϕ2),{displaystyle mu _{text{eff}}=mu _{0}left(1+Bphi +B_{1}phi ^{2}right),}

and the coefficient B1{displaystyle B_{1}} is fit from experimental data or approximated from the microscopic theory. In general, however, one should be cautious in applying such simple formulas since non-Newtonian behavior appears in dense suspensions (ϕ≳0.25{displaystyle phi gtrsim 0.25} for spheres),^[44] or in suspensions of elongated or flexible particles.^[46]

A distinction must be made between a suspension of solid particles, described above, and an emulsion. The latter is a suspension of tiny droplets, which themselves may exhibit internal circulation. The presence of internal circulation can noticeably decrease the observed effective viscosity, and different theoretical or semi-empirical models must be used.^[47]

Amorphous materials

Common glass viscosity curves^[48]

In the high and low temperature limits, viscous flow in amorphous materials (e.g. in glasses and melts)^[49][50][51] has the Arrhenius form:

μ=AeQRT,{displaystyle mu =Ae^{frac {Q}{RT}},}

where Q is a relevant activation energy, given in terms of molecular parameters; T is temperature; R is the molar gas constant; and A is approximately a constant. The activation energy Q takes a different value depending on whether the high or low temperature limit is being considered: it changes from a high value Q_H at low temperatures (in the glassy state) to a low value Q_L at high temperatures (in the liquid state).

Common logarithm of viscosity against temperature for B₂O₃, showing two regimes

For intermediate temperatures, Q{displaystyle Q} varies nontrivially with temperature and the simple Arrhenius form fails. On the other hand, the two-exponential equation

μ=ATexp⁡(BRT)[1+Cexp⁡(DRT)],{displaystyle mu =ATexp left({frac {B}{RT}}right)left[1+Cexp left({frac {D}{RT}}right)right],}

where A{displaystyle A}, B{displaystyle B}, C{displaystyle C}, D{displaystyle D} are all constants, provides a good fit to experimental data over the entire range of temperatures, while at the same time reducing to the correct Arrhenius form in the low and high temperature limits. Besides being a convenient fit to data, the expression can also be derived from various theoretical models of amorphous materials at the atomic level.^[51]

Eddy viscosity

In the study of turbulence in fluids, a common practical strategy for calculation is to ignore the small-scale vortices (or eddies) in the motion and to calculate a large-scale motion with an eddy viscosity that characterizes the transport and dissipation of energy in the smaller-scale flow (see large eddy simulation).^[52][53] Values of eddy viscosity used in modeling ocean circulation may be from 7004500000000000000♠5×10⁴ to 7006100000000000000♠1×10⁶ Pa·s depending upon the resolution of the numerical grid.^{[citation needed]}

Selected substances

In the University of Queensland pitch drop experiment, pitch has been dripping slowly through a funnel since 1927, at a rate of one drop roughly every decade. In this way the viscosity of pitch has been determined to be approximately 230 billion (7011229999999999999♠2.3×10¹¹) times that of water.^[54]

Observed values of viscosity vary over several orders of magnitude, even for common substances. For instance, a 70% sucrose (sugar) solution has a viscosity over 400 times that of water, and 26000 times that of air.^[55] More dramatically, pitch has been estimated to have a viscosity 230 billion times that of water.^[54]

Water

The viscosity of water is about 0.89 mPa·s at room temperature (25 °C). As a function of temperature, the viscosity can be estimated using the semi-empirical relation:

μ=A×10B/(T−C),{displaystyle mu =Atimes 10^{B/(T-C)},}

where A = 6995241400000000000♠2.414×10⁻⁵ Pa·s, B = 247.8 K, and C = 140 K.^{[citation needed]}

Experimentally determined values of the viscosity at various temperatures are given below.

Air

Under standard atmospheric conditions (25 °C and pressure of 1 bar), the viscosity of air is 18.5 μPa·s, roughly 50 times smaller than the viscosity of water at the same temperature. Except at very high pressure, the viscosity of air depends mostly on the temperature.

Other common substances

Honey being drizzled

Substance	Viscosity (mPa·s)	Temperature (°C)
Whole milk[56]	2.12	20
Olive oil[56]	56.2	26
Honey[57]	≈{displaystyle approx } 2000-10000	20
Ketchup[a][58]	≈{displaystyle approx } 5000-20000	25
Peanut butter[a][59]	≈{displaystyle approx } 104-106
Pitch[54]	7011229999999999999♠2.3×1011	10-30 (variable)

See also

References

Further reading

• 
  Hatschek, Emil (1928). The Viscosity of Liquids. New York: Van Nostrand. OCLC 53438464.
  


• 
  Massey, B. S.; Ward-Smith, A. J. (2011). Mechanics of Fluids (9th ed.). London & New York: Spon Press. ISBN 978-0-415-60259-4. OCLC 690084654.
  

External links

Look up viscosity in Wiktionary, the free dictionary.

Wikisource has the text of The New Student's Reference Work article Viscosity of Liquids.

• 
  Fluid properties - high accuracy calculation of viscosity for frequently encountered pure liquids and gases


• Gas viscosity calculator as function of temperature


• Air viscosity calculator as function of temperature and pressure


• 
  Fluid Characteristics Chart - a table of viscosities and vapor pressures for various fluids


• 
  Gas Dynamics Toolbox - calculate coefficient of viscosity for mixtures of gases


• 
  Glass Viscosity Measurement - viscosity measurement, viscosity units and fixpoints, glass viscosity calculation


• 
  Kinematic Viscosity - conversion between kinematic and dynamic viscosity


• 
  Physical Characteristics of Water - a table of water viscosity as a function of temperature


• Vogel–Tammann–Fulcher Equation Parameters


• Calculation of temperature-dependent dynamic viscosities for some common components


• 
  "Test Procedures for Testing Highway and Nonroad Engines and Omnibus Technical Amendments" - United States Environmental Protection Agency
  


• Artificial viscosity


• Viscosity of Air, Dynamic and Kinematic, Engineers Edge