# Viscosity II: Gas Viscosity

*Michael Fowler*

### Gas Between Parallel Plates

Suppose now we repeat Newton’s
suggested experiment, the two parallel plates with one at rest the other moving
at steady speed, but with *gas* rather
than liquid between the plates.

It is found experimentally that the equation $F/A=\eta v/d$ still describes the force necessary to maintain steady motion, but, not surprisingly, for gases anywhere near atmospheric pressure the coefficient of viscosity is far lower than that for liquids (not counting liquid helium$\u2014$a special case):

Gas | Viscosity in 10^{-6} mPa.s |
---|---|

Air at 100K | 7.1 |

Air at 300K | 18.6 |

Air at 400K | 23.1 |

Hydrogen at 300K | 9.0 |

Helium at 300K | 20.0 |

Oxygen at 300K | 20.8 |

Xenon at 300K | 23.2 |

These values are from the CRC
Handbook, 85^{th} Edition, **6**-201.

Note first that, in contrast to the liquid case, gas
viscosity* increases* with temperature.
Even more surprising, it is found experimentally that over a very wide range of
densities, gas viscosity is* independent
of the density* of the gas!

Returning to the two plates, and picturing the gas between as made up of layers moving at different speeds as before.

The first thing to realize is that at atmospheric pressure
the molecules take up something like *a
thousandth* of the volume of the gas.
The previous picture of crowded molecules jostling each other is completely
irrelevant! As we shall discuss in more
detail later, the molecules of air at room temperature fly around at about 500
meters per second, the molecules have diameter around 0.35 nm, are around 3 or
4 nm apart on average, and travel of order 70 nm between collisions with other
molecules.

So where does the gas viscosity come from? Think of two adjacent layers of gas moving at
different speeds. Molecules from the fast layer fly into the slow layer, where
after a collision or two they are slowed to go along with the rest. At the same time, some slower molecules fly
into the fast layer. Even if we assume
that kinetic energy is conserved in each individual molecular collision (so
we’re ignoring for the moment excitation of internal modes of the molecules)
the *macroscopic* kinetic energy of the
layers of gas decreases overall (see the exercise in the preceding section). How can that be? Isn’t energy conserved? Yes, total energy *is* conserved, what happens is that some of the macroscopic kinetic
energy of the gas moving as a smooth substance has been transferred into the
microscopic kinetic energy of the individual molecules moving in random
directions within the gas, in other words, into the random molecular kinetic
energy we call heat.

### Estimating the Coefficient of Viscosity for a Gas: Momentum Transfer and Mean Free Path

The way to find the viscosity of a gas is to calculate the rate of $z$-direction (downward) transfer of $x$ -momentum, as explained in the previous section but one.

The moving top plate maintains a steady horizontal flow pattern, the $x$-direction speed at height $z$

$v\left(z\right)={v}_{0}z/d.$

As explained above, the moving plate is feeding $x$-direction momentum into the gas at a rate $F/A=\eta {v}_{0}/d,$ this momentum moves down through the gas at a steady rate, and the coefficient of viscosity tells us what this rate of momentum flow is for a given velocity gradient.

In fact, this rate of momentum flow can be calculated from a simple kinetic picture of the gas: remember the molecules have about a thousand times more room than they do in the liquid state, so the molecules go (relatively) a long way between collisions. We shall examine this kinetic picture of a gas in much more detail later in the course, but for now we’ll make the simplifying assumption that the molecules all have speed $u,$ and travel a distance $l$ between collisions. Actually this approximation gets us pretty close to the truth.

We take the density of molecules to be $n,$ the molecular mass $m.$ To begin thinking about $x$-direction momentum moving downwards, imagine some plane parallel to the plates and between them. Molecules from above are shooting through this plane and colliding with molecules in the slower moving gas below, on average transferring a little extra momentum in the $x$-direction to the slower stream. At the same time, some molecules from the slower stream are shooting upwards and will slow down the faster stream.

Let’s consider first the molecules passing through the
imaginary plane from above: we’re only interested in the molecules already moving
*downwards*, that’s half of them, so a
molecular density of $n/2.$ If we
assume for simplicity that all the molecules move at the same speed $u,$ then the average *downward* speed of these molecules $\overline{{u}_{z}}=u/2$ ( ${u}_{z}=u\mathrm{cos}\theta ,$ and the average value of $\mathrm{cos}\theta $ over all downward pointing directions is ${\scriptscriptstyle \frac{1}{2}}.$ )

Thus the number of molecules per second passing through the plane from above is

$$\frac{n}{2}\overline{\left|{u}_{z}\right|}=\frac{nu}{4}$$

and the same number are of course coming *up* from below! (Not shown.) (To see what’s going on, the mean free path
shown here is hugely exaggerated compared with the distance between the plates,
of course.)

But this isn’t quite what we want: we need to know how *efficiently* these molecules crossing the
plane are transferring momentum for the fast moving streams above to the slower
ones below. Consider one particular
molecule going from the faster stream at some downward angle into the slower
stream.

