# 1.7: Stokes’ Law

- Page ID
- 234

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\)

\( \newcommand{\id}{\mathrm{id}}\)

\( \newcommand{\Span}{\mathrm{span}}\)

\( \newcommand{\kernel}{\mathrm{null}\,}\)

\( \newcommand{\range}{\mathrm{range}\,}\)

\( \newcommand{\RealPart}{\mathrm{Re}}\)

\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

\( \newcommand{\Argument}{\mathrm{Arg}}\)

\( \newcommand{\norm}[1]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)

\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)

\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vectorC}[1]{\textbf{#1}} \)

\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)## Dropping the Ball (Slowly)

We’ve seen how viscosity acts as a frictional brake on the rate at which water flows through a pipe, let us now examine its frictional effect on an object falling through a viscous medium. To make it simple, we take a sphere. If we use a very viscous liquid, such as glycerin, and a small sphere, for example a ball bearing of radius a millimeter or so, it turns out experimentally that the liquid flows smoothly around the ball as it falls, with a flow pattern like:

If we knew mathematically precisely how the velocity in this flow pattern varied near the ball, we could find the total viscous force on the ball by finding the velocity gradient near each little area of the ball’s surface, and doing an integral. But actually this is quite difficult. It was done in the 1840’s by Sir George Gabriel Stokes. He found what has become known as *Stokes’ Law*: the drag force \(F\) on a sphere of radius *a* moving through a fluid of viscosity \(\eta\) at speed \(v\) is given by:

\[ F= 6\pi a \eta v\]

This drag force is directly proportional to the radius. That’s not obvious—one might have thought it would be proportional to the cross-section area, which would go as the square of the radius. The drag force is also directly proportional to the speed, not, for example to \(v^2\).

Understanding Stokes’ Law with Dimensional Analysis

Is there some way we could see the drag force must be proportional to the radius, and to the speed, without wading through all of Sir George’s mathematics? The answer is yes—by using **dimensions**.

- First we must ask: what can this drag force depend on? Obviously, it depends on the
*size*of the ball: let’s say the radius is*a*, having dimension*L*. - It must depend on the
*speed*\(v\), which has dimension \( LT^{-1} \) - Finally, it depends on the
*coefficient of viscosity \( \eta \)*which has dimensions \( ML^{-1}T^{-1} \).

The drag force *F* has dimensions \( [F] = MLT^{-2} \) what combination of \( [ \alpha ] = L, \, [\nu] = LT^{-1} \) and \([\eta] = ML^{-1}T^{-1} \) will give \( [F] = MLT^{-2} \)?

It’s easy to see immediately that \(F\) must depend linearly on \(\eta\), that’s the only way to balance the \(M\) term. Now let’s look at \( \dfrac{F}{\eta} \), which can only depend on \(a\) and \(v\).

\[ \left[\dfrac{F}{\eta} \right] = L^2T^{-1}. \]

The only possible way to get a function of \(a\), \(v\) having dimension \( L^2T^{-1} \) is to take the product \(av\).

So, the dimensional analysis establishes that the drag force is given by:

\[ F = C \alpha \eta \nu \]

where \(C\) is a constant that cannot be determined by dimensional considerations.

## Experimental Check

We can check this result by dropping small steel balls through glycerin. We choose glycerin because it has a very high viscosity, so the balls fall slowly enough for us to be able to time them (Figure \(\PageIndex{2}\).

One problem is that the viscosity of glycerin is *very* temperature dependent, being 1.49 Pa.sec at 20^{o} C, and 0.95 Pa.sec at 25^{o}C. We measured the temperature of our glycerin to be 23^{o}C, so we assumed its viscosity was 1.17 Pa.sec., just taking a linear interpolation. We used a ball of radius 1.2mm, weighing 0.05 grams. On dropping it through the glycerin, and allowing some distance for it to reach a steady speed, we found it fell 25 cm. in 11.1 seconds, a speed of 0.022 m sec^{-1}. So the drag force should be:

\[ F = 6 \pi \alpha \eta \nu \approx 18.8 \times 1.2 \times 10^{-3} \times 1.17 \times 2.2 \times 10^{-2} \approx 6 \times 10^{-4} N \]

If the ball is dropping at a steady speed, this force should just balance the weight of the ball. The mass is 0.05grams, which is 5.10^{-4 }N. But we should also have subtracted off a buoyancy force, which would get this closer to 4.10^{-4 }N. Since we think our measurements of radius, mass and speed were fairly accurate, the viscosity was evidently less than we thought. Bearing in mind that it drops by 10% for each one-degree rise in temperature, most likely it was not at a uniform temperature, or our measurement of temperature was inaccurate.

We checked the *dimensional* prediction by dropping a ball of exactly *twice* the radius. It fell in exactly *one-quarter* of the time. This confirms the correctness of the dimensional analysis, because once the ball has reached terminal velocity *v*_{term}, and therefore is no longer accelerating, it must feel zero net force. At this stage, the forces of viscous drag and weight must be in balance:

\[ C \alpha \eta \nu_{term} = \left( \dfrac{4}{3} \right) \pi \alpha^3 \rho g. \]

It follows that for two balls of the same density \(\rho\), after canceling *a* from each side, the ratio of their terminal velocities is the *square* of the ratio of their radii, a ball with radius \(2a\) will fall four times faster than a ball with radius \(a\). This is what we found experimentally.

Example \(\PageIndex{1}\)

Assuming the flow pattern in the diagram above has the same proportions for different radii (so for a larger radius ball it’s the same pattern magnified), how does the fluid *velocity gradient* near the “equator” of the ball change on going from a ball of radius *a* to one of radius 2*a*? (Assume the two balls are falling through the fluid at the same speed.) Argue that most of the viscous drag on the sphere takes place in a band surrounding the equator (so, a band shaped like the tropical zone on the earth). From this, make plausible that the total viscous drag will be proportional to the sphere’s radius, not to the square of the radius.

### Contributors

- Michael Fowler (Beams Professor, Department of Physics, University of Virginia)