Dimensional Analysis – Drag on a Sphere R. Shankar Subramanian
Dimensional Analysis – Drag on a Sphere R. Shankar Subramanian
Dimensional Analysis – Drag on a Sphere R. Shankar Subramanian
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
<str<strong>on</strong>g>Dimensi<strong>on</strong>al</str<strong>on</strong>g> <str<strong>on</strong>g>Analysis</str<strong>on</strong>g> – <str<strong>on</strong>g>Drag</str<strong>on</strong>g> <strong>on</strong> a <strong>Sphere</strong><br />
R. <strong>Shankar</strong> <strong>Subramanian</strong><br />
In a typical experiment, we look for how a dependent parameter varies as we change a variety of<br />
independent parameters. By independent parameters, we mean those that we can vary<br />
c<strong>on</strong>veniently, and the dependent parameter is the <strong>on</strong>e whose behavior we seek to establish.<br />
<str<strong>on</strong>g>Dimensi<strong>on</strong>al</str<strong>on</strong>g> analysis permits us to organize the process by which we vary the independent<br />
parameters. In fact, it helps us identify the true dependent and independent parameters in a<br />
situati<strong>on</strong>; in the process, the number of parameters that we must c<strong>on</strong>sider is minimized. But<br />
dimensi<strong>on</strong>al analysis is not foolproof – we must be careful in c<strong>on</strong>sidering all possible parameters<br />
that can affect the dependent parameter. If we omit a crucial parameter in making the list of<br />
independent parameters, dimensi<strong>on</strong>al analysis cannot help us find it. We always need comm<strong>on</strong><br />
sense and physical intuiti<strong>on</strong> in selecting the lists of parameters in a problem.<br />
Here, by c<strong>on</strong>sidering a simple example, I’ll show you how to use dimensi<strong>on</strong>al analysis. We aim<br />
to develop an organized way of examining how the drag <strong>on</strong> a sphere settling through a fluid<br />
varies with the relevant parameters. The first step is to identify the parameters <strong>on</strong> which the drag<br />
F D is likely to depend. These are the diameter D of the sphere, the density ρ and the viscosity<br />
μ of the fluid, and the settling velocity of the sphere V .<br />
Step 1<br />
Make a list of parameters and identify their dimensi<strong>on</strong>s using the fundamental dimensi<strong>on</strong>s of<br />
Mass ( M ) Length ( L ) and Time ( T ) . From now <strong>on</strong>, we’ll use the term “variables” to<br />
designate these parameters, but recognize that they are parameters in the usual sense of that term.<br />
Variable Symbol Dimensi<strong>on</strong><br />
<strong>Sphere</strong> Diameter<br />
D L<br />
Fluid Density ρ<br />
M<br />
3<br />
L<br />
Settling Velocity<br />
V<br />
L<br />
T<br />
Fluid Viscosity μ<br />
M<br />
LT<br />
<str<strong>on</strong>g>Drag</str<strong>on</strong>g> D F ML<br />
2<br />
T<br />
1
Step 2<br />
Identify the repeat variables. There should be as many repeat variables as there are dimensi<strong>on</strong>s.<br />
These are variables that will, in principle, appear in every dimensi<strong>on</strong>less group that we form.<br />
The requirement is that it should not be possible to form a dimensi<strong>on</strong>less group in the form of a<br />
product of powers of these variables. The simplest way to meet this requirement is to select the<br />
repeat variables so that each has a unique dimensi<strong>on</strong> in it. Also, it is helpful to choose repeat<br />
variables that are as simple as possible in their dimensi<strong>on</strong>s. With this in mind, we choose the<br />
first three, namely, the diameter, density, and velocity as our repeat variables, because there are<br />
three dimensi<strong>on</strong>s in this problem. You can see that density has Mass ( M ) appearing uniquely in<br />
it, and velocity has time ( T ) uniquely appearing in it, and the diameter has neither. Therefore, it<br />
is not possible to use these three variables to form a dimensi<strong>on</strong>less group.<br />
Step 3<br />
The number of dimensi<strong>on</strong>less groups is always equal to the number of variables minus the<br />
number of repeat variables. Therefore, we can expect to form two dimensi<strong>on</strong>less groups in this<br />
problem. The group involving the drag will be the dependent dimensi<strong>on</strong>less group and that<br />
involving the viscosity will be the independent dimensi<strong>on</strong>less group. Each group is obtained by<br />
forming a product of each repeat variable to an unknown power and then multiplying by <strong>on</strong>e of<br />
the remaining variables. If we call these two dimensi<strong>on</strong>less groups Π 1 and Π2<br />
, then we might<br />
define Π as follows.<br />
Π =<br />
1<br />
1<br />
a b c<br />
μ D ρ V<br />
To find the values of ab , ,andcthat<br />
would make this group dimensi<strong>on</strong>less, first we write out the<br />
equati<strong>on</strong> in terms of the dimensi<strong>on</strong>s of each side.<br />
b c<br />
0 0 0 M a ⎛M ⎞ ⎛ L⎞<br />
1+ b −+ 1 a− 3b+ c −− 1 c<br />
= × × × =<br />
3<br />
M LT L ⎜ ⎟ ⎜ ⎟ M L T<br />
LT ⎝ L ⎠ ⎝T⎠ Now, match the powers of M , L, and T<br />
the exp<strong>on</strong>ents ab , , and cas<br />
follows.<br />
1+ b = 0<br />
− 1+ a− 3b+<br />
c=<br />
0<br />
−1− c = 0<br />
<strong>on</strong> both sides. This yields three algebraic equati<strong>on</strong>s for<br />
The soluti<strong>on</strong> of these equati<strong>on</strong>s is straightforward, and the result is a = b= c=−<br />
1.<br />
Therefore,<br />
μ<br />
Π 1 =<br />
DρV 2
You can verify that the group in the right side of the above result is indeed dimensi<strong>on</strong>less. Now,<br />
there is nothing in dimensi<strong>on</strong>al analysis that prevents us from using powers of this group as a<br />
dimensi<strong>on</strong>less group in which these four variables appear. It is c<strong>on</strong>venti<strong>on</strong>al to invert this group<br />
and call it the Reynolds number.<br />
So, let us invent a new dimensi<strong>on</strong>less group to replace Π 1 .<br />
Reynolds Number Re =<br />
D V ρ<br />
μ<br />
We shall have more to say about the Reynolds number and its physical significance later in the<br />
course. Next, let us find the dimensi<strong>on</strong>less group that includes the dependent variable F D .<br />
Π = F D ρ V<br />
2<br />
D<br />
a b c<br />
Again, we write out the equati<strong>on</strong> in terms of the dimensi<strong>on</strong>s of each side.<br />
b c<br />
0 0 0 ML a ⎛M⎞ ⎛L⎞ 1+ b 1+ a− 3b+ c −2−c = × × × =<br />
2 3<br />
M LT L ⎜ ⎟ ⎜ ⎟ M L T<br />
T ⎝ L ⎠ ⎝T ⎠<br />
As before, matching the powers of M , L,andT the exp<strong>on</strong>ents ab , ,andc.<br />
1+ b = 0<br />
1+ a− 3b+<br />
c=<br />
0<br />
−2− c = 0<br />
The soluti<strong>on</strong> is a=− 2, b=− 1, c=−2.<br />
Therefore,<br />
Π =<br />
F<br />
DV ρ<br />
2 2<br />
D<br />
2<br />
<strong>on</strong> both sides yields three algebraic equati<strong>on</strong>s for<br />
In practice, a dimensi<strong>on</strong>less drag known as the <str<strong>on</strong>g>Drag</str<strong>on</strong>g> Coefficient, C D , is defined as follows.<br />
C<br />
D<br />
=<br />
8FD DV<br />
2 2<br />
π ρ<br />
The reas<strong>on</strong> for the multiplicative factor ( 8/π ) is that the drag coefficient is defined as the drag<br />
divided by the product of the projected area of the sphere and the velocity head. Thus,<br />
3
C<br />
D<br />
FD 8FD<br />
= = 2 2<br />
π 2 1 2<br />
D × ρV<br />
π DV ρ<br />
4 2<br />
Therefore, dimensi<strong>on</strong>al analysis tells us that the drag coefficient is a universal functi<strong>on</strong> of<br />
Reynolds number, regardless of the choice of fluid, sphere diameter or the settling velocity.<br />
C = φ<br />
D<br />
( Re)<br />
The nature of the functi<strong>on</strong> φ has to be established by experiments or theory as appropriate. The<br />
universality means that when we plot the drag coefficient against the Reynolds number, we<br />
obtain a single curve, regardless of the choice of sphere diameter, fluid, and the velocity at which<br />
the sphere moves through the fluid. This is a remarkable result, because it permits us to use any<br />
fluid(s) of our choice to carry out the experiments to find the universal relati<strong>on</strong>ship between<br />
these two dimensi<strong>on</strong>less variables, and <strong>on</strong>ce determined, use it to infer the drag <strong>on</strong> any sphere<br />
settling through any fluid at any velocity. You can see how powerful a tool this is.<br />
Figure 7.16 (b) in Fluid Mechanics, Sixth Editi<strong>on</strong>, by F.M. White shows the universal curve that<br />
results for the drag <strong>on</strong> a sphere as a functi<strong>on</strong> of the Reynolds number.<br />
4