Let us assume it travels a distance $l$ from its last collision in the “fast” stream
to its first in the “slow” stream. The
average distance between collisions is called the *mean free path*, here “mean” is used in the sense of “average”, and we'll
write the mean free path* **$l.$** * We're simplifying slightly by taking all
distances between collisions to be *l*,
so we don’t bother with statistical averaging of the distance traveled. This
does not make a big difference. The
distance traveled *in the downward
direction* is $\Delta z=l{u}_{z}/u,$,
so the ( $x$-direction) speed difference between the two
streams is

$$\Delta v=\frac{dv\left(z\right)}{dz}\Delta z=\frac{{v}_{0}}{d}\cdot \frac{l{u}_{z}}{u}.$$

The molecule has mass $m,$ so *on
average* the momentum transferred from the fast stream to the slow stream is
$\Delta p=m\Delta v.$ With our simplifying assumption that all
molecules have the same speed $u,$ all downward values of ${u}_{z}$ between 0 and *u* are equally likely, and the density of downward-moving molecules
is $n/2,$ so the rate of transfer of momentum by
downward-moving molecules through the plane is

$$\frac{n}{2}\overline{{u}_{z}m\Delta v}=\frac{nm}{2}\cdot l\cdot \frac{\overline{{u}_{z}^{2}}}{u}\cdot \frac{{v}_{0}}{d}.$$

At the same time, molecules are moving upwards form the slower streams into the faster ones, and the calculation is exactly the same. These two processes have the same sign: in the first case, the slower stream is gaining forward momentum from the faster, in the second, the faster stream is losing forward momentum, and in both cases total forward momentum is conserved. Therefore, the two processes make the same contribution, and the total momentum flow rate (per unit area) across the plane is

$$\text{momentumtransferrate}=\text{}nm\cdot l\cdot \frac{\overline{{u}_{z}^{2}}}{u}\cdot \frac{{v}_{0}}{d}=\frac{1}{3}nmlu\cdot \frac{{v}_{0}}{d}$$

using ${u}_{z}^{2}={u}^{2}/3.$

Evidently this rate of downward transfer of $x$-direction momentum doesn’t depend on what level between the plates we choose for our imaginary plane, it’s the same momentum flow all the way from the top plate to the bottom plate: so it’s simply the rate at which the moving top plate is supplying $x$-direction momentum to the fluid,

$\text{momentumsupplyrate}=F/A=\eta {v}_{0}/d.$

Since the momentum supplied moves steadily downwards through
the fluid, the supply rate *is* the
transfer rate, the two equations above are for the same thing, and we deduce
that the coefficient of viscosity

$\eta =\frac{1}{3}nmlu.$

It should be noted that the $u$ in the above equation is the mean speed of the
atoms or molecules, which is about 8% greater than the rms speed appearing in
the formula for gas pressure. In deriving
the formula, we *did* make the
simplifying assumption that all the molecules move at the same speed, and had
the same mean free path, but Maxwell showed in 1860 that a more careful
treatment gives essentially the same result.
As discussed below, this result was historically important because when
put together with other measurements, in particular the van der Waals isotherms for real gases, *it led to estimates of the sizes of atoms*. A great deal of effort has been expended
improving Maxwell’s formula. Maxwell
assumed that after a collision, the molecules recoiled in random directions. In fact, a fast-moving molecule will tend to
be deflected less, so fast-moving downward molecules will have effectively
longer mean free paths than slow-moving ones, and this will substantially
increase the viscosity over that predicted using the “slow” mean free
path. Further progress can be made if
the interaction between molecules is represented by a definite potential, so
the scattering can be analyzed in detail.
From measurements of viscosity as a function of density, it is then
possible to learn something about the actual intermolecular potential.

### Why Doesn’t the Viscosity of a Gas Depend on Density?

Imagine we have a gas made up of equal numbers of red and
green molecules, which have the same size, mass, etc. One of the molecules
traveling through will on average have half its collisions with red molecules,
half with green. If now all the red
molecules suddenly disappear, the collision rate for our wandering molecule
will be halved. This means its mean free
path $l$ will double.
So, since the coefficient of viscosity
$\eta ={\scriptscriptstyle \frac{1}{3}}nmlu,$ halving the gas number density $n$ at the same time doubles the mean free path $l,$ so $\eta $ *is
unchanged*. ( $\eta $ does finally begin to drop when there is so
little gas left that the mean free path is of order the size of the container.)

Another way to see this is to think about the molecules shooting down from the faster stream into slower streams in the two-plate scenario. If the density is halved, there will only be half the molecules moving down, but each will deliver the $x$-momentum difference twice as far$\u2014$the further they go, the bigger the $x$-velocity difference between where they begin and where they end, and the more effective they are in transporting $x$-momentum downwards.

### Comparing the Viscosity Formula with Experiment

We’ll repeat the earlier table here for convenience:

Gas | Viscosity in 10^{-6} mPa.s |
---|---|

Air at 100K | 7.1 |

Air at 300K | 18.6 |

Air at 400K | 23.1 |

Hydrogen at 300K | 9.0 |

Helium at 300K | 20.0 |

Oxygen at 300K | 20.8 |

Xenon at 300K | 23.2 |

These values are from the CRC
Handbook, 85^{th} Edition, **6**-201.

It’s easy to see one reason why the viscosity increases with temperature: from $\eta =\left(1/3\right)nmlu,$ $\eta $ is proportional to the average molecular speed $u,$ and since this depends on temperature as ${\scriptscriptstyle \frac{1}{2}}m\overline{{u}^{2}}={\scriptscriptstyle \frac{3}{2}}{k}_{B}T$ (if this is unfamiliar, be assured we’ll be discussing it in detail later), this factor contributes a $\sqrt{T}$ dependence. In fact, though, from the table above, the increase in viscosity with temperature is more rapid than $\sqrt{T}.$. We know the density and mass remain constant (we’re far from relativistic energies!) so if the analysis is correct, the mean free path must also be increasing with temperature. In fact this is what happens$\u2014$many of the changes of molecular direction in flight are not caused by hard collisions with other molecules, but by longer range attractive forces (van der Waals forces) when one molecule simply passes reasonably close to another. Now these attractive forces obviously act for a shorter time on a faster molecule, so it is deflected less. This means that as temperature and molecular speed increase, the molecules get further in approximately the same direction, and therefore transport momentum more effectively.

Comparing viscosities of different gases, at the same temperature and pressure they will have the same number density $n.$ (Recall that a mole of any gas has Avogadro’s number of molecules, $6\times {10}^{23},$ and that at standard temperature and pressure this occupies 22.4 liters, for any gas, if tiny pressure and volume corrections for molecular attraction and molecular volume are ignored.) The molecules will also have the same kinetic energy since they’re at the same temperature (see previous paragraph).

Comparing hydrogen with oxygen, for example, the molecular
mass $m$ is up by a factor 16 but the velocity $u$ decreases by a factor 4 (since the molecular
kinetic energies ${\scriptscriptstyle \frac{1}{2}}m{v}^{2}$ are the same at the same temperature). In fact (see the table above), the viscosity
of oxygen is only *twice* that of
hydrogen, so from $\eta ={\scriptscriptstyle \frac{1}{3}}nmlu$ we conclude that the hydrogen molecule has
twice the mean free path distance $l$ between collisions$\u2014$not
surprising since it is smaller. Also,
helium must have an even longer mean free path, again not surprising for the
smallest molecule in existence. The
quite large difference between nitrogen and oxygen, next to each other in the
periodic table, is because N_{2} has a trivalent bonding, tighter than
the divalent O_{2} bonding, and in fact the internuclear distance in
the N_{2} molecule is 10% less than in O_{2}. Xenon is heavy (atomic weight 131), but its
mean free path is shorter than the others because the atom is substantially larger,
and so an easier target.

Historically, the sizes of many atoms and molecules were first estimated from viscosity measurements using this method. In fact, just such a table of atomic and molecular diameters, calculated on the assumption that the atoms or molecules are hard spheres, can be found in the CRC tables. But the true picture is more complicated: the atoms are not just hard spheres, as mentioned earlier they have van der Waals attractive forces between them, beyond the outermost shell. Actually, the electronic densities of atoms and molecules can now be found fairly precisely by quantum calculations, and the resulting “radii” are in reasonable agreement with those deduced from viscosity (there is no obvious natural definition of radius for the electron cloud).

### More General Laminar Flow Velocity Distributions

We’ve analyzed a particularly simple case: for the fluid
between two parallel plates, the bottom plate at rest and the top moving at
steady speed ${v}_{0},$ the fluid stream velocity increases *linearly* from zero at the bottom to ${v}_{0}$ at the top.
For more realistic laminar flow situations, such as that away from the
banks in a wide river, or flow down a pipe, the rate of velocity increase on
going from the fixed boundary (river bed or pipe surface) into the fluid is no
longer linear, that is, $dv\left(z\right)/dz$ is not the same everywhere.

The key to analyzing these more general laminar flow patterns is to find the forces acting on a small area of one of the “sheets”. (Or a larger area of sheet if the flow is uniform in the appropriate direction.) There will be external forces such as gravity or pressure maintaining the flow, which must balance the viscous drag forces exerted by neighboring sheets if the fluid is not accelerating. The sheet-sheet drag force is equivalent to the force $F/A=\eta {v}_{0}/d$ exerted by the top “sheet” of fluid on the top plate in the previous discussion, except that now the forces are within the fluid. As we discussed earlier, the drag force on the top plate comes from molecules close to the plate bouncing off, gaining momentum, which is subsequently transferred to other molecules a little further away. For liquids, this mechanism involves distances of order a few molecular diameters, for gases a few mean free paths. In either case, the distances are tiny on a macroscopic scale. This means that the appropriate formula for the drag force on a plate, or that between one sheet of fluid on another, is

$$\frac{F}{A}=\eta \frac{dv\left(z\right)}{dz}.$$

The *rate of change of
stream velocity close to the interface* determines the drag force.

We shall show in the next lecture how this formula can be used to determine the flow pattern in a river, and that in a circular pipe.